用Gaussian计算大碳片总是收敛不了？

986191088

我用b3lyp/def2-tzvp基组优化了一个带孔的大碳片，优化了好多天了，又接着在Linux系统上优化了两天，一直不能结束（但是能量不变化了），求助大神，这是什么情况？（gview是我中间的out文件）

End of G2Drv F.D. properties file   788 does not exist.
          IDoAtm=11111111111111111111111111111111111111111111111111
          IDoAtm=111111111111111111111
          Differentiating once with respect to electric field.
                with respect to dipole field.
          Differentiating once with respect to nuclear coordinates.
CalDSu exits because no D1Ps are significant.
          There are   216 degrees of freedom in the 1st order CPHF.  IDoFFX=6 NUNeed=     3.
    216 vectors produced by pass  0 Test12= 5.33D-13 1.00D-09 XBig12= 6.93D+03 3.11D+01.
AX will form    24 AO Fock derivatives at one time.
    216 vectors produced by pass  1 Test12= 5.33D-13 1.00D-09 XBig12= 1.07D+03 2.11D+00.

snljty

不要用Windows版，哪怕虚拟机/WSL里面装个Linux版，也能比Windows版快接近一倍。
用这么大基组肯定会相对慢一些，还是开壳层。
应该已经结束优化，开始算频率了。
输入文件建议直接贴文件不要贴图，别人想要测试一下还得一个字一个字敲进去。
最后这个截图左边纵坐标一个小格对应10 kJ/mol左右，不把右边几个点放大看不出来右边到底有多接近收敛。建议直接文本编辑器打开输出文件，搜最后一个Converged，然后看力和预测位移离收敛限还有多远。

不要把B3LYP-D3(BJ)和B3LYP说混。
这个体系很不适合用B3LYP。看卢老师博文《谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的》http://bbs.keinsci.com/thread-17899-1-1.html

986191088

snljty 发表于 2022-1-11 16:51
不要用Windows版，哪怕虚拟机/WSL里面装个Linux版，也能比Windows版快接近一倍。
用这么大基组肯定会相对 ...

多谢!

sobereva

你这个体系做opt freq用def2-TZVP严重不当
用6-311G*足矣，远远比带f极化函数的def2-TZVP便宜得多得多
注意看
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387（http://bbs.keinsci.com/thread-6600-1-1.html）

986191088

sobereva 发表于 2022-1-12 03:27
你这个体系做opt freq用def2-TZVP严重不当
用6-311G*足矣，远远比带f极化函数的def2-TZVP便宜得多得多
注 ...

多谢老师 我知道了