qe计算有铝掺杂的锰酸锂晶体单点能不收敛

yjw1204338757 · 发表于 Post on 2023-5-16 09:09:30

请问一下各位老师，我做的是有铝掺杂的锰酸锂晶体单点能的计算，在前面k点、rh0测试和vc-relax都能正常完成，就是在取vc-relax最终构型跑单点能的时候跑了200步一直不收敛，不知道是不是晶格改变后更难收敛了还是有其它问题。

WilliamH · 发表于 Post on 2023-5-16 10:12:32

建议你在计算之前先看一下赝势文件开头的说明，里面包含有建议的ecutwfc和ecutrho，当前ecutwfc至少对于Li元素来说是不够的。计算锰酸锂，需要考虑自旋极化，并对磁矩进行合理初猜。

yjw1204338757 · 发表于 Post on 2023-5-16 11:07:24

老师您好，里面ecutwfc和ecutrho是经过rh0测试过的值，您说的ecutwfc不够是指的太小了吗？如果太小的话我给他再放大一点，自旋极化和磁矩是没有考虑到，刚接触量化不久，对这个还不是很了解

WilliamH · 发表于 Post on 2023-5-16 12:54:29

yjw1204338757 发表于 2023-5-16 11:07
老师您好，里面ecutwfc和ecutrho是经过rh0测试过的值，您说的ecutwfc不够是指的太小了吗？如果太小的话我给 ...

不太懂你说的“经过rh0测试”是什么意思。你所用的赝势文件中对Li元素的ecutwfc建议值是103Ry，你可能需要再测试一下，或者更换赝势文件。

likejade · 发表于 Post on 2023-5-16 16:28:39

能量毫无下降收敛的趋势，还需要做很多测试，除了k点、rh0等之外，磁性、mixing_beta等的设置对收敛也有影响的……

Quantum198907 · 发表于 Post on 2023-6-3 18:52:44

（1）自旋极化没有考虑
（2）试试DFU+U，可能还需要合理设置一个B值，手册里面有说明
（3）mixing_beta设置成0.2试试看

yjw1204338757 · 发表于 Post on 2023-6-5 16:22:33

Quantum198907 发表于 2023-6-3 18:52
（1）自旋极化没有考虑
（2）试试DFU+U，可能还需要合理设置一个B值，手册里面有说明
（3）mixing_beta设 ...

老师您好，这个考虑自旋极化是在in文件的哪里能体现出来呢，有没有qe手册可以分享一下呢，初次接触这个还有很多没有搞明白

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Quantum ESPRESSO] qe计算有铝掺杂的锰酸锂晶体单点能不收敛