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

Dit artikel berekent de bindingsenergieën van de aangeslagen toestanden van positronium met behulp van de twee-deeltjes Dirac-vergelijkingen van beperking, vergelijkt niet-perturbatieve en perturbatieve resultaten, en corrigeert fouten in de bestaande literatuur.

De kern

Het Dansende Elektron-Positron Koppel: Een Verkenning van de Positronium-Geheime Codes

Stel je voor dat je twee danspartners hebt: een elektron (negatief geladen) en een positron (het positieve tegenhanger van het elektron). Als ze elkaar vinden, dansen ze niet zomaar rond; ze vormen een tijdelijk koppel dat we Positronium noemen. Ze draaien om elkaar heen, net als de aarde om de zon, maar dan veel sneller en op een heel klein niveau.

Deze dans is echter niet perfect. Ze verliezen energie, en soms slaan ze elkaar dood (annihilatie) en verdwijnen ze als lichtflits. De vraag die natuurkundigen zich stellen is: Hoe precies is de energie van deze dans? Hoe snel draaien ze in de verschillende "stappen" (excited states)?

Dit artikel van Robert Johnson is als een rekenmachine voor danspassen. Hij gebruikt een heel specifieke, oude maar krachtige wiskundige methode (de "Two Body Dirac Equations") om te voorspellen hoe deze dans eruit ziet, en vergelijkt dit met andere manieren om het te berekenen.

Hier zijn de belangrijkste punten, vertaald naar alledaagse taal:

1. De Oude Kaart vs. De Nieuwe GPS

Vroeger gebruikten wetenschappers een methode genaamd nrQED (een soort benadering) om de energie van Positronium te berekenen. Het is alsof je een oude, handgetekende kaart gebruikt om een route te plannen: het werkt goed voor de grote lijnen, maar mist soms de kleine bochten.

Johnson gebruikt een andere methode: de Twee-Lichaam Dirac-vergelijkingen. Dit is alsof je een ultra-precieze GPS gebruikt die rekening houdt met elke kleine helling en bocht in de weg, inclusief de zwaartekracht van de auto zelf. Deze methode is "niet-perturbatief", wat betekent dat hij de hele dans in één keer bekijkt in plaats van hem op te splitsen in kleine stukjes.

2. De "Rekenfout" in de Boekjes

Tijdens het werk ontdekte Johnson een grappig maar belangrijk ding: er staan fouten in de oude wiskundige formules die andere wetenschappers gebruikten.

• De Analogie: Stel je voor dat je een recept voor een taart volgt, maar in het boekje staat per ongeluk "100 gram suiker" in plaats van "50 gram". Als je het recept volgt, wordt je taart te zoet.
• Johnson zag dat als hij de formules letterlijk nam, de resultaten niet klopten met de verwachtingen. Hij ontdekte dat er een factor 2 ontbrak of ergens verkeerd stond. Door dit te corrigeren (de "receptfout" op te lossen), kwamen de oude formules eindelijk overeen met de nieuwe, precieze berekeningen.

3. Het Moeilijke Rekenen: Van Rijen naar Logaritmen

Om deze dans te simuleren, moet je een computer laten rekenen met oneindig veel punten. Maar de ruimte rondom de deeltjes is lastig: heel dicht bij elkaar is het heel druk, en verder weg is het rustig.

• De Analogie: Stel je voor dat je een foto maakt van een stad. Als je alles in één keer vastlegt, zie je de details van de gebouwen niet, maar wel de hele stad. Als je inzoomt op één gebouw, zie je de details, maar dan is de rest van de stad weg.
• Johnson gebruikt slimme coördinaattransformaties. Hij "vervormt" de ruimte in zijn computer zo, dat hij heel veel detail kan zien waar de deeltjes dicht bij elkaar zijn (dicht bij de oorsprong), zonder dat de computer vastloopt. Hij gebruikt hiervoor verschillende "lenzen" (coördinaten genaamd x, y en z) om te kijken welke lens het scherpste beeld geeft.

4. De Resultaten: Het Koppel is in Evenwicht

Na al het rekenen en het corrigeren van de fouten, concludeert Johnson:

• Zijn nieuwe, precieze methode geeft exact dezelfde resultaten als de oude methoden (als je de fouten in de oude formules eerst corrigeert).
• Dit betekent dat de oude theorieën eigenlijk wel goed waren, maar dat we ze moesten "kalibreren".
• Hij heeft ook gekeken naar wat er gebeurt als de deeltjes heel sterk aan elkaar trekken (bijvoorbeeld bij zware deeltjes in deeltjesversnellers), en daar werkt zijn methode ook goed.

5. Waarom is dit belangrijk?

Dit klinkt misschien als droge wiskunde, maar het is cruciaal voor de fundamentele natuurkunde.

• Als we de energie van Positronium niet perfect kunnen voorspellen, kunnen we niet zeker weten of er nieuwe deeltjes of nieuwe krachten in het universum zijn die we nog niet kennen.
• Het is alsof je een weegschaal hebt die tot op de microgram nauwkeurig is. Als je die weegschaal niet goed hebt afgesteld (de rekenfouten), denk je dat er een onbekend gewicht op ligt, terwijl het eigenlijk alleen maar een verkeerde instelling was.

Samenvatting

Robert Johnson heeft een oude, complexe wiskundige methode gebruikt om de dans van een elektron en een positron te simuleren. Hij ontdekte dat er rekenfouten in de bestaande formules stonden die de resultaten verdraaiden. Door deze fouten te repareren en slimme computertrucs te gebruiken om de ruimte te "vervormen", heeft hij bewezen dat de theorieën kloppen. Hij heeft zijn code zelfs gratis beschikbaar gesteld (open source) zodat iedereen het kan controleren en verder kan bouwen.

Het is een mooi voorbeeld van hoe wetenschap werkt: niet alleen door nieuwe dingen te bedenken, maar ook door te kijken of de oude dingen wel echt kloppen, en ze dan netjes op te poetsen.

Technische samenvatting

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

Auteur: Robert W. Johnson (Alphawave Research)
Datum: 14 april 2026 (arXiv:2506.21774)

---

1. Het Probleem

De interpretatie van experimentele waarnemingen van het positronium-spectrum (een gebonden systeem van een elektron en een positron) is momenteel grotendeels afhankelijk van niet-relativistische kwantumelektrodynamica (nrQED). Hoewel nrQED succesvol is, is klassieke elektrodynamica per constructie covariënt, en zou een kwantumformulering dit ook moeten zijn.

Er bestaat een alternatieve theorie, gebaseerd op de Twee-Lichaams Dirac-vergelijkingen van Beperking (Two-Body Dirac Equations of Constraint - TBDE), ontwikkeld door Dirac, Todorov en later Crater et al. Deze theorie biedt een covariante beschrijving van gebonden toestanden zonder de singulariteiten die vaak voorkomen in andere benaderingen. Echter, de toepassing van deze theorie op de aangeslagen toestanden van positronium vereist nauwkeurige numerieke implementatie en verificatie tegen bestaande perturbatieve resultaten. Er zijn ook onzekerheden en mogelijke drukfouten in de literatuur over de analytische formules voor het spectrum.

2. Methodologie

De auteur implementeert de TBDE-theorie voor systemen met gelijke massa (QED) en lost de eigenwaardeproblemen numeriek op. De aanpak omvat de volgende stappen:

• Formulering: De vergelijkingen worden opgesteld voor singlet- en triplettoestanden. Voor gelijke massa's ($m_1 = m_2$) decoupleren bepaalde toestanden (zoals $^1J_J$, $^3J_J$ en $^3P_0$), terwijl andere ($^3L_{L+1}$ en $^3L_{L-1}$) gekoppeld blijven.
• Numerieke Methode:
  • De differentiaalvergelijkingen worden discretiseerd met behulp van finite difference methoden (5e orde voor de eerste afgeleide, 7-punts stencil voor de tweede afgeleide).
  • Om de singulariteit bij de oorsprong ($r=0$) en het gedrag op oneindig goed te behandelen, worden verschillende coördinatentransformaties toegepast:
    • $y = rw/2\alpha$ (lineaire schaling).
    • $z = \log y$ (logaritmische schaling).
    • $x = (1 + s/y)^{-1}$ (gecomprimeerde coördinaat).
  • Randvoorwaarden worden zorgvuldig behandeld met "edge correction" voor de finite difference stencil.
• Oplossing: Het eigenwaardeprobleem wordt opgelost met de ARPACK bibliotheek.
• Iteratieve Verbetering: Omdat de energie $w$ in de potentiaaltermen voorkomt, wordt een iteratief proces gebruikt. Eerst wordt $w$ benaderd als $2m_1$, waarna de oplossing wordt gebruikt om de potentiaal te updaten en de eigenwaarde te verfijnen totdat convergentie is bereikt.
• Validatie: De numerieke resultaten worden vergeleken met perturbatieve analytische formules uit de literatuur (Crater et al., 1992) en met een verbeterde Breit-vergelijking.

3. Belangrijkste Bijdragen

• Open Source Implementatie: De auteur heeft de numerieke code beschikbaar gesteld als een toolbox voor de GNU Octave-omgeving, inclusief afleidingen in wxMaxima-formaat.
• Identificatie en Correctie van Fouten in de Literatuur:
  • Er werd een discrepantie gevonden tussen de gepubliceerde analytische formules in referentie [13] en de numerieke resultaten.
  • De auteur identificeerde een drukfout in de formule voor de parameter $\eta_{\pm}$ voor toestanden met $J = L \geq 1$. Door een factor $1/2$ toe te passen op de uitdrukking voor $\eta_{\pm}$, kwamen de perturbatieve resultaten perfect overeen met de numerieke berekeningen.
  • Er werden ook inconsistenties in tekens en factoren van 2 in eerdere publicaties opgemerkt en gecorrigeerd.
• Vergelijking van Coördinaten: Er is een grondige analyse uitgevoerd van de nauwkeurigheid van verschillende coördinatentransformaties ($x, y, z$) voor het oplossen van de vergelijkingen.

4. Resultaten

• Niet-perturbatieve Berekeningen: De bindingenergieën voor diverse aangeslagen toestanden van positronium (tot en met $n=3$) zijn berekend.
  • De resultaten tonen een uitstekende overeenkomst met de (gecorrigeerde) perturbatieve waarden.
  • De nauwkeurigheidsmeting ($\Delta_{10}$, het logaritmische verschil met de perturbatieve waarde) toont aan dat de logaritmische coördinaat ($z$) en de gecomprimeerde coördinaat ($x$) de meest betrouwbare resultaten leveren. De lineaire coördinaat ($y \propto r$) is minder nauwkeurig, vooral voor toestanden met $L=0$.
• Vergelijking met de Breit-vergelijking: De resultaten van de TBDE worden vergeleken met een "verbeterde Breit-vergelijking".
  • Er werd ontdekt dat een directe oplossing van de afgeleide tweede-orde vergelijking (Eqn. 63) leidt tot fysisch onjuiste resultaten (de grondtoestand ontbreekt, en de bindingenergieën zijn ongeveer de helft).
  • Door een factor $w$ van links naar rechts te verplaatsen (Eqn. 64), worden de correcte bindingenergieën hersteld. Dit suggereert dat de matrixoperaties in de oorspronkelijke vorm een ongewenste wortel van de energie impliceren.
• Groot $\alpha$ Regime: In de appendix worden resultaten getoond voor grote waarden van de koppelingsconstante $\alpha$ (tot 0.74), waarbij de bindingenergie meer dan de helft van de rustmassa bedraagt, wat de stabiliteit van de numerieke methode demonstreert.

5. Betekenis

Dit paper levert een belangrijke bijdrage aan de theoretische fysica van gebonden toestanden in de QED:

1. Verificatie van TBDE: Het bevestigt dat de TBDE-formulering, die volledig covariënt is en geen singulariteiten kent, numeriek stabiel is en nauwkeurige resultaten oplevert voor het positronium-spectrum.
2. Correctie van de Wetenschap: Het corrigeert specifieke fouten in de bestaande literatuur, waardoor toekomstige vergelijkingen tussen theorie en experiment nauwkeuriger kunnen worden uitgevoerd.
3. Methodologische Inzicht: Het biedt inzicht in de numerieke valkuilen bij het oplossen van relativistische twee-lichaamsproblemen, met name de keuze van coördinaten en de behandeling van de energie-afhankelijke potentiaal.
4. Toegankelijkheid: Door de code open source te maken, maakt de auteur deze geavanceerde numerieke methode toegankelijk voor andere onderzoekers om de theorie verder te testen en uit te breiden.

Samenvattend presenteert dit werk een robuuste, niet-perturbatieve berekening van het positronium-spectrum, valideert deze tegenover perturbatieve theorieën na correctie van literatuurfouten, en biedt een gedetailleerde technische handleiding voor de numerieke implementatie van de Two-Body Dirac-vergelijkingen.