【求助】vasp 计算磁矩结果比别人的要大
最近看了两篇关于Co2FeAl0.5Si0.5的文献（J. Phys: Condens. Matter 19(2007) 326216; J. Phys. D: Appl. Phys. 40(2007)1582-1586)，觉得比较有兴趣，便做了一下Co2FeAl0.5Si0.5这种四元合金的计算。结果我发现我得到磁矩为5.75，而这两篇文献的值均为5.5，经过多次计算我得到值都是5.75，我想请教大家究竟是怎么回事。有知道的请指点一下，谢谢！
以下是我的incar文件：
Startparameter for this Run:
SYSTEM = L21 Co2FeAlSi-static-calculation
ENCUT = 500
NWRITE = 2 write-flag & timer
PREC = Accurate
ISTART = 0 job : 0-new 1-cont 2-samecut
ICHARG = 2 charge: 1-file 2-atom 10-const
ISPIN = 2 spin polarized calculation?
MAGMOM = 8*3.0 4*4.0 2*0.0 2*0.0
INIWAV = 1 electr: 0-lowe 1-rand 2-diag
NELM = 60
NELMIN = 6; NELMDL= -5 # of ELM steps
EDIFF = 1E-04 stopping-criterion for ELM
VOSKOWN = 1
NPAR = 1
LORBIT = 11
RWIGS = 1.16 1.17 1.18 1.11
LDAU = .TRUE.
LDAUTYPE = 4
LDAUL = 2 2 1 1
LDAUU = 4.453 4.185 0 0
LDAUJ = 0.655 0.625 0 0
Electronic Relaxation 2
IALGO = 38 algorithm
LDIAG = T sub-space diagonalisation
LREAL = Auto real-space projection
ROPT = 2E-4 2E-4 2E-4 2E-4
ALGO = Fast
DOS related values:
ISMEAR = -5
SIGMA = 0.1 broadening in eV -4-tet -1-fermi 0
NEDOS = 1000
Ionic Relaxation
EDIFFG = 1E-03 stopping-criterion for IOM
NSW = 0 number of steps for IOM
ISYM = 2 0-nonsym 1-usesym 2-use
返回小木虫查看更多
京公网安备 11010802022153号
这个要看你用的参数和文献的是否一致
不好意思，我也不知道为什么，我想请教下，就是这里面关于Co和Fe的U和J值你是怎么得到的,谢谢
看你用的是什么泛函，什么赝势，还有 ISMEAR SIGMA 这个对磁矩影响很大
你得总磁矩从OUTCAR 还是 OSZICAR里得到的？
一般看文献，别人用多少就用多少
，
个人觉得可能是你的赝势选取的问题，你可以换一下赝势试试。。。
5.5和5.75差别很大吗？(5.75-5.5)/5.5=4.5%
各位大哥，我的U和J的值是从文献中找到的。这篇文献中算的J和U的值是用wien2k算的。我计算的时候选的是GGA＋U，和这两篇文献的一致。我采用这两篇文献中的U和J的值的时候，能够得到5.5的磁矩，但少数自旋的带隙（该体系为半金属）只有0.34eV，比别人的带隙（约1ev）要小很多，各位有很什么好建议么？