Benchmarking Open-Access Quantum Hardware Through Many-Body Dynamics in the 1D Tilted-Field Ising Model

Resumen Técnico: Evaluación de Rendimiento de Hardware Cuántico de Acceso Abierto mediante la Dinámica de Muchos Cuerpos en el Modelo de Ising de Campo Inclinado 1D

Planteamiento del Problema
La simulación de Hamiltonianos es una aplicación principal propuesta para las computadoras cuánticas de la era actual (near-term), pero su utilidad práctica depende de si el hardware físico puede preservar la dinámica de muchos cuerpos pretendida frente a los límites de muestreo, el ruido de las puertas, la decoherencia y las restricciones de profundidad de circuito. Mientras que el modelo de Ising de campo transversal es un sistema integrable estándar, la adición de un campo longitudinal crea el modelo de Ising de campo inclinado, el cual rompe la integrabilidad e introduce dinámicas complejas y no triviales, incluyendo el caos cuántico y la termalización. El problema central abordado es determinar hasta qué punto el hardware cuántico de acceso abierto actual puede reproducir estos observables de muchos cuerpos físicamente significativos en comparación con simulaciones ideales sin ruido.

Metodología
El estudio utiliza una cadena de Ising de campo inclinado unidimensional de 15 cúbits como sistema de referencia, ejecutada en el backend de IBM Quantum ibm_marrakesh.

• Inicialización del Sistema: El sistema se inicializa con todos los cúbits en el estado de espín hacia arriba ($|0\rangle$), excepto el cúbit central ($q8$), que se invierte a espín hacia abajo ($|1\rangle$) mediante una puerta X para crear una excitación localizada.
• Hamiltoniano y Parámetros: La dinámica está gobernada por el Hamiltoniano $H = -J \sum Z_i Z_{i+1} - h_x \sum X_i - h_z \sum Z_i$. El estudio compara un caso de referencia integrable ($h_z = 0$) contra un caso de campo inclinado no integrable ($h_z = h_x = 0.35625$) a través de tres regímenes de interacción: $J < h_x$, $J = h_x$, y $J > h_x$.
• Técnica de Simulación: La evolución temporal se aproxima mediante una descomposición de Trotter simétrica de segundo orden con un paso de tiempo $\Delta t = 0.1$. Esta elección equilibra la precisión (el error local escala como $O(\Delta t^3)$) con la profundidad del circuito para que sea factible en hardware de acceso abierto.
• Observables: Se rastrean tres observables clave para evaluar la dinámica:
  1. Magnetización Local: $\langle Z_i(t) \rangle$ para rastrear la orientación del espín individual.
  2. Magnetización Total: $M_z(t) = \sum \langle Z_i(t) \rangle$ para medir la memoria de la polarización global del espín.
  3. Correlaciones de Espín Conectadas: $C_{ij}^{conn}(t) = \langle Z_i Z_j \rangle - \langle Z_i \rangle \langle Z_j \rangle$ para analizar la propagación de correlaciones y las estructuras de cono de luz.
• Ejecución Experimental: Los circuitos se ejecutaron en ibm_marrakesh con 2048 disparos (shots) por circuito. El estudio empleó el nivel de resiliencia 1 con desacoplamiento dinámico (secuencia XpXm), pero no utilizó extrapolación de ruido cero implementada por el usuario ni corrección de errores cuánticos completa. Los resultados del hardware fueron muestreados cada 0.5 unidades de tiempo hasta $t=35$.

Resultados Clave
El estudio compara las simulaciones de Trotter clásicas ideales con los resultados de la ejecución en hardware:

• Magnetización Local: El hardware reproduce la tendencia general de decaimiento de la magnetización en los tiempos tempranos. Sin embargo, los valores del hardware están "comprimidos", mostrando una reducción en la separación entre cúbits en comparación con la simulación ideal. En los casos de campo inclinado ($h_z = h_x$), la simulación ideal mantiene valores de magnetización local más altos y una estructura oscilatoria, mientras que los datos del hardware muestran un colapso más rápido hacia cero y una reducción en la variación sitio a sitio.
• Magnetización Total: De manera similar a la magnetización local, los datos del hardware muestran una caída inicial rápida seguida de un estado estacionario más bajo en comparación con la simulación ideal. En el régimen de campo inclinado, la simulación ideal mantiene una magnetización positiva significativa (memoria del estado inicial), mientras que los datos del hardware muestran una reducción mucho más fuerte, indicando una pérdida de la memoria de la magnetización global.
• Correlaciones Conectadas: La simulación ideal exhibe una estructura clara de "cono de luz" donde las correlaciones se propagan hacia afuera desde la excitación central ($q8$) de manera estructurada. En contraste, los resultados del hardware muestran un patrón significativamente más débil y ruidoso. El frente de propagación distintivo queda en gran medida oscurecido, reemplazado por fluctuaciones dispersas de baja intensidad, particularmente en tiempos más largos y distancias mayores desde el centro.

Contribuciones Clave

• Protocolo de Evaluación de Rendimiento: El artículo establece un benchmark reproducible de nivel de observable para la fidelidad de la propagación de correlaciones en una cadena de espines no integrable en hardware IBM de clase Heron.
• Comparación Cuantitativa: Proporciona una comparación directa y cuantitativa de la dinámica ideal frente a la ruidosa para un sistema de 15 cúbits, aislando específicamente los efectos de romper la integrabilidad en el rendimiento del hardware.
• Caracterización del Ruido: El estudio demuestra que, si bien el hardware de acceso abierto puede capturar características gruesas de la dinámica temprana pretendida, la propagación de correlaciones a largo plazo y las estructuras oscilatorias detalladas están cada vez más oscurecidas por el muestreo finito, los errores de puertas de dos cúbits, los errores de lectura y la decoherencia.

Significancia y Afirmaciones
El artículo afirma modestamente que el estudio sirve como un benchmark de corto plazo en lugar de una simulación totalmente corregida de errores. Su principal significancia radica en demostrar los límites actuales del hardware cuántico de acceso abierto para preservar el comportamiento de muchos cuerpos físicamente significativo.

• Preservación de la Dinámica: El hardware preserva algunas características de tiempo temprano (tendencias generales de magnetización) pero no logra preservar de manera confiable la propagación de correlaciones a largo plazo.
• Limitaciones: Los resultados resaltan que, sin mitigación o corrección de errores, es difícil extraer la estructura detallada de la dinámica no integrable (como patrones de correlación específicos y memoria oscilatoria) de los datos ruidosos.
• Direcciones Futuras: Los autores sugieren que, si bien la mitigación de errores (por ejemplo, mitigación de error de lectura, extrapolación de ruido cero) podría mejorar la precisión, la corrección de errores cuánticos completa sigue siendo una carga significativa para los dispositivos actuales. Proponen como trabajo futuro las comparaciones entre diferentes backends, cadenas de espines más pequeñas para aislar el ruido, y la inclusión de observables más complejos como los correladores fuera del tiempo (OTOCs) y la entropía de entrelazamiento para investigar más a fondo el caos cuántico y la propagación de la información.