유한 크기·기저·k-point·함수 선택의 오차

수치 오차·모델 오차·방법 오차·해석 오차를 분리해 결론의 취약점을 찾기

DFT · PART 5 · 숫자의 정확도와 이론의 적용 범위를 시험하는 동안06 / 08

계산과 실험이 0.2 eV 다를 때 원인은 범함수일 수도, 작은 셀·틀린 구조·기준 분자·온도 누락일 수도 있습니다. 모든 차이를 ‘DFT 오차’라고 부르면 개선 방향을 찾을 수 없습니다.

이번 장의 답수치·유한크기·모델·방법·환경·해석 오차로 분류하고 각 항을 독립적인 test로 추정합니다. 절대값보다 같은 근사에서의 차이가 상쇄될 수 있지만 구조·전자상태가 다르면 상쇄를 가정하지 않습니다.

ERROR TAXONOMY

01수치cutoff·k·SCF

02모델·방법cell·XC·spin

03해석·환경T·reference·observable

$\delta Q=\sum_i\frac{\partial Q}{\partial p_i}\delta p_i+\delta_{model}$

오차 항의 크기와 결론 민감도를 곱해 우선순위를 정해야 계산 자원을 효과적으로 배분할 수 있습니다.

01오차를 여섯 상자로 나누기

각 오차는 서로 결합합니다. 예를 들어 작은 결함 셀은 전하 국소화를 바꾸고, 그 변화가 범함수 민감도까지 바꿀 수 있습니다. 가장 큰 항을 줄인 뒤 다른 항을 다시 평가합니다.

같은 원소·결합·셀·범함수의 두 유사 구조를 빼면 core 기준과 큰 에너지 항이 상쇄되어 에너지 차이가 절대값보다 정확할 수 있습니다. 이것이 DFT가 상대 안정성에 강한 이유입니다.

그러나 분자와 금속 표면, 서로 다른 산화상태, 국소화 정도가 다른 상을 비교하면 범함수 편향이 다르게 작용합니다. 상쇄는 희망이 아니라 유사한 benchmark cycle과 방법 변화에서 확인해야 합니다.

실용적 error triage

모든 항을 같은 수준으로 줄일 필요는 없습니다. 1 eV 경향을 설명하는 연구와 10 meV polymorph 순위를 정하는 연구는 필요한 정확도와 방법이 다릅니다.

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

Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.
Kohn, Nobel Lecture: Electronic structure of matter (1999) — 파동함수에서 밀도 범함수로 이어진 전자구조 이론의 맥락을 정리합니다.