| 24小時(shí)熱門(mén)版塊排行榜 |
| 6 | 1/1 | 返回列表 |
| 查看: 2753 | 回復(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ù)的,不存在斜率的跳躍?梢宰C明,對(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ā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +3 | 煙汐憶海 2026-03-16 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +6 | 呼呼呼。。。。 2026-03-17 | 6/300 |
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 求調(diào)劑,總分315,考的生物醫(yī)藥,一志愿湖南師范大學(xué)。調(diào)劑到任何專(zhuān)業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 中科大材料專(zhuān)碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|