| 2 | 1/1 | 返回列表 |
| 查看: 3784 | 回復(fù): 1 | ||
maoxiao木蟲 (正式寫手)
|
[求助]
求助大神關(guān)于SCF不收斂(Convergence failure -- run terminated)的問題 已有1人參與
|
|
在Gaussion 計(jì)算出現(xiàn)下面問題,求助大神怎么解決? Error on total polarization charges = 0.01305 SCF Done: E(RB3LYP) = -1557.98733960 A.U. after 129 cycles NFock=128 Conv=0.17D-04 -V/T= 2.0514 Convergence failure -- run terminated. Error termination via Lnk1e in D:\Gaussion09\G09W\l502.exe at Wed Sep 11 07:12:33 2024. Job cpu time: 0 days 7 hours 45 minutes 7.0 seconds. File lengths (MBytes): RWF= 499 Int= 0 D2E= 0 Chk= 28 Scr= 1 |
版主 (著名寫手)
|
可以參考我這個(gè)帖子里的教程 2.2節(jié) http://m.gaoyang168.com/t-14242735-1 ###################################################################################### A. 首先檢查輸入結(jié)構(gòu)是不是合理。絕大部分情況下都是因?yàn)榻Y(jié)構(gòu)跟 “真實(shí)” 結(jié)構(gòu)差的很遠(yuǎn),導(dǎo)致的SCF難以收斂。盡量給一個(gè)能給出的最好結(jié)構(gòu) (比如先用UFF粗略優(yōu)化一下,再用PM6 在優(yōu)化一下。然后再用DFT計(jì)算之類的) B. 檢查體系的自旋多重度是否設(shè)置的合理。尤其是過渡金屬體系。如果多重度設(shè)置的偏離 “真實(shí)” 情況很遠(yuǎn),則非常難以收斂 C. 檢查基組是否合理。尤其是包含重原子,是否合理的設(shè)置了混合基組之類的都會(huì)影響。具體參考本文1.1節(jié)設(shè)置混合基組和贗勢(shì)的方法。 D. 檢查當(dāng)前給的初始結(jié)構(gòu)是否真實(shí)存在于 “氣相” 中。因?yàn)橛袝r(shí)候氣相和溶液相結(jié)構(gòu)很不一樣。這就會(huì)導(dǎo)致初始結(jié)構(gòu)以為很準(zhǔn)其實(shí)是很差的情況。如果是溶液中的結(jié)構(gòu)則需要加上溶劑模型SCRF,見本文1.2 .有時(shí)候加了溶劑模型就可以幫助收斂。 E. 以上確定是沒有問題的,這時(shí)候就考慮加入一些關(guān)鍵詞幫助收斂 SCF=(XQC,MAXCONVENTIONALCYCLES=64) 這個(gè)關(guān)鍵詞可以應(yīng)用在絕大部分情形下。但是絕不要一上來就用,一定要檢查之前的步驟一都沒問題了再去使用。G09的話就用SCF=XQC就好。 G16可以寫這一串。 SCF=NoVarAcc 高斯在計(jì)算SCF過程中,最開始的幾步為了加速會(huì)降低積分精度。這樣有時(shí)候就會(huì)讓波函數(shù)不小心進(jìn)入了一個(gè)錯(cuò)誤的區(qū)間。很難再優(yōu)化下去。所以可以用這個(gè)關(guān)鍵詞,取消這種加速降低積分精度的方法。 SCF(VSHIFT=500) 這個(gè)一般可以用在過渡金屬體系。這種體系,可能HOMO-LUMO間距過小導(dǎo)致波函數(shù)不容易收斂。可以計(jì)算中人為太高這個(gè)值。這個(gè)不會(huì)影響結(jié)果只會(huì)影響收斂過程。一般設(shè)定在500上下比較合理。 INT(Acc2E=12, UltraFine) 如果用G09進(jìn)行計(jì)算,也可以考慮這個(gè)關(guān)鍵詞提高了DFT的積分精度和格點(diǎn)精度。這個(gè)僅僅適用在某些DFT計(jì)算中。比如M062X 這類對(duì)于格點(diǎn)要求比較高的泛函中。G16中這個(gè)已經(jīng)是默認(rèn)值不需要寫。 實(shí)際計(jì)算可以用以上提到的一個(gè)或者幾個(gè)關(guān)鍵詞一起幫助收斂。 如果體系實(shí)在是太過巨大,這些都做了還是很難收斂。可以使用 SCF(CONVER=6) 這個(gè)關(guān)鍵詞降低SCF收斂限到10-6, 而默認(rèn)是10-8。但是慎用。這樣的結(jié)構(gòu)可能得到的根本不準(zhǔn)。最好是優(yōu)化完成后,讀取波函數(shù)和結(jié)構(gòu)恢復(fù)到默認(rèn)的收斂限再進(jìn)行一次優(yōu)化。 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 售瀆SCI文章,我:8⊙ 55 1⊙ 54。備注【⊙=0】 +3 | vs9ky3szrj 2026-03-06 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +10 | ounce. 2026-03-04 | 12/600 |
|
|
[考研] 070300調(diào)劑一志愿南京大學(xué)化學(xué)專業(yè),初試273分,有論文有專利有競(jìng)賽,求調(diào)劑B區(qū)211 +4 | 82206202 2026-02-28 | 14/700 |
|
|
[考博] 售瀆SCI文章,我:8⊙ 55 1⊙ 54。備注【⊙=0】 +4 | ipys00zm58 2026-03-05 | 9/450 |
|
|
[基金申請(qǐng)] 進(jìn)入個(gè)人成果庫(kù)好難,一下午都沒進(jìn)去 +5 | mi_dilee 2026-03-05 | 5/250 |
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +4 | T123 tt 2026-03-04 | 4/200 |
|
|
[考研] 085600材料與化工(高分子)290分求調(diào)劑 +9 | wengyujian 2026-03-04 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼?~+123456 2026-03-05 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 4/200 |
|
|
[論文投稿]
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|
|
[基金申請(qǐng)] 沒有青基直接申請(qǐng)面上,感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
[考研] 306求調(diào)劑 +7 | 張張張張oo 2026-03-03 | 7/350 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[基金申請(qǐng)] 請(qǐng)問大家,研究風(fēng)險(xiǎn)與應(yīng)對(duì)措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 化工專碩342,一志愿大連理工大學(xué),求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|