Extracting conserved operators from a projected entangled pair state

Cet article présente une méthode basée sur la géométrie quantique et les facteurs de structure statiques pour extraire des opérateurs conservés géométriquement locaux, tels que des Hamiltoniens parents améliorés, dont un état iPEPS donné est un état propre, permettant notamment d'identifier des Hamiltoniens pour les états RVB à courte portée et les codes toriques déformés.

L'essentiel

Imagine que vous avez une photo très complexe d'une foule de personnes (c'est l'état quantique, ou « PEPS »). Vous savez que cette foule bouge selon certaines règles invisibles, comme si elle dansait sur une musique précise. Mais vous ne connaissez pas la partition de cette musique (l'Hamiltonien, ou la loi physique qui régit le système).

La question est : Comment retrouver la partition de musique en regardant seulement la photo de la foule ?

C'est exactement ce que font les auteurs de cet article. Ils ont développé une méthode pour deviner les règles cachées (les opérateurs conservés, y compris l'énergie) d'un système quantique complexe, simplement en analysant les « motifs » de ses corrélations.

Voici une explication simple, avec des analogies :

1. Le problème : Le puzzle inversé

En physique, on a souvent l'inverse : on connaît la musique (la loi physique) et on essaie de deviner comment la foule va danser (l'état du système). Ici, on veut faire l'inverse. On a l'état (la danse), et on veut retrouver la musique.

Le problème, c'est que les systèmes quantiques en 2D (comme une grille de pixels géante) sont extrêmement complexes. Les calculs exacts sont souvent impossibles, un peu comme essayer de compter chaque grain de sable d'une plage à la main. On utilise donc des approximations (des « esquisses » rapides) qui contiennent de petites erreurs. La grande question était : Peut-on encore retrouver la musique si notre photo est un peu floue ?

2. La solution : La « carte de résonance » (Facteurs de structure)

Les auteurs utilisent une astuce géniale. Ils ne regardent pas juste la photo, ils regardent comment les gens réagissent les uns aux autres.

Imaginez que vous tapez dans vos mains à un rythme précis (une onde). Si la foule résonne parfaitement avec ce rythme, cela signifie qu'elle suit une règle spécifique.

• L'outil : Ils utilisent ce qu'on appelle les « facteurs de structure statiques ». C'est comme une carte qui montre à quel point la foule « vibre » ou « résonne » pour différents rythmes.
• L'analogie : Si vous tapez un rythme qui ne correspond à aucune règle de la danse, la foule bouge de façon chaotique (le signal est fort). Mais si vous tapez le rythme exact de la musique cachée, la foule reste parfaitement synchronisée et ne crée pas de « bruit » (le signal s'annule).

3. La méthode : Le test du « zéro »

Leur méthode consiste à tester des millions de rythmes possibles (des opérateurs) et à chercher ceux qui font disparaître le bruit.

• Mathématiquement, ils cherchent les solutions où le « facteur de structure » est zéro.
• C'est comme chercher la clé parfaite qui ouvre une serrure : quand vous trouvez la bonne clé (le bon opérateur), la porte s'ouvre sans effort (l'énergie de fluctuation est nulle).

Même si leur photo (l'état PEPS) est une approximation imparfaite, leur méthode est si robuste qu'elle arrive à retrouver la musique exacte, ou du moins une très bonne copie, malgré le bruit de fond.

4. Les découvertes surprenantes

En appliquant cette méthode, ils ont trouvé des choses fascinantes :

• La musique cachée du RVB : Ils ont pris un état quantique mystérieux appelé « RVB » (Resonating Valence Bond), qui est une sorte de superposition de liaisons entre atomes. Ils ont réussi à trouver la musique (l'Hamiltonien) qui fait danser cette foule. C'est une musique locale (les voisins parlent entre eux) et elle est plus simple que ce que l'on pensait auparavant.
• Les « cicatrices quantiques » (Quantum Many-Body Scars) : C'est la découverte la plus cool. Ils ont trouvé une musique où l'état de la foule n'est pas le début du concert (l'état fondamental), mais une note spécifique jouée au milieu de la symphonie.
  • L'analogie : Imaginez un orchestre qui joue une symphonie chaotique. Normalement, tout le monde se mélange. Mais ici, il y a un groupe de musiciens qui, même au milieu du chaos, continue de jouer une mélodie parfaite et répétitive. Ces musiciens sont les « cicatrices ». Ils résistent au chaos thermique. C'est crucial pour stocker de l'information quantique sans qu'elle ne s'effondre.

En résumé

Les auteurs ont créé un détective mathématique.

1. Il prend une « photo » approximative d'un système quantique complexe.
2. Il teste des milliers de règles possibles en regardant comment le système « résonne ».
3. Il identifie les règles qui font que le système reste calme (zéro bruit).
4. Il en déduit la loi physique (l'Hamiltonien) qui régit ce système, même si la photo de départ n'était pas parfaite.

C'est une avancée majeure car cela permet de comprendre la physique profonde de matériaux complexes (comme les supraconducteurs) simplement en observant leurs états, sans avoir besoin de connaître la théorie complète au préalable. C'est comme pouvoir deviner la recette d'un gâteau en goûtant une miette, même si la miette est un peu brûlée.

Résumé technique

Titre : Extraction d'opérateurs conservés à partir d'un état de paires intriquées projetées (iPEPS)

1. Problématique

L'article aborde un problème inverse fondamental en physique de la matière condensée et en simulation quantique : comment reconstruire l'hamiltonien (ou d'autres opérateurs conservés) sous-jacent à partir d'un état quantique donné, sans connaître la dynamique qui l'a produit ?

• Contexte : Avec l'avènement des simulateurs quantiques et des dispositifs NISQ, il est crucial de certifier que les plateformes implémentent bien les modèles souhaités. De plus, pour la matière fortement corrélée, il est difficile de dériver des hamiltoniens effectifs à basse énergie ($H_{eff}$) à partir de modèles microscopiques complexes.
• Défi spécifique : Les états de réseaux de tenseurs, en particulier les PEPS (Projected Entangled Pair States) en 2D, sont des ansatz puissants pour décrire les états fondamentaux. Cependant, contrairement aux MPS (Matrix Product States) en 1D où des algorithmes efficaces existent, l'extension à la 2D est bloquée par la complexité du calcul exact (problème #P-complet) et la nécessité d'utiliser des schémas de contraction approximatifs (comme CTMRG). La question centrale est : peut-on extraire des informations fiables sur l'hamiltonien parent malgré les erreurs d'approximation inhérentes aux calculs de PEPS ?

2. Méthodologie

Les auteurs proposent un algorithme classique pour extraire des opérateurs conservés géométriquement $k$-locaux (dont l'hamiltonien) pour lesquels un état PEPS donné $|\Psi\rangle$ est un état propre (exact ou approché).

A. Fondement théorique : Les facteurs de structure statique
La méthode repose sur l'observation que les opérateurs conservés $\hat{H}$ correspondent au noyau (kernel) de la matrice des facteurs de structure statique $S$.
Pour une base d'opérateurs hermitiens locaux $\{\hat{o}^\alpha_x\}$, le facteur de structure est défini comme la transformée de Fourier des fonctions de corrélation connectées :
$$S^{\alpha,\beta}(q) = \sum_x e^{-iq\cdot x} \left( \frac{\langle \Psi | \hat{o}^\alpha_x \hat{o}^\beta_0 | \Psi \rangle}{\langle \Psi | \Psi \rangle} - \frac{\langle \Psi | \hat{o}^\alpha_x | \Psi \rangle \langle \Psi | \hat{o}^\beta_0 | \Psi \rangle}{\langle \Psi | \Psi \rangle^2} \right)$$
Si $|\Psi\rangle$ est un état propre de $\hat{H} = \sum_x e^{iq\cdot x} \hat{h}_x$, alors la fluctuation quantique de $\hat{H}$ dans cet état est nulle. Cela implique que le vecteur de coefficients définissant $\hat{H}$ appartient au noyau de la matrice $S$ (les valeurs propres de $S$ associées sont nulles).

B. Évaluation via la fonction génératrice
Pour calculer ces facteurs de structure sur des iPEPS infinis, les auteurs utilisent la méthode de la fonction génératrice :

1. Ils définissent un état déformé $|\langle G_\alpha(\mu, q)| = \langle \Psi | \prod_x (1 + \mu e^{-iq\cdot x} \hat{o}^\alpha_x)$.
2. Le facteur de structure est obtenu par dérivation par rapport au paramètre scalaire $\mu$ à $\mu=0$ :
   $$S^{\alpha,\beta}(q) = \left. \frac{d}{d\mu} \frac{\langle G_\alpha(\mu, q) | \hat{o}^\beta_0 | \Psi \rangle}{\langle G_\alpha(\mu, q) | \Psi \rangle} \right|_{\mu=0}$$
3. Implémentation numérique :
   • L'opérateur produit est traité comme un iPEPO (Infinite Projected Entangled Pair Operator).
   • La contraction est effectuée via la méthode CTMRG (Corner Transfer Matrix Renormalization Group).
   • Pour éviter l'instabilité de la différenciation automatique sur les opérateurs multi-sites, les auteurs utilisent une différence finie centrale à 5 points pour calculer la dérivée.
   • La complexité est polynomiale en la dimension de liaison du PEPS ($D$).

C. Géométrie quantique
L'article établit un lien théorique : la matrice $S$ est proportionnelle au tenseur métrique quantique de la variété d'états déformés. Les opérateurs conservés correspondent à une susceptibilité de fidélité nulle (vanishing fidelity susceptibility).

3. Résultats Clés et Benchmarks

Les auteurs valident leur méthode sur plusieurs cas d'étude, démontrant sa robustesse face aux erreurs d'approximation.

• État AKLT 2D (État exact) :

  • Sur un réseau carré, la méthode retrouve exactement les 4 valeurs propres nulles de la matrice $S$ (identité et les 3 opérateurs de spin $S^x, S^y, S^z$) pour le cas 1-site.
  • Pour le cas 2-sites, elle identifie les 81 hamiltoniens parents sans frustration (frustration-free) connus, confirmant la précision de la méthode sur des états exacts.

• Modèle XX 2D (État variationnel) :

  • Application à l'état fondamental variationnel d'un modèle XX avec champ magnétique en $Z$ alterné.
  • La méthode réussit à extraire le hamiltonien non-frustré (non-frustration-free) et les générateurs de symétrie (U(1)), même lorsque la symétrie est brisée spontanément (phase ferromagnétique).
  • Elle retrouve le hamiltonien original avec une précision de 4 chiffres significatifs, prouvant qu'elle fonctionne au-delà des constructions standards de hamiltoniens parents sans frustration.

• État RVB à courte portée (État critique) :

  • Application à l'état RVB (Resonating Valence Bond) sur réseau carré, qui possède une longueur de corrélation divergente.
  • Malgré les difficultés de convergence du CTMRG pour les systèmes critiques, la méthode extrait un hamiltonien parent local à 4 sites (plaquette).
  • Ce hamiltonien est une superposition de termes d'échange de Heisenberg (1er et 2ème voisins) et d'interactions à 4 spins.
  • L'énergie de l'état RVB sur ce hamiltonien est très proche de l'énergie fondamentale, confirmant que l'état RVB est un bon état propre approché. Ce résultat est significatif car la construction standard de hamiltoniens parents pour RVB nécessite des termes à 8 sites ou plus.

• Ondes de Ising déformées et Toric Code (États excités) :

  • La méthode découvre un hamiltonien pour lequel les états de Ising déformés (ou les états Toric Code déformés) ne sont pas l'état fondamental, mais des états propres excités situés au milieu du spectre.
  • Ces états, qui satisfont la loi d'aire pour l'entropie d'intrication, sont identifiés comme des cicatrices quantiques à plusieurs corps (Quantum Many-Body Scars).
  • Le hamiltonien trouvé commute avec le terme de déformation, ce qui garantit que la famille d'états déformés reste un état propre pour tout paramètre de déformation.

4. Contributions Majeures

1. Algorithme général : Développement d'une méthode systématique pour extraire des opérateurs conservés à partir de n'importe quel état PEPS, qu'il soit exact, variationnel, gappé ou critique.
2. Robustesse aux erreurs : Démonstration que l'extraction d'hamiltoniens parents est possible même avec les erreurs inévitables des schémas de contraction approchés (CTMRG) et des états variationnels non parfaits.
3. Localité améliorée : La méthode produit des hamiltoniens parents avec une localité bien meilleure (ex: 4 sites pour RVB) que les constructions standards (souvent 8 sites ou plus), facilitant potentiellement leur réalisation expérimentale.
4. Découverte de nouvelles phases : Identification d'hamiltoniens possédant des états à faible intrication comme états excités (cicatrices quantiques), ouvrant la voie à l'étude de la thermalisation et de la localisation à plusieurs corps.
5. Lien géométrique : Établissement d'une connexion formelle entre les facteurs de structure statique, la métrique quantique et la susceptibilité de fidélité.

5. Signification et Perspectives

Ce travail démontre que l'apprentissage d'hamiltoniens (Hamiltonian Learning) à partir d'états de réseaux de tenseurs est réalisable en temps polynomial, contredisant l'idée que la reconstruction est impossible en raison du manque d'information dans l'état fondamental.

• Applications : La méthode offre un outil puissant pour valider les simulateurs quantiques, découvrir de nouveaux modèles effectifs pour la matière condensée, et explorer des phénomènes exotiques comme les cicatrices quantiques.
• Limites et futurs : Les auteurs soulignent la nécessité d'améliorer la précision pour les iPEPS critiques (via des méthodes Monte Carlo ou des analyses d'échelle) et d'explorer l'extraction d'opérateurs conservés à partir d'états de Gibbs (température finie) représentés par des PEPO.

En résumé, cet article fournit un cadre théorique et numérique robuste pour "lire" la physique sous-jacente directement depuis la structure d'intrication d'un état quantique simulé, dépassant les limitations des méthodes traditionnelles de construction de hamiltoniens parents.