在PWSCF计算LDA+U，随着U的值增大Mn 原子磁矩增大，文献中是随着U增大磁矩减小的。我的输入文件下面那样，我实在是找的不到错误，请问各位大神帮帮我看看，问题出现在哪里。感觉自己对对的，结果却是错错的。
&system
ibrav= 0,
nat= 8, ntyp= 2,
ecutwfc = 100.0,
ecutrho = 600,
occupations = 'smearing',
smearing = 'gaussian',
degauss = 1.d-8,
lda_plus_u=.true. ,
Hubbard_U(1)= 5.5,
nspin = 2,
starting_magnetization(1) = 0.4
/
&electrons
conv_thr = 1.d-10
mixing_fixed_ns = 0
/
ATOMIC_SPECIES
Mn 54.938 Mn.pz-spn-kjpaw_psl.0.3.1.UPF
Si 28.086 Si.pz-n-kjpaw.UPF

文献发来看看，还有检查一下结构是不是一致。增加U让Mn的d电子更加局域，你计算的材料当中是什么交换作用？比如RKKY？增加U磁矩减小的原因是什么？

楼主提到的两篇文献（PRB93.195101和PhysicaB430，1420，使用的软件Wien2k和VASP），如果想要重复这两篇文章的结果，在QE中是可以+U的，首先+U对于材料做晶格结构驰豫，ecutwfc=100已经比较高，ecutrho=800也就是ecutwfc的8倍，赝势为什么用paw，用rrkjus支持+U吗。楼主说的lda 0.8，pbe 1.0这种差别考虑了晶格驰豫吗。