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

DevilYuan · 发表于 Post on 2026-4-18 22:59:07

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

haoguoyuerban · 发表于 Post on 2026-6-7 11:10:47

我最近也用DFT扫描金属配合物的PCET反应也遇到了类似问题，请问您得到解决方案了吗？谢谢！

haoguoyuerban · 发表于 Post on 2026-6-7 11:16:10

haoguoyuerban 发表于 2026-6-7 11:10
我最近也用DFT扫描金属配合物的PCET反应也遇到了类似问题，请问您得到解决方案了吗？谢谢！

CDFT(constrained DFT)可能是一种办法。您用的chemshell连接的是什么QM软件?

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

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