Explicit Two-Sided Eigenvalue Bounds for Schrödinger Operators with Singular Potentials via Finite Element Method

Technische Samenvatting: Expliciete tweezijdige eigenwaardegrenzen voor Schrödingervergelijkingen met singuliere potentialen via de Finite Element-methode

Probleemstelling
Het artikel behandelt de berekening van expliciete, berekenbare tweezijdige grenzen voor de eigenwaarden $\lambda_k$ van de Schrödingervergelijking $H = -\Delta + V$ op onbegrensde domeinen $\mathbb{R}^N$ ($N=2, 3$). De primaire uitdaging ligt in het hanteren van singuliere aantrekkende potentialen van het Coulomb-type, specifiek $V(x) = -Z/|x|$ (die het waterstofatoom modelleert) en tweegestelde potentialen $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$ (die het $H_2^+$-moleculaire ion modelleert).

Bestaande numerieke methoden leveren doorgaans bovengrenzen via de conformerende Galerkin (Rayleigh–Ritz) methode. Het verkrijgen van rigoureuze ondergrenzen is aanzienlijk moeilijker, vooral voor onbegrensde domeinen en singuliere potentialen waarbij de potentiaal onbegrensd naar beneden gaat. Klassieke ondergrenstechnieken (bijvoorbeeld Lehmann–Maehly–Goerisch, Kato–Temple) vereisen vaak niet-conformerende methoden of specifieke voorwaarden die falen wanneer de potentiaal $V$ negatief en singulier is. Bovendien vertrouwen standaard finite element-benaderingen voor ondergrenzen vaak op het verschuiven van de potentiaal met de $L^\infty$-norm van het negatieve deel ($\gamma_h = \|c_h^-\|_{L^\infty}$). Voor Coulomb-singulariteiten divergeert $\gamma_h$ naarmate het rooster wordt verfijnd ($O(h^{-1})$), waardoor de ondergrens faalt in te convergeren.

Methodologie
Het voorgestelde algoritme combineert domeintroncatie met een uitbreiding van de Composite Enriched Crouzeix–Raviart (CECR) finite element-methode. De aanpak is gestructureerd in drie hoofdcomponenten:

1. Domeintroncatie: Het onbegrensde domein $\mathbb{R}^N$ wordt getronceerd tot een begrensde domein $D(R)$.

   • Bovengrenzen: Verkregen via Dirichlet-randvoorwaarden ($u=0$ op $\partial D(R)$) met behulp van standaard conformerende $P_1$ Galerkin-elementen.
   • Ondergrenzen: Verkregen via Neumann-randvoorwaarden ($\partial u/\partial n = 0$). Een sleutellemma stelt dat als het minimum van de externe potentiaal $\sigma_{ext}(D(R))$ de $k$-de eigenwaarde overschrijdt, de Neumann-eigenwaarde op $D(R)$ een geldige ondergrens biedt voor de eigenwaarde in de volledige ruimte. Deze voorwaarde wordt a posteriori geverifieerd met behulp van de berekende bovengrens.

2. CECR-raamwerk voor tekenwisselende potentialen:

   • De methode maakt gebruik van een paars-ruimte raamwerk waarbij functies worden opgeheven naar paren $(u, \Pi_{0,h}u)$, waardoor de kinetische energie (gradiënt) wordt gescheiden van de reactie-energie (potentiaal).
   • Om de singulariteit te hanteren zonder divergerende verschuivingen, introduceren de auteurs een elementgewijze reactieconstante $\Gamma_h := \max_K c_{h|K}^- \cdot h_K^2 / \pi^2$. In tegenstelling tot de globale $L^\infty$-verschuiving $\gamma_h$, verdwijnt $\Gamma_h$ ($O(h^2)$ of $O(h)$ afhankelijk van de gradatie) naarmate het rooster wordt verfijnd, zelfs voor Coulomb-singulariteiten.
   • Een vaste verschuiving $\sigma > 0$ (onafhankelijk van de roostergrootte $h$) wordt toegepast op de potentiaal om de coerciviteit van de verschoven bilineaire vorm te waarborgen. Dit voorkomt de noodzaak van rooster-afhankelijke verschuivingen die divergeren.

3. Perturbatieanalyse:

   • De continue potentiaal $V$ wordt benaderd door zijn stuksgewijs constante elementgemiddelde $c_h$.
   • Een perturbatielemma begrenst de fout die wordt geïntroduceerd door $V$ te vervangen door $c_h$. De perturbatieverhouding $\epsilon_h$ blijkt te convergeren naar nul ($O(h^2)$ op gegradueerde roosters) onder specifieke aannames over roosterverfijning in de buurt van de singulariteiten.

Belangrijkste bijdragen
Het artikel claimt de volgende specifieke bijdragen:

• Eerste expliciete convergente ondergrenzen: Het presenteert het eerste numerieke algoritme dat expliciete, convergente tweezijdige eigenwaardegrenzen biedt voor Schrödingervergelijkingen met onbegrensde singuliere potentialen op onbegrensde domeinen.
• Paars-ruimte abstracte stelling: Het ontwikkelt een paars-ruimte raamwerk dat de CECR-analyse opheft naar een tweecomponentenruimte, waardoor een zelfstandig max-min-bewijs mogelijk is voor de ondergrensstelling voor tekenonbepaalde potentialen.
• Elementgewijze reactieconstante ($\Gamma_h$): Het introduceert $\Gamma_h$ om de divergente $L^\infty$-verschuiving te vervangen. Dit stelt de CECR-gap in staat om te sluiten met een snelheid van $O(h^2)$ op gegradueerde roosters, waardoor de noodzaak wordt omzeild voor onberekenbare vorm-normconstanten.
• Vaste-verschuiving tweezijdige grens: Het assembleert een volledig berekenbare omhulling $L_k \le \lambda_k \le U_k$ die geen $L^\infty$-controle vereist van het negatieve deel van de potentiaal.
• Toepassing op $H_2^+$: Het specialiseert de theorie tot de tweegestelde Coulomb-potentiaal, waarbij expliciete formules worden geboden voor toelaatbare verschuivingen en perturbatieparameters in zowel 2D als 3D.

Numerieke resultaten
Het algoritme is gevalideerd op vier testgevallen:

1. 2D Waterstofatoom ($V = -1/|x|$): De methode omvat succesvol de exacte grondtoestand $\lambda_1 = -1$. Het gecertificeerde ondergrensinterval versmalt aanzienlijk bij roosterverfijning, waarbij de optimaal-verschuivingsgrenzen op het fijnste niveau een breedte van $0,002$ bereiken.
2. 2D $H_2^+$-ion: Resultaten voor de tweegestelde potentiaal tonen convergentie, waarbij het gecertificeerde interval voor de grondtoestand versmalt van $[-4,43, -1,317]$ naar $[-2,36, -1,317]$ naarmate het rooster verfijnt.
3. 3D Waterstofatoom: Op een bol-domein omvat de methode de exacte $\lambda_1 = -0,25$. De breedte van het gecertificeerde Neumann-omhulsel convergeert met $O(h^2)$.
4. 3D $H_2^+$-ion: Voor het moleculaire ion op een doos-domein biedt de methode een gecertificeerd interval $[-1,045, -0,544]$ (breedte $0,501$) voor de grondtoestand, dat de referentiewaarde $\approx -0,5513$ bevat.

In alle gevallen wordt waargenomen dat de CECR-ondergrensgap $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ sluit met de voorspelde $O(h^2)$-snelheid op gegradueerde roosters. Het artikel merkt op dat hoewel de 3D-gaps breder zijn dan in 2D vanwege de langzamere afname van de elementgewijze constante op grove tetraëdrische roosters, het convergentiegedrag consistent is met de theorie.

Betekenis en claims
Het artikel stelt dat dit werk het eerste finite element-raamwerk vormt dat in staat is expliciete en convergente ondergrenzen te bieden voor eigenwaarden geassocieerd met onbegrensde singuliere potentialen. De betekenis ligt in het overwinnen van de "coerciviteitsgap" veroorzaakt door aantrekkende singuliere potentialen, wat eerder ertoe leidde dat standaard ondergrenstechnieken divergeerden.

De auteurs benadrukken dat de methode expliciet is: alle constanten (interpolatieconstanten, perturbatiegrenzen, verschuivingsparameters) worden in gesloten vorm afgeleid uit de rooster- en potentiaaldata. Hoewel de huidige implementatie dubbelprecisie-aritmetiek gebruikt, stellen de auteurs expliciet dat de conversie naar volledig geverifieerde (IEEE-754-veilige, interval-rekenkundige) omhullingen een aparte verificatiestap is die voor toekomstig werk wordt achtergelaten. Het artikel merkt ook op dat hoewel de huidige grenzen op projectie gebaseerd zijn, ze dienen als een rigoureuze startpunt ($\rho$) voor scherpere nabewerkingstechnieken zoals Lehmann–Goerisch-grenzen, die de intervallen verder kunnen verstrakken.

Het werk claimt niet het multi-elektronenprobleem direct op te lossen, maar identificeert de uitbreiding naar magnetische Hamiltonianen en multi-elektronensystemen met antisymmetrie-beperkingen als natuurlijke toekomstige richtingen.