[已解决]求助Molpro输入文件中对带电分子的电荷设置问题

秋心

各位老师好，我最近在做硫酸分子与硫酸氢根分子的SAPT计算，两个分子的总电子数都是50，自旋量子数(spin)都是0，硫酸电荷为0，而硫酸氢根电荷为-1，然后我根据手册中的定义方式分别对硫酸分子和硫酸氢根分子设置（如图），但是我看out文件中输出NUMBER OF ELECTRONS:      50+   49-    SPACE SYMMETRY=1    SPIN SYMMETRY: Doublet ，硫酸的电子数设置对了，但是硫酸氢根设置错了(少一个电子)，并且整个体系的自旋量子数是错误的，这样继续跑下去就会遇到open-shell不支持的报错，麻烦各位老师帮我看一下我的输入文件那里有问题吗？感激不尽！


解决办法：换用PSI4。我上面遇到的问题大概率是molpro程序本身的问题，多次测试molpro都没办法正确识别体系的电子数目，电荷，自旋量子数，可能我考虑的还不到位，但是继续测试下去时间成本太大且无意义，果断更换性能更好且好用的PSI4。