강도든 반응성이든, 재료의 성질은 결국 전자가 어떻게 배치되는가에서 비롯됩니다. 그 전자를 양자역학 으로 직접 계산하는 것이 원자 스케일 해석의 출발점입니다.

PART 5 · 전자 한 개에서 구조물까지 계산을 이어 가는 동안02 / 05

서로 다른 원자 배열 두 개 중 어느 구조가 안정한지, 수소가 표면에서 어디에 붙는지, 반응 중 원자에 어떤 힘이 작용하는지 알고 싶습니다. 실험값으로 맞춘 힘장 대신 원자번호와 좌표에서 출발합니다.

수많은 전자의 양자 문제를 실제 재료 계산이 가능한 형태로 어떻게 바꿀까요?

이번 장의 답Born–Oppenheimer 근사로 핵과 전자를 나누고, 다체 파동함수 대신 전자 밀도를 쓰는 Kohn–Sham DFT를 자기일관적으로 풀어 에너지·힘을 얻습니다.

01제일원리와 다체 문제

원리적으로 전자의 상태는 슈뢰딩거 방정식을 풀면 얻습니다.

$$ \hat{H}\,\Psi = E\,\Psi $$

문제는 $\Psi$ 가 모든 전자의 좌표에 동시에 의존하는 다체 파동함수라는 점입니다. 전자가 $N$ 개면 변수가 $3N$ 개로 폭발해, 몇 개만 넘어도 직접 풀 수 없습니다. 실험값에 기대지 않는 이런 계산을 제일원리(ab initio)라 합니다.

먼저 전자보다 훨씬 무거운 원자핵은 전자 운동 중 거의 정지해 있다고 보는 Born–Oppenheimer 근사를 사용합니다. 핵 좌표를 고정한 채 전자 바닥상태 에너지를 구하면, 그 에너지가 핵이 움직이는 포텐셜 에너지면이 됩니다. “제일원리”도 근사가 없다는 뜻이 아니라, 경험적으로 맞춘 재료별 파라미터 없이 기본 방정식과 통제 가능한 근사에서 출발한다는 뜻입니다.

02Hohenberg–Kohn에서 Kohn–Sham으로

밀도범함수이론(DFT)의 핵심 발상은, 다루기 힘든 다체 파동함수 대신 전자 밀도$n(\mathbf{r})$(3차원 함수 하나)로 에너지를 표현할 수 있다는 것입니다.

$$ E = E[\,n(\mathbf{r})\,] $$

Hohenberg–Kohn 정리는 바닥상태 에너지가 밀도의 범함수이고, 올바른 밀도에서 최소가 된다는 토대를 줍니다. 실제 계산은 상호작용하는 전자를 같은 밀도를 만드는 유효한 독립 전자 문제로 바꾼Kohn–Sham 방정식을 풉니다.

$$ \left[-\frac{\hbar^2}{2m_e}\nabla^2+V_{\mathrm{eff}}[n](\mathbf r)\right]\psi_i=\varepsilon_i\psi_i, \qquad n(\mathbf r)=\sum_i f_i|\psi_i(\mathbf r)|^2 $$

$V_{\mathrm{eff}}$에는 핵의 외부 퍼텐셜, 전자 밀도의 고전적 Coulomb 항과 교환–상관 퍼텐셜이 들어갑니다. 마지막 항의 정확한 형태를 모르기 때문에 LDA, GGA, hybrid 같은 근사를 선택하며, 이 선택이 결합에너지·격자상수·밴드갭의 체계적 오차를 좌우합니다.

03자기일관성, 기저와 수렴검사

입력 밀도로 유효 퍼텐셜을 만들고 Kohn–Sham 방정식을 풀어 새 밀도를 얻은 뒤, 입력과 출력이 같아질 때까지 섞어 반복합니다. 이것이 SCF(self-consistent field) 과정입니다. 주기적 고체에서는 파동함수를 평면파로 전개하고, 핵심 전자는 pseudopotential 또는 PAW로 묶어 계산비를 줄이는 경우가 많습니다.

CONVERGENCE BEFORE INTERPRETATION
  • 평면파 cutoff와 k-point 격자를 늘려 에너지·힘이 수렴하는지 확인
  • 초격자 크기로 결함·표면의 주기적 이미지 상호작용을 점검
  • 교환–상관 범함수, 분산 보정과 스핀 상태가 질문에 맞는지 비교
  • 구조 최적화의 힘·응력 기준과 전자 SCF 기준을 별도로 기록

04AIMD — 전자 바닥상태 위에서 움직이는 원자

정적 DFT 구조를 시간에 따라 이어 붙인 것이 AIMD입니다. Born–Oppenheimer AIMD에서는 매 순간 전자 SCF를 충분히 수렴시켜 원자에 작용하는 힘 $\mathbf{F}_I = -\nabla_{\mathbf{R}_I} E$ 를 구하고, 그 힘으로 원자를 조금씩 움직이는 것을 반복합니다. 그러면 온도가 있는 상태에서 원자의 진동·이동·반응을 원자 수준에서 직접볼 수 있습니다. 전자구조를 매 스텝 풀기 때문에 보통 수십~수백 원자와 ps 규모에 머물며, 희귀 사건이나 느린 확산을 직접 관찰하지 못할 수 있습니다.

DFT vs AIMD
  • $\hat H\Psi=E\Psi$ → 다체 문제(직접 풀기 불가)
  • DFT: $E=E[n(\mathbf r)]$ — 밀도로 실용화(코온–샴)
  • AIMD: $\mathbf F_I=-\nabla_{\mathbf R_I}E$ 로 원자를 이동

05무엇을 믿고 무엇을 경계할까

DFT/AIMD 결과는 실험적 포텐셜이나 MLFF의 중요한 기준 데이터가 됩니다. 그러나 교환–상관 근사, 유한 크기, 제한된 샘플링과 전자 바닥상태 가정 때문에 절대적인 “정답”은 아닙니다. 밴드갭, 약한 분산력, 강상관 전자계는 특히 방법 선택에 민감합니다. 크기·시간 한계로 미세구조를 직접 다루기 어렵기 때문에 다음 단계에서 고전 MD와 MLFF로 확장합니다.

06이 글의 용어 사전

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

Born–Oppenheimer 근사Born–Oppenheimer approximation
무거운 원자핵의 운동과 빠른 전자 운동을 분리해 고정된 핵 좌표에서 전자 문제를 푸는 근사입니다.
전자 밀도electron density
위치 r 주변에 전자가 존재하는 분포 n(r)입니다. DFT의 기본 변수로 세 공간 좌표에만 의존합니다.
범함수functional
숫자가 아니라 함수 전체를 입력으로 받아 값을 내는 함수입니다. DFT 에너지는 밀도의 범함수입니다.
Kohn–Sham 오비탈Kohn–Sham orbital
실제 다체 파동함수는 아니지만 정확한 전자 밀도를 재현하도록 도입한 유효 독립전자 함수입니다.
교환–상관 범함수exchange-correlation functional
전자 교환과 상관 효과 중 알려진 항에 포함되지 않은 부분을 근사하는 DFT의 핵심 선택입니다.
SCFself-consistent field
입력 밀도에서 만든 퍼텐셜로 새 밀도를 계산하고 두 밀도가 일치할 때까지 반복하는 절차입니다.
k-pointk-point sampling
주기적 결정의 역공간 Brillouin zone 적분을 수치적으로 표본화하는 점입니다.
AIMDab initio molecular dynamics
각 시간 스텝의 힘을 전자구조 계산으로 얻어 원자핵의 운동방정식을 적분하는 방법입니다.

07참고자료와 더 읽을거리

  1. Sholl & Steckel, Density Functional Theory: A Practical Introduction — Ch. 1–3의 DFT 기초와 수렴, Ch. 9의 AIMD, Ch. 10의 정확도와 한계.
  2. Hohenberg & Kohn, Physical Review 136, B864 (1964) — 전자 밀도에 관한 두 기본 정리.
  3. Kohn & Sham, Physical Review 140, A1133 (1965) — 실용적인 독립전자 방정식의 도입.
NEXT NOTE분자동역학(MD) — 원자를 뉴턴 법칙으로 움직이기 →