Excited States Of Positronium From The Two Body Dirac Equations Of Constraint

Diese Arbeit berechnet die Bindungsenergien angeregter Zustände von Positronium mithilfe der Zwei-Körper-Dirac-Gleichungen mit Nebenbedingungen, vergleicht nichtstörungstheoretische mit störungstheoretischen Ergebnissen und korrigiert dabei Fehler in der bestehenden Literatur.

Das Wesentliche

Das Tanzpaar aus Licht und Materie: Eine Reise durch die Positronium-Welt

Stellen Sie sich vor, Sie beobachten ein sehr spezielles Tanzpaar im Universum. Ein Partner ist ein Elektron (ein winziges Teilchen mit negativem Ladung), und der andere ist ein Positron (sein exaktes Spiegelbild, aber mit positiver Ladung). Wenn sie sich treffen, tanzen sie kurzzeitig zusammen, bevor sie sich gegenseitig auslöschen (annihilieren). Dieses kurzlebige Tanzpaar nennt man Positronium.

Physiker wollen genau wissen: Wie fest halten sie sich aneinander? Wie viel Energie steckt in diesem Tanz? Das ist die Frage, die Robert W. Johnson in diesem Papier beantwortet.

1. Das alte Problem: Der ungenaue Tanzlehrer

Bisher haben die meisten Physiker versucht, diesen Tanz mit einer Methode zu beschreiben, die wie ein Schritt-für-Schritt-Anleitung funktioniert (die sogenannte "Störungstheorie"). Das funktioniert gut, wenn die Tänzer langsam und vorsichtig agieren. Aber im Quantenuniversum tanzen sie manchmal sehr wild und schnell. Da wird die Schritt-für-Schritt-Anleitung ungenau, wie eine Landkarte, die nur für flaches Gelände gemacht ist, aber plötzlich Berge zeigt.

Johnson sagt: "Lass uns nicht nur die Schritte zählen, sondern das ganze Tanzpaar gleichzeitig betrachten." Er nutzt dafür eine alte, aber mächtige Methode, die auf den berühmten Physiker Paul Dirac zurückgeht: Die Zwei-Körper-Dirac-Gleichungen.

2. Die neue Brille: Die "Zwangsgleichungen"

Stellen Sie sich vor, die beiden Tänzer sind an einer unsichtbaren, elastischen Schnur gebunden. Die "Zwangsgleichungen" (Constraint Equations) sind wie ein mathematisches Regelwerk, das beschreibt, wie diese Schnur straff wird und wie sich die Tänzer gegenseitig beeinflussen, ohne dass man sie in kleine Stücke zerlegen muss.

Johnson hat diese Gleichungen in einen Computercode übersetzt (für das Programm GNU Octave), um die Energie des Tanzes direkt zu berechnen, ohne die vereinfachten Annahmen der alten Methode. Es ist, als würde man einen Film in Echtzeit abspielen, anstatt nur die einzelnen Bilder zu zählen.

3. Die Entdeckung: Fehler im Notenblatt

Während Johnson die Musik (die Mathematik) spielte, fiel ihm auf, dass das Notenblatt, das andere Physiker in der Vergangenheit benutzt hatten, Fehler enthielt.

• Die Analogie: Stellen Sie sich vor, Sie spielen ein Klavier, aber in der Partitur steht an einer Stelle ein falsches "Halbe" statt eines "Viertel". Wenn Sie das so spielen, klingt es schief.
• Johnson fand heraus, dass in der Literatur ein paar Faktoren (wie die Zahl 2) falsch geschrieben oder vergessen wurden. Er korrigierte diese "Tippfehler" in den Formeln. Erst danach passte seine Berechnung perfekt zu den theoretischen Erwartungen.

4. Der schwierige Weg: Die Landkarte neu zeichnen

Um die Berechnungen im Computer durchzuführen, musste Johnson die Koordinaten ändern.

• Die Analogie: Stellen Sie sich vor, Sie wollen die Küste eines Landes vermessen. Wenn Sie ein lineares Lineal benutzen (jeder Zentimeter gleich weit), verpassen Sie die winzigen Buchten am Anfang (nahe dem Ursprung), wo die Tänzer sich am engsten umarmen.
• Johnson nutzte daher verschiedene "Landkarten" (mathematische Transformationen). Eine Karte zog die Küste in die Länge, eine andere veränderte den Maßstab logarithmisch. Er stellte fest, dass die "logarithmische Karte" (die er mit 'z' bezeichnet) die genauesten Ergebnisse lieferte, besonders für die komplizierten Tanzschritte, die sehr nah am Zentrum stattfinden.

5. Das Ergebnis: Ein perfekter Tanz

Am Ende hat Johnson die Energie berechnet, die nötig ist, um das Positronium in verschiedenen "angeregten" Zuständen (höhere Tanzschritte) zu halten.

• Er verglich seine Ergebnisse mit den alten, vereinfachten Berechnungen.
• Das Fazit: Seine neuen, direkten Berechnungen stimmen fast perfekt mit den korrigierten alten Formeln überein. Das ist ein großer Erfolg! Es bestätigt, dass die "Zwangsgleichungen" von Dirac ein sehr zuverlässiges Werkzeug sind, um die Quantenwelt zu verstehen.

Zusammenfassung für den Alltag

Dieses Papier ist wie eine Reparatur und Verbesserung eines Navigationsgeräts.

1. Der Autor hat ein altes, komplexes mathematisches Modell (Dirac-Gleichungen) genommen.
2. Er hat Tippfehler in den alten Anleitungen gefunden und korrigiert.
3. Er hat die Berechnungsmethode so angepasst, dass sie auch die schwierigsten Stellen (nahe dem Zentrum) genau erfasst.
4. Das Ergebnis ist eine präzise Vorhersage, wie stark das Positronium zusammenhält.

Das Wichtigste ist: Die Mathematik funktioniert! Wenn man die Gleichungen richtig aufsetzt und die kleinen Fehler in der Literatur entfernt, erhält man ein klares, genaues Bild davon, wie die Materie auf der kleinsten Ebene miteinander "tanzt".

---

Warum ist das wichtig?
Weil Positronium ein perfektes Labor ist, um die Gesetze der Physik zu testen. Wenn unsere Berechnungen hier nicht stimmen, wissen wir, dass unsere ganze Theorie der Quantenelektrodynamik (QED) vielleicht Lücken hat. Johnson hat gezeigt, dass unsere Theorie sehr robust ist, solange wir die Mathematik sorgfältig genug prüfen.

Technische Zusammenfassung

Titel: Angeregte Zustände von Positronium aus den Zweikörper-Dirac-Gleichungen der Nebenbedingungen

Autor: Robert W. Johnson
Quelle: arXiv:2506.21774 (Datiert: 14. April 2026)

---

1. Problemstellung und Motivation

Die Arbeit befasst sich mit der Berechnung der Bindungsenergien angeregter Zustände des Positroniums (ein gebundenes System aus einem Elektron und einem Positron).

• Herausforderung: Bisherige Untersuchungen stützen sich primär auf die nicht-relativistische Quantenelektrodynamik (nrQED). Das Ziel dieses Papers ist es, eine vollständig kovariante, relativistische Beschreibung zu verwenden, die auf den Zweikörper-Dirac-Gleichungen der Nebenbedingungen (Two-Body Dirac Equations of Constraint, TBDE) basiert.
• Lücken in der Literatur: Der Autor identifiziert Inkonsistenzen und vermutlich fehlerhafte Druckfehler (Misprints) in früheren analytischen Formeln und numerischen Ergebnissen (insbesondere in Bezug auf die Arbeit von Crater et al., Ref. [13]), die eine direkte Vergleichbarkeit erschweren.
• Ziel: Entwicklung einer nicht-störungstheoretischen (nonperturbative) numerischen Methode zur Lösung dieser Gleichungen für Systeme gleicher Masse und Vergleich der Ergebnisse mit störungstheoretischen Vorhersagen.

2. Methodik

Theoretischer Rahmen

• TBDE-Formalismus: Die Arbeit nutzt den Formalismus, der von Dirac, Todorov und Crater entwickelt wurde. Für Systeme gleicher Masse ($m_1 = m_2$) entkoppeln bestimmte Zustände (Singulett und Triplett mit $J=L$ sowie der $^3P_0$-Zustand), während andere (Triplett-Zustände mit $L=J\pm1$) gekoppelt bleiben.
• Quasipotential: Die Wechselwirkung wird durch ein effektives Quasipotential $\Phi$ beschrieben, das aus den kinematischen Variablen und den hyperbolischen Funktionen der Potentiale abgeleitet wird. Ein wesentlicher Aspekt ist die Behandlung der Singularitäten und die korrekte Berücksichtigung der Rückstoßeffekte durch die Nebenbedingungen.
• Korrektur von Formeln: Der Autor zeigt auf, dass in der Literatur (Ref. [13]) für den Fall $J=L \ge 1$ ein Faktor $1/2$ in der Formel für $\eta_\pm$ fehlte oder falsch angewendet wurde, was zu Diskrepanzen führte.

Numerische Implementierung

• Diskretisierung: Die Differentialgleichungen werden mittels Finite-Differenzen-Methoden gelöst. Es wird ein Stencil (Gitterstencil) der Ordnung 3 (Länge 7) verwendet.
• Randbedingungen: Besondere Aufmerksamkeit gilt den Randbedingungen bei $r=0$ und $r \to \infty$. Da die Wellenfunktion am Ursprung verschwinden muss, werden spezielle Koeffizienten für die Finite-Differenzen an den Rändern (Edge Correction) berechnet.
• Koordinatentransformationen: Um die numerische Stabilität zu erhöhen und die Singularität am Ursprung zu handhaben, werden verschiedene Koordinatentransformationen getestet:
  • $y = r w / (2\alpha)$ (natürliche Skalierung für die Gleichungen).
  • $z = \log y$ (logarithmische Skala).
  • $x = (1 + s/y)^{-1}$ (kompakte Abbildung auf das Intervall $[0,1]$).
• Eigenwertproblem: Das Problem wird in die Form $H V = D V$ überführt, wobei $D$ eine Diagonalmatrix ist. Zur Lösung wird die ARPACK-Bibliothek verwendet.
• Iterative Verbesserung: Da die Energie $w$ in den Potentialtermen vorkommt, wird ein iteratives Verfahren angewendet: Zuerst wird mit einer Näherung ($w \approx 2m$) gelöst, dann wird der erhaltene Eigenwert in die Hamilton-Matrix zurückgespeist und die Lösung verfeinert, bis Konvergenz erreicht ist.

3. Wichtige Beiträge und Entdeckungen

1. Identifikation und Korrektur von Druckfehlern: Der Autor findet und korrigiert signifikante Fehler in der Literatur (Ref. [13]), insbesondere bei den analytischen Formeln für die Störungstheorie. Ohne die Korrektur eines Faktors $1/2$ in den Formeln für $\eta_\pm$ stimmen die numerischen Ergebnisse nicht mit den veröffentlichten Tabellen überein.
2. Open-Source-Implementierung: Die numerische Methode wurde als Open-Source-Toolbox für die GNU Octave-Umgebung implementiert und ist öffentlich verfügbar. Dies ermöglicht eine Reproduzierbarkeit der Ergebnisse.
3. Vergleich verschiedener Koordinatensysteme: Es wird gezeigt, dass die Wahl der Koordinatentransformation ($x, y, z$) einen signifikanten Einfluss auf die numerische Genauigkeit hat. Die logarithmische Koordinate $z$ und die kompakte Koordinate $x$ liefern robustere Ergebnisse als die direkte radiale Koordinate $y$.
4. Vergleich mit der verbesserten Breit-Gleichung: Die Ergebnisse werden mit einer singulärfreien Breit-Gleichung verglichen. Es wird festgestellt, dass eine bestimmte Umformung der Gleichung (Verschiebung eines Faktors $w$) notwendig ist, um die physikalisch korrekten Bindungsenergien zu erhalten, was auf subtile mathematische Äquivalenzen in der Formulierung hinweist.

4. Ergebnisse

• Bindungsenergien: Die nicht-störungstheoretisch berechneten Bindungsenergien für die angeregten Zustände des Positroniums stimmen hervorragend mit den korrigierten störungstheoretischen Werten überein.
• Genauigkeit: Die Übereinstimmung liegt im Bereich von $10^{-8}$ bis $10^{-10}$ (abhängig vom Zustand und der verwendeten Koordinate), was die Validität des TBDE-Formalismus für dieses System unterstreicht.
• Tabelle VI & VII: Die Tabellen zeigen detaillierte Werte für verschiedene Quantenzahlen ($J, S, L, n$). Die Residuen (Fehlermaße) $\chi_{10}$ und $\Delta_{10}$ belegen die hohe numerische Präzision, insbesondere bei Verwendung der Koordinaten $x$ und $z$.
• Große Kopplungskonstante: Im Anhang B werden Ergebnisse für große Werte der Feinstrukturkonstante $\alpha$ (bis $\alpha \approx 0.74$) präsentiert, wo die Bindungsenergie mehr als die Hälfte der Ruhemasse beträgt. Dies demonstriert die Fähigkeit der Methode, auch im stark gekoppelten Regime zu funktionieren.

5. Bedeutung und Fazit

Dieses Papier leistet einen wesentlichen Beitrag zur Präzisionsphysik des Positroniums, indem es:

• Eine robuste, nicht-störungstheoretische numerische Methode etabliert, die unabhängig von der Störungsreihe ist.
• Versteckte Fehler in der etablierten Literatur aufdeckt und korrigiert, was für zukünftige Vergleiche mit experimentellen Daten (z.B. zur Suche nach Physik jenseits des Standardmodells) essenziell ist.
• Die Äquivalenz und Konsistenz zwischen dem TBDE-Formalismus und der verbesserten Breit-Gleichung bestätigt.
• Durch die Bereitstellung des Codes die wissenschaftliche Gemeinschaft in die Lage versetzt, diese Berechnungen selbstständig zu überprüfen und zu erweitern.

Der Autor schließt, dass die TBDE-Gleichungen eine zuverlässige und kovariante Alternative zu nrQED für die Berechnung von gebundenen Zuständen in der QED darstellen, sofern die numerischen Implementierungen sorgfältig auf Randbedingungen und Koordinatentransformationen optimiert werden.