A near-term quantum simulation of the transverse field Ising model hints at Glassy Dynamics

Resumen Técnico: Una simulación cuántica a corto plazo del Modelo de Ising de Campo Transverso sugiere una dinámica vítrea

Planteamiento del Problema
El artículo aborda el desafío computacional de simular sistemas de vidrio de espín a gran escala, los cuales son críticos para comprender las fases exóticas de la materia y las transiciones vítreas en la física de muchos cuerpos cuánticos. Los enfoques computacionales clásicos, como la dinámica molecular y las simulaciones de Monte Carlo, luchan contra el crecimiento exponencial de las funciones de onda electrónicas en sistemas de muchos electrones y la alta complejidad computacional requerida para modelar sistemas desordenados durante escalas de tiempo extendidas. Específicamente, los sistemas vítreos exhiben paisajes de energía accidentados con numerosos mínimos locales y tiempos de relajación extremadamente largos, lo que los hace difíciles de simular con precisión utilizando la computación de alto rendimiento (HPC) clásica. Los autores pretenden explorar si las técnicas de simulación cuántica digital (DQS), específicamente los Solucionadores Cuánticos de Valores Propios Variacionales (VQE), pueden sondear eficazmente estas configuraciones de espín desordenadas en hardware cuántico de corto plazo.

Metodología
Los autores emplearon un enfoque híbrido cuántico-clásico utilizando el algoritmo de Solucionador Cuántico de Valores Propios Variacionales (VQE) para simular el Modelo de Ising de Campo Transverso (TFIM) con un campo longitudinal añadido. El estudio se centró en dos geometrías específicas:

1. Cadena de Ising 1D: Una cadena de longitud $n=25$ (25 cúbits).
2. Red de Ising 2D: Una red cuadrada de $6 \times 6$ (36 cúbits).

Componentes Técnicos Clave:

• Construcción del Hamiltoniano: El Hamiltoniano del sistema incluye interacciones de vecinos más cercanos a lo largo del eje z ($J \sum \sigma_i^z \sigma_{i+1}^z$), un campo transverso en la dirección x ($-B_x \sum \sigma_i^x$) y un campo longitudinal en la dirección z ($-B_z \sum \sigma_i^z$). Se excluyeron las interacciones diagonales.
• Selección del Ansatz: Los autores utilizaron un ansatz "TwoLocal" compuesto por bloques de rotación $R_x$ y $R_y$. Crucialmente, este ansatz fue diseñado sin entrelazamiento entre los cúbits. Esta elección se realizó para reducir la profundidad del circuito y mitigar la susceptibilidad al ruido en dispositivos de corto plazo, manteniendo al mismo tiempo la capacidad de capturar las fluctuaciones cuánticas necesarias para la estructura del problema.
• Optimización: Se utilizó el optimizador clásico COBYLA para minimizar el valor de la esperanza del Hamiltoniano, con un máximo de 25,000 iteraciones (aunque la mayoría de las ejecuciones convergieron antes).
• Implementación de Software: Los autores desarrollaron un paquete de Python personalizado, Qiskit Glassy Dynamics (posteriormente referido como QuGlassyIsing), para construir modularmente los Hamiltonianos para las redes de Ising 1D y 2D. Este paquete gestiona eficientemente la generación de términos de interacción ($\sigma^z\sigma^z$) y términos de campo ($\sigma^x, \sigma^z$).
• Entorno de Simulación: Las simulaciones se realizaron principalmente utilizando el backend de GPU de Qiskit Aer. También se realizó una demostración de prueba de concepto en el dispositivo real de 7 cúbits IBM Q Oslo utilizando un modelo reducido de 7 espines.

Resultados Clave

• Paisajes de Energía: Las simulaciones revelaron que aumentar el campo longitudinal ($B_z$) reduce la energía del sistema al alinear los espines, mientras que aumentar el campo transverso ($B_x$) induce fluctuaciones cuánticas y giros de espín sin reducir la energía de manera tan significativa.
• Fases Desordenadas: El estudio identificó un régimen específico donde emergen configuraciones de espín desordenadas. Esta "fase desordenada" se caracteriza por una mezcla de dominios ferromagnéticos y antiferromagnéticos.
  • En ausencia de un campo longitudinal ($B_z=0$), el sistema exhibió un orden antiferromagnético perfecto.
  • Se observó un estado desordenado solo para combinaciones específicas de $B_x$ y $B_z$ (por ejemplo, $B_z=1, B_x=2$), donde la competencia entre los campos interrumpe el orden uniforme, creando momentos locales de ferromagnetismo dentro de un fondo antiferromagnético.
  • Los autores señalan que las tres configuraciones de espín más probables en las simulaciones 2D eran casi idénticas, lo que sugiere un orden coherente y estructurado incluso dentro de la fase desordenada.
• Demostración de Hardware: La ejecución en el dispositivo IBM Q Oslo (7 cúbits) mostró un perfil de convergencia de energía con fluctuaciones significativas y un valor de energía final ($-4.0139$) superior al mínimo teórico. Este resultado se atribuye a errores de puerta, errores de lectura y la decoherencia de cúbits inherentes a los dispositivos NISQ (Quantum de Escala Intermedia con Ruido).

Significancia y Reivindicaciones
El artículo afirma proporcionar una demostración de "prueba de principio" de que los simuladores cuánticos de corto plazo pueden detectar fases desordenadas en modelos de Ising que son difíciles de acceder mediante métodos clásicos. Los autores argumentan que sus hallazgos:

• Validan la Utilidad de VQE: Demuestran que VQE, incluso con ansatzes no entrelazados, puede modelar funciones de onda complejas en tiempo polinomial, reduciendo los costos computacionales en comparación con los métodos clásicos.
• Elucidan la Dinámica Vítrea: Los resultados sugieren que la interacción entre los campos transverso y longitudinal crea un paisaje de energía "accidentado" con fases metaestables, lo que refleja el comportamiento de los vidrios de espín y los procesos catalíticos donde el desorden conduce a mecanismos de reacción ineficientes.
• Desarrollo de Herramientas: La introducción del paquete QuGlassyIsing ofrece una interfaz modular y amigable para la construcción de Hamiltonianos para la dinámica vítrea, facilitando la investigación futura en sistemas cuánticos complejos.
• Perspectiva Futura: Los autores declaran modestamente que, si bien las limitaciones actuales de NISQ (tiempo de coherencia, fidelidad de puerta) impiden la escalabilidad hacia sistemas más grandes y complejos, su trabajo sienta las bases para futuras aplicaciones de computación cuántica tolerante a fallos (FTQC). Proponen que, a medida que mejore la corrección de errores cuánticos, estos métodos podrían desbloquear la simulación de comportamientos vítreos en sistemas 3D y geometrías más complejas, informando potencialmente el diseño de nuevos materiales de almacenamiento magnético y sistemas catalíticos.

El artículo concluye que, aunque la implementación actual está limitada por las restricciones del hardware, logra resaltar el potencial de las herramientas de procesamiento de información cuántica para aportar una comprensión matizada de la dinámica y la estructura de los sistemas vítreos.