초고속 빌리언 바디 분자동역학: 72랙 블루진/​P에서 1000억 입자 규모의 초대형 유리 전이 시뮬레이션

초록

본 논문은 새로 개발한 Billions‑Body Molecular Dynamics (BBMD) 코드를 IBM BlueGene/P 슈퍼컴퓨터 72랙(≈2.95×10⁵ 코어)까지 확장하여 1000억 입자 시스템을 89 % 효율로 실행하고, 1 억 입자 규모에서는 0.14 초 이하의 시간당 한 스텝을 달성함을 보고한다. 최적화된 도메인 분할, 양방향 연결 리스트, 6방향 통신 겹침 기법 등 소프트웨어 설계와 BlueGene 전용 sqrt 가속을 상세히 설명하고, 1 억 입자 이진 혼합 소프트 구(soft‑sphere) 시스템의 액체‑유리 전이를 온도 구간 T∈

상세 분석

BBMD 코드는 대규모 구조유리와 비정질 물질의 원자 동역학을 목표로 설계된 전용 Lennard‑Jones/soft‑sphere 시뮬레이터이다. 핵심은 3차원 공간을 r_c×r_c×r_c 크기의 셀로 균등 분할하고, 각 셀에 입자 정보를 담은 양방향 연결 리스트를 두어 메모리 이동 없이 포인터 교환만으로 입자 이동을 구현한다. 이는 O(N) 복잡도의 인접 탐색을 가능하게 하며, 전통적인 neighbor list(목록) 방식을 배제해 메모리 사용량을 최소화한다.
병렬화는 도메인 분할 기반 MPI 구현으로, 각 프로세스는 26개의 최근접 이웃과 통신한다. 그러나 BBMD는 통신 창을 확장해 실제 교환 횟수를 6으로 축소하고, 입자 태깅(질량 부호 반전)과 정렬을 온‑더‑플라이로 수행해 통신 버퍼를 연속 메모리 블록으로 만든다. 이와 동시에 계산 단계와 통신을 겹쳐(overlap) MPI Waitall 시간을 크게 감소시켰다.
BlueGene/P 특화 최적화로는 하드웨어 전용 frsqrte 함수와 두 차례 Newton‑Raphson 반복을 이용해 sqrt 연산을 약 7 % 가속하였다. 이러한 저수준 최적화와 고수준 알고리즘 설계가 결합돼 1 억 입자에서 0.14 초 이하, 100 억 입자에서는 2 초 수준의 스텝 타임을 달성했으며, 72랙 전체에서 89 % 확장 효율을 기록했다.
강화된 약한 스케일링 실험에서는 입자 수를 1 억→10 억→100 억으로 3배씩 증가시켜도 실행 시간은 거의 선형(O(N))으로 증가했으며, 계산 대비 통신 비중이 1 % 이하로 감소해 대규모 시스템에서도 통신 병목이 사라졌다.
과학적 적용 사례로, 1 억 입자 이진 혼합 소프트 구 시스템을 온도 2→10⁻⁴까지 서서히 냉각하며 액체‑유리 전이를 관찰했다. 시스템은 m=1, σ₁₁=1, σ₁₂=0.8, σ₂₂=0.88, ε₁₁=1, ε₁₂=1.5, ε₂₂=0.5의 파라미터를 갖고, 밀도 ρ=1.2, 온도 T=0.5에서 초기 액체 상태를 유지한다. 온도 감소에 따라 동적 구조인자와 평균 제곱 변위가 급격히 정체되는 유리 전이 특성을 보였으며, 이는 기존 10⁷ 입자 규모 시뮬레이션 대비 100배 이상 큰 시뮬레이션 볼륨에서 동일 현상이 재현된 최초 사례이다. 이러한 결과는 BBMD가 실험적으로 접근하기 어려운 마이크로미터 스케일 비정질 물질의 동역학을 직접 탐구할 수 있음을 시사한다.