Устойчивость тёмных солитонов в конденсате Бозе-Эйнштейна пузырьковой геометрии

Основная информация

• ID статьи: 2511.04385
• Название: Stability of dark solitons in a bubble Bose-Einstein condensate
• Авторы: Raphael Wictky Sallatti, Lauro Tomio, Dmitry E. Pelinovsky, Arnaldo Gammal
• Учреждения: Universidade de São Paulo, Universidade Estadual Paulista, McMaster University
• Классификация: cond-mat.quant-gas (физика конденсированного состояния - квантовый газ)
• Дата публикации: 7 ноября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2511.04385v1

Аннотация

В данной работе исследуется динамическая устойчивость тёмных солитонов, захваченных на поверхности двумерного сферического пузыря. В сферической геометрии конденсата Бозе-Эйнштейна (БЭК) при параметре взаимодействия ε ≳ 8,37 тёмные солитоны становятся неустойчивыми из-за змеевидной неустойчивости, вызванной дискретными угловыми модами, что приводит к образованию вихревых диполей. Исследование посредством аналитических и численных методов доказывает, что для каждой угловой моды m ≥ 2 существует ровно одна неустойчивая мода, доминирование которой определяет количество m вихревых диполей. Зависящие от времени моделирования подтверждают образование вихревых диполей.

Исследовательский контекст и мотивация

Исследуемая проблема

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

Значимость проблемы

1. Экспериментальная релевантность: Эксперименты в условиях микрогравитации на Международной космической станции реализовали сверхохлаждённые атомные газы в сферической геометрии; наземные эксперименты также реализовали структуры пузырьковых оболочек через двухкомпонентные системы
2. Фундаментальное физическое значение: Замкнутая сферическая топология предоставляет уникальную платформу для изучения низкомерных квантовых газов, включая фундаментальные физические явления, такие как динамика вихрей и фазовый переход Березинского-Костерлица-Таулеса
3. Геометрические эффекты: Кривизна сферической поверхности и замкнутая топология оказывают существенное влияние на распространение солитонов и неустойчивость, отличаясь от плоских систем

Ограничения существующих исследований

1. Отличие от плоских систем: В плоских квазиодномерных или квазидвумерных БЭК тёмные солитоны подвержены змеевидной неустойчивости и распадаются на вихри, однако замкнутая топология сферической поверхности изменяет путь распада
2. Топологические ограничения: На сферической поверхности не допускаются одиночные вихри; вихри должны появляться парами с нулевым полным зарядом (теорема Пуанкаре-Хопфа)
3. Отсутствие систематического анализа: Недостаёт систематического аналитического и численного исследования устойчивости тёмных солитонов в сферической геометрии

Исследовательская мотивация

Объединяя прогресс в экспериментах БЭК на МКС и теоретические потребности, систематически исследовать спектральную устойчивость тёмных солитонов на поверхности сферического пузыря, особенно то, как дискретные угловые моды вызывают неустойчивость и приводят к образованию вихревых пар.

Основные вклады

1. Определение порога устойчивости: Впервые точно определён порог устойчивости тёмных солитонов в БЭК сферического пузыря как ε ≈ 8,37
2. Аналитическая теоретическая база: Установлена полная схема анализа спектральной устойчивости Боголюбова-де Жена (БдЖ), доказывающая:
   • Для каждой угловой моды m ≥ 2 существует ровно одна неустойчивая мода
   • Приведена асимптотическая формула для порога неустойчивости: εₘᵗʰ = 4m(m-1)
3. Предсказание количества вихревых диполей: Установлена прямая соответствие между доминирующей неустойчивой модой m и количеством образующихся вихрь-антивихревых пар
4. Численная верификация: Полная верификация аналитических предсказаний посредством моделирования временной эволюции, демонстрирующая динамику распада тёмного солитона на m вихревых пар
5. Асимптотические решения в малом и большом пределах: Приведены аналитические выражения для профиля тёмного солитона в пределах ε → 0 и ε → ∞

Подробное описание методов

Определение задачи

Исследование линейной устойчивости решений типа тёмных солитонов в БЭК на двумерной сферической поверхности (радиус R, толщина δR << R), определение характеристик устойчивости различных угловых мод при разных интенсивностях взаимодействия.

Входные данные: Параметр взаимодействия ε = g/(2π) = 4√(2π)aN/(δR), где a — длина рассеяния s-волны, N — количество атомов Выходные данные: Спектр устойчивости, пороги неустойчивости, динамика образования вихревых пар

Математическая модель

1. Двумерное уравнение Гросса-Питаевского

В сферических координатах (θ, φ) безразмерное уравнение ГП имеет вид:

i∂ₜψ = -Δ₂Dψ + g|ψ|²ψ

где оператор Лапласа:

Δ₂D = (1/sinθ)∂θ(sinθ∂θ) + (1/sin²θ)∂²φ

2. Стационарное решение тёмного солитона

Ищется решение вида ψₛ(θ,φ,t) = f(θ)/√(2π) e^(-iμt), где f(θ) удовлетворяет:

-[d²/dθ² + cotθ d/dθ]f(θ) + ε|f(θ)|²f(θ) = μf(θ)

Условие нормировки: ∫₀^π sinθ|f(θ)|²dθ = 1

Характеристика тёмного солитона: f(θ) обращается в нуль при θ=π/2, монотонно убывает

3. Асимптотика малого параметра (ε → 0)

f(θ) = √(3/2)cosθ + εf₁(θ) + O(ε²)
μ = 2 + (9/10)ε + O(ε²)

где f₁(θ) ∝ P₃(cosθ) (полином Лежандра третьего порядка)

4. Асимптотика большого параметра (ε → ∞)

Тёмный солитон концентрируется вблизи экватора; введена переменная z = √(ε/2)(π/2 - θ):

f(θ) = f∞(θ) + f∞⁽¹⁾(θ)/√ε + O(1/ε)
f∞(θ) = (1/√2)tanh[√(ε/2)(π/2 - θ)]
μ = ε/2 + √ε + O(1)

Метод анализа устойчивости

1. Метод Боголюбова-де Жена

Рассматривается малое возмущение стационарного решения:

ψ(θ,φ,t) = ψₛ + (e^(-iμt)/√(2π))[u(θ,φ,t) + iv(θ,φ,t)]

После линеаризации получаются связанные уравнения:

∂ₜu = -Δ₂Dv + εf²v - μv
-∂ₜv = -Δ₂Du + 3εf²u - μu

2. Разложение по угловым модам

Возмущение разлагается по угловым модам: u = ûₘ(θ)e^(i(mφ+ωt)), v = v̂ₘ(θ)e^(i(mφ+ωt))

Получается спектральная задача:

ωûₘ = L⁻ₘv̂ₘ,  L⁻ₘ = -Δₘ + εf² - μ
ωv̂ₘ = L⁺ₘûₘ,  L⁺ₘ = -Δₘ + 3εf² - μ

где Δₘ = d²/dθ² + cotθ d/dθ - m²/sin²θ

3. Критерий устойчивости

• Если все собственные значения L±ₘ строго положительны, то все ω вещественны, тёмный солитон устойчив для моды m
• Если L⁻ₘ имеет отрицательное собственное значение, то появляются мнимые частоты Im(ω) ≠ 0, что вызывает неустойчивость

Технические инновации

1. Применение принципа сравнения

Используются два ключевых неравенства:

L⁺ₘ - L⁻ₘ = 2εf²(θ) ≥ 0
L±ₘ₊₁ - L±ₘ = (2m+1)/sin²θ ≥ 0

Строго доказана монотонность порога неустойчивости: εₘ < εₘ₊₁

2. Анализ классификации мод

• m=0: Существуют отрицательные собственные значения энергии, но численно неустойчивость не наблюдается
• m=1: Все собственные значения вещественны, нет ветвления неустойчивости
• m≥2: Для каждого m существует ровно одна неустойчивая мода

3. Техника асимптотического разложения

Для большого предела m посредством теории возмущений выведена:

εₘ = 4m(m-1) + O(1)

Данная формула имеет ошибку менее 5% при m=2

4. Комбинирование численных методов

• Метод стрельбы: Решение нелинейной краевой задачи для получения профиля тёмного солитона
• Конечные разности: Дискретизация оператора для вычисления собственных значений
• Разделённый метод Фурье: Комбинирование со схемой Крэнка-Николсона для временной эволюции

Экспериментальная установка

Методы численного решения

1. Вычисление профиля тёмного солитона

Применяется метод стрельбы в сочетании с методом секущей:

• Переопределение f̃(θ) = f(θ)√ε
• Граничные условия: f̃'(0) = 0, f̃(π/2) = 0
• Для заданного μ корректируется начальное значение f̃(0) до удовлетворения граничных условий
• Использование метода Рунге-Кутты для решения ОДУ
• Сканирование диапазона μ методом продолжения

2. Вычисление собственных значений

• Дискретизация оператора L±ₘ методом конечных разностей
• Количество узлов сетки: до 800 точек
• Вычисление нескольких наименьших собственных значений для определения устойчивости

3. Моделирование временной эволюции

Пространственная дискретизация:

• Направление θ: 256 узлов сетки, шаг δθ = π/256 ≈ 0,0123
• Направление φ: 256 узлов сетки, шаг δφ = 2π/256 ≈ 0,0245

Временная дискретизация:

• Временной шаг: δt = 10⁻⁵

Численный алгоритм:

• Быстрое преобразование Фурье (БПФ) для обработки направления φ
• Метод конечных разностей Крэнка-Николсона для направления θ
• Техника разделённого оператора для обработки нелинейного члена
• Полная схема: ψ(t+δt) ≃ e^(-ig|ψ|²δt/2) БПФ⁻¹ КН БПФ ψ(t)e^(-ig|ψ|²δt/2)

Установка параметров

Пороги устойчивости

Посредством численных расчётов определены критические значения для каждой моды:

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

Выбраны три типичных случая для демонстрации различных доминирующих мод:

• Доминирование m=2: ε=20
• Доминирование m=3: ε=50
• Доминирование m=4: ε=100

Экспериментальные результаты

Основные результаты

1. Фазовая диаграмма устойчивости

Ключевые находки:

• Область устойчивости: При ε ≲ 8,37 тёмный солитон устойчив для всех угловых мод
• Область неустойчивости: При ε ≳ 8,37 начинают возбуждаться моды m≥2
• Переход мод:
  • 8,37 < ε < 35: доминирует m=2
  • ε ≈ 35: конкуренция мод m=2 и m=3 Im(ω₂) = Im(ω₃)
  • ε > 35: начинает доминировать m=3

2. Эволюция собственных значений

На рисунке 2 представлены ключевые результаты:

• Наименьшее собственное значение L⁻ₘ (пунктирная линия): пересекает нулевую точку при ε=εₘ
• Im(ωₘ) (сплошная линия): начинает расти от нуля, характеризуя скорость роста неустойчивости
• Для заданного ε максимальное Im(ωₘ) соответствует доминирующей неустойчивой моде

3. Верификация теоретических предсказаний

Асимптотическая формула εₘᵗʰ = 4m(m-1) хорошо согласуется с численными результатами:

• m=2: теория 8 против численно 8,367 (ошибка 4,6%)
• m=3: теория 24 против численно 24,402 (ошибка 1,7%)
• m≥4: ошибка менее 1%

Анализ динамической эволюции

Случай 1: Доминирование m=2 (ε=20)

Характеристики временной эволюции:

• t=0: тёмный солитон расположен на экваторе, ширина ~1/√ε ≈ 0,22
• Промежуточный этап: появление змеевидной неустойчивости
• Финальное состояние: образование 2 вихрь-антивихревых пар

Распределение плотности:

• На экваторе плотность |ψ(π/2, φ)|² демонстрирует 2-кратную периодическую модуляцию
• 3D визуализация показывает 2 симметрично расположенные вихревые структуры

Случай 2: Доминирование m=3 (ε=50)

• Ширина тёмного солитона: ~0,14
• После распада образуется 3 вихрь-антивихревые пары
• Модуляция плотности демонстрирует 3-кратную симметрию

Случай 3: Доминирование m=4 (ε=100)

• Ширина тёмного солитона: ~0,10
• После распада образуется 4 вихрь-антивихревые пары
• Модуляция плотности демонстрирует 4-кратную симметрию

Проявление топологических ограничений

Верификация теоремы Пуанкаре-Хопфа:

• На сферической поверхности любое непрерывное касательное векторное поле должно иметь по крайней мере один нулевой точка
• Для сверхтекучести это означает невозможность существования одиночного вихря
• Вихри должны появляться парами с нулевой полной циркуляцией
• При возникновении неустойчивости количество вихрей изменяется на целое число, кратное +2

Характеристики профиля тёмного солитона

Верификация предела малого ε

На рисунке 1(a) верхний график показывает, что f(θ) близко к √(3/2)cosθ, согласуясь с теоретическим предсказанием (формула S10)

Верификация предела большого ε

На рисунке 1(a) верхний график показывает, что f(θ) концентрируется вблизи экватора, приближаясь к гиперболической функции тангенса f∞(θ)

Зависимость химического потенциала

На рисунке 1(b) представлена полная кривая μ(ε):

• Область малого ε: линейный рост μ ≈ 2 + 0,9ε (пунктирная линия)
• Область большого ε: рост по квадратному корню μ ≈ ε/2 + √ε + 1 (окружности)
• Численный результат (сплошная линия) плавно соединяет две асимптотические области на всём диапазоне

Связанные работы

Исследования БЭК в сферической геометрии

1. Экспериментальный прогресс:
   • Aveline и др. (2020): Первое наблюдение орбитального БЭК на МКС
   • Carollo и др. (2022): Наблюдение сверхохлаждённых атомных пузырей в условиях микрогравитации
   • Jia и др. (2022): Реализация пузырьковых БЭК в наземной двухкомпонентной системе
2. Теоретические основы:
   • Zobay & Garraway (2001-2004): Теория захвата атомов в полях, создаваемых адиабатическими потенциалами
   • Tononi & Salasnich (2019): Теория БЭК на сферической поверхности
   • Sun и др. (2018): Статические и динамические свойства конденсатов в оболочках

Динамика вихрей

1. Вихри на сферической поверхности:
   • Padavić и др. (2020): Физика вихрь-антивихревых пар в БЭК оболочек
   • Caracanhas и др. (2022): Динамика сверхтекучих вихрей на эллипсоидальных и поверхностях вращения
   • Tononi и др. (2024): Квантовые вихри в искривлённой геометрии
2. Топологические ограничения:
   • Tononi и др. (2022): Топологические фазовые переходы сверхтекучести в конденсатах, захваченных в пузырьковых ловушках
   • Подчёркивается ограничение теоремы Пуанкаре-Хопфа на конфигурации вихрей

Исследования тёмных солитонов

1. Плоские системы:
   • Burger и др. (1999): Экспериментальное наблюдение тёмных солитонов в БЭК
   • Busch & Anglin (2000): Движение тёмных солитонов в захваченных БЭК
   • Feder и др. (2000): Состояния тёмных солитонов в анизотропных ловушках
2. Змеевидная неустойчивость:
   • Gaidoukov & Anglin (2021): БдЖ теория змеевидной неустойчивости серых солитонов в высших размерностях
   • Anderson и др. (2001): Наблюдение распада тёмного солитона на кольцо вихрей
3. Кольцевая геометрия:
   • Kivshar & Yang (1994): Тёмные солитоны в кольцевой геометрии
   • Theocharis и др. (2003): Кольцевые тёмные солитоны и цепочки вихрей в БЭК

Уникальный вклад данной работы

По сравнению с существующими работами, данная статья:

1. Впервые систематически исследует полный спектр устойчивости тёмных солитонов в БЭК сферического пузыря
2. Устанавливает аналитическую базу для точного предсказания порога неустойчивости каждой угловой моды
3. Количественно связывает доминирующую неустойчивую моду с количеством вихревых пар
4. Полностью верифицирует теоретические предсказания на всём диапазоне параметров от малого ε к большому ε

Заключение и обсуждение

Основные выводы

1. Порог устойчивости: Тёмные солитоны в БЭК сферического пузыря устойчивы при ε ≲ 8,37; превышение этого значения вызывает змеевидную неустойчивость угловых мод m≥2
2. Механизм неустойчивости:
   • Каждая мода m≥2 имеет ровно одно неустойчивое собственное значение
   • Пороги неустойчивости удовлетворяют строгой монотонности εₘ < εₘ₊₁
   • Асимптотическая формула εₘ ≈ 4m(m-1) точно предсказывает поведение при большом m
3. Образование вихревых пар: Доминирующая неустойчивая мода m прямо определяет количество образующихся вихрь-антивихревых пар; моделирование временной эволюции полностью подтверждает это предсказание
4. Топологические ограничения: Замкнутая топология сферической поверхности требует появления вихрей парами; вихревые пары не могут убежать на границу, что приводит к сложной динамике вихрей

Ограничения

1. Теоретические предположения:
   • Применяется двумерное приближение, требующее δR << R
   • Игнорируются радиальные возбуждения (более высокие энергии)
   • Не учитываются эффекты тепловых облаков, вызывающие энергетическую неустойчивость
   • Приближение среднего поля, игнорирование квантовых флуктуаций
2. Мода m=0: Теория предсказывает возможную неустойчивость, но численно это не наблюдается; требуется более глубокое исследование
3. Экспериментальные условия:
   • В реальных экспериментах сложно точно контролировать параметр ε
   • Время проведения экспериментов в условиях микрогравитации ограничено
   • Необходимо учитывать эффекты конечной температуры
4. Долгосрочная динамика:
   • Исследуется только начальная неустойчивость и образование вихревых пар
   • Не проводится глубокий анализ последующей эволюции вихревых пар (аннигиляция, стабильные конфигурации и т.д.)

Направления будущих исследований

1. Экспериментальная верификация:
   • Использование платформы МКС или наземных двухкомпонентных систем для верификации порогов устойчивости
   • Наблюдение соответствия между различными значениями ε и количеством вихревых пар согласно теоретическим предсказаниям
2. Эффекты конечной температуры:
   • Включение влияния теплового облака на устойчивость тёмных солитонов
   • Исследование конкуренции между энергетической и динамической неустойчивостью
3. Динамика вихрей:
   • Глубокое исследование долгосрочной эволюции вихревых пар после их образования
   • Анализ механизмов аннигиляции вихрь-антивихревых пар и формирования стабильных конфигураций
4. Двухкомпонентные системы:
   • Расширение на двухкомпонентные смеси БЭК
   • Исследование устойчивости заполненных солитонов (filled solitons) с усиленной стабильностью
5. Эллипсоидальная геометрия:
   • Обобщение на неферические замкнутые поверхности
   • Исследование влияния геометрической асимметрии на устойчивость

Глубокая оценка

Преимущества

1. Теоретическая строгость

• Полная математическая база: От уравнения ГП через метод БдЖ к полной теории линейной устойчивости
• Строгие доказательства: Использование принципа сравнения и теории Штурма для строгого доказательства существования и единственности неустойчивости
• Асимптотический анализ: Асимптотические разложения в пределах малого и большого ε предоставляют аналитические инсайты

2. Продвинутые численные методы

• Многомасштабная обработка: Метод стрельбы + метод продолжения эффективно решают нелинейные краевые задачи
• Комбинирование спектральных методов: БПФ для периодического направления, конечные разности для полюсных областей
• Высокая точность: Численные результаты согласуются с теоретическими предсказаниями с ошибкой менее 5%

3. Ясная физическая картина

• Классификация мод: Полный анализ трёх классов мод m=0, 1, ≥2
• Топологические ограничения: Ясное объяснение физического смысла теоремы Пуанкаре-Хопфа
• Визуализация: 3D графики плотности наглядно демонстрируют процесс образования вихревых пар

4. Экспериментальная релевантность

• Управляемые параметры: Единственный параметр ε содержит всю физическую информацию (длина рассеяния, количество атомов, толщина)
• Реализуемость: Согласуется с условиями экспериментов на МКС и наземных установках
• Чёткие предсказания: Даёт проверяемые количественные предсказания (пороги, количество вихревых пар)

5. Подробные дополнительные материалы

• 19 страниц дополнительных материалов с полными выводами
• Включены все детали реализации численных методов
• Облегчает понимание и воспроизведение результатов

Недостатки

1. Отсутствие экспериментальной верификации

• Чисто теоретическое исследование без сравнения с экспериментальными данными
• Не обсуждаются конкретные схемы и вызовы экспериментальной реализации
• Диапазон управляемости параметра ε в экспериментах не уточнен

2. Неполная долгосрочная динамика

• Демонстрируется только этап до образования вихревых пар (короткое время)
• Последующая эволюция вихревых пар (аннигиляция, перестройка) не исследуется глубоко
• Финальные стабильные конфигурации не определены

3. Эффекты конечной температуры

• Приближение нулевой температуры, игнорирование тепловых возбуждений
• В реальных экспериментах эффекты температуры могут быть значительными
• Связь между энергетической и динамической неустойчивостью не обсуждается

4. Неполный анализ моды m=0

• Теория предсказывает возможную неустойчивость, но численно не наблюдается
• Отсутствует глубокое объяснение этого противоречия
• Может потребоваться более тонкое численное вычисление или теоретическая коррекция

5. Разрыв с экспериментальными условиями

• Предположение идеальной сферической формы; реально может быть эллиптичность
• Игнорирование потерь атомов и эффектов нагрева
• Влияние внешних возмущений (остаточная гравитация) не рассматривается

Оценка влияния

Академическая ценность

• Теоретическая полнота: Предоставляет базовую теорию для устойчивости солитонов в БЭК сферической геометрии
• Методологический вклад: Систематическое применение метода БдЖ в искривлённом пространстве
• Междисциплинарность: Связывает теорию нелинейных волн, топологию и физику холодных атомов

Практическая ценность

• Экспериментальное руководство: Даёт чёткие параметрические рекомендации для экспериментов на МКС и наземных установках
• Предсказание явлений: Количественное предсказание количества вихревых пар может быть прямо проверено
• Техническая справка: Численные методы применимы к другим задачам в искривлённой геометрии

Воспроизводимость

• Высокая: Дополнительные материалы содержат полные детали алгоритмов
• Численные методы стандартные (БПФ, Крэнка-Николсона, метод стрельбы)
• Параметры установки чётко определены, легко воспроизводятся

Потенциальное влияние

1. Краткосрочное: Стимулирование экспериментальных работ по верификации
2. Среднесрочное: Стимулирование расширенных исследований (двухкомпонентные системы, конечная температура)
3. Долгосрочное: Предоставление новой платформы для квантового моделирования и топологической физики

Применимые сценарии

Прямое применение

1. Эксперименты БЭК в микрогравитации: Лаборатория холодных атомов на МКС
2. Наземные пузырьковые БЭК: Двухкомпонентные смешанные системы
3. Численное моделирование: Теоретические исследования БЭК в сферической геометрии

Сценарии обобщения

1. Другие замкнутые поверхности: Эллипсоиды, торы и т.д.
2. Двухкомпонентные системы: Заполненные тёмные солитоны, спин-орбитальная связь
3. Другие квантовые жидкости: Сверхтекучий гелий, экситон-поляритоны

Условия ограничения

• Действительность двумерного приближения: δR/R < 0,1
• Приближение нулевой температуры: T << Tc
• Действительность среднего поля: Na³ >> 1
• Сферическая симметрия

Ключевые ссылки

1. Экспериментальные основы:
   • Carollo et al., Nature 606, 281 (2022) - Наблюдение пузырьковых БЭК на МКС
   • Jia et al., PRL 129, 243402 (2022) - Наземные пузырьковые БЭК
2. Теоретическая база:
   • Geyer & Pelinovsky, AMS Mathematical Surveys 288 (2025) - Устойчивость нелинейных волн в гамильтоновых системах
   • Andriati et al., PRA 104, 033318 (2021) - Устойчивость смесей в пузырьковых ловушках
3. Физика тёмных солитонов:
   • Burger et al., PRL 83, 5198 (1999) - Первое наблюдение тёмных солитонов в БЭК
   • Anderson et al., PRL 86, 2926 (2001) - Распад тёмного солитона на кольцо вихрей
4. БЭК в сферической геометрии:
   • Tononi & Salasnich, PRL 123, 160403 (2019) - Теория БЭК на сферической поверхности
   • Dubessy & Perrin, AVS Quantum Sci. 7, 010501 (2025) - Обзор квантовых газов в пузырьковых ловушках

Общая оценка: Это высокого качества теоретическая работа с строгой математической базой, продвинутыми численными методами и ясной физической картиной. Посредством идеального сочетания аналитических и численных методов систематически решена фундаментальная проблема устойчивости тёмных солитонов в БЭК сферической геометрии. Асимптотическая формула εₘ ≈ 4m(m-1) элегантна и проста; соответствие между доминирующей модой и количеством вихревых пар имеет чёткую экспериментальную проверяемость. Работа предоставляет важное теоретическое руководство для текущих экспериментов БЭК в условиях микрогравитации и методологическую справку для исследований нелинейных волн в искривлённой геометрии. Основные недостатки заключаются в отсутствии экспериментальной верификации и рассмотрения эффектов конечной температуры, однако это естественные направления для последующих исследований.