量子化學第七章-自洽場分子軌道理論

《量子化學》第七章自洽場分子軌道理論Chapter7Self-consistentFieldMolecularOrbitalTheory17.1分子體系7.2單粒子模型7.3哈特利－?？?Hartree-Fock)－羅湯方程27.1分子體系

1．分子體系的薛定諤方程分子由若干原子核和電子組成的多粒子體系。

要嚴格地寫出多電子體系地哈密頓算符是很困難的,因為多電子體系中的相互作用項很多,如:(1)核與電子間的吸引作用能;(2)電子與電子間的排斥作用能;(3)核與核的排斥能3(4)交換能,電子具有全同粒子特性,又得滿足保里原理。則下列兩種狀態(tài)的波函數(shù)不同，能量也有所差異。(a)非對稱波函數(shù)(b)反對稱波函數(shù)例：Li+的某一激發(fā)態(tài)1s12s1,假設電子排布狀態(tài)為1s2s4(5)自旋—軌道偶合作用能;(8)其它效應有關的作用能。(7)軌道—軌道偶合作用能;(6)自旋—自旋偶合作用能;在上述作用項中,前三項是主要的.分子體系的哈密頓算符應當包括所有電子和原子核的動能和勢能,即:5電子動能算符核動能算符核-核排斥能算符電子-電子排斥能算符核-電子作用能算符式中n―電子數(shù)目，N―核的數(shù)目。96事實上，這并不是分子嚴格的哈密頓算符，它僅僅考慮了分子中各個電荷間的庫侖相互作用，而沒有考慮電荷之間一般的電磁作用，也沒有考慮自旋與自旋、自旋與軌道、軌道與軌道之間的相互作用，更沒有考慮電子運動的相對論效應。上述幾項為哈密頓算符的主要項，通常能獲得較好的結果。7在原子單位制中，me

=1，e

=1，a0=1，=27.21165eV=2624.54kJ/mol1原子單位長度=a0

（玻爾半經(jīng)，0.529?)1原子單位質量=me（電子的質量9.1*10-31Kg）1原子單位電量=e（電子的電量1.6*10-19C）1原子單位能量=1個hartree能量原子單位制，其基本物理量有四個：68應用原子單位制，(1)變成：通過求解體系的薛定諤方程，獲得體系的各種性質。92.玻恩－奧潘海默(Born-oppenheimer)近似玻恩－奧潘海默(Born-oppenheimer)近似，也稱核固定近似，即假定核固定不動，來研究電子的運動。對于復雜的多粒子體系，上述(2)的哈密頓算符構成的薛定諤方程是很難求解的，為此，需要采取某些合理的近似。這一假定的依據(jù)是核運動速度遠低于電子的運動速度，前者約103m/s,后者約106~7m/s。10在玻恩－奧潘海默(Born-oppenheimer)近似下，上述(2)中的核的動能算符可以去掉，核與核的排斥能可視為常數(shù)，以I表示體系的薛定諤方程變成：設:1311則：只需求解電子的薛定諤方程(3)這樣，研究一個分子內部運動的問題，就變?yōu)橛懻搉個電子在固定核場中運動的問題。而電子是全同粒子，故分子結構問題的研究就轉化為n個全同粒子體系的問題。其中為體系中電子的哈密頓，為電子的能量，體系的能量為目錄13127.2單粒子模型

1．平均場的概念包含著不能分離的兩個電子的坐標，給方程的求解帶來了困難。方程(3)依然是難于求解的，因為存在所以對多電子體系，除了前面提到的核固定近似外,還引入了單電子近似，也稱獨立電子運動模型、單粒子模型近似。1113其基本思想是：認為每個電子均在諸原子核和其它電子構成的有效平均場中獨立地運動。

電子i在平均場中具有的勢能為：由定態(tài)薛定諤方程確定的定態(tài)波函數(shù)將給出電子在空間的穩(wěn)定分布，從而在整個分子空間出現(xiàn)一個穩(wěn)定的電荷分布，由此將產生一個穩(wěn)定不變的電場－靜電場。分子中是否存在著某種平均場？142．單粒子模型的哈密頓基于平均場的概念，體系中電子的哈密頓可寫成：由于波函數(shù)具有統(tǒng)計意義，所以該靜電場是一個平均場?，F(xiàn)代分子軌道計算中所采用的物理模型建立在三大近似基礎上:核固定近似、非相對論近似、單電子近似.2915令則體系中電子的哈密頓:16

由于，故體系中電子的薛定諤方程就可以用微擾理論來處理。電子瞬時相互作用與平均相互作用之差的哈密頓之和。n個彼此獨立的電子的哈密頓之和，可以作為體系的零級近似哈密頓，式中，i電子的哈密頓173．零級哈密頓的本征方程為n個電子的哈密頓之和，所以上述方程的解可以分解成：18則可通過求解多電子體系中電子i的Schr?dinger方程（如下）獲得。(6)綜上分析可知，在單電子近似下，求解n個粒子體系的薛定諤方程的問題，就歸結為求解一個單粒子的薛定諤方程的問題。31194．單電子薛定諤方程的解求解單電子薛定諤方程(6),可得一系列單電子的能量和相應的本征函數(shù)(一般說來應當有無限多個)這些單電子的本征函數(shù)(即軌道)描述著電子在分子平均場中的運動狀態(tài).20相應的分子零級近似能量為式中為電子i的空間坐標.假定體系中電子

i在軌道中運動,能量為,則整個分子的零級波函數(shù)為:21以上討論的分子軌道沒有考慮電子的自旋,實際上,電子的自旋在分子結構中起著重要的作用,所以分子軌道應包括電子的自旋態(tài).電子是自旋量子數(shù)的粒子,其自旋狀態(tài)由自旋磁量子數(shù)ms來表征.代表取向上自旋態(tài)

,則表示取向下自旋態(tài).22qi―電子i的空間坐標和自旋坐標.考慮電子同時作軌道運動和自旋運動,電子i的完全波函數(shù)為軌道波函數(shù)與自旋波函數(shù)的乘積,即綜上,考慮自旋后整個體系的零級波函數(shù)為:235.能量本征值的交換簡并電荷、質量、自旋完全相同的粒子,稱全同粒子.如:電子。全同粒子體系的完全波函數(shù)必須是對稱的或反對稱的.其對稱性具體取決于體系中粒子的自旋特性.對于自旋量子數(shù)的電子體系的完全波函數(shù)具有反對稱性,這是保里原理的要求.24反對稱性的波函數(shù)可用Slater行列式構建:電子的單電子完全波函數(shù)。256.零級能量本征值的自旋簡并對于含n(偶數(shù))個電子的閉殼層體系,每個軌道上排兩個自旋相反的電子,則n個電子占有n/2個軌道,體系的零級波函數(shù)為:對于單電子能量對應的軌道,上面可能有兩個自旋不同的電子其完全波函數(shù)分別26上式可簡寫成:27含n(奇數(shù))電子體系的基態(tài)開殼層,存在兩種簡并態(tài)完全波函數(shù)分別為:其它開殼層體系可作類似處理.則體系的零級波函數(shù)可寫成兩者的線性組合,即:287.求薛定諤方程的微擾方法前面討論了單粒子模型下分子體系的薛定諤方程的零級近似解.當考慮交換簡并時,零級波函數(shù)是一個Slater行列式,當考慮自旋簡并時,零級波函數(shù)一般寫成若干個Slater行列式的線性組合.基于多電子體系中電子的哈密頓可以分解成兩部分為n個電子的哈密頓之和,是電子的瞬時相互作用能和平均相互作用能之差.1529,可以用微擾法求解

微擾方法的效能取決于微擾哈密頓的大小.越小,微擾方法效能得到較好的結果,雖然這取決于平均場的選擇是否得當.如果適當?shù)剡x擇,使得電子在平均場中勢能逼近電子的瞬時相互作用能,就會很小.目錄307.3哈特利－福克(Hartree-Fock)－羅湯方程采用獨立電子運動模型處理多電子體系，必須求出單電子的軌道，才能構建體系的零級波函數(shù)。采用變分法能有效求解單粒子薛定諤方程(6)。1．變分原理由n個粒子構成的體系,設其哈密頓為.1931E0為體系的基態(tài)能量,這就是變分原理.變分過程如下:不斷改變試探波函數(shù),計算相應的值,直到不再減小,可以認為逼近了體系真實的基態(tài)能量,此時對應的試探波函數(shù)可近似認為體系的基態(tài)波函數(shù).

設

(1,2,…,n)為一品優(yōu)波函數(shù)，則在狀態(tài)下，體系的能量的平均值為：322.自洽場(Self-ConsistentField)迭代法對如下單電子薛定諤方程，通常采取自洽迭代的方法進行求解.實際計算時,除了考慮能量收斂外,也考慮波函數(shù)的收斂情況,一般說來后者更為嚴謹.33自洽場模型方法(Self-ConsistentField,簡稱SCF),是Hartre于1928年提出的.在該模型中,其它電子對i電子的排斥作用寫成i-e34則i電子的薛定諤方程可寫成:求解i

電子的Schr?dinger方程方程的前提是要知道，事實上，也是未知的。為此，Hartree提出采用自洽迭代法求解。35①先假設一套初始波函數(shù)……自洽迭代過程如下：36如此,由零級(初始)波函數(shù)一級波函數(shù)②同理，由一級波函數(shù)二級波函數(shù)，如此反復

現(xiàn)代分子軌道計算，多采用SCF法來求解多電子體系的Schr?dinger方程。③收斂的判據(jù)：能量判據(jù)波函數(shù)判據(jù)373.LCAO-MO和羅湯方程薛定諤方程經(jīng)變分法得到H-F方程，嚴格求解是有困難的,即使采取迭代自洽的方法進行求解,也是相當繁復的,所以人們就設計了若干近似方法求解。20世紀50年代,在H-F方程處理原子結構的基礎上,Roothaan提出,將分子軌道按某個基組集合展開,有效地解決了這個問題。38Roothaan提出,將分子軌道按某個基組集合展開,用有限展開項,按一定精度逼近分子軌道。這樣,對分子軌道的變分就轉變?yōu)閷φ归_系數(shù)的變分。H-F方程就從一組非線性的微積分方程轉化為一組數(shù)目有限的代數(shù)方程,只需迭代求解分子軌道組合系數(shù),這就是Hartree-Fock-Roothaan方程,它是分子軌道計算中的基本方程。39Roothaan的“LCAO-MO”在變分法中假設試探波函數(shù)為:40即:改變使得盡可能的低.41久期方程(代數(shù)方程)Fock算符矩陣Fuv單電子能量重疊積分：原子軌道重疊程度待定系數(shù)Hartree-Fock-Roothaanequation42密度矩陣元雙電子積分單電子積分單電子積分：單電子哈密頓算符43雙電子積分m=占有軌道數(shù)目44Hartree-Fock-Roothaan方程也可寫成其中：S為軌道重疊矩陣；C為軌道組合系數(shù)；

ε為能量本征值上述方程為非線性代數(shù)方程，只能采用SCF迭代來求解?，F(xiàn)代分子軌道計算解的真是Hartree-Fock-Roothaanequation，而非Schr?ding