计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6795|回复 Reply: 6

[量化理论] MOLCAS XMS-CASPT2测试

[复制链接 Copy URL]

597

帖子

20

威望

4153

eV
积分
5150

Level 6 (一方通行)

发表于 Post on 2018-12-29 14:29:08 | 显示全部楼层 Show all |阅读模式 Reading model
本帖最后由 beefly 于 2018-12-29 17:19 编辑

在很多近简并态的情况下,过去的多态多参考二阶微扰理论(主要是MCQDPT2和MS-CASPT2)可能会导致势能面存在个别不连续点。针对此问题,Firefly(PC-GAMESS)作者Granovsky在2011年提出了XMCQDPT2(J. Chem. Phys. 134, 214113, 2011),随后Molpro实现了XMS-CASPT2(J. Chem. Phys. 135, 081106, 2011)。MOLCAS(OpenMOLCAS)虽然几年前也加入了XMS-CASPT2,但是迟迟未写进使用手册。以下是OpenMOLCAS的XMS-CASPT2和MS-CASPT2测试,计算ScO分子最低十个2Π态的势能曲线(局部),两种方法使用相同的CASSCF波函,可见XMS-CASPT2仍然有问题。也许这是MOLCAS未正式公布XMS-CASPT2的原因吧。
这个例子并不说明MS-CASPT2结果好。MS-CASPT2在< 1.9 Angstrom问题更多。
Graph1.png

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
sobereva + 5

查看全部评分 View all ratings

2481

帖子

11

威望

6623

eV
积分
9324

Level 6 (一方通行)

发表于 Post on 2018-12-29 15:26:38 | 显示全部楼层 Show all
why is the left plot blurred?

597

帖子

20

威望

4153

eV
积分
5150

Level 6 (一方通行)

 楼主 Author| 发表于 Post on 2018-12-29 17:16:59 | 显示全部楼层 Show all
我本是个娃娃 发表于 2018-12-29 15:26
why is the left plot blurred?

有几个可能原因:程序有bug;XMS-CASPT2有问题;XMS-CASPT2、XMCQDPT2都有问题。

214

帖子

0

威望

4840

eV
积分
5054

Level 6 (一方通行)

跳跳猪

发表于 Post on 2018-12-30 00:50:25 | 显示全部楼层 Show all
可以试试molpro或者Bagel的MS-CASPT2和XMS-CASPT2。
算这么多态可能需要更大的active space,态多的时候有可能CASSCF就没收敛想要的态上。
我用Gamess,Bagel和Molpro的时候没发现XMS和MS差别这么大(当然没考虑这么多态),原理上XMS比MS精度要好一些,不过没试过Molcas的
流年似水,浮生如梦。

597

帖子

20

威望

4153

eV
积分
5150

Level 6 (一方通行)

 楼主 Author| 发表于 Post on 2018-12-30 09:24:45 | 显示全部楼层 Show all
Mikasa 发表于 2018-12-30 00:50
可以试试molpro或者Bagel的MS-CASPT2和XMS-CASPT2。
算这么多态可能需要更大的active space,态多的时候有 ...

因为后面要做旋轨耦合,能用xms矫正后的casscf波函做soc的程序,目前只有molcas

597

帖子

20

威望

4153

eV
积分
5150

Level 6 (一方通行)

 楼主 Author| 发表于 Post on 2019-3-2 09:51:43 | 显示全部楼层 Show all
本帖最后由 beefly 于 2019-3-2 10:20 编辑

补记:虽然molpro的xms-caspt2非常慢,但还是试了一下,得到的势能曲线更奇怪。一个可能是xms理论本身有问题(或者只能用于pc-gamess的mcqdpt2),另一个可能是molpro、molcas的xms代码有问题。
molpro的ms-caspt2结果正常。前面说molcas的ms-caspt2在1.9 Angstrom之前的结果有问题。检查后发现是molcas的casscf收敛存在问题,与mx-caspt2无关。

XMS-molpro

XMS-molpro

43

帖子

1

威望

476

eV
积分
539

Level 4 (黑子)

发表于 Post on 2019-3-23 22:26:40 | 显示全部楼层 Show all
本帖最后由 I10140317 于 2019-3-25 09:26 编辑
beefly 发表于 2019-3-2 09:51
补记:虽然molpro的xms-caspt2非常慢,但还是试了一下,得到的势能曲线更奇怪。一个可能是xms理论本身有问 ...

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2023-2-7 03:09 , Processed in 0.423415 second(s), 25 queries .

快速回复 返回顶部 返回列表 Return to list