 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3656 | 回復(fù): 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īng)該是沒有收斂，但是應(yīng)該如何優(yōu)化或者修改？救救孩子吧。。。。@jiewei |
回復(fù)此樓» 猜你喜歡
cif求助下載
已經(jīng)有1人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有5人回復(fù)
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有58人回復(fù)
-
基金正文30頁指的是報告正文還是整個申請書
已經(jīng)有9人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有2人回復(fù)
-
化學化工學院 招收化工、化學、材料等相關(guān)方向研究生（學碩、專碩都有調(diào)劑名額）
已經(jīng)有1人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
溫州大學新能源材料微觀溯源課題組招收2026級博士研究生（材料、化學方向）
已經(jīng)有6人回復(fù)
paramecium86
版主 (著名寫手)
- QC強帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26380.1
- 紅花: 349
- 帖子: 2025
- 在線: 1442小時
- 蟲號: 2466443
- 注冊: 2013-05-16
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 量子化學
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
導致scf不收斂的原因很多。通常默認的128步是夠用的。如果不收斂。可能有其他一些原因。直接增加步數(shù)用處可能會很小。有一些比較通用的做法可以看我這個教程里第2.2節(jié)scf不收斂的應(yīng)對方法。 http://m.gaoyang168.com/t-14242735-1 因為沒法發(fā)截圖了 那就把內(nèi)容復(fù)制到這里，不過排版就沒有那么好了。 如下： a. 首先檢查輸入結(jié)構(gòu)是不是合理。絕大部分情況下都是因為結(jié)構(gòu)跟 “真實” 結(jié)構(gòu)差的很遠，導致的scf難以收斂。盡量給一個能給出的最好結(jié)構(gòu) （比如先用uff粗略優(yōu)化一下，再用pm6 在優(yōu)化一下。然后再用dft計算之類的） b. 檢查體系的自旋多重度是否設(shè)置的合理。尤其是過渡金屬體系。如果多重度設(shè)置的偏離 “真實” 情況很遠，則非常難以收斂 c. 檢查基組是否合理。尤其是包含重原子，是否合理的設(shè)置了混合基組之類的都會影響。具體參考本文1.1節(jié)設(shè)置混合基組和贗勢的方法。這里你可以多看看 我看你用的lanl2dz如果你這個體系不是一個重金屬的團簇 對整個體系都用這個贗勢基組未必是好選擇 d. 檢查當前給的初始結(jié)構(gòu)是否真實存在于 “氣相” 中。因為有時候氣相和溶液相結(jié)構(gòu)很不一樣。這就會導致初始結(jié)構(gòu)以為很準其實是很差的情況。如果是溶液中的結(jié)構(gòu)則需要加上溶劑模型scrf，見本文1.2 .有時候加了溶劑模型就可以幫助收斂。 e. 以上確定是沒有問題的，這時候就考慮加入一些關(guān)鍵詞幫助收斂 scf=(xqc,maxconventionalcycles=64) 這個關(guān)鍵詞可以應(yīng)用在絕大部分情形下都有效。但是絕不要一上來就用，一定要檢查之前的步驟一都沒問題了再去使用。g09的話就用scf=xqc就好。 g16可以寫這一串。 scf=novaracc 高斯在計算scf過程中，最開始的幾步為了加速會降低積分精度。這樣有時候就會讓波函數(shù)不小心進入了一個錯誤的區(qū)間。很難再優(yōu)化下去。所以可以用這個關(guān)鍵詞，取消這種加速降低積分精度的方法。 scf(vshift=500) 這個一般可以用在過渡金屬體系。這種體系，可能homo-lumo間距過小導致波函數(shù)不容易收斂�？梢杂嬎阒腥藶樘哌@個值。這個不會影響結(jié)果只會影響收斂過程。一般設(shè)定在500上下比較合理。 int(acc2e=12, ultrafine) 如果用g09進行計算，也可以考慮這個關(guān)鍵詞提高了dft的積分精度和格點精度。這個僅僅適用在某些dft計算中。比如m062x 這類對于格點要求比較高的泛函中。g16中這個已經(jīng)是默認值不需要寫。 實際計算可以用以上提到的一個或者幾個關(guān)鍵詞一起幫助收斂。 如果體系實在是太過巨大，這些都做了還是很難收斂�？梢允褂� scf(conver=6) 這個關(guān)鍵詞降低scf收斂限到10^-6， 而默認是10^-8。但是慎用。這樣的結(jié)構(gòu)可能得到的根本不準。最好是優(yōu)化完成后，讀取波函數(shù)和結(jié)構(gòu)恢復(fù)到默認的收斂限再進行一次優(yōu)化。 |
» 本帖已獲得的紅花（最新10朵）
5樓2021-06-22 05:39:13
DX0929
至尊木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 17385.4
- 散金: 100
- 紅花: 2
- 帖子: 2223
- 在線: 339.6小時
- 蟲號: 1644074
- 注冊: 2012-02-26
- 專業(yè): 金屬有機化學
2樓2021-06-21 16:38:33
sky423
新蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1595.1
- 散金: 352
- 紅花: 8
- 帖子: 1048
- 在線: 103.6小時
- 蟲號: 5800207
- 注冊: 2017-03-01
- 專業(yè): 有機合成
送紅花一朵 |
#p b3lyp/lanl2dz opt pop=(full,npa) SCF=(maxcyc=200) SCF=QC 類似于這樣嗎？謝謝，我太菜了，謝謝您 發(fā)自小木蟲Android客戶端 |
3樓2021-06-21 17:15:48
paramecium86
版主 (著名寫手)
- QC強帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26380.1
- 紅花: 349
- 帖子: 2025
- 在線: 1442小時
- 蟲號: 2466443
- 注冊: 2013-05-16
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 量子化學
4樓2021-06-22 05:30:15
sky423
新蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1595.1
- 散金: 352
- 紅花: 8
- 帖子: 1048
- 在線: 103.6小時
- 蟲號: 5800207
- 注冊: 2017-03-01
- 專業(yè): 有機合成
6樓2021-06-22 18:38:48
sky423
新蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1595.1
- 散金: 352
- 紅花: 8
- 帖子: 1048
- 在線: 103.6小時
- 蟲號: 5800207
- 注冊: 2017-03-01
- 專業(yè): 有機合成
7樓2021-06-22 18:39:53
sky423
新蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1595.1
- 散金: 352
- 紅花: 8
- 帖子: 1048
- 在線: 103.6小時
- 蟲號: 5800207
- 注冊: 2017-03-01
- 專業(yè): 有機合成
8樓2021-06-22 18:40:46
sky423
新蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1595.1
- 散金: 352
- 紅花: 8
- 帖子: 1048
- 在線: 103.6小時
- 蟲號: 5800207
- 注冊: 2017-03-01
- 專業(yè): 有機合成
|
仔細檢查完之后發(fā)現(xiàn)結(jié)果的確有問題，所以重新構(gòu)建了分子，但是新的分子卻壓根無法導入gaussian進行計算，報錯，如下圖分別是03和09的報錯 發(fā)自小木蟲Android客戶端 |
9樓2021-06-22 21:33:03
paramecium86
版主 (著名寫手)
- QC強帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26380.1
- 紅花: 349
- 帖子: 2025
- 在線: 1442小時
- 蟲號: 2466443
- 注冊: 2013-05-16
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 量子化學
10樓2021-06-23 01:09:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 sky423 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿江南大學085701環(huán)境工程專碩總分287求調(diào)劑 +5 | 18266118446 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿河海大學085900土木水利專碩279求調(diào)劑不挑專業(yè) +3 | SunWwWwWw 2026-03-10 | 3/150 |
|
|
|
[考研] 材料工程307，求調(diào)劑 +7 | 我要燃燒你的夢 2026-03-08 | 7/350 |
|
|
|
[考研] 312求調(diào)劑 +5 | 陌宸希 2026-03-10 | 5/250 |
|
|
|
[考研] 標題:撿漏預(yù)警|08工科/09農(nóng)學調(diào)劑!英語要求低，過線即有機會! +7 | 馬超放煙花 2026-03-07 | 12/600 |
|
|
|
[考研] 298求調(diào)劑 +3 | Vv呀！ 2026-03-10 | 3/150 |
|
|
|
[考研] 313分生物學求調(diào)劑 +5 | Yyt楊1 2026-03-09 | 7/350 |
|
|
|
[考研] 0856材料與化工353分求調(diào)劑 +11 | NIFFFfff 2026-03-09 | 11/550 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑 +11 | 與冬清寧 2026-03-07 | 14/700 |
|
|
|
[考研] 材料與化工，291，求調(diào)劑 +12 | 咕嚕咕嚕123123 2026-03-05 | 13/650 |
|
|
|
[碩博家園] 木蟲好像不熱鬧了，是不是？ +4 | 偏振片 2026-03-10 | 4/200 |
|
|
|
[考研] 310 070300化學求調(diào)劑 +4 | 撲風鈴的貓 2026-03-08 | 5/250 |
|
|
|
[考研] 320求調(diào)劑 +4 | 魏zy 2026-03-08 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
|
[考研] 不限學校專業(yè)的調(diào)劑同學看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
|
[考研] 334求調(diào)劑 +6 | Trying] 2026-03-05 | 8/400 |
|