← Back
npj Quantum Materials

Kitaevモデル候補RuBr$_3$の磁気励起

New study reveals RuBr3's magnetic interactions push it from ideal spin liquid state, offering clues for quantum computing materials.

研究分野 Quantum Information
Article Type Research analysis
Authors Nawa et al.
Original Paper Published 2026-03-11
ISOM Posted 2026-04-07 16:00 UTC
Read Time 40M
Open PDF Open DOI Open Source Page

Editorial Disclosure

ISOM follows an editorial workflow that structures the source paper into a readable analysis, then publishes the summary, source links, and metadata shown on this page so readers can verify the original work.

The goal of this page is to help readers understand the paper's core question, method, evidence, and implications before opening the original publication.

背景と学術的系譜

起源と学術的系譜

本論文で取り上げられている問題は、凝縮系物理学の魅力的な分野、特に量子スピン液体（QSL）状態の探索と特徴付けに端を発する。この探求は、2006年にアレクセイ・キタイエフが理論的に提案したキタイエフモデル[1]から本格的に始まった。ハニカム格子上の特異な結合方向依存性を持つイジング相互作用によって定義されるこのモデルは、厳密に解くことができ、QSLの基底状態を明らかにする点で注目に値する。このエキゾチックな状態では、スピンは絶対零度でも秩序化せず、代わりに絡み合ったまま揺らぎ、マヨラナフェルミオンとして知られる分数化励起[7]を生じさせる。これらの特性は、耐故障性量子計算に計り知れない可能性を秘めている。

次に課せられた課題は、これらの捉えどころのないキタイエフ相互作用をホストしうる実在材料を見つけることであった。学術分野は、2009年頃、特定の遷移金属化合物、特にIr$^{4+}$やRu$^{3+}$イオンのような$d^5$電子配置を持つものにおける強いスピン軌道結合が、擬スピン1/2状態を介したメカニズムを通じて、これらの結合方向依存性相互作用を効果的に誘起できるという認識とともに、大きなブレークスルーを迎えた[8]。これにより、「キタイエフ材料」が新たなフロンティアとして出現した。これらのうち、$\alpha$-RuCl$_3$は、最も広範に研究された候補として急速に台頭した[9-38]。

しかしながら、$\alpha$-RuCl$_3$に関する従来の取り組みや研究の根本的な限界、すなわち「ペインポイント」は、高温ではキタイエフスピン液体の振る舞いのいくつかの兆候を示すものの、ある遷移温度($T_N$)以下の基底状態は従来のジグザグ反強磁性秩序であることである。これは、理想的なキタイエフスピン液体状態と競合し、最終的に不安定化させる他の磁気相互作用（ハイゼンベルクや対角線外項など）の存在を示唆している[31, 32, 39-43]。これらの競合する相互作用は、望ましいQSL特性の完全な実現を妨げる。本論文の著者らは、これらの相互作用を理解し制御する必要性に動機づけられている。一つの重要な戦略は、磁性原子を取り囲む配位子イオンを変更することである。これらのイオンは、結合方向依存性相互作用を生じさせる軌道交換メカニズムを媒介する上で重要な役割を果たす。本論文は、塩素を臭素に置換することでこれらの磁気相互作用のバランスが変化し、材料が理想的なキタイエフスピン液体状態に近づくか、あるいは遠ざかることが期待される新しい候補としてRuBr$_3$を調査する。

直感的なドメイン用語

  1. キタイエフモデル (Kitaev Model): 隣接するチェッカー間の相互作用（引き合うか反発するか）のルールが、それらを結ぶ線の方向に完全に依存する、特別な種類の磁気チェッカーボードを想像してほしい。例えば、水平に接続されたチェッカーは同じ方向に整列したがるかもしれないが、対角線に接続されたチェッカーは反対方向に整列したがるかもしれない。キタイエフモデルはこのようなシステムを記述する理論的枠組みであり、「量子スピン液体」状態を予測するため有名である。この状態では、磁石は最も低い温度でも固定パターンに凍結することはない。
  2. 量子スピン液体 (Quantum Spin Liquid, QSL): 水分子が絶えず動き回り相互作用している通常の液体を考えてほしい。量子スピン液体も同様だが、水ではなく、電子の微小な内部磁石（スピン）が、絶えず絡み合った「液体のような」状態にある。絶対零度でも、冷蔵庫の磁石のような単純で秩序だった磁気パターンに落ち着くことはない。代わりに、それらの集団的な振る舞いは非常に複雑で量子力学的であり、マヨラナフェルミオンと呼ばれるエキゾチックな「半粒子」を生じさせる。
  3. マヨラナフェルミオン (Majorana Fermions): 電子のような通常の粒子を「一人の人間」と考えるなら、マヨラナフェルミオンは「半人」のようなものである。量子スピン液体の文脈では、電子のスピンは効果的に2つの独立した部分に分裂し、それぞれがマヨラナフェルミオンのように振る舞うことができる。これらの「半粒子」は、それ自身の反粒子であるという点でユニークであり、局所的な摂動に対する堅牢性から、将来の量子コンピュータの有望な構成要素となっている。
  4. 結合方向依存性イジング相互作用 (Bond-dependent Ising Interactions): 微小なコンパス針のグリッドを想像してほしい。各針は厳密に上か下（バイナリスイッチのように）しか向けない。「結合方向依存性」とは、隣接する2つの針が同じ方向か反対方向を向きたがるかどうかは、グリッド上のそれらを結ぶ線の特定の方向に完全に依存することを意味する。例えば、「南北」結合に沿っては平行を好むかもしれないが、「東西」結合に沿っては反平行を好むかもしれない。この方向性への好みは、キタイエフモデルの核心的な考え方である。
  5. 擬スピン1/2状態 (Pseudospin-1/2 State): いくつかの複雑な材料では、電子の実際のスピン（常に1/2である）が原子周りの軌道運動と絡み合っている。この組み合わせにより、単純なスピン1/2粒子であるかのように振る舞う実効的な磁気モーメントが生じるが、その「上向き」と「下向き」の方向は、単純な磁場だけでなく、結晶構造によって決定される。これは、これらの材料における複雑な電子の振る舞いを記述する巧妙な方法である。

記法表

記法 説明
$H$ 磁気系の全ハミルトニアン
$S_i$ サイト$i$におけるスピン演算子
$J_{ij}$ 最近接結合異方性相互作用行列
$J_1$ 等方性ハイゼンベルク交換結合定数
$K$ キタイエフ相互作用項
$\Gamma$ 対称対角線外交換項
$\Gamma'$ 反対称対角線外交換項
$J_2$ 次最近接結合（NNN）等方性磁気相互作用の結合定数
$J_3$ 第三最近接結合（3NN）等方性磁気相互作用の結合定数
$J_p$ 層間等方性磁気相互作用の結合定数
$T_N$ ネール温度（磁気秩序化温度）
$E_i$ 入射中性子エネルギー
$S(\mathbf{Q}, E)$ 動的構造因子（非弾性中性子散乱で測定）
$\mathbf{Q}$ 波ベクトル（運動量移動）
$E$ エネルギー移動
$n(T)$ ボーズ因子
$k_B$ ボルツマン定数
$\theta_{ab}, \theta_c$ ワイス温度（面内および面外）

問題定義と制約

コア問題の定式化とジレンマ

本論文で取り上げられている中心的な問題は、キタイエフ量子スピン液体の候補材料であるRuBr$_3$の磁気相互作用を精密に特徴付け、これらの相互作用がその基底状態を理想的なキタイエフスピン液体からどのように逸脱させるかを理解することである。

入力/現在の状態は、よく研究されているキタイエフ候補$\alpha$-RuCl$_3$と構造的に同一である層状ハニカム構造を持つ新しい多形であるRuBr$_3$という材料である。両材料は、それぞれのネール温度($T_N$)以下のジグザグ反強磁性（AFM）秩序を示す。しかし、RuBr$_3$は、磁化率の抑制や異なるワイス温度など、$\alpha$-RuCl$_3$とは異なる磁気特性を示す。RuBr$_3$に関する以前のab initio計算では、これらの実験的観測と矛盾する、$\alpha$-RuCl$_3$に匹敵する磁気相互作用が示唆されている。理想的な強磁性キタイエフモデルからの理論的予測は、波ベクトル依存性の弱いスピン励起を予測するが、実在材料は競合する相互作用のために、しばしばより複雑な振る舞いを示す。

望ましい終点（出力/目標状態）は、キタイエフ候補材料におけるアニオン置換（特に臭素を塩素に置換）が磁気相互作用に与える影響の明確化された理解である。本論文は、非弾性中性子散乱（INS）を用いてRuBr$_3$のスピンダイナミクスを実験的に調査し、観測されたジグザグ反強磁性秩序を安定化させる特定の磁気相互作用（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）を特定することを目的とする。最終的な目標は、RuBr$_3$の磁気励起と巨視的特性を正確に再現する、より現実的なモデルを提供することであり、それによってその基底状態が理想的な強磁性キタイエフスピン液体からどのように逸脱するかを説明することである。

欠落しているリンクまたは数学的ギャップは、RuBr$_3$の巨視的磁気特性（ワイス温度やキャンティング角など）と微視的スピンダイナミクス（磁気励起スペクトル）の両方を一貫して説明できる交換パラメータの正確なセットを定量的に決定することにある。理論モデル（例：J$_1$-K-$\Gamma$-$\Gamma$'およびJ$_1$-K-J$_2$-J$_3$モデル）は存在するが、本論文は、粉末平均実験データから複数の交換パラメータを正確に推定することの困難さを強調している。ab initio計算からの理論的予測と、RuBr$_3$の観測された実験的磁気挙動との間には明確な不一致があり、本研究はこの不一致をスピンダイナミクスの詳細な実験的特徴付けを通じて橋渡ししようとしている。

過去の研究者を閉じ込めてきた痛みを伴うトレードオフまたはジレンマであり、本研究の中心でもあるのは、実在材料における望ましいキタイエフ相互作用と他の磁気相互作用（ハイゼンベルク、対角線外、さらなる近接結合）との固有の競合である。キタイエフモデルはスピン液体を実現するために理論的に魅力的であるが、これらの追加相互作用の存在は、しばしばRuBr$_3$で観測されるジグザグ反強磁性のような従来の磁気秩序につながる。研究者は、これらの競合する相互作用を、従来の秩序を抑制するか、キタイエフ物理学を強化するかのジレンマに直面している。配位子置換は一つの戦略であるが、その相互作用の繊細なバランスへの正確な影響は複雑で容易に予測できず、しばしば意図しない結果や新たな課題につながる。本論文は、「キタイエフ相互作用だけでなく、ハイゼンベルクおよび対角線外磁気相互作用を制御することも、実在材料の実用的な利用に不可欠である」と明記しており、この根本的なトレードオフを強調している。さらに、粉末平均実験データから複数の交換パラメータを一意に決定することの困難さは、実験的アクセス性（粉末サンプルを使用）の向上は、パラメータ抽出の精度を犠牲にするというコストを伴うことを意味する。なぜなら、異なるパラメータセットが類似のシミュレーションスペクトルを生み出す可能性があるからである。

制約と失敗モード

RuBr$_3$における磁気励起を理解するという問題は、いくつかの厳しい現実的な制約によって非常に困難になっている。

  • 物理的制約:

    • 材料の複雑性: RuBr$_3$は3層ハニカム構造（図1）を持ち、面内相互作用に加えて層間磁気相互作用($J_p$)を導入する。これにより、ハミルトニアンのパラメータ数が増加し、システムをモデル化することがより複雑になる。
    • 競合する磁気相互作用: この材料はジグザグ反強磁性秩序を示し、非キタイエフ相互作用（ハイゼンベルク $J_1, J_2, J_3$ および対角線外 $\Gamma, \Gamma'$）が重要であり、キタイエフ項 $K$ と競合するか、あるいはそれを支配していることを示唆している。この複雑な相互作用により、各相互作用の個々の寄与を分離し定量化することが困難になる。
    • 温度依存性ダイナミクス: RuBr$_3$の磁気励起は、スペクトル重みのシフトやエネルギーギャップの閉鎖を含む複雑な温度依存性を示す。広範囲の温度にわたるこれらの変化を正確に捉えることは重要であるが、実験的および解析的な負担を増大させる。
    • スピン軌道結合: 結合方向依存性相互作用は、Ru$^{3+}$イオンにおける強いスピン軌道結合に由来し、擬スピン1/2記述を必要とする。この量子力学的な複雑さにより、古典的近似（線形スピン波理論など）は本質的に限界がある。

  • 計算上の制約:

    • 線形スピン波理論（LSWT）の限界: スペクトルシミュレーションに使用されるLSWTは、「擬スピン1/2系には完全に適用できない」と論文は指摘している。これは、理論的枠組み自体が磁気励起の量子的な性質を正確に捉える上で固有の限界があり、実験結果との不一致につながる可能性があることを意味する。
    • 高次元パラメータ空間: ハミルトニアンは多数の交換パラメータ（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）を含む。論文は、「利用可能なデータから複数の交換パラメータを同時に推定することは不可能である」と明記している。これにより、研究者は現象論的な調整に頼らざるを得ず、一意で決定的なパラメータセットではなく、「2つの極端な組み合わせ」を提示することになる。
    • 粉末平均: 粉末非弾性中性子散乱データを使用することは、実験的にはよりアクセスしやすいが、本質的にすべての結晶学的配向を平均化する。この「粉末平均スペクトルから交換パラメータを推定することの困難さ」は、異なるモデルを区別するのに役立つ可能性のある重要な異方性情報が失われることを意味する。

  • データ駆動型制約:

    • 不完全なフォノン減算: 磁気信号を抽出するために、フォノン寄与は高温データから推定され、減算される。しかし、「フォノンの非調和性と、温度依存性がほとんどない背景のために、減算は完全ではない」。これはノイズと不確実性を導入し、微妙な磁気特徴を不明瞭にする可能性がある。
    • 巨視的特性推定の感度: ワイス温度という主要な巨視的特性の推定は、「スピン軌道結合励起状態からのヴァン・フローク常磁性様寄与に非常に敏感」である。これにより、精密なパラメータ決定のための信頼性の低い指標となり、理論モデルの検証をさらに複雑にする。
    • 粉末INSの限定された分解能: 粉末INSは、その性質上、励起の運動量平均ビューを提供する。論文は、弱い励起の場合、「波ベクトルとエネルギー依存性は粉末非弾性散乱スペクトルでは明確ではない」と述べており、競合する理論モデルを区別するのに重要なスピンダイナミクスの微細な詳細を解決することを困難にしている。
Figure 1. Magnetic structure of RuBr3. Ru atoms form a three-layered honeycomb structure with a crystallographic unit cell indicated by thin black lines. The magnetic moments of the Ru atoms form a zigzag antiferromagnetic structure with the unit cell indicated by thin red lines. JX , JY , and JZ represent bond-dependent anisotropic nearest-neighbour magnetic interactions. J2, J3, and Jp represent the next nearest neighbour, third nearest neighbour magnetic interactions within the honeycomb plane, and interplane magnetic interactions, respectively

なぜこのアプローチなのか

選択の必然性

率直に言って、従来の「SOTA」手法（標準的なCNN、基本的な拡散モデル、またはTransformerなど）がこの特定の問題に対して検討され、その後却下されたという考えは、科学分野の文脈においてはやや誤解を招く表現である。本論文は凝縮系物理学の分野で展開されており、磁気励起を調査するための「SOTA」手法は、非弾性中性子散乱のような確立された実験技術と線形スピン波理論のような理論的枠組みである。機械学習モデルは、他の分野では強力であるが、RuBr3のような材料の基本的な磁気相互作用やスピンダイナミクスをこのレベルで直接探求するには、単に適用できない。

著者らが非弾性中性子散乱と線形スピン波理論（LSWT）を組み合わせるという選択は、単なる好みではなく、問題の性質を考えると必然であった。中心的な目的は、キタイエフモデルの候補材料であるRuBr3の磁気励起を実験的に観測し、理論的に解釈することである。非弾性中性子散乱は、磁性材料におけるスピン励起（マグノンまたは分数化励起）のエネルギーと運動量依存性を直接測定するための主要な実験技術である。それは、材料内の動的な磁気相関への直接的な窓を提供する。

著者らが従来の（凝縮系物理学の文脈における）「SOTA」手法では単純な解釈には不十分であると認識した正確な瞬間は、おそらくRuBr3の初期特徴付けから来ている。論文の要旨と導入部は、キタイエフモデルが厳密な量子スピン液体基底状態を予測する一方で、RuBr3は「$T_N$以下のジグザグ反強磁性秩序」を示すと強調している。この観測は、スピン液体につながる理想的な、純粋なキタイエフモデルではRuBr3を記述するには不十分であることを直ちに示唆している。磁気秩序の存在は、拡張キタイエフ・ハイゼンベルクハミルトニアンに適用される線形スピン波理論が提供するものである、秩序相における集団的スピン励起（マグノン）を記述できる理論的枠組みを必要とする。したがって、この認識は計算アルゴリズムの失敗に関するものではなく、材料の物理的特性が最も単純な理論的理想から逸脱したため、より包括的な実験的および理論的アプローチを要求したというものであった。

比較優位性

非弾性中性子散乱実験と線形スピン波理論計算を組み合わせた選択されたアプローチは、微視的な磁気相互作用と巨視的に観測可能なスピンダイナミクスとの間に直接的かつ物理的に解釈可能なつながりを提供することにより、質的な優位性を提供する。

  1. スピンダイナミクスの直接プローブ: 非弾性中性子散乱は、スピン相関関数のフーリエ変換に比例する動的構造因子$S(\mathbf{Q}, E)$を直接測定するため、質的に優れている。これは、空間（運動量$\mathbf{Q}$）と時間（エネルギー$E$）におけるスピンの揺らぎを直接プローブすることを意味する。間接的なプローブとは異なり、中性子散乱は、励起分散関係（運動量とともに励起エネルギーがどのように変化するか）やスペクトル重みを含む、磁気励起スペクトルの完全な画像を提供する。この直接性は、バルク磁性材料の研究において他の技術では比類がない。
  2. LSWTを通じた物理的解釈可能性: 線形スピン波理論は、注意深く構築されたハミルトニアン（式2など）に適用されると、基本的な交換パラメータ（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）の抽出を可能にする。これは構造的な利点を提供する。観測された実験スペクトルを、材料の磁性を支配する微視的な相互作用に結びつける。本論文は、2つの異なるパラメータセット（J1-K-Γ-Γ'およびJ1-K-J2-J3モデル）からのシミュレーションスペクトルを実験データ（図3）と比較することで、これを明示的に実証している。J1-K-J2-J3モデルは、「より良い」一致を示しており、RuBr3における相互作用の複雑な相互作用を捉える能力が優れていることを示唆している。これは、$O(N^2)$から$O(N)$へのメモリ複雑性の削減に関するものではなく、実験的現実を正確に記述する堅牢で物理的に根拠のあるモデルを提供することに関するものである。観測された実験的「強く分散する磁気励起」と「分散しない励起」を再現する能力は、主要な質的利点である。

制約との整合性

選択された方法は、RuBr3のようなキタイエフモデル候補、特にジグザグ反強磁性秩序を示すものの研究という暗黙の制約と完全に整合している。

  1. 磁気励起の実験的特徴付け: 主要な制約は、磁気励起を実験的に特徴付けることである。非弾性中性子散乱は、スピン励起のエネルギーと運動量を直接測定するため、これにユニークに適している。粉末サンプル（単結晶とは対照的に）の使用は実用的な制約であり、解析ではシミュレーション中に全立体角を平均化することでこれを考慮している。
  2. 根底にある磁気相互作用の特定: 重要な要件は、材料の磁気挙動を駆動し、観測されたジグザグ反強磁性秩序を安定化させる可能性のある様々な磁気相互作用（キタイエフ、ハイゼンベルク、対角線外）を特定し定量化することである。拡張キタイエフ・ハイゼンベルクハミルトニアンに適用される線形スピン波理論は、まさにこの目的のために設計されている。これにより、著者らは実験データを理論モデルに適合させ、これらの交換パラメータの値を抽出することができる。この「実験的観測と理論的モデリングの結婚」は、配位子置換がこれらの相互作用にどのように影響するかを理解するために不可欠であり、これは論文の中心的なテーマである。「アニオン置換が磁気相互作用に与える影響を明確にする」という論文の目標は、この組み合わせアプローチによって直接満たされる。

代替案の却下

本論文は、単純な理論モデルを暗黙のうちに却下し、他の特徴付け技術の限界を強調しており、無関係な「SOTA」機械学習アプローチを明示的に却下しているわけではない。

  1. 理想キタイエフモデルの却下: 最も重要な暗黙の却下は、RuBr3の完全な記述としての理想的な強磁性キタイエフモデルの却下である。論文は、「理論的には、理想的な強磁性キタイエフモデルの磁気励起は波ベクトル依存性が弱いという特徴がある」と述べているが、彼らの実験結果は「強く分散する磁気モード」と「ジグザグ反強磁性秩序」を示している。この明確な不一致は、純粋なキタイエフモデルでは不十分であることを意味する。したがって、著者らは、観測された磁気秩序を安定化させ、観測された分散励起を再現するために必要な、追加のハイゼンベルク（$J_1, J_2, J_3$）および対角線外（$\Gamma, \Gamma'$）項を含む拡張ハミルトニアン（式2）を採用する。これは、単純な見方を超えた、理論モデルの改良である。
  2. 他の特徴付け手法の限界: 明示的な却下ではないが、論文は、$\alpha$-RuCl$_3$に関する以前の研究を議論する際に、導入部でラマンスペクトル、核磁気共鳴、X線散乱実験などの他の技術に言及している。これらの手法は貴重な補完情報（例：特定のフォノンモード、局所スピン環境、電子構造）を提供するが、非弾性中性子散乱が提供するような、スピン励起に関する直接的で運動量分解された情報を提供するものではない。例えば、ラマンスペクトルは2マグノン様励起を検出できるが、中性子散乱は完全な分散を提供する。
  3. 機械学習SOTAの無関係性: 前述のように、敵対的生成ネットワーク（GAN）や拡散モデルのようなアプローチは、画像生成、データ合成、または大規模データセットにおける複雑なパターン認識などのタスクを目的としている。それらは、特定の物理的プローブと理論的枠組みを必要とする量子材料におけるスピンダイナミクスの直接的な実験測定や、微視的な交換パラメータの導出には、根本的に不向きである。この文脈でのそれらの適用は、失敗した代替案ではなく、カテゴリエラーであろう。
Figure 5. Inelastic neutron scattering spectrum simulated from (a) the J1–K–Γ–Γ′ and (b) the J1– K–J2–J3 models by using linear spin wave theory. The parameters used for the simulations are (a) J1 = −1.8, K = −7.2, Γ = 10.5, Γ′ = −2.5 meV and (b) J1 = 1.5, K = −8.1, J2 = 0.8, J3 = 5.8, and Γ′ = −0.16 meV. Interplane interactions of Jp = 0.15 meV are adopted in both models

数学的・論理的メカニズム

マスター方程式

RuBr$_3$における磁気励起の分析の根底にある基本的な数学的枠組みは、材料内の様々な磁気相互作用を記述するハミルトニアンである。非弾性中性子散乱データを解釈するために重要な線形スピン波理論計算は、このハミルトニアンに基づいている。中心的な方程式は次のように提示されている。

$$ H = \sum_{NN} S_i J_{ij} S_j + J_2 \sum_{NNN} S_i \cdot S_j + J_3 \sum_{3NN} S_i \cdot S_j + J_p \sum_{interplane} S_i \cdot S_j \quad (2) $$

このハミルトニアンの重要な構成要素は、最近接結合異方性相互作用、$J_{ij}$であり、これは行列で表される。$z$-結合の場合、この行列は次の形式をとる。

$$ J_z = \begin{pmatrix} J_1 & \Gamma & \Gamma' \\ \Gamma & J_1 & \Gamma' \\ \Gamma' & \Gamma' & J_1 + K \end{pmatrix} \quad (3) $$

同様の行列、$J_x$および$J_y$は、対角線要素を循環的に置換して、$J_1+K$項を$x$または$y$方向のいずれかに配置することによって構築される。これは、キタイエフ相互作用の結合方向依存性を反映している。

項ごとの解剖

これらの各方程式の要素を分解して、その役割を理解しよう。

式(2)より – 全ハミルトニアン:

  • $H$: これは磁気系の全ハミルトニアンを表す記号である。

    • 数学的定義: 量子力学的な演算子であり、系の全エネルギーに対応する。その固有値は可能なエネルギー状態を与え、その固有状態は対応する物理的配置を記述する。
    • 物理的/論理的役割: これはRuBr$_3$の磁気挙動を支配する中心的な方程式である。その励起を解くことで、著者らは材料の磁気スペクトルを予測できる。
    • なぜ加算なのか: 系全体のエネルギーは、すべての個々の相互作用エネルギーの合計である。ハミルトニアンの各項は、異なる種類の磁気結合を表し、それらの複合効果が全体の磁気特性を決定する。

  • $\sum_{NN}$: これはすべての最近接結合（NN）スピンペアに対する総和演算子である。

    • 数学的定義: ハニカム格子上の直接隣接するスピン$S_i, S_j$の各ペアに対して、後続の項を合計するように指示する。
    • 物理的/論理的役割: 最も直接的で、しばしば最も強い磁気相互作用が材料全体で考慮されることを保証する。

  • $S_i$: これはサイト$i$におけるスピン演算子を表す。

    • 数学的定義: 特定の原子サイト（この場合はRu$^{3+}$の擬スピン1/2）における電子の量子力学的なスピン角運動量を表すベクトル演算子。$S_i^x, S_i^y, S_i^z$の成分を持つ。
    • 物理的/論理的役割: これらは基本的な磁気自由度である。これらのスピン間の相互作用が、観測される磁気励起を生み出す原因である。

  • $J_{ij}$: これはスピン$S_i$と$S_j$間の最近接結合異方性相互作用行列である。

    • 数学的定義: スピン$S_i$と$S_j$の$x, y, z$成分間の結合の強さと性質を定義する要素を持つ3x3行列（式3に示すように）。
    • 物理的/論理的役割: この項はキタイエフモデルの核心であり、複雑で、結合方向依存性があり、異方性のある交換相互作用を捉え、隣接するRuスピン間で作用する。エキゾチックなスピン液体物理学を理解する上で不可欠である。相互作用は単純な内積ではなく、スピン成分と結合方向の特定の関数であるため、行列形式が必要である。

  • $S_i \cdot S_j$: これはスピン$S_i$と$S_j$間の等方性内積相互作用を表す。

    • 数学的定義: 2つのスピンベクトル、$S_i^x S_j^x + S_i^y S_j^y + S_i^z S_j^z$の内積。
    • 物理的/論理的役割: これは等方性（方向非依存性）の従来のハイゼンベルク型相互作用を記述する。より離れたスピン間の相互作用に使用され、通常、結合依存性はそれほど重要ではない。

  • $J_2$: これは次最近接結合（NNN）等方性磁気相互作用の結合定数である。

    • 数学的定義: NNNペアに対する$S_i \cdot S_j$相互作用の強度をスケーリングするスカラー係数。
    • 物理的/論理的役割: 2サイト離れたスピン間の等方性磁気結合の強さを定量化する。この項は、ジグザグ反強磁性秩序のような特定の磁気秩序を安定化させる役割を果たす。

  • $\sum_{NNN}$: これはすべての次最近接結合ペアに対する総和演算子である。

    • 数学的定義: 2サイト離れたスピン$S_i, S_j$の各ペアに対して、後続の項を合計する。
    • 物理的/論理的役割: すべてのNNN等方性相互作用がエネルギー計算に含まれることを保証する。

  • $J_3$: これは第三最近接結合（3NN）等方性磁気相互作用の結合定数である。

    • 数学的定義: 3NNペアに対する$S_i \cdot S_j$相互作用の強度をスケーリングするスカラー係数。
    • 物理的/論理的役割: 3サイト離れたスピン間の等方性磁気結合の強さを定量化する。本論文は、観測されたジグザグ反強磁性秩序を安定化させる上でその重要性を強調している。

  • $\sum_{3NN}$: これはすべての第三最近接結合ペアに対する総和演算子である。

    • 数学的定義: 3サイト離れたスピン$S_i, S_j$の各ペアに対して、後続の項を合計する。
    • 物理的/論理的役割: すべての3NN等方性相互作用が含まれることを保証する。

  • $J_p$: これは層間等方性磁気相互作用の結合定数である。

    • 数学的定義: 異なる層にあるスピン間の$S_i \cdot S_j$相互作用の強度をスケーリングするスカラー係数。
    • 物理的/論理的役割: 材料の3次元性を考慮して、異なる層にあるスピン間の等方性磁気結合の強さを定量化する。

  • $\sum_{interplane}$: これはすべての層間ペアに対する総和演算子である。

    • 数学的定義: 異なる層にあるスピン$S_i, S_j$の各ペアに対して、後続の項を合計する。
    • 物理的/論理的役割: すべての層間等方性相互作用が含まれることを保証する。

式(3)より – 最近接結合相互作用行列$J_z$:

  • $J_1$: これは等方性ハイゼンベルク交換結合定数である。

    • 数学的定義: $J_{ij}$行列の対角線上に現れるスカラー係数で、$S_x S_x$, $S_y S_y$, $S_z S_z$相互作用の強度を表す（キタイエフ項を除く）。
    • 物理的/論理的役割: 最近接結合交換相互作用の従来の、方向非依存性の部分を表す。磁気相互作用の基本的な構成要素である。

  • $K$: これはキタイエフ相互作用項である。

    • 数学的定義: $J_{ij}$行列の対角線要素の一つ、特に$z$-結合に対する$S_z S_z$成分に加算されるスカラー係数。
    • 物理的/論理的役割: キタイエフモデルの定義的な特徴であり、高度に異方性があり、結合方向依存性のあるイジング様相互作用を表す。エキゾチックなスピン液体物理学の原因である。著者らは、キタイエフ相互作用が特定の結合方向沿いの等方性ハイゼンベルク項に対する追加の異方性成分として作用するため、加算を使用している。

  • $\Gamma$: これは対称対角線外交換項である。

    • 数学的定義: $J_{ij}$行列の対角線外の位置に現れるスカラー係数（例：$S_x S_y$と$S_y S_x$を結合する）。
    • 物理的/論理的役割: スピン成分の交換に対して対称な異方性交換相互作用を表す。スピン軌道結合に由来し、キタイエフ相互作用と競合または強化する可能性があり、磁気基底状態に大きく影響する。

  • $\Gamma'$: これは反対称対角線外交換項である。

    • 数学的定義: $J_{ij}$行列の他の対角線外の位置に現れるスカラー係数（例：$S_x S_z$と$S_z S_x$を結合する）。
    • 物理的/論理的役割: 反対称的な性質を持つ別の種類の異方性交換を表す。$\Gamma$と同様に、スピン軌道結合に由来し、磁気異方性および全体の磁気相図の形成に役割を果たす。

ステップバイステップの流れ

単一の抽象的なデータポイント、すなわちRuBr$_3$材料の磁気状態を表すものを、一連の計算段階を経て移動すると想像してほしい。

  1. 構造ブループリント入力: まず、RuBr$_3$の物理的構造がシステムに入力される。これには、層状ハニカム格子におけるRu原子の正確な配置が含まれ、最近接結合、次最近接結合、第三最近接結合、および層間接続が定義される。各Ru原子には、擬スピン1/2演算子$S_i$が割り当てられる。

  2. ハミルトニアン組み立て: 次に、全ハミルトニアン（式2）が構築される。格子内のすべての結合に対して、適切な相互作用項が選択される。最近接結合の場合、結合方向に基づいて、適切な$J_{ij}$行列（式3のような$J_z$、またはその$J_x, J_y$に対応するもの）が選択され、$J_1, K, \Gamma, \Gamma'$の値が組み込まれる。さらなる近接結合および層間結合の場合、スカラー項$J_2, J_3, J_p$が適用される。このステップは、すべての磁気相互作用を記述する包括的なエネルギー関数を構築する。

  3. 古典的基底状態決定: 量子揺らぎを研究する前に、古典的基底状態、すなわち最低エネルギーのスピン配置が、Luttinger-Tisza法を用いて特定される。ハミルトニアンは古典的に扱われ、システムは各$S_i$のスピン大きさを一定に保ちながら、全エネルギー（式2）を最小化するスピン配置を探索する。これにより、ジグザグ反強磁性秩序のような平衡スピン配置が得られ、その後の量子計算の参照点となる。

  4. 運動量空間変換: 次に、実空間ハミルトニアンがフーリエ変換を用いて運動量空間に変換される。これにより、サイト固有のスピン演算子が運動量依存演算子に変換され、相互作用項が運動量依存結合関数（式4の$J_{AAk}, J_{ABk}$など）になる。この変換は、格子の並進対称性を利用して問題を単純化する。

  5. Holstein-Primakoff近似: 量子励起を記述するために、スピン演算子はボソン生成・消滅演算子を用いて近似される。これはHolstein-Primakoff変換であり、古典的基底状態からの小さな偏差（マグノン）に対してスピンダイナミクスを線形化する。これにより、複雑なスピンハミルトニアンがこれらのボソン演算子に関する二次形式に効果的に変換される。

  6. Bogoliubov-de Gennes行列形成: ジグザグ反強磁性秩序の存在により、システムは実質的に4つの磁気サブ格子によって記述される。前の変換の後、二次ハミルトニアンは、各波ベクトル$k$に対してBogoliubov-de Gennes形式の8x8行列にキャストされる。この行列は、ボソン励起の結合ダイナミクスを数学的に記述する。

  7. 固有値抽出（励起スペクトル）: この8x8 Bogoliubov-de Gennes行列が、各波ベクトル$k$に対して対角化される。得られた固有値は、励起エネルギー（マグノン分散関係）に直接対応し、固有ベクトルはこれらのマグノンモードの性質を記述する。このステップは、運動量に対するエネルギーの変化を示す理論的な非弾性中性子散乱スペクトルを生成する。

  8. 動的構造因子計算: 固有値と固有ベクトルから、動的構造因子が計算される。この量は、非弾性中性子散乱実験で測定される強度に直接比例し、理論と実験との直接的なつながりを提供する。

  9. 粉末平均と分解能畳み込み: 最後に、計算されたスペクトルは、運動量空間のすべての可能な方向について平均化され（粉末非弾性中性子散乱実験をシミュレートするため）、その後、実験的な波ベクトルとエネルギー分解能で畳み込まれる。これにより、図に示される実験データと直接比較できる理論スペクトルが生成される。

最適化ダイナミクス

このメカニズムにおける「学習」または「収束」は、主に機械学習に見られるような自律的な反復最適化アルゴリズムではなく、実験データに対するパラメータ適合とモデル検証のプロセスである。その展開は以下の通りである。

  1. 初期パラメータ推測: プロセスは、交換パラメータ（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）の初期セットから始まる。これらの値は、ab initio計算、類似材料に関する以前の研究、または教育的な推測から得られる可能性がある。

  2. 古典的基底状態最小化: 各パラメータセットに対して、Luttinger-Tisza法が古典的基底状態を見つけるために使用される。これはそれ自体が最適化ステップである。この方法は、スピン magnitudesの固定という制約を遵守しながら、全エネルギー（式2）を最小化するスピン配置を反復的に探索する。これは実質的にエネルギーランドスケープをナビゲートして大域的最小値を見つける。

  3. フォワードモデル計算: 古典的基底状態が確立されると、線形スピン波理論が適用され、マグノン励起スペクトルと動的構造因子が決定論的に計算される。これはモデルの「フォワードパス」であり、現在のパラメータに基づいた理論的予測を生成する。

  4. 比較と不一致評価: 計算された理論スペクトルは、実験的な非弾性中性子散乱データ（例：図2の強度マップまたは図3と4の積分強度）と比較される。 「損失」または「不一致」は、理論的予測と実験的観測との差である。著者らは、ピーク位置、分散、および全体的なスペクトル重みのような特徴において、質的および量的な一致を探している。

  5. 現象論的パラメータ調整: 次に、著者らは交換パラメータを現象論的に調整する。これは典型的な意味での勾配ベースの最適化ではなく、反復的で人間主導のプロセスである。彼らは、特定のジグザグ反強磁性秩序を安定化させるために$J_3$を増加させたり、スピン液体特性に影響を与えるために$K$を調整したりするかもしれない。例えば、分散型および分散しない励起のような特定の特徴を再現するために、2つの「極端な組み合わせ」のパラメータ（$\Gamma$優位の1つと$J_3$優位のもう1つ）を試したと述べている。

  6. 多制約検証: 非弾性中性子散乱スペクトルを超えて、選択されたパラメータは他の巨視的特性とも検証される。例えば、ワイス温度（式1）と磁気モーメントのキャンティング角は、パラメータから計算され、実験値と比較される。これらは、パラメータが満たす必要がある追加の「損失関数」または制約として機能し、実行可能なパラメータ空間を絞り込むのに役立つ。著者らはワイス温度推定の感度を指摘しており、これがフィッティングの繊細な部分であることを示唆している。

  7. 「最良適合」への収束: すべての利用可能な実験データにわたって全体的な一致を最もよく提供するパラメータセットが見つかると、プロセスは収束する。この「最良適合」は、観測された現象を説明する根底にある磁気相互作用に関する著者らの解釈を表す。したがって、最適化は、理論的マグノン分散と強度を実験的観測に一致するように形成するパラメータの探索であり、反復的な比較と調整を通じて根底にある磁気相互作用を効果的に「学習」する。

結果、限界、結論

実験設計とベースライン

RuBr$_3$における磁気励起を厳密に調査するために、研究者らはスピンダイナミクスをプローブする強力な技術である粉末非弾性中性子散乱を採用した。研究対象の材料はRuBr$_3$の新しい多形であり、R3空間群に結晶化し、$\alpha$-RuCl$_3$の低温相と構造的に同一の3層ハニカム構造を形成した。RuBr$_3$を研究する動機は、配位子イオン（臭素を塩素に置換）がキタイエフ、ハイゼンベルク、および対角線外磁気相互作用の複雑なバランスを実験的に調整できるという考えに由来する。以前の研究では、RuBr$_3$が、より小さなバンドギャップ、低減された抵抗、抑制された磁化率、異なるワイス温度、およびより大きな磁気モーメント傾斜角など、$\alpha$-RuCl$_3$とは異なる磁気特性を示すことがすでに示唆されていた。

実験は、さまざまなエネルギーと運動量の範囲にわたって包括的なビューを確保するために、J-PARCのAMATERAS、ANSTOのPELICAN、およびJRR-3のGPTASの3つの異なる中性子分光器を使用して実施された。多結晶RuBr$_3$サンプルは、立方アンビル高圧装置を使用して合成された。9.5gのサンプルを円筒形（直径15mm、高さ15mm）に成形し、ヘリウム交換ガスとともにアルミニウム缶に封入した。非弾性中性子散乱スペクトルは、さまざまな温度（10K、25K、45K、100K、300K）で収集された。入射中性子エネルギー（$E_i$）は20.95、9.70、5.57、3.61 meVに設定され、さらに42.17、15.19、7.75 meVで高エネルギー分解能測定が行われた。

データ解析の重要なステップは、磁気寄与を分離することであった。これは、主に300Kデータから得られたフォノン寄与を推定・減算し、温度依存性ボーズ因子、$1 + n(T) = (1 - e^{-E/(k_B T)})^{-1}$を使用して強度を補正することによって達成された。Utsusemiソフトウェアスイートがこの解析に使用された。

この研究における「犠牲者」またはベースラインモデルは、主に理論的なものであった。すなわち、キタイエフ・ハイゼンベルクハミルトニアンにおける交換パラメータの異なる組み合わせを表すJ1-K-Γ-Γ'およびJ1-K-J2-J3モデルである。これらのモデルは、期待される非弾性中性子散乱スペクトルをシミュレートするために使用され、実験データとの直接比較を可能にした。さらに、RuBr$_3$の実験結果は、$\alpha$-RuCl$_3$やNa$_2$IrO$_3$のような他のキタイエフ候補材料の既知の挙動と、暗黙的および明示的に比較され、RuBr$_3$のユニークな特性を強調するための経験的ベースラインとして機能した。

証拠が証明すること

RuBr$_3$の非弾性中性子散乱実験は、反強磁性ゾーン中心に集中した強く分散する磁気励起として現れる強い反強磁性相互作用の存在を決定的に証明した。ネール温度($T_N = 34$ K、中性子回折で確認)以下、特に10Kおよび25Kでは、顕著な分散励起が観測された。これらの励起は、波ベクトル0.60 Å$^{-1}$および1.55 Å$^{-1}$でピークを示し、磁気反射と一致しており、それらの磁気起源を強く示唆している。10Kで1.5 meVのエネルギーギャップが特定されたが、これは25K以上ではほぼゼロに減少した。

温度が上昇するにつれて、これらの励起の強い波ベクトル依存性は45Kまで持続した。しかし、100K（約$3T_N$）でスペクトル重みがゼロ波ベクトル（$\Gamma$点）にシフトし、その後300Kで完全に消失するという顕著なシフトが発生した。$\Gamma$点へのこのシフトは、システム内に強磁性相互作用が存在することの明確な証拠として解釈される。逆に、ブリルアンゾーン境界付近で観測される磁気励起の堅牢性は、観測されたジグザグ反強磁性秩序を安定化させるために不可欠な強い反強磁性相互作用を示唆している。

他のキタイエフ候補と比較して、RuBr$_3$の磁気励起の温度堅牢性は、$\alpha$-RuCl$_3$よりもNa$_2$IrO$_3$に近いことがわかった。$\alpha$-RuCl$_3$の励起は$T_N$をわずかに超えると$\Gamma$点にシフトするが、RuBr$_3$のシフトは約$3T_N$というより高い温度で発生する。これは、$\alpha$-RuCl$_3$とは異なる相互作用のバランスを示唆している。

実験スペクトルは、2つの線形スピン波理論モデル（J1-K-Γ-Γ'およびJ1-K-J2-J3）からのシミュレーションと比較された。両モデルとも、分散型およびほぼ分散しないモードの両方を含む、スピン励起の特性を再現できた。しかし、大きな第三最近接結合等方性磁気相互作用（$J_3$）を組み込んだJ1-K-J2-J3モデルは、実験的な波ベクトルとエネルギー移動依存性との一致がより良好であった（図3(a)および3(b)を参照）。これは、ジグザグ反強磁性秩序を安定化させ、観測された強い分散を説明するために、大きな最近接結合対称対角線外（$\Gamma$）または第三最近接結合等方性磁気相互作用（$J_3$）が不可欠であることを強く示唆している。

実験的なワイス温度と磁気モーメントのキャンティング角は、モデル予測との間にいくつかの不一致を示したが、論文は最終的に、J1-K-J2-J3モデルがより大きな$J_3$項を持つ方がより現実的であると結論付けている。この結論は、粉末平均データから交換パラメータを正確に推定することの難しさや、ワイス温度計算のさまざまな寄与に対する感度にもかかわらず、主に非弾性中性子散乱スペクトルへのより良い適合に基づいている。

最後に、この研究は、臭素置換が磁気相互作用に著しく影響を与えるという確固たる証拠を提供している。波解析は、強磁性キタイエフ項が持続するが、反強磁性相互作用（特に$\Gamma$または$J_3$）が増強されることを示唆している。この増強は、RuBr$_3$の基底状態をジグザグ反強磁性相にさらに深く押し込み、理想的な強磁性キタイエフスピン液体状態から遠ざける。この効果は、より強いd-p混成と配位子原子を介した間接ホッピングの増強に起因すると考えられる。

限界と今後の方向性

提示された説得力のある証拠にもかかわらず、この研究は、実験的アプローチと理論的モデリングに固有のいくつかの限界を認識している。主な限界は、粉末平均非弾性中性子散乱スペクトルの使用に由来する。この平均化により、交換パラメータの正確な推定が本質的に困難になり、2つの異なる理論モデル（J1-K-Γ-Γ'およびJ1-K-J2-J3）によって再現された類似のスペクトル特徴によって証拠立てられている。モデル化に使用される線形スピン波近似は、擬スピン1/2系には完全に適用できない可能性があるが、磁気相互作用の概算値を提供する。さらに、実験的に導出されたワイス温度とモデルによって予測されたものとの間に顕著な不一致があり、著者らは推定の感度と、スピン軌道結合励起状態からの潜在的なヴァン・フローク常磁性様寄与に起因すると考えている。同様に、モデルで予測されたキャンティング角は実験値と完全に一致しなかったが、著者らは$J_3$におけるさらなる異方性を組み込むことでこれが改善される可能性があると示唆している。決定的に、論文は「その波ベクトルとエネルギー依存性が粉末非弾性散乱スペクトルでは明確ではないため、弱い励起がキタイエフ相互作用から誘起されると結論付けることはできなかった」と述べており、キタイエフ寄与を決定的に特定することにおける重大なギャップを示唆している。さらに、Br置換による反強磁性相互作用増強の正確な起源は「不明」であり、根本的な電子メカニズムに関するさらなる調査が必要である。

将来を見据えると、これらの発見は、将来の研究と開発のためのいくつかのエキサイティングな道を開く。

  • 決定的な証拠のための単結晶研究: 最も直接的かつ影響力のある将来の方向性は、RuBr$_3$の高品質単結晶を用いた非弾性中性子散乱実験を実施することであろう。これにより、粉末平均の曖昧さが解消され、交換パラメータのより正確な決定と、すべての励起、特に捉えどころのないキタイエフ寄与の波ベクトルとエネルギー依存性についてのより明確な理解が可能になる。これにより、キタイエフ相互作用を完全に特徴付けるために必要な、決定的な、否定できない証拠が得られる可能性がある。

  • 高度な理論モデリング: LSWTの限界に対処するために、将来の研究は、より洗練された理論モデルの開発と適用に焦点を当てるべきである。擬スピン1/2系および強相関材料により適した、厳密対角化、密度行列繰り込み群（DMRG）、または量子モンテカルロシミュレーションのような技術は、スピンダイナミクス、キャンティング角、およびワイス温度のより正確な記述を提供し、理論と実験の間のギャップを埋めることができる。

  • 配位子置換メカニズムの解明: 本論文はBr置換の重大な影響を強調しているが、反強磁性相互作用を増強する正確なメカニズムは不明であると指摘している。将来の研究では、高度なab initio計算とターゲットを絞った実験プローブ（例：共鳴非弾性X線散乱、X線吸収分光法）を組み合わせることで、配位子イオンの変化がd-p混成、軌道交換経路、そして最終的には特定の磁気相互作用にどのように影響するかを正確にマッピングできる可能性がある。これにより、望ましい特性を持つキタイエフ材料を設計するための予測フレームワークが得られる可能性がある。

  • RuX$_3$ファミリーと多形の探索: RuX$_3$（X = Br, I）は複数の多形を形成するため、RuBr$_3$およびRuI$_3$の他の構造相の体系的な調査は、新しい磁気基底状態や、キタイエフと他の相互作用の異なるバランスを明らかにする可能性がある。この比較アプローチは、エキゾチックなキタイエフスピン液体状態を支持する特定の構造モチーフまたは化学環境を特定できる可能性がある。

  • 外部摂動と相チューニング: 磁場、圧力、またはひずみなどの外部摂動を適用することは、RuBr$_3$における競合する相互作用のバランスを調整するための強力な方法である。これらの条件下でのスピンダイナミクスの調査は、相転移を明らかにし、競合する相互作用の性質に関するより深い洞察を提供する可能性がある。

  • 分数化励起の探索: 現在の研究はマグノン様励起を強調しているが、より高いエネルギー（例：12および15 meV）での弱い広がりを持つ励起の存在は、キタイエフスピン液体に特徴的な分数化励起の連続を示唆している可能性がある。将来の高分解能実験、特に単結晶での実験は、基底状態が秩序化している場合でも、マヨラナフェルミオンの兆候を特に探すべきである。

  • 層間結合と3D効果: ハミルトニアンには層間相互作用（$J_p$）が含まれているが、観測されたダイナミクスへのその役割は広く議論されていない。3層ハニカム構造を考慮すると、層間結合が磁気励起と全体的な基底状態に与える影響の詳細な調査は、特に純粋な2Dキタイエフ物理学からの逸脱を理解する上で価値がある可能性がある。

実験的改良から高度な理論モデリング、材料設計に至るまで、これらの多様な視点は、RuBr$_3$における相互作用の複雑な相互作用を完全に理解し、キタイエフスピン液体の探求を進める上で不可欠である。

Figure 2. a–e) Inelastic neutron scattering spectrum measured by using AMATERAS with an incident neutron energy of 20.95 meV at a) 10, b) 25, c) 45, d) 100 K and e) 300 K. Dispersive spin excitations were observed at 0.60 and 1.55˚A−1 up to the energy transfer of 15 meV. f–i) Two-dimensional colour maps of the magnetic contributions at f) 10, g) 25, h) 45 and i) 100 K estimated by subtracting the phononic contributions estimated from the 300 K data. Intensities are corrected by the temperature- dependent factor 1 + n(T ) = (1 − e−E/(kBT ))−1, where n(T ) represents a Bose factor Figure 3. Integrated scattering intensities at 10 K plotted as a function of the wavevector or energy transfer after the subtraction of the phonon contribution estimated from the 300 K data. Dashed and solid curves represent the simulated curves based on the J1–K–Γ-Γ′ and J1–K–J2–J3 models, respectively (see text for details). (a) Wavevector dependence of the intensities with integration ranges of [1.0, 2.0], [2.0, 3.0], [3.0, 4.0], [4.0, 5.0], [5.0, 6.0], and [6.0, 7.0] meV. The intensities are shifted for clarity. (b) Energy transfer dependence of the intensities with an integration range of [1.44, 1.74]˚A−1