Laguerre-Sobolev orthogonal Polynomials and Boundary Value Problems on a semi-infinite domain

Technische Samenvatting: Laguerre–Sobolev Orthogonale Polynomen en Randwaardeproblemen op een Semi-Oneindig Domein

Probleemstelling
Het artikel behandelt de numerieke oplossing van tweede-orde randwaardeproblemen (BVP's) op het semi-oneindige interval $(0, +\infty)$, met specifieke aandacht voor stationaire Schrödinger-vergelijkingen met singuliere potentialen, zoals de inverse-afstand potentiaal $V(x) = 1/x$. Hoewel klassieke spectrale methoden gebaseerd op Laguerre-polynomen goed gevestigd zijn voor gladde oplossingen op semi-oneindige domeinen, kunnen zij inefficiënt worden bij het omgaan met randvoorwaarden in oneindigheid of singuliere coëfficiënten in de differentiaaloperator. De auteurs streven ernaar een raamwerk te ontwikkelen dat afgeleide informatie en randgedrag direct in de benaderingsruimte incorporeert om deze singulariteiten effectiever aan te pakken.

Methodologie
De auteurs introduceren een familie van Laguerre–Sobolev orthogonale polynomen, aangeduid als $\{S_n(x)\}_{n \geq 0}$, geassocieerd met een specifieke Sobolev-inwendig product afgeleid van de variationele formulering van het doel-BVP. Het inwendig product wordt gedefinieerd als:
$$\langle u, v \rangle_\lambda = \lambda \int_0^{+\infty} \frac{u(x)v(x)}{x} dx + \int_0^{+\infty} u'(x)v'(x) dx$$
De testfuncties worden gekozen uit de ruimte $\mathcal{P} \times x e^{-x/2}$, wat garandeert dat zij in beide eindpunten ($0$ en $+\infty$) verdwijnen.

De kern van de methodologie omvat:

1. Constructie: Het genereren van de orthogonale basis $\{S_n(x)\}$ via het Gram-Schmidt proces toegepast op de sequentie $\{x^k x e^{-x/2}\}$ met betrekking tot het Sobolev-inwendig product.
2. Analytische Afleiding: Het vaststellen van expliciete connectieformules tussen de nieuwe Sobolev-polynomen $S_n(x)$ en de klassieke Laguerre-polynomen $L_n^{(1)}(x)$. Dit omvat het afleiden van recursie-relaties voor de connectiecoëfficiënten en het analyseren van hun asymptotische gedrag met behulp van verfijnde ratio-asymptotica voor Laguerre-polynomen.
3. Genererende Functie: Het afleiden van een gesloten vorm van een genererende functie voor de Sobolev-polynomen die betrokken zijn bij Bessel-functies.
4. Spectraal Algoritme: Het ontwikkelen van een volledig gedegenereerde spectrale methode. Door gebruik te maken van de connectieformules worden de Fourier–Sobolev-coëfficiënten van de oplossing recursief berekend. Deze aanpak vermijdt het oplossen van lineaire stelsels, aangezien het resulterende algoritme enkel de evaluatie van integralen betreffende klassieke Laguerre-polynomen vereist.

Belangrijkste Bijdragen

• Nieuwe Orthogonale Familie: De definitie en analyse van een specifieke familie van Laguerre–Sobolev orthogonale polynomen die op maat is gemaakt voor BVP's met singuliere potentialen op $(0, +\infty)$.
• Connectieformules: De afleiding van expliciete relaties die $S_n(x)$ koppelen aan $L_n^{(1)}(x)$, inclusief een recursie-relatie voor de connectiecoëfficiënten $a_n$ en hun asymptotische expansie ($a_n \sim 1 - \sqrt{4\lambda/n}$).
• Genererende Functie: De voorziening van een genererende functie voor de Sobolev-polynomen uitgedrukt in termen van Bessel-functies van de eerste soort.
• Gedegenereerd Algoritme: De formulering van een recursieve spectrale methode die de expansiecoëfficiënten van de oplossing berekent zonder lineaire stelsels op te lossen, waardoor de computationele complexiteit wordt verminderd terwijl de stabiliteit voor de specifieke klasse van singuliere problemen behouden blijft.

Resultaten
Het artikel presenteert numerieke experimenten voor het oplossen van een niet-homogene Schrödinger-vergelijking met een inverse-afstand potentiaal ($V(x)=1/x$) en Dirichlet-randvoorwaarden.

• Exponentiële Convergentie: Voor een oplossing die exponentiële afname vertoont ($u(x) = x \cos x e^{-x}$), demonstreert de methode spectrale nauwkeurigheid en exponentiële convergentie. De fout in de Sobolev-norm neemt snel af naarmate het aantal termen toeneemt.
• Beperkingen: De auteurs merken op dat voor oplossingen die geen exponentiële afname vertonen (bijv. rationele afname), de prestaties van de benaderingen slechter zijn, zoals verwacht voor spectrale methoden op semi-oneindige domeinen.
• Efficiëntie: De recursieve aard van het algoritme maakt computationele kosten mogelijk die vergelijkbaar zijn met standaard Laguerre-methoden, terwijl het een verbeterde stabiliteit biedt voor de specifieke klasse van de behandelde singuliere problemen.

Betekenis
Het artikel beweert dat de voorgestelde Laguerre–Sobolev spectrale methode een waardevol alternatief biedt voor klassieke Laguerre-methoden voor problemen met singuliere coëfficiënten en randvoorwaarden in oneindigheid. Door afgeleide-informatie en randgedrag direct in de benaderingsruimte in te bedden, bereikt de methode een volledig gedegenereerd stelsel. Dit elimineert de noodzaak om lineaire stelsels op te lossen, wat de computationele complexiteit aanzienlijk vermindert. Het werk biedt een rigoureuze analytische fundering (connectieformules, asymptotica, genererende functies) die de praktische implementatie van deze methoden voor Schrödinger-type problemen met singuliere potentialen mogelijk maakt.