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

Die Studie untersucht mittels des Algebraischen Bethe-Ansatzes die off-diagonalen Matrixelemente lokaler Operatoren in der integrablen XXX-Spin-Kette und zeigt, dass diese für Eigenzustände innerhalb desselben thermodynamischen Makrozustands exponentiell mit der Systemgröße abfallen und Gumbel-Verteilungen folgen, während sie für Zustände unterschiedlicher Makrozustände schneller abklingen.

Das Wesentliche

Das große Rätsel: Warum wird das Chaos ruhig?

Stell dir vor, du hast einen riesigen, isolierten Raum voller winziger, tanzender Quanten-Teilchen (wie eine Kugelbahn, auf der Millionen Kugeln wild durcheinanderrollen). Normalerweise erwarten wir, dass dieses Chaos sich nie beruhigt. Aber in der echten Welt sehen wir, dass sich Systeme oft "beruhigen" und eine Art Gleichgewichtszustand (Thermodynamik) erreichen, als hätten sie einen Thermostat.

Die Frage ist: Wie passiert das, wenn die Gesetze der Quantenmechanik eigentlich streng reversibel sind und nichts "vergessen" wird?

Die Antwort darauf ist die Eigenstate Thermalization Hypothesis (ETH). Sie besagt im Grunde: "Jeder einzelne Tanzschritt (Zustand) in diesem Chaos enthält bereits die Information über das gesamte Chaos." Wenn du nur einen kleinen Teil des Raums ansiehst, sieht er so aus, als wäre er im Gleichgewicht, weil der Rest des Raums wie ein "Bad" wirkt, das ihn thermalisiert.

Der spezielle Fall: Die perfekten Tänzer (Integrierbare Systeme)

Die ETH funktioniert super für chaotische Systeme (wie ein wilder Moloch). Aber es gibt eine Ausnahme: Integrierbare Systeme. Das sind wie ein Orchester, bei dem jeder Musiker eine eigene, unveränderliche Melodie spielt und sich nie mit den anderen vermischt. In der Physik nennen wir diese Systeme "integrierbar" (wie die hier untersuchte Heisenberg-Kette).

Früher dachte man: "Da gibt es keine Thermalisierung, das Orchester spielt ewig weiter." Aber die Autoren dieser Studie haben genauer hingeschaut. Sie wollten wissen: Was passiert mit den Übergängen zwischen den verschiedenen Tanzschritten?

Die Metapher: Der geheime Code (Matrixelemente)

Stell dir vor, jedes Teilchen im System hat einen "Pass" (seinen Energiezustand). Um von einem Zustand A zu einem Zustand B zu wechseln, muss ein "Geheimsignal" (ein sogenanntes Matrixelement) existieren.

• Bei chaotischen Systemen: Diese Signale sind wie zufälliges Rauschen. Sie sind sehr klein und fallen mit der Größe des Systems exponentiell ab.
• Bei diesem neuen Studium: Die Autoren haben untersucht, wie diese Signale in den "perfekten" Orchestern (integrierbaren Systemen) aussehen.

Sie haben zwei Szenarien getestet:

1. Das "Familientreffen" (Zustände im selben Makrozustand)

Stell dir vor, du vergleichst zwei Tänzer, die beide zur gleichen "Partei" gehören (gleiche Temperatur, gleiche Energie).

• Ergebnis: Die Signale zwischen ihnen sind klein, aber nicht winzig. Sie fallen exponentiell ab, wenn das System größer wird.
• Die Überraschung: Wenn man die Verteilung dieser Signale aufzeichnet, sieht sie nicht wie eine normale Glockenkurve (Gauß) aus, wie man es von Zufall erwartet. Stattdessen folgt sie einer Gumbel-Verteilung.
• Die Analogie: Stell dir vor, du wirfst Tausende von Würfeln. Bei einem chaotischen System würdest du eine perfekte Glockenkurve der Ergebnisse bekommen. Bei diesem System würdest du jedoch feststellen, dass die Würfel eine seltsame Vorliebe für bestimmte extreme Werte haben – wie ein Würfel, der öfter eine 6 wirft als statistisch erwartet, aber nur in einer sehr spezifischen, mathematisch beschreibbaren Weise.

2. Das "Fremden-Treffen" (Zustände in verschiedenen Makrozuständen)

Jetzt vergleichen wir zwei Tänzer, die zu völlig verschiedenen Partys gehören (z.B. einer bei 0 Grad Kelvin, der andere bei unendlich heiß).

• Ergebnis: Hier sind die Signale noch viel kleiner. Sie fallen so schnell ab, dass sie fast sofort verschwinden.
• Die Analogie: Es ist, als würde man versuchen, ein Gespräch zwischen zwei Leuten aufzuzeichnen, die auf gegenüberliegenden Seiten der Erde stehen und durch einen dicken Betonwall getrennt sind. Die Wahrscheinlichkeit, dass sie sich verstehen, ist so gering, dass sie praktisch null ist.

Warum ist das wichtig? (Die "Strings" und die Magie)

Das Besondere an diesem Papier ist, dass es ein System untersucht, in dem die Teilchen "gebundene Zustände" bilden können. Man nennt diese Strings.

• Die Analogie: Stell dir vor, die Teilchen sind nicht nur einzelne Kugeln, sondern manchmal schnüren sie sich zu kleinen Gruppen zusammen (wie Perlen auf einer Schnur). Diese Gruppen verhalten sich wie ein einziges Teilchen.
• Früher war es unmöglich, die Mathematik für diese Gruppen zu lösen, weil die Formeln "zerbrachen" (Singularitäten). Die Autoren haben jedoch eine hochmoderne mathematische Technik (Algebraischer Bethe-Ansatz) verwendet, um diese "zerbrochenen" Formeln zu reparieren und die Zahlen für riesige Systeme (bis zu 240 Teilchen) zu berechnen. Das ist wie der Versuch, ein Puzzle mit 240 Teilen zu lösen, während andere nur 20 schaffen.

Das Fazit in einem Satz

Selbst in einem perfekten, geordneten Quantensystem (dem Orchester), das nicht chaotisch ist, gibt es eine Art "Thermalisierung", aber sie funktioniert anders als im Chaos: Die Übergänge zwischen Zuständen sind extrem klein und folgen einer sehr spezifischen, nicht-zufälligen Statistik (Gumbel), die zeigt, dass die "Strings" (die gebundenen Teilchen) das Verhalten des Systems tiefgreifend verändern.

Zusammenfassend: Die Autoren haben bewiesen, dass selbst in der "strikten Ordnung" der Quantenwelt ein verstecktes Muster existiert, das erklärt, wie sich diese Systeme verhalten, wenn man sie beobachtet – und dieses Muster ist so elegant und seltsam wie ein Gumbel-Würfelwurf.

Technische Zusammenfassung

Titel: Eigenstate Thermalization Hypothesis (ETH) für nicht-diagonale Matrixelemente in integrablen Spin-Ketten

Autoren: Federico Rottoli und Vincenzo Alba (Universität Pisa)

---

1. Problemstellung und Motivation

Ein zentrales Problem der Quanten-Vielteilchenphysik ist das Verständnis, wie isolierte Systeme durch unitäre Evolution lokal ins Gleichgewicht kommen und thermalisieren. Für chaotische (nicht-integrable) Systeme wird dies durch die Eigenstate Thermalization Hypothesis (ETH) beschrieben. Die ETH postuliert, dass Matrixelemente eines lokalen Operators $O$ zwischen Energieeigenzuständen $|E_i\rangle$ und $|E_j\rangle$ die Form
$$\langle E_i | O | E_j \rangle = O(E) \delta_{ij} + f_O(\omega, E) e^{-S(E)/2} R_{ij}$$
annehmen, wobei $R_{ij}$ eine Gauß-verteilte Zufallsvariable ist.

In integrablen Systemen (wie dem Heisenberg-Modell) führt die Existenz einer extensiven Anzahl an Erhaltungsgrößen typischerweise zu Verletzungen der Standard-ETH. Während das Verhalten der diagonalen Matrixelemente in integrablen Systemen bereits untersucht wurde (Abklingen als $L^{-1}$ statt exponentiell), war das Verhalten der nicht-diagonalen Matrixelemente in Gittermodellen (Spin-Ketten) mit gebundenen Zuständen („Strings") noch unklar.

Bisherige Arbeiten (z. B. Ref. [26]) untersuchten das Lieb-Liniger-Modell (ein kontinuierliches Feldtheorie-Modell ohne gebundene Zustände im repulsiven Fall) und fanden:

• Für Eigenzustände im selben makroskopischen Zustand: Exponentielles Abklingen mit logarithmischen Korrekturen ($\sim e^{-L \ln L}$).
• Für Eigenzustände in verschiedenen makroskopischen Zuständen: Schnelleres Abklingen ($\sim e^{-L^2}$).
• Die Statistik der Matrixelemente folgte einer Fréchet-Verteilung.

Die vorliegende Arbeit untersucht, ob dieses Szenario auf integrable Gittermodelle (insbesondere die isotrope Spin-1/2-Heisenberg-Kette, XXX-Modell) übertragbar ist, die gebundene Zustände (Strings) aufweisen. Die Berechnung von Matrixelementen in solchen Systemen ist aufgrund von Singularitäten in den Formeln der Algebraischen Bethe-Ansatz (ABA) extrem schwierig.

---

2. Methodik

Die Autoren verwenden Zustandsformeln des Algebraischen Bethe-Ansatzes (ABA), um Matrixelemente lokaler Operatoren effizient zu berechnen.

• Modell: Isotrope Spin-1/2-Heisenberg-Kette (XXX-Modell) mit Hamiltonian $H = \frac{J}{4} \sum \vec{\sigma}_i \cdot \vec{\sigma}_{i+1}$.
• Operatoren:
  • Ein-Spin-Operatoren (z. B. $E^{(11)}_i$).
  • Zwei-Spin-Operatoren (z. B. $E^{(11)}_i E^{(22)}_{i+1}$).
• Herausforderung: Die ABA-Formeln enthalten fiktive Singularitäten, die durch den homogenen Limes ($\xi_i \to \eta/2$) und das Vorhandensein von Strings (gebundene Zustände) entstehen.
• Lösungsstrategie:
  • Für Ein-Spin-Operatoren können die Singularitäten analytisch entfernt werden, auch bei Vorhandensein von Strings. Dies ermöglicht die Untersuchung von Systemgrößen bis $L \approx 500$.
  • Für Zwei-Spin-Operatoren ist die gleichzeitige Behandlung des homogenen Limes und des Limes verschwindender String-Abweichungen sehr komplex. Daher beschränken sich die Autoren hier auf Eigenzustände ohne gebundene Zustände (keine Strings) und Systemgrößen bis $L \approx 60$.
• Stichprobenziehung: Um thermische makroskopische Zustände zu simulieren, wird ein Metropolis-Algorithmus verwendet, um Eigenzustände gemäß der Gibbs-Verteilung $e^{-\beta E}$ zu sampeln.
• Analyse: Untersucht wird die Verteilungsfunktion (PDF) des Logarithmus der Matrixelemente: $M_{ij} = \ln |\langle E_i | O | E_j \rangle|^2$.

---

3. Wichtige Ergebnisse

A. Skalierung und Abklingen (Same Macrostate)

Für Eigenzustände, die demselben thermischen makroskopischen Zustand angehören (z. B. unendliche Temperatur):

• Die nicht-diagonalen Matrixelemente skalieren wie $|M_{ij}| \sim L$.
• Das Abklingen erfolgt exponentiell mit der Systemgröße $L$, jedoch mit einer logarithmischen Korrektur, analog zu Ref. [26]:
  $$|M_{ij}| \approx c L \ln L + c_0 L$$
• Dies zeigt, dass das Vorhandensein von gebundenen Zuständen (Strings) das qualitative Skalierungsverhalten nicht verändert.

B. Skalierung und Abklingen (Different Macrostates)

Für Eigenzustände aus verschiedenen makroskopischen Zuständen (z. B. Grundzustand vs. unendliche Temperatur):

• Die Matrixelemente zeigen ein deutlich schnelleres Abklingen:
  $$|M_{ij}| \sim e^{-d L^2}$$
• Dies bestätigt das in Ref. [26] für das Lieb-Liniger-Modell gefundene Verhalten auch für Spin-Ketten.

C. Statistische Verteilung (Der Kernbeitrag)

Ein entscheidender Unterschied zu chaotischen Systemen und früheren Ergebnissen im Lieb-Liniger-Modell wurde in der Statistik der Matrixelemente gefunden:

• Gumbel-Verteilung: Die Wahrscheinlichkeitsdichtefunktion (PDF) von $M_{ij}$ (bzw. des verschobenen Werts $M_{ij}/L - c \ln L$) wird hervorragend durch eine Gumbel-Verteilung beschrieben.
  • Im Gegensatz dazu zeigte das Lieb-Liniger-Modell eine Fréchet-Verteilung.
  • Im Gegensatz zur ETH für chaotische Systeme (wo $R_{ij}$ gaußförmig ist) ist die Verteilung hier nicht gaußförmig.
• Nachweis der Nicht-Gauß-Förmigkeit: Die Autoren berechnen das Verhältnis $\Gamma = \overline{|O_{ij}|^2} / |\overline{O_{ij}}|^2$. Für gaußförmige Variablen wäre $\Gamma = \pi/2$. Die numerischen Daten zeigen jedoch $\Gamma \sim 10^5$, was auf eine Divergenz im thermodynamischen Limes hindeutet und die Gumbel-Statistik bestätigt.

D. Ein-Spin vs. Zwei-Spin Operatoren

• Die Ergebnisse für Ein-Spin-Operatoren (mit Strings) und Zwei-Spin-Operatoren (ohne Strings) sind qualitativ identisch. Dies unterstreicht die Robustheit des gefundenen Szenarios.

---

4. Bedeutung und Schlussfolgerungen

1. Verallgemeinerung der ETH für Integrable Systeme: Die Arbeit zeigt, dass das Skalierungsverhalten der nicht-diagonalen Matrixelemente in integrablen Gittermodellen (mit Strings) qualitativ dem von kontinuierlichen Modellen (ohne Strings) entspricht. Das Abklingen ist exponentiell in $L$ (mit Log-Korrektur) für gleiche Makrozustände und wie $e^{-L^2}$ für verschiedene Makrozustände.
2. Neue Statistik: Der wichtigste theoretische Beitrag ist die Identifizierung der Gumbel-Verteilung für die Logarithmen der Matrixelemente in Spin-Ketten. Dies steht im Kontrast zur Fréchet-Verteilung im Lieb-Liniger-Modell und zur Gauß-Verteilung in chaotischen Systemen. Dies deutet darauf hin, dass die Statistik der Matrixelemente empfindlich von den spezifischen Details des Modells (z. B. Vorhandensein von Strings, Gitterstruktur) abhängt.
3. Methodischer Fortschritt: Die Arbeit demonstriert die Machbarkeit der numerischen Berechnung von Matrixelementen für große Systemgrößen ($L \sim 500$) in integrablen Systemen unter Verwendung von ABA-Formeln, unter Umgehung der Probleme durch Singularitäten.
4. Zukunftsaussichten: Die Ergebnisse eröffnen neue Wege, um die Beziehung zwischen der Statistik von Matrixelementen und der Freien Wahrscheinlichkeitstheorie (Free Probability Theory) zu untersuchen, insbesondere im Kontext integrabler Modelle. Zudem könnte dies helfen, die Relaxationsdynamik nach Quanten-Quenches besser zu verstehen.

Zusammenfassend liefert das Paper eine umfassende numerische Charakterisierung der ETH-Verletzung in integrablen Spin-Ketten und etabliert die Gumbel-Statistik als charakteristisches Merkmal für die nicht-diagonalen Matrixelemente in diesem Kontext.