| 5 | 1/1 | 返回列表 |
| 查看: 2756 | 回復(fù): 5 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
伊點若天金蟲 (小有名氣)
|
[求助]
Dmol3中參數(shù)的設(shè)置問題。
|
|
|
在對過渡金屬氧化物進行結(jié)構(gòu)優(yōu)化時,成功優(yōu)化了其中一個表面,但是對于另外一個表面,總是出現(xiàn)以下錯誤: Message: SCF not converging. Choose "Use Smearing" on DMol3 SCF panel or set "Occupation Thermal" in the input file You may also need to change spin or use symmetry Resubmit DMol3 Message: DMol3 job failed Error: DMol3 exiting Use Smearing 選項也有勾選,Occupation Thermal值設(shè)置為0.005,use symmetry也有勾選,請教下Spin值的調(diào)整問題,怎樣才能進行合理的參數(shù)設(shè)置來進行優(yōu)化?或者是Spin值設(shè)定外的其他問題。 (第一個表面spin設(shè)定為0.5成功優(yōu)化,第二個表面spin值也是0.5,出現(xiàn)上述問題,但不確定是否為spin值的原因)。 |
金蟲 (小有名氣)
金蟲 (正式寫手)

金蟲 (正式寫手)

金蟲 (正式寫手)
|
密度泛函理論發(fā)展到現(xiàn)在,已經(jīng)能夠相當好的解決很多實際的問題。特別是“B3LYP”由于其在有機化學(xué)中的巨大成功幾乎已經(jīng)成為DFT的代名詞。DFT的核心就是“交換相關(guān)泛函”。如果這個東西的精確形式被確定,那么,也許從此以后量子化學(xué)家就會失業(yè)了。當然,短時間內(nèi)這種可能性似乎不存在。這里通過一個有趣的現(xiàn)象:“分數(shù)占據(jù)數(shù)”來討論一些交換相關(guān)泛函的一種系統(tǒng)誤差。 我們知道,基于波函數(shù)的從頭算方法是在3N維Hilbert空間中進行的,而DFT則是在3維的Euclid(真實)空間中進行的,在真實空間中,我們可以考察密度的變化對能量的影響。對于真實的N電子體系,我們能夠確定其能量。那么,如果是“分數(shù)”電子呢?比如,H2+在兩個H無窮遠的時候,考慮對稱性,每個核實際上是帶1/2電子的,那么這時的能量如何計算呢? 利用量子統(tǒng)計的方法,Perdew等人證明,對于電子數(shù)N+delta (0 因此,對于精確的密度泛函理論,或者說,精確的交換相關(guān)泛函,能量E對電子數(shù)N的曲線應(yīng)為下圖1的紅色線。很明顯,在整數(shù)電子處,能量對電子數(shù)的一階微分的不連續(xù)的(左,右導(dǎo)數(shù)不相等)!這是一個很奇異的性質(zhì)。也就是說,下面這個關(guān)系成立: 并且,由于線性插值性: 其中I和A分別為電離能和電子親和能。而且,我們發(fā)現(xiàn),E對N是凸的,用化學(xué)語言就是說,第一電離能小于第二電離能,第二電離能小于第三電離能,等等。這個結(jié)論到目前為止無法從數(shù)學(xué)上證明,但是在實驗和數(shù)值計算上沒有發(fā)現(xiàn)例外。 那么近似的泛函滿足這個條件嗎?考慮一個模型,假設(shè)H2電子數(shù)為n(不一定是整數(shù)N),則相應(yīng)的電子密度為: 用這個密度做X-alpha計算,可以得到n belongsto [0,2]的能量曲線,即圖1中藍色線: 圖1 很明顯,Xα泛函失敗了,在N=1處,它的一階微分是連續(xù)的,不存在斜率的跳躍?梢宰C明,對于包括B3LYP在內(nèi)的所有現(xiàn)在的泛函,沒有一個泛函能滿足不連續(xù)微分條件!這意味著,對于H2+,有 這是不可能的。這個結(jié)果是不符合實際的,這是現(xiàn)代所有交換相關(guān)泛函的系統(tǒng)誤差之一:它們?nèi)狈υ谡麛?shù)電子數(shù)處的不可微性。這就是分數(shù)占據(jù)數(shù)(fractional occupation)誤差。這使得DFT容易高估離域分子的穩(wěn)定性。TDDFT常常在計算大共軛體系時出現(xiàn)嚴重的偏差,原因可能也在于此。 http://prb.aps.org/abstract/PRB/v57/i3/p1427_1 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 331求調(diào)劑(0703有機化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 0703,333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123! 2026-03-11 | 4/200 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|