DFT+U的自洽结果与初始的 occupation matrix有关，如果不设，很容易收敛到亚稳态，文献[1]中对算UO2有详细的讨论。文献[2][3]讨论了怎么设置occupation matrix能算出低能态。
Zunger 18年计算NiO, CoO的文章[4]引了以上的文献，并按照以上文献设置occupation matrix，以防算出亚稳态。
我自己算LixCoO2的时候，LMAXMIX设不设，自洽后的局域磁矩都会变化，虽然VASP wiki上说LMAXMIX参数只是加快收敛。甚至按说明书d电子体系设LMAXMIX=4，反而自洽出一个亚稳态，不设LMAXMIX时，自洽结果的能量却更低，应该也是跟初始占据有关。

1. PHYSICAL REVIEW B 79, 235125 (2009)
2. PHYSICAL REVIEW B 82, 195128 (2010)
3. Phys. Chem. Chem. Phys., 2014, 16, 21016—21031
4. PHYSICAL REVIEW B 97, 035107 (2018)

希望大佬解答