跳轉到

1927年 · Hartree 自洽場：平均場臥底案#

卷宗編號：Case-1927-SCF｜偵探：Douglas Hartree｜線人：Bohr、Slater（遠端）｜嫌犯： \(3N\) 維波函數

案發現場#

薛丁格1926年立下大功，也留下大坑：多電子波函數 \(\psi(r_1,\dots,r_N)\) 住在 \(3N\) 維，氦算得動，鈉以上全軍覆沒。光譜學家手握大把原子數據，理論卻交不出能量——案子在1927年堆到 Hartree 桌上。

Hartree 是劍橋的數值獵人，早年玩過火炮彈道差分機，深知硬算 \(3N\) 維等於自殺。案發現場的證物很單純：

  1. 哈密頓量裡的電子排斥項 \(1/r_{ij}\) 讓變數糾纏，方程不可分離。
  2. Bohr 殼層模型暗示電子大體各行其道，彼此只感到平均雲霧。
  3. 實驗游離能與 X 射線譜暗示每個電子彷彿活在某個有效中心勢裡。

委託很明確：把不可解的 \(N\) 體糾纏，降級成 \(N\) 個可解的一體問題，還不能太離譜。

偵查過程#

線索一：臥底計畫——人人只見平均場#

Hartree 的臥底方案是：假裝每個電子看不見其他電子的瞬時位置，只看見它們的平均電荷雲。設第 \(i\) 個電子軌道為 \(\phi_i(r)\) ，則它感受到的有效勢為核吸引加電子平均排斥：

\[ V^{SCF}_i(r) = -\frac{Z}{r} + \sum_{j \ne i} \int \frac{|\phi_j(r')|^2}{|r - r'|} dr' \]

其中 \(Z\) 為核電荷， \(r'\) 為其他電子的座標變數。上式 правой 部分全由軌道密度 \(|\phi_j|^2\) 構成——場由軌道定，軌道又由場定，迴圈於焉形成。

單電子方程於是降維：

\[ \left[ -\frac{1}{2}\Delta + V^{SCF}_i(r) \right] \phi_i = \varepsilon_i \phi_i \]

原子單位下動能為 \(-\Delta/2\) ， \(\varepsilon_i\) 為軌道能量。 \(3N\) 維慘案化為 \(N\) 個三維案件，算力需求從指數掉成線性。

線索二：自洽迴圈——跟自己對口供#

既然場與軌道互為因果，Hartree 發明物理史上第一個迭代演算法：

步驟 動作 收斂判據
\(1\) 猜初始軌道 \(\phi_j\) 經驗屏蔽常數
\(2\) 構造 \(V^{SCF}_i\) Poisson 積分
\(3\) 解徑向方程得新 \(\phi_i\) 打靶法＋節點計數
\(4\) 比較新舊場 差 \(< 10^{-4}\) 即自洽
\(5\) 未收斂則混合後回步驟 \(2\) 阻尼避免振盪

這就是「自洽場」（Self-Consistent Field）。電子們像一屋子互相打聽的嫌疑犯，每輪口供更新，直到沒人再改口。Hartree 用機械微分分析儀（analog computer） millimeters 級地手搖積分，算出 Na 到 Hg 的原子結構。

線索三：總能量不能連加#

新手警員常犯的錯是把軌道能量 \(\varepsilon_i\) 直接加總。Hartree 糾正：排斥能被算了兩次，須扣回。總能量為：

\[ E = \sum_i \varepsilon_i - \frac{1}{2}\sum_{i \ne j} J_{ij} \]

其中 \(J_{ij}\) 為庫侖積分，定義為：

\[ J_{ij} = \iint \frac{|\phi_i(r)|^2 |\phi_j(r')|^2}{|r - r'|} dr dr' \]

係數 \(1/2\) 正是去重複計算的修正。這條會計規則後世 Hartree-Fock、DFT 全盤繼承。

線索四：缺席的證人——交換#

本案有個致命破綻：Hartree 把多體波函數寫成軌道乘積 \(\Phi = \phi_1 \phi_2 \cdots \phi_N\) ，既不反對稱，也沒有交換項。Pauli 原理只靠「每軌道至多兩電子」手工維持，交換能 \(K_{ij}\) 全缺。

後果是算出的激發態、三重態分裂系統性跑偏。Fock 看出問題，Slater 備好行列式，1930年就要來補破網。但1927年當下，平均場臥底已是巨大勝利：重原子第一次能算了。

結案報告#

Hartree 自洽場案形式上告破：用 \(V^{SCF}\) 加迭代，把量子多體變成可算流程。鈉、汞、X射線能階的計算與實驗大致相符，Bohr 模型從定性升級為定量。

遺產清單：

  • 演算法遺產：自洽迭代（猜場→解軌道→更新場）是 SCF、Kohn-Sham、DMFT 的共同祖先。
  • 物理遺產：平均場、屏蔽、有效核電荷概念，化學家至今掛在嘴邊。
  • 計算遺產：證明原子結構是「可算」的，吸引 Slater、Fock、Roothaan 入場。
  • 未竟線索：無交換、無相關、非變分嚴格（乘積函數不滿足反對稱），三宗罪移交1930年 Hartree-Fock 案。

平均場是向維度災難交的第一筆贖金，不是全款。

證據與工具#

  • 有效勢： \(V^{SCF}_i\) 由核項加平均庫侖排斥構成，是平均場的數學化身。
  • 單電子方程： \(\left[ -\Delta/2 + V^{SCF}_i \right] \phi_i = \varepsilon_i \phi_i\) 。
  • 庫侖積分： \(J_{ij}\) 描述電荷雲排斥，總能量須扣 \(1/2\) 重複計算。
  • 迭代工具：打靶法解徑向方程、 Poisson 求解器構造勢、阻尼混合穩定迴圈。
  • 可重演實驗：取氦 \(Z = 2\) ，猜類氫軌道，迭代三輪 \(V^{SCF}\) ，觀察 \(\varepsilon_{1s}\) 如何收斂。
  • 延伸閱讀：Hartree 1928 Proc. Cambridge；後續案件直通 Hartree-Fock，因為缺的交換項已有 Fock 與 Slater 在門外等候。