 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 11 | 1/1 | 返回列表 |
| 查看: 1698 | 回復(fù): 10 | |||||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)8次，作者zzgyb增加金幣 8 個(gè) | |||||
| 當(dāng)前主題已經(jīng)存檔。 | |||||
zzgyb榮譽(yù)版主 (文壇精英)
|
[資源]
自旋污染
|
||||
|
什么是自旋污染？ 關(guān)于 Hartree-Fock計(jì)算（通常采用Rootaan的自洽場(chǎng)（SCF）方式）介紹性的描述集中在單重態(tài)體系，在這樣體系中所有電子自旋都已經(jīng)配對(duì)。通過(guò)假設(shè)計(jì)算是限制在，每個(gè)占據(jù)軌道有兩個(gè)電子，的計(jì)算，計(jì)算可以相對(duì)地容易進(jìn)行。這通常是指自旋限制的Hartree-Fock計(jì)算或RHF。 對(duì)于自旋多重度不為一的體系，是不可能象自旋多重度為一的體系那樣用RHF方式的。 在一個(gè)非限制性的計(jì)算中，有兩組完全的軌道軌道，一組是alpha電子的一組是beta電子的。 通常這兩組軌道用同樣的基函數(shù)但是不同的分子軌道參數(shù)。 非限制性計(jì)算的優(yōu)點(diǎn)是他們計(jì)算的效率非常高。 缺點(diǎn)是波函數(shù)不再是總自旋的本征函數(shù)，因此計(jì)算中可能會(huì)導(dǎo)入一些誤差。 這種誤差成為自旋污染。 自選污染如何影響計(jì)算結(jié)果 自旋污染導(dǎo)致在看來(lái)是我們期望的波函數(shù)中，混合進(jìn)了一點(diǎn)其他的自旋態(tài)。這種情況下由于有了較多的變分自由度，有時(shí)候會(huì)導(dǎo)致輕微的降低計(jì)算的總能量。更多的情況是因?yàn)檩^高能量態(tài)的混入導(dǎo)致總能量的輕微升高。然而這些變化是不正確的波函數(shù)的結(jié)果。因?yàn)檫@不是一種系統(tǒng)誤差，不同態(tài)間的能量差可能受到相反的影響。高自旋的污染可以影響稽核和布局?jǐn)?shù)分析，和顯著地影響自旋密度。 作為自旋污染出現(xiàn)的一種檢查，大多ab initio程序會(huì)打印出期望的總自旋值，. 如果沒(méi)有自旋污染，它會(huì)等于s(s+1)，這里s等于1/2的未成對(duì)電子數(shù)倍。由有機(jī)分子計(jì)算經(jīng)驗(yàn)得到的一個(gè)簡(jiǎn)單規(guī)律是，當(dāng) 和s(s+1) 的差別小雨10%時(shí)，自旋污染的影響可以忽略。雖然這個(gè)規(guī)則提供了一種快速檢驗(yàn)的方法，通常用實(shí)證據(jù)或更嚴(yán)格的計(jì)算來(lái)檢查你的結(jié)果是明智的。 自旋污染通常在非限制的 Hartree-Fock (UHF)和分限制性的Møller-Plesset (UMP2, UMP3, UMP4) 計(jì)算中出現(xiàn)。在DFT計(jì)算中很少會(huì)發(fā)現(xiàn)顯著的自旋污染，甚至在采用非限制性的Kohn-Sham 軌道的情況。 非限制性計(jì)算通常結(jié)合一步自旋湮滅步驟，這個(gè)步驟從計(jì)算的波函數(shù)中除去很大百分比的自旋污染。它有助于最小化自旋污染但是無(wú)法完全阻止自旋污染。的最后值是對(duì)出現(xiàn)的自旋污染的最好的檢查。 在Gaussian中，選項(xiàng)"iop(5/14=2)" 告訴程序用湮滅后的波函數(shù)產(chǎn)生布局?jǐn)?shù)分析。我不知道是否有程序在執(zhí)行幾何優(yōu)化過(guò)程中使用湮滅的波函數(shù)。 限制性開(kāi)殼成計(jì)算 有可能進(jìn)行自旋限制性開(kāi)殼層計(jì)算（ROHF）。這種方法的好處是沒(méi)有自旋污染。缺點(diǎn)是會(huì)需要額外的CPU時(shí)間在對(duì)單占據(jù)的和雙占據(jù)的軌道都進(jìn)行正確處理以及它們之間的相互作用上。作為采用的數(shù)學(xué)算法的一個(gè)結(jié)果，ROHF計(jì)算給出好的總能量和波函數(shù)但是單占據(jù)的軌道不嚴(yán)格服從Koopman定理。 當(dāng)自旋污染導(dǎo)致的誤差不可接受時(shí)，先執(zhí)行開(kāi)殼層計(jì)算是給出可靠波函數(shù)的最好方法。 在Gaussian程序中，限制性開(kāi)殼層計(jì)算可以執(zhí)行Hartree-Fock，密度泛函理論，MP2 和一些半經(jīng)驗(yàn)波函數(shù)計(jì)算。但是ROMP2方式不知吃解析梯度，因此最快的方式是在其他方式優(yōu)化的幾何上進(jìn)行單點(diǎn)能計(jì)算。如果幾何優(yōu)化一定要在這樣的理論水平下進(jìn)行，必須采用一種非基于梯度的方式，如Fletcher-Powell優(yōu)化。（注意在G94手冊(cè)中對(duì)所有情況都還沒(méi)有這種功能） 自旋投影方式 另一種近似是先運(yùn)行一種非限制性計(jì)算然后在得到的波函數(shù)上投影掉自旋污染（PUHF，PUMP2）。 自旋投影結(jié)果并不給出用限制性開(kāi)殼層計(jì)算得到的能量。 這是因?yàn)榉窍拗菩缘能壍辣粌?yōu)化來(lái)描述被污染的自旋態(tài)而不是被優(yōu)化來(lái)描述自旋投影的態(tài)。 用自旋強(qiáng)制的UHF方式（SUHF）可以得到類似的效果。在這種方式中，通過(guò)一個(gè)拉格朗日乘因子法，UHF波函數(shù)中的自旋污染誤差被限制。當(dāng)乘的因子趨于無(wú)窮時(shí)，這種方法完全地除去了自旋污染。實(shí)際上小的正值就可以除去大多數(shù)的自旋污染。 半電子近似 對(duì)自由基計(jì)算，半經(jīng)驗(yàn)程序常用半電子近似。半電子方式是RHF計(jì)算中處理一個(gè)單占據(jù)軌道的數(shù)學(xué)技術(shù)。它在擁有一個(gè)近似的波函數(shù)和軌道能量的情況下能得到一致性的總能量。因?yàn)椴捎玫氖菃涡辛惺接?jì)算，因此沒(méi)有自旋污染。 一致性的總能量使得這種方法可以計(jì)算單態(tài)-三態(tài)能隙：對(duì)單重態(tài)用RHF，對(duì)三重態(tài)用半電子近似。半電子近似計(jì)算不服從Koopman定理。也不能得到自旋密度。Mulliken布居分析通常是相當(dāng)合理的。 補(bǔ)充信息 一些討論和結(jié)果在： W. J. Hehre, L. Radom, P. v.R. Schleyer, J. A. Pople "Ab Initio Molecular Orbital Theory" Wiley (1986) An article that compares unrestricted, restricted and projected results is M. W. Wong, L. Radom J. Phys. Chem. 99, 8582 (1995) Some specific examples and a discussion of the half electron method are given in T. Clark "A Handbook of Computational Chemistry" Wiley (1985) A more mathematical treatment can be found in the paper J. S. Andrews, D. Jayatilake, R. G. A. Bone, N. C. Handy, R. D. Amos Chem. Phys. Lett. 183, 423 (1991) SUHF results are examined in P. K. Nandi, T. Kar, A. B. Sannigrahi Journal of Molecular Structure (Theochem) 362, 69 (1996) P.S 使用 CASSCF，CASPT2，MRMP2，MRCI等multi-reference方法可以比較有效的消除spin contamination的問(wèn)題。根據(jù)個(gè)人體會(huì)使用PUMP2(projected UMP2)似乎也是不錯(cuò)的選擇，而且速度比基于multi-reference的方法更快（CASSCF除外）。ab initio的計(jì)算主要基于error cancelling 比如反應(yīng)熱，反應(yīng)勢(shì)壘等，沒(méi)有人對(duì)絕對(duì)能量感興趣。anyway，雖然UHF有spin contamination的問(wèn)題，但是使用UHF計(jì)算的ground state的能量一般比ROHF更低，因?yàn)榍罢咧蠰CAO的自由度更多。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化計(jì)算及模擬 | 寶貝 |
» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有178人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
2樓2008-04-05 05:24:16
3樓2008-04-05 10:07:07
4樓2008-04-05 13:45:33
5樓2008-04-12 11:15:02
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zzgyb 的主題更新
| 11 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無(wú)極 2026-03-20 | 5/250 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|