计算化学公社

标题: 关于双自由基算法求助 [打印本页]

作者
Author: kevin 时间: 2016-11-22 16:32
标题: 关于双自由基算法求助
各位老师，
   在计算双自由基时，我有一些疑惑，我看过了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
    <分子坐标>


   谢谢各位老师！

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

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

作者
Author: kevin 时间: 2016-12-14 16:04

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看看和文献里是不是相同

作者
Author: yezhonghua 时间: 2016-12-14 21:39

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)