격자상수·상대 안정성·상태방정식

에너지-부피 곡선에서 평형부피·bulk modulus·압력에 따른 상 경쟁 읽기

DFT · PART 4 · 계산 출력에서 재료의 구조와 물성을 읽어 내는 동안03 / 08

구조 최적화 한 번으로 격자상수는 얻을 수 있지만, 부피 주변의 에너지 곡선을 만들면 평형점의 안정성과 압축에 대한 저항까지 함께 읽을 수 있습니다.

이번 장의 답평형부피 주변 여러 부피에서 내부좌표와 필요 시 셀 모양을 완화해 $E(V)$를 만들고 EOS에 적합합니다. 압력에서는 에너지 대신 $H=E+PV$를 비교하며, 수렴·진동·온도 효과를 분리합니다.

ENERGY CURVATURE AS COMPRESSIBILITY

01부피 scan$E(V_i)$

02EOS fitting$V_0,B_0,B'_0$

03압력 안정성$H=E+PV$

$B_0=V_0\left.\frac{d^2E}{dV^2}\right|_{V_0}$

에너지 최소 위치는 평형부피, 곡률은 압축 저항, 접선 기울기는 압력을 나타냅니다.

01평형점 주변의 E-V 데이터를 만들기

예상 평형부피의 대략 ±5~10% 범위에서 여러 부피를 선택하고 각 부피에서 내부 원자좌표를 충분히 완화합니다. 비입방 구조는 부피만 등방 축척할지 c/a와 셀 모양도 완화할지 질문에 따라 결정합니다.

각 점은 같은 cutoff·k-point 밀도·smearing·자성 상태를 사용해야 합니다. 부피에 따라 자성 또는 전자 상태가 바뀌면 하나의 매끄러운 EOS로 억지로 맞추지 말고 상태별 branch를 분리합니다.

Birch-Murnaghan 같은 상태방정식은 $E(V)$ 데이터를 평형부피 $V_0$, bulk modulus $B_0$, 압력 미분 $B'_0$로 요약합니다. $B_0$가 클수록 부피 압축에 대한 에너지 곡률이 큽니다.

격자상수를 실험과 비교할 때는 실험 온도, 영점진동과 열팽창, 시편 조성 및 범함수의 체계적 편향을 고려합니다. 0 K 정적 DFT와 실온 XRD의 차이를 모두 수치 오차로 부르면 안 됩니다.

외부 압력 P에서 안정상은 $H=E+PV$가 낮은 상입니다. 두 상의 enthalpy 곡선이 만나는 압력이 정적 상전이 압력 후보가 됩니다. 온도까지 포함하려면 Gibbs 자유에너지 $G=F+PV$를 비교합니다.

$$ \Delta H(P)=H_B(P)-H_A(P)=\Delta E(P)+P\Delta V(P) $$

고압에서는 작은 부피의 상이 $P\Delta V$ 항으로 유리해질 수 있습니다.

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

Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.
Kohn & Sham, Self-Consistent Equations Including Exchange and Correlation Effects (1965) — 상호작용 전자계를 계산 가능한 유효 독립전자 방정식으로 바꾼 원 논문입니다.