铁磁性计算 (Ferromagnetism - BCC Fe)
这是最基础的入门部分，我们以体心立方铁（BCC Fe）为例，学习如何通过 nsppol 和 spinat 变量开启自旋极化，并让系统产生净磁矩。
输入文件设置 (tspin_1.abi)
这个输入文件设计得很巧妙，在一次运行中执行两个 Datasets：
- Dataset 1: 不考虑自旋极化（非磁性）。
- Dataset 2: 考虑自旋极化（铁磁性）。
结果分析
运行 abinit tspin_1.abi > tspin_1.log 后，我们可以从输出文件中观察到几个关键点：
总磁矩 (Total Magnetization)
在 Dataset 2 的结果中，搜索 "Total magnetization"，可以看到如下数据：
Total magnetization (exact up - dn): 1.822358这表示整个晶胞的完整磁矩约为 1.82 $\mu_B$。磁矩的来源本质上是自旋向上和自旋向下的电子数量不平衡： $$ \text{Total Magnetization} = \text{Spin Up Density} - \text{Spin Down Density} \approx 8.9112 - 7.0888 = 1.8224 $$
能量最小化原理
对比非磁性和铁磁性的总能量 (etotal)：
- Dataset 1 (非磁性):
-123.5378...Ha - Dataset 2 (铁磁性):
-123.5463...Ha
显然，铁磁态的能量更低。根据能量最小化原理，这说明铁在自然界中的基态确实是铁磁性的（Ferromagnetic）。
占据数与 DOS
观察本征值和占据数 (occupation numbers)，我们会发现：
- Spin Up (多数自旋): 能带填充较满，能级向低能方向移动。
- Spin Down (少数自旋): 对应能带几乎是空的。
这种能级移动是由交换作用 (Exchange Interaction) 引起的，也是铁磁态稳定的微观原因。
反铁磁性计算 (Antiferromagnetism - FCC Fe)
接下来我们挑战面心立方铁（FCC Fe）的反铁磁性计算（相邻原子磁矩相反，净磁矩为零）。
利用 Shubnikov 对称性加速计算 (tspin_2.abi)
在这个计算中，我们使用了一个特殊的技巧：
nsppol = 1: 只计算一套波函数。nspden = 2: 分别存储自旋向上和向下的密度。
原理：在完美的反铁磁晶体中，原子 A 的“向上自旋”环境和原子 B 的“向下自旋”环境是完全对称的（Shubnikov 对称性）。Abinit 利用这一点，只计算一半波函数，通过对称变换得到另一半，从而将计算速度提高约 2倍。
结果解读
在 tspin_2.log 中，你会发现所有的占据数 (occupation numbers) 都是 2.00000。这看起来像非磁性，但实际上是因为程序利用对称性简化了显示。
真正的磁性证据在于 Magnetization 数据：
- Maximum (Atom A):
3.1883E-01(正磁性) - Minimum (Atom B):
-3.1883E-01(负磁性) - Integrated (Total):
0.0000E+00(总磁矩为零)
这完美地展示了反铁磁性的特征：局部磁矩大小相等、方向相反，总磁矩相互抵消。
不利用对称性的反铁磁计算 (General Method)
为了验证上述结果，或者是处理结构不那么完美的材料，我们可以使用通用的方法：显式指定每个原子的自旋方向。
输入设置 (tspin_3.abi)
nsppol = 2(全自旋计算)spinat: 显式设置原子 1 为正，原子 2 为负（例如0 0 4和0 0 -4）。
结果对比
- 总能量 (
etotal):tspin_3的计算结果与利用技巧的tspin_2结果惊人地一致（小数点后多位相同）。这说明“条条大路通罗马”，只要物理状态相同，基态能量就是确定的。 - 局部磁矩:
- Atom 1:
1.137602 - Atom 2:
-1.135050 - Total:
0.002512($\approx 0$) 再次证实了反铁磁基态。
- Atom 1:
- DOS 对比: 对比两个原子的态密度，会发现 Atom 1 的 Spin Up 曲线和 Atom 2 的 Spin Down 曲线几乎完全重合（互为镜像）。
自旋轨道耦合 (Spin-Orbit Coupling, SOC)
最后是进阶内容：自旋轨道耦合。这在重原子或强磁性材料中非常关键。
能级分裂 (tspin_6.abi)
运行包含 SOC 的计算后，我们可以观察到原本简并的能级发生了分裂。以 5d 轨道 为例：
- 低能级组:
-0.91225Ha (4个状态) - 高能级组:
-0.80130Ha (6个状态)
能量换算
计算分裂能 $\Delta E$： $$ \Delta E_{SOC} = 0.11095 \text{ Ha} \approx 3.02 \text{ eV} $$
同样，对于 6p 轨道（分裂为 4 个和 2 个状态）： $$ \Delta E_{SOC} \approx 0.070 \text{ Ha} \approx 1.90 \text{ eV} $$
这些巨大的能量分裂表明，在重元素中，自旋轨道耦合效应大到足以显著影响材料的化学成键和物理性质。
总结
通过这个教程，我们掌握了 Abinit 中处理自旋的三种核心模式：
- 铁磁性:
nsppol=2+spinat。 - 反铁磁性 (技巧):
nsppol=1,nspden=2利用对称性。 - SOC: 处理能级分裂和非共线磁性。
本练习的计算文件可以从链接中下载：点击这里