VASP 磁性入门：磁矩与不同磁态比较
用同一铁晶胞比较自旋极化、反平行初始磁矩和非磁态，理解初态、最终磁矩及能量比较。
阅读本篇目录
研究磁性时，第一步通常是找到合适的磁态，而不是直接去算居里温度。设置初始磁矩、检查最终磁矩、比较匹配条件下的能量，是后续SOC、交换参数和磁激发的基础。
1. 从什么结构开始#
本例使用晶格常数2.87 Å、含两个Fe原子的bcc常规晶胞。同一晶胞可以分别给两个原子同向和反向初始磁矩。三个计算保持结构、PBE泛函、550 eV截断、13×13×13 Gamma网格与EDIFF=1E-6一致，采用静态计算。
这是一组固定几何下的磁态比较。如果要比较各磁态的平衡晶格，应分别优化后再作相容比较；也应确认超胞能够容纳要研究的磁周期。
2. 先准备三套输入#
| 目录 | 关键设置 | 代表什么初始条件 |
|---|---|---|
| fe-fm | ISPIN=2，MAGMOM=2.5 2.5 | 两个原子同向 |
| fe-afm | ISPIN=2，MAGMOM=2.5 -2.5 | 两个原子反向 |
| fe-nm | ISPIN=1 | 非自旋极化计算 |
ISPIN = 2
MAGMOM = 2.5 -2.5
ISYM = 0
LORBIT = 11
准备各目录POTCAR后，分别执行 sbatch run.sh。这组首次计算不继承其他磁态的WAVECAR/CHGCAR，以免初态设置与实际读取的数据不一致。实际续算时要理解重启文件中已经包含的自旋信息。
3. 算完看哪几个位置#
- 检查电子迭代达到EDIFF，没有撞到NELM上限。
- 看OSZICAR最后一行的总磁矩；总磁矩接近零不一定是非磁，反铁磁也可能相互抵消。
- 看OUTCAR的局域磁矩表，确认两个原子的最终符号和大小。
- 用同一定义的最终能量，统一到每原子后比较。
- 将结果按最终状态命名，不能仅凭目录名说“得到了反铁磁”。
| 目录 | E / eV·cell⁻¹ | 相对能量 / meV·Fe⁻¹ | 总磁矩 / μB·cell⁻¹ | 局域Fe磁矩 / μB |
|---|---|---|---|---|
| fe-fm | -16.4598066 | 0.000 | 4.4686 | +2.241，+2.241 |
| fe-afm | -15.5400940 | 459.856 | 0.0004 | +1.703，−1.703 |
| fe-nm | -15.2921546 | 583.826 | 0.0000 | 非自旋极化 |
MAGMOM给的是初始猜测，并不强制最终保留该磁态。不同初值可能落入不同局域解；复杂氧化物还可能需要多种磁序、局域畸变或合理的DFT+U设置。
4. 总磁矩、原子磁矩和自旋密度有何区别#
总磁矩是整个晶胞的自旋不平衡。OUTCAR中的原子磁矩来自相应投影区域，不一定严格相加等于总磁矩；差额可能涉及间隙区域与投影定义。自旋密度则描述空间中ρ↑−ρ↓的分布，适合观察磁性主要位于哪些原子或键区。
不要把Bader电荷与局域磁矩混为一谈；一个分析电子数分配，另一个分析自旋不平衡。
5. 下一步能不能算交换参数和温度#
多个磁构型的能量可以映射到明确约定的自旋模型，例如Heisenberg模型。必须写清自旋长度、J的正负号、键数和求和约定；一组FM/AFM能差只约束少量模型参数，不能自动得到完整交换网络。
磁振子、居里/奈尔温度通常还需要自旋波或统计模拟，并检查模型范围。单次自旋极化计算的总磁矩不等于居里温度；同样，看到某个磁态能量更低，也不足以说明实验温度下必然保持有序。
本例VASP 6.6.1 GPU构建采用NVHPC 25.3、CUDA 12.8与配套MPI。涉及非共线/SOC时使用对应 vasp_ncl，不要继续调用 vasp_std。
下载输入、实算数据与图片（ZIP） · 下载44项结果表（CSV）。公开包不含POTCAR和WAVECAR。