求助能垒计算时能垒差过大

TanJing · 发表于 Post on yesterday 17:12

各位老师好！我在做反应能垒计算时遇到一个问题想请教大家：
我初始结构优化 + 频率计算（opt+freq）采用的是 6-31g * 基组，直接读取该输出文件的 Hartree 数值计算能垒，得到能垒过小，不可靠。

opt freq.zip (3.91 MB, 下载次数 Times of downloads: 4)
后续我改用 cc-pvdz 高级基组做单点能计算，算出来的能量差值为 115.54 a.u，和初始结果严重偏离，数值异常不合理。

sp.zip (44.05 KB, 下载次数 Times of downloads: 3)
想请教各位老师这种情况是什么原因导致的、该如何解决？我自己考虑了两种方案：
选用更高水平基组重新计算单点能，结合热力学矫正后再重新计算能垒；
放弃低级别基组结果，直接采用更高精度基组重新做 opt+freq 全优化，直接取用其输出能量进行能垒分析。
麻烦各位老师帮忙指点一下哪种方案更合适，或者有没有其他更好的处理方法，谢谢各位老师！

Stardust0831 · 发表于 Post on yesterday 23:29

1.用类似“stable=opt“之类的关键词，检查波函数稳定性。
2.cc-pvdz这种基组通常不是给B3LYP这种普通杂化泛函用的（通常是给双杂化泛函或者post-HF用的），考虑换类似def2-TZVP这种基组，高斯里关键词直接写”def2TZVP“即可

sobereva · 发表于 Post on 2 hour ago

产生115.54 a.u这种数量级的能量差，必然犯了低级错误，最典型的就是化学组成不对应、用的基组/赝势不对应、数据读取/处理错误

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 求助能垒计算时能垒差过大

浏览过的版块