哈特里 福克方程，又稱為 HF 方程，是一個應用變分法計算波函式的方程式，是量子物理、凝聚態物理學、量子化學中最重要的方程之一。HF 方程形式上是單電子本徵方程，求得的本徵態是單電子波函式，即分子軌道。以 HF 方程為核心的數值計算方法稱為 哈特里 福克方法。 基於分子軌道理論的所有量子化學計算方法都是以 HF 方法為基礎的。鑑於分子軌道理論在現代量子化學中的廣泛應用，HF 方程被視為現代量子化學的基石。 哈特里 福克近似 哈特里 福克近似也稱為分子軌道近似或單行列式近似 ，認為多電子體系波函式可以由體系分子軌道波函式構造的單個斯萊特行列式表示： Psifrac 1sqrt Nbegin vmatrix chi_1 1chi_2 1 cdots chi_N 1 chi_1 2chi_2 2 cdots chi_N 2 vdotsvdotsddotsvdots chi_1 Nchi_2 N cdots chi_N N end vmatrixchi_1chi_2cdotschi_Nrangle 這裡 N 為電子數，chi_i x 是能量最低的 N 個分子軌道波函式。而分子軌道波函式是由體系中所有原子軌道波函式經過線性組合構成的， chi_a xsum_i C_aiphi_i x 通常原子軌道波函式 phi_i x 由計算時指定的基組給出。 在此近似下，由變分法，體系的多電子薛定諤方程可化為單電子方程 f_1chi_a1varepsilon_achi_a1 其中 f_1 為福克算符，chi_a 1 為體系的分子軌道，varepsilon_a 是分子軌道 chi_a 1 的軌道能。這就是哈特里 福克方程。求解該方程，即求解使方程成立的一組分子軌道；這組分子軌道使得在 HF 近似下的體系總能量最低。 發展簡史 1927 年，物理學家瓦爾特海特勒和弗裡茨倫敦完成了氫氣分子的量子力學計算之後，開啟了量子化學的時代。從此，人們便開始嘗試使用量子力學理論來解釋化學物質結構和化學現象。 為了解決多電子體系薛定諤方程近似求解的問題，量子化學家道格拉斯哈特里在 1928 年提出了哈特里假設，將每個電子看做是在其他所有電子構成的平均勢場中運動的粒子，並且首先提出了迭代法的思路。哈特里根據他的假設，將體系電子哈密頓運算元分解為若干個單電子哈密頓運算元的簡單代數和，每個單電子哈密頓運算元中只包含一個電子的座標，因而體系多電子波函式可以表示為單電子波函式的簡單乘積，這就是哈特里方程。但是由於哈特里沒有考慮電子波函式的反對稱要求，他的哈特里方程實際上是非常不成功的。 1930 年，哈特里的學生弗拉基米爾福克和分別提出了考慮泡利原理的自洽場迭代方程和單行列式型多電子體系波函式，這就是今天的哈特里 福克方程。但是由於計算上的困難，HF 方程誕生後整整沉寂了二十年，在 1950 年，量子化學家想到將分子軌道用原子軌道的線性組合來近似展開，而得到了結構的羅特漢方程。 1953 年，美國的、和英國的約翰波普爾花費兩年時間使用手搖計算器分別獨立地實現了對氮氣分子的 RHF 自洽場計算，這是人類首次透過求解 HF 方程獲得對化學結構的量子力學解釋，也是量子化學計算方法第一次實際完成。在第一次成功之後，伴隨著電腦技術的迅猛發展，HF 方程與量子化學一道獲得長足發展，在 HF 方程的基礎上，人們發展出了高階量子化學計算方法，使得計算精度進一步提高，透過對 HF 方程電子積分的簡化和引數化，人們大大縮減了量子化學的計算量，使得對超過 1000 個原子的中等大小分子的計算成為可能。 方程的匯出 哈特里 福克方程源於對多電子體系電子波函式的變分法處理。在玻恩 奧本海默近似條件下，一個多電子體系的電子運動與能量可以與原子核的運動和能量相互分離，這樣利用電子哈密頓算符和多電子波函式便可以計算體系的電子能量。其能量的表示式為 E_0langlePsi_0H_elePsi_0rangle 式中 E_0 表示體系基態電子能量；H_ele 表示體系的電子哈密頓算符， H_elesum_iN frac 12nabla_i2 sum_i1Nsum_a1M frac Z_aR_iasum_i1Nsum_jiN frac 1r_ij 根據作用方式，可以將 H_ele 分解為兩部分 H_eleO_1O_2 其中 O_1 為單電子算符 O_1sum_iN frac 12nabla_i2 sum_i1Nsum_a1M frac Z_aR_iasum_iN left frac 12nabla_i2 sum_a1M frac Z_aR_ia rightsum_iN h_i 描述單個電子的動能和原子核吸引勢能；而 O_2 為二電子算符 O_2sum_i1Nsum_jiN frac 1r_ij 描述電子間相互作用。 Psi_0 代表基態多電子波函式，是體系單電子分子軌道波函式為基函式組建的斯萊特行列式。構建 Psi_0 的各個分子軌道相互之間是正交歸一的，即約束條件為 langlechi_achi_brangledelta_ab 考慮此約束條件，應用拉格朗日乘數法對函式 LE_0 chi_isum_abN varepsilon_ableft langle abrangledelta_abright 變分求極值。式中 varepsilon_ab 是拉格朗日待定因子，langle abrangle 是 langlechi_achi_brangle 的簡寫。 變分法的處理過程如下：令 delta Ldelta E_0 chi_isum_abNvarepsilon_abdelta chi_achi_b0 其中 begin aligndelta E_0sum_aN delta chi_ahchi_asum_aN chi_ahdelta chi_a frac 12sum_abN left deltachi_achi_achi_bchi_bchi_adeltachi_achi_bchi_bchi_achi_adeltachi_bchi_bchi_achi_achi_bdeltachi_bright frac 12sum_abN left deltachi_achi_bchi_bchi_achi_adeltachi_bchi_bchi_achi_achi_bdeltachi_bchi_achi_achi_bchi_bdeltachi_arightend align 這裡定義單電子積分記號 chi_achi_blanglechi_achi_brangleint dxchi_axchi_b x 以及二電子積分記號 chi_achi_bchi_cchi_dint dx_1dx_2chi_ax_1chi_b x_1 r_121chi_cx_2chi_d x_2 考慮到流動座標的不可區分性，可以簡化為： begin aligndelta E_0sum_aN delta chi_ahchi_asum_abN left deltachi_achi_achi_bchi_bdeltachi_achi_bchi_bchi_aright left sum_aN delta chi_ahchi_asum_abN left deltachi_achi_achi_bchi_bdeltachi_achi_bchi_bchi_arightrightend align 同理，delta L 中的 sum_abNvarepsilon_badelta ab 項有： sum_abNvarepsilon_badelta absum_abNvarepsilon_badeltachi_achi_bleft sum_abNvarepsilon_badeltachi_achi_bright 將兩項相加，delta L 表示為： begin aligndelta Lsum_aN delta chi_ahchi_asum_abN left deltachi_achi_achi_bchi_bdeltachi_achi_bchi_bchi_arightsum_abNvarepsilon_badeltachi_achi_b left sum_aN delta chi_ahchi_asum_abN left deltachi_achi_achi_bchi_bdeltachi_achi_bchi_bchi_arightsum_abNvarepsilon_badeltachi_achi_bright 0end align 若 L 函式處於極值點，則變數 leftchi_iright 向各個方向的微小變化都應該有 delta L0。可以取 leftdelta chi_irightleftidelta chi_iright 沿虛軸變分，則在 delta L 表示式中，第一項前會產生一個 i 的係數，第二項複共軛會產生一個 i 係數： begin aligndelta Lisum_aN delta chi_ahchi_asum_abN left deltachi_achi_achi_bchi_bdeltachi_achi_bchi_bchi_arightsum_abNvarepsilon_badeltachi_achi_b ileft sum_aN delta chi_ahchi_asum_abN left deltachi_achi_achi_bchi_bdeltachi_achi_bchi_bchi_arightsum_abNvarepsilon_badeltachi_achi_bright 0end align 消去虛數單位，兩式相加，可以消去表示式中的複共軛項： delta Lsum_aN delta chi_ahchi_asum_abN left deltachi_achi_achi_bchi_bdeltachi_achi_bchi_bchi_arightsum_abNvarepsilon_badeltachi_achi_b0 在引入庫侖算符和交換算符的概念之後，上述表示式可以改寫為： sum_aN int dx_1 delta chi_a1left h_1chi_a1sum_bNleft J_b 1K_b 1rightchi_a1sum_bnvarepsilon_bachi_b1right0 對任意 leftdeltachi_iright 上述等式均應成立，因而必須有 h_1chi_a1sum_bNleft J_b 1K_b 1rightchi_a1sum_bnvarepsilon_bachi_b10 整理得到： left h_1sum_bNleft J_b 1K_b 1rightrightchi_a1sum_bNvarepsilon_bachi_b1 定義福克算符 f_1h_1sum_bNleft J_b 1K_b 1right 方程可以表達為 f_1chi_a1sum_bNvarepsilon_bachi_b1 即哈特里 福克方程。為了求解，透過對分子軌道波函式進行酉變換處理，使得軌道能矩陣 varepsilon_ba 對角化，將一般的不可解的 HF 方程轉化為正則哈特里 福克方程： f_1chi_a1epsilon_achi_a1 此方程形式上為本徵方程，但是福克算符中的庫侖算符和交換算符都與分子軌道有關，因此只能夠透過自洽迭代的方法近似求解，即哈特里 福克自洽場（HFSCF）方法。HFSCF 方法是組態相互作用方法、多體微擾理論、等現代量子化學計算方法的基礎。 方程的解與自洽場方法 正則哈特里 福克方程雖然具有簡單的本徵方程形式，但福克運算元中的庫侖運算元和交換運算元中含有所有 chi_a 的表示式，因而方程的實際形式非常複雜，無法求得精確的解析解，只能使用迭代法求解，即量子化學中的自洽場方法。 在實際操作中，人們會首先給定基組，將本徵方程轉化為矩陣方程。這一轉變使分子軌道可以表示為 chi_a xsum_i C_aiphi_i x，可以用一個變化的 K 維向量 mathbb C_a 來代表分子軌道，K 為基組中基函式的數量。同理，可以將福克算符轉化為福克矩陣。然而由於基函式不一定正交，存在重疊積分 S_ijint dxphi_ixphi_j x 最終 HF 方程的形式轉化為 mathbb Fmathbb Cmathbb varepsilonmathbb Smathbb C 此方程稱為哈特里 福克 羅特漢方程，或羅特漢方程。這樣，相對複雜的本徵方程就轉化為只需要進行簡單代數計算就可以求解的矩陣本徵方程，而原方程中複雜的積分則在上述轉化過程中一次性完成了。 求解時首先需要給定一個分子軌道的初始猜測 mathbb C，用這套分子軌道計算福克矩陣 mathbb F，再求解哈特里 福克 羅特漢方程，得到一組新的分子軌道 mathbb C 和相應的本徵值矩陣 mathbb varepsilon，再使用新軌道 mathbb C 計算新的福克矩陣 mathbb F 重複上述過程，直到 mathbb C 的各個矩陣元數值不再有明顯的變化，稱作自洽迭代收斂，就得到了 HF 方程的解。 得到收斂的 mathbb C 矩陣後，將這些係數代入 chi_asum_i C_iaphi_i 與基函式結合，便獲得了最終的分子軌道波函式以及體系電子總能量等各種性質。 方程的應用 HF 方程在量子化學中有著廣泛的應用，所有分子軌道理論的量子化學計算都是以 HF 方程為基礎的。 組態相互作用方法（CI）：在 CI 方法中，透過 HF 方程解得的一系列分子軌道用於構建多電子基函式集，在構建了多電子基函式集後再透過變分法處理獲得 CI 能量的最低點，因而進行 CI 計算必須首先完成 HF 方程的求解。 多體微擾理論方法（MPn）：MPn 計算將體系所有福克算符的代數和定義為哈特里 福克哈密頓算符，將電子間相互作用所產生的能量項看作是對哈特里 福克哈密頓算符的微擾，經過微擾處理後可以獲得體系能量的近似值。進行多體微擾計算也需要首先進行 HF 方程的求解，以獲得需要的分子軌道波函式形式和分子軌道能量。 ：半經驗量子化學計算是對 HF 方程求解過程的簡化。在 HF 方程的求解中，絕大部分計算量都分佈在由正則 HF 方程向矩陣本徵方程形式轉變的過程中，如果將這一過程中大量的電子積分用經驗數值代替，便可以極大地縮短 HF 方法的計算時間。為此，針對不同的研究體系，量子化學家開發了不同的積分經驗常數，與之相應地產生了各種半經驗量子化學計算方法。本質上講，半經驗計算仍然是透過自洽場方法求解 HF 方程的過程。 參見 量子化學 羅特漢方程 參考文獻 外部連結 A Numerical HartreeFock Program H H Category 方程 Category 電子結構方法