如何计算两种单线态氧

amagi

其中一种似乎可以这样算，另一种需要怎么算呢，nroot应该设成2还是3呢
#p CASSCF(8, 6)/6-31g* opt freq

O2_CASSCF

0 1
O                 -1.90697672    2.06976752    0.00000000
O                 -3.06857672    2.06976752    0.00000000

sobereva

nroot是总共求解的根数，默认是1，显然只能给你0 1设定下的能量最低态。你要算两个必须nroot=2。算多个根的时候一般结合态平均使用，怎么用用首页google框搜索论坛便知

amagi

sobereva 发表于 2018-9-29 10:49
nroot是总共求解的根数，默认是1，显然只能给你0 1设定下的能量最低态。你要算两个必须nroot=2。算多个根的 ...

感谢大神回复！
我先是计算了基态0 3设定下 #p CASSCF(8,6)/6-31g* opt freq ，得到零点能（-149.691123）；
然后计算了0 1设定下 #p CASSCF(8,6，nroot=2)/6-31g* opt freq ，用得到的零点能（-149.658621）减去0 3基态的零点能换算得到0.8844eV，和文献上计算的1Δg能量（0.890eV）相近；
而计算0 1设定下 #p CASSCF(8,6，nroot=3)/6-31g* opt freq，用得到的零点能（-149.632151）减去0 3基态的零点能换算得到1.6047eV，和文献上计算的1Σg+能量（1.6047eV）相近。
最后我再次计算基态0 3设定下#p CASSCF(8,6,nroot=1)/6-31g* opt freq，算出与前面基态不写nroot=1一样的零点能。
是否可以认为确认nroot数值需要不管自旋多重度差异直接叠加呢，比如我想算1Σg+（第二单重态），nroot应该写3（多算个基态三重态）而不是2（不算基态三重态）。
...................
其实我还顺便算了nroot=3     0 3报了l501错
                          nroot=2     0 3 算出-149.554710很高的值
                          nroot=1     0 1 算出-149.658621和nroot=2一模一样
觉得很疑惑

sobereva

amagi 发表于 2018-9-29 11:38
感谢大神回复！
我先是计算了基态0 3设定下 #p CASSCF(8,6)/6-31g* opt freq ，得到零点能（-149.691123 ...

算出的态的自旋多重度是单重态还是三重态是通过0 1、0 3来控制的，跟nroot没关系。nroot设定的是对于你指定的自旋多重度算出多少个态来

fzdzz

amagi 发表于 2018-9-29 11:38
感谢大神回复！
我先是计算了基态0 3设定下 #p CASSCF(8,6)/6-31g* opt freq ，得到零点能（-149.691123 ...

你好，请问关于你用高斯计算那个两种单线态氧的参考文献还能找到么？能发给我吗？