파르포옴: 트위스트 2차원 재료의 모아레 포논 계산을 위한 병렬 툴킷

초록

**
본 논문은 트위스트된 2차원 재료(모아레 시스템)의 대규모 원자 배열에 대해 포논 계산을 효율적으로 수행할 수 있는 오픈소스 패키지 PARPHOM을 소개한다. LAMMPS 기반의 고전력장으로부터 힘 상수를 병렬로 생성하고, ScaLAPACK을 이용해 동적 행렬을 대규모 병렬 대각화한다. 또한 포논 밀도 상태, 밴드 구조, 온도 의존성, 그리고 모드의 챈리티까지 분석할 수 있는 전용 유틸리티를 제공한다.

**

상세 분석

**
PARPHOM은 트위스트 2D 재료에서 발생하는 수천~수만 원자 규모의 모아레 초격자에 대한 포논 계산의 근본적인 병목을 해소한다. 핵심 아이디어는 두 단계로 나뉜다. 첫 번째는 LAMMPS의 클래식 포스 필드를 이용해 원자들을 작은 변위(Δr = 10⁻⁴ Å)로 이동시키고, 중앙 차분을 통해 힘 상수 Φᵅᵝ(jl, j′l′)를 직접 계산한다. 이 과정은 각 원자·방향마다 독립적인 힘 계산을 요구하므로 “embarrassingly parallel” 구조를 갖는다. 구현에서는 MPI와 parallel‑h5py를 이용해 원자 집합을 Nₚ 프로세스에 균등 분배하고, 각 프로세스가 자신에게 할당된 원자들의 6 n(3 방향 × 2 정방향) 변위를 수행한다. 기존 코드가 대칭성을 활용해 변위 수를 감소시키는 반면, PARPHOM은 대규모 병렬성을 활용해 모든 변위를 수행하고, 이후 반복적인 대칭 보정(전치 대칭, 음향 합 규칙 등)을 통해 힘 상수 행렬을 정규화한다. 이는 대칭 연산 자체가 전체 변위보다 비용이 크지 않은 경우에 특히 효율적이다.

두 번째 단계는 생성된 힘 상수 행렬을 이용해 동적 행렬 Dᵅᵝ(jj′, q)를 구성하고, ScaLAPACK의 PZHEEVX 루틴으로 블록‑사이클 방식으로 분산된 행렬을 직접 대각화한다. 여기서 중요한 점은 포논 모드의 그룹 속도(v_qν) 계산을 위해 동적 행렬의 q‑미분을 두 가지 방법(해석적 도함수와 중앙 차분)으로 제공한다는 것이다. 특히 밴드 교차점에서의 속도 계산이 불안정해지는 문제를 해소하기 위해, 파라미터화된 형태의 ⟨Ψ|∂D/∂q|Ψ⟩를 이용한 표현식을 구현하였다.

또한 PARPHOM은 온도 의존적인 포논 밀도 상태(DOS)를 모드‑투사된 속도 자동상관 함수(VACF)를 통해 비교적 정확하게 추정한다. 이때 사용자는 Gaussian 혹은 Lorentzian 스미어링, 혹은 선형 삼각법을 선택해 DOS를 출력할 수 있다. 챈리티 분석 기능은 각 모드의 회전 대칭성을 평가해, 트위스트 각도에 따른 비대칭 모드의 출현을 정량화한다.

성능 측면에서 저자들은 PHONOPY와 PhonoLAMMPS와 비교해, 힘 상수 생성 단계에서 10 ~ 100배 이상의 스케일업을 달성했으며, 전체 파이프라인(힘 상수 → 동적 행렬 → 대각화)까지도 수백 코어 규모에서 효율적으로 동작함을 시연했다. 메모리 사용량 역시 HDF5 기반의 병렬 입출력으로 관리되어, 수십 GB 규모의 힘 상수 행렬도 일반 HPC 클러스터에서 처리 가능하도록 설계되었다.

마지막으로, PARPHOM은 GitHub에 오픈소스로 공개돼 GPLv3 라이선스를 채택하고 있으며, Python 인터페이스(pyphutil)를 통해 후처리(밴드 플롯, DOS, 속도 시각화)를 간편히 수행할 수 있다. 설치 스크립트와 예제 Makefile(intel, aocc)도 제공돼, 다양한 컴파일러와 MPI 환경에 손쉽게 적용할 수 있다.

**