算了3天,一個比較大的結(jié)構(gòu)(分子量2300)基于S0的DFT計算報錯,求助有什么辦法應(yīng)對?
gjf文件:
%chk=1.chk
%mem=10GB
%nproc=32
# B3LYP/genecp opt scf(maxcyc=200)
s0 opt
0 1
。。。
log文件最后一輪計算(未全部收斂):
Item Value Threshold Converged?
Maximum Force 0.023821 0.000450 NO
RMS Force 0.001642 0.000300 NO
Maximum Displacement 0.004522 0.001800 NO
RMS Displacement 0.000251 0.001200 YES
Predicted change in Energy=-5.391250D-04
。。。
log文件最后的報錯:
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GEDIIS/GDIIS optimizer.
Bend failed for angle 51 - 179 - 178
Tors failed for dihedral 50 - 51 - 179 - 178
Tors failed for dihedral 129 - 178 - 179 - 51
Tors failed for dihedral 180 - 178 - 179 - 51
Tors failed for dihedral 181 - 178 - 179 - 51
FormBX had a problem.
Error termination via Lnk1e in /opt/g09/l103.exe at Thu Jul 23 13:22:55 2020.
Job cpu time: 54 days 19 hours 25 minutes 7.5 seconds.
File lengths (MBytes): RWF= 3449 Int= 0 D2E= 0 Chk= 381 Scr= 1 @jiewei
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號
內(nèi)容已刪除
謝謝版主回復(fù)!十分感謝!
當前報錯的最后一步的結(jié)構(gòu)是log文件里的哪個部分。
我傳兩張圖,一張是開始寫的gjf,一張是log報錯前的原子表格(是這個表格嗎?如果是的話,是直接按列復(fù)制到gjf就行還是得把原子名稱寫出來)
gjf.jpg
log.jpg
,
你要不確定就用GaussView打開這個log文件。顯示的最后一幀就是最后算到的那個結(jié)構(gòu)。