Numerical evidence for the non-Abelian eigenstate thermalization hypothesis

Cet article fournit des preuves numériques soutenant l'hypothèse de la thermalisation des états propres (ETH) non abélienne à travers des simulations d'une chaîne de Heisenberg unidimensionnelle et propose une preuve analytique de son auto-cohérence, établissant ainsi un cadre pour la compréhension de la thermalisation dans les systèmes quantiques possédant des quantités conservées non commutatives.

L'essentiel

Imaginez que vous possédez une machine géante et complexe composée de 18 petits interrupteurs (qubits) qui sont tous connectés entre eux. Dans le monde de la physique quantique, ces interrupteurs ne se contentent pas de s'allumer ou de s'éteindre ; ils tournent dans différentes directions, et la direction d'un interrupteur affecte ses voisins.

Pendant longtemps, les physiciens ont cru que si l'on laissait une telle machine fonctionner pendant longtemps, elle finirait par se « stabiliser » dans un état moyen prévisible, un peu comme une tasse de café chaud qui refroidit jusqu'à atteindre la température ambiante. Cette idée est appelée l'Hypothèse de Thermalisation des États Propres (ETH - Eigenstate Thermalization Hypothesis). Elle suggère que, peu importe la manière dont vous démarrez la machine, les parties locales finiront par se comporter comme si elles étaient dans un équilibre « thermique » (chaud/froid) standard.

Le Problème : L'énigme du « Non-Commutatif »
Cependant, il y a un piège : dans cette machine spécifique, les interrupteurs sont régis par une règle spéciale appelée symétrie non-abélienne. Considérez cela comme un jeu de directions :

• Si vous tournez une boussole vers le Nord, puis vers l'Est, vous finissez à un endroit différent de si vous aviez tourné vers l'Est, puis vers le Nord.
• En termes quantiques, ces « directions » (charges) ne s'entendent pas ; elles ne commutent pas. Elles interfèrent les unes avec les autres.

À cause de cette interférence, les anciennes règles (l'ETH standard) s'effondrent. La machine ne devrait pas se stabiliser de la manière habituelle. Mais une nouvelle théorie, appelée l'ETH Non-Abélienne, a été proposée pour expliquer comment ce type de machine finit quand même par se stabiliser, mais avec un ensemble de règles différent.

Ce que cet article a fait
Les auteurs de cet article ont agi comme des détectives testant une nouvelle théorie. Ils ont construit une simulation informatique de cette machine à 18 interrupteurs pour voir si la nouvelle théorie de l'« ETH Non-Abélienne » était vraie.

Voici ce qu'ils ont trouvé, en utilisant des analogies simples :

1. Le motif « lissé » : Ils ont observé le calcul interne de la machine. La théorie prédisait que si vous traciez le comportement de la machine en fonction de son énergie, les points devraient former une courbe fluide et continue (comme une colline douce) plutôt qu'un amas irrégulier et aléatoire. Résultat : Les données ont formé de magnifiques bandes lisses, exactement comme la théorie le prédisait.
2. Le « bruit aléatoire » : La théorie affirmait également que si l'on regardait de près les infimes différences entre les points, elles devraient ressembler à de la neige statique sur un vieux téléviseur (distribution gaussienne). Résultat : En zoomant, le « bruit » ressemblait exactement à de la neige statique aléatoire.
3. La vérification du « volume » : La théorie prédisait une relation spécifique entre la « puissance » du bruit et le nombre d'états différents dans lesquels la machine peut se trouver (la densité d'états). C'est comme dire : « Si la pièce devient plus grande, l'écho doit devenir plus faible d'une manière très spécifique. » Résultat : L'écho est devenu plus faible exactement au rythme prévu par la théorie.
4. Le test du « ratio » : Ils ont comparé le bruit à l'intérieur des groupes principaux de la machine par rapport au bruit entre les groupes. La théorie disait que ce ratio devrait être exactement de 2. Résultat : Leurs mesures sont arrivées à 1,99, ce qui est pratiquement 2.

La preuve de l'« auto-cohérence »
Au-delà de la simulation informatique, les auteurs ont également réalisé une preuve mathématique. Ils ont démontré que la nouvelle théorie ne se contredit pas elle-même. Ils ont dû ajuster légèrement une définition de l'« entropie » (une mesure du désordre) — en soustrayant une petite quantité liée au spin de la machine — pour que les mathématiques fonctionnent parfaitement. Une fois ce petit ajustement effectué, la théorie tenait bon sans aucune faille logique.

L'essentiel
Cet article fournit la première preuve numérique solide que l'ETH Non-Abélienne est réelle. Il confirme que même lorsque les particules quantiques ont des règles « conflictuelles » (des charges non-commutatives) qui les empêchent de se comporter normalement, elles trouvent quand même un moyen de se thermaliser, mais elles suivent de nouvelles instructions, légèrement plus complexes que ce que nous pensions auparavant.

Les auteurs n'ont pas prétendu que cela mènerait à de nouveaux traitements médicaux ou à une technologie immédiate. Au lieu de cela, ils ont réussi à prouver que ce cadre théorique spécifique sur la façon dont les systèmes quantiques se stabilisent est mathématiquement sain et correspond à leurs modèles informatiques.

Résumé technique

Résumé Technique : Preuves Numériques de l'Hypothèse de Thermalisation des États Propres Non-Abélienne

Énoncé du Problème
L'Hypothèse de Thermalisation des États Propres (ETH) explique comment les systèmes quantiques à plusieurs corps isolés et génériques se thermalisent de manière interne, postulant que les valeurs d'attente temporelles moyennes des opérateurs locaux sont égales à leurs valeurs d'attente thermiques. Cependant, des travaux récents par Murthy et al. ont démontré que l'ETH standard échoue en présence de symétries non-abéliennes. Dans de tels systèmes, les quantités conservées (charges) ne commutent pas entre elles, ce qui conduit à des dégénérescences et à l'absence de bases propres communes pour toutes les charges. Cela entre en conflit avec la structure standard de l'ETH et le théorème de Wigner–Eckart, qui dicte la structure des fréquences de transition et des éléments de matrice dans les systèmes dotés de symétries continues. Bien que Murthy et al. aient proposé une « NA-ETH » (ETH non-abélienne) pour résoudre ces conflits, et que Noh ait fourni un soutien numérique préliminaire dans un modèle 2D, une vérification numérique complète dans un système présentant un minimum de symétries extrinsèques faisait défaut.

Méthodologie
Les auteurs modélisent une chaîne unidimensionnelle (1D) de $N=18$ qubits soumise à des conditions aux limites ouvertes. Le système évolue sous un Hamiltonien de Heisenberg non intégrable présentant des interactions de plus proche et de plus proche voisin. L'Hamiltonien est conçu pour présenter une symétrie $SU(2)$ (conservant les composantes du spin total $S^x, S^y, S^z$) tout en brisant explicitement les symétries de translation et d'inversion spatiale pour garantir la non-intégrabilité.

L'étude utilise la diagonalisation exacte pour calculer la base propre jointe $\{|\alpha, m\rangle\}$ partagée par l'Hamiltonien $H$, l'opérateur du spin total au carré $\vec{S}^2$, et la composante $z$ du spin $S^z$. Ici, $\alpha$ étiquette l'énergie $E_\alpha$ et le nombre quantique du spin total $s_\alpha$, tandis que $m$ est le nombre quantique magnétique. Les opérateurs locaux sont représentés par des opérateurs tensoriels sphériques $T^{(k)}_q$. Les auteurs analysent les éléments de matrice réduits $\langle \alpha || T^{(k)} || \alpha' \rangle$, qui sont indépendants des nombres quantiques magnétiques en vertu du théorème de Wigner–Eckart.

L'analyse se concentre sur la vérification de l'ansatz proposé par Murthy et al., qui postule que les éléments de matrice réduits prennent la forme :
$$\langle \alpha || T^{(k)} || \alpha' \rangle = T^{(k)}(E, S) \delta_{\alpha\alpha'} + e^{-S_{th}(E, S)/2} f^{(T)}_\nu(E, S, \omega) R^{(T)}_{\alpha\alpha'}$$
où $E$ et $S$ sont l'énergie et le spin moyens, $\omega$ et $\nu$ sont leurs différences, $R^{(T)}$ représente les fluctuations erratiques, et $S_{th}$ est une entropie thermodynamique modifiée définie comme l'entropie standard moins $\ln(2S+1)$ pour tenir compte de la dégénérescence magnétique.

Contributions Clés et Résultats
L'article fournit des preuves numériques étendues soutenant la NA-ETH à travers sept prédictions spécifiques, parallèlement à une preuve analytique d'auto-cohérence :

1. Auto-cohérence Analytique : Les auteurs prouvent que la NA-ETH est auto-cohérente (au sens de Srednicki) uniquement si l'entropie thermodynamique $S_{th}$ est définie pour exclure le facteur de dégénérescence magnétique $\ln(2S+1)$. Cette définition garantit que l'échelle des éléments de matrice dérivée directement de l'ansatz correspond à l'échelle dérivée indirectement via la composition d'opérateurs (en utilisant les coefficients de Clebsch–Gordan).

2. Vérification de la Non-Intégrabilité : Les statistiques des écarts énergétiques au sein des sous-espaces propres conjoints de $\vec{S}^2$ et $S^z$ suivent la distribution de l'Ensemble Orthogonal Gaussien (GOE), confirmant que le système est non intégrable et approprié pour les tests de l'ETH.

3. Lissé de la Diagonale de Bloc : Les éléments de matrice réduits diagonaux de bloc $\langle \alpha || T^{(1)} || \alpha \rangle$ présentent une dépendance lisse de la densité d'énergie ($E/N$) et de la densité de spin ($s_\alpha/N$) dans la limite thermodynamique, ce qui est cohérent avec la fonction lisse $T^{(k)}(E, S)$ de l'ansatz.

4. Fluctuations Gaussiennes : Après soustraction de la moyenne lisse, les fluctuations des éléments de matrice réduits, tant diagonaux de bloc que hors-diagonaux de bloc, suivent des distributions gaussiennes, soutenant le terme de matrice aléatoire $R^{(T)}_{\alpha\alpha'}$.

5. Structure Hors-Bloc : Les éléments hors-diagonaux de bloc présentent également des fluctuations gaussiennes et dépendent d'une fonction lisse $f^{(T)}_\nu(E, S, \omega)$ qui décroît rapidement lorsque la différence d'énergie $|\omega|$ augmente.

6. Ratio de Variance : Le ratio de la variance des éléments diagonaux de bloc par rapport aux éléments hors-diagonaux de bloc est en moyenne de $2$, une prédiction dérivée de la théorie des matrices aléatoires.

7. Échelle de Variance : Les variances des éléments diagonaux de bloc et hors-diagonaux de bloc décroissent algébriquement avec la densité d'états (DOS). Les pentes observées dans les graphiques log-log ($\approx -0,83$ et $\approx -0,95$) sont proches de la pente prédite de $-1$, les écarts étant attribués à des effets de taille finie.

8. Indépendance de $q$ : La fonction $f^{(T)}_\nu$ est montrée comme étant indépendante du nombre quantique magnétique $q$, ce qui est cohérent avec le théorème de Wigner–Eckart.

Signification
Les auteurs affirment que ce travail initie l'observation et l'application de l'ETH non-abélienne. Il fournit le premier soutien numérique approfondi pour l'hypothèse, s'étendant au-delà des études antérieures en 2D à un système 1D avec un minimum de symétries extrinsèques. En validant la NA-ETH, ce travail comble le fossé entre la thermodynamique quantique (spécifiquement l'étude des charges non-commutatives) et la physique à plusieurs corps. Les auteurs suggèrent que ces résultats permettent de futures investigations sur la manière dont les symétries non-abéliennes modifient la thermalisation, la rétention d'information dans les mémoires quantiques et les relations de dissipation-fluctuation, bien qu'ils ne proposent pas de nouvelles expériences ou applications spécifiques au-delà de ces voies théoriques.