 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3594 | 回復(fù): 19 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者sky423將贈(zèng)送您 50 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
sky423新蟲(chóng) (著名寫(xiě)手)
|
[求助]
高斯不收斂，如何優(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 請(qǐng)問(wèn)高手，我這個(gè)從報(bào)錯(cuò)看應(yīng)該是沒(méi)有收斂，但是應(yīng)該如何優(yōu)化或者修改？救救孩子吧。。。。@jiewei |
回復(fù)此樓» 猜你喜歡
福州大學(xué)梁宇航副教授招收2026年申核制博士研究生/碩士研究生（理論計(jì)算方向）
已經(jīng)有20人回復(fù)
-
深圳理工大學(xué)招聘AI4Science與分子模擬方向科研人員
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有237人回復(fù)
-
中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事
已經(jīng)有0人回復(fù)
-
全獎(jiǎng)博士 英國(guó)利物浦大學(xué) × 臺(tái)灣清華大學(xué) 聯(lián)合培養(yǎng)
已經(jīng)有0人回復(fù)
-
尋博導(dǎo)
已經(jīng)有13人回復(fù)
-
基金正文30頁(yè)指的是報(bào)告正文還是整個(gè)申請(qǐng)書(shū)
已經(jīng)有9人回復(fù)
paramecium86
版主 (著名寫(xiě)手)
- QC強(qiáng)帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26375.6
- 紅花: 349
- 帖子: 2025
- 在線: 1441.9小時(shí)
- 蟲(chóng)號(hào): 2466443
- 注冊(cè): 2013-05-16
- 性別: GG
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
導(dǎo)致scf不收斂的原因很多。通常默認(rèn)的128步是夠用的。如果不收斂�？赡苡衅渌恍┰颉Ｖ苯釉黾硬綌�(shù)用處可能會(huì)很小。有一些比較通用的做法可以看我這個(gè)教程里第2.2節(jié)scf不收斂的應(yīng)對(duì)方法。 http://m.gaoyang168.com/t-14242735-1 因?yàn)闆](méi)法發(fā)截圖了 那就把內(nèi)容復(fù)制到這里，不過(guò)排版就沒(méi)有那么好了。 如下： a. 首先檢查輸入結(jié)構(gòu)是不是合理。絕大部分情況下都是因?yàn)榻Y(jié)構(gòu)跟 “真實(shí)” 結(jié)構(gòu)差的很遠(yuǎn)，導(dǎo)致的scf難以收斂。盡量給一個(gè)能給出的最好結(jié)構(gòu) （比如先用uff粗略?xún)?yōu)化一下，再用pm6 在優(yōu)化一下。然后再用dft計(jì)算之類(lèi)的） b. 檢查體系的自旋多重度是否設(shè)置的合理。尤其是過(guò)渡金屬體系。如果多重度設(shè)置的偏離 “真實(shí)” 情況很遠(yuǎn)，則非常難以收斂 c. 檢查基組是否合理。尤其是包含重原子，是否合理的設(shè)置了混合基組之類(lèi)的都會(huì)影響。具體參考本文1.1節(jié)設(shè)置混合基組和贗勢(shì)的方法。這里你可以多看看 我看你用的lanl2dz如果你這個(gè)體系不是一個(gè)重金屬的團(tuán)簇 對(duì)整個(gè)體系都用這個(gè)贗勢(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，見(jiàn)本文1.2 .有時(shí)候加了溶劑模型就可以幫助收斂。 e. 以上確定是沒(méi)有問(wèn)題的，這時(shí)候就考慮加入一些關(guān)鍵詞幫助收斂 scf=(xqc,maxconventionalcycles=64) 這個(gè)關(guān)鍵詞可以應(yīng)用在絕大部分情形下都有效。但是絕不要一上來(lái)就用，一定要檢查之前的步驟一都沒(méi)問(wèn)題了再去使用。g09的話就用scf=xqc就好。 g16可以寫(xiě)這一串。 scf=novaracc 高斯在計(jì)算scf過(guò)程中，最開(kāi)始的幾步為了加速會(huì)降低積分精度。這樣有時(shí)候就會(huì)讓波函數(shù)不小心進(jìn)入了一個(gè)錯(cuò)誤的區(qū)間。很難再優(yōu)化下去。所以可以用這個(gè)關(guān)鍵詞，取消這種加速降低積分精度的方法。 scf(vshift=500) 這個(gè)一般可以用在過(guò)渡金屬體系。這種體系，可能homo-lumo間距過(guò)小導(dǎo)致波函數(shù)不容易收斂�？梢杂�(jì)算中人為太高這個(gè)值。這個(gè)不會(huì)影響結(jié)果只會(huì)影響收斂過(guò)程。一般設(shè)定在500上下比較合理。 int(acc2e=12, ultrafine) 如果用g09進(jìn)行計(jì)算，也可以考慮這個(gè)關(guān)鍵詞提高了dft的積分精度和格點(diǎn)精度。這個(gè)僅僅適用在某些dft計(jì)算中。比如m062x 這類(lèi)對(duì)于格點(diǎn)要求比較高的泛函中。g16中這個(gè)已經(jīng)是默認(rèn)值不需要寫(xiě)。 實(shí)際計(jì)算可以用以上提到的一個(gè)或者幾個(gè)關(guān)鍵詞一起幫助收斂。 如果體系實(shí)在是太過(guò)巨大，這些都做了還是很難收斂�？梢允褂� 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)化。 |
5樓2021-06-22 05:39:13
DX0929
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 17370.4
- 散金: 100
- 紅花: 2
- 帖子: 2220
- 在線: 339.6小時(shí)
- 蟲(chóng)號(hào): 1644074
- 注冊(cè): 2012-02-26
- 專(zhuān)業(yè): 金屬有機(jī)化學(xué)
2樓2021-06-21 16:38:33
sky423
新蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1595.1
- 散金: 352
- 紅花: 8
- 帖子: 1048
- 在線: 103.6小時(shí)
- 蟲(chóng)號(hào): 5800207
- 注冊(cè): 2017-03-01
- 專(zhuān)業(yè): 有機(jī)合成
送紅花一朵 |
#p b3lyp/lanl2dz opt pop=(full,npa) SCF=(maxcyc=200) SCF=QC 類(lèi)似于這樣嗎？謝謝，我太菜了，謝謝您 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
3樓2021-06-21 17:15:48
paramecium86
版主 (著名寫(xiě)手)
- QC強(qiáng)帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26375.6
- 紅花: 349
- 帖子: 2025
- 在線: 1441.9小時(shí)
- 蟲(chóng)號(hào): 2466443
- 注冊(cè): 2013-05-16
- 性別: GG
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
4樓2021-06-22 05:30:15
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工，270求調(diào)劑 +13 | YXCT 2026-03-01 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+5
|
15779376950 2026-03-01 | 7/350 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +7 | 土土小蟲(chóng) 2026-03-03 | 7/350 |
|
|
|
[考研] 一志愿西工大材料學(xué)，英一數(shù)二總分321分，求調(diào)劑。 +3 | zz05zz 2026-03-02 | 5/250 |
|
|
|
[考研] 理學(xué)，工學(xué)，農(nóng)學(xué)調(diào)劑，少走彎路，這里歡迎您！ +8 | likeihood 2026-03-02 | 11/550 |
|
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
|
[考研] 一志愿東北大學(xué)化學(xué)314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-02 | 6/300 |
|
|
|
[考研] 282求調(diào)劑 +4 | 2103240126 2026-03-02 | 7/350 |
|
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +8 | llx0610 2026-03-02 | 9/450 |
|
|
|
[考研] 085600材料工程一志愿中科大總分312求調(diào)劑 +9 | 吃宵夜1 2026-02-28 | 11/550 |
|
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 哈工大計(jì)算機(jī)劉劼團(tuán)隊(duì)招生 +4 | hit_aiot 2026-03-01 | 6/300 |
|
|
|
[考研] 322求調(diào)劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-02-28 | 5/250 |
|