计算化学公社

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

作者
Author: 秋心 时间: 2024-8-3 12:11
标题: [已解决]求助Molpro输入文件中对带电分子的电荷设置问题
本帖最后由秋心于 2024-8-5 11:05 编辑

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

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

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