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

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

• ID статьи: 2511.00842
• Название: Kostant relation in filtered randomized benchmarking for passive bosonic devices
• Автор: David Amaro-Alcalá (Institute of Physics, Slovak Academy of Sciences)
• Классификация: quant-ph (квантовая физика)
• Дата публикации: 4 ноября 2025
• Ссылка на статью: https://arxiv.org/abs/2511.00842

Аннотация

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

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

Проблема, которую необходимо решить

Характеризация пассивных бозонных устройств (passive bosonic devices) является важным этапом развития квантовых вычислений с непрерывными переменными. Существующие схемы рандомизированного бенчмаркинга бозонных систем, хотя и наследуют преимущества стандартного рандомизированного бенчмаркинга (такие как устойчивость к ошибкам подготовки и измерения), имеют два критических недостатка:

1. Вычислительная сложность: требуется вычисление перманентов матриц (matrix permanents), что является вычислительно сложной задачей (проблема #P-полноты). Более того, необходимые перманенты должны определяться индивидуально посредством сложного разложения через коэффициенты Клебша-Гордана.
2. Экспериментальная сложность: требуется подготовка состояний Фока и использование детекторов с разрешением по числу фотонов (photon-number-resolving detectors), что представляет значительный вызов для большинства лабораторий.

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

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

• контроля качества квантовых компьютеров;
• понимания и смягчения влияния шума;
• проверки квантового преимущества.

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

Основные ограничения исходной схемы 11 включают:

• необходимость вычисления нескольких перманентов матриц (количество не менее ♯λ - 1 + dλ);
• необходимость вычисления коэффициентов Клебша-Гордана для SU(m);
• требование сложного экспериментального оборудования (подготовка состояний Фока и разрешение по числу фотонов);
• сложный процесс анализа данных, зависящий от конкретных начальных состояний и измерений.

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

Данная работа направлена на сохранение преимуществ исходной схемы при одновременном фундаментальном упрощении вычислительных и экспериментальных требований путем использования соотношения Костанта (математического соотношения, связывающего имманты с состояниями нулевого веса).

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

1. Предложение новых фильтрующих функций на основе имманантов: посредством соотношения Костанта процесс фильтрации выражается как вычисление имманантов, полностью избегая вычисления коэффициентов Клебша-Гордана.
2. Снижение вычислительной сложности:
   • сокращение количества необходимых перманентов с ♯λ - 1 + dλ до всего 1;
   • количество необходимых имманантов равно ♯λ - 1 (обычно значительно меньше, чем в исходной схеме);
   • исключение необходимости вычисления коэффициентов Клебша-Гордана.
3. Упрощение экспериментальных требований: доказано, что слабые когерентные состояния (weak coherent states) и измерения интенсивности достаточны для полной характеризации без необходимости в детекторах с разрешением по числу фотонов или сложной подготовке состояний Фока.
4. Теоретические гарантии: доказано, что новые фильтрующие функции по-прежнему обеспечивают однопоказательное затухание и позволяют точно оценить показатель верности F(E).
5. Повышенная устойчивость: метод обладает естественной устойчивостью к потерям фотонов и ошибкам усиления, фильтрующее выражение может быть определено заранее и не зависит от выбора начального состояния или измерения.

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

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

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

• n фотонов, входящих в m-портовый интерферометр;
• K случайно выбранных последовательностей унитарных вентилей Ug_s глубины g;
• начальное состояние ϱ и измерение E (могут быть слабые когерентные состояния и измерения интенсивности).

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

• показатель верности шумового канала E: F(E) = d_λ^(-2) Σ_μ d_μ p_μ(E)

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

• предполагается марковская модель шума, независимая по вентилям и времени: Ũ = E ∘ U

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

1. Математическая основа

Гильбертово пространство системы H_m^n представляет пространство состояний n неразличимых фотонов в m модах. Унитарная операция U ∈ SU(m) действует следующим образом:

U(|n⟩) = (∏_k 1/√(n_k!) (U(a_k†))^(n_k)) |0⟩

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

Γ: U ↦ U ⊗ Ū

2. Разложение по неприводимым представлениям

Тензорное произведение симметричного неприводимого представления λ = (n, 0, ..., 0) и его дуального λ* разлагается как:

Γ := λ ⊗ λ* ≅ ⊕_μ μ

Представления μ вычисляются через формулу Пьери: добавление n клеток в различные столбцы диаграммы Юнга для λ*.

3. Ключевое нововведение: фильтрующие функции на основе имманантов

Теорема 2 (фильтрующие функции на основе имманантов): Определим фильтрующую функцию как:

f_(g,s)_Imm,μ(U^(g,s)) := Imm_μ(U^(g,s))

Тогда:

Φ_g^(f) := E_s [f_μ^(g,s)(U^(g,s)) d_ϱ,E^(g,s)(Ũ^(g,s))] = κ p_μ^(g-1)

где κ — константа, независимая от характеризации.

4. Применение соотношения Костанта

Теорема 1 (соотношение Костанта):

Σ_(|ζ_κ⟩∈Z_κ) ⟨ζ_κ|Γ(U)|ζ_κ⟩ = Imm_κ(U)

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

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

1. Исключение проекционных операторов

Фильтрующая функция исходной схемы имеет вид:

f_λ,orig := ⟨ϱ|P_λ S^+ Γ(U^g_s)† |E⟩

требует вычисления проекционного оператора P_μ на подпространство неприводимого представления μ, что требует коэффициентов Клебша-Гордана.

Предложенная схема использует имманты напрямую:

f_Imm,μ := Imm_μ(U^(g,s))

полностью избегая проекционных операторов и коэффициентов Клебша-Гордана.

2. Использование структуры состояний нулевого веса

Посредством соотношения Костанта имманты могут быть выражены как сумма по состояниям нулевого веса:

Imm_μ(U) = Σ_(|ζ_μ⟩∈Z_μ) ⟨ζ_μ|Γ(U)|ζ_μ⟩

Количество состояний нулевого веса значительно меньше полного базиса, и их структура проста, что облегчает вычисления.

3. Гарантия однопоказательного затухания

Ключевые этапы доказательства:

• усреднение по состояниям нулевого веса;
• использование свойств операции twirling TE;
• в силу ортогональности различных неприводимых представлений каждый p_μ затухает независимо.

Математически:

⟨ζ_μ^(i)|S_Immμ^(i) T[E]^g |ϱ̃⟩ = (Σ_j s_(i,j)⟨ζ_μ^(j)|ϱ̃⟩) p_μ^g

Процесс обработки данных

1. Сбор данных: построение матрицы данных D, где D_(g,s) = ⟨Ẽ|Γ(Ũ^(g,s))|ϱ̃⟩
2. Вычисление имманантов: построение матрицы фильтрации F_μ, где (F_μ)_(g,s) = Imm_μ(U^(g,s))
3. Произведение Адамара: вычисление Φ_g = Σ_s (F_μ ⊙ D)_(g,s)
4. Экспоненциальная аппроксимация: аппроксимация {g, Φ_g} экспоненциальной функцией для извлечения p_μ
5. Вычисление верности: использование F(E) = d_λ^(-2) Σ_μ d_μ p_μ

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

Теоретическая проверка

Данная работа является преимущественно теоретической и методологической, с проверкой корректности метода посредством математических доказательств и конкретных примеров.

Конкретный пример: проверка соотношения Костанта для SU(3)

Параметры системы:

• унитарная группа: SU(3);
• неприводимое представление: μ = (2,1,0);
• состояния нулевого веса:

|ζ_(2,1,0)^(0)⟩ = |2 1 0; 1 1; 1⟩
|ζ_(2,1,0)^(1)⟩ = |2 1 0; 2 0; 1⟩

Процесс проверки:

1. Вычисление D-функций с использованием схем Гельфанда-Цетлина
2. Вычисление диагональных матричных элементов состояний нулевого веса (уравнения A2a и A3)
3. Вычисление имманта с использованием таблицы характеров (уравнение A4)
4. Проверка: Σ_(i=0,1) ⟨ζ^(i)|Γ(U)|ζ^(i)⟩ = Imm_(2,1,0)(U)

Результаты показывают полное совпадение двух методов вычисления, что подтверждает соотношение Костанта.

Осуществимость слабых когерентных состояний и измерений интенсивности

Расширенное гильбертово пространство: Рассмотрим (H_e)m^n := ⊕(n≥0) H_m^n, допускающее изменение числа фотонов.

Ключевое наблюдение: Для входа в виде слабого когерентного состояния:

ϱ = |0,0⟩ + α|1,0⟩

Разложение по неприводимым представлениям в расширенном пространстве:

(0 ⊕ λ) ⊗ (0 ⊕ λ*) = 0 ⊕ λ ⊕ λ ⊕ λ* ⊕ Γ

Целевое неприводимое представление по-прежнему появляется только один раз, поэтому процесс фильтрации остается эффективным.

Экспериментальные преимущества:

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

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

Сравнение вычислительной сложности

Таблица II - итоговое сравнение:

Подробный анализ:

1. Вычисление перманентов: сокращение с нескольких (не менее ♯λ - 1 + d_λ) до всего 1
2. Вычисление имманантов: хотя требуется ♯λ - 1 имманантов:
   • вычислительная сложность имманантов ниже, чем перманентов;
   • могут быть вычислены непосредственно с использованием существующих пакетов программного обеспечения (например, Wolfram);
   • структура проще (обычно одностолбцовые или близкие к одностолбцовым диаграммы Юнга).
3. Коэффициенты Клебша-Гордана: полностью исключены, что является наиболее значительным упрощением.

Численные примеры

Случай SU(3):

• для n=2 фотонов, m=3 модов;
• количество неприводимых представлений ♯λ = 3;
• размерность d_λ = 3.

Требования исходной схемы:

• перманенты: не менее 3 - 1 + 3 = 5;
• коэффициенты Клебша-Гордана: требуется полное вычисление.

Требования предложенного метода:

• перманенты: 1;
• имманты: 3 - 1 = 2;
• коэффициенты Клебша-Гордана: 0.

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

Посредством примера SU(3) в Приложении A конкретно проверено:

1. Правильность идентификации состояний нулевого веса
2. Точность вычисления D-функций
3. Эквивалентность имманта и следа состояний нулевого веса
4. Теоретическая гарантия однопоказательного затухания

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

Рандомизированный бенчмаркинг

1. Стандартный RB: Emerson и др. 5, Magesan и др. 6, Knill и др. 7 установили основу рандомизированного бенчмаркинга для квантовых вентилей конечной размерности
2. Расширение на кудиты: Amaro-Alcalá и др. 8, Jafarzadeh и др. 9 расширили RB на системы высокой размерности
3. Общая основа: Helsen и др. 10 предоставили унифицированную теоретическую основу

Бенчмаркинг бозонных систем

1. Исходная схема: Arienzo и др. 11 впервые предложили схему бозонного RB
2. Параллельные работы: Wilkens и др. 12 исследовали бенчмаркинг бозонной и фермионной динамики
3. Вклад данной работы: значительное упрощение вычислительных и экспериментальных требований 11

Инструменты теории представлений

1. Теория имманантов: классическая работа Littlewood 23
2. Соотношение Костанта: Kostant 15 установил связь между имманантами и состояниями нулевого веса
3. Приложения: de Guise и др. 25,26 применили это в квантовой оптике

Обработка шума

1. Проблема калибровочной свободы: Proctor и др. 19 указали на проблему калибровочной свободы в RB
2. Случайная компиляция: Wallman и др. 20,21 доказали, что случайная компиляция может преобразовать когерентный шум в некогерентный
3. Позиция данной работы: комбинирование RC и RB может эффективно решить проблему калибровки

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

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

1. Теоретический вклад: посредством соотношения Костанта установлены новые фильтрующие функции на основе имманантов с теоретической гарантией однопоказательного затухания и точной оценки верности.
2. Вычислительное упрощение:
   • исключение вычисления коэффициентов Клебша-Гордана;
   • сокращение количества перманентов с O(♯λ + d_λ) до 1;
   • количество имманантов O(♯λ) с простой структурой.
3. Экспериментальное упрощение:
   • замена состояний Фока слабыми когерентными состояниями;
   • замена разрешения по числу фотонов измерениями интенсивности;
   • устойчивость к потерям и ошибкам усиления.
4. Повышенная практичность:
   • фильтрующие функции могут быть определены заранее;
   • упрощение процесса анализа данных;
   • применимость к большему числу экспериментальных платформ.

Ограничения

1. Ограничение пассивными преобразованиями: текущая схема применима только к пассивным бозонным преобразованиям (линейная оптика), не включая активные преобразования (такие как операции сжатия).
2. Предположения о модели шума: предполагается марковский шум, независимый по вентилям и времени; реальные системы могут быть более сложными.
3. Вызовы расширения: расширение на активные бозонные преобразования сталкивается с теоретическими трудностями, поскольку соответствующая группа преобразований некомпактна.
4. Отсутствие экспериментальной проверки: работа является преимущественно теоретической и лишена экспериментальных данных.
5. Приближение слабых когерентных состояний: при использовании слабых когерентных состояний требуется, чтобы α было достаточно малым для обеспечения действительности приближения.

Будущие направления

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

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

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

1. Теоретическая строгость:
   • основана на прочной математической основе (соотношение Костанта);
   • полные доказательства (теорема 2);
   • конкретные примеры для проверки (случай SU(3)).
2. Значительная инновационность:
   • первое применение соотношения Костанта к квантовому бенчмаркингу;
   • умелое избежание коэффициентов Клебша-Гордана как главного узкого места;
   • фундаментальное упрощение структуры проблемы.
3. Высокая практическая ценность:
   • значительное снижение вычислительных затрат (исключение компонентов экспоненциальной сложности);
   • значительное упрощение экспериментальных требований (слабые когерентные состояния + измерения интенсивности);
   • прямое применение к существующим данным.
4. Универсальность метода:
   • применимость к системам с произвольным числом фотонов n и модов m;
   • устойчивость к потерям и ошибкам усиления;
   • фильтрующие функции определяются заранее и не зависят от конкретной реализации.
5. Ясность изложения:
   • логичная структура с четкой логикой;
   • подробные математические выводы;
   • конкретные примеры и приложения.

Недостатки

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

Влияние

1. Вклад в область:
   • предоставление практического инструмента характеризации для квантовых вычислений с непрерывными переменными;
   • демонстрация новых приложений классической математики (теория представлений) в квантовой информации;
   • возможное вдохновение методов бенчмаркинга для других квантовых систем.
2. Практическая ценность:
   • значительное снижение барьера к внедрению;
   • возможность для большего числа лабораторий проводить характеризацию квантовых устройств;
   • возможное ускорение развития квантовых вычислений с непрерывными переменными.
3. Воспроизводимость:
   • полные теоретические выводы;
   • четкое описание метода;
   • возможность реализации с использованием существующих программных инструментов;
   • требуется экспериментальная проверка для подтверждения практичности.
4. Потенциальное влияние:
   • возможность стать стандартным методом бенчмаркинга для бозонных систем;
   • предоставление методологических уроков для других задач с непрерывными переменными;
   • продвижение применения инструментов теории представлений в квантовой информации.

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

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

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

11 M. Arienzo et al., "Bosonic randomized benchmarking with passive transformations," PRX Quantum 6, 020305 (2025) - исходная схема, улучшенная в данной работе

15 B. Kostant, "Immanant inequalities and 0-weight spaces," J. Am. Math. Soc. 8, 181 (1995) - исходная статья о соотношении Костанта

25 H. de Guise et al., "D-functions and immanants of unitary matrices and submatrices," J. Phys. A 49, 09LT01 (2016) - применение соотношения Костанта в квантовой оптике

20 J. J. Wallman and J. Emerson, "Noise tailoring for scalable quantum computation via randomized compiling," Phys. Rev. A 94, 052325 (2016) - метод случайной компиляции

Общая оценка: Это высококачественная теоретическая методологическая работа, которая посредством умелого применения соотношения Костанта значительно упрощает рандомизированный бенчмаркинг бозонных систем. Теория строга, инновационность значительна, практическая ценность высока. Основной недостаток — отсутствие экспериментальной проверки и ограничение области применения. Если последующие эксперименты подтвердят его эффективность, это может стать важной работой в данной области. Рекомендуется читателям обратить внимание на последующие экспериментальные работы по проверке и возможные исследования по расширению метода.