计算化学公社

标题: 求助QE计算SOC的FeGe自洽无法收敛 [打印本页]

作者
Author: zss2022 时间: 2024-3-11 11:35
标题: 求助QE计算SOC的FeGe自洽无法收敛
首先，感谢各位老师百忙中点进来关注这个问题！请教各位老师，我用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

作者
Author: Exodus 时间: 2024-3-11 15:03
初始磁矩确定全是spins up吗？

作者
Author: zss2022 时间: 2024-3-11 15:45

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

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

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

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