궤적을 생성하는 것에서 ensemble 평균과 유효 표본 수를 가진 측정으로 넘어가기
원자 궤적 영상은 그럴듯해 보여도 thermostat이 진동을 왜곡하거나 평형화 구간이 평균에 섞이면 정량 물성이 틀립니다. 연속 frame은 독립 표본도 아닙니다.
짧고 비싼 AIMD 궤적에서 올바른 ensemble 평균과 신뢰구간을 어떻게 얻을까요?
이번 장의 답질문에 맞는 ensemble과 thermostat/barostat의 결합시간을 선택하고 가장 빠른 진동을 충분히 분해하는 timestep을 씁니다. 평형화와 production을 분리하고 autocorrelation time보다 긴 block으로 오차와 유효 표본 수를 평가합니다.
01NVE·NVT·NPT 중 무엇을 모사하는가
| 앙상블 | 고정량 | 주요 용도 |
|---|---|---|
| NVE | N,V,E | 적분·SCF energy conservation 검사, dynamics |
| NVT | N,V,T | 고정 밀도의 온도 평형 구조·액체 |
| NPT | N,P,T | 밀도·부피 평형, 높은 비용·cell fluctuation |
thermostat은 열 reservoir와의 결합을 모사합니다. 결합이 너무 강하면 실제 진동 spectrum과 확산을 왜곡하고, 너무 약하면 짧은 AIMD 안에 목표 온도에 도달하지 못합니다. 구조 평균과 동역학 물성의 설정 요구가 다릅니다.
02가장 빠른 진동과 timestep
time step은 가장 높은 진동주기보다 충분히 작아야 합니다. H-X stretch가 있으면 heavy-atom 고체보다 더 작은 step이 필요합니다. step을 절반으로 줄여 NVE energy drift·온도·구조 분포가 유지되는지 검사합니다.
SCF를 느슨하게 하면 힘에 잡음과 비보존 성분이 생겨 thermostat이 이를 열처럼 흡수할 수 있습니다. NVE 짧은 test로 timestep과 SCF 오차를 분리합니다.
03평형화·autocorrelation·block average
초기 구조와 속도에서 ensemble 정상분포로 가는 구간은 equilibration으로 버리고, 이후 production에서 평균을 계산합니다. 에너지·온도·밀도뿐 아니라 관심 구조 order parameter의 시간 series가 stationary한지 봅니다.
block 길이를 늘렸을 때 block mean 분산이 plateau에 도달하는지 확인합니다. 서로 다른 초기속도의 독립 trajectory는 느린 mode와 trapped state를 감지하는 데 한 긴 궤적보다 유용할 수 있습니다.
04이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- 앙상블ensemble
- 정해진 거시 제약을 공유하는 가능한 미시상태의 확률 집합입니다.
- thermostat
- 계와 열 reservoir의 에너지 교환을 모사해 목표 온도 분포를 만드는 알고리즘입니다.
- 평형화equilibration
- 초기조건의 기억이 사라지고 목표 ensemble의 정상분포에 도달하는 과정입니다.
- production run
- 평형화 뒤 실제 통계량을 수집하는 궤적 구간입니다.
- autocorrelation time
- 시계열 값이 자신의 과거와 유의하게 상관되는 시간척도입니다.
- block average
- 상관된 시계열을 긴 구간별 평균으로 묶어 평균의 통계오차를 추정하는 방법입니다.
05참고자료와 더 읽을거리
- Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
- Car & Parrinello, Unified Approach for Molecular Dynamics and Density-Functional Theory (1985) — 전자구조 계산과 원자동역학을 결합한 ab initio molecular dynamics의 고전적 원 논문입니다.
- Togo & Tanaka, First principles phonon calculations in materials science (2015) — 힘상수·phonon·열역학량 계산을 실무 관점에서 정리한 참고문헌입니다.