| 24小時(shí)熱門(mén)版塊排行榜 |
| 6 | 1/1 | 返回列表 |
| 查看: 2755 | 回復(fù): 5 | |||
伊點(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ě)手)
|
可以嘗試使用 Occupation Thermal 0.005 , 如果計(jì)算能夠收斂的話(huà),再把此值減小,直到,采用feimi占據(jù),因?yàn),要?jì)算體系磁性,盡量使用Fermi占據(jù),不然的話(huà),會(huì)導(dǎo)致分?jǐn)?shù)占據(jù)的。 |

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

金蟲(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ì)稱(chēng)性,每個(gè)核實(shí)際上是帶1/2電子的,那么這時(shí)的能量如何計(jì)算呢? 利用量子統(tǒng)計(jì)的方法,Perdew等人證明,對(duì)于電子數(shù)N+delta (0 因此,對(duì)于精確的密度泛函理論,或者說(shuō),精確的交換相關(guān)泛函,能量E對(duì)電子數(shù)N的曲線(xiàn)應(yīng)為下圖1的紅色線(xiàn)。很明顯,在整數(shù)電子處,能量對(duì)電子數(shù)的一階微分的不連續(xù)的(左,右導(dǎo)數(shù)不相等)!這是一個(gè)很奇異的性質(zhì)。也就是說(shuō),下面這個(gè)關(guān)系成立: 并且,由于線(xià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)例外。 那么近似的泛函滿(mǎn)足這個(gè)條件嗎?考慮一個(gè)模型,假設(shè)H2電子數(shù)為n(不一定是整數(shù)N),則相應(yīng)的電子密度為: 用這個(gè)密度做X-alpha計(jì)算,可以得到n belongsto [0,2]的能量曲線(xiàn),即圖1中藍(lán)色線(xiàn): 圖1 很明顯,Xα泛函失敗了,在N=1處,它的一階微分是連續(xù)的,不存在斜率的跳躍。可以證明,對(duì)于包括B3LYP在內(nèi)的所有現(xiàn)在的泛函,沒(méi)有一個(gè)泛函能滿(mǎn)足不連續(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) (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +9 | 小熊joy 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +9 | 冬十三 2026-03-15 | 10/500 |
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯! 2026-03-15 | 14/700 |
|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
[考研] 材料專(zhuān)碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 機(jī)械專(zhuān)碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專(zhuān)家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂(lè)的香蕉 2026-03-11 | 4/200 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|