【界面体系电子结构难以收敛问题】

Alexander

计算一个包含大概20 A真空层的体系的电子自洽，发现用pulay电荷混合方法，降低混合参数beta到0.01，增加混合维度到30，体系电荷密度依然很难收敛到1e-06，在4e-05附近反复震荡，但是此时体系能量差已经收敛到1e-10. 请问是什么回事呢？初步有以下猜测：
1、目前界面体系不好收敛，需要优化算法；
2、可以增加能量阈值判断收敛条件，当能量差值小于1e-10（或者其它）也可以判定自洽结束；
3、降低电荷密度收敛标准到4e-05，但是不能保证每次能量收敛都非常好。

附件是我的计算输入参数以及log。

zdy-AISI

smearing_sigma                 0.001 #energy range for smearing
out_bandgap                    1

Alexander

尝试了一下你给的参数，确实可以收敛，谢谢。