含钴(Co)体系在Gaussian的过渡态计算求助

Solitude198

使用的条件   #p opt=(calcall,ts,noeigen) freq genecp temperature=313.15 ub3lyp guess=mix em=gd3bj ,基组是6-31g(d,p) Co用lanl2dz
在尝试不同自旋多重度基本确定5（能量最低），然而在用0  1，0  3相继算出过渡态后，使用0    5计算总是算不出来，于是退而求其次不用D3校正，居然成功算出，于是想读取它的chk文件继续用原来的算，但还是没有结果。
目前还尝试了 opt=(calcall,ts,noeigen,maxstep=5,notrust) opt=(calcall,ts,noeigen,maxstep=5,notrust,gdiis) opt=(recalc=3,ts,noeigen,maxstep=5,notrust,gdiis)等，都无法收敛


目前在用PBE0-D3(BJ)尝试（似乎早就应该这么干），但不知道结果如何，求各位老师赐教。

zjxitcc

针对你说的“算不出来”、“都无法收敛”，如果是SCF不收敛，看这篇《解决SCF不收敛问题的方法》http://sobereva.com/61
如果是几何结构优化不收敛，可以看这篇《量子化学计算中帮助几何优化收敛的常用方法》http://sobereva.com/164
你的体系无论是SCF还是几何结构优化 都必可收敛。calcall太贵了，没有必要。在2024年来看，LANL2DZ可能精度太低了，有很多好得多的选择，例如LANL08(f)或def2TZVP，Co原子只有1个，可以用好一点的基组，不要过于小气而丢失了起码的精度。

Solitude198

zjxitcc 发表于 2024-10-8 11:31
针对你说的“算不出来”、“都无法收敛”，如果是SCF不收敛，看这篇《解决SCF不收敛问题的方法》http://sob ...

非常感谢 ORz
不好意思，上文的calcall应该都改成calcfc，我在这里打错了抱歉

sob老师的几篇文章都有看，所以尝试了许多不同方法

基组确实有点老套，我再尝试尝试