求助：Molpro-CASPT2计算过程中计算了同一个态

胡说 · 发表于 Post on 2019-10-8 13:22:40

我用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是不简并的。
输入输出见附件。
不知道是怎么回事，重新算了还是一样的结果。
另外有个分子也发现了同样的情况。

beefly · 发表于 Post on 2019-10-8 18:38:35

如果ci/pt2步骤的根数少于前面mcscf步骤的根数，经常会出现这个问题。

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

胡说 · 发表于 Post on 2019-10-9 09:59:25

本帖最后由胡说于 2019-10-9 13:15 编辑

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

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

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

beefly · 发表于 Post on 2019-10-9 17:50:25

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

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

胡说 · 发表于 Post on 2019-10-10 15:07:40

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

好的多谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Molpro] 求助：Molpro-CASPT2计算过程中计算了同一个态

评分 Rate

浏览过的版块