计算化学公社

标题: 求助双自由基体系用CCSD(T)计算能量比闭壳层单重态能量高 [打印本页]

作者
Author: MageArcane 时间: 2019-5-9 15:33
标题: 求助双自由基体系用CCSD(T)计算能量比闭壳层单重态能量高
计算一个有机分子，用opt=stable进行波函数稳定性检测，得到了双自由基单重态波函数。之后直接用geom=read,guess=check 进行CCSD(T)能量计算发现其能量高于闭壳层的单重态，进行T1诊断值为0.032左右，是不是说明CCSD(T)的单参考态不准确，怎么解决呢？是否需要用CASSCF，CASPT2等方法？

作者
Author: sobereva 时间: 2019-5-9 21:28
看看对称破缺基态波函数的UHF能量和RHF能量相差多少

作者
Author: MageArcane 时间: 2019-5-10 10:09

sobereva 发表于 2019-5-9 21:28
看看对称破缺基态波函数的UHF能量和RHF能量相差多少

SCF Done: E(RB3LYP) = -412.576855235 A.U. after 1 cycles CCSD(T)=-411.4136384
SCF Done: E(UB3LYP) = -412.579847764 A.U. after 1 cycles CCSD(T)=-411.4061666
UHF比RHF稳定1.8kcal/mol 左右，但是CCSD(T)结果是RHF的稳定4.7kcal/mol。

作者
Author: sobereva 时间: 2019-5-10 10:56

MageArcane 发表于 2019-5-10 10:09
SCF Done: E(RB3LYP) = -412.576855235 A.U. after 1 cycles CCSD(T)=-411.4136384
SCF ...

你给的是B3LYP的能量，而不是HF的能量

作者
Author: MageArcane 时间: 2019-5-10 11:23

sobereva 发表于 2019-5-10 10:56
你给的是B3LYP的能量，而不是HF的能量

您指的是CCSD(T)中HF的能量吗？
分别是 SCF Done: E(RHF) = -409.796234785 A.U. after 16 cycles
SCF Done: E(UHF) = -409.838536530 A.U. after 1 cycles

作者
Author: sobereva 时间: 2019-5-10 16:30

MageArcane 发表于 2019-5-10 11:23
您指的是CCSD(T)中HF的能量吗？
分别是 SCF Done: E(RHF) = -409.796234785 A.U. after 16 cycl ...

是

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