求助：模仿Molecular Tailoring Approach方法估计分步开环反应的键解离能

815268308

博文http://sobereva.com/522中提到对A、B、C、D四种状态都做了单点计算（不要分别做优化）。其中A的结构就是对本文体系优化后的结构，B、C、D是把A的结构的基团在GaussView里替换为氢的情况，键轴方向和原先一致，替换后C-H键长度是GaussView默认的1.07埃。

以下是学生估计的分布开环反应的方法：
-E1+EH1=B-A，因此EH1=B-1+E1
-E2+EH2=C-A,因此EH2=C-A+E2,E1、E2是开环反应的能量
E3=-EH1-EH2,EH1、EH2是加氢形成饱和物的能量,E3是脱氢形成两个自由基的能量
由于E=F+G-A(E是最终需要得到的键解离能),F+G=D+E+E3,推导出E=（A-B-C+D+E）/2
问题：
1、请问各位老师这个推导是否合理？
2、我估计的是键解离能，是否需要对先ABCDE做freq（A结构是优化后的。方法等同之前的计算 M062X/6-311g++ d,p），再通过更高精度（def2-TZVP）计算单点能，用焓（焓矫正量+单点能）来表示上述计算中的能量？
3、由于双环结构再GaussView中直接将A结构的基团加氢会改变原本存在的C-H键的键轴方向，并且开环处的C-H键有重叠的现象，学生手动调整H原子的位置是否会影响结果？

sobereva

这种体系直接用拉普拉斯键级估计键能就完了，省事得多，下文讲拉普拉斯键级部分说了
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471（http://bbs.keinsci.com/thread-12488-1-1.html）

815268308

sobereva 发表于 2023-7-30 01:07
这种体系直接用拉普拉斯键级估计键能就完了，省事得多，下文讲拉普拉斯键级部分说了
Multiwfn支持的分析化 ...

感谢sob老师的建议，这就去学习