请问计算能垒时，反应物应该用何种能量

littleear · 发表于 Post on 2024-5-7 16:03:16

各位老师同学们好，最近我在计算一个反应的能垒。我发现如果把过渡态的能量减去所有反应物单独的能量之和，得到的能垒非常大，大约有上千kcal/mol。但是如果通过IRC得到所有反应物在一起的初始体系，则计算出来的能垒较低，约几十kcal/mol。我想请问一下，计算所有反应物的能量时，到底应该是所有反应物单独计算出能量后再求和，还是应该放在同一个体系中计算能量呢？
有此一问的主要原因是，通过IRC得到所有反应物在一起的初始体系后，再用相同水平进行opt+freq时，很容易不收敛，尝试了很多方法也不奏效。但是如果单独对每个反应物进行opt+freq，就不会出现不收敛的问题。
谢谢各位老师同学们。

wzkchem5 · 发表于 Post on 2024-5-7 16:56:00

过渡态的能量减去所有反应物单独的能量之和不应该有这么大，应该是反应物算错了，例如少画了氢原子，电荷/自旋多重度设错了，波函数收敛到了错误的地方，等等。
先解决这个问题，然后再把反应物形成反应复合物也看作一步反应，套用http://sobereva.com/506的结论判断应该用什么能量减什么能量

dddeng · 发表于 Post on 2024-5-7 17:02:07

IRC算出来的相当于是反应前络合物和反应后络合物，使用独立的反应物和反应前络合物这两种处理方式，应该都不会导致能量有如此大的差异，需要检查输入文件。反应前络合不收敛有很多原因，可能收敛标准设置太严格，也可能结构不太对，需要调整一下结构， sobereva老师的博文有关于解决结构不收敛的办法，可以去查看一下。

littleear · 发表于 Post on 2024-5-7 20:46:43

wzkchem5 发表于 2024-5-7 16:56
过渡态的能量减去所有反应物单独的能量之和不应该有这么大，应该是反应物算错了，例如少画了氢原子，电荷/ ...

收到，感谢老师，我去看看

littleear · 发表于 Post on 2024-5-7 20:47:10

dddeng 发表于 2024-5-7 17:02
IRC算出来的相当于是反应前络合物和反应后络合物，使用独立的反应物和反应前络合物这两种处理方式，应该都 ...

好滴好滴，非常感谢，我去看看

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

[Gaussian/gview] 请问计算能垒时，反应物应该用何种能量

评分 Rate