计算化学公社

标题: 使用shermo结合表格计算甲烷氢抽取反应的吉布斯自由能变化 [打印本页]

作者
Author: 糖糖糖豆9988 时间: 2021-12-3 19:37
标题: 使用shermo结合表格计算甲烷氢抽取反应的吉布斯自由能变化
本帖最后由糖糖糖豆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左右。请大家帮帮忙。

作者
Author: wzkchem5 时间: 2021-12-3 19:54
检验了过渡态的虚频吗？有没有做过IRC？

作者
Author: 糖糖糖豆9988 时间: 2021-12-3 20:21

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

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

作者
Author: 糖糖糖豆9988 时间: 2021-12-3 20:23
这是做的irc。

作者
Author: 糖糖糖豆9988 时间: 2021-12-3 20:54
本帖最后由糖糖糖豆9988 于 2021-12-4 08:35 编辑

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

作者
Author: sobereva 时间: 2021-12-4 01:55
哪来的什么thermo
不知名程序就给出链接，知名程序就把名字拼对了

作者
Author: 糖糖糖豆9988 时间: 2021-12-4 08:37

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

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

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