Stationary perturbation theory without sums over intermediate states: Supersymmetric Expansion Algorithm

Technisches Resümee: Stationäre Störungstheorie ohne Summen über Zwischenzustände: Supersymmetrischer Expansionsalgorithmus

Problemstellung
Die Arbeit adressiert eine grundlegende Einschränkung der konventionellen Rayleigh-Schrödinger-Störungstheorie (RSPT). Während die RSPT einen systematischen Rahmen zur Berechnung von Korrekturen für Energieeigenwerte und Eigenzustände bietet, stützt sie sich auf Summen über Zwischenzustände. Die Korrektur erster Ordnung für Eigenzustände erfordert bereits eine „umständliche Summe“, und Korrekturen höherer Ordnung erfordern in jeder Ordnung zunehmend komplexe Summen über Zwischenzustände. Dieser Formalismus wird rechentechnisch mühsam und ist oft jenseits der niedrigsten Ordnungen unpraktikabel. Vorherige Alternativen, wie etwa Summenregeln oder die logarithmische Formulierung der Schrödinger-Gleichung, die von Aharonov und Au vorgeschlagen wurde, konnten die Berechnungen des Grundzustands für eindimensionale Systeme erfolgreich auf Quadraturformen (Integrale) reduzieren, bleiben jedoch aufgrund der Anwesenheit von Knoten für angeregte Zustände weniger handlich.

Methodik: Der Supersymmetrische Expansionsalgorithmus (SEA)
Die Autoren schlagen den Supersymmetrischen Expansionsalgorithmus (SEA) vor, um perturbative Korrekturen beliebiger Ordnung für alle Eigenzustände (Grund- und angeregte Zustände) zu erhalten, ohne Summen über Zwischenzustände zu verwenden. Die Methode nutzt die verborgene Supersymmetrie, die in nicht-relativistischen quantenmechanischen Problemen inhärent ist.

Das Kernverfahren umfasst die folgenden Schritte:

1. Logarithmische Formulierung: Die Schrödinger-Gleichung wird in eine nicht-lineare logarithmische Form gebracht, die eine Superpotenzial $W(x, \lambda) = -\frac{d}{dx} \ln u(x, \lambda)$ beinhaltet, wobei $u$ die Wellenfunktion ist.
2. Perturbative Expansion: Das Superpotenzial, das Potenzial und die Energie werden in Potenzen eines kleinen Parameters $\lambda$ expandiert. Dies transformiert die nicht-lineare Differentialgleichung in eine hierarchische Kaskade gekoppelter erster Ordnung linearer Differentialgleichungen für die Expansionskoeffizienten.
3. Knotenlose „Randzustände“ (Edge States): Der Algorithmus löst zuerst den knotenlosen Eigenzustand des Hamilton-Operators eines gegebenen supersymmetrischen Niveaus. Für eindimensionale Systeme ist dies der Grundzustand; für höherdimensionale Systeme, die auf effektive eindimensionale Potenziale reduzierbar sind (wie z. B. Zentralpotentiale), sind dies die Zustände mit dem höchsten Drehimpuls $l$ für eine gegebene Hauptquantenzahl $n$.
4. Rekursive Lösung mittels Supersymmetrie:
   • Die knotenlose Lösung des $r$-ten supersymmetrischen Partner-Hamilton-Operators $H_r$ wird mittels der logarithmischen Expansion gelöst.
   • Der führende Ordnungsterm für die knotenlose Lösung von $H_r$ wird aus den bekannten ungestörten Lösungen des ursprünglichen Hamilton-Operators $H_0$ abgeleitet, indem die inversen supersymmetrischen Hebeoperatoren ($a^\dagger$) rekursiv angewendet werden.
   • Sobald die knotenlose Lösung für $H_r$ etabliert ist, werden die perturbativen Korrekturen für dessen Energie und Superpotenzial unter Verwendung von Integralen berechnet, die durch die ungestörte Wahrscheinlichkeitsdichte dieses spezifischen knotenlosen Zustands gewichtet sind.
   • Die Eigenzustände des ursprünglichen Hamilton-Operators $H_0$ mit $r$ Knoten werden dann durch Anwendung der Kette von Hebeoperatoren ($a^\dagger_0, a^\dagger_1, \dots$) auf die knotenlose Lösung von $H_r$ rekonstruiert.

Zentrale Beiträge und Ergebnisse

• Eliminierung von Zwischenzustandssummen: Der primäre Beitrag ist die Reduktion perturbativer Korrekturen auf „Quadraturformen“. Energiekorrekturen und Wellenfunktionskorrekturen werden als Integrale ausgedrückt, die durch die Wahrscheinlichkeitsdichten der relevanten Randzustände gewichtet sind, wodurch die in der RSPT erforderliche Summation über Zwischenzustände vollständig vermieden wird.
• Beliebige Ordnung und allgemeine Anwendbarkeit: Das Formalismus ist für beliebige Ordnungen der Störungstheorie anwendbar und erstreckt sich auf alle Eigenzustände, einschließlich angeregter Zustände mit Knoten. Es funktioniert sowohl für polynomielle als auch für nicht-polynomielle Potenziale.
• Analytische Lösungen für spezifische Potenziale: Die Autoren demonstrieren die Effizienz der Methode durch die Rekonstruktion bekannter analytischer Lösungen für die Yukawa-, Hulthén-, Anharmonischen- und Cornell-Potenziale, bei denen die Koeffizienten des Superpotenzials Polynome sind.
• Nicht-polynomiales Beispiel: Das Paper wendet die Methode auf ein Teilchen in einer Box mit einer harmonischen Oszillator-Störung an. Dieser Fall beinhaltet nicht-polynomielle Superpotenziale. Die Autoren leiten erfolgreich die Energiekorrekturen bis zur Ordnung $\lambda^3$ sowohl für den Grundzustand als auch für den ersten angeregten Zustand ab und liefern explizite analytische Ausdrücke für die Energieeigenwerte.

Bedeutung und Behauptungen
Das Paper behauptet, dass der SEA eine „vollständige analytische Lösung“ zu Perturbationsproblemen in Form einer Potenzreihe in dem Expansionsparameter liefert. Durch die Nutzung der verborgenen Supersymmetrie quantenmechanischer Systeme bietet die Methode eine transparente und systematische Alternative zur RSPT. Die Autoren betonen, dass während frühere Methoden (wie die von Aharonov und Au) durch die Handhabung angeregter Zustände begrenzt waren, der SEA die Lösung angeregter Zustände natürlich mit der supersymmetrischen Hierarchie verbindet, was die Berechnung von Korrekturen für alle Zustände einfach macht. Die Arbeit behauptet, dass dieser Ansatz besonders leistungsfähig für Systeme ist, bei denen die Koeffizienten der Potenzialexpansion keine Polynome sind, da er nicht auf die Reduktion des Problems auf algebraische Rekurrenzrelationen angewiesen ist, sondern die Differentialhierarchie direkt löst.