VASP磁矩计算：自旋极化设置与磁结构的正确处理

磁性材料的计算在VASP中有一套特有的参数逻辑，不同于一般的非磁体系。很多人第一次做磁性体系的计算，要么忘了开自旋极化，要么开了却给了不合理的初始磁矩，导致计算收敛到错误的磁态。问题出来了也不容易发现——能量数值看起来合理，磁矩输出却是零，OUTCAR里悄悄告诉你这个体系被当成非磁处理了。

ISPIN：自旋极化的开关

VASP中控制自旋极化的参数是ISPIN：

• ISPIN=1：非自旋极化计算（默认值），α和β自旋共用同一套轨道
• ISPIN=2：自旋极化计算，上下自旋独立处理

对于磁性体系（Fe、Co、Ni、Mn、Cr等过渡金属及其化合物），必须设置ISPIN=2。非磁体系（Cu、Al等）用ISPIN=1即可，强行设置ISPIN=2不会出错，但会让计算时间加倍而毫无意义。

判断一个体系是否需要开自旋极化，可以先查Materials Project数据库，确认该化合物是否有磁基态（磁矩非零）。对于新体系，最保险的做法是分别用ISPIN=1和ISPIN=2跑一遍，比较总能——如果ISPIN=2的能量明显更低（差0.05 eV/atom以上），说明磁性对该体系确实重要。

MAGMOM：初始磁矩的正确设置

ISPIN=2只是打开了自旋极化的可能性，初始磁矩MAGMOM决定了体系在哪个磁态上收敛。MAGMOM的格式是按POSCAR中的原子顺序给出每个原子的初始自旋磁矩，单位μB。

铁磁构型（所有磁性原子磁矩同向）：

MAGMOM = 4 4 4 4  # 4个Fe原子各2 μB

反铁磁构型（交替排列）：

MAGMOM = 4 -4 4 -4  # 正负交替

对于非磁性辅助原子（O、S等）可以设为0：

MAGMOM = 4 0 0 4 0 0  # Fe-O-O-Fe-O-O序列中O为0

初始磁矩的量级选取有一定经验规律：Fe通常给2-4 μB，Co和Ni给1-2 μB，Mn视化合物价态给2-5 μB，稀土元素可以给到7-8 μB。具体数值会在自洽过程中调整到收敛值，初始值的意义在于引导体系进入正确的磁态区域，而不是精确预测最终磁矩。

如果MAGMOM设为全0，体系极有可能收敛到非磁态——即使该体系真实存在磁序，DFT的自洽也可能停留在能量较高的非磁局域极小值而不自动跳出。

磁结构的选择：哪种磁序能量最低

对于未知磁结构的体系，需要比较多种磁序候选态的能量，选取最低能量的作为基态磁结构。通常比较的构型包括：铁磁（FM）、不同类型的反铁磁（AFM-I、AFM-II等）和亚铁磁。

确定磁结构时可以参考实验中子散射数据、文献磁性研究或同类化合物的规律。对于强关联体系（如氧化物、磁性半导体），不同磁序之间的能量差可能只有几meV/atom，这时DFT计算精度本身也接近这个量级，结论需要谨慎。

LORBIT：输出分原子磁矩

设置LORBIT=11（推荐）或LORBIT=10，VASP会在OUTCAR中输出每个原子的分波磁矩（s、p、d分量及总量）。从OUTCAR中查找磁矩信息：

grep "magnetization (x)" OUTCAR  # 磁矩方向（非共线需要）
grep "tot" OUTCAR | tail -n [原子数+3]  # 分原子磁矩表格

或直接用vaspkit → 512（磁矩提取）一步输出，格式更清晰。

对于铁磁体系，总磁矩（自旋向上电子数 – 自旋向下电子数）应为接近整数的值（对于整数磁矩的绝缘体），或非整数（对于金属）。如果总磁矩为0但ISPIN=2，说明体系收敛到了顺磁或非磁态，需要检查MAGMOM设置。

非共线磁性：LNONCOLLINEAR的特殊处理

对于自旋轨道耦合（SOC）主导的体系或磁学拓扑材料（如磁性外尔半金属、斯格明子系统），磁矩方向不能简单地用标量描述，需要开启非共线磁性计算：LNONCOLLINEAR=TRUE，并同时设置LSORBIT=TRUE（包含自旋轨道耦合）。

非共线计算中MAGMOM的格式变为每个原子的三维矢量（MAGMOM = 0 0 4 0 0 4 …），计算时间相比共线磁性大幅增加（通常4-8倍），对体系的物理图像也需要更仔细的分析。

科研学术网提供磁性材料全套计算服务，包括磁矩、磁各向异性和交换耦合常数，欢迎访问 https://www.keyanxueshu.com 了解详情。