Магнитные возбуждения кандидата в китаевские материалы RuBr3

Предыстория и академическая родословная

Происхождение и академическая родословная

Проблема, рассматриваемая в данной статье, берет свое начало в увлекательной области физики конденсированного состояния, в частности, в поиске и характеристике состояний квантового спинового льда (QSL). Этот поиск начался всерьез с теоретического предложения модели Китаева Алексеем Китаевым в 2006 году [1]. Эта модель, определяемая уникальными зависимыми от связи изинговскими взаимодействиями на сотовом (гексагональном) латтисе, примечательна тем, что ее можно точно решить, выявив основное состояние QSL. В этом экзотическом состоянии спины не упорядочиваются даже при абсолютном нуле температуры; вместо этого они остаются запутанными и флуктуирующими, порождая фрактализованные возбуждения, известные как майорановские фермионы [7]. Эти свойства обладают огромным потенциалом для отказоустойчивых квантовых вычислений.

Затем возникла задача найти реальные материалы, которые могли бы поддерживать эти неуловимые китаевские взаимодействия. Академическая область увидела значительный прорыв около 2009 года с осознанием того, что сильное спин-орбитальное взаимодействие в определенных соединениях переходных металлов, особенно с $d^5$ электронной конфигурацией, таких как ионы Ir$^{4+}$ и Ru$^{3+}$, может эффективно индуцировать эти зависимые от связи взаимодействия посредством механизма, включающего состояния псевдоспина-1/2 [8]. Это привело к появлению "китаевских материалов" как нового рубежа. Среди них $\alpha$-RuCl$_3$ быстро стал наиболее интенсивно изучаемым кандидатом [9-38].

Однако фундаментальным ограничением или "болевой точкой" предыдущих подходов и исследований $\alpha$-RuCl$_3$ является то, что, хотя он и демонстрирует некоторые признаки поведения квантового спинового льда Китаева при более высоких температурах, его основное состояние ниже определенной температуры перехода ($T_N$) представляет собой обычный зигзагообразный антиферромагнитный порядок. Это указывает на наличие других магнитных взаимодействий (таких как Гейзенберга и недиагональные члены), которые конкурируют с идеальным состоянием квантового спинового льда Китаева и в конечном итоге дестабилизируют его [31, 32, 39-43]. Эти конкурирующие взаимодействия препятствуют полной реализации желаемых свойств QSL. Авторы данной статьи мотивированы необходимостью понять и контролировать эти взаимодействия. Одной из ключевых стратегий является модификация лигандных ионов, окружающих магнитные атомы, поскольку эти ионы играют решающую роль в опосредовании орбитального обменного механизма, который порождает зависимые от связи взаимодействия. Данная статья исследует RuBr$_3$ как новый кандидат, где замена хлора на бром, как ожидается, изменит баланс этих магнитных взаимодействий, потенциально приближая или отдаляя материал от идеального состояния квантового спинового льда Китаева.

Интуитивные термины предметной области

1. Модель Китаева: Представьте себе особый вид магнитной шашечной доски, где правила взаимодействия двух соседних шашек (притяжение или отталкивание) полностью зависят от направления линии, соединяющей их. Например, шашки, соединенные горизонтально, могут стремиться выстроиться в одном направлении, в то время как те, что соединены по диагонали, могут стремиться выстроиться в противоположных направлениях. Модель Китаева — это теоретическая основа, описывающая такую систему, и она знаменита тем, что предсказывает состояние "квантового спинового льда", в котором магниты никогда не замерзают в фиксированном узоре, даже при самых низких температурах.
2. Квантовый спиновый лед (QSL): Подумайте об обычной жидкости, где молекулы воды постоянно движутся и взаимодействуют. Квантовый спиновый лед похож, но вместо воды это крошечные внутренние магниты (спины) электронов находятся в постоянном, запутанном, "жидкостном" состоянии. Даже при абсолютном нуле они не образуют простой, упорядоченной магнитной структуры, как магнит на холодильнике. Вместо этого их коллективное поведение является высокосложным и квантово-механическим, что приводит к экзотическим "получастицам", называемым майорановскими фермионами.
3. Майорановские фермионы: Если представить обычную частицу, такую как электрон, как целого человека, то майорановский фермион будет подобен "половине человека". В контексте квантовых спиновых льдов спин электрона может эффективно расщепляться на две отдельные, независимые части, каждая из которых ведет себя как майорановский фермион. Эти "получастицы" уникальны тем, что они являются своими собственными античастицами, а их устойчивость к локальным возмущениям делает их перспективным строительным блоком для будущих квантовых компьютеров.
4. Зависимые от связи изинговские взаимодействия: Представьте сетку крошечных стрелок компаса, каждая из которых может указывать строго вверх или вниз (как бинарный переключатель). "Зависимые от связи" означает, что предпочтение двух соседних стрелок указывать в одном или противоположных направлениях полностью зависит от конкретного направления линии, соединяющей их на сетке. Например, вдоль связи "север-юг" они могут предпочитать быть параллельными, а вдоль связи "восток-запад" — антипараллельными. Это направленное предпочтение — основная идея модели Китаева.
5. Состояние псевдоспина-1/2: В некоторых сложных материалах фактический спин электрона (который всегда равен 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$, материале-кандидате на роль квантового спинового льда Китаева, и понимании того, как эти взаимодействия отклоняют его основное состояние от идеального квантового спинового льда Китаева.

Входное/текущее состояние: Материал RuBr$_3$, представляющий собой новый полиморф со слоистой сотовой структурой, изоструктурный хорошо изученному кандидату в китаевские материалы $\alpha$-RuCl$_3$. Оба материала демонстрируют зигзагообразный антиферромагнитный (AFM) порядок ниже своих соответствующих температур Нееля ($T_N$). Однако RuBr$_3$ проявляет отличные от $\alpha$-RuCl$_3$ магнитные свойства, такие как подавленная магнитная восприимчивость и другая температура Вайса. Предыдущие ab initio расчеты для RuBr$_3$ предполагали магнитные взаимодействия, сравнимые с $\alpha$-RuCl$_3$, что противоречит этим экспериментальным наблюдениям. Теоретическое понимание идеальной ферромагнитной модели Китаева предсказывает спиновые возбуждения со слабой зависимостью от волнового вектора, но реальные материалы часто демонстрируют более сложное поведение из-за конкурирующих взаимодействий.

Желаемый конечный результат (выходное/целевое состояние): Четкое понимание влияния замещения анионов (в частности, брома на хлор) на магнитные взаимодействия в кандидатах в китаевские материалы. Статья направлена на экспериментальное исследование спиновой динамики RuBr$_3$ с использованием неупругого рассеяния нейтронов (INS) для идентификации конкретных магнитных взаимодействий ($J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$), которые стабилизируют наблюдаемый зигзагообразный AFM порядок. В конечном итоге цель состоит в том, чтобы предоставить более реалистичную модель для RuBr$_3$, которая точно воспроизводит его магнитные возбуждения и макроскопические свойства, тем самым объясняя, как его основное состояние отклоняется от идеального ферромагнитного квантового спинового льда Китаева.

Отсутствующее звено или математический пробел: Количественное определение точного набора обменных параметров для RuBr$_3$, который мог бы последовательно объяснить как его макроскопические магнитные свойства (например, температуру Вайса и угол наклона), так и его микроскопическую спиновую динамику (спектры магнитных возбуждений). Хотя существуют теоретические модели (например, модели J$_1$-K-$\Gamma$-$\Gamma$' и J$_1$-K-J$_2$-J$_3$), статья подчеркивает трудность точной оценки этих параметров по усредненным по порошку экспериментальным данным. Существует явное расхождение между теоретическими предсказаниями ab initio расчетов и наблюдаемым экспериментальным магнитным поведением RuBr$_3$, которое данное исследование пытается преодолеть посредством детальной экспериментальной характеристики спиновой динамики.

Болезненный компромисс или дилемма: Которая поставила в тупик предыдущих исследователей и является центральной для данной работы, — это присущая конкуренция между желаемыми китаевскими взаимодействиями и другими магнитными взаимодействиями (Гейзенберга, недиагональными, дальних соседей) в реальных материалах. Хотя модель Китаева теоретически привлекательна для реализации спиновых льдов, наличие этих дополнительных взаимодействий часто приводит к обычным магнитным порядкам, таким как зигзагообразный AFM, наблюдаемый в RuBr$_3$. Исследователи сталкиваются с дилеммой, как регулировать эти конкурирующие взаимодействия, чтобы либо подавить обычный порядок, либо усилить китаевскую физику. Замещение лиганда является стратегией, но его точное влияние на тонкий баланс взаимодействий сложно и нелегко предсказуемо, часто приводя к непреднамеренным последствиям или новым проблемам. В статье прямо указано, что "контроль не только китаевских взаимодействий, но и гейзенберговских и недиагональных магнитных взаимодействий имеет важное значение для практического использования реальных материалов", подчеркивая этот фундаментальный компромисс. Кроме того, трудность однозначного определения нескольких обменных параметров по усредненным по порошку экспериментальным данным означает, что улучшение экспериментальной доступности (использование порошковых образцов) достигается за счет точности извлечения параметров, поскольку разные наборы параметров могут давать схожие смоделированные спектры.

Ограничения и режимы отказа

Проблема понимания магнитных возбуждений в RuBr$_3$ чрезвычайно усложняется несколькими суровыми, реалистичными ограничениями:

• Физические ограничения:

  • Сложность материала: RuBr$_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$). В статье прямо указано, что "невозможно одновременно оценить несколько обменных параметров" из имеющихся данных. Это вынуждает исследователей полагаться на феноменологические корректировки и представлять "две крайние комбинации", а не уникальный, окончательный набор параметров.
  • Усреднение по порошку: Использование порошковых данных неупругого рассеяния нейтронов, хотя и более доступно экспериментально, по своей сути усредняет все кристаллографические ориентации. Эта "трудность в оценке обменных параметров по усредненному по порошку спектру" означает потерю критически важной анизотропной информации, которая могла бы помочь различить различные модели.

• Ограничения, основанные на данных:

  • Неполное вычитание фононов: Для выделения магнитного сигнала фононные вклады оцениваются по данным при высоких температурах и вычитаются. Однако "вычитание не является полным из-за ангармонизма фононов и фона, который мало зависит от температуры". Это вносит шум и неопределенность, потенциально скрывая тонкие магнитные особенности.
  • Чувствительность оценки макроскопических свойств: Оценка температуры Вайса, ключевого макроскопического свойства, "очень чувствительна к вкладам парамагнитного типа Ван Флека от возбужденных состояний со спин-орбитальным взаимодействием". Это делает ее ненадежным показателем для точного определения параметров, что еще больше усложняет проверку теоретических моделей.
  • Ограниченное разрешение порошкового 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, базовые диффузионные модели или трансформеры, которые рассматриваются, а затем отвергаются для данной конкретной проблемы, является несколько искаженным представлением научной области. Данная статья относится к области физики конденсированного состояния, где "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$). Это обеспечивает структурное преимущество: оно связывает наблюдаемые экспериментальные спектры с лежащими в основе микроскопическими взаимодействиями, управляющими магнетизмом материала. Статья явно демонстрирует это, сравнивая смоделированные спектры из двух различных наборов параметров (модели 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, особенно того, который демонстрирует зигзагообразный антиферромагнитный порядок.

1. Экспериментальная характеристика магнитных возбуждений: Основное ограничение заключается в экспериментальной характеристике магнитных возбуждений. Неупругое рассеяние нейтронов уникально подходит для этой цели, поскольку оно напрямую измеряет энергию и импульс спиновых возбуждений. Использование порошковых образцов (в отличие от монокристаллов) является практическим ограничением, и анализ учитывает это путем усреднения по всему телесному углу при моделировании.
2. Идентификация лежащих в основе магнитных взаимодействий: Важным требованием является идентификация и количественная оценка различных магнитных взаимодействий (китаевских, гейзенберговских, недиагональных), которые обуславливают магнитное поведение материала и потенциально стабилизируют наблюдаемый зигзагообразный антиферромагнитный порядок. Линейная спин-волновая теория, примененная к расширенному гамильтониану Китаева-Гейзенберга, предназначена именно для этой цели. Она позволяет авторам подгонять экспериментальные данные к теоретической модели, тем самым извлекая значения этих обменных параметров. Этот "брак" между экспериментальным наблюдением и теоретическим моделированием необходим для понимания того, как замещение лиганда влияет на эти взаимодействия, что является центральной темой статьи. Цель статьи "прояснить влияние замещения анионов на магнитные взаимодействия" напрямую достигается этим комбинированным подходом.

Отклонение альтернатив

Статья неявно отвергает более простые теоретические модели и подчеркивает ограничения других методов характеристики, вместо того чтобы явно отвергать несвязанные "SOTA" подходы машинного обучения.

1. Отвержение идеальной модели Китаева: Наиболее значимым неявным отвержением является отвержение идеальной ферромагнитной модели Китаева как полного описания RuBr3. В статье говорится, что "Теоретически, магнитные возбуждения идеальной ферромагнитной модели Китаева характеризуются слабой зависимостью от волнового вектора", но их экспериментальные результаты показывают "сильно диспергирующие магнитные моды" и "зигзагообразный антиферромагнитный порядок". Это явное расхождение означает, что чистая модель Китаева недостаточна. Поэтому авторы принимают расширенный гамильтониан (Уравнение 2), который включает дополнительные члены Гейзенберга ($J_1, J_2, J_3$) и недиагональные ($\Gamma, \Gamma'$) члены, необходимые для стабилизации наблюдаемого магнитного порядка и воспроизведения диспергирующих возбуждений. Это уточнение теоретической модели, выходящее за рамки упрощенного взгляда.
2. Ограничения других методов характеристики: Хотя это и не прямое отвержение, в статье упоминаются другие методы, такие как рамановская спектроскопия, ядерный магнитный резонанс и рентгеновские эксперименты, во введении при обсуждении предыдущих работ по $\alpha$-RuCl$_3$. Эти методы предоставляют ценную дополнительную информацию (например, специфические фононные моды, локальные спиновые среды, электронную структуру), но они не предоставляют прямого, разрешенного по импульсу информации о спиновых возбуждениях, которую дает неупругое рассеяние нейтронов. Например, рамановская спектроскопия может обнаруживать двухмагноноподобные возбуждения, но рассеяние нейтронов дает полную дисперсию.
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$.

  • Математическое определение: Векторный оператор, представляющий квантово-механический спиновый угловой момент электрона на определенном атомном узле (в данном случае, псевдоспин-1/2 для Ru$^{3+}$). Он имеет компоненты $S_i^x, S_i^y, S_i^z$.
  • Физическая/логическая роль: Это фундаментальные магнитные степени свободы. Взаимодействия между этими спинами порождают наблюдаемые магнитные возбуждения.

• $J_{ij}$: Это матрица анизотропного взаимодействия ближайших соседей между спинами $S_i$ и $S_j$.

  • Математическое определение: Матрица 3x3 (как показано в Уравнении 3), элементы которой определяют силу и характер сопряжения между $x, y, z$ компонентами спинов $S_i$ и $S_j$.
  • Физическая/логическая роль: Этот член является ядром модели Китаева, охватывая сложные, зависимые от связи и анизотропные обменные взаимодействия между соседними спинами 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$.
  • Физическая/логическая роль: Это описывает обычные взаимодействия типа Гейзенберга, которые изотропны (независимы от направления). Он используется для взаимодействий между спинами, находящимися на большем расстоянии, где зависимость от связи обычно менее значима.

• $J_2$: Это константа связи для изотропных магнитных взаимодействий следующего ближайшего соседа (NNN).

  • Математическое определение: Скалярный коэффициент, масштабирующий силу взаимодействия $S_i \cdot S_j$ для пар NNN.
  • Физическая/логическая роль: Количественно определяет силу изотропного магнитного сопряжения между спинами, разделенными двумя узлами. Этот член играет роль в стабилизации определенных магнитных порядков, таких как зигзагообразный антиферромагнитный порядок.

• $\sum_{NNN}$: Это оператор суммирования по всем парам следующего ближайшего соседа.

  • Математическое определение: Суммирует последующий член для каждой пары спинов $S_i, S_j$, находящихся на расстоянии двух узлов.
  • Физическая/логическая роль: Гарантирует включение всех изотропных взаимодействий NNN в расчет полной энергии.

• $J_3$: Это константа связи для изотропных магнитных взаимодействий третьего ближайшего соседа (3NN).

  • Математическое определение: Скалярный коэффициент, масштабирующий силу взаимодействия $S_i \cdot S_j$ для пар 3NN.
  • Физическая/логическая роль: Количественно определяет силу изотропного магнитного сопряжения между спинами, разделенными тремя узлами. В статье подчеркивается его важность для стабилизации наблюдаемого зигзагообразного антиферромагнитного порядка.

• $\sum_{3NN}$: Это оператор суммирования по всем парам третьего ближайшего соседа.

  • Математическое определение: Суммирует последующий член для каждой пары спинов $S_i, S_j$, находящихся на расстоянии трех узлов.
  • Физическая/логическая роль: Гарантирует включение всех изотропных взаимодействий 3NN.

• $J_p$: Это константа связи для межплоскостных изотропных магнитных взаимодействий.

  • Математическое определение: Скалярный коэффициент, масштабирующий силу взаимодействия $S_i \cdot S_j$ для спинов в разных слоях.
  • Физическая/логическая роль: Количественно определяет силу изотропного магнитного сопряжения между спинами, расположенными на разных сотовых слоях, учитывая трехмерную природу материала.

• $\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}$, в частности, к компоненте $S_z S_z$ для $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}$ (например, $J_z$ в Уравнении 3 или ее аналоги $J_x, J_y$) в зависимости от направления связи, включая значения $J_1, K, \Gamma, \Gamma'$. Для связей дальних соседей и межплоскостных связей применяются скалярные члены $J_2, J_3, J_p$. Этот шаг создает всеобъемлющую энергетическую функцию, описывающую все магнитные взаимодействия.

3. Определение классического основного состояния: Прежде чем изучать квантовые флуктуации, определяется классическое основное состояние — конфигурация спинов с наименьшей энергией — с использованием метода Люттингера-Тизы. Гамильтониан рассматривается классически, и система ищет такое расположение спинов, которое минимизирует полную энергию (Уравнение 2) при сохранении постоянной величины спина для каждого $S_i$. Это дает равновесную конфигурацию спинов, такую как зигзагообразный антиферромагнитный порядок, которая служит точкой отсчета для последующих квантовых расчетов.

4. Преобразование в импульсное пространство: Гамильтониан в реальном пространстве затем преобразуется в импульсное пространство с помощью Фурье-преобразования. Это преобразует операторы спина, зависящие от узла, в операторы, зависящие от импульса, а члены взаимодействия становятся функциями связи, зависящими от импульса (например, $J_{AAk}, J_{ABk}$ в Уравнении 4). Это преобразование упрощает задачу, используя трансляционную симметрию решетки.

5. Приближение Холстейна-Примакова: Для описания квантовых возбуждений спиновые операторы аппроксимируются бозонными операторами рождения и уничтожения. Это преобразование Холстейна-Примакова, которое линеаризует спиновую динамику для малых отклонений (магнонов) от классического основного состояния. Это эффективно преобразует сложный спиновый гамильтониан в квадратичную форму относительно этих бозонных операторов.

6. Формирование матрицы Боголюбова-де Женера: Учитывая зигзагообразный антиферромагнитный порядок, система эффективно описывается четырьмя магнитными подрешетками. Квадратичный гамильтониан, после предыдущих преобразований, затем приводится к матрице 8x8 в форме Боголюбова-де Женера для каждого волнового вектора $k$. Эта матрица математически описывает связанные динамики бозонных возбуждений.

7. Извлечение собственных значений (спектр возбуждений): Эта матрица 8x8 Боголюбова-де Женера диагонализируется для каждого волнового вектора $k$. Полученные собственные значения напрямую соответствуют энергиям возбуждений (дисперсионные соотношения магнонов), а собственные векторы описывают природу этих магнонных мод. Этот шаг генерирует теоретический спектр неупругого рассеяния нейтронов, показывающий, как энергия зависит от импульса.

8. Расчет динамического фактора структуры: Из собственных значений и собственных векторов рассчитывается динамический фактор структуры. Эта величина прямо пропорциональна интенсивности, измеряемой в экспериментах по неупругому рассеянию нейтронов, обеспечивая прямую связь между теорией и экспериментом.

9. Усреднение по порошку и свертка с разрешением: Наконец, рассчитанный спектр усредняется по всем возможным направлениям в импульсном пространстве (для моделирования порошковых экспериментов по неупругому рассеянию нейтронов), а затем свертывается с экспериментальным разрешением по волновому вектору и энергии. Это дает теоретический спектр, который можно напрямую сравнить с экспериментальными данными, показанными на рисунках.

Динамика оптимизации

"Обучение" или "сходимость" в этом механизме — это в первую очередь процесс подгонки параметров и проверки модели по экспериментальным данным, а не автономный итеративный алгоритм оптимизации, как в машинном обучении. Вот как это происходит:

1. Начальная оценка параметров: Процесс начинается с начального набора обменных параметров ($J_1, K, \Gamma, \Gamma', J_2, J_3, J_p$). Эти значения могут быть получены из ab initio расчетов, предыдущих исследований аналогичных материалов или обоснованных предположений.

2. Минимизация классического основного состояния: Для каждого набора параметров используется метод Люттингера-Тизы для нахождения классического основного состояния. Это само по себе является этапом оптимизации: метод итеративно ищет конфигурацию спинов, которая минимизирует полную энергию (Уравнение 2), соблюдая при этом ограничение на фиксированные величины спинов. Это эффективно навигирует по энергетическому ландшафту для поиска глобального минимума.

3. Расчет прямой модели: После установления классического основного состояния применяется линейная спин-волновая теория для детерминированного расчета спектра магнонных возбуждений и динамического фактора структуры. Это "прямой проход" модели, генерирующий теоретический прогноз на основе текущих параметров.

4. Сравнение и оценка расхождений: Рассчитанный теоретический спектр затем сравнивается с экспериментальными данными неупругого рассеяния нейтронов (например, карты интенсивности на Рис. 2 или интегрированные интенсивности на Рис. 3 и 4). "Потеря" или "расхождение" — это разница между теоретическим прогнозом и экспериментальным наблюдением. Авторы ищут качественное и количественное совпадение таких признаков, как положения пиков, дисперсия и общий спектральный вес.

5. Феноменологическая корректировка параметров: Затем авторы феноменологически корректируют обменные параметры. Это не градиентная оптимизация в обычном смысле, а скорее итеративный, управляемый человеком процесс. Они могут увеличить $J_3$ для усиления антиферромагнитного порядка или скорректировать $K$ для влияния на характеристики спинового льда. Например, они упоминают попытку двух "крайних комбинаций" параметров (одна с преобладающим $\Gamma$, другая с преобладающим $J_3$) для воспроизведения специфических признаков, таких как диспергирующие и дисперсионно-независимые возбуждения.

6. Проверка по нескольким ограничениям: Помимо спектра неупругого рассеяния нейтронов, выбранные параметры также проверяются по другим макроскопическим свойствам. Например, температура Вайса (Уравнение 1) и угол наклона магнитного момента рассчитываются по параметрам и сравниваются с экспериментальными значениями. Они служат дополнительными "функциями потерь" или ограничениями, которым должны удовлетворять параметры, помогая сузить допустимое пространство параметров. Авторы отмечают чувствительность оценки температуры Вайса, указывая, что это деликатная часть подгонки.

7. Сходимость к "наилучшей подгонке": Процесс сходится, когда найден набор параметров, обеспечивающий наилучшее общее соответствие всем доступным экспериментальным данным. Эта "наилучшая подгонка" представляет собой интерпретацию авторами лежащих в основе магнитных взаимодействий, которые объясняют наблюдаемые явления в RuBr$_3$. Таким образом, оптимизация представляет собой поиск параметров, которые формируют теоретическую дисперсию магнонов и интенсивность, чтобы соответствовать экспериментальным наблюдениям, эффективно "изучая" лежащие в основе магнитные взаимодействия посредством итеративного сравнения и корректировки.

Результаты, ограничения и заключение

Экспериментальный дизайн и базовые уровни

Для тщательного исследования магнитных возбуждений в RuBr3 исследователи использовали порошковое неупругое рассеяние нейтронов — мощный метод зондирования спиновой динамики. Исследуемым материалом был новый полиморф RuBr3, который кристаллизуется в пространственной группе R3, образуя трехслойную сотовую структуру, структурно идентичную низкотемпературной фазе $\alpha$-RuCl3. Мотивация изучения RuBr3 проистекает из идеи, что замещение лигандных ионов (брома на хлор) может экспериментально регулировать сложный баланс между китаевскими, гейзенберговскими и недиагональными магнитными взаимодействиями. Предыдущие исследования уже показали, что RuBr3 проявляет отличные от $\alpha$-RuCl3 магнитные свойства, такие как меньшая запрещенная зона, сниженное сопротивление, подавленная магнитная восприимчивость, другая температура Вайса и больший угол наклона магнитного момента.

Эксперименты проводились с использованием трех различных спектрометров нейтронов: AMATERAS в J-PARC, PELICAN в ANSTO и GPTAS в JRR-3, что обеспечило всесторонний обзор в различных диапазонах энергии и импульса. Поликристаллический образец RuBr3 был синтезирован с использованием аппарата высокого давления с кубическим наковальней. Образец массой 9,5 г был сформирован в цилиндр (диаметр 15 мм, высота 15 мм) и запаян в алюминиевую банку с гелиевым теплообменным газом. Спектры неупругого рассеяния нейтронов были собраны при различных температурах: 10 К, 25 К, 45 К, 100 К и 300 К. Энергии падающих нейтронов ($E_i$) были установлены на 20,95, 9,70, 5,57 и 3,61 МэВ, с дополнительными измерениями с высоким энергетическим разрешением при 42,17, 15,19 и 7,75 МэВ.

Важным этапом анализа данных было выделение магнитных вкладов. Это было достигнуто путем оценки и вычитания фононных вкладов, в основном полученных из данных при 300 К, и коррекции интенсивностей с использованием температурно-зависимого бозевского фактора $1 + n(T) = (1 - e^{-E/(k_B T)})^{-1}$. Для этого анализа использовался пакет программ Utsusemi.

"Жертвами" или базовыми моделями в данном исследовании были в основном теоретические: модели линейной спин-волновой теории (LSWT), в частности модели J1-K-Γ-Γ' и J1-K-J2-J3, которые представляют собой различные комбинации обменных параметров в гамильтониане Китаева-Гейзенберга. Эти модели использовались для моделирования ожидаемых спектров неупругого рассеяния нейтронов, что позволило провести прямое сравнение с экспериментальными данными. Кроме того, экспериментальные результаты для RuBr3 косвенно и явно сравнивались с известным поведением других кандидатов в китаевские материалы, таких как $\alpha$-RuCl3 и Na2IrO3, служа эмпирическими базовыми уровнями для выделения уникальных характеристик RuBr3.

Что доказывают свидетельства

Эксперименты по неупругому рассеянию нейтронов на RuBr3 предоставили убедительные доказательства наличия сильных антиферромагнитных взаимодействий, проявляющихся в виде сильно диспергирующих магнитных возбуждений, центрированных в антиферромагнитном центре зоны Бриллюэна. Ниже температуры Нееля ($T_N = 34$ К, подтверждено дифракцией нейтронов), в частности при 10 К и 25 К, наблюдались выраженные диспергирующие возбуждения. Эти возбуждения показывали пики при волновых векторах 0,60 Å$^{-1}$ и 1,55 Å$^{-1}$, что соответствует магнитным отражениям, убедительно свидетельствуя об их магнитном происхождении. При 10 К была обнаружена энергетическая щель в 1,5 МэВ, которая затем уменьшилась почти до нуля выше 25 К.

С повышением температуры сильная зависимость этих возбуждений от волнового вектора сохранялась до 45 К. Однако при 100 К (примерно $3T_N$) произошло значительное смещение, когда спектральная плотность сместилась к нулевому волновому вектору (точка $\Gamma$), прежде чем полностью исчезнуть при 300 К. Это смещение к точке $\Gamma$ интерпретируется как явное свидетельство наличия ферромагнитных взаимодействий в системе. И наоборот, устойчивость магнитных возбуждений, наблюдаемая вблизи границы зоны Бриллюэна, указывает на сильные антиферромагнитные взаимодействия, которые имеют решающее значение для стабилизации наблюдаемого зигзагообразного антиферромагнитного порядка.

По сравнению с другими кандидатами в китаевские материалы, температурная устойчивость магнитных возбуждений RuBr3 оказалась ближе к таковой у Na2IrO3, чем у $\alpha$-RuCl3. В то время как возбуждения $\alpha$-RuCl3 смещаются к точке $\Gamma$ сразу выше $T_N$, смещение RuBr3 происходит при более высокой температуре, около $3T_N$. Это предполагает иной баланс взаимодействий по сравнению с $\alpha$-RuCl3.

Экспериментальные спектры затем были строго сопоставлены с симуляциями из двух моделей линейной спин-волновой теории: J1-K-Γ-Γ' и J1-K-J2-J3. Обе модели смогли воспроизвести характерные особенности спиновых возбуждений, включая как диспергирующие, так и почти дисперсионно-независимые моды. Однако модель J1-K-J2-J3, которая включает большое третье соседнее изотропное магнитное взаимодействие ($J_3$), показала лучшее согласие с экспериментальной зависимостью от волнового вектора и энергии (как показано на Рис. 3(a) и 3(b)). Это убедительно свидетельствует о том, что сильные симметричные недиагональные ($\Gamma$) взаимодействия ближайших соседей или третьего ближайшего соседа ($J_3$) необходимы для стабилизации зигзагообразного антиферромагнитного порядка и объяснения наблюдаемой сильной дисперсии.

Хотя экспериментальные температуры Вайса и углы наклона магнитного момента показали некоторые расхождения с предсказаниями моделей, статья в конечном итоге приходит к выводу, что модель J1-K-J2-J3 с большим членом $J_3$ является более реалистичной. Этот вывод основан в первую очередь на лучшей подгонке к спектрам неупругого рассеяния нейтронов, несмотря на трудности с точной оценкой обменных параметров по усредненным по порошку данным и чувствительность расчетов температуры Вайса к различным вкладам.

Наконец, исследование предоставляет неопровержимые доказательства того, что замещение брома значительно влияет на магнитные взаимодействия. Анализ волновых векторов предполагает, что ферромагнитные китаевские члены сохраняются, но антиферромагнитные взаимодействия (в частности, $\Gamma$ или $J_3$) усиливаются. Это усиление смещает основное состояние RuBr3 глубже в зигзагообразную антиферромагнитную фазу, отдаляя его от идеального ферромагнитного квантового спинового льда Китаева. Этот эффект приписывается более сильной d-p гибридизации и усиленному косвенному переносу через лигандные атомы.

Ограничения и будущие направления

Несмотря на убедительные представленные доказательства, исследование признает несколько ограничений, присущих экспериментальному подходу и теоретическому моделированию. Основное ограничение связано с использованием усредненных по порошку спектров неупругого рассеяния нейтронов. Это усреднение затрудняет точную оценку обменных параметров, как показывают схожие спектральные особенности, воспроизведенные двумя различными теоретическими моделями (J1-K-Γ-Γ' и J1-K-J2-J3). Применяемые для моделирования приближения линейной спин-волновой теории, хотя и полезны, могут не идеально подходить для систем псевдоспина-1/2 и поэтому дают лишь грубую оценку магнитных взаимодействий. Кроме того, наблюдалось заметное расхождение между экспериментально полученными температурами Вайса и теми, что предсказывались моделями, что авторы объясняют чувствительностью оценки и возможными вкладами парамагнитного типа Ван Флека. Аналогично, предсказанные моделью углы наклона не полностью совпадали с экспериментальными значениями, хотя авторы предполагают, что включение дальнейшей анизотропии в $J_3$ может улучшить это. Важно отметить, что статья заявляет, что "не удалось заключить, что слабые возбуждения индуцированы китаевскими взаимодействиями, поскольку их зависимости от волнового вектора и энергии не ясны в спектре порошкового неупругого рассеяния", что указывает на значительный пробел в однозначной идентификации китаевских вкладов. Наконец, точное происхождение усиленных антиферромагнитных взаимодействий из-за замещения Br остается "неясным", требуя дальнейшего изучения лежащих в основе электронных механизмов.

Заглядывая вперед, эти результаты открывают несколько захватывающих направлений для будущих исследований и разработок:

• Исследования монокристаллов для окончательных доказательств: Наиболее непосредственным и значимым будущим направлением будет проведение экспериментов по неупругому рассеянию нейтронов на высококачественных монокристаллах RuBr3. Это устранит неоднозначность усреднения по порошку, позволяя гораздо более точное определение обменных параметров и более четкое понимание зависимостей от волнового вектора и энергии всех возбуждений, включая неуловимые китаевские вклады. Это может предоставить окончательные, неоспоримые доказательства, необходимые для полной характеристики китаевских взаимодействий.

• Продвинутое теоретическое моделирование: Для устранения ограничений LSWT будущая работа должна быть сосредоточена на разработке и применении более сложных теоретических моделей. Такие методы, как точная диагонализация, метод ренормализационной группы матрицы плотности (DMRG) или квантовые Монте-Карло симуляции, которые лучше подходят для систем псевдоспина-1/2 и сильно коррелированных материалов, могут обеспечить более точное описание спиновой динамики, углов наклона и температур Вайса, сокращая разрыв между теорией и экспериментом.

• Прояснение механизмов замещения лиганда: В статье подчеркивается значительное влияние замещения Br, но отмечается, что точный механизм усиления антиферромагнитных взаимодействий неясен. Будущие исследования могли бы сочетать продвинутые ab initio расчеты с целенаправленными экспериментальными зондами (например, резонансным неупругим рассеянием рентгеновских лучей, рентгеновской абсорбционной спектроскопией) для точного картирования того, как изменения лигандных ионов влияют на d-p гибридизацию, орбитальные обменные пути и, в конечном итоге, на специфические магнитные взаимодействия. Это могло бы привести к созданию предсказательной основы для разработки китаевских материалов с желаемыми свойствами.

• Исследование семейства RuX3 и полиморфов: Учитывая, что RuX3 (X = Br, I) образует несколько полиморфов, систематическое исследование других структурных фаз RuBr3 и RuI3 может выявить новые магнитные основные состояния или иные балансы китаевских и других взаимодействий. Этот сравнительный подход может выявить специфические структурные мотивы или химические среды, которые способствуют неуловимому состоянию квантового спинового льда Китаева.

• Внешние возмущения и настройка фазы: Применение внешних возмущений, таких как магнитные поля, давление или деформация, может быть мощным способом регулирования баланса конкурирующих взаимодействий в RuBr3. Исследование спиновой динамики в этих условиях может выявить фазовые переходы, потенциально приводя к стабилизации фазы спинового льда или предоставляя более глубокое понимание природы конкурирующих взаимодействий.

• Поиск фрактализованных возбуждений: Хотя текущее исследование подчеркивает магноноподобные возбуждения, наличие слабых, широких возбуждений при более высоких энергиях (например, 12 и 15 МэВ) может указывать на континуум фрактализованных возбуждений, характерный для квантового спинового льда Китаева. Будущие эксперименты с высоким разрешением, особенно на монокристаллах, должны специально искать признаки майорановских фермионов, даже если основное состояние упорядочено.

• Межплоскостное сопряжение и 3D-эффекты: Гамильтониан включает межплоскостные взаимодействия ($J_p$), но их роль в наблюдаемой динамике не обсуждается подробно. Учитывая трехслойную сотовую структуру, детальное исследование влияния межплоскостного сопряжения на магнитные возбуждения и общее основное состояние может быть ценным, особенно для понимания отклонений от чисто 2D китаевской физики.

Эти разнообразные перспективы, от экспериментального уточнения до продвинутого теоретического моделирования и проектирования материалов, имеют решающее значение для полного понимания сложного взаимодействия взаимодействий в RuBr3 и для продвижения поиска квантовых спиновых льдов Китаева.

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