计算化学公社

标题: 含钴(Co)体系在Gaussian的过渡态计算求助 [打印本页]

作者
Author: Solitude198 时间: 2024-10-8 11:24
标题: 含钴(Co)体系在Gaussian的过渡态计算求助
使用的条件 #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)尝试（似乎早就应该这么干），但不知道结果如何，求各位老师赐教。
(, 下载次数 Times of downloads: 14)

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

作者
Author: Solitude198 时间: 2024-10-8 11:39

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

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

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)