1. DFT磁性计算工程导论

在普通非自旋密度泛函理论（DFT）计算中，电子密度 n(r) 是唯一的基本变量，体系能量对自旋不敏感。然而，铁、钴、镍及其合金、氧化物、二维磁性材料等体系的基态与自旋取向密切相关，必须引入自旋自由度。自旋极化DFT（Spin-Polarized DFT，SP-DFT）把电子密度拆分为自旋向上 n_↑(r) 与自旋向下 n_↓(r) 两个分量，总密度 n = n_↑ + n_↓，自旋磁密度 m = n_↑ − n_↓。交换关联能 E_xc[n_↑, n_↓] 同时依赖两个分量，从而在能量泛函层面允许自旋劈裂。

工程上最常用的实现是局域自旋密度近似（LSDA）与自旋极化广义梯度近似（GGA-PBE、GGA-PBEsol）以及杂化泛函（HSE06）。LSDA 往往高估磁矩，GGA 对金属磁矩较可靠，HSE06 对氧化物带隙与局域磁矩更准确但成本高。计算中需要显式设置自旋初态，例如 VASP 中 ISPIN=2、MAGMOM 给定各原子初始磁矩，LORBIT=11 输出分波投影；Quantum ESPRESSO 中 nspin=2、starting_magnetization 设定初始自旋。自旋极化计算通常比非自旋计算贵约 1.5～2 倍，且可能收敛到亚稳磁态，因此初态设置与收敛判据同样重要。

磁性材料可按磁有序类型分为：铁磁（FM，如 bcc Fe、hcp Co、fcc Ni）、反铁磁（AFM，如 NiO、Cr、MnF₂）、亚铁磁（FiM，如 Fe₃O₄）、顺磁（PM，高温或强阻挫体系）以及螺旋/非共线磁结构（如 γ-Fe、MnSi）。按电子机制又可分为巡游磁性（Fe、Co、Ni，Stoner 判据 I·N(E_F) > 1）、局域磁性（稀土 4f、氧化物中过渡金属 3d 强关联）与二维范德华磁性（CrI₃、Fe₃GeTe₂）。不同类别对泛函、U 值、自旋轨道耦合（SOC）与超胞设置的要求差异很大。

在材料工程中，自旋计算承担三类定位：一是解释与预测磁矩、磁有序温度与磁各向异性；二是为自旋电子学器件（MTJ、自旋轨道转矩、拓扑磁结构）提供参数；三是与实验（SQUID、中子衍射、X 射线磁圆二色）互证。典型工作流为：结构弛豫（非自旋或低精度自旋）→ 自旋极化静态计算 → 磁态枚举与能量比较 → 磁各向异性（SOC）→ 交换耦合常数 J_ij 提取 → 蒙特卡洛或海森堡模型估计 T_C/T_N。每一步的输入、判据与输出都需明确，否则容易得到“数值收敛但物理错误”的结果。

# INCAR 关键行
ISPIN   = 2          # 开启自旋极化
MAGMOM  = 5.0 5.0    # bcc Fe 两个原子初始磁矩（μB）
ENCUT   = 520        # eV
EDIFF   = 1E-6       # eV
ISMEAR  = 1          # 金属体系 Methfessel-Paxton
SIGMA   = 0.05
LORBIT  = 11         # 输出分波投影磁矩
LMAXMIX = 4          # 含 d/f 电子时提高混合
# 若为 NiO：需加 LDAU=.TRUE., LDAUTYPE=2, LDAUU=5.0, LDAUJ=1.0