exaPD 초고속 다원소 상도 구축을 위한 대규모 병렬 워크플로

초록

exaPD는 LAMMPS 기반 MD·MC 시뮬레이션을 활용해 온도·조성 메쉬 전역에서 자유에너지를 동시에 계산하고, Parsl으로 전 세계적인 작업을 관리해 거의 이상적인 확장성을 제공한다. 계산된 자유에너지는 PYCALPHAD와 연계해 CALPHAD 모델링으로 전이상도까지 자동 생성한다.

상세 분석

exaPD는 전통적인 열역학적 적분(TI) 방법을 기반으로 하면서, 고체·액체·고체용액 등 다양한 상에 대해 각각 최적화된 기준계(Einstein crystal, Uhlenbeck‑Ford model 등)를 자동 선택한다. 특히 고체용액의 경우 반그랜드캐노니컬 집합(SGCE)을 구현해 원자 교환 MC와 MD를 교차 수행함으로써 구성 엔트로피를 정확히 샘플링한다. 워크플로의 핵심은 Parsl 기반의 글로벌 컨트롤러로, 수천에서 수만 개에 이르는 독립 MD/MC 작업을 비동기식으로 스케줄링하고, 입력·출력 파일을 JSON 형식으로 정의해 사용자 정의가 용이하도록 설계되었다. 이 구조는 작업 간 의존성을 최소화해 통신 오버헤드가 거의 없으며, 실험적으로 10 000 코어 이상에서 90 % 이상 효율을 유지하는 ‘거의 이상적인 확장성’을 입증한다. 또한 exaPD는 전통적인 EAM, FS 포텐셜뿐 아니라 신경망 포텐셜(NNP)까지 지원해, 근사 ab‑initio 정확도를 대규모 시뮬레이션에 적용할 수 있다. 검증 사례로 Cu‑Zr 시스템에서 EAM‑FS 포텐셜을, La‑Si‑P 시스템에서 최신 NNP를 사용해 자유에너지 곡선을 재현했으며, 실험 데이터와의 차이가 0.02 eV 이하로 매우 작은 오차를 보였다. 마지막 단계에서는 PYCALPHAD를 통해 얻은 자유에너지 데이터를 CALPHAD 모델에 피팅하고, 공통 접선(congruent tangent) 방법으로 상평형 조성을 자동으로 도출한다. 전체 파이프라인이 하나의 파이썬 스크립트로 통합돼, 사용자는 포텐셜 파일과 메쉬 정의만 제공하면 복잡한 다원소 상도를 몇 시간 내에 생성할 수 있다.