RuBr3：Kitaev模型候选材料的磁激发

背景与学术渊源

起源与学术渊源

本文所研究的问题源于凝聚态物理学一个引人入胜的领域，特别是对量子自旋液体（QSL）态的探索与表征。这项探索始于2006年Alexei Kitaev提出的Kitaev模型理论预测[1]。该模型在蜂窝状晶格上定义了独特的键依赖的Ising相互作用，其显著之处在于可以精确求解，揭示出QSL基态。在这种奇异状态下，自旋即使在绝对零度下也不会有序排列；相反，它们保持纠缠并持续涨落，从而产生被称为Majorana费米子的分数化激发[7]。这些特性为容错量子计算带来了巨大的希望。

随后的挑战在于寻找能够承载这些难以捉摸的Kitaev相互作用的真实材料。大约在2009年，随着人们认识到某些过渡金属化合物（特别是具有$d^5$电子构型如Ir$^{4+}$和Ru$^{3+}$离子的化合物）中的强自旋-轨道耦合，可以通过涉及赝自旋-1/2态的机制有效地诱导这些键依赖的相互作用[8]，学术界取得了重大突破。这催生了“Kitaev材料”这一新前沿。在这些材料中，$\alpha$-RuCl$_3$迅速成为研究最广泛的候选材料[9-38]。

然而，先前对$\alpha$-RuCl$_3$的研究和方法存在一个根本性的局限性，即“痛点”：尽管该材料在较高温度下表现出一些Kitaev自旋液体行为的迹象，但其在某个临界温度（$T_N$）以下的基态是一种传统的锯齿形反铁磁有序。这表明存在其他磁相互作用（如Heisenberg和非对角项）与理想的Kitaev自旋液体态竞争并最终破坏其稳定性[31, 32, 39-43]。这些竞争相互作用阻碍了所需QSL性质的完全实现。本文的作者们受到理解和控制这些相互作用的需求的驱动。一个关键策略是改变磁性原子周围的配体离子，因为这些离子在介导产生键依赖相互作用的轨道交换机制中起着至关重要的作用。本文研究了RuBr$_3$作为一种新候选材料，其中用溴取代氯预计会改变这些磁相互作用的平衡，可能使材料更接近或更远离理想的Kitaev自旋液体态。

直观的领域术语

1. Kitaev模型 (Kitaev Model): 想象一种特殊的磁性棋盘格，其中两个相邻棋子之间的相互作用规则（吸引或排斥）完全取决于连接它们的直线的方向。例如，水平连接的棋子可能倾向于沿相同方向排列，而对角连接的棋子可能倾向于沿相反方向排列。Kitaev模型是一个描述这种系统的理论框架，它之所以出名，是因为它预测了一种“量子自旋液体”状态，在这种状态下，即使在最低温度下，磁体也不会冻结成固定的图案。
2. 量子自旋液体 (Quantum Spin Liquid, QSL): 想象一种普通液体，其中的水分子不断运动和相互作用。量子自旋液体与之类似，但不是水，而是电子微小内磁体（自旋）处于一种持续的、纠缠的、“液体状”状态。即使在绝对零度下，它们也不会沉淀成简单的、有序的磁图案，就像冰箱磁铁一样。相反，它们的集体行为是高度复杂的、量子力学的，并产生奇异的“半粒子”，称为Majorana费米子。
3. Majorana费米子 (Majorana Fermions): 如果将电子这样的普通粒子看作一个完整的人，那么Majorana费米子就像“半个人”。在量子自旋液体的背景下，电子的自旋可以有效地分裂成两个独立的、分离的部分，每一部分都表现得像一个Majorana费米子。这些“半粒子”之所以独特，是因为它们是自身的反粒子，并且它们对局部干扰的鲁棒性使其成为未来量子计算机的有希望的构建模块。
4. 键依赖的Ising相互作用 (Bond-dependent Ising Interactions): 想象一个由微小指南针组成的网格，每个指南针只能严格指向向上或向下（就像二进制开关一样）。“键依赖”意味着两个相邻指南针倾向于指向相同方向还是相反方向，完全取决于它们在网格上连接线的特定方向。例如，沿“南北”键，它们可能倾向于平行；但沿“东西”键，它们可能倾向于反平行。这种方向性偏好是Kitaev模型背后的核心思想。
5. 赝自旋-1/2态 (Pseudospin-1/2 State): 在某些复杂材料中，电子的实际自旋（始终为1/2）与其围绕原子的轨道运动纠缠在一起。这种组合产生了一个有效的磁矩，其行为如同一个简单的自旋-1/2粒子，但其“向上”和“向下”方向由晶体的结构决定，而不仅仅是简单的磁场。这是一种简化描述这些材料中复杂电子行为的巧妙方法。

符号表

符号	描述
$H$	磁系统的总哈密顿量
$S_i$	位点 $i$ 处的自旋算符
$J_{ij}$	最近邻各向异性相互作用矩阵
$J_1$	各向同性Heisenberg交换耦合常数
$K$	Kitaev相互作用项
$\Gamma$	对称的非对角交换项
$\Gamma'$	反对称的非对角交换项
$J_2$	次近邻（NNN）各向同性磁相互作用的耦合常数
$J_3$	第三近邻（3NN）各向同性磁相互作用的耦合常数
$J_p$	晶面间各向同性磁相互作用的耦合常数
$T_N$	Néel温度（磁有序温度）
$E_i$	入射中子能量
$S(\mathbf{Q}, E)$	动态结构因子（通过非弹性中子散射测量）
$\mathbf{Q}$	波矢（动量转移）
$E$	能量转移
$n(T)$	Bose因子
$k_B$	Boltzmann常数
$\theta_{ab}, \theta_c$	Weiss温度（面内和面外）

问题定义与约束

核心问题表述与困境

本文所解决的核心问题是精确表征RuBr$_3$（一种Kitaev量子自旋液体候选材料）中的磁相互作用，并理解这些相互作用如何驱动其基态偏离理想的Kitaev自旋液体。

输入/当前状态是RuBr$_3$这种材料，它是一种具有层状蜂窝状结构的新多晶型物，与广泛研究的Kitaev候选材料$\alpha$-RuCl$_3$同构。两种材料在各自的Néel温度（$T_N$）以下都表现出锯齿形反铁磁（AFM）有序。然而，RuBr$_3$与$\alpha$-RuCl$_3$相比表现出不同的磁性质，例如磁化率的抑制和不同的Weiss温度。先前对RuBr$_3$的从头算计算表明其磁相互作用与$\alpha$-RuCl$_3$相当，这与实验观测相矛盾。理想铁磁Kitaev模型的理论预测自旋激发具有较弱的波矢依赖性，但由于竞争相互作用，真实材料通常表现出更复杂的行为。

期望终点（输出/目标状态）是澄清阴离子取代（特别是用溴取代氯）对Kitaev候选材料中磁相互作用的影响。本文旨在通过非弹性中子散射（INS）实验研究RuBr$_3$的自旋动力学，以识别稳定观测到的锯齿形AFM有序的特定磁相互作用（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）。最终目标是为RuBr$_3$提供一个更现实的模型，该模型能准确重现其磁激发和宏观性质，从而解释其基态如何偏离理想的铁磁Kitaev自旋液体。

缺失环节或数学鸿沟在于定量确定RuBr$_3$的精确交换参数集，该参数集能够一致地解释其宏观磁性质（如Weiss温度和倾斜角）及其微观自旋动力学（磁激发谱）。虽然存在理论模型（例如，J$_1$-K-$\Gamma$-$\Gamma$'和J$_1$-K-J$_2$-J$_3$模型），但本文强调了从粉末平均实验数据中准确估计这些参数的困难。RuBr$_3$的理论预测与实验观测到的磁行为之间存在明显的差异，本研究试图通过详细的自旋动力学实验表征来弥合这一差距。

困扰先前研究者并成为本研究核心的痛苦权衡或困境是真实材料中期望的Kitaev相互作用与其他磁相互作用（Heisenberg、非对角、更远邻）之间固有的竞争。虽然Kitaev模型在实现自旋液体方面具有理论吸引力，但这些额外相互作用的存在通常会导致传统的磁有序，如RuBr$_3$中观察到的锯齿形AFM。研究人员面临着如何调节这些竞争相互作用以抑制传统有序或增强Kitaev物理学的困境。配体取代是一种策略，但其对相互作用精细平衡的确切影响是复杂的，不易预测，常常导致意想不到的后果或新的挑战。本文明确指出，“控制不仅Kitaev相互作用，而且Heisenberg和非对角磁相互作用对于实际材料的实际利用至关重要”，这强调了这一根本性权衡。此外，从粉末平均实验数据中唯一确定多个交换参数的困难意味着，通过提高实验可及性（使用粉末样品）是以牺牲参数提取的精度为代价的，因为不同的参数集可以产生相似的模拟谱。

约束与失效模式

理解RuBr$_3$中磁激发的难题因几个严峻的现实约束而变得异常困难：

• 物理约束:

  • 材料复杂性: RuBr$_3$具有三层蜂窝状结构（图1），除了面内相互作用外，还引入了晶面间磁相互作用（$J_p$）。这增加了哈密顿量中参数的数量，使得系统更难建模。
  • 竞争磁相互作用: 该材料表现出锯齿形反铁磁有序，表明非Kitaev相互作用（Heisenberg $J_1, J_2, J_3$ 和非对角 $\Gamma, \Gamma'$）很重要，并且与Kitaev项 $K$ 竞争甚至占主导地位。这种复杂的相互作用使得分离和量化每种相互作用的个体贡献变得具有挑战性。
  • 温度依赖动力学: RuBr$_3$中的磁激发表现出复杂的温度依赖性，包括谱权重转移和能量隙的闭合。在宽温度范围内准确捕捉这些变化至关重要，但这增加了实验和分析的负担。
  • 自旋-轨道耦合: 键依赖相互作用源于Ru$^{3+}$离子中的强自旋-轨道耦合，这需要赝自旋-1/2描述。这种量子力学复杂性使得经典近似（如线性自旋波理论）本身具有固有的局限性。

• 计算约束:

  • 线性自旋波理论（LSWT）的局限性: 本文指出，用于模拟谱的LSWT“不完全适用于赝自旋-1/2系统”。这意味着理论框架本身在准确捕捉磁激发量子性质方面存在固有局限性，可能导致与实验结果的差异。
  • 高维参数空间: 哈密顿量涉及多个交换参数（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）。本文明确指出“无法同时从可用数据中估计多个交换参数”。这迫使研究人员依赖现象学调整，并呈现“两个极端组合”，而不是一个独特、确定的参数集。
  • 粉末平均: 使用粉末非弹性中子散射数据虽然在实验上更容易获得，但本质上是对所有晶体学取向进行了平均。这种“从粉末平均谱中估计交换参数的困难”意味着关键的各向异性信息（有助于区分不同模型）丢失了。

• 数据驱动约束:

  • 声子减法不完整: 为了提取磁信号，声子贡献是根据高温数据估算并减去的。然而，“由于声子的非谐性和几乎不随温度变化的背景，减法并不完全。”这引入了噪声和不确定性，可能掩盖了细微的磁特征。
  • 宏观性质估算的敏感性: Weiss温度（关键宏观性质）的估算“对自旋-轨道耦合激发态的范德瓦尔斯顺磁类贡献非常敏感”。这使得它成为精确参数确定的不可靠指标，进一步复杂化了理论模型的验证。
  • 粉末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”（State-of-the-Art）方法，如标准CNN、基础扩散模型或Transformer，视为并随后拒绝用于此特定问题，是对科学领域的一种误解。该论文属于凝聚态物理学范畴，其中研究磁激发的“SOTA”方法是成熟的实验技术，如非弹性中子散射，以及理论框架，如线性自旋波理论。机器学习模型虽然在其他领域功能强大，但对于直接探测RuBr3等材料的根本磁相互作用和自旋动力学，在这一层面是不可适用的。

作者选择非弹性中子散射结合线性自旋波理论（LSWT）并非仅仅是偏好，而是由于问题的性质所决定的必然选择。核心目标是实验上观测并理论上解释RuBr3（一种Kitaev模型候选材料）中的磁激发。非弹性中子散射是直接测量磁性材料中自旋激发（磁振子或分数化激发）的能量和动量依赖性的首要实验技术。它为材料内部的动态磁关联提供了直接的窗口。

作者意识到传统（在凝聚态物理学背景下）“SOTA”方法不足以进行简单解释的时刻，很可能来自于对RuBr3的初步表征。论文的摘要和引言强调，尽管Kitaev模型预测了精确的量子自旋液体基态，但RuBr3表现出“$T_N$以下的锯齿形反铁磁有序”。这一观察立即表明，理想的、纯粹的Kitaev模型（将导致自旋液体）不足以描述RuBr3。磁有序的存在需要一个能够描述有序相中集体自旋激发（磁振子）的理论框架，而这正是线性自旋波理论应用于扩展Kitaev-Heisenberg哈密顿量时所提供的。因此，这种认识并非关于计算算法的失败，而是关于材料的物理性质偏离了最简单的理论理想，从而需要更全面的实验和理论方法。

比较优势

所选方法结合了非弹性中子散射实验和线性自旋波理论计算，通过提供微观磁相互作用与宏观可观测自旋动力学之间直接且物理上可解释的联系，展现出定性优势。

1. 直接探测自旋动力学: 非弹性中子散射在定性上更优越，因为它直接测量动态结构因子 $S(\mathbf{Q}, E)$，该因子与自旋-自旋关联函数的傅里叶变换成正比。这意味着它直接探测自旋在空间（动量$\mathbf{Q}$）和时间（能量$E$）上的涨落。与间接探测相比，中子散射提供了磁激发谱的完整图景，包括色散关系（激发能量如何随动量变化）和谱权重。这种直接性是其他研究块体磁性材料的技术所无法比拟的。
2. 通过LSWT实现物理可解释性: 线性自旋波理论应用于精心构建的哈密顿量（如方程2所示）时，可以提取基本的交换参数（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）。这提供了一个结构优势：它将观测到的实验谱与控制材料磁性的基本微观相互作用联系起来。本文通过比较两个不同参数集（J$_1$-K-$\Gamma$-$\Gamma$'和J$_1$-K-J$_2$-J$_3$模型）的模拟谱与实验数据（图3）明确展示了这一点。J$_1$-K-J$_2$-J$_3$模型显示出“更好”的一致性，表明其在捕捉RuBr3中相互作用的复杂相互作用方面具有更强的能力。这并非关于将内存复杂度从$O(N^2)$降低到$O(N)$，而是关于提供一个稳健的、物理上合理的模型，该模型能准确描述实验现实。能够重现实验中观测到的“强色散磁激发”和“无色散激发”是关键的定性优势。

与约束的对齐

所选方法与研究RuBr3等Kitaev模型候选材料（特别是表现出锯齿形反铁磁有序的材料）的隐含约束完美契合。

1. 磁激发实验表征: 主要约束是实验表征磁激发。非弹性中子散射因其直接测量自旋激发能量和动量的能力而特别适合此任务。使用粉末样品（而非单晶）是一个实际约束，分析中考虑了这一点，方法是在模拟过程中对整个实空间进行平均。
2. 识别潜在磁相互作用: 一个关键要求是识别和量化驱动材料磁行为并可能稳定观测到的锯齿形反铁磁有序的各种磁相互作用（Kitaev、Heisenberg、非对角）。应用于扩展Kitaev-Heisenberg哈密顿量的线性自旋波理论正是为此目的而设计的。它允许作者将实验数据拟合到理论模型，从而提取这些交换参数的值。这种实验观测与理论建模的“结合”对于理解配体取代如何影响这些相互作用至关重要，这是本文的核心主题。本文旨在“阐明阴离子取代对磁相互作用的影响”的目标直接通过这种组合方法得到满足。

替代方案的拒绝

本文隐含地拒绝了更简单的理论模型，并强调了其他表征技术的局限性，而不是直接否定不相关的“SOTA”机器学习方法。

1. 拒绝理想Kitaev模型: 最显著的隐含拒绝是拒绝将理想铁磁Kitaev模型作为RuBr3的完整描述。本文指出，“理论上，理想铁磁Kitaev模型的磁激发具有较弱的波矢依赖性”，但其实验结果显示出“强色散磁模式”和“锯齿形反铁磁有序”。这种明显的差异意味着纯Kitaev模型不足以描述。因此，作者采用了扩展哈密顿量（方程2），该哈密顿量包含额外的Heisenberg（$J_1, J_2, J_3$）和非对角（$\Gamma, \Gamma'$）项，这些项对于稳定观测到的磁有序和重现色散激发是必需的。这是对理论模型的改进，超越了简单的观点。
2. 其他表征方法的局限性: 虽然不是直接拒绝，但本文在引言中讨论$\alpha$-RuCl$_3$的先前工作时，提到了拉曼散射、核磁共振和X射线散射实验等其他技术。这些方法提供了有价值的补充信息（例如，特定的声子模式、局部自旋环境、电子结构），但它们没有提供非弹性中子散射所提供的关于自旋激发的直接、动量分辨信息。例如，拉曼散射可以检测双磁振子类激发，但中子散射提供完整的色散。
3. 机器学习SOTA的无关性: 如前所述，生成对抗网络（GANs）或扩散模型等方法旨在执行图像生成、数据合成或大型数据集中的复杂模式识别等任务。它们在直接实验测量自旋动力学或量子材料中导出微观交换参数方面存在根本性不适用性，这需要特定的物理探针和理论框架。它们在此背景下的应用将是类别错误，而不是失败的替代方案。

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_1+K$项放置在$x$或$y$方向，可以构建类似的矩阵$J_x$和$J_y$，这反映了Kitaev相互作用的键依赖性质。

项对项解剖

让我们逐一分解这些方程的每个元素以理解其作用：

来自方程（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$之间的最近邻各向异性相互作用矩阵。

  • 数学定义: 一个3x3矩阵（如方程3所示），其元素定义了自旋$S_i$和$S_j$的$x, y, z$分量之间耦合的强度和性质。
  • 物理/逻辑作用: 这一项是Kitaev模型的关键特征，捕捉了相邻Ru自旋之间复杂、键依赖和各向异性的交换相互作用。它对于理解奇异的自旋液体物理至关重要。矩阵形式是必需的，因为相互作用不是简单的点积，而是取决于特定的自旋分量和键的方向。

• $S_i \cdot S_j$: 这表示自旋$S_i$和$S_j$之间的各向同性点积相互作用。

  • 数学定义: 两个自旋向量的点积，$S_i^x S_j^x + S_i^y S_j^y + S_i^z S_j^z$。
  • 物理/逻辑作用: 这描述了传统的Heisenberg型相互作用，它们是各向同性的（方向无关）。它用于距离更远的自旋之间的相互作用，其中键依赖性通常不那么显著。

• $J_2$: 这是次近邻（NNN）各向同性磁相互作用的耦合常数。

  • 数学定义: 一个标量系数，缩放了NNN对的$S_i \cdot S_j$相互作用的强度。
  • 物理/逻辑作用: 它量化了相隔两个位点的自旋之间的各向同性磁耦合强度。该项在稳定特定的磁有序（如锯齿形反铁磁有序）中起作用。

• $\sum_{NNN}$: 这是所有次近邻（NNN）对的求和算符。

  • 数学定义: 对相隔两个位点的自旋对 $S_i, S_j$ 求和后续项。
  • 物理/逻辑作用: 确保所有NNN各向同性相互作用都包含在总能量计算中。

• $J_3$: 这是第三近邻（3NN）各向同性磁相互作用的耦合常数。

  • 数学定义: 一个标量系数，缩放了3NN对的$S_i \cdot S_j$相互作用的强度。
  • 物理/逻辑作用: 它量化了相隔三个位点的自旋之间的各向同性磁耦合强度。本文强调了它在稳定观测到的锯齿形反铁磁有序中的重要性。

• $\sum_{3NN}$: 这是所有第三近邻（3NN）对的求和算符。

  • 数学定义: 对相隔三个位点的自旋对 $S_i, S_j$ 求和后续项。
  • 物理/逻辑作用: 确保所有3NN各向同性相互作用都包含在内。

• $J_p$: 这是晶面间各向同性磁相互作用的耦合常数。

  • 数学定义: 一个标量系数，缩放了不同层中自旋对的$S_i \cdot S_j$相互作用的强度。
  • 物理/逻辑作用: 它量化了位于不同层中的自旋之间的各向同性磁耦合强度，考虑了材料的三维性质。

• $\sum_{interplane}$: 这是所有晶面间对的求和算符。

  • 数学定义: 对位于不同层中的自旋对 $S_i, S_j$ 求和后续项。
  • 物理/逻辑作用: 确保所有晶面间各向同性相互作用都包含在内。

来自方程（3）——最近邻相互作用矩阵 $J_z$：

• $J_1$: 这是各向同性Heisenberg交换耦合常数。

  • 数学定义: 出现在$J_{ij}$矩阵对角线上的标量系数，代表$S_x S_x$, $S_y S_y$, 和$S_z S_z$相互作用（不包括Kitaev项）的强度。
  • 物理/逻辑作用: 代表传统的、方向无关的最近邻交换相互作用部分。它是磁相互作用的基本组成部分。

• $K$: 这是Kitaev相互作用项。

  • 数学定义: 一个标量系数，添加到$J_{ij}$矩阵的对角线元素之一，特别是对于$z$键的$S_z S_z$分量。
  • 物理/逻辑作用: 这是Kitaev模型的定义特征，代表高度各向异性、键依赖的Ising类相互作用。它负责奇异的自旋液体物理。作者使用加法是因为Kitaev相互作用沿特定键方向作为各向同性Heisenberg项的附加各向异性分量。

• $\Gamma$: 这是对称的非对角交换项。

  • 数学定义: 出现在$J_{ij}$矩阵非对角位置的标量系数（例如，耦合$S_x S_y$和$S_y S_x$）。
  • 物理/逻辑作用: 代表一种各向异性交换相互作用，它在自旋分量交换下是对称的。它源于自旋-轨道耦合，可以与Kitaev相互作用竞争或增强，显著影响磁基态。

• $\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方法识别经典基态——最低能量自旋构型。哈密顿量被视为经典，系统搜索最小化总能量（方程2）同时保持每个$S_i$恒定自旋幅度的自旋排列。这产生了平衡自旋构型，例如锯齿形反铁磁有序，它作为后续量子计算的参考点。

4. 动量空间变换: 然后使用傅里叶变换将实空间哈密顿量转换为动量空间。这会将位点特定的自旋算符转换为动量依赖的算符，并将相互作用项转换为动量依赖的耦合函数（例如，方程4中的$J_{AAk}, J_{ABk}$）。这种变换通过利用晶格的平移对称性来简化问题。

5. Holstein-Primakoff近似: 为了描述量子激发，使用玻色子产生和湮灭算符来近似自旋算符。这是Holstein-Primakoff变换，它将小偏差（磁振子）从经典基态的自旋动力学线性化。这有效地将复杂的自旋哈密顿量转换为这些玻色子算符的二次形式。

6. Bogoliubov-de Gennes矩阵形成: 考虑到锯齿形反铁磁有序，系统由四个磁子晶格有效描述。二次哈密顿量在经过上述变换后，被转化为每个波矢$k$的Bogoliubov-de Gennes形式下的8x8矩阵。该矩阵在数学上描述了玻色子激发的耦合动力学。

7. 本征值提取（激发谱）: 对每个波矢$k$，对这个8x8的Bogoliubov-de Gennes矩阵进行对角化。得到的本征值直接对应于激发能量（磁振子色散关系），本征向量描述了这些磁振子模式的性质。此步骤生成非弹性中子散射的理论谱，显示能量如何随动量变化。

8. 动态结构因子计算: 从本征值和本征向量中，计算动态结构因子。该量与非弹性中子散射实验中测量的强度成正比，提供了理论与实验之间的直接联系。

9. 粉末平均与分辨率卷积: 最后，将计算出的谱在动量空间的所有可能方向上进行平均（以模拟粉末非弹性中子散射实验），然后与实验的波矢和能量分辨率进行卷积。这会产生一个理论谱，可以直接与图中所示的实验数据进行比较。

优化动力学

该机制中的“学习”或“收敛”主要是参数拟合和模型验证与实验数据对比的过程，而不是像机器学习中那样自主的迭代优化算法。其展开方式如下：

1. 初始参数猜测: 该过程从一组初始交换参数（$J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$）开始。这些值可能来自从头算计算、先前关于类似材料的研究或有根据的猜测。

2. 经典基态最小化: 对于每组参数，采用Luttinger-Tisza方法找到经典基态。这本身就是一个优化步骤：该方法迭代搜索最小化总能量（方程2）同时遵守固定自旋幅度的约束的自旋构型。这有效地在能量景观中导航以找到全局最小值。

3. 正向模型计算: 确定经典基态后，应用线性自旋波理论，通过确定性地计算磁振子激发谱和动态结构因子。这是模型的“前向传播”，基于当前参数生成理论预测。

4. 比较与差异评估: 然后将计算出的理论谱与实验非弹性中子散射数据（例如，图2中的强度图或图3和4中的积分强度）进行比较。“损失”或“差异”是理论预测与实验观测之间的差异。作者正在寻找在峰位、色散和整体谱权重等特征上定性和定量匹配。

5. 现象学参数调整: 然后作者现象学地调整交换参数。这并非典型的基于梯度的优化，而是一个迭代的、由人类指导的过程。他们可能会增加$J_3$以增强反铁磁有序，或调整$K$以影响自旋液体特性。例如，他们提到尝试了两种参数的“极端组合”（一种是主要的$\Gamma$，另一种是主要的$J_3$）来重现特定的特征，如色散和无色散激发。

6. 多约束验证: 除了非弹性中子散射谱之外，所选参数还通过其他宏观性质进行验证。例如，从参数计算Weiss温度（方程1）和磁矩倾斜角，并与实验值进行比较。这些作为参数必须满足的附加“损失函数”或约束，有助于缩小可行参数空间。作者注意到Weiss温度估算的敏感性，这表明这是拟合中一个微妙的部分。

7. 收敛到“最佳拟合”: 当找到一组参数能够跨所有可用实验数据提供最佳整体一致性时，该过程就收敛了。这个“最佳拟合”代表了作者对解释RuBr$_3$中观测现象的潜在磁相互作用的解读。因此，优化是一个搜索参数的过程，这些参数塑造理论磁振子色散和强度以匹配实验观测，通过迭代比较和调整有效地“学习”潜在的磁相互作用。

结果、局限性与结论

实验设计与基线

为了严格研究RuBr$_3$中的磁激发，研究人员采用了粉末非弹性中子散射这一探测自旋动力学的强大技术。研究的材料是RuBr$_3$的一种新多晶型物，它在R3空间群中结晶，形成三层蜂窝状结构，在结构上与$\alpha$-RuCl$_3$的低温相相同。研究RuBr$_3$的动机源于这样一个想法：取代配体离子（用溴取代氯）可以实验性地调节Kitaev、Heisenberg和非对角磁相互作用之间复杂的平衡。先前的研究已经表明，RuBr$_3$表现出与$\alpha$-RuCl$_3$不同的磁性质，例如较小的带隙、较低的电阻、抑制的磁化率、不同的Weiss温度以及较大的磁矩倾斜角。

实验使用三种不同的中子光谱仪进行：J-PARC的AMATERAS、ANSTO的PELICAN和JRR-3的GPTAS，确保了在各种能量和动量范围内的全面视图。多晶RuBr$_3$样品使用立方砧高压装置合成。将9.5克的样品制成圆柱体（直径15毫米，高15毫米），并用氦交换气密封在铝罐中。在不同温度下收集了非弹性中子散射谱：10 K、25 K、45 K、100 K和300 K。入射中子能量（$E_i$）设置为20.95、9.70、5.57和3.61 meV，并在42.17、15.19和7.75 meV下进行了额外的全能量分辨率测量。

数据分析的一个关键步骤是分离磁贡献。这是通过估算并减去声子贡献（主要来自300 K数据）以及使用温度依赖的Bose因子$1 + n(T) = (1 - e^{-E/(k_B T)})^{-1}$校正强度来实现的。Utsusemi软件套件用于此分析。

本研究中的“受害者”或基线模型主要是理论模型：线性自旋波理论（LSWT）模型，特别是J1-K-Γ-Γ'和J1-K-J2-J3模型，它们代表了Kitaev-Heisenberg哈密顿量中交换参数的不同组合。这些模型用于模拟预期的非弹性中子散射谱，从而可以直接与实验数据进行比较。此外，RuBr$_3$的实验结果与$\alpha$-RuCl$_3$和Na$_2$IrO$_3$等其他Kitaev候选材料的已知行为进行了隐含和明确的比较，作为经验基线以突出RuBr$_3$的独特特征。

证据证明的内容

RuBr$_3$的非弹性中子散射实验提供了确凿的证据，表明存在强反铁磁相互作用，这表现为以反铁磁区域中心为中心的强色散磁激发。在Néel温度（$T_N = 34$ K，通过中子衍射确认）以下，特别是在10 K和25 K时，观察到了明显的色散激发。这些激发在波矢为0.60 Å$^{-1}$和1.55 Å$^{-1}$处显示出峰值，这与磁反射一致，强烈表明其磁起源。在10 K时识别出1.5 meV的能量隙，该隙在25 K以上降至接近零。

随着温度升高，这些激发的强波矢依赖性一直持续到45 K。然而，在100 K（约$3T_N$）时发生了一个显著的变化，谱权重转移到零波矢（$\Gamma$点），然后在300 K时完全消失。这种转移到$\Gamma$点被解释为系统中存在铁磁相互作用的明确证据。相反，在布里渊区边界附近观察到的磁激发的鲁棒性表明存在强反铁磁相互作用，这对于稳定观测到的锯齿形反铁磁有序至关重要。

与RuBr$_3$与其他Kitaev候选材料相比，其磁激发的温度鲁棒性比$\alpha$-RuCl$_3$更接近Na$_2$IrO$_3$。尽管$\alpha$-RuCl$_3$的激发在$T_N$以上立即转移到$\Gamma$点，但RuBr$_3$的转移发生在更高的温度，约$3T_N$。这表明其相互作用的平衡与$\alpha$-RuCl$_3$不同。

实验谱随后与两个线性自旋波理论模型（J1-K-Γ-Γ'和J1-K-J2-J3）的模拟进行了严格比较。两种模型都能重现自旋激发的特征，包括色散和近乎无色散的模式。然而，J1-K-J2-J3模型（包含大的第三近邻各向同性磁相互作用$J_3$）与实验的波矢和能量转移依赖性（如图3(a)和3(b)所示）更一致。这强烈表明，大的最近邻对称非对角（$\Gamma$）或第三近邻各向同性磁相互作用（$J_3$）对于稳定锯齿形反铁磁有序和解释观测到的强色散至关重要。

尽管实验Weiss温度和磁矩倾斜角与模型预测存在一些差异，但本文最终得出结论，具有大$J_3$项的J1-K-J2-J3模型更为现实。这一结论主要基于其与非弹性中子散射谱的更好拟合，尽管从粉末平均数据中精确估计交换参数存在挑战，并且Weiss温度计算对各种贡献敏感。

最后，该研究提供了硬性证据，表明溴取代显著影响了磁相互作用。波矢分析表明，铁磁Kitaev项仍然存在，但反铁磁相互作用（特别是$\Gamma$或$J_3$）得到了增强。这种增强将RuBr$_3$的基态推向了更深的锯齿形反铁磁相，使其远离了理想的铁磁Kitaev自旋液体态。这种效应归因于更强的d-p杂化和通过配体原子增强的间接跳跃。

局限性与未来方向

尽管本文提出了令人信服的证据，但该研究承认了实验方法和理论建模固有的几项局限性。主要局限性源于粉末平均非弹性中子散射谱的使用。这种平均使得精确估计交换参数变得困难，正如两种不同理论模型（J1-K-Γ-Γ'和J1-K-J2-J3）重现的相似谱特征所证明的那样。用于建模的线性自旋波近似虽然有用，但可能不完全适用于赝自旋-1/2系统，因此仅提供磁相互作用的粗略估计。此外，实验得出的Weiss温度与模型预测之间存在显著差异，作者将其归因于估算的敏感性和潜在的范德瓦尔斯顺磁类贡献。同样，模型预测的倾斜角与实验值并不完全一致，尽管作者认为包含$J_3$的进一步各向异性可能会改善这一点。至关重要的是，本文指出“由于其波矢和能量依赖性在粉末非弹性散射谱中不清楚，因此无法断定弱激发是由Kitaev相互作用引起的”，这表明在明确识别Kitaev贡献方面存在重大差距。最后，Br取代导致反铁磁相互作用增强的确切起源“尚不清楚”，需要进一步研究潜在的电子机制。

展望未来，这些发现为未来的研究和发展开辟了几个激动人心的方向：

• 用于确定性证据的单晶研究: 最直接和影响深远的方向将是对高质量RuBr$_3$单晶进行非弹性中子散射实验。这将消除粉末平均的模糊性，从而能够更精确地确定交换参数，并更清晰地理解所有激发（包括难以捉摸的Kitaev贡献）的波矢和能量依赖性。这可以提供完全表征Kitaev相互作用所需的确定性、无可辩驳的证据。

• 先进的理论建模: 为了解决LSWT的局限性，未来的工作应侧重于开发和应用更先进的理论模型。像精确对角化、密度矩阵重整化群（DMRG）或量子蒙特卡洛模拟等技术，它们更适合赝自旋-1/2系统和强关联材料，可以更准确地描述自旋动力学、倾斜角和Weiss温度，从而弥合理论与实验之间的差距。

• 阐明配体取代机制: 本文强调了Br取代的显著影响，但指出增强反铁磁相互作用的确切机制尚不清楚。未来的研究可以结合先进的从头算计算和有针对性的实验探测（例如，共振非弹性X射线散射、X射线吸收光谱），以精确绘制配体离子变化如何影响d-p杂化、轨道交换途径，并最终影响特定的磁相互作用。这可能为设计具有所需性质的Kitaev材料提供一个预测框架。

• 探索RuX$_3$家族和多晶型物: 鉴于RuX$_3$（X = Br, I）形成多种多晶型物，系统地研究RuBr$_3$和RuI$_3$的其他结构相可能会揭示新的磁基态或Kitaev与其他相互作用的不同平衡。这种比较方法可以识别有利于难以捉摸的Kitaev自旋液体态的特定结构基序或化学环境。

• 外部扰动与相调控: 应用磁场、压力或应变等外部扰动可能是调控RuBr$_3$中竞争相互作用平衡的有力方法。在这些条件下研究自旋动力学可能会揭示相变，从而可能稳定自旋液体相或提供对竞争相互作用性质的更深入见解。

• 寻找分数化激发: 尽管当前研究强调类磁振子激发，但高能量（例如12和15 meV）处弱而宽的激发可能暗示着分数化激发连续体，这是Kitaev自旋液体的特征。未来的高分辨率实验，特别是在单晶上进行的实验，应该专门寻找Majorana费米子的迹象，即使基态是有序的。

• 晶面间耦合与三维效应: 哈密顿量包含晶面间相互作用（$J_p$），但其在观测动力学中的作用并未得到广泛讨论。考虑到三层蜂窝状结构，对晶面间耦合对磁激发和整体基态影响的更详细研究可能是有价值的，尤其是在理解偏离纯二维Kitaev物理学方面。

这些多样化的视角，从实验改进到先进的理论建模和材料设计，对于全面理解RuBr$_3$中相互作用的复杂相互作用以及推进对Kitaev自旋液体的探索至关重要。

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

与其他领域的同构性

结构骨架

本文提出了一个框架，用于模拟和分析结构化晶格上相互作用元素的集体动力学行为，其中相互作用是各向异性的且依赖于键。