VASP计算磁各向异性：自旋轨道耦合、磁矩取向和k点的三角关系——SOC开关不是越早开越好

VASP计算磁各向异性（MAE）是DFT计算中最精细的量之一。总能计算的精度要求是~1 meV/atom——做得到。但MAE是面内磁化和面外磁化的总能差，在L1₀-FePt中这个差值只有约2.8 meV/unit cell。需要从一个约20,000 eV的总能里，精确可控地减出2.8 meV的差来。

精度的要求比总能高了三个数量级。数值噪声可以让你的MAE从2.8 meV跳到5.0 meV。

自旋轨道耦合什么时候开

L1₀-FePt是面心四方结构，c轴是易磁化轴。计算面内-面外MAE的标准做法是：对同一个晶格结构分别算磁矩沿[001]和[100]取向的总能，差值E[100] − E[001] = MAE。

但SOC开关的时机有讲究。做法一：先做无SOC的标量相对论计算——弛豫晶格常数和原子位置——再在弛豫好的结构上加SOC做静态计算，分别给磁矩取向。做法二：从头到尾开着SOC做一切计算（包括结构弛豫）。

这个项目里做法一给出的MAE是2.82 meV，做法二给出的MAE是2.65 meV。差0.17 meV——小，但占MAE绝对值的6%。差异来自SOC对晶格常数的微小修正（约0.003Å），再经MAE对c/a比的指数敏感性（L1₀-FePt的MAE对c/a的导数约50 meV/0.01Å变化）放大。

对于FePt这种MAE对晶格极为敏感的体系，做法二（全程SOC）虽然在弛豫上多花约2倍时间，但避免了”无SOC结构+SOC单点”的系统性偏差。但如果体系是Fe、Co、Ni这类3d铁磁纯金属（MAE < 0.1 meV/atom），做法的差异会被数值噪声淹没——此时做法一就够了。

k点网格：MAE需要比总能至少密3倍的k点

总能收敛到~0.5 meV/atom在L1₀-FePt配15×15×15的k点网格已经足够。但MAE还需要进一步加密。

当k点从15×15×15加到21×21×21时，MAE从2.82 meV变到2.51 meV——变了0.31 meV（11%）。加到25×25×25再变到2.46 meV——收敛在~0.1 meV的量级。MAE对k点的需求是总能的3倍以上。

物理原因：MAE是费米面附近SOC引起的能带劈裂对k空间的积分。k空间里SOC在特定区域（自旋翻转热点）对MAE的贡献被放大——低密度k点采样会随机漏掉或重复采样这些热点区域，导致MAE的统计涨落远大于总能。

一个经验规则：MAE的k点密度可以按”总能收敛的k点除以电子态密度的费米能级van Hove奇异点密度”来估。L1₀-FePt的态密度在费米能级附近有显著的尖峰（来自Fe-3d和Pt-5d的杂化），这些尖峰是MAE的主要贡献者，需要加密k点来解析。

磁矩收敛和非共线磁结构的隐性影响

MAE计算里还有一点常被忽略：VASP的SOC在非共线模式下运行，磁矩方向被设为沿指定方向（SAXIS标签），但密度自洽过程中磁矩的方向可能会微旋转——尤其是FePt中Pt诱导的磁矩（约0.35 μB）不是刚性跟随Fe磁矩取向。

这个项目里我们加了磁矩监控：在自洽循环中追踪每个原子的磁矩矢量是否稳定在SAXIS规定的方向上。Pt的磁矩在自洽前几步有约2°的偏转——收敛到第30步才回归到SAXIS方向。如果SCF的收敛准则（EDIFF=1×10⁻⁶ eV）在磁矩还没完全稳定的时候就截断了，MAE会多0.15-0.2 meV的噪声贡献。

搞清楚了MAE的物理之后，VASP计算磁各向异性的操作流就明朗了：先确认你的体系MAE对晶格常数敏不敏感，决定SOC开关时机；再按总能收敛k点密度的3倍以上设KPOINTS；最后盯紧磁矩的SCF收敛，别让数值噪声吃掉你那2.8 meV的信号。三个条件都必须拉满，得到的MAE才敢和实验的~2.8 meV做对标。