Novel method for evaluating the eigenvalues of the Heun differential equation with an application to the Breit equation

Technische Samenvatting: Een Nieuwe Methode voor het Evalueren van de Eigenwaarden van de Heun-differentiaalvergelijking met een Toepassing op de Breit-vergelijking

Probleemstelling
De Breit-vergelijking, die de interactie tussen twee fermionen (bijv. een elektron en een proton) via een statisch Coulomb-potentiaal beschrijft, leidt tot een tweede-orde differentiaalvergelijking met vier singuliere punten. Dit classificeert de vergelijking als een Heun-differentiaalvergelijking, die meer singulier is dan de standaard geconfluente hypergeometrische vergelijkingen die gewoonlijk in de wiskundige natuurkunde worden aangetroffen. Hoewel eerdere studies (bijv. Yamaguchi; Tursunov et al.) de Breit-vergelijking hebben geanalyseerd met behulp van variatiemethoden of door series-oplossingen te matchen, bestaat er momenteel geen bepaling van energieniveaus die direct gebaseerd is op de tweede-orde Heun-vergelijking afgeleid van de Breit-vergelijking. Bovendien bestaat er geen expliciete uitdrukking voor de expansie van energieniveaus in machten van de fijnstructuurconstante gekwadrateerd ($\alpha^2$) voorbij de eerste orde.

Methodologie
De auteurs ontwikkelen een theoretisch kader om de Breit-vergelijking op te lossen voor spin-singlet ($S=0$) en specifieke spin-triplet ($S=1$) toestanden. De methodologie verloopt via twee afzonderlijke, wederzijds validerende benaderingen:

1. Continued Fraction (CF) Benadering:

   • De radiale Breit-vergelijking wordt getransformeerd naar een standaardvorm van de geconfluente Heun-vergelijking.
   • De eigenfuncties worden geëxpandeerd in termen van geassocieerde Laguerre-polynomen.
   • Deze expansie levert een drie-term recursierelatie voor de expansiecoëfficiënten op.
   • De recursierelatie wordt omgezet in een continued fraction (doorlopende breuk). De eigenwaarden worden bepaald door de nulpunten van deze continued fraction te vinden.
   • Numerieke evaluatie maakt gebruik van Miller's backward recurrence algoritme voor de normalisatie van de golffunctie en Peter Wynn's $\rho$-algoritme (een niet-lineaire sequentie-versneller) om de trage convergentie van de CF nabij eigenwaarden te behandelen.

2. Green Function en Semi-Infinite Determinant Benadering:

   • De differentiaalvergelijking wordt geherformuleerd als een homogene Fredholm integraalvergelijking van de tweede soort met behulp van een Green-functie, geconstrueerd uit de eigenfuncties van de bijbehorende Laguerre-operator.
   • Het expanderen van de oplossing van deze integraalvergelijking in geassocieerde Laguerre-polynomen leidt tot een conditie op een semi-oneindig determinant ($\det \Delta_{m,n} = 0$).
   • De auteurs introduceren specifieke speciale functies ($I_{m,n}$, $S^m_n$, $\hat{S}^m_n$) om de determinantelementen te vereenvoudigen.
   • Om numerieke stabiliteit en convergentie te garanderen, wordt een "convergerende factor" geïntroduceerd om de diagonale elementen van de determinant te modificeren.
   • De determinant wordt numeriek opgelost met behulp van drie methoden: directe berekening (met behulp van Dodgson's condensatiemethode), forward recursion, en de CF-methode voor kruisverificatie.

Kernbijdragen

• Afleiding van de Heun-vergelijking: Het artikel leidt expliciet de tweede-orde Heun-differentiaalvergelijking af uit de Breit-vergelijking voor het spin-singlet geval ($m=M$) en generaliseert de parametrisering voor diverse spin-triplet toestanden.
• Nieuwe Determinant-methode: De auteurs presenteren een nieuwe methode voor het berekenen van eigenwaarden gebaseerd op een semi-oneindig determinant afgeleid van een Green-functie integraalvergelijking. Dit biedt een alternatief voor de continued fraction methode.
• Hoge-precisie Numerieke Resultaten: De studie berekent energieniveaus voor de eerste vijf hoofdkwantumgetallen ($N=1$ tot $5$) voor spin-singlet en spin-triplet toestanden ($J=0$ en $J=L$). De resultaten bereiken een nauwkeurigheid tot 25–30 decimalen.
• Consistentiecontrole: De twee verschillende wiskundige benaderingen (CF en Determinant) blijken consistent binnen 25 decimalen van nauwkeurigheid, wat de numerieke methoden valideert.
• Reproductie van Bekende Fysica: De methoden reproduceren succesvol bekende resultaten uit de literatuur (bijv. Scott et al., Yamaguchi et al.) en identificeren correct het "spurious state" gedrag nabij $\rho \sim 1/y$ dat eerder door Patterson werd besproken.
• Dirac-limiet: Het kader toont aan dat de Dirac-vergelijking oplossing kan worden teruggewonnen uit de determinantconditie via een geschikte mapping van variabelen, wat inzicht geeft in de structuur van de determinant.

Resultaten

• Eigenwaarden: Het artikel geeft getabelde waarden voor de eigenwaarden ($z_n$) en energieniveaus ($E_n$) voor de eerste vijf niveaus. Voor het spin-singlet geval komen de resultaten overeen met literatuurwaarden tot 18–25 decimalen.
• Golffuncties: De radiale golffuncties ($F, K, G$) worden efficiënt berekend met behulp van Miller's backward recurrence algoritme. Het artikel visualiseert deze componenten, wat de verwachte asymptotische gedraging bevestigt.
• Convergentieanalyse: De studie analyseert de convergentie van de series-expansie, bewijst absolute convergentie voor $\rho \ge 0$. Het behandelt ook de numerieke uitdagingen van de determinant-methode voor grote fysieke waarden van de parameter $y$ (gerelateerd aan de fijnstructuurconstante), waarbij wordt aangetoond dat forward recursion superieur is aan directe berekening voor hoge precisie in dit regime.
• Analytische Structuur: De analyse van de determinant suggereert dat de series-expansie van eigenwaarden niet-analytische termen kan bevatten, specifiek logaritmen van de fijnstructuurconstante ($\ln(1/\alpha)$), wat contrasteert met sommige bestaande literatuur maar overeenkomt met recente hogere-orde berekeningen voor positronium.

Betekenis en Claims
De auteurs beweren dat hun werk eerdere berekeningen met enkele decimalen extra nauwkeurigheid heeft uitgebreid en een robuust, alternatief theoretisch kader biedt voor het oplossen van de Breit-vergelijking. De primaire betekenis ligt in:

1. Het bieden van een directe methode om eigenwaarden te berekenen op basis van de Heun-vergelijking zonder te vertrouven op het matchen van series-oplossingen bij grenzen.
2. Het aanbieden van een semi-oneindig determinant-conditie die dient als een rigoureuze consistentiecontrole voor de continued fraction methode.
3. Het leggen van de theoretische grondslag voor toekomstig werk om expliciete expansies van energieniveaus in machten van $\alpha^2$ (voorbij de eerste orde) te verkrijgen en potentieel de volledige Breit-interactie (fotonuitwisseling) in te sluiten om de ordening van toestanden op te lossen (bijv. de omkering van $2^1P_1$ en $2^1S_0$ toestanden).

Het artikel blijft bescheiden over de analytische expansie en merkt op dat hoewel de determinant-structuur de aanwezigheid van $\ln(1/\alpha)$ termen suggereert, een expliciete expansie in $\alpha^2$ nog niet succesvol is afgeleid in dit werk, maar een doel is voor toekomstig onderzoek.