过渡态能垒为负值怎么解决

vixtrunt · 发表于 Post on 2026-5-15 15:46:38

老师们好，本人最近在计算一个金属催化氢转移的体系，算出过渡态之后，跑完IRC之后，选取第一帧和最后一帧优化出了前体和后体，发现前体自由能比过渡态自由能低了20kcal/mol，然后m062x方法重新算了一遍，发现能垒变成-30kcal/mol，且EE值也是-27kcal/mol，这种情况怎么解决。结构优化和过渡态优化关键词分别是
#p opt=noeigen freq m062x/6-31g(d) scrf=(iefpcm,solvent=acetone) 5d em=gd3
#p opt=(calcfc,ts,noeigen,maxstep=3,notrust) freq m062x/6-31g(d,p) scrf=(iefpcm,solvent=acetone) 5d em=gd3

非常非常感谢！

wzkchem5 · 发表于 Post on 2026-5-15 16:25:57

基组必须统一

睡觉不打鼾 · 发表于 Post on 2026-5-15 18:42:35

这不好说。个人经验，多手调，尝试各种可能，只靠程序很容易就从极小值跳出来了。

cym · 发表于 Post on 2026-5-16 09:29:12

二楼说的有道理。基组得统一。
一个是6-31g(d)，一个是6-31g(d,p) ，这肯定是不行的。

J.Liu · 发表于 Post on 2026-5-16 10:35:48

计算水平要统一

阿司匹林泡腾片 · 发表于 Post on 2026-5-18 08:31:02

之前用vasp算也遇到过类似问题，我的经验是前后统一多做调试

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

[Gaussian/gview] 过渡态能垒为负值怎么解决

浏览过的版块