求助，QM/MM计算蛋白PCET过程遇到能量曲线骤降的问题

DevilYuan

各位大佬，我现在在用chemshell的QM/MM计算蛋白的pcet过程，用反应路径扫描的方法，限制的是质子转移的两个键长之差，用的泛函是wb97xd，基组是6-31G**。反应是Y373+ W322·⁺ +D321- → Y373· +W322+D321·H。扫描结果如图所示。最高点计算了频率发现只有一个虚频。过了最高点后，能量和自旋分布会发生骤变，最高点前后的构象会有突然的变化，导致曲线不平滑。应该是发生了电子态切换，SSCF 跟到了另一支解，导致能量、构型和自旋分布一起突变。用B3LYP的话，初始态的两个残基的自旋密度分布不对，用M062X和CAM-B3LYP也是这个情况。我现在应该怎么做才能使得能量曲线和构象变化变得平滑，得到正确的反应路径能量曲线。