Fast design and scaling of multi-qubit gates in large-scale trapped-ion quantum computers

Resumen Técnico: Diseño Rápido y Escalamiento de Puertas Multiqubit en Computadoras Cuánticas de Iones Atrapados a Gran Escala

Declaración del Problema
Las computadoras cuánticas de iones atrapados ofrecen conectividad de largo alcance y tiempos de coherencia prolongados, lo que las convierte en una plataforma líder para la computación cuántica. Sin embargo, escalar estos sistemas a cientos de iones presenta un desafío conceptual y tecnológico significativo. Aunque los iones poseen conectividad de todos contra todos a través de sus modos de movimiento colectivos (bus de fonones), diseñar puertas de entrelazamiento multiqubit programables y de alta fidelidad para cristales grandes es computacionalmente prohibitivo.

La dificultad central reside en el problema del diseño de la puerta: encontrar un conjunto de impulsos de control que satisfagan simultáneamente $N(N-1)/2$ restricciones cuadráticas (fases de entrelazamiento objetivo entre todos los pares de qubits) mientras se minimiza la potencia del impulso y se asegura que los modos de movimiento regresen al estado fundamental. Esto se formula como un problema de optimización con restricciones cuadráticas que escala como $N^2$ con el número de iones $N$, lo que lo convierte en un problema NP-duro. Los enfoques existentes suelen restringir el problema a subconjuntos pequeños de iones o dependen de algoritmos de minimización basados en gradientes que tienen dificultades con la complejidad y las restricciones de sistemas a gran escala.

Metodología: El Método de Gran Escala Rápido (LSF)
Los autores introducen un marco metodológico denominado "Gran Escala Rápido" (LSF, por sus siglas en inglés) para mitigar la naturaleza NP-dura del problema de diseño de puertas, permitiendo el diseño de puertas robustas y programables para cristales de hasta cientos de iones.

1. Formulación en el Dominio Espectral: El diseño de la puerta se formula en el dominio espectral, representando los campos de impulso como una suma de componentes de Fourier desintonizados alrededor de la frecuencia de transición del qubit. El problema se reduce a encontrar un vector $R$ (amplitudes de los tonos de impulso) que minimice la norma $\|R\|$ sujeto a restricciones lineales (desplazamientos de movimiento nulos) y restricciones cuadráticas (fases de entrelazamiento objetivo $\phi_{nn'}$).
2. Semilla de Solución de Fase Cero (ZPS): El método LSF se basa en la generación de una solución "semilla" conocida como Solución de Fase Cero (ZPS). Una ZPS es una solución no trivial al problema homogéneo donde todas las fases de entrelazamiento objetivo son cero ($\phi^{(0)}_s = 0$). Estas soluciones se acumulan antes de resolver para puertas específicas.
3. Conversión y Optimización: El algoritmo procede en dos pasos principales:
   • Conversión: Una ZPS se escala por un factor grande $\lambda$ y se perturba con un vector pequeño $D$ para satisfacer las restricciones cuadráticas objetivo. Al linealizar las restricciones alrededor de la ZPS escalada, los autores derivan un sistema lineal para resolver $D$, convirtiendo la ZPS en una solución aproximada $R_{conv}$ que satisface las fases objetivo.
   • Reducción de Norma: Un proceso iterativo minimiza la norma del vector de solución mientras se adhiere estrictamente a las restricciones. Esto implica alternar entre pasos de "reducción de norma" (moverse a lo largo de la variedad de la solución para reducir la potencia) y pasos de "reducción de error" (corregir desviaciones de las fases objetivo). Los tamaños de paso se calculan analíticamente utilizando métodos de búsqueda de raíces polinómicas (métodos de Cardano y Ferrari) para asegurar que el error permanezca dentro de una tolerancia especificada (típicamente $10^{-4}$).
4. Agregación: El método agrega múltiples ZPSs para asegurar un espacio de búsqueda diverso, permitiendo que el algoritmo encuentre soluciones eficientes para diversos mapas de entrelazamiento.

Contribuciones Clave y Resultados

• Escalamiento de Tiempo Polinomial: El método LSF transforma el problema de diseño NP-duro en una tarea de generar ZPSs (que pueden hacerse eficientemente para simetrías específicas) seguida de una optimización eficiente de tiempo polinomial. Las pruebas de referencia muestran que LSF logra un tiempo de solución aproximadamente un orden de magnitud más rápido que los optimizadores estándar de región de confianza y de gradiente conjugado, particularmente para sistemas con $N > 20$ iones.
• Estimación de la Potencia de Impulso: Los autores derivan un estimador heurístico para la frecuencia de Rabi total ($\Omega_{nuc}$) requerida para impulsar una puerta. Este estimador depende de la norma nuclear de la matriz de entrelazamiento objetivo $\|\phi\|_{nuc}$ y escala linealmente con el número de iones $N$. Sorprendentemente, los resultados de la simulación muestran que la potencia de impulso requerida es ampliamente independiente de la estructura de modos específica o el espaciamiento de frecuencia, colapsando en una curva universal cuando se normaliza por la escala de tiempo de la densidad de modos promedio $T_\nu$.
• Tiempo de Puerta Mínimo ($T_{min}$): El estudio identifica un límite inferior fundamental en la duración de la puerta. Por debajo de un tiempo crítico $T_{min} \approx 2.9 \cdot T_\nu$ (donde $T_\nu$ es proporcional al inverso del espaciamiento de modo promedio), no se pueden encontrar soluciones de alta fidelidad. Crucialmente, $T_{min}$ escala linealmente con el número de iones $N$. Esto implica que, si bien el número de operaciones de entrelazamiento simultáneas escala como $N^2$, la duración de la puerta solo escala como $N$, proporcionando una ventaja de escalamiento sobre las implementaciones de puertas secuenciales.
• Análisis de Escalamiento de Errores: El artículo cuantifica el impacto de tres fuentes de error dominantes en la fidelidad de la puerta:
  • Derivas de la frecuencia de los modos de movimiento: Los errores escalan cuadráticamente con la deriva de frecuencia y el número de acoplamientos ($N^2$).
  • Fluctuaciones de la amplitud del impulso: Los errores escalan con la varianza del ruido y el número de acoplamientos.
  • Calentamiento de los modos de movimiento: El escalamiento de los errores inducidos por el calentamiento depende de la correlación espacial del ruido. Para ruido no correlacionado, la fidelidad decae cuadráticamente con $N$; para ruido fuertemente correlacionado, decae linealmente.
    El marco LSF permite la inclusión de restricciones lineales para robustecer las puertas contra estos errores, reduciendo significativamente la sensibilidad a las derivas y al calentamiento.
• Ejemplo de Código de Superficie: El método se demuestra en un cristal de 49 iones para implementar una puerta de estabilizador de código de superficie (entrelazando 33 iones en un solo pulso). La solución muestra una distribución de potencia eficiente, asignando amplitudes de impulso más altas a los iones "nodo" y amplitudes más bajas a los iones "borde", y suprime eficazmente el entrelazamiento cruzado no deseado entre pares desconectados.

Significado y Reivindicaciones
El artículo afirma delinear un camino hacia el escalamiento de las computadoras cuánticas de iones atrapados a cientos de qubits al superar la barrera conceptual de diseñar puertas multiqubit. Al proporcionar un algoritmo rápido y robusto (LSF), los autores permiten la compilación offline de circuitos que involucran decenas de iones y puertas multiqubit complejas, una tarea anteriormente considerada formidable debido a la naturaleza NP-dura de la optimización.

El trabajo establece que las puertas multiqubit de alta fidelidad son físicamente factibles para cristales grandes, siempre que el tiempo de la puerta exceda el umbral de escalamiento lineal $T_{min}$. Además, las relaciones derivadas entre la potencia de impulso, el tiempo de puerta y las fuentes de error proporcionan directrices de diseño prácticas y tolerancias para los experimentales. Los autores señalan que, si bien su método maneja los mecanismos de error dominantes dentro de la aproximación de Lamb-Dicke, las interacciones con los portadores de qubits y los términos de orden superior en la expansión de Magnus para soluciones de alta amplitud requieren mayor verificación. El marco se presenta como una herramienta fundacional para realizar el potencial de la conectividad de largo alcance en arquitecturas de iones atrapados para aplicaciones de corrección de errores cuánticos y optimización.