Eigenstate Thermalization Hypothesis (ETH) for off-diagonal matrix elements in integrable spin chains

En utilisant des résultats récents de l'Ansatz de Bethe algébrique, cette étude démontre que dans la chaîne de Heisenberg isotrope, les éléments de matrice hors-diagonale d'opérateurs locaux décroissent exponentiellement avec la taille du système pour des états appartenant au même macro-état thermique (suivant une loi de Gumbel), tandis que leur décroissance est beaucoup plus rapide pour des états de macro-états différents.

L'essentiel

🎻 L'Orchestre Quantique : Quand la musique s'arrête-t-elle ?

Imaginez un gigantesque orchestre composé de millions d'instruments (des atomes) qui jouent tous en même temps. En physique quantique, cet orchestre est un système isolé : personne ne l'écoute, personne ne le dirige, il joue tout seul selon les règles strictes de la mécanique quantique.

Le grand mystère que ces chercheurs (Federico Rottoli et Vincenzo Alba) tentent de résoudre est le suivant : Comment un tel orchestre, qui ne fait que jouer, finit-il par sembler "calme" et "chaud" (comme un objet qui a atteint l'équilibre thermique) ?

Pour comprendre cela, ils étudient comment les différents instruments (les états d'énergie) "parlent" entre eux. C'est là qu'intervient l'hypothèse de l'ETH (Eigenstate Thermalization Hypothesis), qui est comme une règle d'or pour les systèmes chaotiques (imprévisibles).

1. La Règle du Chaos vs. La Règle de l'Intégrabilité

• Les systèmes chaotiques (Le Jazz) : Dans un système désordonné, les instruments sont si bien connectés que si vous écoutez deux notes très proches, elles se mélangent parfaitement. La règle dit que les connexions entre les notes sont faibles et suivent une distribution normale (comme une cloche de Gauss). C'est la "santé" de la thermalisation.
• Les systèmes intégrables (La Musique Classique Stricte) : Ici, comme dans la chaîne de spins étudiée (le modèle XXX), il y a trop de règles conservées. C'est comme un orchestre où chaque musicien a une partition rigide et ne peut pas improviser. On pensait que cela empêchait la thermalisation.

Les chercheurs se sont demandé : "Même dans cette musique très stricte, comment les notes se parlent-elles ?"

2. L'Expérience : Mesurer les "Chuchotements"

Pour répondre, ils ont utilisé une technique de pointe appelée Ansatz de Bethe (une sorte de super-calculatrice mathématique) pour calculer comment deux états d'énergie différents interagissent via un petit opérateur (une petite touche sur un instrument).

Ils ont regardé deux situations :

A. Les musiciens du même groupe (Même macro-état)
Imaginez deux musiciens qui jouent dans le même style, à la même température.

• Ce qu'ils ont trouvé : Les "chuchotements" (les éléments de matrice hors-diagonale) entre eux sont très faibles. Ils décroissent exponentiellement quand l'orchestre grandit. C'est comme si le bruit de fond devenait si fort que deux musiciens ne peuvent plus se parler distinctement.
• La surprise : La distribution de ces chuchotements ne ressemble pas à une cloche de Gauss (comme on le pensait pour le chaos). Elle ressemble plutôt à une distribution de Gumbel.
  • L'analogie : Imaginez que vous lancez des dés. Dans un système chaotique, les résultats forment une courbe en cloche. Dans ce système intégrable, c'est comme si vous regardiez les records de vitesse du vent : il y a beaucoup de vents faibles, mais quelques rafales extrêmes qui tirent la moyenne. C'est une statistique différente, plus "sauvage".

B. Les musiciens de groupes opposés (Différents macro-états)
Imaginez maintenant un musicien jouant une symphonie lente (température zéro) et un autre jouant du rock rapide (température infinie).

• Ce qu'ils ont trouvé : Leurs "chuchotements" sont encore plus faibles ! Ils disparaissent beaucoup plus vite (comme $e^{-L^2}$). C'est comme si les deux groupes étaient dans des univers parallèles qui ne se touchent presque plus.

3. Pourquoi est-ce important ?

C'est comme si on découvrait que même dans un système très ordonné et prévisible (la chaîne de spins), la nature trouve un moyen de "thermaliser" (de se mélanger) d'une manière très subtile.

• Le résultat clé : Même si le système est "intégrable" (trop de règles), les éléments de matrice hors-diagonale suivent une loi de décroissance similaire aux systèmes chaotiques, mais avec une statistique différente (Gumbel au lieu de Gaussien).
• L'implication : Cela signifie que la frontière entre "système ordonné" et "système chaotique" n'est pas aussi nette qu'on le pensait. La thermalisation peut survenir même dans des systèmes très structurés, mais elle laisse une "signature" statistique particulière dans la façon dont les états interagissent.

En résumé

Cette étude est comme un examen de routine pour un orchestre quantique. Les chercheurs ont écouté comment les notes se parlent entre elles. Ils ont découvert que :

1. Même dans un système très ordonné, les interactions entre états voisins s'effacent très vite (ce qui permet une forme de thermalisation).
2. La façon dont ces interactions sont réparties (les "records" de force des chuchotements) est unique et différente de celle des systèmes chaotiques.

C'est une preuve que la nature est subtile : même avec des règles strictes, le chaos et l'équilibre thermique peuvent émerger, mais avec une "signature" mathématique bien spécifique.

Résumé technique

1. Problématique et Contexte

L'étude de la thermalisation dans les systèmes quantiques isolés hors équilibre repose fondamentalement sur l'Hypothèse de Thermalisation des États Propres (ETH). Pour les systèmes chaotiques (non intégrables), l'ETH postule que les éléments de matrice d'un opérateur local $O$ entre deux états propres d'énergie $|E_i\rangle$ et $|E_j\rangle$ suivent une forme spécifique :
$$\langle E_i| O |E_j\rangle = O(\bar{E}) \delta_{ij} + f_O(\omega, E) e^{-S(E)/2} R_{ij}$$
où $R_{ij}$ est une variable aléatoire gaussienne de moyenne nulle et de variance unité. Cela implique que les éléments diagonaux convergent vers une fonction lisse de l'énergie, tandis que les éléments hors-diagonaux sont exponentiellement supprimés avec la taille du système $L$.

Cependant, dans les systèmes intégrables (comme la chaîne de Heisenberg XXX), la présence d'un nombre extensif de charges conservées viole l'ETH standard. Une question ouverte majeure est de comprendre comment les éléments de matrice hors-diagonaux se comportent dans ces systèmes, en particulier la dépendance en taille finie et la statistique de ces éléments. Des travaux récents sur le modèle de Lieb-Liniger (théorie des champs 1D) ont suggéré des comportements spécifiques, mais l'extension à des modèles de réseaux intégrables (chaînes de spins) avec des états liés (strings) restait à explorer.

2. Méthodologie

Les auteurs se concentrent sur la chaîne de Heisenberg isotrope spin-1/2 (modèle XXX). Pour surmonter les limitations des méthodes de diagonalisation exacte (limitées à de petites tailles de système), ils utilisent des résultats avancés de l'Ansatz de Bethe Algébrique (ABA).

• Opérateurs étudiés : Ils analysent les éléments de matrice d'opérateurs locaux agissant sur un site (opérateurs à un spin, ex: $E^{(11)}_i$) et sur deux sites voisins (opérateurs à deux spins, ex: $E^{(11)}_i E^{(22)}_{i+1}$).
• Échantillonnage des états : Ils considèrent des états propres extraits de macro-états thermodynamiques, notamment l'ensemble de Gibbs à température infinie et à température finie, ainsi que l'état fondamental (température nulle).
• Traitement des singularités : Le calcul des éléments de matrice via l'ABA pour des états contenant des "strings" (états liés de quasi-particules) est complexe en raison de singularités fictives.
  • Pour les opérateurs à un spin, les auteurs utilisent une procédure analytique pour supprimer les singularités liées aux déviations de string, permettant d'étudier des chaînes jusqu'à $L \approx 500$.
  • Pour les opérateurs à deux spins, la complexité est plus élevée. Ils se restreignent donc aux états sans strings (sans états liés) pour atteindre $L \approx 60$.
• Quantité mesurée : Ils analysent la distribution de probabilité de $M_{ij} = \ln |\langle E_i| O |E_j\rangle|^2$.

3. Contributions Clés et Résultats

A. Éléments de matrice entre états du même macro-état

Lorsque les états $|E_i\rangle$ et $|E_j\rangle$ appartiennent au même macro-état thermique (même densité d'énergie et de magnétisation) :

1. Décroissance exponentielle : Les éléments de matrice hors-diagonaux décroissent exponentiellement avec la taille du système $L$. La forme est compatible avec :
   $$|\langle E_i| O |E_j\rangle| \propto \exp\left( -c_O L \ln L - f_{M_{ij}} L \right)$$
   Cela confirme que la présence de strings (états liés) dans le modèle XXX ne modifie pas qualitativement le scénario observé dans le modèle de Lieb-Liniger.
2. Statistique (Distribution de Gumbel) : Contrairement aux systèmes chaotiques où les fluctuations sont gaussiennes, la distribution de probabilité de $M_{ij}$ (ou de sa partie résiduelle après soustraction de la tendance moyenne) est parfaitement décrite par une distribution de Gumbel.
   • Les paramètres de cette distribution sont indépendants de $L$.
   • Cela implique que les éléments de matrice $\langle E_i| O |E_j\rangle$ ne suivent pas une loi gaussienne, ce qui est une signature forte de l'intégrabilité.
3. Ratio $\Gamma$ : Les auteurs calculent le ratio $\Gamma = \overline{|O_{ij}|^2} / (\overline{|O_{ij}|})^2$. Pour une distribution gaussienne, $\Gamma = \pi/2$. Ils trouvent que $\Gamma$ diverge exponentiellement avec $L$ (atteignant $\sim 10^5$), confirmant la non-gaussianité et la nature "lourde" de la queue de la distribution.

B. Éléments de matrice entre états de macro-états différents

Lorsque les états appartiennent à des macro-états thermodynamiques différents (par exemple, état fondamental vs état à température infinie) :

1. Décroissance plus rapide : Les éléments de matrice décroissent beaucoup plus vite, selon une loi :
   $$|\langle E_i| O |E_j\rangle| \propto \exp(-d_O L^2)$$
2. Statistique : La distribution des logarithmes de ces éléments suit également une loi de Gumbel, mais avec une échelle de décroissance quadratique en $L$.

C. Opérateurs à deux spins

Malgré la limitation de taille de système ($L \le 60$) due à la complexité computationnelle, les résultats pour les opérateurs à deux spins confirment le même scénario que pour les opérateurs à un spin : décroissance exponentielle en $L$ (avec correction logarithmique) et statistique de Gumbel pour les états du même macro-état.

4. Signification et Implications

• Validation du scénario d'intégrabilité : L'étude confirme que le comportement des éléments de matrice hors-diagonaux dans les modèles de spins intégrables (XXX) est universel par rapport aux modèles de champs intégrables (Lieb-Liniger), malgré la présence de structures complexes comme les strings.
• Distinction Intégrable vs Chaotique : La découverte que la statistique des éléments de matrice suit une loi de Gumbel (et non gaussienne) fournit un critère robuste pour distinguer les systèmes intégrables des systèmes chaotiques, au-delà de la simple vérification de la décroissance exponentielle.
• Théorie des Probabilités Libres : Les résultats ouvrent la voie à l'application de la Théorie des Probabilités Libres (Free Probability Theory) aux modèles intégrables, suggérant que les structures algébriques sous-jacentes (comme l'ABA) dictent des statistiques de valeurs extrêmes (Gumbel) plutôt que centrales (Gaussiennes).
• Dynamique hors équilibre : Comprendre la statistique de ces éléments de matrice est crucial pour reconstruire la dynamique temporelle complète des observables après un quench quantique, reliant ainsi les propriétés spectrales statiques à la relaxation dynamique.

En résumé, cet article établit une caractérisation précise et quantitative de l'ETH dans les chaînes de spins intégrables, démontrant que l'intégrabilité se manifeste non seulement par la violation de la thermalisation standard, mais aussi par une statistique spécifique (Gumbel) des fluctuations des éléments de matrice hors-diagonaux.