An explicit formula for perturbation theory at any order with infinitely many perturbations

Die Autoren stellen eine systematische Formel vor, die auf ganzzahligen Partitionen basiert und eine explizite Störungstheorie beliebiger Ordnung für unendlich viele Störungen mittels einer einzigen Matrixgleichung ermöglicht, die sowohl Eigenwert- als auch Eigenvektorkorrekturen umfasst und die Berechnungen gegenüber traditionellen Methoden vereinfacht.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, das Wetter vorherzusagen. Sie haben eine perfekte, einfache Formel für einen sonnigen Tag ohne Wind (das ist Ihr ungestörtes System). Aber die Realität ist chaotisch: Es gibt plötzlich Windböen, Regen, Hagel und sogar einen kleinen Tornado. Jeder dieser Störfaktoren verändert das Wetter ein wenig.

In der Physik nennt man das Störungstheorie. Man versucht, die komplexe Realität durch das Hinzufügen von immer kleineren „Korrektur-Schichten" zu beschreiben. Das Problem ist: Je genauer man sein will (also je mehr Störfaktoren man berücksichtigt), desto unübersichtlicher wird die Mathematik. Die Formeln werden so lang und kompliziert, dass sie kaum noch zu lesen sind – wie ein Knoten, den niemand mehr entwirren kann.

Dieses Papier von Joseph Jones und M. W. Long bietet nun eine neue, clevere Anleitung, um diesen Knoten zu entwirren. Hier ist die Erklärung in einfachen Worten:

1. Das Problem: Der endlose Knoten

Normalerweise betrachtet Physiker nur einen oder zwei Störfaktoren (z. B. nur den Wind). Aber in der modernen Physik (besonders bei Quantenmechanik) gibt es oft unendlich viele kleine Störungen gleichzeitig. Wenn man versucht, die Auswirkungen dieser unendlichen Menge bis zur 10. oder 20. Genauigkeit zu berechnen, explodiert die Komplexität. Die alten Methoden sind wie ein Werkzeugkasten, der nur einen kleinen Schraubenzieher hat – für einen riesigen Schraubenberg ungeeignet.

2. Die Lösung: Der „Zahlensalat" (Ganzzahl-Partitionen)

Die Autoren haben eine brillante Idee: Sie nutzen die Mathematik der Ganzzahl-Partitionen.
Stellen Sie sich vor, Sie haben eine Zahl, sagen wir 4. Wie können Sie diese Zahl als Summe kleinerer Zahlen schreiben, wobei die Reihenfolge wichtig ist?

• 4
• 3 + 1
• 1 + 3
• 2 + 2
• 2 + 1 + 1
• 1 + 2 + 1
• 1 + 1 + 2
• 1 + 1 + 1 + 1

Das sind alle Möglichkeiten, die Zahl 4 zu „zerlegen". Die Autoren sagen: „Jede dieser Zerlegungen entspricht genau einem Schritt in unserer physikalischen Berechnung."

Statt sich in endlosen Formeln zu verlieren, sagen sie einfach: „Nimm die Zahl der gewünschten Genauigkeit (z. B. 4), zerlege sie in alle möglichen Reihenfolgen, und jede Zerlegung ist eine Bauregel für einen Term in der Formel."

3. Der magische Bauplan (Die Matrix)

Die Autoren haben eine Art Master-Formel (eine Matrix-Gleichung) entwickelt.

• Die alte Methode: Man musste für jede neue Stufe (z. B. von Stufe 3 auf 4) die ganze Formel neu herleiten und anpassen. Das war wie jedes Mal ein neues Haus zu bauen, statt nur einen Stockwerk hinzuzufügen.
• Die neue Methode: Sie haben einen einzigen, universellen Bauplan. Wenn Sie wissen wollen, wie sich das System bei der 100. Störung verhält, müssen Sie nicht neu erfinden, wie man rechnet. Sie nutzen einfach die „Zerlegungs-Liste" für die Zahl 100 und setzen die Bausteine (die Störungen) in die Lücken ein.

Es ist, als hätten sie einen 3D-Drucker für physikalische Formeln gebaut. Sie geben nur die Zahl ein (die gewünschte Genauigkeit), und der Drucker spuckt die perfekte, korrekte Formel aus – egal, ob es nur einen Störfaktor gibt oder unendlich viele.

4. Warum ist das revolutionär?

• Unendliche Störungen: Früher war es fast unmöglich, Formeln für unendlich viele Störungen aufzustellen. Jetzt ist es ein Standard-Verfahren.
• Ein Werkzeug für alles: Die gleiche Formel berechnet sowohl die Energieänderung (wie stark sich das Wetter ändert) als auch die Zustandsänderung (wie sich die Wolken bewegen). Früher brauchte man dafür zwei verschiedene, komplizierte Wege.
• Einfachheit: Was früher 10 Seiten an Algebra brauchte, passt jetzt auf eine Seite, wenn man die „Zerlegungs-Symbole" benutzt.

Zusammenfassung mit einer Analogie

Stellen Sie sich vor, Sie wollen ein riesiges Puzzle zusammenbauen.

• Früher: Sie mussten für jedes Puzzleteil eine eigene Anleitung schreiben. Wenn das Puzzle 1000 Teile hatte, waren Sie 1000 Seiten lang in Anleitungstexten verloren.
• Jetzt: Die Autoren haben ein System entwickelt. Sie sagen: „Nimm die Zahl 1000. Zerlege sie in alle möglichen Summen. Jede Summe sagt dir, welches Puzzleteil wohin gehört."
  • Die „Zerlegung" ist der Bauplan.
  • Die „Störungen" sind die Puzzleteile.
  • Das Ergebnis ist das fertige Bild (die physikalische Lösung).

Fazit:
Dieses Papier ist wie ein Universal-Schlüssel für die Physik. Es verwandelt ein chaotisches, unübersichtliches mathematisches Problem in ein systematisches, fast spielerisches Sortier-Verfahren. Es ermöglicht Wissenschaftlern, Berechnungen durchzuführen, die bisher als zu schwierig galten, und macht die Mathematik hinter den komplexesten Phänomenen des Universums endlich verständlich und handhabbar.

Technische Zusammenfassung

Titel: Eine explizite Formel für die Störungstheorie beliebiger Ordnung mit unendlich vielen Störungen

Autoren: Joseph M. Jones und M. W. Long (Universität Birmingham)

---

1. Problemstellung

Die Störungstheorie ist ein Grundpfeiler der theoretischen Physik zur Berechnung näherungsweiser Lösungen, wenn exakte Ergebnisse nicht verfügbar sind. Trotz ihrer weit verbreiteten Anwendung beschränken sich praktische Berechnungen meist auf niedrige Ordnungen (oft bis zur zweiten Ordnung).

• Herausforderung: Explizite Formeln für Korrekturen beliebiger Ordnung sind aufgrund der algebraischen Komplexität höherer Terme selten und schwer zu handhaben.
• Lücke: Bisherige Ansätze (wie Rayleigh-Schrödinger oder Brillouin-Wigner) behandeln typischerweise nur einen einzelnen Störoperator oder erfordern iterative, selbstkonsistente Gleichungen. Es fehlte eine systematische, geschlossene Formel, die sowohl für eine endliche als auch für eine unendliche Menge von Stör-Hamiltonianern gilt und gleichzeitig Eigenwerte und Eigenvektoren aus einer einzigen Matrixgleichung ableitet.
• Motivation: Ein physikalischer Auslöser war die Notwendigkeit, Störungstheorie auf die unendliche Reihenentwicklung der Baker-Campbell-Hausdorff-Formel anzuwenden, bei der $A$ als exakt lösbarer Referenzpunkt und $B$ als unendliche Reihe von Störungen fungiert.

2. Methodik

Die Autoren entwickeln einen neuen formalen Rahmen, der auf geordneten ganzzahligen Partitionen (im Englischen oft als "Compositions" bezeichnet) basiert.

• Generierende Funktionen:
  Es werden generierende Funktionen für den Hamiltonian $H(z)$, die Energie $E(z)$ und den Eigenzustand $\psi(z)$ definiert, wobei $z$ die Ordnungsparameter steuert.
  $$H(z) = H_0 + \sum_{n=1}^{\infty} H^{(n)} z^n$$
• Definition der zentralen Größe $\Delta H(n)$:
  Um die Struktur zu vereinfachen, führen die Autoren eine modifizierte Störgröße ein:
  $$\Delta H(n) \equiv H^{(n)} - E^{(n)}Q$$
  wobei $Q = 1 - P$ der Projektor auf den Unterraum orthogonal zum ungestörten Grundzustand $|0\rangle$ ist. Diese Definition fasst Energie- und Zustandskorrekturen in einem einzigen Operator zusammen und eliminiert "disconnected" (nicht-extensive) Beiträge.
• Die Resolvente $\Gamma$:
  Es wird die Diagonal-Matrix $\Gamma \equiv Q(\epsilon_0 - H_0)^{-1}Q$ eingeführt, die als Nennermatrix in der Störungstheorie fungiert.
• Rekursive Matrixgleichung:
  Die Kernidee ist die Einführung einer Matrix $M(z)$, die durch eine rekursive Beziehung definiert ist:
  $$M(z) = \Delta H(z) + \Delta H(z) \Gamma M(z)$$
  Dies lässt sich explizit als geometrische Reihe auflösen:
  $$M(z) = (1 - \Delta H(z)\Gamma)^{-1} \Delta H(z)$$

3. Hauptergebnisse und Formeln

Das Paper liefert eine explizite, geschlossene Formel für die Korrekturen $N$-ter Ordnung, ausgedrückt als Summe über alle geordneten Partitionen der Zahl $N$.

• Allgemeine Lösung für $N$-te Ordnung:
  Die Korrekturen für den Eigenwert $E^{(N)}$ und den Eigenvektor $|\psi^{(N)}\rangle$ lauten:
  $$E^{(N)} = \langle 0 | M^{(N)} | 0 \rangle$$
  $$|\psi^{(N)}\rangle = \Gamma M^{(N)} | 0 \rangle$$
  wobei $M^{(N)}$ gegeben ist durch:
  $$M^{(N)} = \sum_{p=1}^{N} \sum_{n_1, \dots, n_p} \delta_{n_1 + \dots + n_p, N} \, \Delta H^{(n_1)} \Gamma \Delta H^{(n_2)} \dots \Gamma \Delta H^{(n_p)}$$
  Hierdurch werden alle möglichen geordneten Zerlegungen (Partitionen) von $N$ in positive ganze Zahlen $n_i$ abgedeckt.

• Notation und Vereinfachung:
  Um die Formeln übersichtlich zu halten, führen die Autoren eine Kurzschreibweise ein, z.B. $(n_1 n_2 \dots n_p)$ für den Ausdruck $\Delta H^{(n_1)} \Gamma \dots \Gamma \Delta H^{(n_p)}$.

  • Beispiel für $N=4$: $M^{(4)} = (4) + (13) + (22) + (31) + (112) + (121) + (211) + (1111)$.
  • Die Anzahl der Terme entspricht $2^{N-1}$ (die Anzahl der geordneten Partitionen von $N$).

• Einzelne Störung (Grenzfall):
  Im Fall einer einzigen Störung ($H^{(n)} = 0$ für $n > 1$) reduziert sich die Formel auf den bekannten Rayleigh-Schrödinger-Fall. Dabei gilt $\Delta H^{(n)} \to -E^{(n)}Q$ für $n \ge 2$. Die Formel reproduziert automatisch die bekannten Ergebnisse bis zur zehnten Ordnung, jedoch in einer kompakteren, algorithmisch generierbaren Form.

4. Technische Details und Implementierung

• Algorithmische Handhabbarkeit: Die Autoren stellen Code (in Mathematica) bereit, der die Terme für beliebige Ordnungen automatisch generiert.
• Speicherbedarf: Für viele Körperprobleme (z.B. $2^N \times 2^N$ Matrizen) ist es effizienter, nur den Vektor $M(z)|0\rangle$ zu berechnen, anstatt die volle Matrix $M(z)$ zu speichern.
• Rekursion: Die Berechnung von $E^{(N)}$ erfordert nur Korrekturen bis zur Ordnung $N-2$ (da Terme der Ordnung $N-1$ an den Rändern durch Projektion auf $|0\rangle$ verschwinden), während $|\psi^{(N)}\rangle$ die Ordnung $N-1$ benötigt. Dies vermeidet selbstkonsistente Gleichungen.

5. Bedeutung und Vergleich mit bestehender Literatur

• Vergleich mit Kato, Löwdin, Soliverez und BP:
  • Während Kato die Verbindung zu Partitionen für Konvergenzanalysen nutzte, liefert dieses Paper eine direkte, explizite Formel für die Korrekturen.
  • Im Gegensatz zu Löwdins Ansatz, der oft auf Brillouin-Wigner (selbstkonsistent) basiert, liefert diese Methode eine direkte rekursive Lösung ohne Selbstkonsistenzprobleme.
  • Die Formeln von Soliverez und Bracci/Picasso (BP) sind äquivalent, aber algebraisch umständlicher und weniger kompakt als die hier vorgestellte Summe über Partitionen.
• Innovation:
  1. Unendliche Störungen: Der Ansatz ist der erste, der eine unendliche Menge von Stör-Hamiltonianern systematisch behandelt.
  2. Einheitliche Formel: Eigenwerte und Eigenvektoren werden aus derselben Matrixgleichung $M^{(N)}$ abgeleitet.
  3. Vereinfachung: Die Verwendung von $\Delta H(n)$ und der Partitionssumme eliminiert die Notwendigkeit, Terme manuell zu gruppieren oder zu kürzen.

6. Fazit und Ausblick

Die Autoren haben eine systematische, auf ganzzahligen Partitionen basierende Formel entwickelt, die die Störungstheorie beliebiger Ordnung für eine beliebige Anzahl von Stör-Hamiltonianern explizit und effizient berechnet.

• Anwendungsgebiete: Die Methode ist besonders relevant für Probleme, die durch die Baker-Campbell-Hausdorff-Formel motiviert sind (z.B. in der Quantenfeldtheorie oder statistischen Mechanik).
• Zukünftige Arbeiten: Die Autoren planen, diese Methode auf spezifische physikalische Modelle anzuwenden, bei denen die Störungen durch die BCH-Reihe definiert sind, was zu dramatischen Vereinfachungen führen könnte, etwa bei der Untersuchung von Phasenübergängen mittels Transfermatrix-Methoden.

Zusammenfassend bietet das Paper ein mächtiges Werkzeug, das die algebraische Komplexität höherer Ordnungen der Störungstheorie durch eine elegante kombinatorische Struktur beherrschbar macht.