计算化学公社

标题: 优化过程位移与受力波动导致优化无法成功 [打印本页]

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

作者
Author: ldatea 时间: 2020-11-6 14:44
本帖最后由 ldatea 于 2020-11-6 15:15 编辑

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

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

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

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

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