원자번호와 좌표에서 에너지·힘·전자 밀도를 얻어 재료의 질문에 답하는 첫 지도

DFT · PART 1 · 안정한 원자 배열을 찾기 위해 전자 문제의 문 앞에 서는 동안01 / 05

FCC와 BCC 두 구조의 금속이 있습니다. 실험을 하기 전 어느 구조가 더 안정한지, 표면의 어느 자리에 분자가 붙는지 알고 싶습니다. 출발점으로 주어진 것은 원자의 종류와 좌표뿐입니다.

DFT 계산은 어떤 입력을 받아 무엇을 예측하며, 그 숫자는 재료의 질문과 어떻게 연결될까요?

이번 장의 답DFT는 고정된 핵 배치에서 전자 바닥상태를 구해 총에너지·전자 밀도·힘을 제공합니다. 서로 일관된 모델의 에너지 차이와 에너지의 미분을 해석하면 안정성·구조·반응을 비교할 수 있습니다.

01재료의 질문을 전자 문제로 바꾸기

재료 연구에서 DFT가 강한 이유는 실험을 대체해서가 아니라, 실험에서 분리하기 어려운 원자 수준의 원인을 통제된 모델로 비교할 수 있기 때문입니다. 결정 구조 A와 B의 에너지 차이, 결함이 없는 셀과 공공이 있는 셀의 차이, 깨끗한 표면과 흡착된 표면의 차이를 계산하면 관찰된 현상의 가능한 원인을 좁힐 수 있습니다.

질문은 반드시 계산 가능한 비교로 번역해야 합니다. “이 촉매가 좋은가?”는 너무 큽니다. “같은 범함수와 수치 조건에서 N$_2$ 해리 장벽이 어느 표면 자리에서 낮은가?”처럼 모델·비교 대상·관심량을 명시하면 DFT가 답할 수 있는 문제가 됩니다.

DFT가 자주 답하는 네 가지 질문
  • 어느 원자 배열의 총에너지가 더 낮은가
  • 원자를 조금 움직였을 때 힘과 응력은 어느 방향인가
  • 전자 밀도와 상태밀도에서 결합·전도·자성의 단서는 무엇인가
  • 초기상태에서 최종상태로 가는 에너지 장벽은 얼마인가

02핵심 출력 세 가지: 에너지·밀도·힘

첫 번째 출력은 총에너지 $E_0$입니다. 절대값 하나보다 같은 기준으로 계산한 차이가 중요합니다. 원자 수와 조성이 같으면 구조 간 에너지 차이를 비교할 수 있고, 조성이 다르면 화학퍼텐셜을 포함한 열역학적 순환을 세워야 합니다.

$$ \Delta E = E_{\mathrm{model\,B}}-E_{\mathrm{model\,A}} $$
$\Delta E<0$이면 선택한 계산 조건 안에서 B가 A보다 낮은 전자 바닥상태 에너지를 갖습니다.

두 번째는 전자 밀도 $n(\mathbf r)$입니다. 전하가 공간에 어떻게 분포하는지 보여 주며, 차분 밀도·국소 상태밀도·자화 밀도 같은 파생량으로 결합과 전하 재배치를 해석합니다. 세 번째는 힘 $\mathbf F_I=-\partial E/\partial\mathbf R_I$입니다. 힘이 0에 가까워질 때까지 원자를 움직이면 주어진 모델 안의 안정 구조를 찾을 수 있습니다.

출력직접 의미재료 질문으로 번역
총에너지주어진 핵 배치의 전자 바닥상태 에너지상대 안정성·결합에너지·반응에너지
전자 밀도공간별 전자 분포결합·전하이동·전자구조
힘·응력에너지의 좌표·변형률 미분구조 최적화·격자·탄성·동역학

03첫날부터 알아야 할 경계

보통의 Kohn-Sham DFT는 바닥상태 이론입니다. 유한온도의 엔트로피, 긴 시간의 확산, 실제 시편의 결함 분포가 자동으로 포함되지 않습니다. 필요한 물리를 자유에너지, phonon, AIMD, 통계역학 또는 상위 스케일 모델로 덧붙여야 합니다.

또한 결과는 범함수·기저·k-point·초격자·스핀 상태에 의존합니다. 따라서 좋은 DFT 연구는 계산을 많이 한 연구가 아니라, 질문에 맞는 모델을 만들고 오차가 결론보다 작은지 보인 연구입니다. 이 과목 전체는 바로 그 판단 능력을 만드는 과정입니다.

04이 글의 용어 사전

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

제일원리first principles / ab initio
재료별 실험값으로 맞춘 힘장 대신 원자번호와 기본 양자역학에서 출발하는 계산 접근입니다.
바닥상태ground state
주어진 외부 퍼텐셜과 입자 수에서 허용되는 가장 낮은 에너지 상태입니다.
총에너지total energy
전자 운동·전자-핵·전자-전자·핵-핵 상호작용을 일관되게 포함한 비교용 에너지입니다.
전자 밀도electron density
위치 r 주변의 전자 분포를 나타내는 세 공간 좌표의 함수입니다.
atomic force
총에너지를 원자 좌표로 미분한 음의 기울기로, 구조가 이동하려는 방향을 나타냅니다.
계산 모델computational model
원자 종류·좌표·격자·경계조건·전하·스핀을 포함해 실제 계를 제한적으로 표현한 대상입니다.

05참고자료와 더 읽을거리

  1. Sholl & Steckel, Density Functional Theory: A Practical Introduction평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
  2. Kohn, Nobel Lecture: Electronic structure of matter (1999)파동함수에서 밀도 범함수로 이어진 전자구조 이론의 맥락을 정리합니다.
  3. Hohenberg & Kohn, Inhomogeneous Electron Gas (1964)바닥상태 전자 밀도에 관한 두 정리로 DFT의 존재 근거를 세운 원 논문입니다.
NEXT DFT NOTE파동함수·연산자·고유값: 양자역학의 최소 언어