2025-11-18T17:19:13.403952

Digital quantum simulation of many-body localization crossover in a disordered kicked Ising model

Hayata, Seki, Yunoki
Simulating nonequilibrium dynamics of quantum many-body systems is one of the most promising applications of quantum computers. However, a faithful digital quantum simulation of the Hamiltonian evolution is very challenging in the present noisy quantum devices. Instead, nonequilibrium dynamics under the Floquet evolution realized by the Trotter decomposition of the Hamiltonian evolution with a large Trotter step size is considered to be a suitable problem for simulating in the present or near-term quantum devices. In this work, we propose simulating the many-body localization crossover as such a nonequilibrium problem in the disordered Floquet many-body systems. As a demonstration, we simulate the many-body localization crossover in a disordered kicked Ising model on a heavy-hex lattice using $60$ qubits from $156$ qubits available in the IBM Heron r2 superconducting qubit device named ibm\_fez. We compute out-of-time-ordered correlators as an indicator of the many-body localization crossover. From the late-time behavior of out-of-time-ordered correlators, we locate the quantum chaotic and many-body localized regimes as a function of the disorder strength. The validity of the results is confirmed by comparing two independent error mitigation methods, that is, the operator renormalization method and zero-noise extrapolation.
academic

Цифровое квантовое моделирование кроссовера многотельной локализации в неупорядоченной периодически возбуждаемой модели Изинга

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

  • ID статьи: 2510.01983
  • Название: Digital quantum simulation of many-body localization crossover in a disordered kicked Ising model
  • Авторы: Tomoya Hayata, Kazuhiro Seki, Seiji Yunoki
  • Классификация: quant-ph cond-mat.stat-mech hep-lat hep-th
  • Дата публикации: 15 октября 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2510.01983

Аннотация

В данной работе предложен метод моделирования явления кроссовера многотельной локализации (Many-body localization, MBL) на современных квантовых устройствах. Исследователи использовали 60 квантовых битов на сверхпроводящем квантовом процессоре IBM Heron r2 ibm_fez для моделирования кроссовера MBL в неупорядоченной периодически возбуждаемой модели Изинга на тяжелой гексагональной решётке. Путём расчёта временных неупорядоченных корреляционных функций (Out-of-time-ordered correlators, OTOCs) в качестве индикатора кроссовера MBL авторы выявили области квантового хаоса и многотельной локализации из поведения на поздних временах. Результаты верифицированы двумя независимыми методами смягчения ошибок: переормализацией операторов и экстраполяцией с нулевым шумом.

Научный контекст и мотивация

Центральные проблемы

  1. Моделирование неравновесной динамики квантовых многотельных систем: это одно из наиболее перспективных приложений квантовых компьютеров, однако точное моделирование эволюции гамильтониана на современных шумных квантовых устройствах представляет значительные трудности.
  2. Явление многотельной локализации: MBL — это механизм, препятствующий термализации неупорядоченных квантовых многотельных систем, имеющий важное значение для понимания нарушения эргодичности в квантовых системах.
  3. Системы Флоке: периодически возбуждаемые квантовые системы могут демонстрировать экзотические фазы, недостижимые в равновесных системах, такие как временные кристаллы.

Научная мотивация

  • Технологическая мотивация: шум в современных квантовых устройствах ограничивает точное моделирование эволюции гамильтониана, тогда как эволюция Флоке, реализуемая через разложение Троттера, лучше подходит для современных квантовых устройств
  • Научная мотивация: переход MBL включает неравновесную динамику и полный спектр собственных состояний, что затруднительно для классических компьютеров при работе с большими системами; квантовое моделирование имеет естественное преимущество
  • Практическая мотивация: исследование новых направлений для приложений современных квантовых устройств до реализации отказоустойчивых квантовых компьютеров

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

  1. Первая реализация цифрового квантового моделирования кроссовера MBL на 60-кубитной системе, превосходящей возможности классических методов точной диагонализации
  2. Предложен метод идентификации кроссовера MBL на основе OTOCs, локализирующий области квантового хаоса и многотельной локализации через поведение на поздних временах
  3. Разработана неупорядоченная периодически возбуждаемая модель Изинга для тяжёлой гексагональной решётки, полностью использующая встроенную связность и операции вентилей квантовых устройств IBM
  4. Установлены два независимых метода смягчения ошибок: переормализация операторов и экстраполяция с нулевым шумом для взаимной верификации надёжности результатов
  5. Доказана осуществимость моделирования неравновесной динамики многотельных систем Флоке на современных квантовых устройствах

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

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

Исследование явления кроссовера многотельной локализации в неупорядоченной периодически возбуждаемой модели Изинга на тяжёлой гексагональной решётке путём варьирования интенсивности беспорядка W и идентификации границы между областью квантового хаоса (W < Wc) и областью многотельной локализации (W > Wc).

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

Неупорядоченная периодически возбуждаемая модель Изинга

Зависящий от времени гамильтониан определяется как:

H_X = \sum_{i=1}^N B_X(i)X_i, & t \in [0,T/2) \\ H_Z = J\sum_{(i,j)} Z_iZ_j + B_Z\sum_{i=1}^N Z_i, & t \in [T/2,T) \end{cases}$$ где: - $X_i$, $Z_i$ — операторы Паули для i-го узла решётки - $B_X(i)$ — случайное поле, равномерно распределённое на $[B_{X0}-W, B_{X0}+W]$ - $J$, $B_Z$, $B_{X0}$ — фиксированные параметры - $W$ — параметр интенсивности беспорядка #### Оператор эволюции Флоке Оператор эволюции за один период возбуждения: $$U_F = e^{-iH_ZT/2}e^{-iH_XT/2}$$ Параметры: $JT = π/2$, $B_ZT = 1.3$, $B_{X0}T = π/2$, обеспечивающие максимальное хаотическое состояние системы при $W=0$. #### Временные неупорядоченные корреляционные функции (OTOCs) Определяются как: $$\text{OTOC} = \langle 0^N|(U_F)^{\dagger n}X_1(U_F)^nZ_m(U_F)^{\dagger n}X_1(U_F)^n|0^N\rangle$$ Нормализованная форма: $$\text{OTOC} = \frac{\langle 0^N|(U_F)^{\dagger n}X_1(U_F)^nZ_m(U_F)^{\dagger n}X_1(U_F)^n|0^N\rangle}{\langle 0^N|(U_F)^{\dagger n}I_1(U_F)^nZ_m(U_F)^{\dagger n}I_1(U_F)^n|0^N\rangle}$$ ### Технические инновации 1. **Оптимизация схем**: ручное удаление избыточных вентилей вне причинного конуса, значительно улучшающее качество смягчения ошибок 2. **Эффективный квантовый объём**: использование знаменателя для оценки роста причинного конуса как дополнительного индикатора MBL 3. **Двойное смягчение ошибок**: комбинирование переормализации операторов и метода ZNE для взаимной верификации 4. **Анализ в зависимости от расстояния**: идентификация характеристик распространения волнового фронта через поведение OTOC на различных расстояниях x ## Экспериментальная установка ### Квантовое устройство - **Устройство**: сверхпроводящий квантовый процессор IBM Heron r2 ibm_fez (156 квантовых битов) - **Используемый масштаб**: 60 квантовых битов, образующих тяжёлую гексагональную решётку 2×3 - **Частоты ошибок вентилей**: вентиль CZ ~3.7×10⁻³, вентиль SX ~2.93×10⁻⁴, ошибка считывания ~1.76×10⁻² ### Экспериментальные параметры - **Количество реализаций беспорядка**: 25 - **Количество шагов Троттера**: 10 - **Общее количество схем**: 250 (25×10) - **Количество измерений на схему**: 16 000 - **Количество схем обмотки Паули**: 32 ### Метрики оценки 1. **Нормализованный OTOC**: основной индикатор, отражающий распространение возмущения информации 2. **Эффективный квантовый объём**: $V_{eff} = \log F / \log(1-p)$, характеризующий размер причинного конуса 3. **Локализация точки кроссовера**: $W_c = \arg\max_W \partial\langle\langle\text{OTOC}\rangle\rangle/\partial W$ ### Методы смягчения ошибок 1. **Переормализация операторов**: использование знаменателя для коррекции влияния шума 2. **Экстраполяция с нулевым шумом (ZNE)**: экстраполяция через коэффициенты усиления шума f=1.0, 1.5 ## Экспериментальные результаты ### Основные результаты #### Характеристики поведения OTOC - **Область слабого беспорядка (W≤0.18)**: OTOC резко снижается на расстоянии x=n, сильно затухает на поздних временах, проявляя поведение квантового хаоса - **Область сильного беспорядка (W≥0.18)**: OTOC остаётся близким к 1 или медленно затухает, демонстрируя характеристики многотельной локализации - **Точка кроссовера**: $W_c \sim 0.18$, результаты обоих методов смягчения ошибок согласуются #### Анализ распространения волнового фронта - **W=0.05**: чёткий волновой фронт, распространяющийся вблизи x=n - **W=0.5**: волновой фронт исчезает, распространение информации ограничено #### Эффективный квантовый объём - **Слабый беспорядок**: сверхлинейный рост на поздних временах, отражающий быстрое расширение причинного конуса - **Сильный беспорядок**: линейный рост, указывающий на прекращение расширения причинного конуса ### Верификация смягчения ошибок Два метода смягчения ошибок дают согласованные результаты в областях хаоса и MBL; область кроссовера (0.14≤W≤0.2) требует осторожного анализа. Метод ZNE делает поведение фазового перехода более чётким, но имеет относительно большие полосы ошибок. ### Экспериментальные находки 1. **Зависимость от расстояния**: OTOC при x=n/2=5 наиболее подходит в качестве индикатора кроссовера 2. **Эффекты размера системы**: 60-кубитная система демонстрирует масштаб, недостижимый для классических вычислений 3. **Робастность к шуму**: нормализованный OTOC обладает определённой устойчивостью к шуму ## Связанные работы ### Теоретические исследования MBL - Теоретическая база, основанная на нарушении гипотезы термализации собственных состояний - Логарифмическое поведение светового конуса в одномерных системах - Дискуссии о стабильности MBL в многомерных системах ### Экспериментальные работы по квантовому моделированию - Реализация временных кристаллов на сверхпроводящих кубитах - Идентификация локальных интегралов движения - Применение вариационного квантового решателя для возбуждённых состояний ### Системы Флоке - Экзотические фазы периодически возбуждаемых систем - Эквивалентность динамики Троттера и систем Флоке - Квантовое моделирование явлений предварительного нагрева ## Заключение и обсуждение ### Основные выводы 1. **Успешная реализация моделирования кроссовера MBL на 60 кубитах**, локализирующая точку кроссовера $W_c \sim 0.18$ 2. **Верификация эффективности OTOC в качестве индикатора MBL**, различающего области хаоса и локализации 3. **Доказательство практических возможностей современных квантовых устройств**, открывающее новые направления для приложений до эры отказоустойчивых компьютеров ### Ограничения 1. **Ограничение количества шагов Троттера**: шум ограничивает наблюдение асимптотического поведения на длительных временах 2. **Единственный размер системы**: отсутствие анализа конечномерного масштабирования затрудняет определение наличия истинного фазового перехода 3. **Стабильность двумерного MBL**: теоретически стабильность фазы MBL в многомерных системах остаётся спорной 4. **Ограничение причинного конуса**: невозможность расширения схемы при покрытии причинным конусом всех кубитов ### Будущие направления 1. **Устройства с меньшим шумом**: исследование логарифмического поведения светового конуса на более длительных временах 2. **Анализ конечномерного масштабирования**: исследование систем различных размеров для определения природы фазового перехода 3. **Другие неравновесные фазы**: расширение на квазикристаллические временные структуры и другие экзотические фазы 4. **Классически неразрешимые задачи**: поиск дополнительных приложений с квантовым преимуществом ## Глубокая оценка ### Преимущества 1. **Технологическая инновативность**: первая реализация моделирования MBL на крупномасштабном квантовом устройстве; методы оптимизации схем и смягчения ошибок обладают универсальностью 2. **Экспериментальная полнота**: 25 реализаций беспорядка, двойное смягчение ошибок, анализ на множественных расстояниях обеспечивают надёжность результатов 3. **Физическое значение**: успешная идентификация кроссовера MBL предоставляет новую экспериментальную платформу для квантовой многотельной физики 4. **Методологический вклад**: установление парадигмы исследования неравновесной динамики на шумных квантовых устройствах ### Недостатки 1. **Теоретическая глубина**: отсутствие теоретических предсказаний и сравнения численных значений точки кроссовера 2. **Анализ масштабирования**: единственный размер системы ограничивает понимание природы фазового перехода 3. **Модель шума**: предположение о деполяризирующем канале может быть чрезмерно упрощённым 4. **Универсальность**: применимость методов к другим квантовым устройствам и моделям требует верификации ### Влияние 1. **Академическая ценность**: предоставление важного экспериментального эталона для пересечения квантовой многотельной физики и квантовых вычислений 2. **Технологический прогресс**: методы смягчения ошибок и оптимизации схем могут быть распространены на другие задачи квантового моделирования 3. **Перспективы приложений**: доказательство потенциала современных квантовых устройств для научных вычислений 4. **Воспроизводимость**: детальные экспериментальные параметры и открытые данные поддерживают повторение результатов ### Применимые области 1. **Исследования квантовой многотельной физики**: неравновесная динамика, фазовые переходы, критические явления 2. **Разработка квантовых алгоритмов**: методы смягчения ошибок и оптимизации схем 3. **Оценка квантовых устройств**: тестирование производительности и характеризация 4. **Образование и обучение**: примеры для преподавания квантового моделирования и многотельной физики ## Библиография Данная работа цитирует 48 важных источников, охватывающих: - Теоретические основы MBL [10-17, 25, 29-33] - Экспериментальные работы по квантовому моделированию [19, 26-28, 35] - Системы Флоке [1-6, 15-17] - Методы смягчения ошибок [22, 34, 35, 44] - Теорию квантового хаоса [36-39] Ключевые источники включают обзор MBL Nandkishore & Huse, доказательство квантового преимущества командой IBM и теоретические основы OTOC. --- **Общая оценка**: Это важная экспериментальная работа на пересечении квантовой многотельной физики и квантовых вычислений. Авторы успешно реализовали цифровое моделирование кроссовера MBL на крупномасштабном квантовом устройстве с инновационными техническими методами и надёжными результатами. Несмотря на некоторые ограничения в теоретической глубине и размере системы, работа открывает новые направления для научных приложений современных квантовых устройств и имеет значительную академическую ценность и технологическое значение.