计算化学公社

标题: Gaussian过渡态优化，Maximum Force不收敛 [打印本页]

作者
Author: 王一伊 时间: 2021-5-5 20:28
标题: Gaussian过渡态优化，Maximum Force不收敛
各位老师好，我在找过渡态的过程中 Maximum Displacement，RMS Displacement 两项都已经收敛，RMS Force 数值刚好是0.000300卡在收敛限上，而 Maximum Force 一直无法收敛，且在后面的几十步的计算中各项数值都是不变的。目前已经尝试了调整结构、使用不同的关键词进行优化，所得结果都相同。也尝试了在优化过程中降低收敛限（opt=loose），频率计算震动方向正确，但是接下来的IRC只能往一个方向跑。不知道自己这种解决办法是否合理，请问这种情况下该怎么办呢？
图片是优化最后一步（当然前面几十步也相同）的数据

(, 下载次数 Times of downloads: 32)

(, 下载次数 Times of downloads: 35)

附件是开始的时候使用的输入文件

(, 下载次数 Times of downloads: 147)

作者
Author: wzkchem5 时间: 2021-5-5 21:24
以当前结构优化的最后一步作为初始结构重新优化。
如果还不行，试试提高积分精度。有可能是积分精度太低，能量和梯度不匹配。

作者
Author: 王一伊 时间: 2021-5-5 22:21

wzkchem5 发表于 2021-5-5 21:24
以当前结构优化的最后一步作为初始结构重新优化。
如果还不行，试试提高积分精度。有可能是积分精度太低， ...

谢谢老师回复我试一下~

作者
Author: sobereva 时间: 2021-5-6 02:37
体系这么小，显然应当考虑用recalc或calcall帮助收敛
而且通过设小maxstep帮助收敛是极为常见的做法
没必要写cartesian，只会令收敛更费劲
你当前给的能量变化图什么有用的也看不出来，显然应当把中、后期部分放大了看

作者
Author: 王一伊 时间: 2021-5-6 14:12

sobereva 发表于 2021-5-6 02:37
体系这么小，显然应当考虑用recalc或calcall帮助收敛
而且通过设小maxstep帮助收敛是极为常见的做法
没必 ...

谢谢老师回复，按您给的方法解决啦！能量图的问题我以后注意~

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