Digital Quantum Simulation of Nonequilibrium Dynamics in the Schwinger Model under a Strong External Electric Field

Resumen Técnico: Simulación Cuántica Digital de la Dinámica Fuera del Equilibrio en el Modelo de Schwinger bajo un Campo Eléctrico Externo Fuerte

Planteamiento del Problema
El artículo aborda el desafío de simular la dinámica fuera del equilibrio en la electrodinámica cuántica (QED) de campo fuerte, específicamente el efecto Schwinger, donde un campo eléctrico fuerte induce la decadencia del vacío y la producción de pares partícula-antipartícula. Los métodos perturbativos tradicionales fallan al describir este fenómeno no perturbativo en tiempo real de manera fiable. Además, los métodos estándar de Monte Carlo de integral de camino euclidiano encuentran el "problema del signo" al tratar con la evolución en tiempo real o pesos complejos, lo que conduce a un crecimiento rápido de los errores estadísticos. Si bien el formalismo hamiltoniano evita el problema del signo, requiere el manejo de un espacio de Hilbert exponencialmente grande. Los autores proponen el uso de la simulación cuántica digital para superar estas limitaciones mediante la representación directa de estados cuánticos de muchos cuerpos y la evolución unitaria en el tiempo sobre cúbits. El enfoque específico es el modelo de Schwinger en (1+1) dimensiones bajo un campo eléctrico externo constante, investigando la dinámica fuera del equilibrio de tamaño finito tras un quench (cambio brusco) de campo eléctrico.

Metodología
El estudio emplea un protocolo de simulación cuántica digital basado en el modelo de Schwinger en red de (1+1) dimensiones con condiciones de contorno abiertas. La metodología consta de tres etapas principales:

1. Formulación del Modelo y Mapeo a Cúbits:

   • La teoría continua se discretiza utilizando el formalismo de fermiones escalonados de Kogut-Susskind.
   • Se utiliza la ley de Gauss para eliminar los grados de libertad de calibre, expresando el campo eléctrico en cada enlace como una función de la distribución de carga de materia y del campo de frontera. Esto reduce el modelo a un Hamiltoniano puramente fermiónico con interacciones de largo alcance.
   • Los operadores fermiónicos se mapean a operadores de Pauli mediante la transformación de Jordan-Wigner, resultando en un Hamiltoniano de cúbits compuesto por términos de energía cinética, de masa y de campo eléctrico.

2. Preparación del Estado (VQE):

   • El estado inicial para la dinámica es el estado de vacío de campo cero ($\epsilon=0$).
   • Este estado se prepara utilizando el Algoritmo Cuántico Variacional de Autovalores (VQE) con un ansatz de RealAmplitudes eficiente para el hardware y el optimizador clásico SLSQP.
   • Para estudiar las propiedades estáticas del sistema bajo campos externos variables, los autores también utilizan VQE para encontrar el estado fundamental y la Deflación Cuántica Variacional (VQD) para estimar la energía del primer estado excitado.

3. Evolución en Tiempo Real (Trotter-Suzuki):

   • El sistema se somete a un "quench" repentino donde se activa un campo eléctrico externo en $t=0$.
   • La evolución en tiempo real se simula utilizando una descomposición de Trotter-Suzuki de segundo orden para aproximar el operador de evolución temporal $e^{-iHt}$.
   • La simulación se realiza en computadoras clásicas utilizando el marco de trabajo Qiskit, sin incluir ruido cuántico, para comparar el algoritmo con la Diagonalización Exacta (ED).

Contribuciones Clave y Resultados
El artículo presenta resultados numéricos para tamaños de red $N=8$ y $N=12$ (con unidades de red $a=m=g=1$), comparando el protocolo de simulación cuántica contra los puntos de referencia de la Diagonalización Exacta (ED).

• Propiedades Estáticas del Estado Fundamental:

  • El estado fundamental de campo cero preparado mediante VQE logra una fidelidad alta ($>99.99\%$) en comparación con los resultados de ED.
  • Los escaneos sobre la intensidad del campo externo $\epsilon$ revelan un comportamiento de tipo escalón en el condensado quiral y una brecha de energía que desaparece entre el estado fundamental y el primer estado excitado en intensidades de campo críticas específicas. Esto indica una reconstrucción del campo inducida por el campo de la estructura del vacío, que corresponde a una transición de fase de primer orden en volúmenes finitos.
  • La intensidad de campo crítica $\epsilon_c$ exhibe una clara dependencia de tamaño finito sobre $1/N$. La extrapolación lineal arroja un campo crítico de volumen infinito de $\epsilon_c(\infty) \approx 0.469$, el cual los autores atribuyen a la influencia de las condiciones de contorno abiertas al reducir la brecha de energía.

• Dinámica Fuera del Equilibrio:

  • Leyes de Conservación: La evolución de Trotter de segundo orden preserva con éxito la carga $U(1)$ total y la energía total (dentro de la precisión numérica), validando la estabilidad del algoritmo.
  • Separación de Carga: Tras activar el campo, la distribución espacial de carga muestra la producción y aniquilación periódica de pares carga-anticarga cerca de las fronteras abiertas. La amplitud y el periodo de estas oscilaciones aumentan con la intensidad del campo externo.
  • Redistribución de Energía: La energía total del campo eléctrico $H_E(t)$ exhibe un comportamiento cuasiperiódico, reflejando el intercambio de energía entre el término del campo eléctrico y el sector fermiónico. Esto es distinto de la periodicidad más simple de la distribución de carga, lo que sugiere contribuciones complejas de modos de mayor energía.
  • Fidelidad del Vacío: La fidelidad del estado de vacío $P_{vac}(t)$ decae a medida que el sistema es impulsado lejos del vacío inicial. En la ventana de tiempo temprano ($0 \le t/a \le 1.0$), esta decadencia sigue una ley exponencial efectiva, $P_{vac}(t) \sim e^{-\gamma_{eff} t}$. La tasa de decadencia efectiva $\gamma_{eff}$ aumenta monótonamente con la intensidad del campo externo.
  • Efectos de Tamaño Finito: Las comparaciones entre los sistemas $N=8$ y $N=12$ confirman que las características dinámicas observadas (separación de carga en la frontera, decadencia de la fidelidad del vacío) no son artefactos numéricos de tamaños de sistema pequeños, sino que persisten en redes más grandes, donde los sistemas mayores muestran una desviación más rápida del vacío inicial.

Significancia y Reivindicaciones
Los autores afirman que su trabajo demuestra la viabilidad de combinar la preparación de estados basada en VQE con la evolución digital en tiempo real (vía descomposición de Trotter-Suzuki) para estudiar la dinámica fuera del equilibrio en teorías de calibre en red de campo fuerte.

• Validación: El protocolo reproduce con éxito las principales características de los fenómenos inducidos por el campo —específicamente la separación de carga en la frontera, la decadencia de la fidelidad del estado de vacío y la redistribución cuasiperiódica de la energía— cuando se compara con la diagonalización exacta.
• Utilidad Metodológica: El estudio sirve como un punto de referencia para la simulación cuántica digital en física de altas energías, mostrando que la evolución digital explícita de Trotter puede capturar eficazmente dinámicas complejas de campo fuerte en sistemas de red finitos.
• Perspectiva Futura: Los autores sugieren modestamente que este enfoque proporciona una base para futuras extensiones a tamaños de sistema más grandes, ansatzes variacionales mejorados, implementación en hardware cuántico ruidoso y el estudio de campos externos dependientes del tiempo. No pretenden una realización experimental inmediata en el hardware ruidoso actual, sino que enfatizan la verificación algorítmica en simuladores clásicos como un paso necesario.