VASP计算磁各向异性：从SOC开关到MAE数值的全流程

L1₀有序相的FePt是已知磁晶各向异性能（MAE）最高的材料之一——实验报道的单轴各向异性常数K_u高达7×10⁶ J/m³，是Nd₂Fe₁₄B的两倍以上。这个变态级各向异性使得FePt成为下一代超高密度磁记录介质和微机电系统（MEMS）永磁体的候选材料。

但MAE这个量在DFT计算里是最难算准的几个物理量之一。原因很简单：MAE的大小通常在meV/atom量级，而DFT总能的收敛精度通常在meV量级。你要求两个总能之间的差值比单次总能计算的数值噪声还小——这不是测不准，而是信噪比问题。

MAE的物理定义与计算路径

磁晶各向异性能是体系总能量随磁矩方向变化的函数：

$$MAE=E(\text{hard axis})-E(\text{easy axis})$$

对于L1₀-FePt，易磁化轴是c轴（[001]方向），难磁化轴在ab面内（比如[100]方向）。MAE就是这两个方向的能量差——正值表示易轴在[001]，负值表示在面内。

VASP中做MAE计算的标准路径是三层嵌套：

第一层，做非磁（或collinear磁性，但不开SOC）的结构优化。FePt的c/a比接近0.96（略扁的四方晶胞），这个c/a比对MAE的符号和大小都敏感——如果结构没弛豫到底，MAE的符号都可能算错。

第二层，在优化好的结构上，做collinear磁性计算（ISYM=0，开LSORBIT=FALSE），确认磁矩方向和大小。Fe的磁矩约2.9 μ_B，Pt感应出约0.3 μ_B，整体与文献一致。

第三层，开启SOC（LSORBIT=TRUE，设置SAXIS指定磁矩方向），分别算[001]和[100]两个方向的总能，取差值。关键参数：LORBMOM=FALSE——不预设轨道磁矩，让SOC自然感应。

数值收敛的三重检验

MAE对k点网格极其敏感。在5×5×5的k网格下，L1₀-FePt的MAE算出来是1.8 meV/f.u.（公式单元，含1个Fe和1个Pt）。加密到11×11×11，MAE涨到2.7 meV/f.u.；进一步到15×15×15，稳定在2.8 meV/f.u.左右。

这个k点敏感性源自SOC矩阵元在布里渊区某些高对称点附近对k的快速变化。FePt的SOC强度主要来自Pt原子——Pt的5d轨道SOC劈裂约0.5-0.6 eV，在DFT能带中表现为一些能带交叉点附近的SOC带隙开合，而这些交叉点恰好落在布里渊区的高对称线上。

平面波截断能（ENCUT）的收敛也需要验证。在500 eV的默认值下MAE为2.5 meV/f.u.，提到700 eV后收敛到2.8 meV/f.u.。建议ENCUT不低于600 eV。

最终的收敛值2.8 meV/f.u.，换算成K_u = MAE × 原子数密度 = 2.8 meV × (2/30.6 ų) ≈ 6.5 MJ/m³，与实验值7 MJ/m³的偏差不到10%。

Force Theorem vs 自洽SOC

DFT社区做MAE有两种主流方法。一种是前面走的标准自洽SOC路径——两个磁矩方向分别做完全自洽的SOC计算，直接取总能差。另一种是Magnetic Force Theorem（MFT）——只做一次scalar-relativistic的自洽计算，然后把SOC作为微扰加到单电子能量上：

MAE≈∑i[εiSOC(hard)−εiSOC(easy)]MAE≈i∑​[εiSOC​(hard)−εiSOC​(easy)]

MFT避免了两次全自洽SOC计算的成本，对FePt来说可以节省约一半机时。两种方法在FePt体系中的差异约0.15 meV/f.u.，在实际精度需求下几乎等价。MFT的唯一风险是当SOC导致费米面附近的能带发生较大重排时（比如某些重费米子体系），微扰假设可能不成立——但FePt不在此列。

应变工程：c/a比对MAE的调控

FePt薄膜在生长中经常受到衬底的应变。做了一组c/a比扫描（从0.90到1.05），发现MAE对c/a比的依赖相当陡峭：c/a=0.94时MAE只有1.5 meV/f.u.，c/a=0.98时达到2.8 meV/f.u.，c/a超过1.0后急剧下降到接近零。

这个趋势的物理来源是：c/a比改变了Fe 3d轨道的晶体场劈裂。四方晶场下，d_xy和d_{x²-y²}的能级差随c/a变化，而SOC耦合的矩阵元恰好在这两个轨道之间。c/a偏离最优值时，这两个能级的间距增大，SOC的二阶微扰贡献减小，MAE自然跌落。

工程意义上的价值

2.8 meV/f.u.这个数本身不是新发现。真正有用的是计算框架的可迁移性——一旦在FePt上验证了SOC计算参数（k点15×15×15、ENCUT 700 eV、c/a弛豫到位），就可以把它搬到FePt基合金体系（比如掺Co、掺Ni、或者跟Pt-Mn形成多层膜）中做MAE的高通量筛选。

对于磁记录介质和永磁体材料设计来说，实验合成一种新组分可能需要数周，DFT跑一周可以做上百个候选组分。VASP计算磁各向异性这件事的终极价值，就在这里。