跳轉到

1930年 · Hartree-Fock 與 Slater:行列式結盟案#

卷宗編號:Case-1930-HF|偵探:Vladimir Fock、John Slater|前任:Hartree|伏筆:Roothaan

案發現場#

1927年 Hartree 自洽場破了重原子可算案,卻留下具無名屍:交換能。氦的單重—三重分裂算不準,分子鍵能系統偏差,審稿人頻頻退件。Pauli 原理像幽靈在走廊遊蕩——每軌道至多兩電子,可乘積波函數 \(\phi_1 \phi_2 \cdots\) 交換兩電子居然不變號,這在量子法條裡是非法的。

1929至1930年,兩位偵探從不同方向抵達同一現場:列寧格勒的 Fock 懂變分與算符,哈佛的 Slater 懂光譜與行列式。他們的共同委託是:給平均場補上反對稱,讓多電子波函數合法化,同時保留可算性。

案發現場的白板上寫著三行血字:波函數須反對稱、能量須變分最低、計算須可自洽。

偵查過程#

線索一:Slater 行列式——合法身份證#

Slater 1929年的靈感是:把軌道排成行列式,交換兩電子即交換兩行,行列式自動變號。 \(N\) 電子 Slater 行列式為:

\[ \Phi = \frac{1}{\sqrt{N!}} \det[\phi_1 \phi_2 \cdots \phi_N] \]

其中行標是電子座標 \(r_1,\dots,r_N\) ,列標是自旋軌道。上式展開即 \(N!\) 項的反對稱化和,每項都是軌道乘積,但整體嚴格滿足 Pauli 原理。兩電子同自旋同軌道時兩列相同,行列式為零——Pauli 禁戒變成線性代數。

波函數形式 是否反對稱 是否含交換
Hartree 乘積 否 否
Slater 行列式 是 是
精確 \(3N\) 維 \(\psi\) 是 全相關

一張表看清本案在階梯上的位置:比 Hartree 合法,比精確解便宜。

線索二:Fock 算符——交換現形#

Fock 1930年用變分法對行列式求極值,證明最優軌道滿足本徵方程。定義 Fock 算符為:

\[ \hat F = \hat h + \hat J - \hat K \]

其中 \(\hat h\) 為單電子核動能項, \(\hat J\) 為庫侖算符, \(\hat K\) 為交換算符。對軌道 \(\phi_i\) 有:

\[ \hat F \phi_i = \varepsilon_i \phi_i \]

庫侖項 \(\hat J\) 是經典電荷雲排斥,交換項 \(\hat K\) 則純量子、無經典對應:它只作用於同自旋電子,效果是挖一個「Fermi 洞」,讓同自旋電子自動迴避。這正是 Hartree 缺席的證人。

線索三:能量帳本與變分保證#

Hartree-Fock 總能量可寫成軌道期望加雙電子修正:

\[ E_{HF} = \sum_i h_i + \frac{1}{2}\sum_{ij} (J_{ij} - K_{ij}) \]

其中 \(h_i\) 為單電子積分, \(J_{ij}\) 為庫侖積分, \(K_{ij}\) 為交換積分,係數 \(1/2\) 照例去重。關鍵是變分原理保證:

\[ E_{HF} \ge E_{精確} \]

任何試探行列式的能量都是真基態的上界。偵探終於有法官:能量越低,波函數越好。自洽迴圈(猜軌道→建 \(\hat F\) →對角化→更新)於是有了明確的收斂 target。

線索四:相關能——新懸案的誕生#

Löwdin 後來把差值命名為相關能:

\[ E_{corr} = E_{精確} - E_{HF} < 0 \]

Hartree-Fock 抓到了交換,卻放走了「電子瞬時躲避」的動態相關。 He 原子 \(E_{HF}\) 誤差約 \(0.04\) Hartree(約 \(1\) 電子伏特),化學鍵離解、范德華力、過渡金屬多重態全是重災區。平均場讓每顆電子只見雲霧,看不見對方此刻的眼神——這筆欠款要等 CI、耦合簇、DFT、QMC 來還。

伏筆:Roothaan 在門外#

1930年方程是微分—積分方程,手算僅限原子。Roothaan(1951)把軌道展成基函數 \(\phi_i = \sum_\mu C_{\mu i} \chi_\mu\) ,化為矩陣廣義本徵問題 \(FC = SC\varepsilon\) ,分子計算才真正開張。本案結尾,Roothaan 還是研究生,但分子軌道法的子彈已上膛。

結案報告#

1930年行列式結盟案告破:Slater 給身份(反對稱),Fock 給武器(交換算符),Hartree 方法升級為 Hartree-Fock。自洽場從此合法,量子化學有了第一套可變分、可系統改進的流程。

遺產清單:

  • 正軌確立:單行列式+自洽+變分,後世半經驗、從頭算、DFT 軌道圖像全沿此軌。
  • 基準地位: \(E_{HF}\) 是 \(100\%\) 交換、 \(0\%\) 相關的參照點,相關能由此定義。
  • 計算遺產:四指標雙電子積分 \((\mu\nu|\lambda\sigma)\) 成為算力黑洞,直接催生高斯基組、直接 SCF、密度擬合。
  • 未竟線索:相關能、相對論、強關聯,移交戰後 CI、CC、MCSCF、DFT、QMC。

第一幕在此落幕:世界已寫成方程,方程已有近似解法,下一幕,電腦要把人力換成電子。

證據與工具#

  • 行列式: \(\Phi = \det[\phi_i]/\sqrt{N!}\) ,反對稱的最小合法結構。
  • Fock 算符: \(\hat F = \hat h + \hat J - \hat K\) ,自洽本徵方程 \(\hat F \phi_i = \varepsilon_i \phi_i\) 。
  • 能量公式: \(E_{HF} = \sum_i h_i + \frac{1}{2}\sum_{ij}(J_{ij} - K_{ij})\) ,滿足 \(E_{HF} \ge E_{精確}\) 。
  • 相關能定義: \(E_{corr} = E_{精確} - E_{HF}\) ,恆為負值。
  • 可重演實驗:用最小基算 He 原子 \(E_{HF}\) ,與 Hylleraas 精確值 \(-2.9037\) Hartree 對比,親眼見證相關能缺口。
  • 延伸閱讀:Fock 1930 Z. Physik、Slater 1929 Phys. Rev.;後續直通 Roothaan 方程與戰後 ENIAC 暴力計算。