Abinit 自旋计算教程指导

铁磁性计算 (Ferromagnetism - BCC Fe)

这是最基础的入门部分，我们以体心立方铁（BCC Fe）为例，学习如何通过 nsppol 和 spinat 变量开启自旋极化，并让系统产生净磁矩。

输入文件设置 (tspin_1.abi)

这个输入文件设计得很巧妙，在一次运行中执行两个 Datasets：

1. Dataset 1: 不考虑自旋极化（非磁性）。
2. 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）。

结果对比

1. 总能量 (etotal): tspin_3 的计算结果与利用技巧的 tspin_2 结果惊人地一致（小数点后多位相同）。这说明“条条大路通罗马”，只要物理状态相同，基态能量就是确定的。
2. 局部磁矩:
   • Atom 1: 1.137602
   • Atom 2: -1.135050
   • Total: 0.002512 ($\approx 0$) 再次证实了反铁磁基态。
3. DOS 对比: 对比两个原子的态密度，会发现 Atom 1 的 Spin Up 曲线和 Atom 2 的 Spin Down 曲线几乎完全重合（互为镜像）。

自旋轨道耦合 (Spin-Orbit Coupling, SOC)

最后是进阶内容：自旋轨道耦合。这在重原子或强磁性材料中非常关键。

能级分裂 (tspin_6.abi)

运行包含 SOC 的计算后，我们可以观察到原本简并的能级发生了分裂。以 5d 轨道 为例：

• 低能级组: -0.91225 Ha (4个状态)
• 高能级组: -0.80130 Ha (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 中处理自旋的三种核心模式：

1. 铁磁性: nsppol=2 + spinat。
2. 反铁磁性 (技巧): nsppol=1, nspden=2 利用对称性。
3. SOC: 处理能级分裂和非共线磁性。

本练习的计算文件可以从链接中下载：点击这里