Fe2O3 计算DOS 对称性问题请教

wprison

请教各位老师，VASP计算Fe2O3的DOS图时，文献中给出的结果都呈现自旋对称Total DOS，例如

但是自己计算的时候，为什么结果显示出不对称？是INCAR哪个参数设置问题吗？谢谢！

INCAR
DFT+U Calculation
LDAU   = .TRUE.        (Activate DFT+U)
LDATYPE=  2            (Dudarev; only U-J matters)
LDAUL  =  2 -1         (Orbitals for each species)
LDAUU  =  5.3  0         (U for each species)
LDAUJ  =  0  0         (J for each species)
LMAXMIX=  4            (Mixing cut-off; 4-d, 6-f)

Global Parameters
ISTART =  1           (Read existing wavefunction; if there)
ISPIN =  2           (Spin polarised DFT)
ICHARG =  11         (Non-self-consistent: GGA/LDA band structures)
LREAL  = Auto          (Projection operators: automatic)
ENCUT  =  520      (Cut-off energy for plane wave basis set, in eV)
PREC   =  Normal       (Precision level)
LWAVE  = .TRUE.        (Write WAVECAR or not)
LCHARG = .TRUE.        (Write CHGCAR or not)

ISIF   =  3            (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions)
Electronic Relaxation
ISMEAR =  0            (Gaussian smearing; metals:1)
SIGMA  =  0.05         (Smearing value in eV; metals:0.2)
NELM   =  60           (Max electronic SCF steps)
NELMIN =  6            (Min electronic SCF steps)
EDIFF  =  1E-05        (SCF energy convergence; in eV)
# GGA  =  PS           (PBEsol exchange-correlation)

Ionic Relaxation
NSW    =  0          (Max electronic SCF steps)
IBRION =  -1            (Algorithm: 0-MD; 1-Quasi-New; 2-CG)
EDIFFG = -0.03        (Ionic convergence; eV/AA)
# ISM =  2            (Symmetry: 0=none; 2=GGA; 3=hybrids)

MAGMOM = 12*5.0 18*0.6

卡开发发

Fe2O3应该是AFM的磁序，MAGMOM你应该重新调整，有些Fe是正的有些是负的，不应该全是正的。

wprison

卡开发发 发表于 2021-12-22 16:50
Fe2O3应该是AFM的磁序，MAGMOM你应该重新调整，有些Fe是正的有些是负的，不应该全是正的。

谢谢老师，按照这个设置MAGMOM，DOS就是自旋对称的了。十分感谢！

八十番

wprison 发表于 2021-12-22 23:23
谢谢老师，按照这个设置MAGMOM，DOS就是自旋对称的了。十分感谢！

你这个原子的磁矩分布图是在哪里找的

sigma-jlu

八十番 发表于 2022-6-10 18:38
你这个原子的磁矩分布图是在哪里找的

PhysRevB.96.094426论文给的Morin相变前后的反铁磁组态图比较清晰。

sigma-jlu

wprison 发表于 2021-12-22 23:23
谢谢老师，按照这个设置MAGMOM，DOS就是自旋对称的了。十分感谢！

哦，这个赤铁矿低温单胞内反铁磁示意图的出处： https://www.mdpi.com/2075-163X/4/1/89

sigma-jlu

sigma-jlu 发表于 2023-4-23 19:25
哦，这个赤铁矿低温单胞内反铁磁示意图的出处： https://www.mdpi.com/2075-163X/4/1/89

不过这个图用于吸附，恐怕不现实，因为这中c向磁矩一般在Morin转变，即250K以下才存在，室温反铁磁几乎是ab面内取向。