ハートリー＝フォック方程式
ハートリー=フォック方程式(ハートリーフォックほうていしき、Hartree–Fock equation)は、多電子系を表すハミルトニアンの固有関数(波動関数)を一個のスレーター行列式で近似(ハートリー=フォック近似)した場合に、それが基底状態に対する最良の近似となるような(スピンを含む)1電子分子軌道の組を探し出すための方程式である。ウラジミール・フォックによって導かれた。分子軌道法の基本となる方程式である。
ハートリー=フォック方程式
は、\{ \varphi_i \}の近似的な解が与えられた場合、方程式中の\{ \varphi_i \}置換することで方程式
が誘導される。すなわちこの方程式の\hat{F}には固有関数\psiは含まれず、普通の固有値方程式として解くことが出来る。
これにより得られた解を近似解として適用し再帰的に解く事で、多電子系のフェルミ粒子(この場合は電子)全体の作る平均場と、その中で一粒子運動をするフェルミ粒子の波動関数を自己無撞着に決定することができる(SCF法)。
ハートリー=フォック法では大きく4つの近似をする。
ボルン–オッペンハイマー近似を適用する。すなわち、本来ならば分子の原子核と電子それぞれの座標についての関数である分子全体の波動関数を、原子核座標は不変とし電子のみの座標の関数とみなす。
相対論の効果は無視し、運動量演算子は非相対論的なものと仮定する。
それぞれのエネルギー固有関数(定常状態のシュレーディンガー方程式の解)は1つのスレイター行列式で記述できると仮定する。
平均場近似を適用する。つまりある1つの電子が受ける相互作用の大きさは、その電子の位置のみに依存し、他の電子の位置には依存しない。この仮定から外れることによる効果、つまり電子相関は、反平行スピンどうしの電子では無視されるが、平行スピンどうしの電子では考慮される。 (電子相関と「電子の交換」を混同しないように。電子の交換相互作用はハートリーフォック法で考慮されている。)
最後の2つの近似を仮定しないことによって、多くのポスト-ハートリー-フォック法が作られている。
N個のフェルミオン系を考える。分子全体の波動関数を1つのスレーター行列式とし、時間依存しないシュレーディンガー方程式に代入すると、エネルギーは次のように書ける。
E= \sum_{i=1}^{N} \int \mathrm dx \varphi_i^{*}(x) \left(-\frac{1}{2m} \nabla^{2}\right) \varphi_i(x)+ \frac{1}{2} \sum_{i,k=1}^{N}\int \mathrm dx \mathrm dy \varphi_i^{*}(x)\varphi_k^{*}(y)v(x-y)\left(\varphi_i(x)\varphi_k(y)-\varphi_i(y)\varphi_k(x)\right)
これを一粒子波動関数\varphi_i^{*}(x)で変分する。つまり\varphi_i^{*}(x)が規格直交化されており、かつ最低のエネルギーをとるようものを、ラグランジュの未定乗数法などで探すことで、以下のハートリー=フォック方程式を得る。
-\frac{1}{2m} \nabla^{2} \varphi_i(x)+V_H(x)\varphi_i(x)-\int \mathrm dy V_E(x,y)\varphi_i(y) = \epsilon_i \varphi_i(x)
ここで、
\varphi_i (x) \  :一粒子波動関数
V_H(x)\equiv \int \mathrm dy \sum_{k=1}^{N} v(x-y) |\varphi_k(y)|^2 :ハートリーポテンシャル
V_E(x,y) \equiv \sum_{k=1}^{N} v(x-y) \varphi_k(x)\varphi_k^{*}(y) :フォックポテンシャル
フォックポテンシャルは、波動関数の反対称化が必要なフェルミオン多体系に特有のものであり、ボソン多体系の平均場を求める方程式(グロス・ピタエフスキー方程式と呼ばれている)には存在しない。
ハートリー=フォック方程式の解をユニタリ変換したものも、ハートリー=フォック方程式の解になっている。よってユニタリ変換をどのように選ぶかによって、いろいろな解の表現の仕方がある。そこで、ユニタリ変換後のハートリー-フォック方程式の未定定数が対角形になるようなユニタリ変換を選んで表したものを正準ハートリー=フォック方程式 (canonical Hartree-Fock equation) と呼ぶ。
正準ハートリー=フォック方程式は、フォック演算子の固有値方程式である。つまり固有値としてスピン軌道エネルギー εi、それに属する固有関数としてスピン軌道 \varphi_iをもつ固有値方程式である。
\epsilon_i \varphi_i
\hat{F} \equiv \hat{h} + \sum_{j=1}^n (\hat{J}_j - \hat{K}_j)
核–一電子ハミルトニアン
\hat{h} \equiv -\frac{1}{2}\nabla_i^2 - \sum_{A=1}^N \frac{Z_A}{r_{iA}
註)第一項は i 番目の電子の運動エネルギー、第二項は原子核-電子間の引力のポテンシャルエネルギーを表す。
\hat{J}_j(1) \varphi_i(1) \equiv \int \varphi_j^{*}(2) \varphi_j(2) \frac{1}{r_{12} \varphi_i(1) \,\mathrm d \boldsymbol{x}_2
註)位置 x2 にある一個の電子が χj で表される一個の電子から感じる平均的なポテンシャルを表す。
\hat{K}_j(1) \varphi_i(1) \equiv \int \varphi_j^{*}(2) \varphi_i(2) \frac{1}{r_{12} \varphi_j(1) \,\mathrm d\boldsymbol{x}_2
交換演算子は古典的解釈のできない演算子であり、単にスレーター行列式のアーティファクトである。つまり、区別が付かない電子に番号を付けたため、パウリの原理が要請する反対称性を満たす波動関数としてスレーター行列式を導入したことが原因で生じている。なお、クーロン項(ハートリーポテンシャル)中に存在する電子とそれ自身との「自己相互作用」は交換(フォック)ポテンシャル中にも存在するため打ち消される。
ここで、x は電子の空間座標 r とスピン座標 ω をまとめた空間スピン座標、
\varphi_j(\boldsymbol{x}_1)=\varphi_j(\boldsymbol{r}_1,\omega_1)
である。
ハートリー=フォック方程式はこのままの形では解くことが難しい。そこで通常は求めるスピン軌道を既知の基底関数の組で展開し行列方程式の形へ変換して解く。
いずれにしろ、フォック演算子のうちクーロン演算子と交換演算子が求めようとしているスピン軌道を含むため、つじつまの合った場の方法(自己無撞着場の方法あるいはSCF法とも呼ばれる)によって解く。
電子の出入りによって分子軌道が変化しないと仮定する。ハートリー=フォックエネルギーE[N]から、k番目の電子が抜き取られた後のN−1電子系のエネルギーE[N-1]を引くと\varepsilon_kとなり、-\varepsilon_kは電子を抜き取るために必要なエネルギー、つまりイオン化エネルギーの意味を持つ(クープマンズの定理)。したがって、ハートリー=フォック方程式の未定乗数εiは分子軌道エネルギーと解釈することができる。しかし、ハートリー=フォックエネルギーはεiの総和ではないことに留意すべきである。