使用shermo结合表格计算甲烷氢抽取反应的吉布斯自由能变化

糖糖糖豆9988 · 发表于 Post on 2021-12-3 19:37:55

本帖最后由糖糖糖豆9988 于 2022-10-25 19:14 编辑

刚开始接触计算，想从简单的反应开始再结合老师的表格尝试计算速率常数。使用高斯得到的甲烷氢抽取反应的过渡态，结合thermo得到过渡态的G（T），甲烷和H原子的G同理，计算的甲烷氢抽取反应的吉布斯自由能变化和帖子里发布的kisthelp计算出来差异很大？大家可以帮助找找原因吗？使用的理论和基组：b3lyp/6-311g(d,p)，甲烷：-40.5069616 hartree氢原子：-0.5128102 hartree，过渡态=G(298.15 K)= -40.9988028 hartree，△G=55.1541095KJ/mol，例子里面是自由能变化是80左右。请大家帮帮忙。

wzkchem5 · 发表于 Post on 2021-12-3 19:54:08

检验了过渡态的虚频吗？有没有做过IRC？

糖糖糖豆9988 · 发表于 Post on 2021-12-3 20:21:17

wzkchem5 发表于 2021-12-3 19:54
检验了过渡态的虚频吗？有没有做过IRC？

老师您好，做了IRC，虚频是-1108.95cm^-1.

糖糖糖豆9988 · 发表于 Post on 2021-12-3 20:23:57

这是做的irc。

糖糖糖豆9988 · 发表于 Post on 2021-12-3 20:54:51

本帖最后由糖糖糖豆9988 于 2021-12-4 08:35 编辑

我换成了例子中使用的M602x泛函，采用相同的基组，频率校正因子也一样0.97，高斯计算的虚频首先就变大了，是1500多。这回通过shermo计算的自由能变是77.6026912，与80接近了。但是我发现不论是高斯还是shermo都和例子里面kisthelp计算的G（一系列热力学数据）有差异。然后借助表格计算的化学反应速率常数是2.193E-20,例子里面是2.1713，这样应该还可以吧？另外请问表格里的旋转对称数如果是利用shermo计算的，就默认值1，对吗？

sobereva · 发表于 Post on 2021-12-4 01:55:42

哪来的什么thermo
不知名程序就给出链接，知名程序就把名字拼对了

糖糖糖豆9988 · 发表于 Post on 2021-12-4 08:37:01

sobereva 发表于 2021-12-4 01:55
哪来的什么thermo
不知名程序就给出链接，知名程序就把名字拼对了

很抱歉老师，shermo，非常抱歉使用人家的成果还给说错了。

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

[新手求助] 使用shermo结合表格计算甲烷氢抽取反应的吉布斯自由能变化