计算化学公社

标题: 明尼苏达泛函算SP第一次SCF不收敛报错中止但第二次却正常收敛，求问为何 [打印本页]

作者
Author: shenchaoren 时间: 2020-4-5 18:50
标题: 明尼苏达泛函算SP第一次SCF不收敛报错中止但第二次却正常收敛，求问为何
本帖最后由 shenchaoren 于 2020-4-5 18:55 编辑

软件用G16，因为主要关注反应能垒和反应热，所以采用明尼苏达系列泛函（一般是M06或M06L，偶尔用MN15L）算有机金属络合物的体系，小基组（一般是6-31G(d,p)-SDD 或者 def2SVP）优化结构，大基组（一般是def2TZVP或者ma-def2TZVP甚至更大的基组）算单点能。已经多次出现如下状况：算SP的时候，第一次提交作业，128圈SCF迭代后，不收敛报错中止。重新提交作业，作业里的关键字未做任何变动（积分格点精度是G16默认的ultrafine），SCF顺利收敛，任务完成。也用过B3LYP或B3PW91做类似的流程，但在算SP时从没有遇到过类似情况。请教各位高人，这是什么原因？我隐约感觉是明尼苏达系列泛函的问题，但是却不知道是什么原因，求解释

作者
Author: liyuanhe211 时间: 2020-4-5 19:59
除非关键词里有涉及读取初猜的情况，否则发生这种情况的概率很低。输出文件对比看看就知道了

作者
Author: sobereva 时间: 2020-4-6 19:53
本来明尼苏达系列泛函SCF收敛难度整体就比其它泛函要大，而且M06L之类纯泛函普遍也不如杂化泛函容易收敛，所以碰见SCF不收敛很正常。
如果完全是同一个输入文件，一次收敛一次不收敛，有一定数值巧合性（并行运算导致的），也可能关键词使用有问题，如2L所述。

作者
Author: shenchaoren 时间: 2020-4-18 19:39
本帖最后由 shenchaoren 于 2020-4-18 19:42 编辑

sobereva 发表于 2020-4-6 19:53
本来明尼苏达系列泛函SCF收敛难度整体就比其它泛函要大，而且M06L之类纯泛函普遍也不如杂化泛函容易收敛， ...

谢谢卢老师的解答，我这里给了一个具体例子，同一个优化过的结构，两次算单点能，结构完全一样，关键词没有任何变化，但是计算出的EE相差很多

作者
Author: sobereva 时间: 2020-4-19 00:31

shenchaoren 发表于 2020-4-18 19:39
谢谢卢老师的解答，我这里给了一个具体例子，同一个优化过的结构，两次算单点能，结构完全一样，关键词没 ...

这个属于极度罕见的巧合，两次结果不同和并行计算有关，见
数值误差对计算化学结果重现性的影响
http://sobereva.com/88

用能量更低的那个

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