Study of multi-particle states with tensor renormalization group method

Die Studie untersucht die Vielteilchenzustände des (1+1)-dimensionalen Ising-Modells mittels der Tensor-Renormierungsgruppen-Methode, identifiziert die Eigenzustände durch Symmetrien und Interpolationsoperatoren und berechnet die Zwei-Teilchen-Streuphasenverschiebung unter Verwendung von Lüschers Formel sowie der Wellenfunktionsmethode.

Das Wesentliche

Titel: Wie man unsichtbare Partikel-Partys mit einem digitalen Lupe-Verfahren findet

Stellen Sie sich vor, Sie sind ein Detektiv in einer winzigen, zweidimensionalen Welt (wie auf einem Schachbrett), in der kleine Magnete (Spin-Teilchen) hin und her wackeln. Ihr Ziel ist es herauszufinden, was passiert, wenn diese Magnete nicht nur einzeln, sondern in Gruppen zusammenarbeiten – also wenn sie „Partys" feiern. Die Wissenschaftler in diesem Papier haben eine neue Methode entwickelt, um diese Partys zu zählen und zu verstehen, ohne dabei den Überblick zu verlieren.

Hier ist die Geschichte, wie sie das gemacht haben, einfach erklärt:

1. Das Problem: Der Lärm im Keller

Normalerweise versuchen Physiker, solche Teilchen-Partys zu beobachten, indem sie riesige Simulationen auf Computern laufen lassen. Das ist wie der Versuch, ein leises Flüstern in einem lauten Stadion zu hören. Je höher die Energie (je mehr Gäste auf der Party), desto lauter wird der Lärm (statistisches Rauschen), und desto schwieriger ist es, die einzelnen Stimmen zu unterscheiden.

2. Die Lösung: Ein smarter „Kochtopf" (Tensor-Netzwerk)

Die Forscher haben eine andere Methode gewählt, die sie Tensor-Renormierungsgruppe nennen. Stellen Sie sich das vor wie einen riesigen, komplexen Kochtopf, in dem alle möglichen Zustände der Magnete enthalten sind.

• Das alte Rezept: Bisher haben sie den Topf immer in quadratischen Formen gekocht. Das funktionierte gut für die ersten, einfachen Gerichte (die Grundzustände), aber sobald sie komplexere Gerichte (höhere angeregte Zustände) zubereiten wollten, wurde das Essen verbrannt (die Fehler wurden zu groß).
• Das neue Rezept: In dieser Arbeit haben sie den Topf anders geschnitten. Anstatt ihn quadratisch zu kochen, haben sie ihn langgestreckt gemacht (wie einen langen Streifen). Durch dieses „Rechtwinklig-Schneiden" konnten sie die komplexeren Gerichte viel klarer schmecken, ohne dass sie verbrannten.

3. Die Detektivarbeit: Wer ist auf der Party?

Sobald sie die Energie der Zustände berechnet hatten, mussten sie herausfinden: Wie viele Gäste sind auf dieser Party?

• Der Zähler: Sie benutzten eine Art „Zauberspiegel" (einen mathematischen Operator), der ihnen sagte, ob eine Gruppe von Teilchen eine gerade oder ungerade Anzahl hat (eine Eigenschaft namens Quantenzahl).
• Die Größe der Party: Um zu wissen, ob es eine Ein-Personen-Party, eine Zwei-Personen-Party oder eine Drei-Personen-Party ist, haben sie die Größe des Raumes (des Schachbretts) verändert.
  • Die Analogie: Wenn Sie eine einzelne Person in einem Raum haben, ändert sich ihre Energie kaum, wenn Sie den Raum vergrößern. Aber wenn zwei oder drei Personen miteinander tanzen (wechselwirken), ändert sich ihre Energie deutlich, je größer der Tanzboden wird. Indem sie den Raum immer größer machten, konnten sie genau zählen: „Aha, diese Energiekurve passt zu einer Zwei-Personen-Party!"

4. Der Tanzschritt: Wellenfunktionen und Streuung

Nicht nur die Anzahl der Gäste war wichtig, sondern auch, wie sie sich bewegen.

• Die Wellenfunktion: Das ist wie eine Karte, die zeigt, wo die Gäste sich wahrscheinlich aufhalten. Die Forscher haben diese Karte für die Zwei-Personen-Partys erstellt. Sie sahen, wie die beiden Teilchen sich anziehen oder abstoßen, als würden sie auf einem Tanzboden tanzen.
• Der Streuungs-Winkel: Am Ende wollten sie wissen: Wenn zwei Teilchen aufeinanderprallen, wie stark werden sie abgelenkt? Das nennen sie den „Streuungs-Winkel" (Phasenverschiebung).
  • Sie haben dies auf zwei Arten berechnet: einmal über die Energie der Party (Lüscher-Formel) und einmal über die Tanzkarte (Wellenfunktion).
  • Das Ergebnis: Beide Methoden lieferten das exakt gleiche Ergebnis, das auch mit der perfekten, theoretischen Vorhersage übereinstimmte.

Zusammenfassung

Die Forscher haben bewiesen, dass man durch eine clevere Umstellung der Rechenmethode (den „langgestreckten Topf") nicht nur die einfachen, sondern auch die sehr komplexen Teilchen-Partys in einer einfachen Welt (dem Ising-Modell) genau beobachten kann. Sie haben die Anzahl der Teilchen gezählt, ihre Tanzbewegungen kartiert und bewiesen, dass ihre neue Methode so präzise ist wie ein Maßband.

Dies ist ein großer Schritt, um in Zukunft auch noch komplexere Universen und Teilchenwechselwirkungen zu verstehen, ohne vom Lärm der Simulation überwältigt zu werden.

Technische Zusammenfassung

Titel: Untersuchung von Vielteilchenzuständen mit der Tensor-Renormierungsgruppen-Methode

Autoren: Fathiyya Izzatun Az-zahra, Shinji Takeda, Takeshi Yamazaki
Kontext: 42. Internationales Symposium für Gitter-Feldtheorie (Lattice Field Theory), Mumbai, Indien, November 2025.

---

1. Problemstellung und Motivation

Die Untersuchung von Vielteilchenwechselwirkungen in der Gitterfeldtheorie ist eine zentrale Herausforderung. Während die Monte-Carlo-Simulationen ein leistungsfähiges Werkzeug darstellen, leiden sie unter spezifischen technischen Limitationen:

• Statistisches Rauschen: Die Extraktion höherer angeregter Zustände ist durch großes statistisches Rauschen belastet.
• Zeitliche Ausdehnung: Um Kontaminationen durch angeregte Zustände zu unterdrücken, sind Gitter mit sehr großer zeitlicher Ausdehnung ($L_t$) erforderlich.

Alternativansätze basierend auf Tensor-Netzwerken (speziell für den Lagrange-Formalismus) wurden eingeführt. Diese sind deterministisch und benötigen nur eine kleine zeitliche Ausdehnung, um das Spektrum zu extrahieren. Allerdings zeigen bisherige Methoden, die quadratische Tensor-Netzwerke ($L_t = L_s$) verwenden, drastisch anwachsende Fehler für höher angeregte Zustände. Dies liegt an der starken Entartung der Eigenwerte des Transferoperators bei kleinen $L_t$ und der Notwendigkeit vieler Koarse-Graining-Iterationen bei großen $L_t$.

Ziel der Arbeit: Entwicklung und Anwendung einer modifizierten Koarse-Graining-Strategie, um auch höher angeregte Zustände (die Vielteilchenzuständen entsprechen) mit hoher Präzision aus dem Tensor-Netzwerk zu extrahieren. Als Testfall dient das (1+1)-dimensionale Ising-Modell.

---

2. Methodik

A. Spektralanalyse mit Transfermatrix

Die Autoren nutzen die Transfermatrix-Methode ($T$), wobei die Partitionfunktion als Spur über $T^{L_t}$ geschrieben wird. Die Energieeigenwerte $E_a$ ergeben sich aus den Eigenwerten $\lambda_a$ der Transfermatrix via $\lambda_a = e^{-E_a}$.

• Identifikation von Quantenzahlen: Da die Quantenzahlen der Eigenzustände nicht a priori bekannt sind, werden Matrixelemente eines geeigneten interpolierenden Operators $\hat{O}_q$ berechnet. Unter Ausnutzung der $Z_2$-Symmetrie des Ising-Modells werden Zustände nach ihrer Parität ($q = \pm 1$) klassifiziert.
• Impuritäts-Tensor-Netzwerk: Um Matrixelemente wie $\langle \Omega | \hat{O}_q | a \rangle$ zu berechnen, wird ein "Impurity"-Tensor-Netzwerk eingeführt, bei dem der Operator an einer bestimmten Stelle in das Netzwerk eingefügt wird.

B. Modifizierte Tensor-Renormierungsgruppe (HOTRG)

Der Kern der methodischen Innovation liegt in der Anpassung des Koarse-Graining-Prozesses für Tensor-Netzwerke mit $1 < L_t < L_s$:

1. Zweidimensionales Koarse-Graining: Zuerst werden $p$ Iterationen des Higher-Order Tensor Renormalization Group (HOTRG) Algorithmus in beiden Richtungen (Raum und Zeit) durchgeführt.
2. Eindimensionales Koarse-Graining: Anschließend werden $(n - p - 1)$ Iterationen nur in der Raumrichtung durchgeführt.
3. Exakte Kontraktion: Die verbleibenden Tensoren werden exakt kontrahiert, um eine effektive Transfermatrix zu erhalten, deren Eigenwerte das Energiespektrum approximieren.

Diese Strategie ($L_t = 8, L_s = 64$ im Beispiel) minimiert die Entartung der Singulärwerte und reduziert die Anzahl der numerisch fehleranfälligen Iterationen im Vergleich zu quadratischen Gittern.

C. Wellenfunktions- und Streuphase-Bestimmung

• Wellenfunktion: Die Wellenfunktion von Zwei-Teilchen-Zuständen wird über Matrixelemente eines Operators $\hat{O}_2(x)$ extrahiert.
• Streuphase: Die Streuphase $\delta(k)$ wird auf zwei Wegen berechnet und verglichen:
  1. Über die endliche Volumen-Energie mittels Lüscher-Formel.
  2. Durch Anpassung (Fitting) der extrahierten Wellenfunktion an eine freie Wellenfunktion außerhalb des effektiven Wechselwirkungsbereichs $R$.

---

3. Wichtige Beiträge und Ergebnisse

A. Verbesserte Genauigkeit des Energiespektrums

Die Simulationen am Ising-Modell bei $T=2.44$ (symmetrische Phase) zeigen, dass die Wahl von $L_t = 8$ (bei $L_s = 64$) im Vergleich zu $L_t = 2, 4, 16, 64$ die geringsten relativen Fehler im Energiespektrum liefert.

• Dies ermöglicht die zuverlässige Berechnung von Energien bis hin zu sehr hoch angeregten Zuständen (bis $a=42$), was mit früheren Methoden nicht möglich war.
• Die Quantenzahlen ($q=\pm 1$) wurden erfolgreich durch Analyse der Matrixelemente identifiziert und stimmen exakt mit der analytischen Vorhersage überein.

B. Identifikation von Vielteilchenzuständen

Durch die Analyse der Systemgrößenabhängigkeit ($L_s$) der Energieniveaus konnten die Autoren die Teilchenzahl der Zustände bestimmen:

• $q = -1$ Sektor: Der tiefste Zustand entspricht einem Einteilchenzustand (Masse $m$). Höhere Zustände nähern sich $3m$ an und wurden als Dreiteilchenzustände identifiziert.
• $q = +1$ Sektor: Die drei tiefsten Niveaus nähern sich $2m$ an und wurden als Zweiteilchenzustände (Grundzustand und zwei angeregte Zustände) identifiziert. Ein weiterer Zustand (violette Dreiecke in Fig. 4b) konnte aufgrund seiner Wellenfunktionsform und der nur bei kleinen Systemgrößen verfügbaren Daten nicht eindeutig klassifiziert werden.

C. Wellenfunktionen und Streuphase

• Die Wellenfunktionen für die Zweiteilchenzustände wurden erfolgreich extrahiert. Die Analyse zeigt, dass der effektive Wechselwirkungsbereich bei $T=2.44$ etwa $R \approx 40$ beträgt.
• Die Streuphase $\delta(k)$ wurde sowohl über die Lüscher-Formel als auch über das Wellenfunktions-Fitting berechnet.
• Ergebnis: Beide Methoden liefern konsistente Ergebnisse, die exzellent mit der analytischen Vorhersage $\delta(k) = -\pi/2$ übereinstimmen. Dies validiert die neue Wellenfunktions-Methode als alternative und zuverlässige Technik zur Bestimmung von Streuphasen.

---

4. Bedeutung und Ausblick

Diese Arbeit demonstriert einen signifikanten Fortschritt in der Anwendung von Tensor-Netzwerk-Methoden auf Vielteilchenprobleme in der Quantenfeldtheorie:

1. Überwindung von Limitierungen: Die modifizierte Koarse-Graining-Strategie ($L_t \ll L_s$) löst das Problem der stark wachsenden Fehler bei hoch angeregten Zuständen, das bei früheren quadratischen Ansätzen bestand.
2. Deterministische Vielteilchenphysik: Die Methode ermöglicht die präzise Untersuchung von Vielteilchenzuständen (bis zu drei Teilchen) ohne statistisches Rauschen, was für Monte-Carlo-Simulationen oft prohibitiv ist.
3. Validierung der Wellenfunktions-Methode: Die erfolgreiche Extraktion der Streuphase aus der Wellenfunktion bestätigt, dass Tensor-Netzwerke nicht nur Energien, sondern auch räumliche Korrelationsfunktionen und Streudaten hochpräzise liefern können.

Zukünftige Arbeiten: Die Autoren planen, die identifizierten Dreiteilchenzustände detaillierter zu analysieren und das Schema auf andere Quantenfeldtheorien anzuwenden.

Fazit: Das vorgestellte Spektrum-Verfahren auf Basis des Tensor-Netzwerks ist ein leistungsfähiges, deterministisches Werkzeug zur Untersuchung der Vielteilchen-Dynamik in niedrigdimensionalen Modellen und bietet eine vielversprechende Alternative oder Ergänzung zu herkömmlichen Gitter-QCD-Methoden.