| 8 | 1/1 | 返回列表 |
| 查看: 2753 | 回復(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é)果(當(dāng)然也有麻煩的,看運氣和經(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=當(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ā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +9 | zjl的號 2026-03-16 | 14/700 |
|
|---|---|---|---|---|
|
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 3/150 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯! 2026-03-15 | 14/700 |
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|