求助QE计算SOC的FeGe自洽无法收敛

zss2022

首先，感谢各位老师百忙中点进来关注这个问题！请教各位老师，我用QE计算加了SOC的FeGe始终无法收敛能量，里面是不是有啥微妙之处啊？感觉QE的SOC经常很难收敛。。。scf.in贴在下方了。
#-------------------------------------------------------------------------------------------------------
&control
                calculation  = 'scf'
                    prefix = 'FeGe',
                 pseudo_dir = 'pp/',
                    outdir = './tmp'
                  verbosity = 'high',
/
&system
                    ibrav  = 0,
                       nat = 8,
                      ntyp = 2,  
                   ecutwfc = 80,   
                   ecutrho = 600,               
              occupations  = 'smearing',
                 smearing  = 'mp',
                   degauss = 0.01
                lspinorb   = .true.
                noncolin   = .true.
starting_magnetization(1) = 1.0
              angle1(1)    = 54.7356
              angle2(1)    = 45.0
/
&electrons
                  conv_thr  = 1.0d-7
/

ATOMIC_SPECIES
    Fe 55.84500  Fe.rel-pbe-spn-kjpaw_psl.1.0.0.UPF
    Ge 72.64     ge_pbe_v1.4.uspp.F.UPF

ATOMIC_POSITIONS (crystal)
Fe      0.1352200059999973  0.1352200059999973  0.1352200059999973
Fe      0.3647800090000004  0.8647800090000004  0.6352199909999996
Fe      0.8647800090000004  0.6352199909999996  0.3647800090000004
Fe      0.6352199909999996  0.3647800090000004  0.8647800090000004
Ge      0.8418599959999966  0.8418599959999966  0.8418599959999966
Ge      0.6581400040000034  0.1581400040000034  0.3418599959999966
Ge      0.1581400040000034  0.3418599959999966  0.6581400040000034
Ge      0.3418599959999966  0.6581400040000034  0.1581400040000034


K_POINTS {automatic}
    9 9 9 0 0 0

CELL_PARAMETERS {angstrom}
     4.6995000839000003   -0.0000000000000000    0.0000000000000000
    -0.0000000000000000    4.6995000839000003   -0.0000000000000000
    -0.0000000000000000    0.0000000000000000    4.6995000839000003

HUBBARD (ortho-atomic)
    U Fe-3d 2.0

Exodus

初始磁矩确定全是spins up吗？

zss2022

Exodus 发表于 2024-3-11 15:03
初始磁矩确定全是spins up吗？

应该是吧，FeGe是常见的FM金属。而且VASP这样计算是没问题的。

zss2022

唉，实在搞不定，不知道应该调啥。尝试了增大K点，能量截断，重新优化了结构，改变occupation，degauss，换了rel-pbe赝势，都无法收敛。