优化过程位移与受力波动导致优化无法成功

王寓于

我优化的体系有344个原子（C,H,O,N,S,Br),用的Gaussian09，wb97xd/6-31g**.最开始优化的时候最大位移和最大受力都趋向收敛，但一段时间后又开始震荡，之后提取了结构最好的一帧重新优化发现最大受力曲线平坦，但位移曲线始终波动，我参考了  量子化学计算中帮助几何优化收敛的常用方法（http://sobereva.com/164）添加了opt=maxstep=5，开始几天优化曲线还好，但今天又开始有了震荡的趋势（如下图），请问我还是继续等吗？谢谢。

ldatea

opt=maxstep=3甚至更低都是可以考虑的。

可以考虑opt=modredundant，手动加几个冗余内坐标，算几个冗余内坐标的数值Hessian。但这个我没操作过，只是手册上看到过。

如果计算量和内存承受得起，建议用opt=recalc=x，效果是每x步计算一次Hessian,一般3-5比较合适。
再降低就是opt=calcfc只算一次Hessian。不过我感觉300多原子的Hessian计算量会很大很大，慎用。

sobereva

光写maxstep根本没用，必须再加上notrust
另外，势能面平坦的大体系建议加上gdiis
明显http://sobereva.com/164文中的知识你还没充分掌握，再一个字一个字反复多读几遍此文