MPStab: an hybrid stabilizers tensor-network quantum circuit simulator

Resumen Técnico de "MPStab: un simulador de circuitos cuánticos de redes de tensores y estabilizadores híbridos"

Planteamiento del Problema
La simulación de sistemas cuánticos en computadoras clásicas enfrenta un cuello de botella fundamental: el crecimiento exponencial de la dimensión del espacio de estados con el número de cúbits. Aunque las simulaciones de vector de estado son exactas, se vuelven intratables más allá de unas pocas docenas de cúbits. Los algoritmos estructurados existentes intentan sortear esto explotando propiedades específicas del circuito, pero enfrentan limitaciones distintas:

• Redes de Tensores (p. ej., MPS): Eficientes para sistemas con entrelazamiento acotado (interacciones locales, dinámica de tiempo corto), pero fallan cuando el entrelazamiento crece rápidamente, requiriendo un aumento exponencial de la dimensión de enlace (bond dimension).
• Métodos de Estabilizadores: Eficientes para circuitos de Clifford (tiempo polinomial) mediante el seguimiento de estabilizadores de Pauli, pero carecen de universalidad; no pueden manejar eficientemente puertas no-Clifford (p. ej., puertas $T$, rotaciones arbitrarias), las cuales introducen "no-estabilizabilidad" (magic).

El desafío es simular circuitos que contienen tanto un entrelazamiento significativo (que requiere redes de tensores) como operaciones no-Clifford (que requieren el seguimiento de estabilizadores), regímenes donde ningún método por sí solo es eficiente.

Metodología: El Formalismo HSMPO
El artículo presenta mpstab, una biblioteca de Python que implementa el formalismo de Operador de Producto Matricial de Estabilizador-Híbrido (HSMPO). La metodología central se basa en la descomposición de un circuito cuántico universal $U$ en bloques alternos de unitarias de Clifford ($C_j$) y rotaciones de un solo cúbit de tipo "magia" ($R_j$):
$$U = C_0 R_1 C_1 R_2 \dots R_N C_N$$

El algoritmo explota la propiedad de que la conjugación de Clifford mapea cadenas de Pauli en cadenas de Pauli. Al "empujar" las puertas de Clifford a través de las rotaciones de magia mediante conjugación ($C_j R_j C_j^\dagger = \tilde{R}_j$), el circuito se reescribe como una secuencia de rotaciones "vestidas" $\tilde{R}_j$ seguidas de una evolución de Clifford residual.

• Rotaciones Vestidas: Cada $\tilde{R}_j$ es una rotación alrededor de una cadena de Pauli (potencialmente no local). Crucialmente, cualquier rotación de este tipo admite una representación exacta de Operador de Producto Matricial (MPO) con una dimensión de enlace fija de 2, independientemente de la no-localidad de la cadena de Pauli.
• Evolución Híbrida:
  1. Retropropagación (Imagen de Heisenberg): Las operaciones de Clifford se aplican a los observables utilizando aritmética de tablas (tableau) eficiente (formalismo de estabilizadores). Esto es exacto y escala polinomialmente.
  2. Evolución hacia Adelante (Imagen de Schrödinger): Las rotaciones vestidas $\tilde{R}_j$ se aplican al estado cuántico (representado como un MPS) como MPOs. El estado se comprime después de cada paso para mantener una dimensión de enlace máxima $\chi$.
• Valores de Expectación: El valor de expectación final se calcula mediante la contracción del MPS evolucionado con el observable retropropagado.

Contribuciones Clave

1. Biblioteca mpstab: Los autores presentan un paquete de Python de código abierto y modular que implementa el marco de trabajo HSMPO. Se integra con el ecosistema Qibo, permitiendo a los usuarios definir circuitos usando sintaxis estándar mientras aprovechan la simulación híbrida.
2. Arquitectura Modular: El paquete separa la simulación en motores independientes:
   • Motores de Estabilizadores: Maneja la retropropagación de Clifford (p. ej., vía Stim).
   • Motores de Redes de Tensores: Maneja la evolución MPS/MPO (p. ej., vía Quimb con backend de Jax).
     Esta modularidad permite a los usuarios intercambiar backends, habilitar la diferenciación automática o integrar motores personalizados.
3. Funcionalidades de Mitigación de Errores: La biblioteca incluye funcionalidad para generar circuitos "sustitutos" reemplazando probabilísticamente puertas de magia por puertas de Clifford, una característica útil para técnicas de mitigación de errores basadas en datos, como la Regresión de Datos de Clifford.
4. Análisis de Complejidad: Los autores derivan un modelo de costo computacional comparando HSMPO con un MPS estándar. Demuestran que HSMPO ofrece una ventaja de tiempo polinomial para circuitos dominados por puertas de Clifford con contenido no-Clifford moderado, siempre que la dimensión de enlace requerida para una fidelidad objetivo sea significativamente menor en el enfoque híbrido que en un MPS puro.

Resultados y Benchmarks
El artículo evalúa el rendimiento de mpstab frente a un simulador de red de tensores puro (Quimb) utilizando un ansatz personalizado que controla independientemente la cantidad de entrelazamiento (vía profundidad de la escalera CNOT) y el contenido de "magia" (vía la fracción de rotaciones no-Clifford, $\beta$).

• Validación de Costo: Las mediciones de tiempo de ejecución (wall-clock) confirman el modelo de costo analítico. Existe un punto de cruce donde mpstab se vuelve más eficiente que un MPS puro; esto ocurre cuando el circuito tiene una alta proporción de puertas de Clifford respecto a las no-Clifford y la dimensión de enlace requerida es grande.
• Ventaja de Fidelidad: A una dimensión de enlace $\chi$ fija, mpstab logra consistentemente una fidelidad mayor que un MPS puro, particularmente a medida que aumenta la profundidad del circuito. El componente de estabilizador absorbe exactamente el entrelazamiento generado por las puertas de Clifford, dejando que el MPS represente únicamente el entrelazamiento proveniente de las rotaciones no-Clifford vestidas.
• Escalabilidad: La ventaja de mpstab es más pronunciada en regímenes con recuentos bajos de puertas no-Clifford ($N_{nC}$) y altos de puertas de Clifford ($N_C$). A medida que el tamaño del sistema aumenta, la ventaja persiste solo si la retención de la magia $\beta$ se ajusta de tal manera que el número total de puertas de magia no escale de forma prohibitiva.

Significancia
El artículo posiciona a mpstab como una herramienta práctica para extender los límites de la simulación cuántica clásica. No pretende simular computadoras cuánticas de gran escala arbitrarias, sino que apunta a regímenes específicos donde la computación clásica sigue siendo viable: circuitos con una estructura de Clifford sustancial e inyección moderada de elementos no-Clifford.

La significancia radica en:

1. Cerrar la Brecha: Combina efectivamente la eficiencia de los métodos de estabilizadores para la generación de entrelazamiento con la universalidad de las redes de tensores para las operaciones no-Clifford.
2. Utilidad Práctica: Al integrarse con Qibo y ofrecer motores modulares, facilita el uso de la simulación híbrida en algoritmos variacionales y flujos de trabajo de mitigación de errores.
3. Identificación de Regímenes: Los autores proporcionan un "mapa de fases" concreto (basado en el recuento de puertas de Clifford vs. no-Clifford y dimensiones de enlace) que guía a los usuarios sobre cuándo emplear mpstab sobre otros simuladores estándar.

El trabajo concluye que mpstab es el simulador preferido para circuitos dominados por Clifford con contenido no-Clifford moderado, donde logra una mayor precisión con un menor costo computacional en comparación con los enfoques de redes de tensores puros.