计算化学公社

标题: 求助用scan刚性扫描两个键角的变化出现l502收敛失败的错误，该怎么解决？ [打印本页]

作者
Author: leee 时间: 2023-7-3 12:49
标题: 求助用scan刚性扫描两个键角的变化出现l502收敛失败的错误，该怎么解决？
本帖最后由 leee 于 2023-7-3 13:29 编辑

各位老师，您好，我在做C3O2-与CO分子的范德华力构型的势能面扫描时，设置了A2（50.0 50 1.0 ）和A5（144.0 1 1.0）两个角度同时变化，当扫描到A5=144, A2=98时显示收敛失败，但是当我设置缩小角度变化范围A2（97.0 3 1.0），A5（143.0 3 1.0）就能正常结束了，直接只算这一个单点能也能正常结束，请问这是什么原因导致的？谢谢各位老师！

作者
Author: wzkchem5 时间: 2023-7-3 15:30
势能面扫描的时候，程序读取上一步的结构的波函数作为初猜。显然，两个结构越接近，它们的收敛的波函数一般而言也就越接近。所以当步长足够小的时候，上一步收敛的波函数放在这一步也几乎是收敛的，所以收敛很容易；当步长很大的时候，上一步收敛的波函数放在这一步离收敛很远，所以收敛很难甚至不收敛。

作者
Author: leee 时间: 2023-7-3 16:48

wzkchem5 发表于 2023-7-3 15:30
势能面扫描的时候，程序读取上一步的结构的波函数作为初猜。显然，两个结构越接近，它们的收敛的波函数一般 ...

谢谢老师的回复！您说的话我认真的读了好几遍有点明白这其中的道理了
但是我还有一个疑问，为什么我把扫描的两个角度范围变小，步长不变，它就能收敛呢?

投稿之后我尝试了2种方法：
第一这期间我也尝试了步长改为0.5，还是不收敛。可能还要往小了改吧。
第二我加了guess(always,mix)（参考了http://sobereva.com/474）,这会倒是收敛了。

作者
Author: wzkchem5 时间: 2023-7-3 19:36

leee 发表于 2023-7-3 09:48
谢谢老师的回复！您说的话我认真的读了好几遍有点明白这其中的道理了
但是我还有一个疑问，为什么我把扫 ...

可能是偶然因素，你这个体系恰好处在能收敛和不能收敛的边缘上，所以任何偶然扰动都可以把收敛变成不收敛，或反之

作者
Author: leee 时间: 2023-7-3 19:47

wzkchem5 发表于 2023-7-3 19:36
可能是偶然因素，你这个体系恰好处在能收敛和不能收敛的边缘上，所以任何偶然扰动都可以把收敛变成不收敛 ...

好的，十分感谢老师

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