ms2 2.0: 하이브리드 병렬화와 질량포텐셜 파생량 구현을 통한 열역학·전송 특성 확장

초록

ms2 시뮬레이션 툴의 2.0 버전은 MPI와 OpenMP를 결합한 하이브리드 병렬화로 확장성을 높였으며, Lustig의 질량포텐셜(Massieu) 파생량 공식을 도입해 단일 NVT 시뮬레이션에서 다양한 열역학 파생량을 동시에 샘플링한다. 또한 Green‑Kubo 형식을 확장해 전기 전도도와 체류 시간을 계산하고, Ewald summation을 구현해 전하를 가진 시스템도 다룰 수 있게 하였다. 방사형 분포 함수(RDF) 샘플링 기능도 추가되었다.

상세 분석

본 논문은 기존 ms2(2011)에서 제시된 분자 동역학(MD) 기반 열역학 특성 계산 프레임워크를 크게 세 가지 측면에서 개선하였다. 첫째, MPI만 사용하던 병렬 구조를 OpenMP와 결합한 하이브리드 형태로 전환함으로써, 1 000 ~ 2 000개의 입자를 다루는 경우에도 2 048코어 규모에서 20 % 정도의 추가 속도 향상을 달성하였다. 이는 MPI 통신이 병목이 되는 구간에서도 각 프로세스가 충분한 연산량을 보유하도록 스레드 수준의 병렬화를 적용한 결과이다. 특히, 힘 계산에서 발생하는 레이스 컨디션을 회피하기 위해 상호작용별 힘을 임시 리스트에 저장한 뒤 후처리 단계에서 합산하는 전략을 사용했으며, 이는 원자적 업데이트나 크리티컬 섹션 사용에 비해 오버헤드를 크게 감소시켰다.

둘째, Lustig이 제시한 Massieu 포텐셜 F/T(N,V,1/T)의 파생량 A_mn을 직접 샘플링하는 기능을 구현하였다. 이는 A_10, A_01, A_20, A_11, A_02 등 최대 2차까지의 파생량을 계산할 수 있게 하며, 각 파생량은 내부 에너지 U에 대한 부피 미분(∂U/∂V, ∂²U/∂V²)과 연계된다. 논문은 Lennard‑Jones 및 Coulomb 포텐셜에 대해 이러한 미분식을 사전에 분석하고, 장거리 보정(LRC)을 포함한 정확한 표현식을 제공한다. 결과적으로, 단일 NVT 시뮬레이션에서 자유에너지, 엔탈피, 열용량, 압축성 등 다양한 열역학 특성을 동시에 얻을 수 있어, 실험 데이터와 결합한 고성능 상태 방정식 개발에 유리하다.

셋째, 전송 특성 계산을 위한 Green‑Kubo 공식이 전기 전도도와 체류 시간(residence time)까지 확장되었다. 전기 전도도 σ는 전류 플럭스 j_e(t)의 자기상관함수 적분으로 정의되며, 이때 전하를 가진 입자만을 대상으로 하여 계산 효율을 높였다. 체류 시간은 특정 거리 r_ij 내에 머무는 시간의 평균값으로, Heaviside 함수와 입자 간 거리 판정을 통해 자동으로 샘플링한다. 또한, 기존에 전기 중성 시스템에만 적용 가능했던 반응장 방법을 대체해 Ewald summation을 도입함으로써 이온 용액 및 전하를 가진 분자 시스템도 정확히 시뮬레이션할 수 있게 되었다. 다만, 현재 Ewald과 Massieu 파생량, 하이브리드 병렬화는 동시에 사용할 수 없는 제한이 존재한다는 점을 명시한다.

성능 평가에서는 동일 조건 하에 GROMACS와 비교했을 때, ms2는 반응장(RF) 모드에서는 비슷하거나 약간 뒤처지지만, Ewald(EW) 모드에서는 2배 이상 빠른 속도를 보였다. 또한, 메모리 사용량과 입자 수에 따른 스케일링 그래프를 제시해 하이브리드 병렬화의 효용성을 실증하였다.

전반적으로 ms2 2.0은 고성능 컴퓨팅 환경에서 중소 규모(10³ ~ 10⁴ 입자) 시스템의 열역학·전송 특성을 포괄적으로 다루는 도구로서, 오픈소스 배포와 상세 매뉴얼 제공을 통해 학계·산업계의 활용 가능성을 크게 확대하였다.