计算化学公社

标题: 【burai】请教一下优化切面以后的ZnFe2O4结构收敛很慢的原因 [打印本页]

作者
Author: 一束光 时间: 2021-5-24 12:54
标题: 【burai】请教一下优化切面以后的ZnFe2O4结构收敛很慢的原因
各位大佬好，求助一个问题，我在优化切面以后的ZnFe2O4结构时（mixing beta设置的0.04），发现结构收敛非常慢，计算了8个多小时，走了130多步，但是能量还在一点点下降。收敛的非常慢，请问这种情况正常吗？如果不正常应该采取什么方法加速收敛啊？谢谢

作者
Author: 卡开发发 时间: 2021-5-24 13:02

发现结构收敛非常慢

到底是构型优化收敛慢还是自洽场迭代收敛慢？

作者
Author: 喵星大佬 时间: 2021-5-24 13:50
说结构收敛非常慢，然后贴一个SCF收敛慢的图上来，是什么意思

作者
Author: 一束光 时间: 2021-5-24 14:26

卡开发发发表于 2021-5-24 13:02
到底是构型优化收敛慢还是自洽场迭代收敛慢？

不好意思哈老师，我用BURAI选择optimize，然后结果里就是在跑这个收敛，然后收敛以后会重新再跑一次SCF。这个图第一个，应该是构型优化。然后第二次应该是自洽场迭代。

作者
Author: 一束光 时间: 2021-5-24 14:27

喵星大佬发表于 2021-5-24 13:50
说结构收敛非常慢，然后贴一个SCF收敛慢的图上来，是什么意思

感谢您的回复哈，不好意思。我用的是burai计算，，选择的optimize，他一共会进行两次scf计算，第一次应该是构型优化，收敛以后会再跑一次（这个应该是自洽场迭代计算）。这个是第一个图。

作者
Author: 卡开发发 时间: 2021-5-24 14:38

一束光发表于 2021-5-24 14:26
不好意思哈老师，我用BURAI选择optimize，然后结果里就是在跑这个收敛，然后收敛以后会重新再跑一次SCF。 ...

不如直接把out文件发上来算了

作者
Author: 喵星大佬 时间: 2021-5-24 15:04
我觉得你应该是很多常识没搞清楚

作者
Author: 一束光 时间: 2021-5-24 15:10

卡开发发发表于 2021-5-24 14:38
不如直接把out文件发上来算了

好的，谢谢老师，1号是我为了看能不能收敛把convergence threshold设成了0.0001的结果，也是跑了130多步。第二个是我改变了输入结构异构的输出文件，每次计算到63步就会又重新计算。可能因为我有很多基础东西没弄懂，问的题没表达清楚，给您添麻烦了。

作者
Author: 一束光 时间: 2021-5-24 15:12

喵星大佬发表于 2021-5-24 15:04
我觉得你应该是很多常识没搞清楚

嗯嗯，谢谢您的指导。因为还是个菜鸟，刚起步。很多问题可能不专业，请您多包涵。请问您什么好的补充基础类的帖子和视频链接吗？非常感谢

作者
Author: 卡开发发 时间: 2021-5-24 15:32

一束光发表于 2021-5-24 15:10
好的，谢谢老师，1号是我为了看能不能收敛把convergence threshold设成了0.0001的结果，也是跑了130多步 ...

其实我还是不知道你指的130步到底是自洽场的步数还是构型优化的步数，你要是找得到out文件的话再一块研究。

作者
Author: 一束光 时间: 2021-5-24 16:03

卡开发发发表于 2021-5-24 15:32
其实我还是不知道你指的130步到底是自洽场的步数还是构型优化的步数，你要是找得到out文件的话再一块研究 ...

老师，我在问题下面修改了回答，贴上了两个输出文件，您看是吗？

作者
Author: 卡开发发 时间: 2021-5-24 16:50
本帖最后由卡开发发于 2021-5-24 16:52 编辑

一束光发表于 2021-5-24 15:10
好的，谢谢老师，1号是我为了看能不能收敛把convergence threshold设成了0.0001的结果，也是跑了130多步 ...

给的两个文件看上去是SCF收敛比较慢，构型优化那个没跑完。SCF收敛缓慢对于QE主要能够调整的是mixing beta和mixing mode以及DIIS，目前先建议调整mixing mode，以及DIIS也可以适当放大些。另外算力有限的话可以考虑适当减少点真空，当然真空大小最好好好测试下多少才使得能量偏差比较小。如果能够把c轴做成垂直于ab的模型，可以考虑做assume_isolated='esm'和esm='bc1'，这样应该可以使用稍微小一些的真空层。

基本的概念可以慢慢补充，但只要是要继续做计算就应该多了解一些。构型优化和自洽计算的关系，一般而言一步自洽计算时求解确定结构的KS方程得到波函数、电子密度、能量、力等信息，然后构型优化则是由n个自洽计算组成的，通过自洽计算得到的能量（函数）和力（导数）来找能量为极小（未必是最小）的原子移动方向，然后移动原子产生新的结构，继续通过自洽如此往复直至力和能量满足收敛标准。

另外材料的结构和磁性我没办法帮你把握，我不知道你具体做的是ZnFeO的那种材料，这个可能需要你自己了解更多。

作者
Author: 一束光 时间: 2021-5-24 17:49

卡开发发发表于 2021-5-24 16:50
给的两个文件看上去是SCF收敛比较慢，构型优化那个没跑完。SCF收敛缓慢对于QE主要能够调整的是mixing bet ...

非常感谢老师您的耐心解答，对我帮助很大。

作者
Author: 杨凤磊 时间: 2021-10-21 16:55

一束光发表于 2021-5-24 17:49
非常感谢老师您的耐心解答，对我帮助很大。

您好，我刚开始尝试用QE算吸附，我能请教一下您如何建模吗？

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