同一个晶体结构，为什么用不同 MAGMOM 起算会得到不同能量、不同总磁矩，甚至不同 PDOS？原因在于 VASP 求的是电子自洽解，而磁性体系的自旋分布可能有多个局域极小点。非磁态、铁磁态和反铁磁态不是三个随手命名的标签，而是三种不同的自旋约束和初始有序方式。同一个几何结构可以对应多个磁性解；比较之前，先要弄清楚自己让体系从哪一种自旋排列出发。
对过渡金属氧化物、硫化物、M-N4 单原子位点、缺陷表面和吸附中间体来说，这一点尤其重要。一个结构在非磁计算里能收敛，不代表它的真实电子态就是非磁；一个铁磁初值能收敛，也不代表它一定是能量最低的磁性排列。若不把磁性解单独比较，后面的吸附能、Bader 电荷、局域磁矩和反应自由能都可能建立在不同电子态上。
非磁态不是“没有设置 MAGMOM”这么简单
非磁态通常对应 ISPIN = 1，体系不区分自旋向上和自旋向下电子。它适合闭壳层分子、很多主族半导体和明确没有局域磁矩的金属体系。但对含 Fe、Co、Ni、Mn、Cr、稀土元素或氧空位的体系，直接用非磁态可能把本来应该出现的局域自旋平均掉。此时得到的总能量和键合强度可能只是一个人为压平后的结果。
铁磁态则让多个磁性原子的初始磁矩同向排列，反铁磁态让相邻或指定子晶格磁矩反向排列。反铁磁态往往需要足够大的超胞才能表达，例如一个原胞只有一个金属位点，就无法写出相邻金属正负交替的排列。若模型尺寸不够，所谓“反铁磁比较”其实没有真正建立。
图1 MPX3 中 AFM-Neel、AFM-stripy、AFM-zigzag 自旋排列与晶体结构示意。DOI：10.1021/acs.jpcc.2c00646
还有一种常被忽略的亚铁磁态
如果不同子晶格的磁矩反向排列，但大小不完全抵消，总磁矩不为零，这类状态更接近亚铁磁态。许多掺杂、合金、缺陷和表面模型会自然出现这种情况：不是所有金属原子的磁矩都一样，也不是所有反向磁矩都刚好抵消。写结果时不能只看总磁矩为不为零，还要看每个磁性原子的局域磁矩方向和大小。
Mtot = Σi mi
这里的 Mtot 是体系总磁矩，mi 是第 i 个磁性原子的局域磁矩。铁磁态通常表现为多个 mi 同号；反铁磁态表现为相邻 mi 异号且总和接近零；亚铁磁态则是异号但不能完全抵消。这个区分会影响后面对能量排序和反应活性的解释。
交换作用和配位环境会竞争
磁性来自未成对电子，自旋之间的相互作用可以偏向同向，也可以偏向反向。金属-金属直接交换、金属-配体-金属超交换、配体场分裂和 d 轨道占据共同决定哪种排列更低能。改变一个键角、一个配体、一个吸附物或一个缺陷，就可能让铁磁和反铁磁排序发生切换。磁性解不是只由元素决定，而是由局域配位和电子占据共同决定。
图2 NiPS3 中直接反铁磁交换与铁磁超交换过程示意。DOI：10.1021/acs.jpcc.2c00646
DFT+U 会让这种差异更明显。U 值增强局域 d/f 电子的排斥，可能让电子更局域、磁矩更清楚，也可能把不同磁性态之间的能量差拉大。若一篇工作比较了多个 U 值，却只保留一个磁性初值，得到的趋势未必完整。更合适的做法是固定 U 值后先比较磁性排列，再在必要时讨论 U 值对排序的影响。
MAGMOM 给的是入口，不是最终答案
MAGMOM 主要决定初始自旋分布，并影响对称性处理。它不是强行规定最终磁矩必须等于输入值。比如 Fe 初始给 4 μB，自洽后可能变成 3.5 μB；给 -4 μB，最终也可能翻转或被相邻原子部分抵消。真正要比较的是收敛后的总能、总磁矩、各原子局域磁矩和电子结构，而不是只看 INCAR 里写了什么。
SYSTEM = magnetic_ordering_check
ENCUT = 520
EDIFF = 1E-6
ISPIN = 2
IBRION = -1
NSW = 0
LORBIT = 11
LCHARG = .TRUE.
LASPH = .TRUE.
LREAL = .FALSE.
# FM: MAGMOM = 4*4.0 32*0
# AFM: MAGMOM = 2*4.0 2*-4.0 32*0
这段片段对应同一几何结构下的磁性单点比较。关键不是参数多，而是 FM 和 AFM 任务除 MAGMOM 排列外保持一致。若要让 AFM 排列不被晶体对称性平均掉，还应检查结构是否允许正负磁性原子成为不同位置，必要时使用更大的超胞或关闭过高对称性。
先比同一几何，再决定是否分别优化
最清楚的比较顺序，是先在同一几何结构上做非磁、铁磁、反铁磁和可能的亚铁磁单点，判断电子态能量差。若某几种磁性态能量相差很小，再分别做结构优化，查看键长、层间距、吸附构型和磁矩是否随磁性态改变。这样可以区分“电子态差异”与“几何弛豫差异”。
ΔEmag = (Estate – Elowest) / Nmag
这里 Nmag 是磁性原子数。把能量差除以磁性原子数，可以避免大超胞和小超胞直接比较时产生误读。若 ΔEmag 只有几个 meV，温度、U 值、k 点和收敛精度都可能改变排序；若差到几十或上百 meV，磁性基态判断就更有把握。
图3 掺杂 MPX3 中 FM、AFM-N、AFM-s、AFM-z 不同磁性构型的能量排序。DOI：10.1021/acs.jpcc.2c00646
电子结构也要跟着磁性态一起读
不同磁性解可能有不同带隙、不同费米能级附近 d 态占据，也可能改变吸附物 p 态与金属 d 态的杂化。若只用最低能磁性态做总能，再用另一种磁性态的 PDOS 做解释，就会把电子结构和能量口径拆开。磁性体系的 DOS、Bader、差分电荷密度和吸附能最好来自同一轮自洽结果。
图4 二维亚铁磁合金的能带、PDOS 与局域自旋构型示意。DOI：10.1021/acs.jpcc.2c00646
对催化体系，还要注意吸附物本身可能改变磁性排序。干净表面是 AFM，吸附 *O 后可能局域磁矩下降；干净单原子位点是高自旋，吸附 CO 后可能进入低自旋。此时吸附能应写明使用哪一种磁性态作为反应起点。
Eads = Eslab+X,state – Eslab,state – EX
不要把“收敛到一个解”写成“只有一个解”
一个磁性初值收敛了，只说明从这个入口可以到达一个自洽解。更完整的写法，是说明比较了 NM、FM、AFM 或若干 AFM 排列，最低能态是哪一个，能量差是多少，最终总磁矩和局域磁矩如何分布。若不同磁性态能量接近，应把它们作为可能竞争的电子态保留讨论。
图5 多种 FM、AFM、FiM 排列的能量排序对比，说明磁性构型需要单独比较。DOI：10.1073/pnas.2422973122
用于吸附和催化时，推荐把结论写成两层：先说干净表面的磁性基态，再说关键吸附态是否保持同一磁性排序。若吸附前后磁性改变明显，就不要把吸附能变化全部归因于电荷转移或成键增强，还应说明自旋态重排参与了能量变化。
要点一：非磁态、铁磁态、反铁磁态对应不同自旋处理方式，不是同一个结果的不同命名。
要点二：同一几何结构可能有多个自洽磁性解，MAGMOM 决定入口，最终判断要看收敛后的能量和磁矩。
要点三：反铁磁和亚铁磁常需要足够大的超胞表达，原胞不够时不能强行比较。
要点四：吸附能、PDOS、Bader 和差分电荷密度应来自同一磁性态，避免电子态口径错位。
要点五：若多个磁性态能量接近，应把它们作为竞争电子态讨论，而不是只保留最顺眼的一组结果。
