하트리-폭 방법

편집 역사 토론
양자화학 양자역학 계산화학 계산물리 마지막 수정: 2026-07-05 04:13:05
하트리-폭 방법
Hartree-Fock method
약칭HF
분야양자화학 × 전자구조 이론
핵심 아이디어평균장 근사 + 슬레이터 행렬식
푸는 방식자기무모순장(SCF) 반복
대표 프로그램Gaussian, GAMESS, PySCF, Psi4

1. 개요[편집]

하트리-폭 방법(Hartree-Fock method, HF)은 다전자 시스템의 슈뢰딩거 방정식을 근사적으로 풀기 위해, 전자 하나가 나머지 전자 전체가 만드는 평균장(mean field) 속에서 움직인다고 가정하는 전자구조 계산법이다. 즉 “전자들끼리 매 순간 서로를 밀쳐내며 춤추는” 진짜 상관관계를 포기하는 대신, “다른 전자들의 평균적인 전하 분포”만 느끼게 만들어 N체 문제를 N개의 1체 문제로 쪼개는 방법이다.1

다전자 슈뢰딩거 방정식은 전자가 3개만 넘어가도 해석적으로 풀리지 않는다. 전자-전자 반발항 1/rij1/r_{ij}이 모든 전자 쌍을 얽어놓기 때문이다. HF는 이 얽힘을 평균으로 뭉개버림으로써 계산 가능한 형태로 만든 최초의 실용적 방법이자, 오늘날 거의 모든 고급 양자화학 방법(밀도범함수이론의 하이브리드 범함수, MP2, CCSD 등)이 딛고 서는 출발점이다.

2. 다전자 문제와 평균장 근사[편집]

분자의 전자 해밀토니안(원자단위, 보른-오펜하이머 근사)은 대략 이렇게 생겼다.

H^=−∑i12∇i2−∑i,AZAriA+∑i<j1rij\hat{H} = -\sum_i \frac{1}{2}\nabla_i^2 - \sum_{i,A}\frac{Z_A}{r_{iA}} + \sum_{i<j}\frac{1}{r_{ij}}

앞의 두 항(운동에너지, 핵-전자 인력)은 전자별로 분리되지만, 마지막 전자-전자 반발항이 문제의 원흉이다. 이 항 때문에 파동함수를 전자별로 깔끔하게 분리할 수 없다.

하트리(D. R. Hartree)의 초기 발상은 이 반발을 각 전자가 느끼는 평균 퍼텐셜로 바꾸는 것이었다. 그런데 하트리의 원래 곱함수(Hartree product)는 전자가 페르미온이라는 사실, 즉 **반대칭성(antisymmetry)**을 무시했다. 두 전자를 맞바꾸면 파동함수의 부호가 뒤집혀야 하는데 단순 곱은 그러지 못한다. 여기서 폭(V. Fock)과 슬레이터(J. C. Slater)가 등판한다.

3. 슬레이터 행렬식[편집]

반대칭성을 강제하는 가장 우아한 장치가 **슬레이터 행렬식(Slater determinant)**이다. NN개의 스핀오비탈 χi\chi_i로 파동함수를 행렬식 형태로 쓴다.

Ψ=1N!∣χ1(x1)χ2(x1)⋯χN(x1)χ1(x2)χ2(x2)⋯χN(x2)⋮⋮⋱⋮χ1(xN)χ2(xN)⋯χN(xN)∣\Psi = \frac{1}{\sqrt{N!}} \begin{vmatrix} \chi_1(\mathbf{x}_1) & \chi_2(\mathbf{x}_1) & \cdots & \chi_N(\mathbf{x}_1) \\ \chi_1(\mathbf{x}_2) & \chi_2(\mathbf{x}_2) & \cdots & \chi_N(\mathbf{x}_2) \\ \vdots & \vdots & \ddots & \vdots \\ \chi_1(\mathbf{x}_N) & \chi_2(\mathbf{x}_N) & \cdots & \chi_N(\mathbf{x}_N) \end{vmatrix}

행렬식은 두 행(두 전자)을 바꾸면 부호가 반전되므로 반대칭성이 공짜로 딸려온다. 게다가 두 전자가 같은 스핀오비탈을 차지하면 두 열이 같아져 행렬식이 0이 되는데, 이게 바로 파울리 배타 원리다. 수학이 물리를 이렇게 얌전히 챙겨주는 경우도 드물다.2

4. 자기무모순장(SCF)[편집]

HF의 진짜 골때리는 지점은 닭이 먼저냐 달걀이 먼저냐 문제다. 각 전자가 느끼는 평균장을 계산하려면 다른 전자들의 오비탈을 알아야 하는데, 그 오비탈은 또 각자의 평균장을 풀어야 나온다. 서로가 서로를 정의한다.

그래서 반복법으로 돌린다. 이게 자기무모순장(Self-Consistent Field, SCF) 절차다.

  1. 초기 오비탈(추정치)로 전자 밀도를 만든다.
  2. 그 밀도로 포크 연산자 F^\hat{F}를 구성한다.
  3. 포크 방정식 F^χi=εiχi\hat{F}\chi_i = \varepsilon_i \chi_i을 풀어 새 오비탈을 얻는다.
  4. 새 오비탈로 밀도를 다시 만든다.
  5. 밀도(또는 에너지)가 이전 반복과 거의 안 변할 때까지 2~4를 반복한다.

실무에서 기저함수(basis set)를 도입해 포크 방정식을 행렬 형태로 바꾼 것이 그 유명한 루탄 방정식(Roothaan equation) FC=SCE\mathbf{FC} = \mathbf{SCE}이다. 이 일반화 고윳값 문제를 뉴턴-랩슨법 계열의 가속 기법(DIIS 등)과 함께 반복해서 푼다. SCF가 수렴 안 하면? 초기 추정을 바꾸거나, 감쇠(damping)를 걸거나, 결국 신에게 맡긴다.3

5. 교환에너지와 상관 누락[편집]

HF의 포크 연산자에는 고전적인 쿨롱 반발(Coulomb, J^\hat{J}) 외에 순수 양자역학적 항인 교환(exchange, K^\hat{K}) 항이 등장한다. 교환에너지는 같은 스핀 전자들이 반대칭성 때문에 서로를 자연스럽게 회피하면서 생기는 에너지 이득으로, 고전물리에는 대응물이 없는 순수 양자 효과다. HF는 이 교환을 정확히(exactly) 다룬다는 점이 큰 장점이다.4

그러나 여기 치명적 약점이 있다. HF는 반대칭성에서 오는 같은 스핀 전자의 회피는 챙기지만, 스핀이 다른 전자들이 순간순간 서로를 밀쳐내며 피하는 **동적 상관(dynamic correlation)**은 평균장으로 뭉개버린다. 이 놓친 부분을 **상관에너지(correlation energy)**라 부르며, 정의상

Ecorr=Eexact−EHFE_\text{corr} = E_\text{exact} - E_\text{HF}

이다. 크기는 전체 에너지의 1% 남짓으로 작아 보이지만, 화학 반응의 에너지 차이(결합 에너지, 반응 장벽)는 딱 이 스케일에서 결정되기 때문에 상관에너지를 놓치면 정량적 화학은 물 건너간다. 그래서 HF 위에 상관을 얹는 포스트-HF 방법(MP2, CI, Coupled Cluster)이 발달했고, 다른 노선인 밀도범함수이론은 상관을 교환-상관 범함수로 근사해 훨씬 싼 값에 챙긴다.

6. 장단점 요약[편집]

항목HF의 태도
반대칭성/교환정확히 처리 (자랑거리)
전자 상관동적 상관 통째로 누락
계산 비용형식상 기저함수 4제곱 스케일
위치거의 모든 고급 방법의 출발점

7. 여담[편집]

  • HF로 계산한 원자화 에너지, 반응열은 실험과 수십 kcal/mol씩 어긋나는 일이 흔하다. “정성적으로는 맞는데 정량적으로는 틀린” 방법의 대명사.
  • 그럼에도 HF를 배우지 않고 양자화학을 논할 수는 없다. 모든 교재의 1장이 HF인 데는 이유가 있다.
  • 현대 계산화학의 주력은 밀도범함수이론이지만, 그 하이브리드 범함수(B3LYP 등)는 HF 교환을 일정 비율 섞어 쓴다. HF는 죽지 않고 부품으로 살아남았다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 그래서 HF를 “독립 입자 근사(independent particle approximation)“라고도 부른다. 전자들이 서로 눈치는 보되(평균장), 실시간으로 손잡고 피하지는 않는(상관 없음) 세계관이다. ↩

  2. 슬레이터 행렬식 한 개만으로는 사실 여전히 근사다. 정확한 파동함수는 무한히 많은 행렬식의 선형결합(Full CI)으로 표현되며, HF는 그중 딱 하나의 최적 행렬식이다. ↩

  3. SCF 발산은 특히 전이금속 착물, 라디칼, 늘어난 결합 같은 “전자가 헷갈리는” 시스템에서 잘 터진다. DIIS(Direct Inversion of the Iterative Subspace)를 켜는 것이 국룰이지만, 그마저도 안 되면 초기 게스를 손으로 만져야 한다. ↩

  4. 정확히는 “HF 근사의 틀 안에서 정확히”다. 이 교환을 근사(local exchange)로 바꾼 것이 초기 DFT의 슬레이터 교환이고, HF 교환을 그대로 가져다 섞은 것이 하이브리드 DFT다. ↩