The emergence of (3+1)-dimensional expanding spacetime from complex Langevin simulations of the Lorentzian type IIB matrix model with deformations

1. Problemstellung und Motivation

Das Lorentzsche Typ-IIB-Matrixmodell (auch IKKT-Modell) gilt als vielversprechende nichtstörungstheoretische Formulierung der Superstringtheorie. In diesem Modell entsteht die Raumzeit dynamisch aus den Eigenwertverteilungen der $N \times N$-Boson-Matrizen $A_\mu$ ($\mu = 0, \dots, 9$) im Grenzwert großer $N$.

Das Hauptproblem bei der Untersuchung dieses Modells ist das Vorzeichenproblem (Sign Problem) in der Pfadintegral-Formulierung, verursacht durch den oszillierenden Phasenfaktor $e^{iS}$.

• Äquivalenz zum euklidischen Modell: Ohne Modifikationen ist das Lorentzsche Modell äquivalent zum euklidischen Modell. Dies führt dazu, dass die emergente Raumzeit komplex und euklidisch ist (mit imaginärer Zeit), was nicht der physikalischen Realität entspricht.
• Singularitäten und Drift: Bei der Anwendung der Komplexen Langevin-Methode (CLM), um das Vorzeichenproblem zu umgehen, tritt ein „singuläres Drift-Problem" auf. Dies wird durch die Pfaffian-Determinante (entstanden durch Integration über Fermionen) verursacht, deren Eigenwerte nahe Null liegen und zu divergierenden Drift-Termen führen, was die Konvergenz der Simulation verhindert.
• Frühere Ergebnisse: Bisherige Studien zeigten zwar expandierende Raumzeiten, litten aber oft an singulären Strukturen oder Artefakten durch Lorentz-Boosts, die die Interpretation der Raumzeit-Evolution erschwerten.

2. Methodik

Die Autoren verwenden die Komplexe Langevin-Methode (CLM) mit Matrixgrößen bis zu $N=128$, um das Vorzeichenproblem zu lösen. Um die oben genannten Schwierigkeiten zu überwinden, wurden folgende technische Anpassungen vorgenommen:

• Einführung eines lorentzinvarianten Massenterms: Um die Äquivalenz zum euklidischen Modell zu brechen und eine reelle, lorentzsche Raumzeit zu ermöglichen, wird ein Massenterm $S_\gamma = -\frac{N}{2}\gamma \text{Tr}(A_\mu A^\mu)$ mit $\gamma > 0$ hinzugefügt. Dies ändert die Dynamik so, dass klassische Lösungen eine expandierende Raumzeit darstellen.
• Deformation zur Vermeidung singulärer Drifts: Um das Problem der singulären Drifts durch die Fermionen-Pfaffian zu lösen, wird das Modell deformiert:
  1. Fermionen-Massenterm: Ein Term $S_{mf}$ wird hinzugefügt, der die Eigenwerte der Matrix $M$ (die den Pfaffian bestimmt) vom Ursprung wegschiebt. Dies stabilisiert die CLM, unterdrückt aber gleichzeitig den Beitrag der Fermionen.
  2. Bosonische Anisotropie (SUSY-inspirierte Deformation): Um den Verlust der Fermionenbeiträge zu kompensieren und Supersymmetrie (SUSY) zu approximieren, wird der Massenterm für die Bosonen anisotrop gestaltet. Die 9 Raumdimensionen werden in $\tilde{d}$ und $9-\tilde{d}$ unterteilt, wobei letztere durch einen Parameter $\xi$ unterdrückt werden. Dies ähnelt der „polarisierten" Deformation im euklidischen Fall.
• Entfernung von Lorentz-Boosts: Da das Modell Lorentz-invariant ist, neigen die Konfigurationen dazu, durch zufällige Boosts verzerrt zu werden. Die Autoren entwickeln ein Verfahren, bei dem für jede Konfiguration eine Lorentz-Transformation angewendet wird, um die „Schwerpunktsbewegung" zu minimieren und die Artefakte zu entfernen.
• Dynamische Stabilisierung: Eine Technik, bei der nach jedem Langevin-Schritt die Matrizen leicht in Richtung ihrer hermiteschen Komponente verschoben werden ($A_i \to A_i + \eta A_i^\dagger$), wird angewendet, um die Stabilität zu erhöhen.

3. Wichtige Beiträge

• Durchbruch bei der CLM-Anwendung: Die erfolgreiche Anwendung der CLM auf das Lorentzsche Typ-IIB-Modell mit Fermionenbeiträgen bei großen $N$ ($N=128$), was zuvor aufgrund des singulären Drift-Problems nicht möglich war.
• Neue Deformationsstrategie: Die Kombination aus einem fermionischen Massenterm (zur Stabilisierung) und einer anisotropen bosonischen Deformation (zur Wiederherstellung der SUSY-Effekte) als praktikabler Weg, um die physikalische Dynamik zu untersuchen.
• Entwicklung einer Methode zur Boost-Entfernung: Ein systematischer Ansatz, um Lorentz-Boost-Artefakte aus den numerischen Daten zu filtern, um die echte Expansion der Raumzeit sichtbar zu machen.

4. Ergebnisse

Die numerischen Simulationen führten zu folgenden Schlüsselergebnissen:

• Emergenz einer (3+1)-dimensionalen Raumzeit: Im deformierten Modell mit Fermionenbeiträgen (insbesondere bei kleinem $m_f = 6$) wurde eine Phase beobachtet, in der sich die Raumzeit dynamisch entwickelt.
  • Die Zeit (Eigenwerte von $A_0$) wird reell und glatt.
  • Die Raumdimensionen zeigen eine spontane Symmetriebrechung (SSB) von $SO(9)$ (bzw. $SO(\tilde{d})$) zu $SO(3)$.
  • Drei der neun Raumrichtungen beginnen zu einem bestimmten Zeitpunkt zu expandieren, während die anderen sechs klein bleiben.
• Glätte und Realismus: Im Gegensatz zu früheren Studien, die singuläre Strukturen (beschrieben durch Pauli-Matrizen) zeigten, sind die emergenten Raum und Zeit in diesem Modell glatt und reell. Die Eigenwertverteilungen der Raummatrizen sind dicht verteilt, was auf eine kontinuierliche Raumzeit hindeutet.
• Rolle der Fermionen: Ohne Fermionen (reines Bosonen-Modell) bleibt die $SO(9)$-Symmetrie erhalten, und es tritt keine spontane Dimensionenreduktion auf. Die Fermionenbeiträge (moduliert durch $m_f$) sind entscheidend für die Brechung der Symmetrie hin zu einer (3+1)-dimensionalen Struktur.
• Phasenübergang: Es wurde ein Phasenübergang beobachtet; bei zu großem $m_f$ (zu starke Unterdrückung der Fermionen) kollabiert die (3+1)-dimensionale Struktur wieder zu einem isotropen Zustand.

5. Bedeutung und Ausblick

Diese Arbeit liefert starke numerische Belege dafür, dass das Lorentzsche Typ-IIB-Matrixmodell in der Lage ist, eine physikalisch plausible (3+1)-dimensionale expandierende Raumzeit aus reinen quantenmechanischen Prinzipien zu erzeugen, ohne zusätzliche Parameter für die Kompaktifizierung zu benötigen.

• Herausforderungen: Die Studie nutzt noch immer Näherungen (z. B. die manuelle Entfernung von Boosts, die die Holomorphie der CLM verletzt, und die Verwendung von Massentermen, die nicht exakt der SUSY-Deformation entsprechen).
• Zukünftige Arbeiten: Die Autoren planen, die Lorentz-Symmetrie nicht-störungstheoretisch mittels Faddeev-Popov-Verfahren zu fixieren und die exakte SUSY-deformierte Version des Modells zu simulieren. Zudem wird die Lefschetz-Thimble-Methode als alternative Methode zur Untersuchung des dominanten Sattelpunkts diskutiert, um die Ergebnisse der CLM zu validieren.

Zusammenfassend demonstriert das Paper, dass nicht-störungstheoretische Simulationen der Stringtheorie in der Lage sind, die Entstehung unserer beobachtbaren Raumzeit-Dimensionen zu reproduzieren, was einen bedeutenden Schritt in Richtung einer fundamentalen Theorie des Universums darstellt.