| 10 | 1/1 | 返回列表 |
| 查看: 2718 | 回復(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ì)于這類(lèi)體系,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 的值,然后逐步降低至零,類(lèi)似于一個(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ā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +18 | rare12345 2026-03-18 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 0817調(diào)劑 +3 | 沒(méi)有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 材料專(zhuān)碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽(yáng)CF 2026-03-18 | 4/200 |
|
|
[考研] 085410人工智能專(zhuān)碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 085601專(zhuān)碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 材料專(zhuān)碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|