빠른 전자와 느린 핵을 분리해 구조·힘·반응 경로의 공통 지도를 만들기
원자핵과 전자를 모두 동시에 양자적으로 추적하면 문제는 더 커집니다. 그러나 양성자는 전자보다 약 1836배 무겁고, 같은 에너지 척도에서 훨씬 느리게 움직입니다.
핵과 전자의 시간척도 차이를 이용하면 재료 계산은 어떻게 단순해질까요?
이번 장의 답핵 좌표를 순간적으로 고정한 채 전자 바닥상태를 풀면 에너지가 핵 좌표의 함수가 됩니다. 그 포텐셜 에너지면의 기울기는 힘이고, 최소점은 안정 구조이며, 안장점은 전이상태입니다.
01질량 차이가 만드는 시간척도 분리
핵은 전자보다 훨씬 무겁기 때문에 전자가 핵 배치의 변화에 빠르게 적응한다고 가정할 수 있습니다. 전체 파동함수를 전자 부분과 핵 부분의 곱으로 근사하고, 먼저 고정된 핵 좌표에서 전자 고유문제를 풉니다.
이것이 Born-Oppenheimer 근사입니다. 제일원리 계산도 근사를 사용합니다. 중요한 것은 근사의 물리적 조건을 알고, 전자 여기와 핵 운동이 강하게 결합하는 상황에서는 한계를 의식하는 것입니다.
02포텐셜 에너지면에서 구조와 반응 읽기
핵 좌표가 많으면 $E_0(\{\mathbf R_I\})$는 고차원 지형이 됩니다. 골짜기의 국소 최소점은 준안정 구조, 가장 낮은 최소점은 주어진 조건의 바닥상태 구조 후보입니다. 두 골짜기 사이 가장 낮은 고개는 1차 안장점이며 전이상태와 연결됩니다.
- 기울기를 따라 내려가 구조 최적화
- 셀 변형률 미분으로 응력과 격자 최적화
- 최소에너지경로를 찾아 NEB와 반응 장벽
- 힘으로 뉴턴 방정식을 적분해 Born-Oppenheimer AIMD
DFT 총에너지는 0 K의 전자 바닥상태 지형을 주는 경우가 많습니다. 실제 온도에서 안정성을 논하려면 진동·배치·자기 엔트로피와 압력 항을 포함한 자유에너지로 확장해야 합니다.
03근사가 흔들리는 경우
금속에서 빠른 비평형 전자와 강하게 결합하는 원자 운동, 광여기 상태, 원뿔교차 부근처럼 여러 전자상태가 가까우면 전자가 항상 하나의 바닥상태를 즉시 따른다는 가정이 깨질 수 있습니다. 이런 경우 비단열 결합이나 시간의존 전자구조 방법이 필요합니다.
그러나 대부분의 평형 구조·결합에너지·표면·결함과 열적 원자운동에는 Born-Oppenheimer 그림이 매우 유용합니다. 이후 모든 계산에서 ‘고정된 핵 배치의 전자 문제를 풀고 핵을 갱신한다’는 공통 구조를 계속 사용합니다.
04이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- Born-Oppenheimer 근사
- 핵과 전자의 질량·시간척도 차이를 이용해 전자 문제를 고정된 핵 좌표에서 먼저 푸는 근사입니다.
- 포텐셜 에너지면potential energy surface
- 모든 핵 좌표에 대해 전자 바닥상태 에너지를 대응시킨 고차원 함수입니다.
- 국소 최소점local minimum
- 작은 모든 변위에 대해 에너지가 증가하지만 전체에서 가장 낮다고 보장되지는 않는 구조입니다.
- 안장점saddle point
- 반응좌표 한 방향으로는 최대이고 나머지 방향으로는 최소인 전이상태 후보입니다.
- 단열adiabatic
- 핵이 움직이는 동안 전자가 선택한 전자상태를 계속 따라간다는 분리 그림입니다.
- 비단열 결합nonadiabatic coupling
- 핵 운동과 서로 다른 전자상태 사이 전이가 무시될 수 없을 때의 결합입니다.
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의 고전적 원 논문입니다.
- Kohn, Nobel Lecture: Electronic structure of matter (1999) — 파동함수에서 밀도 범함수로 이어진 전자구조 이론의 맥락을 정리합니다.