计算化学公社

标题: 求助：Molpro-CASPT2计算过程中计算了同一个态 [打印本页]

作者
Author: 胡说 时间: 2019-10-8 13:22
标题: 求助：Molpro-CASPT2计算过程中计算了同一个态
我用molpro2012计算一个分子的激发态能量，先是SA(4)-CASSCF(12,9),然后计算MS-CASPT2.输入参照molpro手册：“Performing SS-SR-CASPT2 calculations”这一节.
但是发现CASPT2计算过程中在算Root2和Root3的时候参考态都是CASSCF里的Root3，导致CASPT2算得的Root2和Root3能量完全一致，也就是算的同一个态。但CASSCF的结果Root2和Root3是不简并的。
输入输出见附件。
不知道是怎么回事，重新算了还是一样的结果。
另外有个分子也发现了同样的情况。

作者
Author: beefly 时间: 2019-10-8 18:38
如果ci/pt2步骤的根数少于前面mcscf步骤的根数，经常会出现这个问题。

在ci/pt2中加入以下选项，其中xxx是mcscf中的根数
option,nstati=xxx;

作者
Author: 胡说 时间: 2019-10-9 09:59
本帖最后由胡说于 2019-10-9 13:15 编辑

beefly 发表于 2019-10-8 18:38
如果ci/pt2步骤的根数少于前面mcscf步骤的根数，经常会出现这个问题。

在ci/pt2中加入以下选项，其中xxx ...

多谢老师，确实有效，解决了问题。
但是发现加了nstati后，耗时会增加，有些增加得还蛮多。
这个选项是否可以只加在程序识别参考态错误的那个态上呢，其他比如root1不加。根据我的算的一些分子，发现CASSCF阶段如果有两个根能量很接近时会出现这种情况。
从能量偶极矩上看，“option，nstati”这个选项加与不加，差别微乎其微。

作者
Author: beefly 时间: 2019-10-9 17:50

胡说发表于 2019-10-9 09:59
多谢老师，确实有效，解决了问题。
但是发现加了nstati后，耗时会增加，有些增加得还蛮多。
这个选项是 ...

如果默认能算对就不用加了

作者
Author: 胡说 时间: 2019-10-10 15:07

beefly 发表于 2019-10-9 17:50
如果默认能算对就不用加了

好的多谢老师

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