关于双自由基算法求助

kevin

各位老师，
       在计算双自由基时，我有一些疑惑，我看过了sob老师的两篇博文《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82和《CASSCF计算双自由基以及双自由基特征的计算》http://sobereva.com/264。其中均以C4H8(.CH2-Ch2-CH2-CH2.)为例，我的理解这种双自由基两个单电子之间的作用较小，处于自旋极化单重态，最简单直接的算法是：
     #p  opt  ub3lyp/6-31g*  guess=mix
      
      0 1
      <分子坐标>
     我最近阅读的文献Energy Fuels 2011, 25, 3739–3746，也计算了类似的双自由基，结构上，两个单电子距离很远。不同之处是作者计算了单重态和三重态，结果显示三重态比单重态能量还低，有点疑惑。
我的问题是：
     1. 我对计算双自由基的理解是否正确?
     2. 文献中这种出现三重态比单重态能量更低是什么原因？
     3. 如果计算三重态的双自由基，计算文件是否如下
     #p  opt  ub3lyp/6-31g*  guess=mix
      
      0 3
      <分子坐标>

      
       谢谢各位老师！

sobereva

那篇文献我暂时没时间看，对于第三点，用不着写guess=mix，直接0 3就够了
你要是怀疑三重态能量更低，就三重态和自旋极化单重态都算算

yezhonghua

您好，我在计算的BDE（键解离能）的时候，环键断裂的情况产生的双自由基，opt freq ub3lyp/6-31G（d），0，3。计算的结果和文献差异较大。
请问您有接触过么？谢谢

kevin

yezhonghua 发表于 2016-12-14 11:30
您好，我在计算的BDE（键解离能）的时候，环键断裂的情况产生的双自由基，opt freq ub3lyp/6-31G（d），0， ...

按照我的理解，你这种关键词写法，计算的是三重态的，首先你得确定是三重态。如果是三重态的话，可能有其他的问题，这你需要再请教下sob老师。
如果不确定的话，可以试试单重态的 opt freq ub3lyp/6-31G(d) guess=mix  0 1看看和文献里是不是相同

yezhonghua

kevin 发表于 2016-12-14 16:04
按照我的理解，你这种关键词写法，计算的是三重态的，首先你得确定是三重态。如果是三重态的话，可能有其 ...

对，可能是单线态的能量，谢谢我算一下