| 8 | 1/1 | 返回列表 |
| 查看: 2748 | 回復(fù): 7 | |||||
亍寞金蟲 (小有名氣)
|
[求助]
opt=TS 循環(huán)步驟問題
|
|
命令這么寫的,但是走完第一圈就9999出錯了,請高人指教! %mem=15000MB %nprocshared=8 #p opt(maxcyc=300,TS,calcfc) freq b3lyp/6-311++g(2d,2p) scf(maxcyc=300,tight) 該怎么弄呢? |
科研資料 | computational chemistry (gaussian) | Gaussian |
榮譽版主 (著名寫手)
|
看到好幾個帖子說L999問題了,我把量化版的FAQ上面的貼上來 L9999錯誤是屬于構(gòu)型優(yōu)化沒有找到最后結(jié)果,即在規(guī)定的優(yōu)化步驟內(nèi)沒有找到最終的合適的結(jié)構(gòu)。(L9999錯誤)。對于優(yōu)化不收斂,即L9999錯誤,實際上是在規(guī)定的步數(shù)內(nèi)沒有完成優(yōu)化,即還沒有找到極小值點。(或者對于過渡態(tài)優(yōu)化,還沒有找到過渡態(tài)),可以增加優(yōu)化步驟,或是拿最后的結(jié)構(gòu) 接著進(jìn)行優(yōu)化。此外這有幾種可能性:1. 看一下能量的收斂的情況,可能正在單調(diào)減小,眼看有收斂的趨勢,這樣的情況下,只要加大循環(huán)的步數(shù)opt(maxcycle=200),可能就可以解決問題了。2. 加大循環(huán)步數(shù)還不能解決的(循環(huán)步數(shù)有人說超過200再不收斂,再加也不會有用了,這雖然不一定絕對正確,但200步應(yīng)該也差不多了),有兩種可能。一是查看能量,發(fā)現(xiàn)能量在振蕩了,且變化已經(jīng)很小了,這時可能重新算一下,或者構(gòu)型稍微變一下,繼續(xù)優(yōu)化,就可以得到收斂的結(jié)果(當(dāng)然也有麻煩的,看運氣和經(jīng)驗了);二是構(gòu)型變化太大,和你預(yù)計的差別過大,這很可能是你的初始構(gòu)型太差了,優(yōu)化不知道到哪里去了,這時最好檢查一下初始構(gòu)型,再從頭優(yōu)化。3. 對于L9999快達(dá)到收斂時,考慮減小優(yōu)化步長有時對于能量振蕩的情況也是有用的,opt(maxstep=1).一個建議是,對于大體系,難收斂體系,先用小基組,低精度算法優(yōu)化一下,以得到較好的初始構(gòu)型,再用高精度的計算接著算。如果前面的方法保留了chk文件,重新計算時需要使用 geom=allcheck 讀入構(gòu)型(就不必麻煩地寫構(gòu)型了), guess=read(讀入初始波函數(shù),可以加快第一步SCF收斂)。 |
鐵桿木蟲 (著名寫手)
金蟲 (小有名氣)
金蟲 (小有名氣)
L9999的錯誤的引發(fā)一般是因為Optimization stopped.這個可以從輸出文件中看到,這個產(chǎn)生的原因有很多,其中一個是-- Number of steps exceeded, NStep= N (N=當(dāng)前優(yōu)化的步數(shù)) 另一個是-- Wrong number of Negative eigenvalues: Desired= 1 Actual= N (N=當(dāng)前Hessian矩陣的本征值數(shù),找過渡態(tài)需要該值為1,通過添加noeigentest參數(shù)可以取消該檢查,雖然有時也能找到過渡態(tài),但一般情況下這說明過渡態(tài)的初始猜測并不是很合理)還有一個是 Optimization aborted. -- No acceptable step. 這個很明顯,構(gòu)象優(yōu)化出錯了,解決方法調(diào)構(gòu)象或者更改計算方法或條件…… 還有很多其他的原因,我這就不一一列舉了。所以解決l9999的錯誤最好的方法是在你輸出文件中仔細(xì)找一下到底是什么原因?qū)е掠嬎憬K止,然后對癥下藥 盲目的加maxcyc=N opt=loose opt=gdiis等等是沒什么意義的 |


金蟲 (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +3 | 簡單點0 2026-03-17 | 3/150 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|