Моделирование гамильтониана с явными формулами для цифро-аналогового квантовых вычислений

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

• ID статьи: 2511.11404
• Название: Hamiltonian simulation with explicit formulas for Digital-Analog Quantum Computing
• Авторы: Mikel Garcia de Andoin (Университет Страны Басков), Thorge Müller (Немецкий центр аэрокосмических исследований), Gonzalo Camacho (Немецкий центр аэрокосмических исследований)
• Категория: quant-ph (квантовая физика), math-ph (математическая физика), math.MP (математическая физика)
• Дата публикации: 14 ноября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2511.11404v1

Аннотация

В данной работе предложен новый протокол моделирования гамильтониана для парадигмы цифро-аналогового квантового вычисления (Digital-Analog Quantum Computing, DAQC). DAQC использует естественный гамильтониан взаимодействия системы в качестве ресурса запутанности в сочетании с однокубитными вентилями для реализации универсальных квантовых операций. Традиционные методы проектирования оптимальных DAQC-цепей требуют экспоненциального времени классических вычислений. В статье предложено субоптимальное решение, которое сводит задачу к полиномиально разрешимой, а именно: путём собственного разложения матрицы связи размером 3N×3N (где N — число кубитов) генерируются эффективные DAQC-цепи за время O(N³) с количеством цифро-аналоговых блоков не более 12N².

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

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

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

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

• Фундаментальная потребность квантового моделирования: моделирование гамильтониана является одним из основных приложений квантовых вычислений с широкими перспективами применения в химии и физике конденсированного состояния
• Преимущества DAQC: по сравнению с чистыми цифровыми квантовыми вычислениями, DAQC объединяет универсальность цифровой парадигмы с устойчивостью к шумам аналоговых вычислений
• Узкое место масштабируемости: существующие методы требуют экспоненциальных вычислительных ресурсов на этапе предварительной обработки, что серьёзно ограничивает применение к крупномасштабным квантовым системам

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

• Экспоненциальная сложность: поиск оптимальной DAQC-цепи считается NP-трудной задачей, требующей экспоненциальных ресурсов
• Зависимость от численной оптимизации: существующие протоколы требуют численной оптимизации в большом пространстве параметров с высокой вычислительной стоимостью
• Ограничения гамильтониана источника: ранние методы накладывают строгие требования на топологическую структуру гамильтониана источника (требуется hᵢⱼˢ ≠ 0 тогда и только тогда, когда hᵢⱼᴾ ≠ 0)
• Ограничения MPS-прокси: использование матричного произведения состояний (MPS) для оптимизации эффективно только в системах, удовлетворяющих закону площади

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

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

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

1. Алгоритм построения за полиномиальное время: предложен аналитический метод построения DAQC-цепей за время O(N³), избегающий экспоненциальной численной оптимизации
2. Явные формулы разложения (Результат 1): доказано, что произвольный двухчастичный гамильтониан может быть представлен как сумма не более O(N²) локальных унитарных преобразований, действующих на гамильтониан Изинга типа ZZ: $TH_P = \sum_{q=1}^{O(N^2)} t_q U_q H_S U_q^\dagger, \quad t_q > 0$
3. Метод собственного разложения: задача преобразуется в собственное разложение положительно полуопределённой матрицы B размером 3N×3N с использованием стратегии «разделяй и властвуй» для построения эффективного разложения каждого собственного вектора
4. Анализ сложности цепи: генерируемая DAQC-цепь содержит не более 12N² цифро-аналоговых блоков, что находится на одном уровне с предыдущими методами, требующими 9N(N-1)/2 блоков
5. Границы времени моделирования: предоставлены верхние границы для общего времени моделирования: $t_A \leq \sum_k \lambda_k = 3N|\tilde{\lambda}_{min}|$

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

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

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

• Проблемный гамильтониан: $H_P = \sum_{i<j,\mu,\nu} h_{ij}^{\mu\nu} \sigma_i^\mu \sigma_j^\nu$ (произвольный двухчастичный гамильтониан)
• Гамильтониан источника: $H_S = \sum_{i<j} h_{ij}^{zz} \sigma_i^z \sigma_j^z$ (гамильтониан Изинга типа ZZ)
• Время эволюции: T

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

• Параметры DAQC-цепи: времена моделирования блоков {tₖ} и однокубитные вентили {Uₖ}
• Удовлетворяющие: $e^{-iTH_P} \approx \prod_q U_q e^{-it_q H_S} U_q^\dagger$

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

• tₖ > 0 (физическая реализуемость)
• однокубитные вентили должны удовлетворять условиям нормализации

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

1. Общая структура

DAQC-цепь состоит из чередующихся цифровых блоков (однокубитные вентили) и аналоговых блоков (свободная эволюция):

[Цифровой блок U₁] → [Аналоговая эволюция t₁Hₛ] → [Цифровой блок U₁†] → [Цифровой блок U₂] → [Аналоговая эволюция t₂Hₛ] → ...

2. Построение матрицы задачи

Коэффициенты связи гамильтониана организуются в матрицу B размером 3N×3N: $B_{3i+\mu,3j+\nu} = Th_{ij}^{\mu\nu}/h_{ij}^{zz}, \quad \forall i \neq j$

где диагональные блоки (i=j) имеют неопределённые элементы, которые можно свободно выбирать для удовлетворения положительной полуопределённости.

3. Обработка положительной полуопределённости

• Установить неопределённые диагональные элементы в нуль и вычислить минимальное собственное значение $\tilde{\lambda}_{min}$
• Установить диагональные блоки равными $-\tilde{\lambda}_{min}$, чтобы сделать B положительно полуопределённой матрицей
• Выполнить собственное разложение: $B = U^\dagger \lambda U$

4. Разложение собственных векторов (ключевое нововведение)

Для каждого собственного вектора $\vec{v}_k$ (соответствующего собственному значению λₖ) построить 2N пар векторов: $\vec{\gamma}_{+k}^{(\ell)}, \vec{\gamma}_{-k}^{(\ell)}, \quad \ell = 1, ..., 2N$

Каждый вектор с i-м трёхмерным блоком определяется как: $\gamma_{ik}^{(\ell)} = \frac{v_{ik} \pm \epsilon_{ik}^{(\ell)}}{\sqrt{\|v_{ik}\|^2 + \|\epsilon_{ik}^{(\ell)}\|^2}}$

где вектор возмущения: $\epsilon_{ik}^{(\ell)} = \cos\theta_{ik}^{(\ell)} \eta_{ik} + \sin\theta_{ik}^{(\ell)} \xi_{ik}$

параметры углов: $\theta_{ik}^{(\ell)} = \frac{\pi(i-1)(\ell-1)}{N}$

здесь $v_{ik} \perp \eta_{ik} \perp \xi_{ik}$, и $\|\eta_{ik}\|^2 = \|\xi_{ik}\|^2 = \max_i\|v_{ik}\|^2 - \|v_{ik}\|^2$

5. Полная формула разложения

$B = \sum_{k=1}^{3N} \lambda_k \vec{v}_k \vec{v}_k^\dagger = \sum_{k=1}^{3N} t_k \sum_{\ell=1}^{2N} \left(\vec{\gamma}_{+k}^{(\ell)}\vec{\gamma}_{+k}^{(\ell)\dagger} + \vec{\gamma}_{-k}^{(\ell)}\vec{\gamma}_{-k}^{(\ell)\dagger}\right)$

где время моделирования блока: $t_k = \frac{\lambda_k \max_i\|v_{ik}\|^2}{4N}$

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

1. Преобразование условий нормализации

Основная сложность: собственные векторы удовлетворяют глобальной нормализации $\|\vec{v}_k\| = 1$, но DAQC требует локальной нормализации $\|\gamma_{ik}\| = 1, \forall i$

Инновационное решение: путём добавления ортогональных возмущений $\epsilon_{ik}^{(\ell)}$ каждый трёхмерный блок нормализуется независимо, сохраняя при этом соотношение суммы внешних произведений

2. Стратегия «разделяй и властвуй»

Вместо прямого решения задачи оптимизации с 9NK переменными (K — количество блоков):

• Разложить B на 3N матриц ранга 1 (внешние произведения собственных векторов)
• Для каждой матрицы ранга 1 независимо построить разложение в 2N шагов
• Общая сложность снижается с экспоненциальной на полиномиальную

3. Аналитическое построение ортогональных возмущений

Путём параметризации векторов возмущений тригонометрическими функциями, используя условия ортогональности: $\sum_{\ell=1}^{2N} \epsilon_{ik}^{(\ell)}\epsilon_{jk}^{(\ell)\dagger} = 0, \quad \forall i \neq j$

это эквивалентно ортогональности дискретного преобразования Фурье и имеет аналитическое решение.

4. Различия с базовыми методами

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

Набор данных

Генерация случайных задач:

• Прямая генерация матрицы B размером 3N×3N (а не конкретного физического гамильтониана)
• Элементы выбираются из равномерного распределения U-1,1
• Нормализация: $\max|B_{3i+\mu,3j+\nu}| = 1$
• Масштаб системы: N = 1–50 кубитов
• Для каждого значения N генерируется 10⁴ случайных экземпляров

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

Общее время моделирования: $t_A = \sum_{q=1}^{12N^2} t_q$

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

• глубину цепи
• накопление ошибок декогеренции
• фактическое время выполнения

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

• Теоретическая верхняя граница44: $t_A \leq 3N|\tilde{\lambda}_{min}|$
• Сравнение с методами из литературы15-17, требующими 9N(N-1)/2 блоков по количеству блоков

Детали реализации

• Собственное разложение использует стандартные библиотеки численной линейной алгебры
• Для каждого собственного значения λₖ < ε (порог) соответствующие члены могут быть отброшены для уменьшения размера цепи
• Параметризация однокубитных вентилей: $R(\theta, \hat{n}) = e^{-i\frac{\theta}{2}(n_x\sigma^x + n_y\sigma^y + n_z\sigma^z)}$

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

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

1. Изменение времени моделирования в зависимости от масштаба системы

Как показано на рисунке 2:

• Сплошная линия: среднее значение tₐ для 10⁴ запусков
• Пунктирная линия: теоретическая верхняя граница 3N|λ̃ₘᵢₙ|
• Цветные области: диапазон между максимальным и минимальным значениями

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

• При условии нормализации $\max|B_{3i+\mu,3j+\nu}| = 1$ tₐ остаётся почти постоянным при увеличении N
• Теоретическая верхняя граница растёт линейно с N, но фактические значения значительно ниже верхней границы
• Дисперсия (ширина цветной области) немного увеличивается с N, но остаётся управляемой

2. Сравнение с теоретической верхней границей

• Скорость роста теоретической верхней границы: O(N)
• Скорость роста фактического среднего значения: ~O(1) (приблизительно постоянная)
• это показывает, что метод в практических приложениях значительно более эффективен, чем оценка наихудшего случая

3. Проверка масштабируемости

• Успешная проверка системы с N=50 (матрица размером 3×50=150)
• Время вычисления остаётся полиномиальным (O(N³))
• По сравнению с методами, требующими экспоненциальных ресурсов, демонстрирует значительное преимущество в масштабируемости

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

1. Линейная зависимость

$t_A \sim T \max|h_{ij}^{\mu\nu} / h_{ij}^{zz}|$

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

2. Влияние нормализации

Когда максимальное значение элементов B фиксировано, tₐ не растёт с масштабом системы, что соответствует интуиции:

• более крупные системы имеют больше степеней свободы
• распределение собственных значений более рассеянное
• средний вклад каждого собственного вектора уменьшается

3. Фактические требования к количеству блоков

Хотя теоретически требуется 12N² блоков, путём отбрасывания членов, соответствующих малым собственным значениям, фактически требуемое количество блоков может быть значительно уменьшено.

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

1. Оптимизация цифровых квантовых вычислительных цепей

• Универсальные наборы вентилей1,2: произвольные однокубитные вентили + двухкубитные вентили
• Оптимизация компиляции19-21: полиномиальные субоптимальные алгоритмы (например, разбиение матроида, AlphaTensor)
• Теория сложности18: точная оптимизация квантовых цепей является co-NQP-трудной

2. Парадигма DAQC

• Фундаментальная теория3,4: использование естественного гамильтониана + SQG для реализации универсальных квантовых вычислений
• Преимущества устойчивости к шумам5: лучшая устойчивость к шумам по сравнению с чистыми цифровыми методами
• Ранние протоколы14,15: требуют экспоненциальных ресурсов или эвристических алгоритмов

3. Временно-оптимальные многокубитные вентили

• Метод группы Клиффорда16: использование конечного набора SQG, ослабление требований к гамильтониану источника
• Оптимизация MPS-прокси17: использование произвольного SQG, но требует моделирования MPS, применимо только к системам с низкой запутанностью
• Исследования границ времени15: исследование сложности и временных границ многокубитных вентилей

4. Квантовое моделирование

• Теоретические основы22-24: универсальный квантовый симулятор Ллойда, разложение Троттера
• Экспериментальный прогресс25-30: эксперименты по квантовому моделированию с использованием ионных ловушек, ультрахолодных атомов, массивов атомов Ридберга
• Области применения37-39: квантовое моделирование в химии и физике конденсированного состояния

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

• Вычислительная эффективность: O(N³) против экспоненциальной
• Аналитичность: явные формулы против численной оптимизации
• Универсальность: применимо к произвольным двухчастичным гамильтонианам
• Реализуемость: совместимо с существующими квантовыми платформами (сверхпроводящие, ионные ловушки, нейтральные атомы)

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

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

1. Полиномиальная разрешимость: доказано существование полиномиального времени (O(N³)) субоптимального решения для задачи построения DAQC-цепи
2. Аналитический метод построения: предоставлена явная формула на основе собственного разложения матрицы, не требующая численной оптимизации
3. Практическая производительность: при типичном распределении задач общее время моделирования tₐ не растёт с масштабом системы
4. Масштабируемость: делает возможным проектирование DAQC-цепей для крупномасштабных квантовых систем (N≥50)

Ограничения

1. Ограничения гамильтониана источника

• Текущие требования: применимо только к гамильтониану Изинга типа ZZ
• Расширяемость: может быть обобщено на симметричные члены (XX, YY), но универсальный гамильтониан источника требует вложенного применения, приводящего к O(N⁴) блокам

2. Субоптимальность

• Не гарантирует минимизацию tₐ, только предоставляет субоптимальное решение
• Количество блоков 12N² находится на одном уровне, но может быть не минимальным

3. Ошибка Троттера

• Из-за того, что эффективные гамильтонианы не коммутируют, существует ошибка Trotterization
• Требуется разделение эволюции на nₜ шагов для контроля ошибки, увеличивая сложность цепи

4. Проблемы физической реализации

• Требуются произвольные однокубитные вентили, что предъявляет высокие требования к точности оборудования
• 12N² блоков для крупных систем остаётся значительной глубиной цепи

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

1. Универсальный гамильтониан источника: поиск эффективных методов использования произвольных двухчастичных гамильтонианов источника, избегая O(N⁴) затрат
2. Анализ ошибки Троттера: систематическое исследование связи между ошибкой и количеством шагов nₜ, оптимизация компромисса между точностью и эффективностью
3. Оптимизация количества блоков: исследование возможности существования методов построения с меньшим количеством блоков
4. Экспериментальная проверка: проверка данного протокола на платформах сверхпроводящих цепей, ионных ловушек и других
5. Оптимизация для конкретных задач: разработка специализированных методов оптимизации для конкретных классов гамильтонианов в химии и материаловедении
6. Расширение на отказоустойчивые вычисления: расширение данного метода на рамки отказоустойчивых квантовых вычислений

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

Достоинства

1. Инновационность метода ⭐⭐⭐⭐⭐

• Теоретический прорыв: преобразование NP-трудной задачи в полиномиально разрешимую, новаторский подход
• Математическая элегантность: хитроумное построение с использованием собственного разложения и ортогональных возмущений
• Стратегия разделяй и властвуй: независимая обработка каждого собственного вектора, снижение сложности

2. Достаточность экспериментов ⭐⭐⭐⭐

• Крупномасштабная проверка: тестирование N=1–50, 10⁴ экземпляров для каждого масштаба
• Статистический анализ: предоставлены среднее значение, максимум/минимум, сравнение с теоретической верхней границей
• Ограничение: отсутствие сравнения с конкретными физическими гамильтонианами (например, молекулярными гамильтонианами)

3. Убедительность результатов ⭐⭐⭐⭐⭐

• Доказательство масштабируемости: сложность O(N³) подтверждена численно
• Демонстрация практичности: находка tₐ~O(1) имеет большое значение для практических приложений
• Теоретические гарантии: предоставлены строгие математические доказательства (приложение B)

4. Ясность изложения ⭐⭐⭐⭐⭐

• Чёткая структура: логическая последовательность от определения задачи → метод → эксперименты → обсуждение
• Стандартизированная нотация: последовательное использование математических символов
• Вспомогательные графики: рисунок 1 (схема цепи) и рисунок 2 (результаты) эффективно передают информацию

Недостатки

1. Ограничения метода

• Ограничение гамильтониана источника: ограничение типом ZZ ограничивает область применения
• Отсутствие количественной оценки ошибки Троттера: недостаёт количественного соотношения между ошибкой и nₜ
• Неоптимальное количество блоков: 12N² может быть дополнительно сжато

2. Недостатки экспериментальной установки

• Синтетические данные: использование только случайных матриц B без тестирования реальных физических задач
• Отсутствие сравнительных экспериментов: отсутствие прямого сравнения tₐ с методами из литературы15-17 на одинаковых задачах
• Отсутствие анализа ошибок: отсутствие исследования влияния ошибки Троттера на верность моделирования

3. Недостаточный анализ

• Теоретическое объяснение tₐ~O(1): отсутствует строгое доказательство того, почему tₐ не растёт при условии нормализации
• Нижняя граница оптимального количества блоков: отсутствует обсуждение существования нижней границы ω(N²)
• Влияние шумов оборудования: отсутствует анализ влияния шумов реального квантового оборудования на данный протокол

Влияние

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

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

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

• Высокая: применимо к системам с 50+ кубитами
• Средняя: требует высокоточных произвольных SQG, предъявляет высокие требования к оборудованию
• Требует проверки: производительность на реальных физических задачах требует экспериментального подтверждения

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

• Полнота алгоритма: предоставлены полные математические формулы и шаги построения
• Потенциал открытого кода: метод основан на стандартной линейной алгебре, легко реализуется
• Явные параметры: все гиперпараметры (например, способ нормализации) имеют пояснения

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

Наиболее подходящие сценарии

1. Крупномасштабное квантовое моделирование: N>20, где вычислительная стоимость традиционных методов чрезмерна
2. Платформы с гамильтонианом источника типа ZZ: сверхпроводящие кубиты, ионные ловушки с естественным взаимодействием
3. Ограниченные ресурсы предварительной обработки: невозможность проведения крупномасштабной численной оптимизации

Менее подходящие сценарии

1. Малые системы: N<10, где численная оптимизация может быть более эффективной
2. Гамильтониан источника не типа ZZ: требуется вложенное применение, снижение эффективности
3. Требования экстремальной точности: ошибка Троттера может требовать большого количества nₜ шагов

Потенциальные области применения

• Квантовая химия: моделирование молекулярных гамильтонианов
• Физика конденсированного состояния: моделирование спиновых систем, модели Хаббарда
• Квантовое машинное обучение: проектирование цепей для вариационных квантовых алгоритмов
• Квантовая оптимизация: эффективная реализация алгоритмов типа QAOA

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

1. Lloyd (1996): Universal quantum simulators - фундаментальная основа теории квантового моделирования
2. Barenco et al. (1995): Elementary gates for quantum computation - теория универсальных наборов вентилей
3. Parra-Rodriguez et al. (2020): Digital-analog quantum computation - введение парадигмы DAQC
4. Dodd et al. (2002): Universal quantum computation using any entangling Hamiltonian - доказательство универсальности DAQC
5. Garcia-de-Andoin et al. (2024): Digital-analog quantum computation with arbitrary two-body Hamiltonians - предыдущие работы
6. Baßler et al. (2023, 2024): Time-optimal multi-qubit gates - связанные методы оптимизации

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