求助降低MM对QM的梯度的办法

科研狗

我在使用COBRAMM程序作为接口连接Amber跟Gaussian去做QM/MM模拟，想冻结高层去优化中层。在脚本中找到了command[130]=1可以冻结高层原子，但是当我输入这个命令并且提交任务（用的方法是pbe0/6-31G(d)）后，任务结束 ，发现高层原子构型发生了变化，并没有冻结，

我与开发者联系是因为高层跟中层挨得太近，相互作用太强，导致梯度太大。我试着先用gaussian的半经验方法跟mopac的半经验方法去做预优化，结果与前面的相似，只是构型没有像前面的那样发生特别大的变化，
我想问一下还有什么方法可以使这种梯度降低？
谢谢您的解答。