Operator delocalization in disordered spin chains via exact MPO marginals

En introduisant la notion de « longueur d'opérateur » comme mesure complémentaire de la complexité, cette étude démontre que le cadre des opérateurs de produit matriciel (MPO) permet de calculer exactement la localisation des opérateurs dans les chaînes de spins désordonnées, révélant une saturation rapide dans le cas d'Anderson et une croissance logarithmique robuste dans le régime de localisation à plusieurs corps (MBL).

L'essentiel

🧠 Le Grand Voyage de l'Information : Quand le Chaos Rencontre le Désordre

Imaginez que vous avez une pièce de monnaie (une information) posée sur la première case d'un très long tapis roulant composé de milliers de cases (une chaîne d'atomes ou de spins).

Dans un monde normal et ordonné, si vous lancez cette pièce, elle va rouler, rebondir et se mélanger à tout le tapis très vite. C'est ce qu'on appelle la thermalisation : l'information se disperse partout, devenant impossible à retrouver à sa place d'origine. C'est comme une goutte d'encre qui se diffuse rapidement dans un verre d'eau.

Mais, que se passe-t-il si le tapis est très désordonné (rempli de nids-de-poule, de bosses, de trous aléatoires) ?

1. Les Deux Types de "Blocages"

Les chercheurs étudient deux scénarios où l'information a du mal à bouger :

• Le Cas "Libre" (Anderson) : Imaginez un tapis désordonné mais sans frottement (pas d'interactions entre les particules). Si vous lancez la pièce, elle va vite se coincer dans un trou et rester là. Elle ne bouge plus. C'est ce qu'on appelle la localisation d'Anderson. L'information est prisonnière.
• Le Cas "Interagissant" (MBL - Localisation à Many-Body) : C'est là que ça devient intéressant. Maintenant, imaginez que les particules sur le tapis se parlent entre elles (elles interagissent). Même avec un tapis très désordonné, on pensait que cette conversation ferait bouger la pièce.
  • La découverte : Même avec des interactions, si le désordre est assez fort, la pièce ne bouge pas vite. Elle avance, mais extrêmement lentement, comme une tortue qui grimpe une montagne de sable. C'est la Localisation à Many-Body (MBL).

2. Les Nouveaux Outils de Mesure : La "Masse" et la "Longueur"

Pour mesurer comment cette information se déplace, les auteurs ont inventé deux nouvelles règles du jeu, basées sur la façon dont l'information est "écrite" en langage quantique (le langage des "Pauli strings", qui sont comme des mots composés de lettres spéciales).

• La Masse de l'Opérateur (Operator Mass) :
  Imaginez que votre information est un mot. La "masse", c'est le nombre de lettres dans ce mot.

  • Au début, le mot est court : "A". (Masse = 1).
  • En se déplaçant, il se transforme en "ABCDEF...". (Masse augmente).
  • Plus la masse est grande, plus l'information est complexe et étalée sur beaucoup de particules.

• La Longueur de l'Opérateur (Operator Length) :
  C'est la distance entre le début et la fin du mot.

  • Si le mot est "A...Z" (avec des lettres vides au milieu), la longueur est la distance jusqu'à la lettre "Z".
  • Cela nous dit jusqu'où l'information a physiquement voyagé le long du tapis.

3. La Révolution : Regarder sans "Deviner"

Avant, pour mesurer ces choses, il fallait faire des millions de devinettes (échantillonnage aléatoire), un peu comme essayer de deviner la couleur d'un océan en regardant une seule goutte d'eau à la fois. C'était long et imprécis.

Ces chercheurs ont développé une méthode (basée sur les "réseaux de tenseurs" ou MPO) qui permet de voir tout le tapis d'un coup, exactement et sans deviner. C'est comme passer d'une lampe torche à un projecteur géant qui éclaire tout instantanément.

4. Les Résultats : La Croissance Logarithmique

Voici ce qu'ils ont découvert en regardant le tapis se dérouler dans le temps :

• Sans interactions (Cas libre) : L'information se coince immédiatement. La masse et la longueur s'arrêtent de croître très vite. C'est un mur.
• Avec interactions (Cas MBL) : L'information commence à bouger, mais très bizarrement. Elle ne suit pas une ligne droite (comme une voiture qui accélère). Elle suit une courbe très plate : elle grandit comme le logarithme du temps.
  • L'analogie : Imaginez que pour doubler la distance parcourue, vous devez attendre non pas deux fois plus de temps, mais beaucoup plus de temps. C'est une progression ultra-lente, mais qui ne s'arrête jamais vraiment. C'est ce qu'on appelle le "cône de lumière logarithmique".

5. Pourquoi est-ce important ?

• Pour la science : Cela confirme que même dans un système désordonné et interactif, l'information ne disparaît pas, elle se "fige" très lentement. Cela aide à comprendre pourquoi certains matériaux ne deviennent pas chauds (ne thermalisent pas) même après des milliards d'années.
• Pour le futur (Ordinateurs Quantiques) : Si vous voulez construire un ordinateur quantique, vous voulez que l'information reste stable et ne se mélange pas trop vite avec l'environnement. Ce phénomène de "localisation lente" pourrait être un moyen de protéger l'information quantique contre les erreurs.
• Pour les expériences : Les auteurs montrent aussi comment mesurer tout cela dans de vrais laboratoires (avec des atomes froids ou des circuits supraconducteurs) en utilisant des techniques de "ombres classiques" (une méthode de photo rapide et intelligente).

En Résumé

C'est comme si vous observiez une foule dans une pièce remplie d'obstacles.

• Si les gens ne se parlent pas, ils restent bloqués à leur place.
• S'ils se parlent, ils commencent à bouger, mais très lentement, comme s'ils avançaient dans du miel épais.
• Les chercheurs ont inventé une caméra spéciale pour compter exactement combien de personnes ont bougé et jusqu'où, prouvant que cette progression lente est une règle fondamentale de la nature quantique désordonnée.

Résumé technique

1. Problématique et Contexte

L'article s'intéresse à la dynamique de la délocalisation des opérateurs quantiques dans les systèmes de spins désordonnés en une dimension, en particulier dans le contexte de la localisation à plusieurs corps (MBL - Many-Body Localization).

• Le défi : Comprendre comment l'information quantique se propage (scrambling) et comment les opérateurs locaux deviennent non locaux au cours du temps. Dans les systèmes chaotiques, cette propagation est rapide (balistique ou diffusive), tandis que dans les systèmes MBL, elle est extrêmement lente.
• Limites des méthodes existantes : La mesure standard du scrambling, le corrélateur hors ordre temporel (OTOC), fournit des informations spatio-temporelles mais ne dépend d'aucune base d'opérateurs spécifique. D'autres approches basées sur la base de Pauli (comme l'entropie de Rényi de stabilisateur ou SRE) nécessitent un échantillonnage stochastique coûteux, dont le nombre d'échantillons requis peut croître exponentiellement avec la taille du système.
• Objectif : Développer des mesures complémentaires de la complexité des opérateurs (la masse et la longueur de l'opérateur) qui puissent être calculées exactement et efficacement sans échantillonnage stochastique, en utilisant des réseaux de tenseurs.

2. Méthodologie

Les auteurs proposent un cadre théorique et numérique basé sur l'utilisation d'opérateurs sous forme de produit matriciel (MPO - Matrix Product Operators) dans l'espace de Hilbert des opérateurs.

A. Définition des observables clés

En développant un opérateur évolué $\hat{O}(t)$ dans la base des chaînes de Pauli normalisées $\{\hat{P}_\mu\}$, les auteurs définissent deux grandeurs :

1. Masse de l'opérateur ($m(t)$) : Le nombre moyen de sites sur lesquels l'opérateur agit de manière non triviale (c'est-à-dire le nombre de matrices de Pauli différentes de l'identité dans la chaîne). Elle quantifie la densité de l'opérateur.
2. Longueur de l'opérateur ($h(t)$) : La position du site le plus à droite où l'opérateur agit de manière non triviale. Elle caractérise l'étendue spatiale du support de l'opérateur.

Contrairement à la masse, qui compte les sites actifs, la longueur mesure la propagation du "front" de l'opérateur.

B. Implémentation via MPO et Réseaux de Tenseurs

• Représentation Choi-Liouville : L'opérateur $\hat{O}$ est mappé sur un vecteur d'état $|\hat{O}\rangle$ dans un espace de Hilbert doublé (système + ancilla).
• Calcul exact des distributions : En représentant $|\hat{O}\rangle$ sous forme d'un MPS (Matrix Product State), les auteurs montrent que les distributions de probabilité complètes de la masse et de la longueur peuvent être calculées déterministiquement.
  • La longueur est obtenue via des relations de récurrence impliquant les entropies de Rényi d'ordre 2 tronquées.
  • La masse est calculée comme la valeur moyenne d'un super-opérateur diagonal agissant sur l'état MPS.
• Avantage computationnel : Cette méthode évite le bruit statistique de l'échantillonnage. Dans le régime MBL, la dimension de liaison (bond dimension) nécessaire pour représenter l'opérateur ne croît que linéairement avec le temps (car l'entropie d'intrication de l'opérateur croît logarithmiquement), rendant les simulations à long temps et sur de grands systèmes (jusqu'à $L=32$) possibles.

C. Modèle Physique

Les auteurs étudient la chaîne de spins-1/2 XXZ désordonnée avec :

• Terme de saut $J$ (fixé à 1).
• Interaction de densité-densité $\Delta$.
• Champ désordonné $h_j$ uniformément distribué dans $[-W, W]$.
• Ils comparent le cas non-interactif ($\Delta=0$, modèle de fermions libres localisés d'Anderson) et le cas interactif ($\Delta \neq 0$, régime MBL).

3. Résultats Clés

Les simulations numériques (ED pour $L=12$ et MPO pour $L$ jusqu'à 32) révèlent des comportements distincts selon le régime :

A. Cas Non-Interactif ($\Delta = 0$)

• Comportement : La masse, la longueur et l'entropie d'intrication de l'opérateur saturent rapidement à des valeurs constantes.
• Interprétation : Le système correspond à un modèle d'Anderson localisé. Les états propres à une particule sont localisés, empêchant toute propagation significative de l'opérateur au-delà de la longueur de localisation. Il n'y a pas de scrambling.

B. Cas Interactif ($\Delta \neq 0$, Régime MBL)

• Croissance Logarithmique : Pour un désordre fort ($W \gtrsim 5$), la longueur $h(t)$, la masse $m(t)$ et l'entropie d'intrication de l'opérateur $S_{op}(t)$ présentent une croissance logarithmique en fonction du temps ($h(t) \sim a \ln t + b$).
• Robustesse : Ce comportement persiste même pour des interactions très faibles ($\Delta = 0.2, 0.4$). La transition entre le comportement saturé (Anderson) et la croissance logarithmique (MBL) est non perturbative : une interaction infinitésimale suffit à induire une délocalisation lente mais illimitée.
• Cône de lumière logarithmique : Cette croissance logarithmique confirme l'existence d'un "cône de lumière" logarithmique pour la propagation des corrélations quantiques, caractéristique du régime MBL.
• Efficacité des simulations : La croissance logarithmique de l'entropie d'intrication de l'opérateur garantit que la dimension de liaison du MPO reste gérable, permettant des simulations sur des temps très longs ($t \sim 10^5$).

C. Validation par le Modèle $\ell$-bit

Les résultats numériques sont parfaitement reproduits par un modèle analytique simple basé sur des intégrales de mouvement localisées ($\ell$-bits).

• Dans ce modèle, les interactions entre les $\ell$-bits décroissent exponentiellement avec la distance ($J_{jl} \sim e^{-\kappa |j-l|}$).
• La résolution analytique montre que le temps nécessaire pour qu'un opérateur atteigne une distance $r$ est $t_r \sim e^{r/\xi}$, ce qui implique directement une croissance logarithmique de la distance en fonction du temps : $r(t) \sim \xi \ln t$.

D. Analyse de la taille du système et du désordre

• L'analyse des moyennes temporelles sur une fenêtre proportionnelle à la taille du système $L$ montre un effondrement des données lorsque l'on trace $h/\ln L$ en fonction du désordre $W$. Cela confirme la loi d'échelle $h(t) \sim \xi(W) \ln t$.
• Une dépendance non monotone est observée pour la masse en fonction du désordre, avec un pic à un désordre intermédiaire, probablement lié à la nature intégrable du système à $W=0$.

4. Contributions et Signification

Contributions Techniques

1. Nouvelles Observables : Introduction de la "longueur de l'opérateur" comme mesure complémentaire à la "masse" pour quantifier la délocalisation spatiale.
2. Algorithme Exact : Développement d'une méthode basée sur les MPO pour calculer exactement les distributions de probabilité de ces grandeurs, éliminant le besoin d'échantillonnage stochastique et réduisant la complexité computationnelle de manière drastique.
3. Preuve Numérique de la Dynamique MBL : Fourniture de preuves numériques solides que la croissance logarithmique de la complexité des opérateurs est une signature universelle du régime MBL, même pour des interactions faibles.

Signification Physique

• Distinction Anderson vs MBL : L'article clarifie la différence fondamentale entre la localisation d'Anderson (où la délocalisation est absente) et la localisation MBL (où elle est lente mais présente), soulignant le rôle crucial des interactions pour briser la localisation stricte.
• Lien avec le Scrambling : La croissance logarithmique de la longueur et de la masse corrobore le fait que le scrambling de l'information dans les systèmes MBL est extrêmement lent, cohérent avec la propagation logarithmique des corrélations.

Perspectives Expérimentales

Les auteurs proposent un protocole expérimental réalisable sur les plateformes de simulation quantique actuelles (qubits supraconducteurs, ions piégés, atomes de Rydberg) :

• Utilisation de l'isomorphisme Choi-Jamiołkowski pour préparer un état système-ancilla.
• Évolution temporelle contrôlée (avant et arrière).
• Mesures basées sur les ombres classiques (classical shadows) et des mesures de paires de Bell (réalisables via des portes Clifford locales) pour reconstruire les distributions de masse et de longueur sans tomographie complète.

Conclusion

Cet article établit que la longueur et la masse des opérateurs, calculées via des techniques de réseaux de tenseurs exacts, sont des sondes puissantes et efficaces de la dynamique de délocalisation. Ils confirment que dans les systèmes désordonnés en interaction, la complexité de l'opérateur croît logarithmiquement avec le temps, une signature robuste du régime MBL qui persiste même pour des interactions infinitésimales, contrairement au cas non-interactif où l'opérateur reste strictement localisé.