Оценка математических ожиданий с параметризованными квантовыми схемами

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

• ID статьи: 2407.19499
• Название: Expectation value estimation with parametrized quantum circuits
• Авторы: Bujiao Wu, Lingyu Kong, Yuxuan Yan, Fuchuan Wei, Zhenhuan Liu
• Классификация: quant-ph (квантовая физика)
• Дата публикации: июль 2024 г. (препринт arXiv, версия 2 обновлена 16 октября 2025 г.)
• Ссылка на статью: https://arxiv.org/abs/2407.19499

Аннотация

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

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

Определение проблемы

Оценка линейных свойств квантового состояния Tr(ρH) является центральной задачей квантовой информатики, где ρ — квантовое состояние, а H — наблюдаемая. Такие задачи широко встречаются в:

1. Квантовой химии: расчет энергии основного состояния молекулы
2. Физике многих тел: измерение корреляционных функций
3. Квантовой информации: оценка верности состояния

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

1. Протокол классической тени (Classical Shadow):
   • Локальный протокол CS имеет сложность выборки O(4^k) для k-локальных наблюдаемых
   • Глобальный протокол CS, хотя и достигает сложности O(1), требует схемы логарифмической глубины
   • Оба являются "независимыми от измерения", не используя априорную информацию о целевой наблюдаемой
2. Метод разложения Паули:
   • Ограничен реализацией на схемах Клиффорда
   • Разложение ограничено наблюдаемыми Паули
   • Может требовать глубоких схем или высокой сложности выборки

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

Существующие методы сталкиваются со следующими проблемами на современных квантовых устройствах:

• Ограничение глубины схемы
• Ограничения связности
• Влияние шума
• Недостаточное использование информации о наблюдаемой

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

1. Унифицированная структура: предложена общая структура для оценки линейных свойств с использованием параметризованных квантовых схем, объединяющая существующие протоколы разложения Паули
2. Два алгоритма разложения:
   • Разложение жадной проекции (GPD): применимо к общим гамильтонианам
   • Разложение тензорной сети (TND): применимо к гамильтонианам с компактным представлением тензорной сети
3. Теоретические нижние границы: выведены фундаментальные нижние границы сложности выборки, необходимой для оценки целевой наблюдаемой с использованием заданной неглубокой квантовой схемы
4. Численная верификация: алгоритмические преимущества проверены на разреженных/плотных гамильтонианах и оценке внутреннего произведения определителей Слейтера

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

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

Дано:

• Неизвестное квантовое состояние ρ
• Целевая наблюдаемая H
• Параметризованная квантовая схема глубины L: U_L(θ)
• Требования к точности ε и вероятности успеха 1-δ

Цель: оценить Tr(ρH) с минимизацией сложности выборки

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

Структура разделена на классический и квантовый этапы:

Классический этап: разложение целевой наблюдаемой на $H \approx \sum_{k=1}^K U_L(\theta^{(k)})^\dagger \Lambda_k U_L(\theta^{(k)})$ где Λ_k — вещественные диагональные матрицы

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

• Выборка члена k с вероятностью p_k ∝ ||Λ_k||_2
• Выполнение U_L(θ^{(k)}) и измерение в вычислительном базисе
• Применение метода медианы средних для получения финальной оценки

Алгоритм разложения жадной проекции (GPD)

Основная идея: итеративный поиск оптимального приближения члена U_L(θ)†ΛU_L(θ)

Процедура алгоритма:

1. Инициализация H^{(0)} = H, k = 0
2. Пока ||H^{(k)}||_2 ≥ ε:
   • Решение оптимизационной задачи: θ^{(k)} = argmin_θ ||U_L(θ)H^{(k)}U_L†(θ) - diagU_L(θ)H^{(k)}U_L†(θ)||_F
   • Установка Λ_k = diagU_L(θ^{(k)})H^{(k)}U_L†(θ^{(k)})
   • Обновление H^{(k+1)} = H^{(k)} - U_L†(θ^{(k)})Λ_k U_L(θ^{(k)})
   • k = k + 1

Анализ сложности: время классической обработки составляет O(poly(n)·2^{ωn}), где ω ≈ 2.37 — показатель матричного умножения

Алгоритм разложения тензорной сети (TND)

Область применения: целевой гамильтониан имеет эффективное представление матричного произведения операторов (MPO)

Цель оптимизации: минимизация функции потерь $L = ||H - \sum_k U_L(\theta_k)^\dagger \Λ_k U_L(\theta_k)||_F^2$

Ключевые технологии:

• Представление U_L(θ_k) как унитарной тензорной сети глубины L
• Представление Λ_k в форме MPO
• Использование сжатия тензорной сети для вычисления функции потерь
• Оптимизация параметров {θ^{(k)}, Λ_k} методом градиентного спуска

Анализ сложности выборки

Верхняя граница: алгоритм 1 требует T = O(||Λ||_1^2 log(1/δ)/ε_2^2) выборок, где ||Λ||_1 — сумма всех ||Λ_k||_2

Нижняя граница: любая адаптивная стратегия с одной копией, использующая параметризованную схему U_L(θ), требует $T = Ω\left(\frac{\text{Tr}(H_0^2)^2}{\varepsilon^2 \delta(H_0) 4^n}\right)$ где H_0 — бесследовая часть H, δ(H_0) — квадрат максимального математического ожидания H_0 на множестве достижимых состояний

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

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

1. Оценка энергии основного состояния разреженного гамильтониана: система из 8 кубитов, 64 ненулевых элемента
2. Оценка математического ожидания плотного гамильтониана: 4 кубита, случайная эрмитова матрица
3. Оценка внутреннего произведения определителей Слейтера: система из 3 кубитов, внутреннее произведение τ-определителя Слейтера и чистого состояния

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

• Протокол классической тени: глобальная CS и локальная CS
• Методы разложения Паули: Derandomized, C-LBCS, SG, Adaptive, OGM и др.
• Специализированные методы: классическая тень фермионов (FCS)

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

• Параметризованные вентили: вентили iSWAP + тензорное произведение двух произвольных однокубитовых вентилей
• Алгоритм GPD: глубина L=4, K=20 или 80 членов разложения
• Алгоритм TND: глубина L=1, K=3 члена разложения

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

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

Разреженный гамильтониан (8 кубитов):

• При 25848 выборках ошибка GPD составляет 0.030, значительно лучше, чем лучший метод сравнения OGM с ошибкой 0.097
• С увеличением количества выборок GPD постоянно сохраняет наименьшую ошибку

Плотный гамильтониан (4 кубита):

• При 25848 выборках ошибка GPD составляет 0.046, лучше, чем лучший метод сравнения OGM с ошибкой 0.053
• Преимущество более выраженно при меньшем количестве выборок

Внутреннее произведение определителей Слейтера (3 кубита):

• GPD достигает наименьшей ошибки при всех количествах выборок
• При 25848 выборках ошибка 0.009, лучший метод сравнения — 0.012

Анализ сходимости

Численные результаты показывают:

1. При фиксированном количестве членов разложения K расстояние Фробениуса уменьшается с увеличением глубины схемы L
2. При фиксированной глубине схемы расстояние Фробениуса экспоненциально убывает с увеличением количества членов разложения K

Производительность метода тензорной сети

Для гамильтонианов с низкой размерностью связи:

• Метод TND использует только 3 члена разложения и глубину схемы 1
• При 18000 шагов ошибка 0.050, лучше, чем традиционные методы

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

Обучение квантовых состояний

• Квантовая томография: полная реконструкция квантового состояния, экспоненциальная сложность
• Теневая томография: классическое описание состояния, поддержка оценки множественных свойств

Протоколы случайного измерения

• Локальные измерения: однокубитовая группа Клиффорда, подходит для локальных наблюдаемых
• Глобальные измерения: глобальная группа Клиффорда, требует глубокой схемы
• Неглубокие схемы: компромиссный вариант, но все еще недостаточно использует информацию о наблюдаемой

Методы разложения Паули

• Основаны на разложении Паули наблюдаемой
• Реализованы через схемы Клиффорда и измерение в вычислительном базисе
• Предложенная структура объединяет эти методы

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

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

1. Предложенная структура успешно объединяет существующие протоколы измерения и расширяет их на общие параметризованные схемы
2. Алгоритмы GPD и TND значительно превосходят существующие методы в различных сценариях
3. Установленные теоретические нижние границы раскрывают фундаментальные ограничения неглубоких схем в задачах квантового обучения

Ограничения

1. Алгоритм GPD:
   • Сложность классической оптимизации остается высокой
   • Жадная стратегия не гарантирует глобальный оптимум
   • Теоретический анализ количества членов разложения K затруднен
2. Алгоритм TND:
   • Применим только к гамильтонианам с эффективным представлением MPO
   • Требует дополнительных методов оптимизации тензорной сети
3. Теоретические нижние границы:
   • Могут быть недостаточно точными для низкоранговых наблюдаемых (таких как верность)
   • Зависят от точной оценки параметра способности схемы δ(H_0)

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

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

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

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

1. Теоретические вклады:
   • Предоставляет унифицированную структуру, интегрирующую множество существующих методов
   • Устанавливает важные теоретические нижние границы, углубляя понимание способностей неглубоких схем
2. Методологические инновации:
   • Алгоритм GPD применим к общим гамильтонианам, высокая практичность
   • Алгоритм TND оптимизирован для специфических структур, высокая эффективность
   • Полное использование априорной информации о наблюдаемой
3. Достаточные эксперименты:
   • Охватывают множество сценариев применения (разреженные/плотные гамильтонианы, оценка внутреннего произведения)
   • Сравнение с несколькими основными методами, убедительные результаты
   • Предоставляют анализ сходимости и среднюю производительность

Недостатки

1. Проблемы масштабируемости:
   • Сложность классической оптимизации GPD экспоненциально растет с количеством кубитов
   • Практичность для крупномасштабных систем требует проверки
2. Ограничения экспериментов:
   • Масштаб численных экспериментов относительно небольшой (максимум 8 кубитов)
   • Отсутствие верификации на реальных квантовых устройствах
   • Не учитывается влияние шума на производительность алгоритма
3. Теоретические пробелы:
   • Существует значительный разрыв между верхней и нижней границами
   • Отсутствуют строгие теоретические гарантии скорости сходимости количества членов разложения K

Влияние

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

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

1. Квантовая химия: расчет энергии основного состояния молекулы и свойств
2. Квантовое моделирование: измерение корреляционных функций систем многих тел
3. Квантовое машинное обучение: отображение признаков и методы ядра
4. Квантовая оптимизация: оценка математического ожидания целевой функции
5. Квантовая коррекция ошибок: оценка верности кодовых слов и производительность исправления ошибок

Библиография

Статья цитирует 66 связанных работ, охватывающих важные исследования в основных областях обучения квантовых состояний, случайных измерений, классической тени и разложения Паули, обеспечивая прочную теоретическую основу для исследования.