Столетняя загадка струн Бете наконец решена

Большинство научных изысканий находят свое экспериментальное подтверждение незамедлительно либо по прошествии непродолжительного времени. Тем не менее, научное сообщество прошлых эпох породило колоссальный массив концепций, долгое время остававшихся лишь гипотезами из-за несовершенства технологической базы. Часть из них была опровергнута современными исследователями, иные по сей день остаются спорными, однако некоторые все же дождались своего триумфального подтверждения. К последней категории относится гипотеза, сформулированная физиком Хансом Бете еще в 1931 году и получившая название «струны Бете». Сотрудники Инсбрукского университета (Австрия) сумели искусственно сформировать и зафиксировать эти уникальные квантовые состояния. Каким образом физикам удалось воплотить эксперимент, какими свойствами обладали полученные структуры и какое значение это открытие имеет для современной квантовой механики — рассмотрим в данной публикации.

Основа исследования

«Струны Бете» представляют собой связанные конфигурации элементарных частиц в рамках одномерных интегрируемых квантовых многочастичных систем, играющих ключевую роль в магнетизме, физике ультрахолодных атомов и наноструктур. Они выступают в качестве квантового эквивалента ярких солитонов. Будучи фундаментальными коллективными возбуждениями, согласно теоретическим выкладкам, эти структуры обязаны радикально трансформировать равновесные и динамические характеристики сред. Если ранее короткие единичные струны Бете уже удавалось обнаружить в квантовых магнитах и сверхпроводящих кубитах, то синтез комплексных ансамблей, состоящих из множества разноразмерных (включая протяженные) струн, долгое время казался невыполнимой экспериментальной задачей. В новом исследовании этот барьер был успешно преодолен.

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

В интегрируемых одномерных квантовых моделях силы притяжения вовсе не обязательно приводят к катастрофическому коллапсу. Вместо этого они порождают «струны Бете» — многочастичные связанные состояния, предсказанные около века назад. Они зарождаются как стабильные возбуждения внутри спиновых цепочек, квантовых магнитов и газов холодных атомов. Их стабильность обусловлена характерным для интегрируемых систем отсутствием дифракционного рассеяния, что надежно защищает их от распада в низкоэнергетические состояния. Последние изыскания демонстрируют, что протяженные струны Бете определяют аномальные транспортные характеристики интегрируемых спиновых цепочек и формируют универсальное поведение в рамках класса Кардара-Паризи-Жанга. Помимо фундаментальной физики конденсированных сред, эти объекты находят важнейшее применение в теориях струн и высоких энергий. Глубокое постижение столь широкого спектра проявлений требует безупречного владения методами синтеза и контроля плотных ансамблей струн Бете различного масштаба.

Предполагается, что в системах холодных атомов, описываемых моделью Либа-Линигера, также зарождаются струны Бете. В квазиклассическом пределе они трансформируются в яркие солитоны нелинейного уравнения Шредингера. Солитонные цепочки, сохраняющие стабильность исключительно при слабых взаимодействиях, ранее уже наблюдались в бозонных средах с притяжением вблизи квазиклассического режима. Здесь формирование солитонов интерпретируется как результат модуляционной неустойчивости уравнения Шредингера. В описываемой работе ученые задействовали ультрахолодные атомы цезия с эффектом притяжения в рамках одномерной геометрии, сумев сформировать плотный гибрид струн путем плавной регулировки межчастичного взаимодействия от отталкивания к притяжению через точку нулевого взаимодействия. Исследователи функционировали в квантовой области, лежащей за пределами применимости уравнения Шредингера, при этом стремительный коллапс газовой среды сдерживался благодаря приближенной интегрируемости системы. За счет контролируемого изменения параметров было создано глубоко неравновесное состояние, насыщенное множеством струн Бете, включая масштабные протяженные конфигурации. Их габариты определялись исходными условиями — главным образом плотностью и температурными показателями атомного облака. Фиксация этих объектов осуществлялась посредством сопоставления данных, полученных после деактивации удерживающего потенциала в 1D- и 3D-форматах, а также через оценку импульсных распределений и вычисление параметра контакта Тана. Полученные эмпирические результаты демонстрируют прекрасное совпадение с теоретическими прогнозами обобщенной гидродинамики.

Результаты исследования

Струны Бете проявляют себя как собственные состояния классического гамильтониана Либа-Линигера:

где первый член отражает кинетическую энергию частиц массой m, а второй симулирует парные взаимодействия интенсивностью c̃. Собственные состояния, вычисленные посредством анзаца Бете, задаются набором быстрот θ, играющих роль обобщенных импульсов возникающих стабильных квазичастиц.

Изображение №1

На панели 1a приведена принципиальная схема различных режимов взаимодействия одномерного бозе-газа. При c̃ > 0 преобладают отталкивающие силы, а быстроты принимают чисто действительные значения. В условиях мощного отталкивания бозоны фермионизируются, переходя в состояние Тонкса-Жирардо. Сдвиг в сторону сильного притяжения (c̃ < 0) инициирует возбужденное газовое состояние, известное как супер-Тонкс-Жирардо. Тем не менее, в зоне притяжения быстроты приобретают комплексный характер с ненулевыми мнимыми компонентами. Они концентрируются в кластерах с идентичной вещественной частью λ ≡ Re(θ), симметрично распределенных по мнимой оси с шагом |c̃|. Эти специфические образования и есть струны Бете. Автономная струна Бете ранга n обладает волновой функцией ψ(z₁ < . . . < z_n) ∝ exp(iΣⁿ_j=1θ_jz_j) с быстротами θ_j = λ + |c̃|(n + 1 − 2j). Протяженность такой структуры Δz_n оценивается по затуханию волновой функции как Δz_n = 2[|c̃|(n − 1)]⁻¹, где масштабом выступает величина |c̃|⁻¹. Полная энергия струны E_n складывается из энергии центра масс E_com;n = nℏ²λ²/(2m) и энергии связи E_b;n = ℏ²c̃²n(n² − 1)/(24m), давая в разности E_com;n − E_b;n итоговый показатель. Фундаментальное основное состояние для N частиц представляет собой колоссальный комплекс Бете. Малокомпонентные струны Бете должны рождаться в неравновесной обстановке, которая в данном опыте моделировалась за счет прохождения системы с конечной скоростью сквозь точку нулевого взаимодействия из фазы отталкивания в фазу притяжения. Специфика предстартового состояния, динамика изменения параметров и неизбежная конечная температура определяют итоговый спектр распределения состояний Бете в подобных условиях.

Экспериментальный протокол

Экспериментальный цикл стартовал с получения трехмерного конденсата Бозе-Эйнштейна атомов цезия, помещенного в оптическую дипольную ловушку с пересекающимися лучами, где гравитация компенсировалась магнитным полем. Общее число атомов в конденсате варьировалось от 2 × 10⁴ до 6 × 10⁴ при 10-процентной погрешности, а термодинамическая температура удерживалась в диапазоне 10–35 нК с точностью до 1 нК. Использование широкого резонанса Фешбаха позволило прецизионно управлять длиной s-рассеяния a_3D при помощи подстраиваемого магнитного поля B_z. Первоначальное значение поля выставлялось на отметке 20.8 Гс, что соответствует a_3D = 209 a₀, при этом конденсат пребывал в режиме Томаса-Ферми. Среда адиабатически переносилась в скрещенные под прямым углом лазерные лучи, формирующие двумерную оптическую решетку (1b). Это приводило к образованию массива изолированных 1D-трубок с поперечной и продольной частотами удержания ω_⊥/2π = 10.5 кГц и ω_z/2π = 29.3 Гц соответственно (при глубине решетки 25 E_r, где E_r — энергия отдачи, а d = 532.2 нм — период решетки). Данные параметры соответствовали длинам гармонического осциллятора a_ho = 1.6 мкм (вдоль) и a_⊥ = 83 нм (поперек). Варьируя плотность исходного конденсата, исследователи добивались среднего числа частиц в каждой трубке от 7 до 30, что обеспечивало заполнение от 3000 до 35000 аналогичных каналов.

В рассматриваемой одномерной модели интенсивность межчастичной связи задается соотношением c̃ ≃ 4a_3D/a²_⊥. По завершении этапа загрузки решетки параметр взаимодействия устанавливался на уровне c̃ = 3.4 мкм⁻¹, переведя систему в умеренно взаимодействующий режим Тонкса-Жирардо. Одномерная температура T_1D лежала в пределах 3–35 нК. В данной области параметров система находилась в глубоком квантовом режиме, где приближение Гросса-Питаевского становится несостоятельным, а обобщенная гидродинамика (GHD) выступает единственным адекватным инструментом анализа. Экспериментаторы плавно снижали параметр взаимодействия в сторону отрицательных значений, пересекая точку нулевого взаимодействия с постоянной скоростью ∂c̃/∂t ≈ −0.038 мкм⁻¹мс⁻¹. Ожидалось, что эта процедура трансформирует систему в ансамбль струн Бете по мере нарастания сил притяжения. В частности, при c̃ = −1 мкм⁻¹ струна ранга n = 4 обладала длиной связи порядка 0.67 мкм и энергией связи, эквивалентной k_B × 9 нК. Стоит подчеркнуть, что струны Бете корректно определены лишь тогда, когда их размер Δz_n существенно превышает поперечный масштаб a_⊥. На сверхмалых дистанциях происходят неизбежные потери за счет короткодействующих молекулярных реакций, не охватываемых моделью Либа-Линигера. Отсутствие возбуждений поперечных степеней свободы в методе времени пролета (TOF) подтвердило сохранение строгой одномерности системы.

Для исследования характеристик среды применялся протокол измерения быстрот, задействованный ранее для систем без связанных комплексов. Он включает двухступенчатую процедуру расширения: начальное разлетание внутри трубок на интервал t_1D с последующим свободным трехмерным разлетом в течение t_3D (1c). Время t_1D подбиралось достаточно большим для перехода среды в разреженный режим, где взаимодействия нивелируются. В фазе отталкивания последующее абсорбционное изображение напрямую отражает спектр быстрот газа. Однако в режиме притяжения связь между распределением быстрот и оптическим снимком усложняется. На практике продольное расширение в трубках запускалось одновременным отключением дипольной ловушки и устранением остаточного продольного потенциала решетки с помощью вспомогательного синего луча. На втором этапе снималось поперечное ограничение и одновременно обнулялось a_3D, исключая динамическое влияние взаимодействий на плотность разлетающегося облака.

Измерения энергии высвобождения

Изображение №2

Первоочередной задачей эксперимента являлась фиксация следов «струн Бете» посредством регистрации энергии высвобождения. Величины E_1D и E_3D вычислялись через анализ эволюции второго момента разлетающегося облака на соответствующих этапах протокола. Параметр E_1D определялся варьированием времени t_1D при фиксированном t_3D, тогда как E_3D — наоборот. Графики зависимости этих энергий от параметра c̃ наглядно демонстрируют, что в отталкивающем режиме обе величины совпадают в пределах погрешности, плавно снижаясь от ~h × 100 Гц при сильном отталкивании до ~h × 60 Гц около нулевой точки. Напротив, в зоне притяжения возникает качественное расхождение: если E_1D выходит на плато около h × 35 Гц, то E_3D стремительно возрастает до h × 110 Гц при c̃ = −1.0 мкм⁻¹. Данная аномалия служит неопровержимым доказательством генерации струн Бете.

В ходе 1D-расширения эти связанные комплексы претерпевают соударения, но сохраняют целостность, из-за чего показатель E_1D остается низким. Переход к 3D-разлету и мгновенное снятие потенциала приводят к разрушению связанных состояний, высвобождая энергию связи в виде кинетической энергии (преимущественно вдоль оси z), что и обусловливает всплеск E_3D.

В отталкивающей фазе отсутствие связанных структур гарантирует тождество E_1D = E_3D = U, где U — внутренняя энергия газа. В режиме притяжения величина E_1D отражает кинетику движения центров масс струн (⟨E_com⟩), в то время как E_3D суммирует кинетическую энергию центров масс и энергию связи (⟨E_com⟩ + ⟨E_b⟩).

Результаты численного моделирования в рамках GHD демонстрируют превосходное соглашение с экспериментом. Несмотря на конечность числа частиц, гидродинамический подход обеспечивает исчерпывающее количественное описание. Главным фактором, нарушающим идеальную интегрируемость в реальном эксперименте, выступают трехчастичные потери, не учитываемые теорией GHD. Если при отталкивании они ничтожно малы, то по мере уплотнения струн Бете в зоне притяжения трехчастичные процессы становятся заметными, приводя к потере наиболее крупных и энергоемких струн. Это обстоятельство объясняет небольшие расхождения между моделью и экспериментом.

На вставке к изображению №2 показана зависимость 〈E_b〉 от c̃, полученная из разности E_3D − E_1D. Без привлечения громоздких расчетов GHD эти данные позволяют установить нижнюю границу максимальных размеров струн Бете в системе. Поскольку энергия связи монотонно растет с увеличением числа компонентов в струне, наблюдаемый уровень 〈E_b〉 требует присутствия комплексов размером по меньшей мере n* = 6, хотя детальный анализ GHD предсказывает генерацию еще более массивных структур.

Характеристика одномерного расширения

Изображение №3

Дополнительные доказательства существования струн Бете были получены при сопоставлении импульсных распределений f(p/ℏ) после одномерного расширения с прогнозами GHD (при t_1D = 10 мс и t_3D = 46.3 мс). Импульсные профили извлекались из продольного распределения облака по формуле p = mz/(t_1D + t_3D). На панели 3a отражена эволюция f(p/ℏ) при переходе от отталкивания к притяжению. Изначально функция имеет колоколообразную форму с полушириной (FWHM) ℏ × 6.0 мкм⁻¹. При приближении к нулевому взаимодействию пик заостряется, сужаясь до ℏ × 3.0 мкм⁻¹ при c̃ = 0.1 мкм⁻¹. После пересечения нулевой отметки распределение вновь уширяется, увеличиваясь почти двукратно к концу диапазона при c̃ = −2.0 мкм⁻¹.

Сопоставление с расчетами GHD, где оптимальная начальная температура была оценена как T_1D = 15 нК, выявило отличную корреляцию во всей области параметров. В отталкивающем режиме f(p/ℏ) эквивалентна распределению быстрот, поэтому сужение профиля объясняется замедлением квазичастиц. Напротив, в режиме притяжения квазичастицы объединяются в струны Бете, что приводит к уширению импульсного спектра. Согласно принципу неопределенности, струна с длиной Δz_n генерирует неопределенность импульса Δp_n ~ ℏ/(2Δz_n), причем этот эффект усиливается с ростом размера струн и глубины притяжения.

Влияние температуры на формирование комплексов проиллюстрировано на панели 3b, где представлена зависимость FWHM от c̃ для трех различных значений T_1D. С ростом температуры уширение при удалении от точки c̃ = 0 сглаживается. На панели 3c приведены предсказания GHD относительно доли атомов nP_n, задействованных в цепочках ранга n. Повышение температуры приводит к доминированию кинетической энергии над энергией связи, подавляя рождение массивных кластеров и снижая среднюю энергию связи на одну частицу.

Контакт Тана

Присутствие струн Бете кардинально трансформирует короткодействующие парные корреляции среды G₂ = ⟨Σ_i≠jδ(z_i − z_j)⟩, определяющие вероятность обнаружения двух частиц в одной точке. Данный параметр прямо пропорционален контакту Тана: C = c̃²G₂, который устанавливает мост между микроскопической структурой и термодинамикой системы. Теорема вириала связывает контакт C с внутренней энергией U и потенциальной энергией ловушки E_V следующим соотношением:

Как выяснилось, данное уравнение корректно и для неравновесных сценариев, описываемых стационарными решениями GHD.

Экспериментальное определение U базировалось на измерениях энергии 3D-высвобождения. Если в отталкивающем режиме полная энергия находилась напрямую из E_3D, то для притяжения использовалось приближение 〈E_com⟩(c̃ < 0) ≃ E_3D(c̃ = 0), а энергия связи оценивалась как разность E_3D − E_3D(c̃ = 0). Потенциальная энергия измерялась по пространственному профилю облака при коротком TOF (1.5 мс) без предварительного 1D-расширения, составляя E_V = (1/2)mω_z²〈z²〉.

Изображение №4

На графике выше представлена зависимость параметра контакта от сил взаимодействия. В отталкивающей зоне контакт Тана ожидаемо снижается по мере уменьшения c̃, однако при переходе в область притяжения он демонстрирует резкий скачок, оказываясь почти в шесть раз выше при |c̃| = 1 мкм⁻¹. Этот эффект обусловлен уплотнением струн Бете и ростом локальной плотности G₂. Теоретические кривые GHD качественно описывают эксперимент, демонстрируя идеальное совпадение при отталкивании и небольшие расхождения в зоне притяжения, обусловленные эффектами конечного размера системы и потерей части атомов (до 20% в сильнейших полях притяжения). Полученные данные могут послужить надежным бенчмарком для верификации альтернативных теоретических подходов.

Для углубленного изучения методологии рекомендуем ознакомиться с оригинальным докладом исследователей и дополнительными материалами.

Эпилог

В представленном исследовании физикам впервые удалось экспериментально сформировать и зафиксировать трудноуловимые квантовые состояния — струны Бете, названные в честь выдающегося физика, предсказавшего их существование еще в 1931 году.

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

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

В перспективе научный коллектив планирует исследовать роль струн Бете в формировании универсального поведения по системе Кардара-Паризи-Жанга, а также изучить их поведение при размерных кроссоверах и их потенциальную связь с эффектами физики Ефимова.

Немного рекламы

Спасибо, что остаетесь с нами. Вам нравятся наши статьи? Хотите видеть больше интересных материалов? Поддержите нас, оформив заказ или порекомендовав знакомым, облачные VPS для разработчиков от $4.99, уникальный аналог entry-level серверов, который был придуман нами для Вас: Вся правда о VPS (KVM) E5-2697 v3 (6 Cores) 10GB DDR4 480GB SSD 1Gbps от $19 или как правильно делить сервер? (доступны варианты с RAID1 и RAID10, до 24 ядер и до 40GB DDR4).

Dell R730xd в 2 раза дешевле в дата-центре Maincubes Tier IV в Амстердаме? Только у нас 2 х Intel TetraDeca-Core Xeon 2x E5-2697v3 2.6GHz 14C 64GB DDR4 4x960GB SSD 1Gbps 100 ТВ от $199 в Нидерландах! Dell R420 — 2x E5-2430 2.2Ghz 6C 128GB DDR3 2x960GB SSD 1Gbps 100TB — от $99! Читайте о том Как построить инфраструктуру корп. класса c применением серверов Dell R730xd Е5-2650 v4 стоимостью 9000 евро за копейки?

 

Источник