计算化学公社

标题: 体系波函数检测稳定后仍然不收敛 [打印本页]

作者
Author: somnus133937 时间: 2021-8-29 10:43
标题: 体系波函数检测稳定后仍然不收敛
通过使用stable=opt关键词，得到稳定波函数后，在其基础上进行单点能的计算，但仍然出现了收敛失败的现象

麻烦问一下各位老师，为什么在稳定波函数基础下计算为什么结果还会出现收敛失败的问题呢？

作者
Author: sobereva 时间: 2021-8-29 11:33
直接把两个任务用的关键词贴出来，最好输出文件压缩后也上传

作者
Author: wzkchem5 时间: 2021-8-29 15:40
先搞清楚一个概念，收敛是否失败和波函数是否稳定没有必然联系。只要单点能读取了已经收敛的波函数（不管该波函数是否稳定），且单点能计算级别和读取的这个波函数用的计算级别严格一致，就必然可以收敛。如果两者的计算级别不一致，那么即使初猜的波函数是稳定的，也可能不收敛，且不收敛的概率可能并不比不稳定波函数低

作者
Author: somnus133937 时间: 2021-8-30 08:22
这是计算的gjf和输出文件

作者
Author: somnus133937 时间: 2021-8-30 08:23
本帖最后由 somnus133937 于 2021-8-30 08:24 编辑

wzkchem5 发表于 2021-8-29 15:40
先搞清楚一个概念，收敛是否失败和波函数是否稳定没有必然联系。只要单点能读取了已经收敛的波函数（不管该 ...

计算级别是保持了一致，但我在进行stable=opt的计算时，加了scf=xqc帮助收敛，是这个原因导致的后续单点能计算存在不收敛的现象吗？

作者
Author: somnus133937 时间: 2021-8-30 08:25

sobereva 发表于 2021-8-29 11:33
直接把两个任务用的关键词贴出来，最好输出文件压缩后也上传

好的老师，文件以上传。

作者
Author: sobereva 时间: 2021-8-30 11:29
你真要算三重态的Li2？
为什么用MN12L？
当前的计算也太奇怪了

作者
Author: somnus133937 时间: 2021-8-30 14:18

sobereva 发表于 2021-8-30 11:29
你真要算三重态的Li2？
为什么用MN12L？
当前的计算也太奇怪了

确实挑的体系和泛函比较特殊，我只是比较好奇为什么会出现波函数稳定后仍然不收敛的情况，

sob老师，请问是因为这个体系和泛函比较特殊的原因吗？

作者
Author: wzkchem5 时间: 2021-8-30 15:58

somnus133937 发表于 2021-8-30 01:23
计算级别是保持了一致，但我在进行stable=opt的计算时，加了scf=xqc帮助收敛，是这个原因导致的后续单点 ...

加一个level shift（比如0.3 au）再试。第二步迭代能量突然上升，且一直降不下来，这个一般是意味着有的轨道的顺序反转了，所以加level shift阻止轨道顺序变化。
如果还不行，增大格点再试，Minnesota泛函对格点比较敏感

作者
Author: somnus133937 时间: 2021-8-30 16:11

wzkchem5 发表于 2021-8-30 15:58
加一个level shift（比如0.3 au）再试。第二步迭代能量突然上升，且一直降不下来，这个一般是意味着有的 ...

好的老师，我去试试，谢谢您提的建议

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