AIMD: DFT 힘으로 원자를 움직이기 | 조우현의 연구 노트

각 시간 step의 전자구조 계산을 핵의 운동방정식과 결합하는 Born-Oppenheimer MD

DFT · PART 6 · 정적인 바닥상태를 온도·반응·시간의 세계로 확장하는 동안05 / 08

phonon은 평형점 주변의 작은 진동을 선형화하지만 액체, 확산, 결합 교환과 큰 구조변화는 하나의 조화 우물로 설명하기 어렵습니다. 원자를 실제 시간 순서로 움직여야 합니다.

이번 장의 답Born-Oppenheimer AIMD는 각 핵 배치에서 전자 SCF를 충분히 수렴해 힘을 얻고 핵의 고전 운동방정식을 적분합니다. 화학 반응을 자연스럽게 포함하지만 수십~수백 원자와 ps 규모의 비용, 전자·시간 step 수렴을 감수합니다.

SCF INSIDE EVERY TIME STEP

01핵 좌표$\mathbf R(t)$

02DFT SCF·힘$\mathbf F=-\nabla E$

03시간 적분$\mathbf R(t+\Delta t)$

$M_I\ddot{\mathbf R}_I=\mathbf F_I^{DFT}$

전자 바닥상태 계산과 핵 운동 적분을 매 step 반복해 반응 가능한 제일원리 궤적을 만듭니다.

01Born-Oppenheimer AIMD의 한 step

time step loop

핵은 보통 고전적으로 다룹니다. H의 tunneling·영점 운동처럼 nuclear quantum effect가 중요하면 path-integral MD 같은 확장이 필요합니다.

Car-Parrinello MD는 Kohn-Sham 오비탈에 가상의 질량을 주고 핵과 함께 동역학적으로 진화시켜 매 step 완전한 SCF를 피합니다. 전자 가상운동과 핵 진동의 시간척도를 분리해야 하며 fictitious kinetic energy를 감시합니다.

접근	전자 처리	핵심 검사
BO-AIMD	매 step SCF ground state	SCF residual·energy drift
Car-Parrinello	가상 전자 동역학	fictitious mass·adiabatic separation

DFT 계산이 매 step 필요해 보통 femtosecond step으로 picosecond~수십 ps를 샘플링합니다. 장벽이 큰 확산·핵생성·느린 재배열은 궤적에서 한 번도 일어나지 않을 수 있습니다. 관찰되지 않았다는 것이 불가능하다는 뜻은 아닙니다.

희귀사건은 metadynamics, umbrella sampling, accelerated dynamics 또는 NEB/TST와 결합합니다. 더 큰 계·긴 시간은 검증된 MLFF로 확장하되 DFT 학습영역 밖을 감지해야 합니다.

본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.

Car & Parrinello, Unified Approach for Molecular Dynamics and Density-Functional Theory (1985) — 전자구조 계산과 원자동역학을 결합한 ab initio molecular dynamics의 고전적 원 논문입니다.
Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
Kohn & Sham, Self-Consistent Equations Including Exchange and Correlation Effects (1965) — 상호작용 전자계를 계산 가능한 유효 독립전자 방정식으로 바꾼 원 논문입니다.