求助无能垒反应过程的能量及结构优化问题

qiubiye

【求助】各位老师好，我最近在做过渡态的尝试，因为刚接触gaussian，问题可能有些初级，希望各位老师不要介意，望各位老师赐教。
1、在反应过程中会有无能垒过程，比如C3H8——CH3.+C2H5.
这个反应过程的能量是用(CH3.单独做结构优化后的能量)+(C2H5.单独做结构优化的能量)-(C3H8单独做结构优化的能量）得到呢？
还是用(CH3.+C2H5放在一个gjf文件里做结构优化的能量）-(C3H8单独做结构优化的能量）得到呢？
我用前者方法得到的能量是92Kcal/mol,文献上计算的83kcal/mol；我用后者方法没做成功，在把CH3.+C2H5.放在一起优化后会得到与丙烷类似的结构，而不是想要的两个自由基的结构，可见图：图左为输入的文件gjf文件,是将C3H8copy过去后断键，再把CH3.和C2H5.拉开较长距离（比C3H8的C-C键长），图右为opt后的log文件。
2、第二个问题和也与上图相关，假如上图的反应有过渡态，在过渡态搜索前那肯定是需要把两个自由基放在一个gjf文件里作为生成物做优化的，那优化遇到上图情况有什么好的解决方法吗？

提前感谢各位老师啦！

qiubiye

granvia 发表于 2019-5-13 17:21
1.第一种方法就行。2. 没懂。既然自由基反应无能垒，干嘛找过渡态？

感谢前辈的帮助，第二个问题可以换种说法，图中的优化过后成为丙烷的这个现象是正常的吗，还是说可以通过一些方法比如说调整结构改变角度等等使其优化成CH3.+C2H5.两个片段，我是想请教下以后在做其他结构优化过程中遇到类似情况能否有较好的解决方法？

qiubiye

sobereva 发表于 2019-5-14 03:10
1
用(CH3.单独做结构优化后的能量)+(C2H5.单独做结构优化的能量)-(C3H8单独做结构优化的能量）

非常感谢sob老师的帮助，我又尝试算了一遍，发现是忘记了做零点能校正，加上之后就正确了。

sobereva

1
用(CH3.单独做结构优化后的能量)+(C2H5.单独做结构优化的能量)-(C3H8单独做结构优化的能量）

计算级别对结果影响极大

2 没过渡态

granvia

1.第一种方法就行。2. 没懂。既然自由基反应无能垒，干嘛找过渡态？