| 10 | 1/1 | 返回列表 |
| 查看: 2717 | 回復(fù): 9 | ||||
xiaoxiao136木蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】結(jié)構(gòu)優(yōu)化不收斂 已有4人參與
|
| 我用Gaussian優(yōu)化好的結(jié)構(gòu)用到adf,但是還是不收斂,該怎么調(diào)節(jié)呢,大家說(shuō)一說(shuō)啊,手冊(cè)上寫(xiě)的看得有點(diǎn)迷茫 |
ADF技術(shù) |
捐助貴賓 (正式寫(xiě)手)
|
以下內(nèi)容是我以前在一個(gè)壇子里看到的,抄到這里,你看看是否有用: 對(duì)于含有單電子的開(kāi)殼層體系,如過(guò)渡金屬或者稀土化合物,往往容易遇到SCF 不收斂的問(wèn)題。對(duì)于這類體系,ADF 軟件的使用手冊(cè)(ADFUsersGuide.pdf)針對(duì)SCF 不收斂問(wèn)題提供了以下幾種建議: 1, New DIIS 使用NEWDIIS 關(guān)鍵字,在SCF 過(guò)程中調(diào)用新的DIIS 方法,這可以解決一部分的不收斂問(wèn)題; 2, Keep Orbitals 對(duì)于某些SCF 特別難收斂的體系,如開(kāi)殼層高自旋體系可以使用關(guān)鍵字OCCUPATIONS 并指定keeporbitals,把軌道順序固定下來(lái)避免能量震蕩。這種方法幾乎一定能讓SCF 收斂,但是軌道的順序是否正確卻無(wú)法保證,除非事先指定軌道排布。圖形界面下ADFinput 可以先列出各個(gè)軌道并分別指定其中的電子數(shù),對(duì)于較小的分子體系不妨一試,利用這種方法甚至可以設(shè)置出一個(gè)激發(fā)態(tài)的電子軌道占據(jù)。 3, Mixing SCF 關(guān)鍵字中有個(gè)選項(xiàng)mix,其意義是根據(jù)占據(jù)軌道算的當(dāng)前勢(shì)的相對(duì)權(quán)重,默認(rèn)值為0.2。對(duì)于需要強(qiáng)阻尼的困難體系,適當(dāng)減少mix 的值有助于SCF 收斂,一般設(shè)置為0.02-0.05 之間比較合適。 4, Level Shift SCF 關(guān)鍵字中還有個(gè)選項(xiàng)能級(jí)移動(dòng)參數(shù)lShift,可以人為地將非占據(jù)軌道的能量拔高( 單位:hatree)。在SCF 迭代中,當(dāng)電荷在能量非常接近且都位于Fermi 能級(jí)附近的 不同軌道之間前后擺動(dòng)時(shí),這可幫助解決收斂問(wèn)題。不過(guò)需要注意,涉及到能級(jí)躍遷的計(jì)算( 如預(yù)測(cè)UV 譜) 使用此參數(shù)將導(dǎo)致結(jié)果出現(xiàn)偏差。 5, Smearing 對(duì)于非常難收斂的體系,還可以使用OCCUPATIONS關(guān)鍵字的smearq 選項(xiàng), 這相當(dāng)于給體系一個(gè)溫度并按照Boltzman 分布排電子, 電子排布不再是aufbau 的。Smearing 選項(xiàng)可以讓SCF 很快收斂,不過(guò)這將會(huì)引入分?jǐn)?shù)形式的軌道占據(jù)數(shù),對(duì)此可以先設(shè)定一個(gè)smearq 的值,然后逐步降低至零,類似于一個(gè)退火的過(guò)程。早期版本中這需要依次修改smearq 的大小然后多次restart, 現(xiàn)在直接在"smearq=" 之后填一串?dāng)?shù)字即可。ADFinput 中的設(shè)置更簡(jiǎn)單,分別指定smearq 的初始值和終值,再給出退火的步數(shù),它將自動(dòng)生成一系列遞減數(shù)值并且步長(zhǎng)也遞減,比手動(dòng)設(shè)置的會(huì)更精細(xì)一些。對(duì)于smearq 的每一個(gè)值,SCF 都將作一次自洽,因此SCF 的最大循環(huán)圈數(shù)需要放寬,一般按每個(gè)值50 圈估算。smearq 的初始值如何設(shè)置,這取決于具體的體系,需要嘗試。 |
木蟲(chóng) (正式寫(xiě)手)

木蟲(chóng) (正式寫(xiě)手)

木蟲(chóng) (正式寫(xiě)手)

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫(xiě)郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考研] 生物學(xué)調(diào)劑招人!! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 321求調(diào)劑 +8 | 何潤(rùn)采123 2026-03-18 | 10/500 |
|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|