Открытие пространства Крылова для доступа ко всей динамике посредством динамических симметрий

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

• ID статьи: 2503.07403
• Название: Opening Krylov space to access all-time dynamics via dynamical symmetries
• Авторы: Nicolas Loizeau, Berislav Buča, Dries Sels
• Учреждения: Niels Bohr Institute (Копенгаген), Université Paris-Saclay, University of Oxford, New York University, Flatiron Institute
• Классификация: quant-ph (квантовая физика)
• Дата публикации: 22 августа 2025 г. (версия v2)
• Ссылка на статью: https://arxiv.org/abs/2503.07403v2

Аннотация

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

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

Основные проблемы

1. Вычислительные трудности квантовой многотельной динамики: Решение краткосрочной и долгосрочной динамики закрытых квантовых многотельных систем является одной из главных проблем атомной и конденсированной материи физики
2. Сложность извлечения динамических симметрий: Хотя теоретически динамика локальных наблюдаемых может быть полностью охарактеризована собственными модами Лиувиллиана (динамическими симметриями), вывод этих операторов оставался невозможным за исключением редких моделей

Важность проблемы

• Полная характеризация динамики: Динамические симметрии обеспечивают полную характеризацию динамики системы; для произвольного оператора O и начального состояния ρ временная эволюция может быть представлена как: Tr(O(t)ρ) = ∑_ω e^(iωt)μ_ω
• Понимание неравновесных явлений: Критично для понимания многотельной локализации, временных кристаллов, квантовых рубцов, фрагментации и других нетривиальных динамических поведений
• Механизмы термализации: В хаотических системах связывает гипотезу роста коррелированных операторов с термализацией и экспоненциальным затуханием наблюдаемых

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

1. Ограничения конечных систем: Спектр Лиувиллиана конечной закрытой системы чисто вещественный, что не позволяет исследовать свойства термодинамического предела
2. Методы разложения по псевдомодам: Предыдущие методы (Refs 96,97) ограничены динамикой вблизи равновесия
3. Неисследованная псевдолокальность: Предыдущие исследования не рассматривали (псевдо)локальность динамических симметрий, что является ключевым для неравновесной физики

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

Предложить более строгий метод путём введения открытых граничных условий на цепи Крылова для реализации разделения система-окружение, позволяющий:

• Локальным величинам поддерживаться в системе
• Нелокальным величинам составлять окружение
• Идентифицировать, какие динамические симметрии являются псевдолокальными и, следовательно, релевантными

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

1. Новый метод открытых граничных условий: Предложено введение диссипативных граничных условий в конкретной позиции L цепи Крылова, избегая нефизического отражения волновой функции, с формой: ∂t φ_L = (b_L + b{L+1})φ_ - 2b_{L+1}φ_L
2. Структура разделения система-окружение: Установлено естественное разложение пространства Крылова, где:
   • Левая сторона (система): простые локальные операторы с точно известными коэффициентами прыжков
   • Правая сторона (окружение): сложные нелокальные операторы, составляющие эффективную ванну
3. Аналитические решения для идеальных случаев: Получены точные решения для двух идеализированных случаев:
   • Линейный рост b_n = n: собственные векторы выражены через полиномы Мейксера
   • Рост по квадратному корню b_n = √n: собственные векторы выражены через полиномы Эрмита
4. Прямая связь между ростом операторов и термализацией: Доказано, что в случае линейного роста (хаотические системы) все наблюдаемые экспоненциально затухают со скоростью -2α, где α — скорость роста коэффициентов Ланцоша
5. Численная верификация и приложения: Метод верифицирован на реальных спиновых системах, таких как XXZ-цепь и трансверсальная цепь Изинга, с успешным извлечением динамических симметрий

Детальное описание метода

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

Входные данные:

• Гамильтониан H
• Начальный оператор O_0 (или начальная матрица плотности ρ)

Выходные данные:

• Динамические симметрии A_ω, удовлетворяющие H, A_ω = -ωA_ω
• Частоты ω (могут быть комплексными): вещественная часть соответствует осцилляциям, мнимая часть — затуханию

Ограничения:

• Рассматриваются только псевдолокальные динамические симметрии (влияющие на локальную физику)
• Применимо к системам с локальными взаимодействиями и начальными состояниями с короткодействующими корреляциями

Архитектура модели

1. Построение пространства Крылова (алгоритм Ланцоша)

Начиная с затравочного оператора O_0, рекурсивно применяется Лиувиллиан L = H, · для построения ортонормированного базиса {O_n}:

O_1 = LO_0/b_1 = [H,O_0]/b_1,  b_1 = ||LO_0||

Для n > 2:
O'_n = LO_{n-1} - b_{n-1}O_{n-2}
O_n = O'_n/b_n
b_n = ||O'_n||

где норма определена как ||O||² = 1/(2N) TrO²

2. Одномерное отображение одной частицы

Временная эволюция в базисе Крылова отображается на задачу одномерного прыгания одной частицы:

O(t) = 1/(2N) ∑_n i^n φ_n(t)O_n

Уравнение эволюции коэффициентов: ∂t φ_n = b_n φ - b_{n+1} φ_{n+1}, φ_n(0) = δ_

3. Открытые граничные условия

Цепь обрезается в позиции L, предполагая, что φ_n достаточно гладкая для локальной аппроксимации линейной функцией: φ_{L+1} ≈ φ_L + (φ_L - φ_)

Получаются диссипативные граничные условия: ∂t φ_L = (b_L + b{L+1})φ_ - 2b_{L+1}φ_L

4. Модифицированная матрица Лиувиллиана

Обрезанный Лиувиллиан имеет трёхдиагональную форму плюс модификация граничных условий:

L̃ = i × [трёхдиагональная матрица, последняя строка имеет специальную форму]
Последняя строка: [0, 0, ..., b_L + b_{L+1}, -2b_{L+1}]

Эта матрица неэрмитова, что позволяет иметь комплексные собственные значения

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

1. Физический смысл граничных условий

• Отсутствие введения параметров: В отличие от Refs 96,97 с внешней диссипацией и экстраполяцией, данный метод не требует дополнительных параметров
• Предположение гладкости: Основано на естественном предположении о непрерывности производной волновой функции Крылова
• Разделение локального и нелокального: Операторы справа от цепи Крылова в основном состоят из k-локальных строк Паули с k > L

2. Отличие от традиционных методов

• Граничные условия Дирихле: φ_{L+1} = 0 приводит к нефизическому отражению
• Данный метод: Позволяет потоку входить в окружение и возвращаться, более реалистично отражая распространение оператора

3. Обоснованность конструкции

В идеальных случаях (b_n = n или √n):

• Динамика при открытых граничных условиях практически неотличима от бесконечной цепи (рис. 2)
• Вещественная часть спектра почти не изменяется, мнимая часть вводит физическое затухание (рис. 3)
• Случай линейного роста: все нетривиальные корни удовлетворяют Im(ω) = -2i (численно верифицировано до n=2^12)

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

Модельные системы

1. XXZ-цепь (интегрируемая модель)

H = ∑_i (s^x_i s^x_{i+1} + s^y_i s^y_{i+1} + Δs^z_i s^z_{i+1} + hs^z_i)
Параметры: Δ = -1/2, h = 2
Начальный оператор: O_0 = Q_3 = ∑_i (s^+_i s^+_{i+1}s^+_{i+2} + s^-_i s^-_{i+1}s^-_{i+2})

Известно существование динамической симметрии с частотой ω = 12

2. Хаотическая цепь

H = ∑_i (s^x_i s^x_{i+1} - 1.05s^z_i + 0.5s^x_i)
Начальный оператор: O_0 = ∑_i (1.05s^x_i s^x_{i+1} + s^z_i)

Выбран оператор без перекрытия с H для исключения тривиальных сохраняющихся величин

Численная реализация

• Инструменты: Пакет Julia PauliStrings.jl (эффективные вычисления вложенных коммутаторов с использованием представления строк Паули)
• Трансляционная инвариантность: Хранится только элементарная ячейка, а не экстенсивные операторы
• Усиленное усечение Паули: На каждом шаге алгоритма Ланцоша сохраняются M строк Паули с наибольшим весом (алгоритм 1)

Методы сравнения

1. Моделирование матричного произведения операторов (MPO): В качестве эталона
2. Стандартные граничные условия Дирихле: Демонстрация нефизического отражения
3. Метод Теретентькова и др.: Линейная экстраполяция + искусственная диссипация (дополнительные материалы)

Метрики оценки

• Точность частоты динамической симметрии: Сравнение с известными результатами (XXZ) или аналитическими предсказаниями (хаотические)
• Верность временной эволюции: TrO_0(0)O_0(t) vs точные/MPO результаты
• Мера локальности: Распределение собственного вектора в базисе Крылова |ψ_0|² (более локальное — более тёмное)

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

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

1. Верификация идеальных случаев (рис. 2-3)

Линейный рост b_n = n:

• Динамика при открытых граничных условиях неотличима от бесконечной цепи
• Все собственные значения удовлетворяют Im(ω) = -2i (кроме двух тривиальных корней ω = -i, -3i)
• Граничные условия Дирихле приводят к значительному нефизическому отражению при t≈5

Рост по квадратному корню b_n = √n:

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

2. Результаты для XXZ-цепи (рис. 4a,c)

• Восстановление динамической симметрии: Успешное извлечение известной динамической симметрии вблизи ω = 12
  • L=26: Re(ω) ≈ 12.0, Im(ω) ≈ -0.5
  • Усиленное усечение Паули (L=50): значительное улучшение сходимости
• Временная эволюция: TrO_0(0)O_0(t) демонстрирует устойчивые осцилляции (совпадает с MPO до t≈3)
• Локальность: Только вещественнозначная динамическая симметрия является локальной (|ψ_0|² сосредоточена на малых n)

3. Результаты для хаотической цепи (рис. 4b,d)

• Отсутствие вечных симметрий: Все моды затухают, нет чисто вещественных частот
• Скопление скоростей затухания: Собственные значения сосредоточены вблизи Im(ω) ≈ -0.72i
• Согласие с теоретическим предсказанием: -0.72 ≈ -2λ, где λ — скорость роста коэффициентов Ланцоша
• Экспоненциальное затухание: TrO_0(0)O_0(t) быстро затухает к нулю (совпадает с MPO)
• Нелокальность: Все моды крайне нелокальны (|ψ_0|² равномерно распределена)

4. Закалка из стационарного состояния (рис. 5)

Закалка из |ψ_0⟩ = 2^{-N/2}(|0⟩+|1⟩)^⊗N:

• Q_1 = ∑_i (s^+_i + s^-_i): быстрое затухание
• Q_5 (3-локальная динамическая симметрия 5-го порядка): долгосрочные осцилляции
• Хорошая сходимость результатов для различных L (24-30)
• Стандартное усечение (серая линия) демонстрирует нефизическое поведение

Абляционные исследования

Влияние усечения Паули (дополнительный рис. 4)

• Без усечения Паули: L_max = 26 (предел точных вычислений)
• С усечением Паули (сохранение 2^22 строк): L_max = 50
  • Частота динамической симметрии: Re(ω) улучшается с ≈12.1 до ≈12.05
  • Im(ω) улучшается с ≈-0.8 до ≈-0.2

Анализ сходимости (дополнительные рис. 2-3)

• Модель XXZ: Re(ω) стабилизируется после L=10, Im(ω) продолжает улучшаться
• Хаотическая цепь: Спектр в основном сходится после L=30, собственные значения при большом L плотно сосредоточены вблизи линии Im(ω)=-2λ

Сравнение с методом Теретентькова (дополнительные рис. 6-7)

Хаотическая цепь:

• Собственные значения метода Теретентькова (γ=0.2-0.4) не сходятся к аналитическому результату
• Данный метод точно воспроизводит скопление Im(ω)=-2λ

XXZ-цепь:

• Обе методы имеют сходную точность для локальной динамической симметрии (ω=12)
• Предсказания для нелокальных мод значительно отличаются

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

1. Количественная связь между ростом операторов и термализацией: В случае линейного роста (хаотические) скорость затухания точно равна -2α, что прямо верифицирует теоретическое предсказание
2. Критерий локальности: |ψ_0|² в представлении Крылова автоматически обеспечивает количественный показатель локальности динамической симметрии
3. Применимость метода:
   • Хаотические системы (линейный рост): отличные результаты, предположение гладкости хорошо выполняется
   • Интегрируемые системы (рост по квадратному корню): возможно извлечение динамических симметрий, но сходимость медленнее
4. Вычислительные узкие места: Основное ограничение — вычисление коэффициентов Ланцоша высокого порядка (операторы становятся крайне нелокальными), усечение Паули обеспечивает эффективное облегчение

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

Методы пространства Крылова

• Сложность операторов: Parker и др. 105 для исследования роста операторов в квантовом хаосе
• Гидродинамика: Refs 97, 117-119 для исследования транспортных свойств
• Системы Флоке: Refs 120-124 применение к периодически управляемым системам

Разложение по псевдомодам

• Проекция Мори: Ранний метод вычисления временных корреляционных функций 92
• Недавние разработки: Refs 96, 97 предлагают добавление диссипации + экстраполяцию или исследование резонансов Рюэля-Поллико 98-104
• Ограничения: Ограничено динамикой вблизи равновесия, не рассматривает псевдолокальность

Теория динамических симметрий

• Вечное равновесие: Теоретическая структура Buča 82, форма состояния ρ(t) = Z^{-1}exp(∑_u μ_u e^{iλ_u t}A_u)
• Спектрогенерирующие алгебры: Refs 33, 34, 90, 91 связанные, но не обязательно псевдолокальные
• Временные кристаллы и рубцы: Refs 12, 29-44 приложения динамических симметрий

Преимущества данной работы

1. Без параметров: Не вводит искусственные скорости диссипации и другие дополнительные параметры
2. Без экстраполяции: Прямое использование известных коэффициентов Ланцоша
3. Псевдолокальность: Явно идентифицирует релевантные псевдолокальные симметрии
4. Вдали от равновесия: Применимо к сильным закалкам и другим неравновесным ситуациям

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

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

1. Эффективность метода: Открытые граничные условия успешно реализуют численное вычисление спектра Лиувиллиана бесконечных закрытых систем
2. Теоретические инсайты:
   • Хаотические системы: линейный рост Ланцоша → все наблюдаемые затухают со скоростью -2α
   • Интегрируемые системы: рост по квадратному корню → существуют вечные осциллирующие моды
3. Практическая ценность: В сочетании с методом строк Паули обеспечивает мощный инструмент для вычисления многотельной динамики

Ограничения

1. Предположение гладкости: Метод зависит от гладкости φ_n в представлении Крылова
   • Хаотические системы: линейный рост гарантирует гладкость
   • Другие случаи: требуется верификация
2. Вычислительная сложность:
   • Операторы высокого порядка O_n становятся крайне нелокальными
   • Модель XXZ может быть точно вычислена только до L=26
   • Усечение Паули вводит приближение, но эффективно
3. Точность на долгих временах:
   • Мнимая часть динамических симметрий не полностью сходится, приводя к ошибкам на долгих временах
   • Требуется больший L или улучшенная стратегия усечения
4. Трудности закалки из тепловых состояний: Закалка из теплового состояния ρ = e^{-βO_0} требует вычисления экспоненты суммы некоммутирующих операторов, численно сложно (реализована только закалка из стационарного состояния)

Будущие направления

1. Квантовые схемы: Обобщение на случаи, когда Лиувиллиан имеет верхнюю форму Хессенберга 107
2. Другие псевдолокальные симметрии:
   • Полулокальные динамические симметрии 137-139
   • Псевдолокальные симметрии квантовых рубцов 30, 31, 36, 82
3. Моделирование ЯМР экспериментов: Прямое вычисление преобразования Фурье Tr(Z_tot(t)Z_tot(0)) 140
4. Улучшение сходимости:
   • Более тонкие стратегии усечения
   • Адаптивное усечение Паули
   • Оптимизация рекурсивного алгоритма (дополнительные материалы)

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

Сильные стороны

1. Методологическая инновативность ⭐⭐⭐⭐⭐

• Оригинальные граничные условия: Основаны на естественном физическом предположении о гладкости волновой функции, отличаются от искусственного добавления диссипации
• Разделение система-окружение: Хитрая структура для разделения локального и нелокального на уровне пространства операторов
• Синтез теории и численности: Аналитические решения для идеальных случаев (полиномы Мейксера/Эрмита) дополняют численную верификацию на реальных системах

2. Теоретическая глубина ⭐⭐⭐⭐⭐

• Связь роста операторов и термализации: Впервые количественно доказана универсальная связь линейный рост → затухание -2α
• Количественная локальность: Через |ψ_0|² автоматически идентифицируются релевантные динамические симметрии
• Интеграция с теорией: Глубокая интеграция с теоретической структурой Buča 82

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

• Множественные модели: Интегрируемая (XXZ) и хаотическая цепь как типичные примеры
• Множественные сравнения: MPO, граничные условия Дирихле, метод Теретентькова
• Анализ сходимости: Детальное исследование зависимости от L и усечения Паули
• Недостаток: Отсутствие тестирования на большем количестве реальных физических систем (например, модель Хаббарда)

4. Техническая реализация ⭐⭐⭐⭐⭐

• Открытый исходный код: Реализация на Julia (GitHub + PauliStrings.jl)
• Оптимизация алгоритма: Трансляционная инвариантность, усечение Паули, рекурсивный алгоритм и другие практические приёмы
• Воспроизводимость: Детальные дополнительные материалы и параметры установки

Недостатки

1. Строгость предположения гладкости

• Отсутствие теоретического доказательства: Когда φ_n достаточно гладкая для линейной аппроксимации?
• Случаи отказа: Может не работать при сильной фрагментации или специальных граничных условиях
• Рекомендация: Требуется более строгое математическое обоснование или критерии применимости

2. Масштабируемость вычислений

• Верхний предел L: XXZ достигает только 26 (50 требует усечения), ограничивает высокоточные приложения
• Зависимость от размера: Не обсуждается зависимость от N (размер системы)
• Обобщение на 2D: Только 1D цепи, высшие размерности вычислительно сложны

3. Ограниченная экспериментальная верификация

• Типы закалок: Только закалка из стационарного состояния, не реализована закалка из теплового состояния
• Временные масштабы: Ограниченная точность на долгих временах (t<5 для XXZ)
• Физические величины: В основном автокорреляции операторов, отсутствуют другие наблюдаемые

4. Недостаточное сравнение с методом Теретентькова

• Только дополнительные материалы: Основной текст не содержит детального обсуждения
• Чувствительность к параметрам: Не систематически исследована зависимость от γ и схемы экстраполяции
• Анализ преимуществ: Когда данный метод превосходит? Когда эквивалентен?

Влияние на область

Вклад в область ⭐⭐⭐⭐⭐

• Парадигмальный сдвиг: Введение идей открытых систем в исследование динамики закрытых систем
• Универсальный инструмент: Применим к широкому спектру квантовых многотельных задач
• Теоретическая основа: Единая структура для понимания термализации, временных кристаллов, рубцов и других явлений

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

• Численный метод: Может быть прямо применён к вычислению динамики реальных систем
• Теоретическое руководство: Поддержка экспериментов ЯМР и других
• Ограничение: Вычислительные затраты всё ещё ограничивают приложения к сверхбольшим системам

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

• Открытый исходный код: Полная реализация на GitHub
• Подробная документация: Все детали алгоритмов в дополнительных материалах
• Прозрачность параметров: Все экспериментальные установки явно указаны

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

Идеальные сценарии

1. 1D системы с локальными взаимодействиями: Спиновые цепи, решёточные бозоны и т.д.
2. Хаотические системы: Линейный рост Ланцоша гарантирует эффективность метода
3. Краткосредняя динамика: t<10τ (характерная временная шкала)
4. Локальные наблюдаемые: Измерение k-локальных операторов

Сценарии, требующие осторожности

1. Сильно фрагментированные системы: Предположение гладкости может не выполняться
2. Дальнодействующие взаимодействия: Структура пространства Крылова может отличаться
3. Очень долгие времена: Неполная сходимость мнимой части приводит к накоплению ошибок
4. Высшие размерности: Вычислительные затраты резко возрастают

Неприменимые сценарии

1. Открытые системы: Существует более подходящая структура Линдблада
2. Системы с сильными измерениями: Требуются другие теоретические инструменты
3. Классический предел: Чрезмерно сложно, классические методы проще

Дополнительные технические детали

Ключевые формулы и выводы

Главное уравнение граничных условий (дополнительные материалы S20)

Точная временная эволюция Ȯ = iH,O, жёсткие граничные условия соответствуют:

Ȯ = iH,O - 1/(2N) Tr(O_L† iH,O)O_L

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

Ȯ = iH,O - 1/(2N) Tr(O_L† iH,O)O_L + ib_{L+1}Tr(O_† O)O_L - 2b_{L+1}Tr(O_L† O)O_L

Демонстрирует, как граничные условия балансируют поток операторов в окружение

Решение через полиномы Мейксера (дополнительные материалы S3-S7)

При линейном росте b_n = n собственные векторы:

φ_n(ω) = i^n M_n(ω)/n!

где M_n удовлетворяет рекуррентности: M_{n+1}(x) = xM_n(x) - n²M_(x)

Граничные условия дают ограничение на частоты: (ω + 2i(L+1))M_L(ω) = (2L+1)LM_(ω)

Тривиальные корни ω = -i, -3i соответствуют φ_n ∝ 1 и φ_n ∝ 2n-1

Численные приёмы

Рекурсивный алгоритм (дополнительные материалы)

Улучшенное извлечение динамических симметрий:

1. Случайная инициализация k-локального O_0
2. Диагонализация L̃, выбор собственной моды с частотой ≈ω, обозначаемой A_ω
3. Усечение A_ω до k-локального, обновление O_0 ← truncate(A_ω)
4. Повтор до сходимости

Для модели XXZ с k=3 сходимость за 4 шага: Re(ω) от 12.14 → 12.02

Избранные ссылки

Методологические основы

• 82 B. Buča, PRX 13, 031013 (2023) - теория вечного равновесия
• 105 D. Parker et al., PRX 9, 041017 (2019) - сложность Крылова
• 96,97 A. Teretenkov, O. Lychkovskiy et al., PRB (2024) - разложение по псевдомодам

Области приложения

• 12 B. Buča et al., Nat. Commun. 10, 1730 (2019) - временные кристаллы
• 29-44 Серия работ по квантовым рубцам
• 119 N. Loizeau et al., SciPost (2025) - PauliStrings.jl

Теоретический фундамент

• 1 L. D'Alessio et al., Adv. Phys. 65, 239 (2016) - обзор термализации
• 86-89 Теория псевдолокальности (Prosen, Doyon и др.)

Общая оценка: Это отличная статья с важными инновациями в области квантовой многотельной динамики (оценка 4.5/5). Основной метод (открытая цепь Крылова) элегантен и имеет ясный физический смысл, теоретический вывод строг (аналитические решения для идеальных случаев), численная верификация полна. Главные вклады заключаются в установлении количественной связи между ростом операторов и термализацией, а также в предоставлении практического численного инструмента. Основные ограничения касаются области применимости предположения гладкости и масштабируемости вычислений. Для исследователей, работающих над неравновесной динамикой квантовых многотельных систем, временными кристаллами, квантовыми рубцами и связанными явлениями, эта статья является обязательной к прочтению.