求助，Molpro是否不支持开壳层系统CASPT2的解析梯度？

commanderwilson · 发表于 Post on 2025-5-28 02:13:37

事情是这样的，我准备在CASPT2(3,6)/def2-TZVPP级别下对[BeH2]+阳离子做几何优化，一开始我想先用CASSCF进行运算，计算体系的³B₂态，输入文件：

1.inp (209 Bytes, 下载次数 Times of downloads: 0)
这一步是可以正常输出的，优化的结果也基本正确：

1.out (17.6 KB, 下载次数 Times of downloads: 0)
接下来我加入rs2输入块，试图调用CASPT2梯度的时候，输入文件：

2.inp (210 Bytes, 下载次数 Times of downloads: 0)
然后出现了这样的报错：

2.out (21.57 KB, 下载次数 Times of downloads: 3)
我还不是很熟悉Molpro，求大佬们讲讲为什么会这么报错

beefly · 发表于 Post on 2025-5-28 19:09:15

本帖最后由 beefly 于 2025-5-28 19:19 编辑

中性分子wf,6,1,0、阴离子wf,7,3,1都可以算，唯独阳离子wf,5,3,1（三个活性电子）出错。说明是程序bug：rs2梯度代码在处理四电子密度矩阵的时候，没有考虑少于四个活性电子的特殊情况。临时的解决方法是人为增加活性电子数，例如把closed空间的2个电子放入活性空间

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

[Molpro] 求助，Molpro是否不支持开壳层系统CASPT2的解析梯度？

评分 Rate

浏览过的版块