计算化学公社

标题: VASP计算稀土掺杂磁性体系时磁矩该如何设置？ [打印本页]

作者
Author: CalebParadoX 时间: 2025-4-27 10:00
标题: VASP计算稀土掺杂磁性体系时磁矩该如何设置？
1，我是想做稀土元素掺杂Ti2CO2，两个La原子替代两个Ti原子形成我这么一个结构，根据文献我想将两个掺杂原子设为FM和AFM构型，这个时候我该如何设置磁矩，以及这个时候我是否要加自选轨道耦合和DFT+U
这是我之前的INCAR，我尝试加了磁矩，但是计算以后OUTCAR显示磁矩为0:
Global Parameters
ISTART =  0          (Read existing wavefunction; if there)
NPAR=2
ISPIN =  2          (Spin polarised DFT)
ICHARG =  2          (Non-self-consistent: GGA/LDA band structures)
LREAL  = Auto       (Projection operators: automatic)
ENCUT  =  410       (Cut-off energy for plane wave basis set, in eV)
PREC =  Fine          (Precision level)
LWAVE  = .FALSE.       (Write WAVECAR or not)
LCHARG = .TRUE.       (Write CHGCAR or not)
ADDGRID= .TRUE.       (Increase grid; helps GGA convergence)
ISYM = 0
MAGMOM = +4 +4 45*1

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  =  0.0001    (SCF energy convergence; in eV)
GGA  =  PE          (PBE exchange-correlation)

Ionic Relaxation
NSW =  200       (Max electronic SCF steps)
IBRION =  2          (Algorithm: 0-MD; 1-Quasi-New; 2-CG)
EDIFFG = -0.02       (Ionic convergence; eV/AA)
VOSKOWN = 1
IVDW = 11
LAECHG=.TRUE.
LORBIT = 11
原子顺序方面，两个掺杂原子顺序是前两个

作者
Author: Stars 时间: 2025-5-9 10:40
vasp的赝势对稀土元素的处理有大核和小核的版本，要根据价层电子数量判定。

作者
Author: wqYang 时间: 2025-10-1 17:41

Stars 发表于 2025-5-9 10:40
vasp的赝势对稀土元素的处理有大核和小核的版本，要根据价层电子数量判定。

您好老师，请问如果是选取的冻结4f轨道的赝势，算吸附能是否需要开磁矩呢？

作者
Author: Stars 时间: 2025-10-13 20:59

wqYang 发表于 2025-10-1 17:41
您好老师，请问如果是选取的冻结4f轨道的赝势，算吸附能是否需要开磁矩呢？

看体系内是否有这样的有磁矩的原子，如有都应当打开

作者
Author: wqYang 时间: 2025-10-16 19:01

Stars 发表于 2025-10-13 20:59
看体系内是否有这样的有磁矩的原子，如有都应当打开

感谢解答！

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