| 8 | 1/1 | 返回列表 |
| 查看: 2757 | 回復(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) 接著進行優(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é)果(當然也有麻煩的,看運氣和經(jīng)驗了);二是構(gòu)型變化太大,和你預(yù)計的差別過大,這很可能是你的初始構(gòu)型太差了,優(yōu)化不知道到哪里去了,這時最好檢查一下初始構(gòu)型,再從頭優(yōu)化。3. 對于L9999快達到收斂時,考慮減小優(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=當前優(yōu)化的步數(shù)) 另一個是-- Wrong number of Negative eigenvalues: Desired= 1 Actual= N (N=當前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的錯誤最好的方法是在你輸出文件中仔細找一下到底是什么原因?qū)е掠嬎憬K止,然后對癥下藥 盲目的加maxcyc=N opt=loose opt=gdiis等等是沒什么意義的 |


金蟲 (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +3 | 想上岸的鯉魚 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[基金申請] 被我言中:新模板不強調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 302求調(diào)劑 +9 | 負心者當誅 2026-03-11 | 9/450 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|