Классическое моделирование схем с реалистичными состояниями Готтесмана-Китаева-Преськилла нечётной размерности

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

• ID статьи: 2412.13136
• Название: Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states
• Авторы: Cameron Calcluth, Oliver Hahn, Juani Bermejo-Vega, Alessandro Ferraro, Giulia Ferrini
• Категория: quant-ph (квантовая физика)
• Дата публикации: 17 декабря 2024 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2412.13136

Аннотация

В данной работе предложен алгоритм классического моделирования для симуляции квантовых схем нечётной размерности на основе кодирования Готтесмана-Китаева-Преськилла (ГКП). Алгоритм специально оптимизирован для наиболее сложного, но практически релевантного сценария — кодовых слов с высокой степенью сжатия, но конечным сжатием. Алгоритм использует функцию Вигнера Зака-Гросса, которая обеспечивает положительное представление кодированных стабилизаторных состояний с бесконечным сжатием. Время выполнения алгоритма пропорционально отрицательности функции Вигнера, что позволяет эффективно моделировать определённые крупномасштабные схемы. Для стабилизаторных ГКП-состояний со сжатием 12 дБ алгоритм может моделировать схемы с числом мод до 1000, требуя количество выборок, менее чем в два раза превышающее количество, необходимое для одной входной моды, что резко контрастирует с существующими симуляторами.

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

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

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

Анализ значимости

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

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

1. Метод разложения по базису Фока: Поскольку кодовые слова ГКП включают бесконечно много энергетических уровней, прямое моделирование становится невозможным даже для нескольких ГКП-кубитов
2. Методы квазивероятностного распределения: Функция Вигнера кодированных ГКП-состояний имеет большую отрицательность, что приводит к экспоненциальному росту времени выполнения
3. Методы разложения гауссовых состояний: Неприменимы к практически релевантным высокосжатым ГКП-состояниям из-за их значительной негауссовости

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

1. Предложена многомодовая функция Вигнера Зака-Гросса (ЗГВ): Обобщена одномодовая функция ЗГВ на n-модовый случай и доказано, что она удовлетворяет модифицированным аксиомам Стратоновича-Вейля
2. Установлена связь между функцией ЗГВ и функцией Вигнера Гросса: Доказано соответствие между функцией ЗГВ ГКП-состояний непрерывной переменной и функцией Вигнера Гросса кодированных логических состояний
3. Разработан эффективный алгоритм моделирования:
   • Предоставлен алгоритм слабого моделирования для идеальных стабилизаторных ГКП-состояний
   • Предоставлен алгоритм вероятностной оценки для ГКП-состояний с конечным сжатием
4. Достигнуто значительное улучшение производительности: Для стабилизаторных ГКП-состояний со сжатием 12 дБ моделирование схемы с 1000 модами требует менее чем в два раза больше выборок, чем для одной моды

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

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

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

• Вход: Произвольные ГКП-кодированные кубиты с параметром антисжатия Δ
• Эволюция: Унитарные операции, описываемые целочисленными симплектическими матрицами, и произвольные смещения (включая все кодированные операции Клиффорда)
• Измерение: Модульные измерения в базисе положения $M_Z(s) = \frac{1}{d\ell}\sum_n e^{-i\ell s \cdot n}T_{(0,n)}$

Определение функции ЗГВ

Определение 1: Многомодовая функция ЗГВ нечётной размерности определяется как: $W_{\hat{\rho}}(\eta) = \text{Tr}(\hat{\rho}\hat{A}_\eta)$

где оператор фазовой точки имеет вид: $\hat{A}_\eta = \frac{1}{(2\pi)^n}\sum_{a \in \mathbb{Z}^{2n}} e^{i\ell[a,\eta]+i\pi a_X^T a_Z}\hat{T}_a$

где $\eta \in [0,d\ell)^{2n}$, $d$ — положительное нечётное число, $\hat{T}_a = e^{i\pi a_X^T a_Z/d}\hat{T}_{a_X}\hat{T}_{a_Z}$.

Ключевые теоретические результаты

Лемма 1: Связь между функцией ЗГВ и функцией Вигнера Гросса $W_{\hat{\rho}}(\eta) = \bar{W}_{\bar{\rho}(u)}(t)$ где $\eta = \ell(u+t)$, $u = \frac{1}{\ell}\eta \bmod 1$, $t = \frac{1}{\ell}\eta - u$.

Теорема 1: Преобразование функции ЗГВ при симплектической эволюции Для гауссовой унитарной эволюции, описываемой целочисленной симплектической матрицей $S$, функция ЗГВ эволюционирующего состояния имеет вид: $W_{\hat{U}_S\hat{\rho}\hat{U}_S^\dagger}(\eta) = W_{\hat{\rho}}(S\eta - t)$

Теорема 2: Вычисление вероятностей измерения Распределение вероятностей измерения логического оператора $\hat{Z}_L^{\otimes n}$ имеет вид: $\text{Tr}(\hat{\rho}\hat{M}_Z(s)) = \int d\eta_X W_{\hat{\rho}}\left(\begin{pmatrix}\eta_X\\s\end{pmatrix}\right)$

Алгоритм моделирования

Слабое моделирование идеальных стабилизаторных ГКП-состояний

1. Выборка вектора $\eta$ из эффективного распределения вероятностей $W_{\hat{\rho}_0}(\eta)$
2. Преобразование под действием симплектической матрицы $S$ и вектора смещения $t$: $S\eta - t$
3. Применение линейного оператора смещения: $S\eta - t + c$
4. Результат измерения определяется второй половиной вектора

Вероятностная оценка для ГКП-состояний с конечным сжатием

Использование метода выборки с квазивероятностным распределением: $\text{Pr}(\eta) = \frac{1}{M_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}}|W_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}(\eta)|$

Требуемое количество выборок: $N = \frac{2}{\epsilon^2}M_{\hat{\rho}_0}^2\log(2/\delta)$

где $M_{\hat{\rho}}$ — отрицательность функции ЗГВ (отрицательный объём).

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

Теоретическая аналитическая база

Работа сосредоточена на теоретическом анализе и численных расчётах, уделяя особое внимание:

1. Вычисление отрицательности функции ЗГВ: Для ГКП-состояний с различными уровнями сжатия
2. Анализ сложности выборки: На основе поведения отрицательности при масштабировании
3. Сравнение с существующими методами: Теоретическое сравнение сложности

Оценочные показатели

• Логарифмическая отрицательность ЗГВ: $\log M_{\hat{\rho}}$, характеризующая сложность моделирования
• Масштабирование количества выборок: Рост количества выборок относительно одномодового случая
• Ёмкость по числу мод: Максимальное число мод, которое можно эффективно моделировать

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

Функция ЗГВ для ГКП-состояний с конечным сжатием

Посредством явного вычисления получена функция ЗГВ реалистичного логического 0-состояния ГКП: $W_{\text{GKP}}\left(\begin{pmatrix}u\\v\end{pmatrix}\right) \propto \vartheta(\Gamma; z)$

где $z = (v/(d\ell), -u/(d\ell), 0, 0)^T$, $\Gamma$ — комплексная матрица, определяемая параметром сжатия $\Delta$.

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

Для $\Delta = 0.25$ (соответствует сжатию 12 дБ):

• Одномодовая отрицательность: $M_{\hat{\rho}_0^\Delta} \approx e^{3 \times 10^{-4}}$
• Моделирование тысячи мод: Благодаря мультипликативному характеру отрицательности моделирование 1000 входных мод требует менее чем в два раза больше выборок, чем для одной моды
• Улучшение производительности: На несколько порядков превосходит существующие методы

Анализ магических состояний

Для π-состояния $|\psi_\pi\rangle = \frac{1}{\sqrt{3}}(|0_L\rangle + |1_L\rangle - |2_L\rangle)$:

• Демонстрирует значительную отрицательность при всех значениях $\Delta$
• Значительные вычислительные затраты при моделировании крупномасштабных схем, что соответствует ожиданиям

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

Развитие алгоритмов классического моделирования

1. Теорема Готтесмана-Нилла: Полиномиальное моделирование схем стабилизаторов
2. Методы квазивероятности: Алгоритмы выборки на основе функции Вигнера Гросса
3. Моделирование бозонных систем: Использование положительности квазивероятностных распределений, тензорных сетей, звёздчатого представления и других методов

Современное состояние моделирования ГКП-кодирования

Существующие методы в основном ограничены:

• Идеальными (бесконечно сжатыми) ГКП-состояниями
• Гауссовыми операциями и гомодинным измерением
• Неспособностью обрабатывать практически релевантные случаи с конечным сжатием

Инновационные аспекты данной работы

По сравнению с существующими работами, данная статья впервые реализует эффективное моделирование ГКП-состояний с конечным сжатием, особенно достигая значительного улучшения производительности в режиме высокого сжатия.

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

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

1. Теоретический вклад: Установлена полная теоретическая база функции ЗГВ, включая многомодовое обобщение и законы эволюции
2. Алгоритмическая инновация: Предложен эффективный алгоритм моделирования, применимый к реалистичным ГКП-состояниям
3. Практическая ценность: Предоставляет мощный инструмент для эталонного тестирования бозонных квантовых процессоров

Ограничения

1. Ограничение нечётной размерностью: Текущий метод применим только к ГКП-кубитам нечётной размерности; расширение на чётные размерности (например, кубиты) остаётся открытой проблемой
2. Ограничение операций: Ограничено операциями, описываемыми целочисленными симплектическими матрицами; не может обрабатывать произвольное одномодовое сжатие
3. Модель шума: Не учитывает влияние реальных шумов, таких как потери, на отрицательность функции ЗГВ

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

1. Расширение на ГКП-кубиты: Разработка методов для двумерного случая
2. Более широкие бозонные коды: Обобщение框架 на другие бозонные коды помимо ГКП
3. Квазивероятностная база: Развитие комплексной квазивероятностной теории, специально разработанной для бозонных кодов

Углубленная оценка

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

1. Теоретическая строгость: Полные математические выводы, детальные доказательства, прочная теоретическая база
2. Высокая практичность: Решает проблему моделирования практически релевантных ГКП-состояний с конечным сжатием
3. Выдающаяся производительность: Достигает улучшения на несколько порядков в режиме высокого сжатия
4. Методологическая инновация: Многомодовое обобщение и применение функции ЗГВ представляет важный теоретический вклад

Недостатки

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

Влияние

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

Сценарии применения

• Эталонное тестирование и верификация бозонных квантовых процессоров
• Теоретический анализ схем отказоустойчивых квантовых вычислений с ГКП-кодированием
• Исследование классического моделирования квантовых вычислений с непрерывными переменными

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

Статья ссылается на 68 соответствующих работ, включая:

• Оригинальные работы по ГКП-кодированию
• Классические работы по алгоритмам квантового моделирования
• Последние достижения в области бозонных квантовых вычислений
• Теоретические основы функций Вигнера и квазивероятностных распределений

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