DFT 힘의 선형응답에서 force constant·dynamical matrix·진동 분산을 얻기
0 K 최적화 구조에서도 원자는 정지한 점입자가 아니라 진동합니다. 작은 변위에 대한 힘의 변화를 알면 결정 전체의 집단 진동과 열역학의 첫 보정을 만들 수 있습니다.
원자별 DFT 힘에서 파동벡터별 phonon 주파수와 동역학적 안정성을 어떻게 얻을까요?
이번 장의 답평형점 주변 에너지를 원자 변위의 2차까지 전개해 force constant를 구하고 질량으로 정규화한 dynamical matrix를 대각화합니다. 고유값은 $\omega^2$, 고유벡터는 진동모드이며 음의 $\omega^2$는 구조 불안정 방향을 뜻합니다.
01평형점 주변의 조화근사
원자 변위 $u$가 작으면 에너지를 Taylor 전개할 수 있습니다. 완화된 평형점에서는 일차항인 힘이 0이고, 이차 force constant 행렬이 복원력을 정합니다.
조화근사는 충분히 낮은 온도와 작은 진동에서 유효합니다. 열팽창·phonon 수명·강한 soft mode에는 3차 이상 비조화 force constant가 필요합니다.
02finite displacement와 DFPT
finite displacement 방법은 초격자의 대칭적으로 독립인 원자를 조금 움직이고 DFT 힘을 계산해 $\Phi=-\partial F/\partial u$를 유한차분합니다. 구현이 직관적이고 다양한 코드에 적용되지만 큰 셀과 엄격한 힘 정확도가 필요합니다.
density-functional perturbation theory는 작은 주기 perturbation에 대한 전자 밀도 응답을 선형으로 풀어 q점별 dynamical matrix를 직접 구합니다. 두 방법 모두 q-grid·supercell·cutoff·k-point를 수렴시켜야 합니다.
| 방법 | 계산 단위 | 주요 수렴 |
|---|---|---|
| Finite displacement | 초격자 변위와 힘 | 변위 크기·초격자·force |
| DFPT | q별 선형응답 | q-grid·전자 SCF·response |
| 둘의 공통 | dynamical matrix | ASR·k-grid·구조 완화 |
03분산·DOS·허수모드 읽기
원시셀에 s개 원자가 있으면 각 q에 3s개 branch가 있고, q→0에서 세 acoustic branch의 주파수가 0으로 갑니다. 나머지는 optical branch입니다. acoustic sum rule은 전체 평행이동에 복원력이 없어야 한다는 조건입니다.
plot에서 음수로 그리는 imaginary frequency는 $\omega^2<0$인 음의 곡률입니다. 수치 잡음·불충분 수렴인지 실제 구조 불안정인지 q-grid와 더 큰 셀, mode 방향 변위 후 재최적화로 구분합니다.
04이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- phonon
- 조화격자 진동의 양자이자 결정의 집단적 정상모드입니다.
- force constant
- 원자 변위에 대한 다른 원자 힘의 선형 변화율인 에너지 2차 미분입니다.
- dynamical matrix
- force constant를 원자 질량으로 정규화하고 q공간으로 Fourier 변환한 행렬입니다.
- acoustic mode
- 장파장 한계에서 원자들이 함께 움직여 주파수가 0으로 가는 진동 branch입니다.
- optical mode
- 원시셀 안 서로 다른 원자가 상대적으로 움직이는 유한주파수 branch입니다.
- 허수모드imaginary mode
- 에너지 곡률이 음수인 구조 불안정 방향으로 $\omega^2<0$인 모드입니다.
05참고자료와 더 읽을거리
- Togo & Tanaka, First principles phonon calculations in materials science (2015) — 힘상수·phonon·열역학량 계산을 실무 관점에서 정리한 참고문헌입니다.
- Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
- Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.