Algebraic Reduction to Improve an Optimally Bounded Quantum State Preparation Algorithm

Resumen Técnico: Reducción Algebraica para Mejorar un Algoritmo de Preparación de Estados Cuánticos con Límite Óptimo

Planteamiento del Problema
La Preparación de Estados Cuánticos (QSP, por sus siglas en inglés) es un subprocedimiento fundamental requerido para numerosos algoritmos cuánticos, incluyendo resolvedores de sistemas lineales, simulación de Hamiltonianos, química cuántica y aprendizaje automático. El desafío principal radica en minimizar la complejidad del circuito (profundidad, recuento de puertas y recuento de CNOT) del operador unitario $U$ que mapea el estado inicial $|0\rangle^{\otimes n}$ hacia un estado objetivo deseado $|\psi\rangle = \sum c_i |i\rangle$. Si bien el algoritmo de Sun et al. (SUN) es conocido por estar óptimamente limitado respecto al compromiso espacio-tiempo asintótico (utilizando cúbits ancilla para reducir la profundidad), su implementación involucra una descomposición que utiliza tres operadores de tipo $\Lambda$ para cada Puerta Controlada Uniformemente (UCG). Este artículo aborda la oportunidad de reducir la profundidad del circuito y los recuentos de puertas del algoritmo SUN mediante una descomposición algebraica más eficiente, particularmente en el rango paramétrico donde el número de cúbits ancilla $m$ satis-faga $2^n \leq m \leq 2^n/n$.

Metodología
Los autores proponen un nuevo algoritmo de QSP, denominado OSUN, que introduce una descomposición algebraica más simple del operador unitario. El cambio metodológico central consiste en separar la preparación del estado deseado en dos bloques unitarios distintos:

1. Parte Real ($U_{mod}$): Un unitario responsable de preparar los módulos reales (amplitudes) de los coeficientes del estado objetivo.
2. Parte Compleja ($D_{ph}$): Un unitario diagonal responsable de aplicar las fases complejas necesarias.

En el marco tradicional de SUN, cada UCG en la estructura de escalera se descompone en tres operadores de tipo $\Lambda$ (que manejan desplazamientos de fase diagonales para los componentes reales y complejos por separado). El algoritmo OSUN propuesto aprovecha el hecho de que $U_{mod}$ solo requiere la preparación de amplitudes reales. Al restringir las UCG en este bloque a matrices de tipo bloque-diagonal $R_y$ y utilizando un cambio específico en el conjunto de puertas ($R_y(\gamma) \equiv SH \cdot R_z(\gamma) \cdot HS^\dagger$), los autores demuestran que cada UCG en la preparación de la parte real puede implementarse utilizando un único operador de tipo $\Lambda$ en lugar de tres. La información de la fase compleja se aplica entonces al final del circuito mediante un único operador $\Lambda$ de máxima dimensión.

Esta descomposición permite la paralelización de las operaciones en el registro ancilla utilizando el mismo método de suma en el espacio de Fourier de 5 etapas establecido en las implementaciones previas de SUN, pero con un menor número de constructores $\Lambda$ por nivel.

Contribuciones Clave
El artículo describe seis contribuciones específicas:

1. Optimización de SUN: Un nuevo algoritmo de QSP que mejora al algoritmo SUN en términos de profundidad y recuento de puertas en todo el primer rango del compromiso espacio-tiempo óptimo definido en Sun et al. [40, Teorema 1].
2. Descomposición Algebraica: Una derivación que muestra cómo la separación de las partes real y compleja permite que cada UCG se implemente con un único constructor de tipo $\Lambda$, reduciendo el recuento de tres a uno por nivel de puerta.
3. Límites Asintóticos Más Estrechos: Un análisis formal que proporciona límites óptimos asintóticos más ajustados para el compromiso espacio-tiempo en el rango paramétrico especificado.
4. Implementación: La implementación completa del algoritmo OSUN utilizando la librería PennyLane, disponible en GitHub.
5. Verificación Exhaustiva: Una batería exhaustiva de pruebas ejecutadas en el sistema HPC de la Universidad de Parma, que cubre estados densos y dispersos, coeficientes aleatorios (reales positivos, reales negativos, complejos) y estados físicamente relevantes (Bell, GHZ, W, Dicke) para recuentos de cúbits de hasta 10.
6. Análisis Comparativo: Una comparación directa del rendimiento con el algoritmo de Möttönen et al. (el estándar para QSP sin ancillas) y la implementación original de SUN, evaluando métricas como profundidad, total de puertas, recuento de CNOT y tipos de puertas.

Resultados
Los resultados de la simulación indican que el algoritmo OSUN logra una reducción sistemática en la profundidad del circuito en comparación con la implementación original de SUN.

• Análisis de Complejidad: El análisis teórico confirma que, si bien la clase de complejidad asintótica permanece en $O(n^2 + 2^n/m)$, los prefactores de los términos dominantes se reducen significamente. Específicamente, el término logarítmico-lineal ($n \log m$) muestra una mejora de un factor de 3, y el término exponencial ($2^n/m$) muestra una mejora de un factor de 2.
• Métricas de Simulación: Para estados aleatorios con coeficientes complejos, OSUN demuestra una profundidad menor comenzando desde los 6 cúbits en comparación con la versión original de SUN y el algoritmo de Möttönen. La curva de crecimiento de la profundidad es más suave que la de la versión original.
• Recuento de Puertas: Aunque OSUN reduce la profundidad, el número total de puertas y CNOT no muestra una reducción en comparación con el algoritmo de Möttönen. Sin embargo, los autores señalan un cambio en la composición del conjunto de puertas: el enfoque de Möttönen utiliza CNOTs y rotaciones, mientras que OSUN utiliza CNOTs y desplazamientos de fase (operadores diagonales).
• Fidelidad: El algoritmo mantiene una fidelidad alta (cercana a 1) y una distancia de traza baja en todos los tipos de estados y recuentos de cúbits probados (2 a 10).

Significancia y Reivindicaciones
El artículo sostiene que la reducción algebraica propuesta ofrece una mejora práctica al algoritmo SUN de vanguardia sin alterar su clase de complejidad asintótica fundamental. Al reducir los factores constantes en la profundidad y el recuento de puertas, OSUN proporciona un circuito más eficiente para la preparación de estados en el régimen donde hay cúbits ancilla disponibles ($m \geq 2^n$). Los autores posicionan este trabajo como un refinamiento del compromiso espacio-tiempo óptimo, destacando que la reducción en el número de operadores $\Lambda$ (de 3 a 1 por UCG) es el motor principal de la ganancia de rendimiento.

Los autores declaran explícitamente que esta optimización se aplica a todo el primer rango paramétrico del compromiso óptimo. Sugieren que trabajos futuros podrían extender estos beneficios a otros puntos del dominio donde los registros ancilla crecen exponencialmente. El artículo concluye señalando que, si bien OSUN mejora la profundidad, investigaciones futuras deberían abordar la conectividad de cúbits en hardware real y la traducción de estos circuitos al conjunto de puertas Clifford+T para evaluar plenamente su complejidad en contextos de tolerancia a fallos. El trabajo no pretende resolver el problema general de QSP para todas las restricciones de hardware, sino que ofrece una mejora específica y verificada al marco de SUN utilizando estrategias de ancilla existentes.