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

Dit onderzoek toont aan dat in de isotrope spin-1/2 Heisenberg-ketting de niet-diagonale matrixelementen van lokale operatoren exponentieel afnemen met de systeemgrootte, waarbij de afname voor toestanden binnen dezelfde macrotoestand lineair is en voor verschillende macrotoestanden kwadratisch verloopt, en dat de verdelingen van deze elementen goed worden beschreven door Gumbel-distributies.

De kern

Stel je voor dat je een gigantisch, ingewikkeld labyrint hebt met miljoenen deuren. Elke deur leidt naar een andere kamer, en elke kamer vertegenwoordigt een specifieke toestand van een quantum-systeem (zoals een keten van magnetische deeltjes).

Dit artikel, geschreven door Federico Rottoli en Vincenzo Alba, onderzoekt hoe deze deuren met elkaar verbonden zijn. Ze kijken naar een heel specifiek type labyrint: de isotrope Heisenberg-keten (een rij van spin-1/2 deeltjes). Dit is een "integreerbaar" systeem, wat betekent dat het heel geordend is en veel wetten heeft die het gedrag voorspelbaar maken, in tegenstelling tot chaotische systemen die als een warboel van chaos lijken.

Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen:

1. Het Grote Doel: Waarom wordt het warm?

In de natuurkunde is een groot raadsel: hoe wordt een koud, geïsoleerd systeem (zoals een magnetische keten) op een dag warm en evenwichtig?

• De oude theorie (ETH): Voor chaotische systemen dachten wetenschappers dat de "deuren" tussen de kamers (de matrix-elementen) heel specifiek gedragen. Ze zeggen: "De meeste deuren zijn gesloten, maar degenen die open zijn, gedragen zich als een willekeurige ruis."
• Het probleem: Voor de geordende systemen (zoals onze Heisenberg-keten) wisten ze niet zeker of deze regels ook opgingen.

2. De Methode: Een magische sleutel

Om dit te onderzoeken, gebruiken de auteurs een geavanceerde wiskundige techniek genaamd Algebraïsche Bethe-Ansatz.

• De analogie: Stel je voor dat je normaal gesproken een heel labyrint moet doorlopen om te zien welke deuren open zijn. Dat zou jaren duren. De auteurs hebben echter een "magische sleutel" (de Bethe-Ansatz formules) die hen direct de plattegrond geeft. Hiermee kunnen ze berekenen hoe sterk twee kamers met elkaar verbonden zijn, zelfs voor systemen die zo groot zijn dat ze ze niet fysiek kunnen bouwen (tot wel 240 deeltjes lang!).

3. De Ontdekkingen: Twee soorten deuren

Ze keken naar twee situaties:

Situatie A: Kamers in dezelfde "buurt" (Zelfde macrotoestand)
Stel je voor dat je twee kamers kiest die beide in een drukke, warme stad liggen (dezelfde temperatuur).

• Wat ze zagen: De verbinding tussen deze deuren wordt kleiner naarmate het systeem groter wordt. Het is alsof de deur een beetje dichtklapt.
• De verrassing: De manier waarop deze deuren dichtklappen, is niet helemaal willekeurig (zoals een muntworp). De statistiek van deze deuren volgt een heel specifiek patroon, genaamd de Gumbel-verdeling.
• Vergelijking: Stel je voor dat je een berg van blokken bouwt. In een chaotisch systeem vallen de blokken willekeurig om. In dit geordende systeem vallen ze om volgens een heel specifiek, voorspelbaar patroon dat lijkt op de vorm van een bepaalde soort golf (de Gumbel-kromme). Dit betekent dat het systeem "anders" is dan de chaotische theorie voorspelde, maar het gedraagt zich toch op een manier die leidt tot thermisch evenwicht.

Situatie B: Kamers in verschillende "werelden" (Verschillende macrotoestanden)
Nu kiezen ze twee kamers die in totaal verschillende werelden liggen: één in een ijskoude pool (0 Kelvin) en één in een gloeiende oven (oneindige temperatuur).

• Wat ze zagen: Hier sluit de deur niet alleen een beetje, maar wordt hij dichtgelast. De verbinding verdwijnt veel sneller naarmate het systeem groter wordt.
• De analogie: Het is alsof je probeert te fluisteren van de ene kant van de oceaan naar de andere. De kans dat iemand je hoort, is niet alleen klein, maar wordt exponentieel kleiner naarmate de oceaan breder wordt.

4. De Belangrijkste Conclusie: Orde binnen de chaos

Het belangrijkste wat deze paper laat zien, is dat zelfs in een heel geordend, "integreerbaar" systeem (waar je zou denken dat er geen thermische wetten gelden), de deuren tussen de toestanden zich gedragen alsof het systeem toch warm wordt.

• Ze ontdekten dat de "ruis" in de verbindingen niet normaal verdeeld is (zoals een klokkromme), maar een Gumbel-verdeling volgt.
• Dit is als het verschil tussen het gooien van een dobbelsteen (chaos) en het gooien van een dobbelsteen die een beetje scheef is (geordend, maar nog steeds voorspelbaar op een andere manier).

Samenvattend in één zin:

De auteurs hebben bewezen dat zelfs in een perfect geordend quantum-systeem, de verbindingen tussen verschillende energietoestanden zo gedragen dat het systeem uiteindelijk "warm" wordt, maar dat deze verbindingen een heel specifiek, niet-willekeurig patroon volgen dat verschilt van wat we in chaotische systemen zien.

Het is alsof ze hebben ontdekt dat zelfs in een perfect georganiseerd orkest, de muziek die uit de verschillende instrumenten komt, op een heel specifieke, niet-toevallige manier samensmelt tot een harmonieus geheel.

Technische samenvatting

1. Probleemstelling en Context

Het fundamentele probleem in de niet-evenwichts kwantumveeldeeltjesfysica is te begrijpen hoe unitaire evolutie in geïsoleerde systemen leidt tot lokale equilibratie en thermalisatie. Voor chaotische (niet-integrabele) systemen wordt dit beschreven door de Eigenstate Thermalization Hypothesis (ETH). Volgens de ETH gedragen de matrixelementen van lokale operatoren $O$ tussen energie-eigenstaten $|E_i\rangle$ en $|E_j\rangle$ zich als:
$$\langle E_i | O | E_j \rangle = O(E) \delta_{ij} + f_O(\omega, E) e^{-S(E)/2} R_{ij}$$
waarbij $R_{ij}$ een willekeurige variabele is met een Gaussische verdeling.

In integrabele systemen (zoals de Heisenberg-spinketen) is de ETH echter niet direct van toepassing vanwege het bestaan van een extensief aantal behouden ladingen. Hoewel recente studies (zoals Ref. [26] voor het Lieb-Liniger-model) de schaling van diagonale en niet-diagonale matrixelementen in integrabele veldtheorieën hebben onderzocht, ontbreekt het aan een vergelijkend onderzoek voor integrabele roostermodellen (spin-ketens) die gebonden toestanden ("strings") bevatten. De aanwezigheid van deze strings maakt de numerieke berekening van matrixelementen aanzienlijk complexer door het optreden van fictieve singulariteiten.

Het doel van dit werk is om de ETH-schaling voor niet-diagonale matrixelementen in de isotrope spin-1/2 Heisenberg-keten (XXX-keten) te onderzoeken, met name de invloed van gebonden toestanden en het verschil tussen eigenstaten binnen dezelfde of verschillende thermodynamische macrotoestanden.

2. Methodologie

De auteurs maken gebruik van geavanceerde resultaten uit de Algebraïsche Bethe-Ansatz (ABA) om matrixelementen efficiënt te berekenen voor systemen die veel groter zijn dan wat met exacte diagonalisatie mogelijk is.

• Systeem: De isotrope spin-1/2 Heisenberg-keten (XXX-model) met lengte $L$ tot $240$ (voor één-spin operatoren) en $L \approx 60$ (voor twee-spin operatoren).
• Operatoren: Lokale operatoren met ondersteuning op één site ($E^{(11)}_i$) en twee naburige sites ($E^{(11)}_i E^{(22)}_{i+1}$).
• Macrotoestanden: Er wordt gewerkt met thermische macrotoestanden (Gibbs-ensemble) bij verschillende temperaturen (inclusief oneindige temperatuur en nul-temperatuur).
• Berekeningsstrategie:
  • Gebruik van de Slavnov-formule voor scalair producten tussen Bethe-toestanden.
  • Toepassing van de string-hypothese om complexe snelle getallen (rapidity) te behandelen als gebonden toestanden.
  • Regularisatie: Een cruciaal onderdeel is het zorgvuldig verwijderen van fictieve singulariteiten die ontstaan door de homogene limiet en de limiet van verdwijnende string-afwijkingen ($\delta \to 0$). Voor één-spin operatoren kan dit analytisch; voor twee-spin operatoren is dit veel complexer, waardoor de auteurs zich beperken tot eigenstaten zonder strings voor de twee-spin case.
  • Monte Carlo Sampling: Een Metropolis-algoritme wordt gebruikt om eigenstaten te bemonsteren volgens het Gibbs-gewicht $e^{-\beta E}$.

3. Belangrijkste Resultaten

A. Matrixelementen binnen dezelfde macrotoestand

Voor eigenstaten die tot dezelfde thermodynamische macrotoestand behoren (bijv. beide bij oneindige temperatuur):

• Schaling: De grootte van de niet-diagonale matrixelementen $|\langle E_i | O | E_j \rangle|$ vertoont een exponentiële afname met de systeemgrootte $L$, vergelijkbaar met de ETH maar met een logaritmische correctie:
  $$|\langle E_i | O | E_j \rangle| \propto \exp\left( -c_O L \ln L - f M_{ij} L \right)$$
  Dit bevestigt dat de aanwezigheid van gebonden toestanden (strings) het kwalitatieve beeld van Ref. [26] niet verandert.
• Verdeling (Statistics): De verdelingsfunctie (PDF) van de logaritme van de matrixelementen, $M_{ij} = \ln |\langle E_i | O | E_j \rangle|^2$, wordt uitstekend beschreven door een Gumbel-verdeling.
  • Dit staat in schril contrast met het Lieb-Liniger-model (Ref. [26]), waar een Fréchet-verdeling werd gevonden.
  • De parameters van de Gumbel-verdeling zijn onafhankelijk van $L$.
• Niet-Gaussisch: De verdeling is niet-Gaussisch. Dit wordt bevestigd door de verhouding $\Gamma = \overline{|O_{ij}|^2} / (\overline{|O_{ij}|})^2$ te analyseren. Voor een Gaussische verdeling zou $\Gamma = \pi/2$ moeten zijn; de auteurs vinden echter dat $\Gamma \sim 10^5$ en divergeert met $L$, wat consistent is met de Gumbel-statistiek.

B. Matrixelementen tussen verschillende macrotoestanden

Voor eigenstaten die tot verschillende macrotoestanden behoren (bijv. grondtoestand vs. oneindige temperatuur):

• Snellere Afname: De matrixelementen vervallen veel sneller dan in het geval van dezelfde macrotoestand:
  $$|\langle E_i | O | E_j \rangle| \propto \exp(-d_O L^2)$$
• Verdeling: Ook hier volgt de verdeling van de logaritme van de matrixelementen een Gumbel-verdeling, hoewel de schaling met $L^2$ verschilt.

C. Vergelijking Eén-spin vs. Twee-spin Operatoren

• De resultaten voor één-spin operatoren (tot $L \approx 240$) en twee-spin operatoren (tot $L \approx 60$) tonen hetzelfde kwalitatieve gedrag.
• De beperking in systeemgrootte voor twee-spin operatoren komt door de exponentiële toename van de rekentijd met de ondersteuning van de operator, maar de data bevestigen dat het scenario robuust is voor lokale operatoren met bredere ondersteuning.

4. Bijdragen en Significatie

• Uitbreiding naar Roostermodellen: Het artikel vult een belangrijke lacune op door de ETH-schaling in integrabele systemen te testen voor roostermodellen (spin-ketens) in plaats van alleen voor continue veldtheorieën (Lieb-Liniger).
• Rol van Gebonden Toestanden: Het werk demonstreert dat de aanwezigheid van "strings" (gebonden toestanden), die experimenteel relevant zijn, de fundamentele schaling van niet-diagonale matrixelementen niet verandert, maar wel de statistische verdeling beïnvloedt.
• Statistische Universiteit: De ontdekking dat de verdeling Gumbel is (in plaats van Fréchet of Gaussisch) voor de XXX-keten suggereert dat de statistiek van matrixelementen in integrabele systemen fundamenteel verschilt van die in chaotische systemen en mogelijk afhankelijk is van het specifieke model.
• Methodologische Doorbraak: De succesvolle toepassing van ABA-formules voor matrixelementen van operatoren met ondersteuning op twee sites in de aanwezigheid van strings (of het zorgvuldig vermijden daarvan) opent de weg voor verdere studies van grootschalige integrabele systemen.

Conclusie

De auteurs concluderen dat hoewel de exponentiële afname van niet-diagonale matrixelementen in integrabele systemen lijkt op die van chaotische systemen (ETH), de statistiek van deze elementen fundamenteel anders is. In plaats van Gaussische fluctuaties, vertonen ze een Gumbel-verdeling. Dit onderstreept dat integrabiliteit niet alleen de dynamica beïnvloedt, maar ook de fundamentele statistische eigenschappen van de eigenstaten in het spectrum. Dit heeft implicaties voor het begrijpen van relaxatiedynamica en de toepasbaarheid van vrije waarschijnlijkheidstheorie (Free Probability Theory) op integrabele modellen.