请问这样进行双自由基优化+能量计算的步骤是否正确？

林守龙

对于长链（22个原子)碳氢燃料分子解离反应，RH=R+H,对于C-H断键后的双自由基体系，进行优化和能量计算，进行了以下操作，不知道是否正确，麻烦老师们帮我看一下：
第一步：获得稳定的对称破缺波函数，为下一步结构优化提供初猜（结果是<S**2>=1.000）
%chk=R2.6.chk
# stable=opt um062x/6-311g(d,p) guess=mix

第二步：做限制性优化（一开始做了普通的优化opt,但是发现在优化后又变成未断键时的键长，请问是什么原因呢？），这时候是在上一步自旋极化波函数的基础上优化结构，得到体系的合理结构



%oldchk=R2.6.chk
# um062x/6-311g(d,p) guess=read opt=modredundant nosymm
坐标
18 20 F

第三步：取第二步优化后的结构，进行能量计算
%oldchk=R2.6.chk
# cbs-qb3

不知道这样的步骤是否正确，恳请各位老师帮忙解答。

sobereva

甭用CBS-QB3
否则又当闭壳层计算了

若断键结构下两个自旋相反单电子位点太近，即便初猜波函数是对称破缺单重态，优化的时候如果不加恰当的约束，也会自发成键

林守龙

sobereva 发表于 2024-4-19 02:58
甭用CBS-QB3
否则又当闭壳层计算了

明白了，谢谢sob老师