| 5 | 1/1 | 返回列表 |
| 查看: 3642 | 回復: 19 | ||
| 【懸賞金幣】回答本帖問題,作者sky423將贈送您 50 個金幣 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
sky423新蟲 (著名寫手)
|
[求助]
高斯不收斂,如何優(yōu)化 已有1人參與
|
|
|
#p b3lyp/lanl2dz opt pop=(full,npa) >>>>>>>>>> Convergence criterion not met. SCF Done: E(UB+HF-LYP) = -1714.45115901 A.U. after 129 cycles Convg = 0.8068D-05 -V/T = 2.0529 S**2 = 0.9172 KE= 1.628346096785D+03 PE=-1.275059909267D+04 EE= 5.005782636124D+03 Annihilation of the first spin contaminant: S**2 before annihilation 0.9172, after 0.7568 Convergence failure -- run terminated. Error termination via Lnk1e in d:program filesg03l502.exe at Fri Jun 18 18:59:12 2021. Job cpu time: 1 days 5 hours 55 minutes 14.0 seconds. File lengths (MBytes): RWF= 206 Int= 0 D2E= 0 Chk= 18 Scr= 1 請問高手,我這個從報錯看應該是沒有收斂,但是應該如何優(yōu)化或者修改?救救孩子吧。。。。@jiewei |
版主 (著名寫手)
至尊木蟲 (著名寫手)

新蟲 (著名寫手)
送紅花一朵 |
#p b3lyp/lanl2dz opt pop=(full,npa) SCF=(maxcyc=200) SCF=QC 類似于這樣嗎?謝謝,我太菜了,謝謝您 發(fā)自小木蟲Android客戶端 |
版主 (著名寫手)
|
導致scf不收斂的原因很多。通常默認的128步是夠用的。如果不收斂?赡苡衅渌恍┰。直接增加步數用處可能會很小。有一些比較通用的做法可以看我這個教程里第2.2節(jié)scf不收斂的應對方法。 http://m.gaoyang168.com/t-14242735-1 因為沒法發(fā)截圖了 那就把內容復制到這里,不過排版就沒有那么好了。 如下: a. 首先檢查輸入結構是不是合理。絕大部分情況下都是因為結構跟 “真實” 結構差的很遠,導致的scf難以收斂。盡量給一個能給出的最好結構 (比如先用uff粗略優(yōu)化一下,再用pm6 在優(yōu)化一下。然后再用dft計算之類的) b. 檢查體系的自旋多重度是否設置的合理。尤其是過渡金屬體系。如果多重度設置的偏離 “真實” 情況很遠,則非常難以收斂 c. 檢查基組是否合理。尤其是包含重原子,是否合理的設置了混合基組之類的都會影響。具體參考本文1.1節(jié)設置混合基組和贗勢的方法。這里你可以多看看 我看你用的lanl2dz如果你這個體系不是一個重金屬的團簇 對整個體系都用這個贗勢基組未必是好選擇 d. 檢查當前給的初始結構是否真實存在于 “氣相” 中。因為有時候氣相和溶液相結構很不一樣。這就會導致初始結構以為很準其實是很差的情況。如果是溶液中的結構則需要加上溶劑模型scrf,見本文1.2 .有時候加了溶劑模型就可以幫助收斂。 e. 以上確定是沒有問題的,這時候就考慮加入一些關鍵詞幫助收斂 scf=(xqc,maxconventionalcycles=64) 這個關鍵詞可以應用在絕大部分情形下都有效。但是絕不要一上來就用,一定要檢查之前的步驟一都沒問題了再去使用。g09的話就用scf=xqc就好。 g16可以寫這一串。 scf=novaracc 高斯在計算scf過程中,最開始的幾步為了加速會降低積分精度。這樣有時候就會讓波函數不小心進入了一個錯誤的區(qū)間。很難再優(yōu)化下去。所以可以用這個關鍵詞,取消這種加速降低積分精度的方法。 scf(vshift=500) 這個一般可以用在過渡金屬體系。這種體系,可能homo-lumo間距過小導致波函數不容易收斂?梢杂嬎阒腥藶樘哌@個值。這個不會影響結果只會影響收斂過程。一般設定在500上下比較合理。 int(acc2e=12, ultrafine) 如果用g09進行計算,也可以考慮這個關鍵詞提高了dft的積分精度和格點精度。這個僅僅適用在某些dft計算中。比如m062x 這類對于格點要求比較高的泛函中。g16中這個已經是默認值不需要寫。 實際計算可以用以上提到的一個或者幾個關鍵詞一起幫助收斂。 如果體系實在是太過巨大,這些都做了還是很難收斂?梢允褂 scf(conver=6) 這個關鍵詞降低scf收斂限到10^-6, 而默認是10^-8。但是慎用。這樣的結構可能得到的根本不準。最好是優(yōu)化完成后,讀取波函數和結構恢復到默認的收斂限再進行一次優(yōu)化。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學,英一數二305分求調劑,四六級已過 +5 | 小小番的茄 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿山東大學,總分327,英語二79,有論文,有競賽,已過四六級 +3 | 木木目目1 2026-03-09 | 3/150 |
|
|
[考研] 一志愿山東大學105500藥學專碩,總分302求調劑 +5 | 五維天空 2026-03-04 | 11/550 |
|
|
[考研] 材料與化工,291分,求調劑 +13 | 咕嚕咕嚕123123 2026-03-04 | 15/750 |
|
|
[考研] 一志愿武理314求調劑 +4 | ( ̄~ ̄;) 2026-03-08 | 5/250 |
|
|
[考研] 安徽農業(yè)大學材料與化學學院0856材料招收調劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
[考研] 085701環(huán)境工程專業(yè),初試305,均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
[考研] 081700學碩一志愿北京化工大學數二英一過六級有競賽求調劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
[考研] 材料化工求調劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 材料科學(0805)338 求調劑 +7 | xiaokang3286 2026-03-07 | 7/350 |
|
|
[考研] 化學290求調劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
[考研] 誠求調劑,323分 有專利/科研/實習經歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
[考研] 308求調劑 +7 | 倘若起風了呢 2026-03-05 | 9/450 |
|
|
[考研] 一志愿中石油,化學工程,301求調劑 +4 | 豚饅枳枳 2026-03-04 | 4/200 |
|
|
[考研] 085600,一志愿鄭州大學,280分求調劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
[論文投稿]
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
[考研] 331求調劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研]
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|
|
[考研] 化工專碩調劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|