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

Technische Samenvatting: Stationaire Perturbatietheorie zonder Sommaties over Tussenliggende Toestanden: Supersymmetrisch Expansie-algoritme

Probleemstelling
Het artikel behandelt een fundamentele beperking in de conventionele Rayleigh-Schrödinger perturbatietheorie (RSPT). Hoewel RSPT een systematisch kader biedt voor het berekenen van correcties aan energie-eigenwaarden en eigenfuncties, is het afhankelijk van sommaties over tussenliggende toestanden. De eerste-orde correctie aan eigenfuncties vereist al een "onhandige som", en hogere-orde correcties vereisen steeds complexere sommaties over tussenliggende toestanden bij elke orde. Deze formalisering wordt computationeel omslachtig en vaak onpraktisch voorbij de laagste orden. Eerdere alternatieven, zoals somregels of de logaritmische formulering van de Schrödingervergelijking voorgesteld door Aharonov en Au, hebben er succesvol in geslaagd om berekeningen voor de grondtoestand te reduceren tot kwadratuurvormen (integralen) voor eendimensionale systemen, maar blijven minder handig voor aangeslagen toestanden vanwege de aanwezigheid van knopen.

Methodologie: Het Supersymmetrisch Expansie-algoritme (SEA)
De auteurs stellen het Supersymmetrisch Expansie-algoritme (SEA) voor om perturbatieve correcties van willekeurige orde voor alle eigenfuncties (grond- en aangeslagen toestanden) te verkrijgen zonder te sommeren over tussenliggende toestanden. De methode maakt gebruik van de verborgen supersymmetrie die inherent is aan niet-relativistische kwantummechanische problemen.

De kernprocedure omvat de volgende stappen:

1. Logaritmische Formulering: De Schrödingervergelijking wordt gegoten in een niet-lineaire logaritmische vorm met een superpotentiaal $W(x, \lambda) = -\frac{d}{dx} \ln u(x, \lambda)$, waarbij $u$ de eigenfunctie is.
2. Perturbatieve Expansie: De superpotentiaal, de potentiaal en de energie worden geëxpandeerd in machten van een kleine parameter $\lambda$. Dit transformeert de niet-lineaire differentiaalvergelijking in een hiërarchische cascade van gekoppelde eerste-orde lineaire differentiaalvergelijkingen voor de expansiecoëfficiënten.
3. Knoploze "Randtoestanden" (Edge States): Het algoritme lost eerst de knoplote eigenfunctie (de "randtoestand") van de Hamiltoniaan op een gegeven supersymmetrisch niveau op. Voor eendimensionale systemen is dit de grondtoestand; voor hogere dimensies die reduceerbaar zijn tot effectieve eendimensionale potentialen (zoals centrale potentialen), zijn dit de toestanden met het hoogste impulsmoment $l$ voor een gegeven hoofdkwantumgetal $n$.
4. Recursieve Oplossing via Supersymmetrie:
   • De knoplote oplossing voor de $r$-de supersymmetrische partner-Hamiltoniaan $H_r$ wordt opgelost met behulp van de logaritmische expansie.
   • De leidende orde term voor de knoplote oplossing van $H_r$ wordt afgeleid van de bekende ongestoorde oplossingen van de oorspronkelijke Hamiltoniaan $H_0$ door recursief de inverse van de supersymmetrische verhogingsoperatoren ($a^\dagger$) toe te passen.
   • Zodra de knoplote oplossing voor $H_r$ is vastgesteld, worden de perturbatieve correcties voor de energie en superpotentiaal berekend met behulp van integralen gewogen door de ongestoorde waarschijnlijkheidsdichtheid van die specifieke knotloze toestand.
   • De eigenfuncties van de oorspronkelijke Hamiltoniaan $H_0$ met $r$ knopen worden vervolgens gereconstrueerd door de keten van verhogingsoperatoren ($a^\dagger_0, a^\dagger_1, \dots$) toe te passen op de knoplote oplossing van $H_r$.

Belangrijkste Bijdragen en Resultaten

• Eliminatie van Tussenliggende Toestandssommen: De primaire bijdrage is de reductie van perturbatieve correcties naar "kwadratuurvormen". Energiecorrecties en eigenfunctiecorrecties worden uitgedrukt als integralen gewogen door de waarschijnlijkheidsdichtheden van de relevante randtoestanden, waarbij de sommatie over tussenliggende toestanden die vereist is in RSPT volledig wordt vermeden.
• Willekeurige Orde en Algemene Toepasbaarheid: De formalisering is toepasbaar op willekeurige orden van de perturbatietheorie en breidt zich uit naar alle eigenfuncties, inclusief aangeslagen toestanden met knopen. Het werkt voor zowel polynomiale als niet-polynomiale potentialen.
• Analytische Oplossingen voor Specifieke Potentialen: De auteurs demonstreren de effectiviteit van de methode door bekende analytische oplossingen te reproduceren voor de Yukawa-, Hulthén-, anharmonische- en Cornell-potentialen, waar de superpotentiaal-coëfficiënten polynomiale functies zijn.
• Niet-Polynomiale Voorbeeld: Het artikel past de methode toe op een deeltje in een doos met een harmonische oscillator-perturbatie. Dit geval omvat niet-polynomiale superpotentialen. De auteurs hebben de methode succesvol toegepast om de energiecorrecties tot orde $\lambda^3$ af te leiden voor zowel de grondtoestand als de eerste aangeslagen toestand, waarbij expliciete analytische expressies voor de energie-eigenwaarden worden geleverd.

Betekenis en Claims
Het artikel claimt dat de SEA een "volledige analytische oplossing" biedt aan perturbatieproblemen in de vorm van een machtsreeks in de expansieparameter. Door gebruik te maken van de verborgen supersymmetrie van kwantumsystemen, biedt de methode een transparant en systematisch alternatief voor RSPT. De auteurs benadrukken dat, hoewel eerdere methoden (zoals die van Aharonov en Au) beperkt waren in het behandelen van aangeslagen toestanden, de SEA de oplossing van aangeslagen toestanden natuurlijk verbindt met de supersymmetrische hiërarchie, waardoor de berekening van correcties voor alle toestanden eenvoudig wordt. Het werk stelt dat deze benadering bijzonder krachtig is voor systemen waarbij de coëfficiënten in de potentiaalexpansie geen polynomialen zijn, omdat het niet vertrouwt op het reduceren van het probleem tot algebraïsche recursierelaties, maar de differentiaalhiërarchie direct oplost.