高斯双自由基计算的焓变计算

yezhonghua

需要通过计算自由基的焓变值从而计算键的解离能；
A进行均裂产生两个自由基或者环断链形成一个双自由基：
A---B*+C* or A---B**（*代表单自由基；**代表双自由基）（如文献所示方法）
通过计算B*+C*的焓变减去A的焓变（ZPE矫正不做）计算BDE，计算结果和文献对应；但是如果计算B**双自由基的结果和文献差距很大（见计算结果统计图）；
目前有几点疑问：
1、一般来说相同泛函和基组B3LYP/6-31G*计算的BDE理论上应该可以作为键能大小的对比依据（比如说通过分子内，不同键计算的BED的结果和实际情况应该有正相关关系）；
2、双自由基计算的结果高斯有时候会自动将两个自由基链接在一起；为什么会这样？
3、双自由基的计算采用（0,1）或者（0,3）计算应该选取能量最低结果作为参考；
4、双自由基的计算如附件M1-3-N-O.log

谢谢大师指教！

sobereva

1 这个级别太烂。强烈建议用M06-2X/def2TZVP计算单点，把B3LYP/6-31G*得到的焓的热校正量加到上面，来作为体系的焓

2 gview纯粹根据距离来判断是否成键，不用管

3 双自由基的正确计算方式看
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

如果计算方式和级别都合理了，结果还是和实验对不上，考虑是否对于环内断键的计算模型和实际没有可比性。

yezhonghua

谢谢sober老师指教，茅塞顿开。