计算化学公社

标题: qe计算有铝掺杂的锰酸锂晶体单点能不收敛 [打印本页]

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

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

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

作者
Author: WilliamH 时间: 2023-5-16 12:54

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

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

作者
Author: likejade 时间: 2023-5-16 16:28
能量毫无下降收敛的趋势，还需要做很多测试，除了k点、rh0等之外，磁性、mixing_beta等的设置对收敛也有影响的……

作者
Author: Quantum198907 时间: 2023-6-3 18:52
（1）自旋极化没有考虑
（2）试试DFU+U，可能还需要合理设置一个B值，手册里面有说明
（3）mixing_beta设置成0.2试试看

作者
Author: yjw1204338757 时间: 2023-6-5 16:22

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

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

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