哈特里-福克方程式(英語:Hartree–Fock equation),又稱為HF方程式,是一個應用變分法計算多電子系統波函數的方程式,是量子物理、凝聚態物理學、量子化學中最重要的方程式之一。HF方程式形式上是單電子本徵方程式,求得的本徵態是單電子波函數,即分子軌態。以HF方程式為核心的數值計算方法稱為「哈特里-福克方法」(Hartree–Fock method)。
基於分子軌態理論的所有量子化學計算方法都是以HF方法為基礎的。鑑於分子軌態理論在現代量子化學中的廣泛應用,HF方程式被視為現代量子化學的基石。
1927年,物理學家瓦爾特·海特勒和弗里茨·倫敦完成了氫氣分子的量子力學計算之後,開啟了量子化學的時代。從此,人們便開始嘗試使用量子力學理論來解釋化學物質結構和化學現象。
為了解決多電子體系薛定諤方程式近似求解的問題,量子化學家道格拉斯·哈特里在1928年提出了哈特里假設,將每個電子看做是在其他所有電子構成的平均勢場中運動的粒子,並且首先提出了迭代法的思路。哈特里根據他的假設,將體系電子哈密頓算子分解為若干個單電子哈密頓算子的簡單代數和,每個單電子哈密頓算子中只包含一個電子的坐標,因而體系多電子波函數可以表示為單電子波函數的簡單乘積,這就是哈特里方程式。但是由於哈特里沒有考慮電子波函數的反對稱要求,他的哈特里方程式實際上是非常不成功的。
1930年,哈特里的學生弗拉基米爾·福克和約翰·斯萊特分別提出了考慮鮑利原理的自洽場迭代方程式和單行列式型多電子體系波函數,這就是今天的哈特里-福克方程式。但是由於計算上的困難,HF方程式誕生後整整沉寂了二十年,在1950年,量子化學家克萊門斯·羅特漢想到將分子軌態用原子軌態的線性組合來近似展開,而得到了閉殼層結構的羅特漢方程式。
1953年,美國的R.帕里瑟、R.帕爾和英國的約翰·波普爾花費兩年時間使用手搖計算器分別獨立地實現了對氮氣分子的RHF自洽場計算,這是人類首次通過求解HF方程式獲得對化學結構的量子力學解釋,也是量子化學計算方法第一次實際完成[2]。在第一次成功之後,伴隨着電腦技術的迅猛發展,HF方程式與量子化學一道獲得長足發展,在HF方程式的基礎上,人們發展出了高級量子化學計算方法,使得計算精度進一步提高,通過對HF方程式電子積分的簡化和參數化,人們大大縮減了量子化學的計算量,使得對超過1000個原子的中等大小分子的計算成為可能。
哈特里-福克方程式源於對多電子體系電子波函數的變分法處理。在玻恩-奧本海默近似法條件下,一個多電子體系的電子運動與能量可以與原子核的運動和能量相互分離,這樣利用電子哈密頓算符和多電子波函數便可以計算體系的電子能量。其能量的表達式為:
![{\displaystyle E_{0}=\langle \Psi _{0}|H_{ele}|\Psi _{0}\rangle }](//wikimedia.org/api/rest_v1/media/math/render/svg/f153a6f7695d68f7454a5dedc48b882decc75bd1)
式中
表示體系基態電子能量;
表示體系的電子哈密頓算符,
![{\displaystyle H_{ele}=\sum _{i}^{N}{\frac {-1}{2}}\nabla _{i}^{2}-\sum _{i=1}^{N}\sum _{a=1}^{M}{\frac {Z_{a}}{R_{ia}}}+\sum _{i=1}^{N}\sum _{j>i}^{N}{\frac {1}{r_{ij}}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/75e2a28a5761c96f9c4675953948c87e09071d3a)
根據作用方式,可以將
分解為兩部分
![{\displaystyle H_{ele}=O_{1}+O_{2}}](//wikimedia.org/api/rest_v1/media/math/render/svg/5146251c5cdd2016539871bd11e903697220632d)
其中
為單電子算符
![{\displaystyle O_{1}=\sum _{i}^{N}{\frac {-1}{2}}\nabla _{i}^{2}-\sum _{i=1}^{N}\sum _{a=1}^{M}{\frac {Z_{a}}{R_{ia}}}=\sum _{i}^{N}-\left({\frac {1}{2}}\nabla _{i}^{2}+\sum _{a=1}^{M}{\frac {Z_{a}}{R_{ia}}}\right)=\sum _{i}^{N}h_{i}}](//wikimedia.org/api/rest_v1/media/math/render/svg/d302b997c96fb8051c1f14c0b0e8b7b60f2a4bd5)
描述單個電子的動能和原子核吸引位能;而
為二電子算符
![{\displaystyle O_{2}=\sum _{i=1}^{N}\sum _{j>i}^{N}{\frac {1}{r_{ij}}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/2aad17980c7905867e75e6d5268a4b9fa0a73fb7)
描述電子間相互作用。
代表基態多電子波函數,是體系單電子分子軌態波函數為基函數組建的斯萊特行列式。構建
的各個分子軌態相互之間是正交歸一的,即約束條件為
![{\displaystyle \langle \chi _{a}|\chi _{b}\rangle =\delta _{ab}}](//wikimedia.org/api/rest_v1/media/math/render/svg/6a0ea991789da7d412082bdb24c77e96edae90db)
考慮此約束條件,應用拉格朗日乘數法對函數
變分求極值。式中
是拉格朗日待定因子,
是
的簡寫。
變分法的處理過程如下:令
![{\displaystyle \delta L=\delta E_{0}[\{\chi _{i}\}]-\sum _{a,b}^{N}\varepsilon _{ab}\delta [\chi _{a}|\chi _{b}]=0}](//wikimedia.org/api/rest_v1/media/math/render/svg/b0bedf76d2098845939029a7329d9d26d596d2f1)
其中
![{\displaystyle {\begin{aligned}\delta E_{0}=&\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a}^{N}[\chi _{a}|h|\delta \chi _{a}]\\&+{\frac {1}{2}}\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]+[\chi _{a}\delta \chi _{a}|\chi _{b}\chi _{b}]+[\chi _{a}\chi _{a}|\delta \chi _{b}\chi _{b}]+[\chi _{a}\chi _{a}|\chi _{b}\delta \chi _{b}]\right)\\&+{\frac {1}{2}}\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]+[\chi _{a}\delta \chi _{b}|\chi _{b}\chi _{a}]+[\chi _{a}\chi _{b}|\delta \chi _{b}\chi _{a}]+([\chi _{a}\chi _{b}|\chi _{b}\delta \chi _{a}]\right)\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/6ccc819e62bd6c09e8e72f5a2244bf067c81e634)
這裏定義單電子積分記號
![{\displaystyle [\chi _{a}|\chi _{b}]=\langle \chi _{a}|\chi _{b}\rangle =\int dx\chi _{a}^{*}(x)\chi _{b}(x)}](//wikimedia.org/api/rest_v1/media/math/render/svg/2a79fc02d4dce206eea49a0f31ae9d65d9836c1e)
以及二電子積分記號
![{\displaystyle [\chi _{a}\chi _{b}|\chi _{c}\chi _{d}]=\int dx_{1}dx_{2}\chi _{a}^{*}(x_{1})\chi _{b}(x_{1})r_{12}^{-1}\chi _{c}^{*}(x_{2})\chi _{d}(x_{2})}](//wikimedia.org/api/rest_v1/media/math/render/svg/13986ee99dfbe01cf535108cda382d9db85ecfec)
考慮到流動坐標的不可區分性,可以簡化為:
![{\displaystyle {\begin{aligned}\delta E_{0}=&\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]-[\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]\right)\\&+\left(\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]-[\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]\right)\right)^{*}\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/bde6fa8ded659c8d11e013a6d20b4801ca742739)
同理,
中的
項有:
![{\displaystyle \sum _{a,b}^{N}\varepsilon _{ba}\delta [a|b]=\sum _{a,b}^{N}\varepsilon _{ba}[\delta \chi _{a}|\chi _{b}]+\left(\sum _{a,b}^{N}\varepsilon _{ba}[\delta \chi _{a}|\chi _{b}]\right)^{*}}](//wikimedia.org/api/rest_v1/media/math/render/svg/d2a8b3ac059c08140650b230b95197f0c6d1c931)
將兩項相加,
表示為:
![{\displaystyle {\begin{aligned}\delta L=&\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]-[\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]\right)+\sum _{a,b}^{N}\varepsilon _{ba}[\delta \chi _{a}|\chi _{b}]\\&+\left(\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]-[\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]\right)+\sum _{a,b}^{N}\varepsilon _{ba}[\delta \chi _{a}|\chi _{b}]\right)^{*}\\=&0\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/dfd6d56db76c7362450648b448a00310bfa95c00)
若
函數處於極值點,則變量
向各個方向的微小變化都應該有
。可以取
沿虛軸變分,則在
表達式中,第一項前會產生一個
的系數,第二項複共軛會產生一個
系數:
![{\displaystyle {\begin{aligned}\delta L=&i\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]-[\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]\right)+\sum _{a,b}^{N}\varepsilon _{ba}[\delta \chi _{a}|\chi _{b}]\\&-i\left(\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]-[\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]\right)+\sum _{a,b}^{N}\varepsilon _{ba}[\delta \chi _{a}|\chi _{b}]\right)^{*}\\=&0\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/b79c7000d1d68d49c6f8a33b6e4e5c17948cf835)
消去虛數單位,兩式相加,可以消去表達式中的複共軛項:
![{\displaystyle \delta L=\sum _{a}^{N}[\delta \chi _{a}|h|\chi _{a}]+\sum _{a,b}^{N}\left([\delta \chi _{a}\chi _{a}|\chi _{b}\chi _{b}]-[\delta \chi _{a}\chi _{b}|\chi _{b}\chi _{a}]\right)+\sum _{a,b}^{N}\varepsilon _{ba}[\delta \chi _{a}|\chi _{b}]=0}](//wikimedia.org/api/rest_v1/media/math/render/svg/3cbc6475fff2f5662e02be454a0e37f822dd12e2)
在引入庫侖算符和交換算符的概念之後,上述表達式可以改寫為:
![{\displaystyle \sum _{a}^{N}\int dx_{1}\delta \chi _{a}^{*}(1)\left[h_{1}\chi _{a}(1)+\sum _{b}^{N}\left(J_{b(1)}-K_{b(1)}\right)\chi _{a}(1)-\sum _{b}^{n}\varepsilon _{ba}\chi _{b}(1)\right]=0}](//wikimedia.org/api/rest_v1/media/math/render/svg/8f732dee57c75a92e553a0b37e8706d94b230d45)
對任意
上述等式均應成立,因而必須有:
![{\displaystyle h_{1}\chi _{a}(1)+\sum _{b}^{N}\left(J_{b(1)}-K_{b(1)}\right)\chi _{a}(1)-\sum _{b}^{n}\varepsilon _{ba}\chi _{b}(1)=0}](//wikimedia.org/api/rest_v1/media/math/render/svg/7373f7b560f70a5aa5d1647870321c3caeb6ce3a)
整理得到:
![{\displaystyle \left[h_{(1)}+\sum _{b}^{N}\left(J_{b(1)}-K_{b(1)}\right)\right]\chi _{a}(1)=\sum _{b}^{N}\varepsilon _{ba}\chi _{b}(1)}](//wikimedia.org/api/rest_v1/media/math/render/svg/75430c638273f985d1700acb098a9b686ae94ea9)
定義福克算符
![{\displaystyle f_{(1)}=h_{(1)}+\sum _{b}^{N}\left(J_{b(1)}-K_{b(1)}\right)}](//wikimedia.org/api/rest_v1/media/math/render/svg/fc95a98380e251fd8c2d5bec02041123d7d01053)
方程式可以表達為
![{\displaystyle f_{(1)}\chi _{a}(1)=\sum _{b}^{N}\varepsilon _{ba}\chi _{b}(1)}](//wikimedia.org/api/rest_v1/media/math/render/svg/3572494743e8b3027bc3f1335095b3bb947e4d0f)
即哈特里-福克方程式。為了求解,通過對分子軌態波函數進行酉轉換處理,使得軌態能矩陣
對角化,將一般的不可解的HF方程式轉化為正則哈特里-福克方程式:
![{\displaystyle f_{(1)}\chi _{a}(1)=\epsilon _{a}\chi _{a}(1)}](//wikimedia.org/api/rest_v1/media/math/render/svg/5ed77949e8a4595ef238a70a42ff830485cac4ab)
此方程式形式上為本徵方程式,但是福克算符中的庫侖算符和交換算符都與分子軌態有關,因此只能夠通過自洽迭代的方法近似求解,即哈特里-福克自洽場(HF-SCF)方法。HF-SCF方法是排佈相互作用方法、多體微擾理論、半經驗量子化學計算等現代量子化學計算方法的基礎。
HF方程式在量子化學中有着廣泛的應用,所有分子軌態理論的量子化學計算都是以HF方程式為基礎的。
- 排佈相互作用方法(CI):在CI方法中,通過HF方程式解得的一系列分子軌態用於構建多電子基函數集,在構建了多電子基函數集後再通過變分法處理獲得CI能量的最低點,因而進行CI計算必須首先完成HF方程式的求解。
- 多體微擾理論方法(MPn):MPn計算將體系所有福克算符的代數和定義為哈特里-福克哈密頓算符,將電子間相互作用所產生的能量項看作是對哈特里-福克哈密頓算符的微擾,經過微擾處理後可以獲得體系能量的近似值。進行多體微擾計算也需要首先進行HF方程式的求解,以獲得需要的分子軌態波函數形式和分子軌態能量。
- 半經驗量子化學計算:半經驗量子化學計算是對HF方程式求解過程的簡化。在HF方程式的求解中,絕大部分計算量都分佈在由正則HF方程式向矩陣本徵方程式形式轉變的過程中,如果將這一過程中大量的電子積分用經驗數值代替,便可以極大地縮短HF方法的計算時間。為此,針對不同的研究體系,量子化學家開發了不同的積分經驗常數,與之相應地產生了各種半經驗量子化學計算方法。本質上講,半經驗計算仍然是通過自洽場方法求解HF方程式的過程。
Pople, JA; Pariser, R; Parr, RG; Koutecky, J; Murrell, JN; Harget, AJ. Semi-Empirical Self-Consistent-Field Molecular Orbital Theory of Molecules. Trans. Faraday Soc. 1953, 49: 1375.