计算化学公社

标题: 如何计算两种单线态氧 [打印本页]

作者
Author:
amagi    时间: 2018-9-29 00:53
标题: 如何计算两种单线态氧
其中一种似乎可以这样算,另一种需要怎么算呢,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





作者
Author:
sobereva    时间: 2018-9-29 10:49
nroot是总共求解的根数,默认是1,显然只能给你0 1设定下的能量最低态。你要算两个必须nroot=2。算多个根的时候一般结合态平均使用,怎么用用首页google框搜索论坛便知
作者
Author:
amagi    时间: 2018-9-29 11:38
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一模一样
觉得很疑惑
作者
Author:
sobereva    时间: 2018-9-29 12:59
amagi 发表于 2018-9-29 11:38
感谢大神回复!
我先是计算了基态0 3设定下 #p CASSCF(8,6)/6-31g* opt freq ,得到零点能(-149.691123 ...

算出的态的自旋多重度是单重态还是三重态是通过0 1、0 3来控制的,跟nroot没关系。nroot设定的是对于你指定的自旋多重度算出多少个态来
作者
Author:
fzdzz    时间: 2021-9-23 10:10
amagi 发表于 2018-9-29 11:38
感谢大神回复!
我先是计算了基态0 3设定下 #p CASSCF(8,6)/6-31g* opt freq ,得到零点能(-149.691123 ...

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3