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$は、キタエフ量子スピン液体状態を実現するための主要な候補物質であるが、その基本的な磁気特性については依然として議論の対象となっている。主な未解決の問題は以下の通りである：

1. バルク構造： 近年の実験では低温度相が$R\bar{3}$空間群に属すること（約150 Kで$C2/m$へ転移）が示唆されているが、$R\bar{3}$構造がエネルギー的に有利な低温度相であることを確認する理論的研究はこれまで存在しなかった。
2. 電子構造： $\alpha$-RuCl$_3$における$J_{\text{eff}} = 1/2$の描像の妥当性、特にネールベクトルに対する角運動量量子化軸の向きに関する検証が必要である。
3. 磁気ハミルトニアン： 従来の研究では、磁気交換パラメータ（例：最近接キタエフ交換$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}}$基底状態に投影することで解析した。極めて重要な点として、$J_{\text{eff}}$の描像をテストするために、角運動量量子化軸をネールベクトルに平行に合わせた。
• 磁気ハミルトニアンの構築： 様々な磁気構成（ジグザグ、ネール、ストライピー、強磁性、および特定の1-by-4構成）からの全エネルギーを用いて、以下の3つのモデルハミルトニアンにフィッティングを行った：
  1. 従来のキタエフ・ハイゼンベルク（KH）モデル。
  2. スピンアイス（SI）モデル。
  3. 異方的な第二および第三近接相互作用を含む、新しい**異方的$J_1$-$J_2$-$J_3$（A-$J_1J_2J_3$）**モデル。
• $g$因子計算： 2つの手法を比較した：
  1. 原子軌道投影法： ハバード多様体にモーメントを投影する。
  2. Wannier補間を用いたモーメント傾斜法： $J_{\text{eff}}$モーメントを1°傾け、並進等価なWannier補間（WannierBerri）を用いて全磁気モーメントの変化を計算する。
• 歪み解析： 完全な八面体からの2つの特定の歪みを分離した：三方歪み（Ru-Cl-Z角の変化）と、相対的なねじれ（XY平面内での上下のCl三角形の回転）。

主要な貢献と結果

1. $R\bar{3}$積層の確認：
   本研究は、$R\bar{3}$構造が異なるHubbard $U$値および緩和タイプにおいて、$C2/m$構造よりも約5–6 meV/Ru低エネルギーであることを示している。これは、実験的に観察されている低温度$R\bar{3}$相に対する初の理論的裏付けとなる。

2. 洗練された$J_{\text{eff}}$の描像：
   伝導帯下端の状態は、角運動量量子化軸をネールベクトルに平行に設定した場合にのみ、ほぼ純粋な$J_{\text{eff}} = 1/2$および$m_{\text{eff}} = -1/2$の性質を持つことが示された。この視点は、これまで無視されていたが、状態の軌道特性を正しく特定するために不可欠である。$J_{\text{eff}}$モーメントの大きさは、混成や歪みにより理想的な$1/2$よりはわずかに小さいものの、各構成間でロバストであることが判明した。

3. 磁気ハミルトニアンと交換パラメータ：

• 従来のKHモデルおよびスピンアイスモデルは、特にネール構成において、DFT全エネルギーを正確に再現することに失敗している。
• 提案された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_{XY} \approx 2.22$および$g_Z \approx 1.84$の$g$因子が得られた。これらの値は小さな異方性を示しており、最近の実験報告（Loidlら）と一致し、大きな異方性を主張する以前の報告とは矛盾する。
• 原子軌道投影法は、$g$因子の大きさを著しく過小評価しており（$g_{XY} \approx 1.54$）、全磁気モーメントの無視できない部分を見落としている。著者らは、正確な軌道磁化の計算にはWannier補間アプローチが必要であると結論付けている。

5. 構造歪みの役割：
   歪みを分離して解析することにより、上下のCl三角形の相対的なねじれ（XY平面内での回転）が、磁気異方性と正しい基底状態のモーメント方向（$\theta \sim 60^\circ$）を決定する支配的な要因であることを突き止めた。三方歪み単独では、実験的な基底状態を再現できない。このねじれ歪みは、外部からの摂動なしにバルク材料内に存在しており、従来の理論モデルでは見落とされてきたものである。

意義
本論文は、$\alpha$-RuCl$_3$の構造および磁気的記述における長年の矛盾を解決することにより、包括的な第一原理的理解を提供することを目的としている。

• $R\bar{3}$低温度構造を理論的に検証した。
• $J_{\text{eff}}$の描像を正しく解釈するためには、量子化軸をネールベクトルに合わせることが不可欠であることを確立した。
• 正確な磁気モデリングには、異方的な第二近接交換相互作用を含める必要があることを示した。
• $g$因子の計算手法を修正し、原子軌道投影法は不十分であり、真の異方性は小さいことを示した。
• Cl三角形の相対的なねじれが、磁気異方性を支配する、これまで見落とされてきた重要な構造パラメータであることを特定し、キタエフ材料における磁性を制御するための新たな道を提示した。

著者らは、本研究が$\alpha$-RuCl$_3$の磁性に対する異なる視点を提供し、今後の理論モデルや高圧実験においては、これらの特定の構造歪みや交換項を取り入れるべきであることを示唆している。