Operator delocalization in disordered spin chains via exact MPO marginals

Titel: Operator-delokalisatie in ongeordende spin-ketens via exacte MPO-marginaalverdelingen

Auteurs: Jonnathan Pineda, Mario Collura, Gianluca Passarelli, Procolo Lucignano, Davide Rossini, en Angelo Russomanno.

1. Probleemstelling en Context

Het artikel onderzoekt hoe kwantuminformatie verspreidt (scrambling) en hoe operatoren delokaliseren in kwantumsystemen, met name in de context van veeldeeltjeslokalisatie (MBL) en Anderson-lokalisatie.

• Achtergrond: In thermiserende systemen worden operatoren gaandeweg de tijd steeds meer gedelokaliseerd, wat leidt tot scrambling. Dit wordt vaak gemeten via Out-of-Time-Order Correlators (OTOC), die een "lichtkegel" van informatieverspreiding tonen. In chaotische systemen is dit ballistisch, terwijl het in geïntegreerde systemen sub-diffusief is.
• Het MBL-fenomeen: In sterk ongeordende, interagerende systemen (MBL) wordt thermalisatie onderdrukt. In plaats van een lineaire lichtkegel, verspreidt informatie zich logaritmisch in de tijd ($t \sim e^r$ of $r \sim \ln t$).
• De uitdaging: Bestaande methoden om operator-complexiteit te meten (zoals de stabilizer Rényi-entropy of operator-magic) vereisen vaak stochastisch sampling (Monte Carlo), wat inefficiënt is en exponentieel schaalt met het systeemgrootte. Er is behoefte aan een methode die de volledige waarschijnlijkheidsverdeling van operator-groei exact en efficiënt kan berekenen zonder sampling.

2. Methodologie

De auteurs introduceren twee nieuwe, complementaire maatstaven voor operator-complexiteit en ontwikkelen een numeriek raamwerk om deze exact te berekenen.

A. Nieuwe Maatstaven: Operator-massa en Operator-lengte

In plaats van alleen te kijken naar de verspreiding van informatie, definiëren de auteurs twee grootheden gebaseerd op de expansie van een tijd-geëvolueerde operator $\hat{O}(t)$ in de Pauli-basis:

1. Operator-massa ($m$): Het aantal sites waarop de operator niet-triviaal werkt (het aantal niet-identiteit Pauli-matrices in een Pauli-string). Dit kwantificeert het aantal betrokken deeltjes.
2. Operator-lengte ($h$): De positie van het meest rechtse niet-identiteit element in de Pauli-string. Dit kwantificeert de ruimtelijke uitbreiding (hoe ver de operator reikt).

Deze grootheden worden berekend als het gemiddelde over de waarschijnlijkheidsverdeling $P(\hat{Q}) = |A_{\hat{Q}}|^2$, waarbij $A_{\hat{Q}}$ de coëfficiënten zijn van de expansie in de genormaliseerde Pauli-basis.

B. Numerieke Implementatie: MPS/MPO Raamwerk

• Operator als Toestand: De auteurs gebruiken de Choi-Jamiołkowski-isomorfie om een operator $\hat{O}$ af te beelden op een vector $|\hat{O}\rangle$ in een verdubbelde Hilbertruimte.
• Matrix Product States (MPS): In plaats van de volledige operator te opslaan, wordt deze gerepresenteerd als een MPS. Dit maakt het mogelijk om de tijdsevolutie van grote systemen (tot $L=32$ sites) te simuleren.
• Exacte Berekening: Een cruciale innovatie is dat de auteurs laten zien dat de verdelingen voor massa en lengte exact en deterministisch kunnen worden berekend binnen het MPS-formaat.
  • De lengte wordt berekend via een recursieve relatie die gebaseerd is op de tweede-orde Rényi-entropy van de gereduceerde dichtheidsmatrix.
  • De massa wordt berekend door een diagonale superoperator te contracteren met de MPS.
• Voordelen: Deze methode vereist geen stochastisch sampling en vermijdt de exponentiële kostprijs van volledige operator-tomografie. In het MBL-regime groeit de benodigde bond-dimension slechts lineair met de tijd, wat lange tijdsimulaties mogelijk maakt.

C. Experimenteel Protocol

Het artikel schetst ook een experimenteel protocol om deze grootheden te meten op kwantumsimulatoren (bijv. supergeleidende qubits of gevangen ionen):

• Gebruik van een verdubbelde systeem-ancilla staat.
• Toepassing van "klassieke schaduwen" (classical shadows) met Bell-paar metingen (geïmplementeerd via lokale Clifford-gates) om de marginaalverdelingen van massa en lengte direct uit experimentele shots te reconstrueren.

3. Belangrijkste Resultaten

De auteurs passen hun methode toe op een ongeordende XXZ spin-1/2 keten met willekeurige velden $h_j$ en interactie $\Delta$.

1. Verschil tussen Anderson en MBL:

   • Niet-interagerend geval ($\Delta = 0$): Het systeem is een Anderson-geïsoleerd vrije fermionenmodel. De operator-massa, -lengte en entanglement-entropy satureren snel tot een constante waarde. Er is geen echte delokalisatie.
   • Interagerend geval ($\Delta \neq 0$): Zelfs voor zeer zwakke interacties, vertonen alle grootheden een robuste logaritmische groei in de tijd ($h(t) \sim \ln t$). Dit bevestigt het bestaan van een logaritmische lichtkegel in het MBL-regime.

2. Koppeling met $\ell$-bit model:

   • De numerieke resultaten komen perfect overeen met een analytisch oplosbaar $\ell$-bit model (lokaal geïntegreerde bewegingsintegralen).
   • In dit model wordt de logaritmische groei verklaard door de exponentiell zwakke interacties tussen de $\ell$-bits, wat leidt tot een tijdschaal $t_r \sim e^{r/\xi}$ en dus $r(t) \sim \xi \ln t$.

3. Systeemgrootte en Disorder:

   • De logaritmische groei is robuust voor verschillende systeemgroottes ($L$) en disorder-strengths ($W$).
   • Bij versterking van de interactie ($\Delta > 0$) verdwijnt het saturatiegedrag van het Anderson-geval abrupt; de systemen vertonen een niet-perturbatieve overgang naar het MBL-gedrag.
   • De tijd-gegemiddelde grootheden, geschaald met $\ln L$, collapse tot een universele kromme, wat een sterke numerieke bevestiging is van de MBL-logaritmische lichtkegel.

4. Bijdragen en Significatie

• Nieuwe Observabelen: De introductie van "operator-lengte" als een complementaire maatstaf voor "operator-massa" biedt een scherpere kijk op de ruimtelijke structuur van operator-groei.
• Efficiëntie en Nauwkeurigheid: De methode biedt een manier om de volledige waarschijnlijkheidsverdeling van operator-complexiteit exact te berekenen zonder sampling, wat een doorbraak is ten opzichte van eerdere benaderingen die afhankelijk waren van stochastische sampling of beperkt waren tot kleine systemen.
• MBL Karakterisering: Het werk levert overtuigend numeriek bewijs dat de logaritmische lichtkegel een fundamenteel kenmerk is van het MBL-fase, zelfs bij zeer zwakke interacties, en onderscheidt dit duidelijk van Anderson-lokalisatie.
• Experimentele Toepasbaarheid: Door het koppelen van de theoretische berekeningen aan een haalbaar experimenteel protocol (via klassieke schaduwen en Bell-metingen), opent dit onderzoek de weg voor directe experimentele verificatie van operator-delokalisatie in huidige kwantumplatforms.

Conclusie:
Dit artikel biedt een krachtig, exact en experimenteel haalbaar raamwerk om de dynamiek van kwantuminformatie in ongeordende systemen te bestuderen. Het bevestigt dat interacties essentieel zijn voor de langzame, logaritmische delokalisatie in MBL-systemen en dat tensor-netwerk methoden (MPS) ideaal zijn om deze complexe dynamiek te simuleren.