Эффективная эмуляция квантового оборудования на нейтральных атомах

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

• ID статьи: 2510.09813
• Название: Efficient Emulation of Neutral Atom Quantum Hardware
• Авторы: Кемаль Бидживев, Стефано Грава, Пабло ле Энафф, Мауро Мендизабал, Эли Мерхей, Антон Кюэль (PASQAL)
• Категория: quant-ph (квантовая физика)
• Дата публикации: 14 октября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2510.09813

Аннотация

Моделирование динамики массивов нейтральных атомов представляет собой сложную вычислительную задачу. Для её решения авторы представляют два эмулятора — emu-sv и emu-mps, разработанные как вычислительные бэкенды пакета Pulser компании Pasqal. Emu-sv предназначен для высокоточного моделирования вектора состояния и может эмулировать системы с до 27 кубитами на GPU A100 40GB, что идеально подходит для случаев, требующих численно точных результатов. В отличие от этого, emu-mps использует представление матричного произведения состояний (MPS) и другие контролируемые приближения для эффективного моделирования больших массивов атомов при сохранении управляемой ошибки. Сравнительные тесты показывают, что оба эмулятора обеспечивают значительное ускорение по сравнению с универсальными решателями, такими как QuTiP. Кроме того, авторы предоставляют практическое руководство по выбору между двумя эмуляторами.

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

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

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

Значимость исследования

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

Технический фундамент

Исследование основано на зависящем от времени уравнении Шрёдингера: $i\frac{d|\psi\rangle}{dt} = \hat{H}(t)|\psi\rangle$

где гамильтониан имеет вид: $\hat{H}(t) = \sum_{i=1}^N \frac{\Omega_i(t)}{2}\hat{\sigma}^x_i - \sum_{i=1}^N \delta_i(t)\hat{n}_i + \sum_{i>j} \frac{C}{|r_{ij}|^6}\hat{n}_i\hat{n}_j$

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

1. Разработка двух высокопроизводительных квантовых эмуляторов:
   • emu-sv: точный эмулятор на основе вектора состояния, поддерживающий до 27 кубитов
   • emu-mps: приближённый эмулятор на основе матричного произведения состояний для систем большего масштаба
2. Достижение значительного повышения производительности: значительное ускорение по сравнению с универсальными решателями, такими как QuTiP
3. Обеспечение полной интеграции программного экосистема:
   • Полная совместимость с пакетом Pulser
   • Поддержка рабочих процессов машинного обучения PyTorch
   • Реализация унифицированного интерфейса API
4. Разработка практического руководства по выбору: эмпирические правила для выбора эмулятора в различных сценариях применения

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

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

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

• Система из N кубитов на нейтральных атомах
• Зависящие от времени параметры управления Ωᵢ(t) (частоты Раби) и δᵢ(t) (расстройка)
• Начальное квантовое состояние |ψ(0)⟩

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

• Квантовое состояние после временной эволюции |ψ(t)⟩
• Ожидаемые значения наблюдаемых

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

• Ограничения по памяти (GPU/CPU)
• Требования к точности
• Ограничения по времени вычисления

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

Архитектура Emu-sv

1. Представление состояния: использование полного вектора состояния размерности 2^N
2. Хранение гамильтониана: использование структурных особенностей гамильтониана Ридберга, хранение только диагональных элементов и значений Ωᵢ
3. Алгоритм временной эволюции:
   • Использование кусочно-постоянного приближения: |ψ(tₖ₊₁)⟩ = e^(-iHₖdt)|ψ(tₖ)⟩
   • Применение алгоритма Ланцоша для вычисления матричной экспоненты
   • Формально представляет собой интегратор второго порядка

Архитектура Emu-mps

1. Представление MPS: представление квантового состояния как матричного произведения состояний $c_{i_1...i_N} = A^{i_1}_{j_1} A^{i_2}_{j_1,j_2} \ldots A^{i_N}_{j_N}$
2. Алгоритм временной эволюции:
   • Использование зависящего от времени вариационного принципа (TDVP)
   • Двухсайтовая схема TDVP второго порядка
   • Обработка моделей шума через квантовые скачки Монте-Карло
3. Оптимизация памяти: контроль использования памяти через размерность связи, достижение экспоненциального сбережения памяти для состояний с низкой запутанностью

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

1. Структурированное хранение гамильтониана:
   • Использование того факта, что σˣᵢ члены производят только недиагональные элементы
   • Реализация матрично-векторного умножения через обмен подвекторов
   • Экономия памяти по сравнению с форматами COO и CSR
2. Адаптивный порядок интегрирования:
   • Алгоритм Ланцоша достигает адаптивного порядка через построение подпространства Крылова
   • Избегает фиксированных верхних пределов методов, таких как ZVODE
   • Лучшая производительность на больших системах с широким спектром гамильтониана
3. Интеллектуальное упорядочение кубитов:
   • Использование алгоритма Катхилла-Макки для оптимизации порядка кубитов в MPS
   • Снижение ошибки алгоритма TDVP в двумерных системах
   • Уменьшение требуемой размерности связи

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

Аппаратное окружение

• GPU: NVIDIA A100 40GB DRAM
• CPU: AMD EPYC 7742, до 1TB RAM
• Версии программного обеспечения:
  • Pulser 1.5 (на основе QuTiP 5.2, NumPy 2.2, SciPy 1.15)
  • Версия эмулятора 2.3 (на основе PyTorch 2.8)

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

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

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

• Стандартный бэкенд Pulser: решатель ODE ZVODE на основе QuTiP
• Различные конфигурации: сравнение с различными временными шагами dt и параметрами точности

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

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

Производительность Emu-sv

1. Сравнение времени выполнения:
   • Начиная с 9 кубитов, emu-sv значительно превосходит Pulser
   • Экспоненциальное масштабирование проявляется позже, но в конечном итоге ограничено памятью
   • Ускорение на GPU значительно на больших системах
2. Потребление памяти:
   • Для 26 кубитов требуется примерно 1GB памяти для хранения вектора состояния
   • При предположении 15 векторов Крылова, общее потребление памяти для 26 кубитов составляет менее 20GB
3. Анализ точности:
   • Основной источник ошибки — схема временной дискретизации
   • При dt=5ns типичная ошибка составляет примерно 10⁻⁵
   • Параметр точности p в основном влияет на сходимость Ланцоша, а не на общую ошибку

Производительность Emu-mps

1. Предсказание памяти:
   • Возможно предварительно оценить требования к памяти для заданного числа кубитов N и максимальной размерности связи χ
   • Масштаб системы, обрабатываемой 40GB GPU, зависит от размерности связи
2. Оценка времени выполнения:
   • Время одного временного шага зависит от N и χ
   • Зависимость числа итераций Ланцоша от диапазона параметров относительно слабая
3. Контроль точности:
   • Точность усечения по умолчанию p=10⁻⁵
   • Граница общей ошибки: p × m × (2N-3), где m — число временных шагов
   • Фактическая ошибка обычно меньше теоретической границы

Руководство по выбору

На основе тестов GPU эмпирические правила:

• ≤27 кубитов: рекомендуется использовать emu-sv
• >27 кубитов: рекомендуется использовать emu-mps

Примеры применения

1. Вычисление динамического структурного фактора: emu-mps успешно применён к системам с более чем 20 кубитами
2. Квантовая оптимизация графов: оба эмулятора применены к задачам квантовой оптимизации графов

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

Программное обеспечение для квантового моделирования

• QuTiP: универсальный квантовый инструментарий, использующий решатель ZVODE для зависящих от времени задач
• Pulser: специализированный пакет Python с открытым исходным кодом для квантовых вычислений на нейтральных атомах
• Другие платформы: инструменты моделирования для различных аппаратных платформ (сверхпроводящие кубиты, ионные ловушки, фотонные устройства)

Методы тензорной сети

• MPS/DMRG: группа перенормировки матрицы плотности и методы матричного произведения состояний
• TDVP: зависящий от времени вариационный принцип в квантовой динамике
• Квантовый Монте-Карло: стохастические методы для обработки динамики открытых систем

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

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

1. Значительное повышение производительности: оба эмулятора обеспечивают значительное улучшение производительности по сравнению с QuTiP
2. Дополняющий дизайн: emu-sv подходит для точных вычислений на малых системах, emu-mps — для приближённых вычислений на больших системах
3. Удобство использования: унифицированный дизайн API упрощает переключение между различными эмуляторами

Ограничения

1. Ограничения памяти Emu-sv: ограничено экспоненциальным ростом вектора состояния, максимум примерно 27 кубитов
2. Контроль точности Emu-mps: требует правильной конфигурации параметров усечения, проблемы упорядочения в двумерных системах
3. Порядок интегрирования: в настоящее время используется интегратор второго порядка, не применяются более высокие методы, такие как экстраполяция Ричардсона

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

1. Дифференцируемость: emu-sv уже поддерживает дифференцируемые вычисления, планируется добавить эту функцию в emu-mps
2. Отслеживание основного состояния: реализация алгоритма DMRG для отслеживания основного состояния при временной эволюции
3. Расширение моделей шума: добавление новых источников шума на основе анализа аппаратной команды
4. Интеграция машинного обучения: поддержка рабочих процессов машинного обучения на основе PyTorch

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

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

1. Высокая практическая ценность: непосредственно решает реальные потребности в области квантовых вычислений на нейтральных атомах
2. Отличная техническая реализация:
   • Полное использование структурных особенностей гамильтониана Ридберга
   • Адаптивные характеристики алгоритма Ланцоша хорошо работают на больших системах
   • Значительная оптимизация памяти метода MPS для состояний с низкой запутанностью
3. Качество программной инженерии:
   • Хорошая интеграция с существующей экосистемой
   • Удобный унифицированный дизайн API
   • Предоставлено практическое руководство по выбору
4. Полное тестирование: подробные сравнения по нескольким измерениям (время выполнения, память, точность)

Недостатки

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

Влияние

1. Вклад в область: предоставляет важный инструмент моделирования для квантовых вычислений на нейтральных атомах
2. Практическая ценность: уже применён и проверен в нескольких исследовательских проектах
3. Вклад в открытый исходный код: как часть экосистемы Pulser способствует развитию области

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

1. Разработка квантовых алгоритмов: проверка алгоритмов перед реализацией на оборудовании
2. Тестирование квантового оборудования: оценка производительности квантового оборудования
3. Исследование квантовой физики: изучение динамики многотельных квантовых систем
4. Образование и подготовка: обучение и исследовательская подготовка в области квантовых вычислений

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

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

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