First-principles study of bulk stacking, $J_{\rm eff}$ picture, magnetic Hamiltonian, $g$ factors, and structural distortions of $α$-RuCl$_3$

技术摘要：关于 $\alpha$-RuCl$_3$ 的体相堆叠、$J_{\text{eff}}$ 图像、磁哈密顿量、$g$ 因子及结构畸变的基态原理研究

问题陈述
$\alpha$-RuCl$_3$ 是实现 Kitaev 量子自旋液体状态的主要候选材料，但其基本磁学性质仍存在争议。关键的未解决问题包括：

1. 体相结构： 虽然近期的实验表明其低温相属于 $R\bar{3}$ 空间群（在 $\sim$150 K 时向 $C2/m$ 转变），但此前尚无理论研究证实 $R\bar{3}$ 结构是能量上有利的低温相。
2. 电子结构： $\alpha$-RuCl$_3$ 中 $J_{\text{eff}} = 1/2$ 图像的有效性需要验证，特别是角动量量子化轴相对于 Néel 矢量的取向问题。
3. 磁哈密顿量： 先前的研究报告了磁交换参数（例如最近邻 Kitaev 交换项 $K_1$ 的符号）存在显著差异，并且往往忽略了第二近邻相互作用。
4. $g$ 因子各向异性： 关于 $g$ 因子各向异性（$g_{XY}$ vs. $g_Z$）的实验报告存在矛盾，范围从微小到巨大的各向异性不等。此外，计算这些因子的理论方法（特别是原子轨道投影法）可能并不充分。
5. 结构畸变： 对 RuCl$_6$ 八面体中特定结构畸变（三方畸变 vs. Cl 三角形的相对扭转）如何影响磁各向异性的研究尚未得到系统性的分离分析。

方法论
作者在 Quantum ESPRESSO 框架内，采用局部密度近似（LDA）和旋转不变 LDA+U 方法，使用了约束密度泛函理论 (cDFT)。

• 结构弛豫： 针对 $C2/m$ 和 $R\bar{3}$ 堆叠构型，在实验、固定胞和可变胞弛豫条件下进行计算，以确定能量稳定性。
• $J_{\text{eff}}$ 分析： 通过将波函数投影到 $J_{\text{eff}}$ 基态上进行电子结构分析。至关重要的是，将角动量量子化轴与 Néel 矢量平行对齐，以测试 $J_{\text{eff}}$ 图像。
• 磁哈密顿量构建： 作者将来自各种磁构型（zigzag、Néel、stripy、铁磁以及特定的 1-by-4 构型）的总能量拟合到三种模型哈密顿量中：
  1. 常规 Kitaev-Heisenberg (KH) 模型。
  2. 自旋冰 (SI) 模型。
  3. 一种包含各向异性第二近邻和第三近邻相互作用的新型各向异性 $J_1$-$J_2$-$J_3$ (A-$J_1J_2J_3$) 模型。
• $g$ 因子计算： 对比了两种方法：
  1. 原子轨道投影法： 将磁矩投影到 Hubbard 流形上。
  2. 利用 Wannier 插值进行的磁矩偏转法： 将 $J_{\text{eff}}$ 磁矩倾斜 $1^\circ$，并使用平移等变 Wannier 插值（WannierBerri）计算总磁矩的变化。
• 畸变分析： 作者隔离了来自完美八面体的两种特定畸变：三方畸变（改变 Ru-Cl-Z 角度）和相对扭转（在 XY 平面内旋转上下两个 Cl 三角形）。

核心贡献与结果

1. 确认 $R\bar{3}$ 堆叠：
   本研究表明，$R\bar{3}$ 结构的能量低于 $C2/m$ 结构（在不同 Hubbard $U$ 值和弛豫类型下，每 Ru 低约 5–6 meV）。这为实验观察到的低温 $R\bar{3}$ 相提供了首个理论支持。

2. 精炼的 $J_{\text{eff}}$ 图像：
   研究显示，导带底态几乎纯粹由 $J_{\text{eff}} = 1/2$ 和 $m_{\text{eff}} = -1/2$ 特征组成，但前提是必须将角动量量子化轴设定为与 Néel 矢量平行。这一此前被忽视的视角对于正确识别态的轨道特征至关重要。尽管由于杂化和畸变导致 $J_{\text{eff}}$ 磁矩略小于理想值，但其量级在不同构型下表现得十分稳健。

3. 磁哈密顿量与交换参数：

• 常规 KH 模型和自旋冰模型无法准确重现 DFT 总能量，特别是在 Néel 构型下。
• 提出的 A-$J_1J_2J_3$ 模型（以及简化的 A-$J_1J_2$-I-$J_3$ 变体）能够准确重现 DFT 能量。
• 第二近邻相互作用至关重要。第二近邻交换的各向异性分量（$K_2, \Gamma_2, \Gamma'_2$）不可忽略，且在量级上大于通常被考虑的最近邻各向异性项（$\Gamma'_1$）。包含这些项使模型相对于常规模型的总能量均方根（RMS）误差降低了 3.6 倍。

4. $g$ 因子各向异性与计算方法：

• 结合 Wannier 插值的磁矩偏转法得出 $g$ 因子为 $g_{XY} \approx 2.22$ 且 $g_Z \approx 1.84$。这些数值表明存在较小的各向异性，符合近期实验报告（如 Loidl 等人），并反驳了早期关于存在大各向异性的说法。
• 原子轨道投影法显著低估了 $g$ 因子的大小（$g_{XY} \approx 1.54$），遗漏了相当一部分总磁矩。作者得出结论，准确计算轨道磁化率需要使用 Wannier 插值法。

5. 结构畸变的作用：
   通过隔离畸变，研究发现上下 Cl 三角形的相对扭转（在 XY 平面内）是决定磁各向异性和正确基态磁矩方向（$\theta \sim 60^\circ$）的主导因素。仅靠三方畸变无法重现实验基态。这种在本体材料中无需外部扰动即存在的扭转畸变，在以往的理论模型中被忽视了。

意义
该论文声称通过解决 $\alpha$-RuCl$_3$ 在结构和磁性描述上的长期矛盾，提供了全面的第一性原理理解。具体而言：

• 它从理论上验证了 $R\bar{3}$ 低温结构。
• 它确立了将量子化轴与 Néel 矢量对齐对于正确解释 $J_{\text{eff}}$ 图像的必要性。
• 它证明了准确的磁建模需要包含各向异性的第二近邻交换相互作用。
• 它纠正了 $g$ 因子的计算方法，表明原子轨道投影是不充分的，且真实的各向异性较小。
• 它识别出 Cl 三角形的相对扭转是决定磁各向异性的关键、此前被忽视的结构参数，这为工程化设计 Kitaev 材料的磁性提供了一条新途径。

作者总结道，他们的工作为 $\alpha$-RuCl$_3$ 的磁性提供了不同的视角，并建议未来的理论模型和高压实验应纳入这些特定的结构畸变和交换项，以实现准确的解释。