| 5 | 1/1 | 返回列表 |
| 查看: 1769 | 回復(fù): 5 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||
[求助]
HF方程 已有1人參與
|
|||||
| 推導(dǎo)閉殼層分子體系的HF方程 |
專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |
|
1、這里我確實(shí)沒說清楚,h1中動能項(xiàng)1是包含對r1的微分,以及后面那一項(xiàng)吸引項(xiàng)也是r1的函數(shù),其實(shí)能夠求單粒子平均是因?yàn)檫@兩個算符都是單粒子的。 2、只是為了簡單的把交換勢導(dǎo)出來,要不然多個粒子推導(dǎo)會有些繁復(fù),確實(shí)加上個kronecker δ(ms,ms)會更嚴(yán)格一些。如你所說,實(shí)際上兩自旋相反的電子這一項(xiàng)是不該存在的。不過這個問題確實(shí)是我忘記說了。 3、因?yàn)槭阶犹L了,再寫上一堆根號可能會有點(diǎn)亂。 寫得有點(diǎn)快,水平有點(diǎn)低,難免會出錯,不過多謝指正。 |

專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +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)的項(xiàng)(增量i(1)無影響的項(xiàng),變分中會被約去) Ω'=<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ù)共軛項(xiàng)<δ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) 其中,第一項(xiàng)包含動能及核對電子吸引hi,第二項(xiàng)為電子-電子平均排斥Ji,第三項(xiàng)為交換勢Ki,總體構(gòu)成Fock算符Fi,右側(cè)為i能級能量Ei。 打字打的比較費(fèi)勁,其實(shí)直接參考一些量子力學(xué)和量子化學(xué)的書就可以了,上面公式是自己推導(dǎo)的,可能個別之處有不嚴(yán)謹(jǐn)甚至是出錯的地方,由于本人數(shù)學(xué)及量子化學(xué)基礎(chǔ)不太好,希望大家能夠批評指正。 |

木蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本2一志愿C9-333分,材料科學(xué)與工程,求調(diào)劑 +9 | 升升不降 2026-03-31 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +6 | 鐵憨憨123425 2026-03-31 | 6/300 |
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +10 | 26要上岸 2026-03-28 | 10/500 |
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
[考研] 347求調(diào)劑 +11 | 山頂見α 2026-03-25 | 11/550 |
|
|
[考研] 272求調(diào)劑,接受跨專業(yè)調(diào)劑! +3 | 閑魚盧 2026-03-31 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-30 | 6/300 |
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | 圖鑒212 2026-03-30 | 4/200 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +12 | Li李魚 2026-03-26 | 12/600 |
|
|
[考研] 291求調(diào)劑 +12 | Y-cap 2026-03-29 | 15/750 |
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 8/400 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|