DFT 에너지와 온도·압력·화학퍼텐셜을 결합해 bulk·표면 상도를 만들기
진공에서 가장 낮은 표면에너지를 가진 구조가 산소 분위기에서도 안정하다고 볼 수 없습니다. 환경과 원자를 교환하면 조성이 다른 구조의 안정 순위가 화학퍼텐셜에 따라 바뀝니다.
0 K의 여러 DFT 총에너지를 실제 온도·분압 조건의 상 안정성 지도로 어떻게 바꿀까요?
이번 장의 답열린 계의 구조마다 Gibbs 자유에너지에서 reservoir 화학퍼텐셜 항을 빼 grand potential을 구성합니다. 기체 화학퍼텐셜의 T,p 의존성과 경쟁상 제약을 적용해 가장 낮은 구조의 영역을 그립니다.
01열린 계의 thermodynamic potential
원자 종류 j를 reservoir와 교환할 수 있는 계에서는 $G-\sum_jN_j\mu_j$를 최소화합니다. 표면은 면적으로 나누어 surface grand potential을 비교하고, bulk 화합물은 조성 제약과 경쟁상 조건을 사용합니다.
고체 $G_i$는 DFT 에너지에 진동·배치·전자 자유에너지를 더해 근사합니다. 많은 표면 계산에서 고체 진동 차이를 생략하지만 그 크기가 안정영역 경계보다 작은지 확인해야 합니다.
02기체 화학퍼텐셜의 T,p 의존성
O 원자 reservoir를 O$_2$ 기체로 잡으면 $\mu_O=\frac12\mu_{O_2}$. DFT O$_2$ 결합에너지 오차, 표준 열화학표의 thermal correction과 기준 일관성을 관리합니다.
허용 chemical potential은 원소 응축을 막는 상한과 host·경쟁상의 안정성을 만족하는 하한 사이에 있습니다. 이 범위를 벗어난 phase diagram 영역은 물리적으로 접근 불가능합니다.
03교차선과 안정영역 해석
각 구조의 $\Omega_i(\mu)$는 조성 차이에 따라 서로 다른 기울기의 선 또는 면이 됩니다. 가장 아래 envelope가 열역학적 안정 구조이며 교차점은 공존 경계입니다.
- chemical potential의 물리적 허용범위
- T,p 변환에 사용한 기체 모델과 표준상태
- 포함한 구조 후보와 누락 가능성
- 진동·배치 엔트로피 포함 여부
- DFT reference 보정과 범함수 민감도
04이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- ab initio thermodynamics
- DFT 에너지와 통계열역학·화학퍼텐셜을 결합해 환경별 안정성을 예측하는 접근입니다.
- grand potential
- 입자 reservoir와 교환하는 계에서 최소화되는 $\Omega=G-\sum N_i\mu_i$ 퍼텐셜입니다.
- 표준상태standard state
- chemical potential의 기준이 되는 정해진 압력·농도 상태입니다.
- 경쟁상competing phase
- 같은 원소 조합에서 host보다 안정해질 수 있어 chemical potential 범위를 제한하는 상입니다.
- 안정영역stability region
- 후보 중 특정 구조의 thermodynamic potential이 가장 낮은 환경 범위입니다.
- lower envelope
- 여러 에너지 선·면 가운데 각 조건에서 가장 낮은 부분을 이은 안정 경계입니다.
05참고자료와 더 읽을거리
- Reuter & Scheffler, Composition, structure, and stability of RuO2(110) (2001) — DFT 에너지와 화학퍼텐셜을 결합한 표면 ab initio thermodynamics의 대표 사례입니다.
- Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
- Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.