4. 自旋计算关键参数设置
自旋极化计算与普通非自旋计算的最大区别，在于电子密度被拆分为自旋向上与自旋向下两套通道，体系能量对磁矩构型产生依赖。因此，参数设置不只是“算得准不准”的问题，更直接决定计算能否收敛到正确的磁态。本节围绕六个关键参数展开：自旋极化开关、初始磁矩、k点网格、截断能、展宽方法、收敛判据，并重点讨论磁性计算特有的收敛陷阱。
一、自旋极化开关：先确认是否真的打开了自旋通道
在 VASP 中，自旋极化由 ISPIN 控制：ISPIN=1 为非自旋极化，ISPIN=2 为自旋极化。对含 Fe、Co、Ni、Mn 等磁性元素的体系，必须使用 ISPIN=2。若忘记打开，即使初始磁矩设置得再合理，计算也会收敛到非磁态，得到的总能、态密度和磁矩全部错误。
打开自旋极化后，OUTCAR 中会出现自旋向上/向下两套本征值，且总磁矩由 MAGMOM 初始值出发，经自洽迭代逐步收敛。判断是否真正进入磁性计算，可检查 OUTCAR 中是否出现 ISPIN = 2，以及是否输出 magnetization (x) 等磁矩分量。
二、初始磁矩设定：给自洽迭代一个合理的起点
MAGMOM 是磁性计算中最容易被忽视、却最影响收敛结果的参数。它按 POSCAR 中原子顺序逐个给出初始磁矩，单位为 μB。例如 Fe 单胞两个原子，可写为：
MAGMOM = 2*5.0
对于 Fe、Co、Ni，常见初始值可取 3～5 μB；Mn 可更高，约 4～5 μB；若为反铁磁，则需按子晶格交替设置正负号，例如：
MAGMOM = 2.0 -2.0 2.0 -2.0
初始磁矩不是最终磁矩，但它决定自洽迭代从哪个磁态附近出发。若初始值过小，可能收敛到非磁或低自旋态；若符号设置错误，可能收敛到错误的磁序。因此，MAGMOM 应结合实验磁矩、文献值或化学直觉给出，而不是随意填 1.0。
三、k点网格、截断能与展宽：精度与稳定性的平衡
k点网格决定布里渊区积分精度。对金属磁性体系，k点过稀会导致费米面附近积分误差大，磁矩振荡甚至不收敛。一般过渡金属体相材料，Monkhorst-Pack 网格可取 9×9×9 至 15×15×15；二维材料在真空方向只需 1 个 k 点，如 12×12×1。判断 k 点是否足够，可比较相邻两套网格的总能与磁矩变化，若磁矩变化小于 0.01 μB，通常可认为收敛。
截断能 ENCUT 决定平面波基组大小。对含磁性 3d 元素的体系，ENCUT 一般取 400～520 eV。若使用 PAW 赝势，可参考 POTCAR 中 ENMAX 的 1.3～1.5 倍。截断能不足会使磁矩偏小、能量偏高，且这种误差不会因自洽收敛而消失。
展宽方法 ISMEAR 与 SIGMA 对金属磁性体系尤为关键。金属推荐 ISMEAR=1（Methfessel-Paxton）或 ISMEAR=2（Fermi 展宽），SIGMA 取 0.05～0.2 eV。若 SIGMA 过大，磁矩会被“抹平”；若过小，自洽过程容易振荡。对绝缘体或半导体，则用 ISMEAR=0 并配合小 SIGMA。
四、收敛判据：能量、力与磁矩三者都要看
自洽收敛由 EDIFF 控制，通常取 1E-5～1E-6 eV。但磁性计算不能只看能量收敛，还要检查磁矩是否稳定。建议同时满足：
- 能量变化小于
EDIFF； - 总磁矩在最后若干离子步中变化小于 0.01 μB；
- 各原子局域磁矩不再明显漂移；
- 若做结构弛豫，还需
EDIFFG控制力收敛，一般取 -0.01～-0.02 eV/Å。
只盯能量收敛，是磁性计算中最常见的误判来源。
五、磁性收敛陷阱与排查路径
磁性计算常见的收敛陷阱包括：收敛到非磁态、磁矩振荡、亚稳态磁序、以及“能量收敛但磁矩未收敛”。下图给出从参数设置到结果判读的排查路径，帮助定位问题出在开关、初始磁矩、k点、展宽还是收敛判据。
排查时建议按“开关→初始磁矩→k点→展宽→收敛判据”的顺序逐项检查。若磁矩始终为零，先确认 ISPIN=2；若磁矩振荡，优先调大 SIGMA 或加密 k 点；若收敛到异常磁矩，尝试改变 MAGMOM 初值并比较总能，判断是否落入亚稳态。
SIGMA（如 0.2 eV）和较稀 k 点快速试探，确认磁矩量级后再逐步收紧参数做高精度计算。这样既能避免一开始就陷入收敛困难，也能节省机时。
EDIFF 而忽略磁矩，可能得到一个能量看似收敛、但磁态错误的解。务必在 OUTCAR 中同时检查 energy(sigma->0) 与 tot 磁矩的迭代历史。
综上，自旋计算参数设置的核心逻辑是：用 ISPIN 打开自旋通道，用 MAGMOM 给出合理起点，用 k 点、截断能和展宽控制精度与稳定性，最后用能量与磁矩双判据确认收敛。只有把这六个参数作为一个整体来调，才能稳定得到物理上正确的磁态。