서로 다른 DFT 출력이 무엇을 보장하고 어디서 과해석되기 쉬운지 구분하기

DFT · PART 2 · 다체 파동함수를 전자 밀도와 유효 독립전자로 바꾸는 동안06 / 06

계산이 끝나면 총에너지, 수백 개의 고유값, 전하 밀도와 여러 파생량이 한꺼번에 나옵니다. 가장 보기 쉬운 숫자를 물리적 관측량으로 곧바로 부르면 오류가 시작됩니다.

DFT 출력 가운데 무엇이 바닥상태 이론의 직접 결과이고 무엇은 보조계의 해석 도구일까요?

이번 장의 답총에너지는 점유 Kohn-Sham 고유값의 단순 합이 아니며 이중계수 보정이 필요합니다. 정확한 DFT의 밀도와 그 범함수는 바닥상태량을 주지만 일반 Kohn-Sham 고유값 전체가 실제 여기 에너지는 아닙니다.

01왜 고유값을 더하면 총에너지가 아닌가

Kohn-Sham 고유값은 운동에너지와 유효 퍼텐셜의 기대값을 포함합니다. 이를 점유 상태에 대해 더하면 Hartree와 교환-상관 퍼텐셜 기여가 이미 들어 있으므로, 에너지 범함수의 항을 그대로 더할 경우 전자 상호작용이 중복 계수됩니다.

$$ E=\sum_i f_i\varepsilon_i-E_H[n]-\int v_{xc}n\,d\mathbf r+E_{xc}[n]+E_{NN} $$
국소 Kohn-Sham 형태의 개념적 총에너지 식입니다. 구현과 PAW 보정에 따라 세부 표현은 달라집니다.

따라서 구조 안정성·결함 형성에너지·흡착에너지는 total energy 차이로 계산해야 합니다. 특정 band eigenvalue의 합이나 DOS 면적을 총에너지 대신 쓰면 안 됩니다.

02고유값은 어디까지 물리적인가

정확한 Kohn-Sham DFT에서 최고 점유 고유값은 전자의 제거에너지와 연결되는 특별한 지위를 갖습니다. 그러나 나머지 고유값은 정확한 다전자 여기 스펙트럼을 일반적으로 보장하지 않습니다. 밴드구조가 유용한 이유와 모든 밴드 에너지가 엄밀한 관측값이라는 주장은 다릅니다.

기본 범함수는 band gap을 과소평가하는 경우가 많습니다. 이는 근사 범함수 오차와 전자 수가 정수를 지날 때 교환-상관 퍼텐셜이 뛰어야 하는 derivative discontinuity의 누락이 함께 관련됩니다.

출력별 안전한 첫 해석
  • 총에너지 차이: 같은 조성·방법·수치 기준에서 비교
  • 힘·응력: 에너지 미분으로 구조와 변형 방향 판단
  • 밀도·자화밀도: 공간 분포와 차분을 비교하되 분할법 의존성 명시
  • Kohn-Sham 밴드: 상태의 분산과 성격을 보는 모델이며 여기 에너지와의 차이 경계

03밀도에서 파생된 양도 정의가 필요하다

총 전자 밀도는 DFT의 기본 변수지만 원자별 전하·결합차수·국소 DOS는 공간이나 기저를 나누는 분석법에 의존합니다. Bader, Mulliken, Löwdin 같은 서로 다른 분할은 같은 밀도에서 다른 원자 전하를 줄 수 있습니다.

따라서 ‘원자에 0.63e가 이동했다’는 숫자는 사용한 분할법·격자·기저와 함께 보고하고, 절대 전하보다 동일한 방법 안의 경향을 비교하는 편이 안전합니다. DFT 출력 해석은 계산 단계가 아니라 모델링의 마지막 단계입니다.

04이 글의 용어 사전

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

이중계수double counting
같은 상호작용 기여를 서로 다른 항에서 두 번 포함하는 문제입니다.
Kohn-Sham 고유값
보조 독립전자 방정식의 고유값으로, 일반적으로 모든 실제 여기 에너지와 같지는 않습니다.
Fermi 준위Fermi level
점유를 정하는 전자 화학퍼텐셜로, 0 K 금속에서는 최고 점유 경계와 연결됩니다.
derivative discontinuity
전자 수가 정수를 지날 때 정확한 교환-상관 퍼텐셜에 나타나는 불연속적 이동입니다.
전하 분할charge partitioning
연속 전자 밀도를 원자 또는 영역별 기여로 나누는 분석 절차입니다.
상태밀도density of states
에너지 구간당 허용된 일전자 상태 수를 나타내는 함수입니다.

05참고자료와 더 읽을거리

  1. Kohn & Sham, Self-Consistent Equations Including Exchange and Correlation Effects (1965) — 상호작용 전자계를 계산 가능한 유효 독립전자 방정식으로 바꾼 원 논문입니다.
  2. Kohn, Nobel Lecture: Electronic structure of matter (1999) — 파동함수에서 밀도 범함수로 이어진 전자구조 이론의 맥락을 정리합니다.
  3. Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
NEXT DFT NOTE주기적 고체와 Bloch 정리 →