算了3天,一个比较大的结构(分子量2300)基于S0的DFT计算报错,求助有什么办法应对?
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
返回小木虫查看更多
京公网安备 11010802022153号
内容已删除
谢谢版主回复!十分感谢!
当前报错的最后一步的结构是log文件里的哪个部分啊?
我传两张图,一张是开始写的gjf,一张是log报错前的原子表格(是这个表格吗?如果是的话,是直接按列复制到gjf就行还是得把原子名称写出来)
gjf.jpg
log.jpg
,
你要不确定就用GaussView打开这个log文件。显示的最后一帧就是最后算到的那个结构。