【burai】请教一下优化切面以后的ZnFe2O4结构收敛很慢的原因

一束光

各位大佬好，求助一个问题，我在优化切面以后的ZnFe2O4结构时（mixing beta设置的0.04），发现结构收敛非常慢，计算了8个多小时，走了130多步，但是能量还在一点点下降。收敛的非常慢，请问这种情况正常吗？如果不正常应该采取什么方法加速收敛啊？谢谢

卡开发发

发现结构收敛非常慢

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

喵星大佬

说结构收敛非常慢，然后贴一个SCF收敛慢的图上来，是什么意思

一束光

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

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

一束光

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

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

卡开发发

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

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

喵星大佬

我觉得你应该是很多常识没搞清楚

一束光

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

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

一束光

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

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

卡开发发

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

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

一束光

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

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

卡开发发

一束光 发表于 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的那种材料，这个可能需要你自己了解更多。

一束光

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

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

杨凤磊

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

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