Плавление доменных стенок и запутанность в цепочках свободных фермионов с зонной структурой

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

• ID статьи: 2507.03780
• Название: Domain-wall melting and entanglement in free-fermion chains with a band structure
• Автор: Виктор Айслер (Технологический университет Граца и Университет Граца)
• Классификация: cond-mat.stat-mech, cond-mat.quant-gas, quant-ph
• Дата публикации: 18 сентября 2025 г. (подача в SciPost Physics)
• Ссылка на статью: https://arxiv.org/abs/2507.03780

Аннотация

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

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

Проблемный контекст

1. Проблема плавления доменных стенок: Плавление доменных стенок является классической установкой для изучения динамики одномерных квантовых цепочек вдали от равновесия путём подготовки ступенчатого начального распределения плотности и наблюдения его временной эволюции для понимания свойств переноса.
2. Ограничения существующих исследований:
   • Большинство исследований сосредоточены на трансляционно-инвариантных однородных цепочках
   • Недостаточно изучены случаи, когда сам оператор временной эволюции нарушает трансляционную инвариантность
   • Отсутствует систематическое понимание распространения запутанности в системах с зонной структурой
3. Исследовательская мотивация:
   • Исследовать влияние периодических амплитуд прыжков на динамику плавления доменных стенок
   • Понять закономерности роста энтропии запутанности в многозонных системах
   • Установить связь между системами с зонной структурой и однородными системами

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

1. Построение теоретической базы: Первое систематическое исследование проблемы плавления доменных стенок в цепочках свободных фермионов с p-точечной трансляционной инвариантностью
2. Закономерности роста запутанности: Доказано, что энтропия запутанности сохраняет логарифмический рост, но предфактор пропорционален числу заполненных квазичастичных зон p
3. Гидродинамический предел: Для димеризованной цепочки (p=2) показано, что распределения плотности и электрического тока могут быть получены из однородного случая простой переномировкой скоростей
4. Разделение универсальных и неуниверсальных вкладов:
   • Универсальный вклад может быть получен аналитически и удвоен по сравнению с однородным случаем
   • Неуниверсальный вклад требует численного расчёта через блочные матрицы Тёплица
5. Анализ осцилляционных явлений: Обнаружены и объяснены осцилляции в росте энтропии, частота которых определяется величиной энергетической щели

Методологические подробности

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

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

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

Гамильтониан димеризованной цепочки с прыжками

Ĥ = -∑_{m=-N/2+1}^{N/2} [(1+δ)/2 c†_{2m-1}c_{2m} + (1-δ)/2 c†_{2m}c_{2m+1} + h.c.]

Процесс диагонализации

1. Преобразование подрешётки: Введение фермионных операторов двух подрешёток a_m = c_{2m-1}, b_m = c_{2m}
2. Преобразование Фурье:

   (a_m, b_m)^T = (1/√N) ∑_k e^{ikm} (a_k, b_k)^T
   

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

   Ĥ = -∑_k ω_k(α†_k α_k - β†_k β_k)
   

   
   где дисперсионное соотношение имеет вид:

   ω_k = √[cos²(k/2) + δ²sin²(k/2)]
   

Расчёт корреляционных функций

Гидродинамический предел

В пределе m,n,t → ∞ при фиксированных m/t, n/t через метод стационарной фазы получено:

C_{m,n}(t) = ∫_{q_-}^{q_+} dq/(2π) e^{iq(2n-2m)} [(1, e^{iφ_q}); (e^{-iφ_q}, 1)]
           + ∫_{-q_+}^{-q_-} dq/(2π) e^{iq(2n-2m)} [(1, -e^{iφ_q}); (-e^{-iφ_q}, 1)]

где точки Ферми q_±(m,n,t) определяются условием v_{q_±} = (m+n-1)/t.

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

1. Представление блочными матрицами: Использование двухточечной трансляционной инвариантности для описания корреляционных функций матрицами 2×2
2. Образ квазичастиц: Два интегральных члена соответствуют независимым вкладам квазичастиц верхней и нижней зон
3. Комбинирование КТП и ОГД: Объединение методов конформной теории поля с обобщённой гидродинамикой для обработки распространения запутанности

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

Параметры численных расчётов

• Длина цепочки: N = 300 (открытые граничные условия для избежания удвоения в передней области)
• Временная эволюция: t = 200, 300
• Параметры димеризации: δ = 0.02, 0.05, 0.1, 0.2, 0.3, 0.5

Подготовка начального состояния

Начальное состояние доменной стенки: левая половина цепочки полностью заполнена, правая половина пуста:

C(0) = [1 0; 0 0]

Методы оценки

1. Плотность и электрический ток: Сравнение с аналитическими результатами гидродинамики
2. Энтропия запутанности: Расчёт через собственные значения матрицы корреляций
3. Поведение масштабирования: Проверка логарифмического закона роста

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

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

1. Распределение плотности и электрического тока

• Плотность: ρ(n,t) = (1/π)arccos(n/(v_max t)), где v_max = 1-δ
• Электрический ток: J(n,t) = (v_max/π)√1-(n/(v_max t))²
• Проверка: Численные результаты идеально совпадают с аналитическими выражениями

2. Рост энтропии запутанности

• Энтропия половины цепочки: S ~ (1/3)ln t + const (удвоена по сравнению с однородным случаем 1/6 ln t)
• Универсальный вклад: S_u = (1/3)ln2v_max t(1-ζ²)
• Неуниверсальный вклад: Требует численного расчёта через блочные матрицы Тёплица

3. Осцилляционные явления

Обнаружены затухающие осцилляции в росте энтропии:

S - (1/3)ln t ≈ α + β cos(2δt + φ)/(δt)

Частота осцилляций точно равна величине энергетической щели 2δ.

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

Предел δ → 1

В пределе полной димеризации аналитически получено:

lim_{δ→1} S_0(ζ) = (1/3)ln√(1-ζ²) + s_0

Трёхзонная модель

Для трёхточечной периодичности проверено:

• Энтропия половины цепочки: S = (1/2)ln t + const
• Распределение энтропии показывает несколько областей масштабирования, отражающих различия в скоростях разных зон

Анализ конкретных случаев

Поведение временного масштабирования

• t ≪ δ^{-1}: Поведение совпадает с однородной цепочкой
• t ≫ δ^{-1}: Переход к удвоенному логарифмическому росту
• Масштаб кроссовера определяется обратной величиной энергетической щели

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

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

1. Однородные цепочки: Первоначальное исследование плавления доменных стенок в цепочке XX проведено Antal и соавторами 10
2. Взаимодействующие системы: Расширение на взаимодействующие системы через метод ОГД 11
3. Исследования запутанности: Понимание распространения запутанности методами КТП 21,22

Неоднородные системы

1. Локальные дефекты: Одиночная примесь приводит к линейному росту энтропии 38-44
2. Периодические модуляции: Первое систематическое исследование влияния периодических амплитуд прыжков

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

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

1. Универсальное масштабирование: Энтропия запутанности сохраняет логарифмический рост, предфактор пропорционален числу зон
2. Переномировка скоростей: Гидродинамические величины могут быть получены простой переномировкой скоростей
3. Эффект энергетической щели: Создаёт характерные осцилляции, частота которых определяется щелью

Ограничения

1. Неуниверсальная часть: Требует численных расчётов, отсутствуют аналитические выражения
2. Структура краёв: Анализ тонкой структуры за пределами метода стационарной фазы затруднён
3. Расширение на взаимодействующие системы: Димеризация нарушает интегрируемость, метод ОГД неприменим

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

1. Произвольное заполнение: Исследование случаев с неполным заполнением
2. Краевая физика: Анализ универсального поведения в области масштабирования t^{1/3}
3. Взаимодействующие системы: Исследование существования эффективного гидродинамического описания

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

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

1. Теоретическая полнота: Полная теоретическая база от диагонализации гамильтониана до расчёта запутанности
2. Методологическая инновативность: Искусное объединение методов КТП, ОГД и матричного анализа
3. Универсальность результатов: Установление общих закономерностей для многозонных систем
4. Численная верификация: Аналитические результаты получили точную численную проверку

Недостатки

1. Технические ограничения: Отсутствуют аналитические результаты для постоянного члена блочной матрицы Тёплица
2. Физическая интерпретация: Глубокий физический механизм осцилляционных явлений требует дальнейшего изучения
3. Экспериментальная доступность: Не обсуждены схемы экспериментальной проверки теоретических результатов

Влияние

1. Теоретический вклад: Новая перспектива на динамику неоднородных квантовых систем
2. Методологическая ценность: Метод блочных матриц Тёплица может быть обобщён на другие системы
3. Перспективы применения: Теоретическое руководство для экспериментов по квантовому моделированию

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

1. Системы холодных атомов: Реализация периодической модуляции прыжков в оптических решётках
2. Физика твёрдого тела: Динамика экситонов в димеризованных материалах
3. Квантовая информация: Понимание распространения запутанности в многозонных системах

Список литературы

Статья цитирует 76 связанных работ, охватывающих:

• Теорию обобщённой гидродинамики 3-8
• Классические работы по плавлению доменных стенок 9-17
• Теорию энтропии запутанности 18-27
• Методы конформной теории поля 47-49
• Теорию матриц Тёплица 53-55,72-76