计算化学公社

标题: Fe2O3 计算DOS 对称性问题请教 [打印本页]

作者
Author: wprison 时间: 2021-12-22 09:42
标题: Fe2O3 计算DOS 对称性问题请教
本帖最后由 wprison 于 2021-12-22 10:33 编辑

请教各位老师，VASP计算Fe2O3的DOS图时，文献中给出的结果都呈现自旋对称Total DOS，例如
(, 下载次数 Times of downloads: 17)
但是自己计算的时候，为什么结果显示出不对称？是INCAR哪个参数设置问题吗？谢谢！
(, 下载次数 Times of downloads: 21)
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

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

作者
Author: wprison 时间: 2021-12-22 23:23

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

谢谢老师，按照这个设置MAGMOM，DOS就是自旋对称的了。十分感谢！ (, 下载次数 Times of downloads: 25) (, 下载次数 Times of downloads: 25)

作者
Author: 八十番 时间: 2022-6-10 18:38

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

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

作者
Author: sigma-jlu 时间: 2023-4-23 19:19

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

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

作者
Author: sigma-jlu 时间: 2023-4-23 19:25

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

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

作者
Author: sigma-jlu 时间: 2023-4-23 19:28

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

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

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