전자 바닥상태 에너지의 기울기를 따라 안정 구조를 찾고 잔류력의 의미를 읽기
실험 구조를 그대로 계산하면 원자마다 힘이 남을 수 있습니다. DFT가 선택한 포텐셜 에너지면에서 원자가 움직이려는 방향을 따라가야 일관된 에너지와 물성을 얻습니다.
힘을 어떻게 계산하고 어떤 조건에서 구조가 충분히 최적화됐다고 판단할까요?
이번 장의 답원자 힘은 자기일관 총에너지의 좌표 미분이며, 완전한 기저에서는 Hellmann-Feynman 형태로 얻습니다. 유한 기저·PAW·SCF 오차를 통제하고 최대 힘·에너지 변화·대칭·불안정 모드를 함께 확인해야 합니다.
01Hellmann-Feynman 힘과 보정
정확한 고유상태에서 Hamiltonian의 매개변수 미분에 대한 에너지 미분은 Hamiltonian 미분의 기대값으로 계산됩니다. 핵 좌표 미분에 적용하면 전자 밀도와 이온 사이 힘을 효율적으로 얻습니다.
평면파 기저는 원자 위치에 직접 붙지 않아 위치 Pulay force가 작다는 장점이 있지만 PAW projector·augmentation와 유한 SCF 오차가 남습니다. 코드가 제공하는 힘 공식을 믿되 cutoff와 SCF에 대해 힘 자체를 검증합니다.
02steepest descent에서 quasi-Newton까지
힘 방향으로 조금씩 이동하는 steepest descent는 안정하지만 좁고 굽은 골짜기에서 느립니다. conjugate gradient는 이전 방향을 활용하고, BFGS 같은 quasi-Newton 방법은 힘 변화로 Hessian의 역을 근사해 효율적인 변위를 예측합니다.
한 번에 너무 큰 원자 이동은 SCF 실패와 원자 충돌을 만들 수 있습니다. 반대로 step이 지나치게 작으면 수치 잡음 아래에서 정체됩니다. preconditioner와 trust radius는 결합 강도와 최대 이동량을 조절합니다.
03힘 0 하나보다 넓은 종료 조건
- 최대 원자 힘과 RMS force가 허용치 아래인지
- 마지막 여러 step의 에너지 변화와 원자 이동이 작은지
- 고정 원자·대칭 제약이 의도한 것인지
- 셀까지 최적화할 경우 잔류 응력이 작은지
- phonon 또는 작은 변위에서 허수 모드가 없는지
가벼운 원자의 진동·전이상태·phonon은 일반 구조에너지보다 엄격한 힘 기준이 필요합니다. 반응물과 생성물은 같은 수준으로 최적화해야 장벽과 반응에너지에 편향이 생기지 않습니다.
04이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- Hellmann-Feynman 정리
- 정확한 고유상태 에너지의 매개변수 미분을 Hamiltonian 미분의 기대값으로 바꾸는 정리입니다.
- Pulay force
- 원자 위치에 따라 변하는 불완전 기저 때문에 에너지 미분에 추가되는 힘 항입니다.
- Hessian
- 에너지를 원자 좌표로 두 번 미분한 곡률 행렬입니다.
- quasi-Newton 방법
- 힘 변화로 Hessian 또는 그 역을 갱신해 최적화 step을 정하는 방법입니다.
- trust radius
- 국소 곡률 근사가 신뢰된다고 보는 한 step의 최대 이동 범위입니다.
- 잔류력residual force
- 최적화 종료 구조에서 아직 0이 아닌 원자별 힘입니다.
05참고자료와 더 읽을거리
- Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
- Blöchl, Projector augmented-wave method (1994) — 평면파 효율성과 핵 주변의 전전자 정보를 잇는 PAW 방법의 원 논문입니다.
- Kohn & Sham, Self-Consistent Equations Including Exchange and Correlation Effects (1965) — 상호작용 전자계를 계산 가능한 유효 독립전자 방정식으로 바꾼 원 논문입니다.