自旋量子数的计算：从第一性原理到磁性材料的自旋极化处理

自旋量子数的计算，是磁性材料DFT研究中最早遇到、也最容易埋坑的地方。很多用户在VASP里直接开ISPIN=2就以为搞定了自旋极化——实际上，初始磁矩的设定（MAGMOM）是否合理，直接决定了你最终收敛到的是物理上正确的自旋态，还是一个能量更低的、但磁学性质完全错误的电子态。

MAGMOM不是随便填的数字

VASP的MAGMOM标签决定了每个原子初始的局域磁矩方向。对于过渡金属氧化物（比如Fe2O3、Co3O4），如果不按离子轮的磁性状态来设初始值，计算经常会掉进自旋污染的陷阱里。

一个实际的经验：对Fe³⁺（高自旋，S=5/2），每个Fe原子设MAGMOM=5是比较合理的起点。但如果体系里有反铁磁耦合，还必须用LSORBIT=.TRUE.加上SAXIS来指定磁轴方向，否则VASP默认把所有磁矩沿z轴排列，直接把反铁磁基态给忽略了。

自旋多重度与总自旋的对应关系

自旋量子数S和自旋多重度（2S+1）的对应，是解读VASP磁学输出的基础。OUTCAR里的magetization (x)给出的是每个原子的局域磁矩（μ_eff，单位μ_B），而总自旋S需要通过公式 μ_eff = g√[S(S+1)] 来反推（g≈2 for 3d体系）。

实践中，我习惯在INCAR里设置LORBIT=11，让PROCAR和DOSCAR输出每个轨道的磁矩分解。这样不仅能看到总自旋，还能判断磁矩主要来自d_x²-y²还是d_xy轨道——这个信息对后续分析磁各向异性和交换耦合常数J是不可或缺的。

DFT+U框架下的自旋修正

对强关联体系（NiO、MnO、La2CuO4等），普通GGA-PBE会系统性低估带隙和磁矩。这时候必须加U（在INCAR里设LDAU=.TRUE. + LDAUTYPE=2 + LDAUU值）。

但U值的选取不能拍脑袋。一个可行的方案是：先做一组不同U值的扫描计算（比如U=2, 3, 4, 5 eV），把得到的磁矩和带隙与实验值对比，取偏差最小的U。更准确但更费力的做法是线性响应方法（Cococcioni & de Gironcoli, PRB 2005），VASP用户可以用LDAUTYPE=3来自动算U的自洽值。

非共线自旋与自旋轨道耦合（SOC）

当体系有重元素（Au、Pt、Bi等）或者研究磁各向异性（MAE）时，必须打开自旋轨道耦合。INCAR里对应的设置是LSORBIT=.TRUE. + SAXIS= 0 0 1（磁轴沿z方向）。

SOC的计算量通常是非SOC计算的3-5倍，因为PAW投影算子要处理自旋oribtal耦合项。一个实际的折中方案：先不加SOC做结构优化和自洽，得到收敛的波函数后，再在ICHARG=11下做SOC的静态计算。这样能节省大约40%的计算时间，同时对MAE的计算精度损失很小（<0.5 meV/atom）。

自旋量子数的计算不是”设个ISPIN=2″就能交差的。初始磁矩的合理设定、U值的系统校准、SOC的按需引入，这三步做好了，DFT给出的磁学性质才能和实验对上。