| 5 | 1/1 | 返回列表 |
| 查看: 1768 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
HF方程 已有1人參與
|
||||
| 推導(dǎo)閉殼層分子體系的HF方程 |
木蟲 (小有名氣)
專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗: +224 |
|
最簡單的閉殼層模型就是兩個電子,方便起見以這個為例子, 其中Hamilton量H=T+V(Ne)+V(ee),動能算符T=-1/2▽^2,核對電子吸引V(Ne)=ΣZ/r1,電子-電子排斥V(ee)=1/r12 根據(jù)Pauli原理,體系波函數(shù)ψ寫成Slater行列式,兩個電子1和2有兩個軌道i和j,分子體系的H不包含自旋變量并不考慮spin-orbital的耦合,可將自旋部分分離不考慮,ψ=1/√2[i(1)j(2)-i(2)j(1)]代入,方便起見先消去1/√2這個系數(shù)。 方便起見,把-1/2▽^2+ΣZ/r1寫成hi算符,<[i(1)j(2)-i(2)j(1)]|h1+h2+V(ee)|[i(1)j(2)-i(2)j(1)]>-E<[i(1)j(2)-i(2)j(1)]|[i(1)j(2)-i(2)j(1)]>=0 h1對r1微分,根據(jù)i(1)j(1)正交及i(2)j(2)正交,<[i(1)j(2)-i(2)j(1)]|h1|[i(1)j(2)-i(2)j(1)]>=<i(1)|h1i(1)>+<j(1)|h(1)|j(1)> h2對r2微分,根據(jù)i(1)j(1)正交及i(2)j(2)正交,<[i(1)j(2)-i(2)j(1)]|h2|[i(1)j(2)-i(2)j(1)]>=<i(2)|h2i(2)>+<j(2)|h(2)|j(2)> <ψ|Vee|ψ>=<[i(1)j(2)-i(2)j(1)]|1/r12|[i(1)j(2)-i(2)j(1)]>=[<i(1)j(2)|1/r12|i(1)j(2)>+<i(2)j(1)|1/r12|i(2)j(1)>]-[<i(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|i(1)j(2)>] 令上式第一個[]=J=∫ρi(1)ρj(2)/r12dr1dr2+∫ρi(2)ρj(1)/r12dr1dr2,第二個[]=K=<i(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|i(1)j(2)> 對上式做i、j歸一化約束條件的變分,構(gòu)造泛函Ω=<i(1)|h1i(1)>+<j(1)|h(1)|j(1)>+<i(2)|h2i(2)>+<j(2)|h(2)|j(2)>+J-K-Ei<i(1)|i(1)>-Ej<j(1)|j(1)>-Ei<i(2)|i(2)>-Ej<j(2)|j(2),其中Ei和Ej為Lagrange乘子。 方便起見,僅挑出對i(1)變分相關(guān)的項(增量i(1)無影響的項,變分中會被約去) Ω'=<i(1)|h1i(1)>+<i(1)j(2)|1/r12|i(1)j(2)>-[<i(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|i(1)j(2)>]-Ei<i(1)|i(1)>=0 對i(1)做作增量δi(1),則Ω‘的增量為 δΩ'=<δi(1)|h1i(1)>+<i(1)|h1δi(1)]>+<δi(1)j(2)|1/r12|i(1)j(2)>+<i(1)j(2)|r12|δi(1)j(2)>+[<δi(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|δi(1)j(2)>]-Ei<δi(1)|i(1)>-Ei<i(1)|δi(1)> 當(dāng)δΩ'=0時,H為極值,挑出復(fù)共軛項<δi(1)| <δi(1)|h1i(1)>+<δi(1)j(2)|1/r12|i(1)j(2)>+<δi(1)j(2)|1/r12|i(2)j(1)>-Ei<δi(1)|i(1)>=0 由于變分過程中δi(1)的任意性,h1i(1)+<j(2)|1/r12|j(2)>i(1)+<j(2)|1/r12|i(2)>j(1)=Eii(1) 其中,第一項包含動能及核對電子吸引hi,第二項為電子-電子平均排斥Ji,第三項為交換勢Ki,總體構(gòu)成Fock算符Fi,右側(cè)為i能級能量Ei。 打字打的比較費(fèi)勁,其實(shí)直接參考一些量子力學(xué)和量子化學(xué)的書就可以了,上面公式是自己推導(dǎo)的,可能個別之處有不嚴(yán)謹(jǐn)甚至是出錯的地方,由于本人數(shù)學(xué)及量子化學(xué)基礎(chǔ)不太好,希望大家能夠批評指正。 |

專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗: +224 |
|
1、這里我確實(shí)沒說清楚,h1中動能項1是包含對r1的微分,以及后面那一項吸引項也是r1的函數(shù),其實(shí)能夠求單粒子平均是因為這兩個算符都是單粒子的。 2、只是為了簡單的把交換勢導(dǎo)出來,要不然多個粒子推導(dǎo)會有些繁復(fù),確實(shí)加上個kronecker δ(ms,ms)會更嚴(yán)格一些。如你所說,實(shí)際上兩自旋相反的電子這一項是不該存在的。不過這個問題確實(shí)是我忘記說了。 3、因為式子太長了,再寫上一堆根號可能會有點(diǎn)亂。 寫得有點(diǎn)快,水平有點(diǎn)低,難免會出錯,不過多謝指正。 |


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +22 | 阿斯蒂芬2004 2026-03-28 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 08工科,295,接受跨專業(yè)調(diào)劑 +3 | lmnlzy 2026-03-31 | 3/150 |
|
|
[考研] 英一數(shù)一總分334求調(diào)劑 +3 | 陳陽坤 2026-03-31 | 3/150 |
|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +5 | 里子木yy 2026-03-29 | 5/250 |
|
|
[考研] 282求調(diào)劑 +3 | 呼吸都是減肥 2026-03-31 | 3/150 |
|
|
[考研] 一志愿西交288化工專碩求調(diào)劑 +4 | 好運(yùn)好運(yùn)接接 2026-03-24 | 4/200 |
|
|
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
[考研] 0703一志愿9,初試成績:338,四六級已過,有科研經(jīng)歷,求調(diào)劑! +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
[基金申請] 面上5B能上會嗎? +7 | redcom 2026-03-29 | 7/350 |
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|