| 5 | 1/1 | 返回列表 |
| 查看: 2760 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
伊點(diǎn)若天金蟲(chóng) (小有名氣)
|
[求助]
Dmol3中參數(shù)的設(shè)置問(wèn)題。
|
||
|
在對(duì)過(guò)渡金屬氧化物進(jìn)行結(jié)構(gòu)優(yōu)化時(shí),成功優(yōu)化了其中一個(gè)表面,但是對(duì)于另外一個(gè)表面,總是出現(xiàn)以下錯(cuò)誤: 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 選項(xiàng)也有勾選,Occupation Thermal值設(shè)置為0.005,use symmetry也有勾選,請(qǐng)教下Spin值的調(diào)整問(wèn)題,怎樣才能進(jìn)行合理的參數(shù)設(shè)置來(lái)進(jìn)行優(yōu)化?或者是Spin值設(shè)定外的其他問(wèn)題。 (第一個(gè)表面spin設(shè)定為0.5成功優(yōu)化,第二個(gè)表面spin值也是0.5,出現(xiàn)上述問(wèn)題,但不確定是否為spin值的原因)。 |
金蟲(chóng) (正式寫(xiě)手)
|
密度泛函理論發(fā)展到現(xiàn)在,已經(jīng)能夠相當(dāng)好的解決很多實(shí)際的問(wèn)題。特別是“B3LYP”由于其在有機(jī)化學(xué)中的巨大成功幾乎已經(jīng)成為DFT的代名詞。DFT的核心就是“交換相關(guān)泛函”。如果這個(gè)東西的精確形式被確定,那么,也許從此以后量子化學(xué)家就會(huì)失業(yè)了。當(dāng)然,短時(shí)間內(nèi)這種可能性似乎不存在。這里通過(guò)一個(gè)有趣的現(xiàn)象:“分?jǐn)?shù)占據(jù)數(shù)”來(lái)討論一些交換相關(guān)泛函的一種系統(tǒng)誤差。 我們知道,基于波函數(shù)的從頭算方法是在3N維Hilbert空間中進(jìn)行的,而DFT則是在3維的Euclid(真實(shí))空間中進(jìn)行的,在真實(shí)空間中,我們可以考察密度的變化對(duì)能量的影響。對(duì)于真實(shí)的N電子體系,我們能夠確定其能量。那么,如果是“分?jǐn)?shù)”電子呢?比如,H2+在兩個(gè)H無(wú)窮遠(yuǎn)的時(shí)候,考慮對(duì)稱性,每個(gè)核實(shí)際上是帶1/2電子的,那么這時(shí)的能量如何計(jì)算呢? 利用量子統(tǒng)計(jì)的方法,Perdew等人證明,對(duì)于電子數(shù)N+delta (0 因此,對(duì)于精確的密度泛函理論,或者說(shuō),精確的交換相關(guān)泛函,能量E對(duì)電子數(shù)N的曲線應(yīng)為下圖1的紅色線。很明顯,在整數(shù)電子處,能量對(duì)電子數(shù)的一階微分的不連續(xù)的(左,右導(dǎo)數(shù)不相等)!這是一個(gè)很奇異的性質(zhì)。也就是說(shuō),下面這個(gè)關(guān)系成立: 并且,由于線性插值性: 其中I和A分別為電離能和電子親和能。而且,我們發(fā)現(xiàn),E對(duì)N是凸的,用化學(xué)語(yǔ)言就是說(shuō),第一電離能小于第二電離能,第二電離能小于第三電離能,等等。這個(gè)結(jié)論到目前為止無(wú)法從數(shù)學(xué)上證明,但是在實(shí)驗(yàn)和數(shù)值計(jì)算上沒(méi)有發(fā)現(xiàn)例外。 那么近似的泛函滿足這個(gè)條件嗎?考慮一個(gè)模型,假設(shè)H2電子數(shù)為n(不一定是整數(shù)N),則相應(yīng)的電子密度為: 用這個(gè)密度做X-alpha計(jì)算,可以得到n belongsto [0,2]的能量曲線,即圖1中藍(lán)色線: 圖1 很明顯,Xα泛函失敗了,在N=1處,它的一階微分是連續(xù)的,不存在斜率的跳躍。可以證明,對(duì)于包括B3LYP在內(nèi)的所有現(xiàn)在的泛函,沒(méi)有一個(gè)泛函能滿足不連續(xù)微分條件!這意味著,對(duì)于H2+,有 這是不可能的。這個(gè)結(jié)果是不符合實(shí)際的,這是現(xiàn)代所有交換相關(guān)泛函的系統(tǒng)誤差之一:它們?nèi)狈υ谡麛?shù)電子數(shù)處的不可微性。這就是分?jǐn)?shù)占據(jù)數(shù)(fractional occupation)誤差。這使得DFT容易高估離域分子的穩(wěn)定性。TDDFT常常在計(jì)算大共軛體系時(shí)出現(xiàn)嚴(yán)重的偏差,原因可能也在于此。 http://prb.aps.org/abstract/PRB/v57/i3/p1427_1 |

金蟲(chóng) (正式寫(xiě)手)
|
可以嘗試使用 Occupation Thermal 0.005 , 如果計(jì)算能夠收斂的話,再把此值減小,直到,采用feimi占據(jù),因?yàn),要?jì)算體系磁性,盡量使用Fermi占據(jù),不然的話,會(huì)導(dǎo)致分?jǐn)?shù)占據(jù)的。 |

金蟲(chóng) (小有名氣)
金蟲(chóng) (正式寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工專(zhuān)碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|