计算化学公社

标题: 用高斯算单点能，收敛限放宽100倍时，能量差值已满足收敛限，但是还是输出收敛失败 [打印本页]

作者
Author: zhongyuabc 时间: 2017-12-18 11:09
标题: 用高斯算单点能，收敛限放宽100倍时，能量差值已满足收敛限，但是还是输出收敛失败
一个问题，望伙伴们解答，谢谢！
（1）用高斯算单点能，把收敛限放宽到conver=6,最后的收敛限达到了，但是还是输出收敛失败。图1是输出文件，图2是部分输出文件，附件是压缩的log文件

作者
Author: liyuanhe211 时间: 2017-12-18 13:10
分明写着convg=0.9E-4，哪儿满足收敛限了

作者
Author: zhongyuabc 时间: 2017-12-18 15:14

liyuanhe211 发表于 2017-12-18 13:10
分明写着convg=0.9E-4，哪儿满足收敛限了

后面几个cycle都基本是DE=xx-08,为什么Convg还是-4勒？

作者
Author: zjxitcc 时间: 2017-12-18 15:18

zhongyuabc 发表于 2017-12-18 15:14
后面几个cycle都基本是DE=xx-08,为什么Convg还是-4勒？

不是看DE。收敛条件还没满足。看>>>>>>之后的Convg

作者
Author: sobereva 时间: 2017-12-18 16:07
翻译：Convergence criterion not met 明显没达到收敛限

MaxDP高达2.83D-02，对照 Requested convergence on MAX density matrix=1.00D-04.，显然不可能判定成收敛

此例充分说明maxcyc=几百纯粹糟蹋计算资源。

作者
Author: zhongyuabc 时间: 2017-12-18 17:07

zjxitcc 发表于 2017-12-18 15:18
不是看DE。收敛条件还没满足。看>>>>>>之后的Convg

好的，谢谢。我最初以为是看DE，但是又看到下面有一个Convg,就有些糊涂他们之间的关系。

作者
Author: zhongyuabc 时间: 2017-12-18 17:11

sobereva 发表于 2017-12-18 16:07
翻译：Convergence criterion not met 明显没达到收敛限

MaxDP高达2.83D-02，对照 Requested convergenc ...

社长大大，那从MaxDP来看，是否意味着增加cycle是无用之举？那如果convg=xxE-5,增加cycle是否就有可能可行呢？

作者
Author: zjxitcc 时间: 2017-12-18 17:22

zhongyuabc 发表于 2017-12-18 17:11
社长大大，那从MaxDP来看，是否意味着增加cycle是无用之举？那如果convg=xxE-5,增加cycle是否就有可能可 ...

只有当观察到迭代最后几十圈的能量是在一直下降时，增加最大圈数才是有意义的。
另外，我感觉你想通过调升收敛标准来使结果收敛，这种做法是掩耳盗铃，使计算结果误差增大或者无意义。有些网站有人建议的这种做法：放松收敛标准，应该予以批判，十分不可取。
而是应该积极改善初猜，检查结构是否合理，寻找SCF不收敛的解决办法，最后无路可走，才考虑选择放松收敛标准，但是，单点不能低于10-6，优化不能低于10-7，否则结果可信度是比较低的。

作者
Author: zhongyuabc 时间: 2017-12-18 20:01

zjxitcc 发表于 2017-12-18 17:22
只有当观察到迭代最后几十圈的能量是在一直下降时，增加最大圈数才是有意义的。
另外，我感觉你想通过调 ...

好的，明白了，非常谢谢您

作者
Author: sobereva 时间: 2017-12-18 20:53

zhongyuabc 发表于 2017-12-18 17:11
社长大大，那从MaxDP来看，是否意味着增加cycle是无用之举？那如果convg=xxE-5,增加cycle是否就有可能可 ...

都不用看这个，99.999%的情况增加Gaussian默认的cyc上限都是无用的，而那0.001%能靠增大maxcyc解决的情况，总有其它更理想的做法解决。

作者
Author: zhongyuabc 时间: 2017-12-18 21:00

sobereva 发表于 2017-12-18 20:53
都不用看这个，99.999%的情况增加Gaussian默认的cyc上限都是无用的，而那0.001%能靠增大maxcyc解决的情 ...

谢谢社长大大

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