VASP SOC与磁各向异性入门
用铁的三个磁化方向说明SOC设置、微小能差和易轴判断的精度要求。
阅读本篇目录
自旋轨道耦合把自旋方向与晶体方向联系起来。磁各向异性通常比较同一结构在不同磁化方向下的能量，往往比一般结构能差小得多。
1. 工作顺序#
- 先完成无SOC的共线自旋SCF，确认总磁矩和局域磁态，并保存CHGCAR。
- 三个SOC目录使用完全相同的晶胞、赝势、k网格、NBANDS与展宽，从同一收敛密度启动。
- 使用vasp_ncl，开启LSORBIT，分别设定SAXIS；同时检查MAGMOM在所用自旋坐标系下的含义。
- 分别提交并检查电子收敛以及最终磁矩方向；若研究中必须固定方向，应按约束磁矩方法显式约束，不能仅凭初始方向认定保持不变。
- 比较同一能量口径的差值，报告每晶胞或每磁性原子的单位。
本例bcc Fe晶胞2个原子、11³ k网格、ISYM=−1、ISMEAR=1、SIGMA=0.1 eV，EDIFF=1E-6。这里演示设置与能差读取。
2. 实际能差#
| 方向 | 总能量 / eV·cell⁻¹ | 相对[001] / meV·Fe⁻¹ |
|---|---|---|
| 001 | -16.47273265 | 0.000000 |
| 100 | -16.47273674 | -0.002045 |
| 111 | -16.47277217 | -0.019760 |
立方对称下[001]与[100]应等价，本例仍相差约0.002 meV/Fe，这是很有用的数值精度参照。[111]的差值也很小，所以本页不据此宣布可靠易轴。
3. 定量研究需要补什么#
收紧电子能量阈值，增加k点，比较展宽和能量口径，再检查磁矩方向。用户日常结构参考采用1E-6没有问题，但它不是所有微电子伏量级磁性结论的保证。
多个FM/AFM能量不能直接给出居里温度。交换参数、DMI、自旋螺旋、磁振子和有限温度磁性还需要明确磁哈密顿量、相互作用范围及后续统计模型。普通能量排序与这些性质应分开处理。
环境与下载#
本页数据在鸢算科技集群实际计算。VASP 6.6.1、PBE PAW版本64、ENCUT=550 eV、EDIFF=1E-6；所用优化结构按自由原子力0.02 eV/Å准备。GPU程序采用NVHPC 25.3、CUDA 12.8、HPC-X 2.22.1和MKL 2023.2。方法专用参数以各目录INCAR为准。
解压后先看README。准备有权使用、与POSCAR元素顺序匹配的POTCAR，按目录run.sh提交Slurm作业；后处理结果须由自己的新计算重新提取。附件中的reference结果是本次参考值。公开包不含POTCAR、WAVECAR及原始私有日志。
需要根据你的体系选择资源？联系讨论计算需求。实际可用软件与权限，以当前账户和运行环境为准。