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

Résumé technique : Simulation quantique numérique de la dynamique hors équilibre du modèle de Schwinger sous un champ électrique externe intense

Énoncé du problème
L'article traite du défi que représente la simulation de la dynamique hors équilibre en électrodynamique quantique (QED) à champ fort, plus précisément l'effet Schwinger, où un champ électrique intense induit la désintégration du vide et la production de paires particules-antiparticules. Les méthodes perturbatives traditionnelles ne parviennent pas à décrire de manière fiable ce phénomène non perturbatif en temps réel. De plus, les méthodes classiques de Monte Carlo par intégrale de chemin euclidienne rencontrent le « problème du signe » lorsqu'elles traitent l'évolution en temps réel ou des poids complexes, ce qui entraîne une croissance rapide des erreurs statistiques. Bien que le formalisme hamiltonien évite le problème du signe, il nécessite de gérer un espace de Hilbert exponentiellement grand. Les auteurs proposent d'utiliser la simulation quantique numérique pour surmonter ces limitations en représentant directement les états quantiques à plusieurs corps et l'évolution unitaire dans le temps sur des qubits. L'accent est mis spécifiquement sur le modèle de Schwinger en (1+1) dimensions sous un champ électrique externe constant, en étudiant la dynamique hors équilibre à taille finie suite à une trempe de champ électrique.

Méthodologie
L'étude emploie un protocole de simulation quantique numérique basé sur le modèle de Schwinger sur réseau en (1+1) dimensions avec des conditions aux limites ouvertes. La méthodologie se compose de trois étapes principales :

1. Formulation du modèle et cartographie sur qubits :

   • La théorie continue est discrétisée en utilisant le formalisme des fermions échelonnés de Kogut-Susskind.
   • La loi de Gauss est utilisée pour éliminer les degrés de liberté de jauge, exprimant le champ électrique sur chaque lien comme une fonction de la distribution de charge de la matière et du champ aux limites. Cela réduit le modèle à un Hamiltonien purement fermionique avec des interactions à longue portée.
   • Les opérateurs fermioniques sont cartographiés sur des opérateurs de Pauli via la transformation de Jordan-Wigner, résultant en un Hamiltonien de qubits composé de termes d'énergie cinétique, de masse et de champ électrique.

2. Préparation de l'état (VQE) :

   • L'état initial pour la dynamique est l'état du vide à champ nul ($\epsilon=0$).
   • Cet état est préparé à l'aide de l'algorithme VQE (Variational Quantum Eigensolver) avec un ansatz RealAmplitudes efficace pour le matériel et l'optimiseur classique SLSQP.
   • Pour étudier les propriétés statiques du système sous des champs externes variables, les auteurs utilisent également le VQE pour trouver l'état fondamental et la déflation quantique variationnelle (VQD) pour estimer l'énergie du premier état excité.

3. Évolution en temps réel (Trotter-Suzuki) :

   • Le système est soumis à une « trempe » soudaine où un champ électrique externe est activé à $t=0$.
   • L'évolution en temps réel est simulée à l'aide d'une décomposition de Trotter-Suzuki de second ordre pour approximer l'opérateur d'évolution temporelle $e^{-iHt}$.
   • La simulation est réalisée sur des ordinateurs classiques à l'aide du framework Qiskit, sans inclure de bruit quantique, afin de comparer l'algorithme à la diagonalisation exacte (ED).

Contributions clés et résultats
L'article présente des résultats numériques pour des tailles de réseau $N=8$ et $N=12$ (avec des unités de réseau $a=m=g=1$), comparant le protocole de simulation quantique aux références de la diagonalisation exacte (ED).

• Propriétés statiques de l'état fondamental :

  • L'état fondamental du vide préparé par VQE atteint une fidélité élevée ($>99,99\%$) par rapport aux résultats de l'ED.
  • Des balayages sur l'intensité du champ externe $\epsilon$ révèlent un comportement par paliers du condensat chiral et un écart d'énergie s'annulant entre l'état fondamental et le premier état excité à des intensités de champ critiques spécifiques. Cela indique une reconstruction du vide induite par le champ, correspondant à une transition de phase de premier ordre dans des volumes finis.
  • L'intensité du champ critique $\epsilon_c$ présente une dépendance claire à la taille finie en $1/N$. Une extrapolation linéaire donne un champ critique en volume infini de $\epsilon_c(\infty) \approx 0,469$, ce que les auteurs attribuent à l'influence des conditions aux limites ouvertes réduisant l'écart d'énergie.

• Dynamique hors équilibre :

  • Lois de conservation : L'évolution de Trotter de second ordre préserve avec succès la charge $U(1)$ totale et l'énergie totale (dans la précision numérique), validant la stabilité de l'algorithme.
  • Séparation de charge : Lors de l'activation du champ, la distribution spatiale de la charge montre une production et une annihilation périodiques de paires charge-anticharge près des limites ouvertes. L'amplitude et la période de ces oscillations augmentent avec l'intensité du champ électrique.
  • Redistribution de l'énergie : L'énergie totale du champ électrique $H_E(t)$ présente un comportement quasi-périodique, reflétant l'échange d'énergie entre le terme du champ électrique et le secteur fermionique. Ceci est distinct de la périodicité plus simple de la distribution de charge, suggérant des contributions complexes de modes de plus haute énergie.
  • Fidélité du vide : La fidélité de l'état du vide $P_{vac}(t)$ décroît à mesure que le système est éloigné de son état initial. Dans la fenêtre de temps précoce ($0 \le t/a \le 1,0$), ce déclin suit une loi exponentielle effective, $P_{vac}(t) \sim e^{-\gamma_{eff} t}$. Le taux de décroissance effectif $\gamma_{eff}$ augmente de manière monotone avec l'intensité du champ électrique.
  • Effets de taille finie : Les comparaisons entre les systèmes $N=8$ et $N=12$ confirment que les caractéristiques dynamiques observées (séparation de charge aux limites, décroissance de la fidélité) ne sont pas des artefacts numériques dus à de petites tailles de système, mais persistent dans des réseaux plus grands, les systèmes plus vastes montissant une déviation plus rapide par rapport au vide initial.

Signification et affirmations
Les auteurs affirment que leur travail démontre la viabilité de la combinaison de la préparation d'état par VQE avec l'évolution numérique en temps réel (via la décomposition de Trotter-Suzuki) pour l'étude de la dynamique hors équilibre dans les théories de jauge sur réseau à champ fort.

• Validation : Le protocole reproduit avec succès les principales caractéristiques des phénomènes induits par le champ — spécifiquement la séparation de charge aux limites, la décroissance de la fidélité de l'état du vide et la redistribution quasi-périodique de l'énergie — lorsqu'il est comparé à la diagonalisation exacte.
• Utilité méthodologique : Cette étude sert de référence pour la simulation quantique numérique en physique des hautes énergies, montrant que l'évolution numérique explicite de Trotter peut capturer efficacement des dynamiques complexes à champ fort dans des systèmes de réseau finis.
• Perspectives : Les auteurs suggèrent modestement que cette approche constitue une base pour de futures extensions vers des tailles de systèmes plus importantes, des ansatz variationnels améliorés, une implémentation sur du matériel quantique bruyant et l'étude de champs externes dépendants du temps. Ils ne revendiquent pas de réalisation expérimentale immédiate sur le matériel bruyant actuel, mais soulignent que la vérification algorithmique sur des simulateurs classiques est une étape nécessaire.