请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 100|回复: 4

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

[复制链接]

25

帖子

0

威望

1189

eV
积分
1214

Level 4 (黑子)

发表于 6 天前 | 显示全部楼层 |阅读模式
我用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是不简并的。
输入输出见附件。
不知道是怎么回事,重新算了还是一样的结果。
另外有个分子也发现了同样的情况。


输入

输入

输出

输出

output.out

65.04 KB, 下载次数: 0

input.com

1.47 KB, 下载次数: 0

411

帖子

16

威望

2768

eV
积分
3499

Level 5 (御坂)

发表于 6 天前 | 显示全部楼层
如果ci/pt2步骤的根数少于前面mcscf步骤的根数,经常会出现这个问题。

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

评分

参与人数 1eV +2 收起 理由
sobereva + 2

查看全部评分

25

帖子

0

威望

1189

eV
积分
1214

Level 4 (黑子)

 楼主| 发表于 5 天前 | 显示全部楼层
本帖最后由 胡说 于 2019-10-9 13:15 编辑
beefly 发表于 2019-10-8 18:38
如果ci/pt2步骤的根数少于前面mcscf步骤的根数,经常会出现这个问题。

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

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

411

帖子

16

威望

2768

eV
积分
3499

Level 5 (御坂)

发表于 5 天前 | 显示全部楼层
胡说 发表于 2019-10-9 09:59
多谢老师,确实有效,解决了问题。
但是发现加了nstati后,耗时会增加,有些增加得还蛮多。
这个选项是 ...

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

25

帖子

0

威望

1189

eV
积分
1214

Level 4 (黑子)

 楼主| 发表于 4 天前 | 显示全部楼层
beefly 发表于 2019-10-9 17:50
如果默认能算对就不用加了

好的多谢老师
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-10-14 22:21 , Processed in 0.162713 second(s), 27 queries .

快速回复 返回顶部 返回列表