计算化学公社

标题: 求助，Molpro是否不支持开壳层系统CASPT2的解析梯度？ [打印本页]

作者
Author: commanderwilson 时间: 2025-5-28 02:13
标题: 求助，Molpro是否不支持开壳层系统CASPT2的解析梯度？
事情是这样的，我准备在CASPT2(3,6)/def2-TZVPP级别下对[BeH2]+阳离子做几何优化，一开始我想先用CASSCF进行运算，计算体系的³B₂态，输入文件：
(, 下载次数 Times of downloads: 0)
这一步是可以正常输出的，优化的结果也基本正确： (, 下载次数 Times of downloads: 0)
接下来我加入rs2输入块，试图调用CASPT2梯度的时候，输入文件： (, 下载次数 Times of downloads: 0)
然后出现了这样的报错： (, 下载次数 Times of downloads: 3)
我还不是很熟悉Molpro，求大佬们讲讲为什么会这么报错

作者
Author: beefly 时间: 2025-5-28 19:09
本帖最后由 beefly 于 2025-5-28 19:19 编辑

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

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