自旋密度是计算磁学材料时最重要的物理量之一，它描述了体系中自旋向上与自旋向下电子的差值分布。在VASP计算中，自旋密度分析能够揭示材料的磁有序、局域磁矩、磁各向异性等关键信息。
自旋密度
定义为自旋向上与自旋向下电子密度之差：
其中：
为自旋向上电子密度
为自旋向下电子密度
根据VASP手册和官方文档，自旋密度也可表达为：
其中
和i分别为自旋和能带索引，
是归一化的单粒子波函数，
表示对所有占据态求和。
对于自旋极化体系，VASP中电荷密度
和磁化（自旋）密度
定义为：
在非自旋极化情况下，
，此时自旋密度为零。
应用领域 |
分析内容 |
典型研究体系 |
|
铁磁材料 |
自旋极化强度 |
Fe、Co、Ni、Gd等 |
|
反铁磁材料 |
磁矩反平行排列 |
MnO、NiO、Cr等 |
|
半金属 |
自旋通道导电性差异 |
Heusler合金、CrO₂ |
|
单分子磁体 |
局域磁矩稳定性 |
Mn₁₂、Fe₈等 |
|
二维磁性材料 |
维度效应下的磁性 |
CrI₃、CrGeTe₃ |
|
稀磁半导体 |
掺杂诱导的磁性 |
(Ga,Mn)As、ZnO:Co |
|
催化磁性 |
磁性对催化活性影响 |
金属氧化物表面 |
# ============ 自旋极化计算核心参数 ============
SYSTEM = Spin Density Analysis
PREC = Accurate # 高精度计算
ENCUT = 520 # 平面波截断能
# ============ 自旋极化控制参数 ============
ISPIN = 2 # 开启自旋极化计算
# 1:非极化 2:极化
MAGMOM = 5.0*1.0 1.0*2.0 # 每个原子的初始磁矩
# 格式：N*value 或 逐个指定
# 默认值为1.0/原子
# ============ 自旋轨道耦合(可选) ============
LNONCOLLINEAR = .TRUE. # 非共线自旋计算
SAXIS = 0 0 1 # 自旋量化轴方向
LSORBIT = .TRUE. # 开启自旋轨道耦合
# ============ 电子步参数 ============
EDIFF = 1E-6 # 电子自洽精度
ISMEAR = -5 # 四面体方法（精确积分）
SIGMA = 0.01
NEDOS = 5000 # DOS点数
# ============ 输出控制 ============
LCHARG = .TRUE. # 保存CHGCAR（包含自旋信息）
LWAVE = .TRUE. # 保存WAVECAR
LORBIT = 11 # 输出局域磁矩到OUTCAR
参数值 |
含义 |
适用场景 |
|
ISPIN = 1 |
非自旋极化 |
非磁性材料、闭壳层体系 |
|
ISPIN = 2 |
自旋极化 |
磁性材料、开壳层体系 |
|
ISPIN = 3 |
非共线自旋 |
磁畴结构、自旋螺旋 |
重要提示：
对于自旋极化计算，必须设置ISPIN=2
MAGMOM参数仅在ISPIN=2时生效
默认初始磁矩为1.0μB/原子
# ============ 非共线自旋计算 ============
LNONCOLLINEAR = .TRUE. # 开启非共线自旋
ISPIN = 2 # 必须设置
LSORBIT = .TRUE. # 自旋轨道耦合
SAXIS = 0 0 1 # 自旋量化轴
# 初始磁矩格式（x,y,z分量）
MAGMOM = 1*1.0 1.0 0.0 # 第一个原子
注意：
非共线计算需要VASP 5.4.1及以上版本
需要读取WAVECAR作为初始波函数（ICHARG=11）
磁矩由矢量(x,y,z)描述而非标量
# ============ DFT+U自旋计算 ============
LDAU = .TRUE. # 开启LDA+U
LDAU_TYPE = 2 # Dudarev方法
LDAUTYPE = 2
LDAUL = 2 2 0 # d轨道U (过渡金属)
LDAUU = 4.0 4.0 0.0 # U值
LDAUJ = 1.0 1.0 0.0 # J值
LDAUPRINT = 0 # 输出级别
ISPIN = 2 # 自旋极化
对于具有局域d或f电子的体系（如过渡金属氧化物、稀土材料），LDA+U能显著改善自旋磁矩的计算精度。
文件 |
内容 |
用途 |
|
OUTCAR |
总磁矩、原子磁矩、自旋密度 |
读取磁矩信息 |
|
CHGCAR |
自旋分辨电荷密度 |
分析自旋分布 |
|
PROCAR |
投影态密度 |
自旋分辨能带 |
|
EIGENVAL |
本征值 |
自旋能级结构 |
# OUTCAR中磁矩相关输出示例
mag = -2.9947154231 mu_B (site 1)
mag = 0.1523412567 mu_B (site 2)
…
total magnetic moment = 5.0000000000 mu_B
# 每个原子的磁矩
atom ion tot-mag up dn
1 1 2.50 1.75 -0.75
2 2 -2.50 -0.75 1.75
在自旋极化计算中，CHGCAR包含以下信息：
第一部分： 自旋向上电荷密度
第二部分： 自旋向下电荷密度
对于非共线计算，需要查看CHGCAR_UP和CHGCAR_DN或使用专门工具解析。
文件结构示例：
CHGCAR
1.0
a1_x a1_y a1_z
…
元素符号
原子数
Direct
原子坐标
…
NX NY NZ # 网格点数
自旋向上密度值… # 第一部分
自旋向下密度值… # 第二部分
自旋密度的物理意义与基本原理：自旋密度的定义、总电荷密度与自旋密度的关系、自旋密度分析的主要应用领域
VASP自旋极化计算参数设置：INCAR关键参数详解、ISPIN参数说明、非共线自旋计算设置
自旋密度输出文件分析：输出文件类型与内容、从OUTCAR读取磁矩信息、CHGCAR文件中的自旋密度
下章将正式引入本次教程的核心—VASP电荷密度计算。我们将从差分电荷密度的物理意义与理论基础、VASP差分电荷密度计算工作流程、差分电荷密度可视化工具方面详细介绍VASP差分电荷密度计算。
