Efficient classical simulation of two-dimensional long-range systems: Rydberg arrays and beyond

Résumé Technique : Simulation Classique Efficace de Systèmes à Longue Portée en Deux Dimensions

Énoncé du Problème
La méthode de Monte Carlo Variationnel (VMC) est un outil de premier plan pour l'étalonnage des simulateurs quantiques, particulièrement pour les systèmes bidimensionnels (2D) à interactions à longue portée (par exemple, les réseaux d'atomes de Rydberg avec des couplages en $1/r^3$ ou $1/r^6$). Cependant, un goulot d'étranglement computationnel important existe lors de l'évaluation de l'énergie locale ($H_{\text{loc}}$) pour des systèmes possédant des interactions de deux corps « tout-vers-tout » arbitraires.

• Dans le cadre du VMC standard, alors que les coûts d'échantillonnage sont de $O(N)$ pour les États de Réseaux de Tenseurs (TNS) et de $O(N^2)$ pour les États Quantiques Neuronaux (NQS), l'évaluation de $H_{\text{loc}}$ pour les couplages hors-diagonaux à longue portée nécessite la sommation de $O(N^2)$ termes non locaux.
• Comme chaque terme nécessite une évaluation complète du réseau prenant un temps de $O(N)$, le coût total par balayage de Monte Carlo (MC) suit une mise à l'échelle en $O(N^3)$. Cette mise à l'échelle cubique rend l'improductive la simulation classique précise de grands systèmes 2D à longue portée (par exemple, des réseaux $10 \times 10$).
• De plus, une pathologie technique survient lors de l'évolution des TNS à partir d'états produits (fréquents dans les expériences de Rydberg) : aux états produits exacts, le Tenseur de Géométrie Quantique (QGT) développe des directions nulles si la dimension de liaison est supérieure à un, ce qui rend les équations de la Principe Variationnel Temporel (TDVP) singulières.

Méthodologie
Les auteurs introduisent un cadre tVMC à base adaptative par réseaux de tenseurs conçu pour réduire la complexité computationnelle de l'évaluation de l'énergie locale de $O(N^3)$ à $O(N)$, tout en conservant la capacité de gérer des couplages inhomogènes généraux et des géométries réalistes.

1. Échantillonnage de Base Adaptative :

   • Au lieu d'échantillonner dans une base unique fixe (par exemple, la base $z$), le Hamiltonien est décomposé en composantes qui sont diagonales dans des bases spécifiques.
   • Pour le modèle XY dipolaire ($H = H_{XY} + H_Z$), les auteurs divisent l'évolution en $H_X$ (diagonale dans la base $x$) et $H_Y + H_Z$ (diagonale dans la base $y$).
   • Deux chaînes de Markov indépendantes sont maintenues : l'une échantillonnant à partir de $|\langle x|\psi_\theta\rangle|^2$ et l'autre à partir de $|\langle y|U|\psi_\theta\rangle|^2$, où $U$ est l'opérateur de changement de base.
   • Puisque les couplages hors-diagonaux à longue portée s'annulent au sein de chaque base spécifique, l'évaluation de l'énergie locale pour chaque composante devient $O(N)$.
   • Les mises à jour des paramètres sont calculées séparément pour chaque base ($\dot{\theta}_X$ et $\dot{\theta}_Y$) et combinées via $\dot{\theta} = \dot{\theta}_X + U^\dagger \dot{\theta}_Y$.

2. Gestion des Conditions Initiales d'État Produit :

   • Pour résoudre la singularité du QGT lorsqu'on part d'un état produit (où la dimension de liaison $D > 1$), les auteurs emploient une approche de tVMC projeté (p-tVMC) pour le premier pas de temps.
   • Au lieu de résoudre l'équation TDVP singulière, le premier pas ($\delta t$) est déterminé en maximisant la fidélité $|\langle \psi_\theta | V | \psi(0) \rangle|^2$, où $V$ est l'expansion du second ordre de l'opérateur d'évolution temporelle.
   • Cette optimisation est efficace car l'état initial est un état produit, rendant la distribution d'échantillonnage $|\langle s|V|\psi(0)\rangle|^2$ computationnellement traitable. Après ce premier pas, le tVMC standard est repris.

3. Exploitation des Symétries :

   • Le cadre impose explicitement une symétrie de réflexion le long des diagonales du réseau dans l'ansatz TNS. Cela réduit le nombre de paramètres variationnels indépendants d'un facteur quatre, allégeant considérablement les exigences de mémoire GPU et réduisant la variance des mises à jour de paramètres.

Résultats Clés
La méthode a été appliquée pour simuler le protocole de dynamique réelle adiabatique d'un modèle XY dipolaire $10 \times 10$ réalisé dans une expérience récente de Rydberg (Chen et al., Nature 2023), un système dont la taille était auparavant hors de portée de la simulation classique.

• Étalonnage : Sur un système $4 \times 4$, les résultats TNS-tVMC (avec des dimensions de liaison $D=4, 5$) convergent systématiquement vers les résultats de la diagonalisation exacte, reproduisant avec précision la dynamique de relaxation pilotée et les oscillations rapides.
• Simulation $10 \times 10$ :
  • La simulation a suivi l'évolution de l'énergie, montrant que le système se relaxe vers un plateau proche de l'état stationnaire, légèrement au-dessus de l'énergie du fondamental.
  • Corrélations Antiferromagnétiques (AFM) : L'étude a traité une divergence présente dans l'expérience originale, où l'état préparé dynamiquement présentait des corrélations AFM à longue portée supprimées par rapport aux calculs de l'état fondamental par MPS (Matrix Product State).
  • La simulation de la dynamique idéale des auteurs (partant d'un état de Néel parfait avec un désaccord plus grand $\delta_0 = 19.48$) démontre que le rampe adiabatique prescrit produit bien des corrélations AFM spatialement étendues, cohérentes avec l'état fondamental variationnel.
  • Cela suggère que la suppression observée dans l'expérience est probablement due à des imperfections expérimentales (telles qu'une puissance laser limitée empêchant l'utilisation du $\delta_0$ optimal sur l'ensemble $10 \times 10$) plutôt qu'à un échec intrinsèque du protocole adiabatique.
• Convergence : Les résultats ont montré une convergence robuste par rapport à la taille de l'échantillon ($N_s = 81\,920$ vs $163\,840$), confirmant la fiabilité de l'approche par base adaptative.

Signification et Revendications
L'article affirme fournir un outil de benchmarking classique scalable pour les plateformes quantiques à longue portée, spécifiquement les simulateurs de Rydberg et de pièges à ions, qui ont actuellement un besoin urgent de tels outils.

• Réussite Technique : Ce travail résout le goulot d'étranglement $O(N^3)$, permettant des simulations précises de la dynamique en temps réel et des états fondamentaux pour des systèmes 2D à longue portée avec un coût de $O(N)$ par balayage MC.
• Résolution de la Pathologie : Il offre une solution robuste pour faire évoluer les TNS à partir d'états produits, un scénario courant dans les protocoles expérimentaux.
• Impact Scientifique : En fournissant une référence classique directe pour l'expérience de Rydberg $10 \times 10$, ce travail clarifie la nature de la divergence entre les résultats expérimentaux et les prédictions théoriques de l'état fondamental. Il démontre que le protocole adiabatique idéal est capable de générer l'ordre à longue portée attendu, déplaçant ainsi l'attention de l'optimisation expérimentale vers les imperfections de réalisation.
• Perspectives : Les auteurs positionnent cette méthode comme un pont pour la collaboration classique-quantique, où les simulations classiques peuvent guider l'identification de configurations physiques intéressantes avant que la croissance de l'intrication ne devienne trop sévère pour les réseaux de tenseurs, moment auquel la dynamique peut être confiée aux simulateurs quantiques. La méthode est également notée comme une nouvelle voie pour l'étude des systèmes 2D à longue portée présentant des problèmes de signe, là où d'autres méthodes numériques peinent.