Guaranteed Lower Eigenvalue Bounds for Spectral Galerkin Methods with Application to Schrödinger Operators

Technische Samenvatting: Gegarandeerde ondergrenzen voor eigenwaarden van spectrale Galerkin-methoden met toepassing op Schrödinger-operatoren

Probleemstelling
Spectrale Galerkin-methoden staan algemeen bekend om hun hoge precisie bij de benadering van eigenwaarden, in het bijzonder voor gladde problemen waar zij spectrale (super-algebraïsche) convergentie vertonen. Echter, een rigoureuze, berekenbare ondergrens voor deze eigenwaarden is binnen het spectrale kader een hardnekkig probleem gebleven. Hoewel klassieke insluitingsmethoden (bijv. Kato, Weinstein–Temple, Lehmann–Goerisch) ondergrenzen kunnen bieden, vereisen zij doorgaans a priori informatie over naburige eigenwaarden of specifieke probleemstructuren. Daartegenover staan bestaande kaders voor gegarandeerde ondergrenzen, zoals die gebaseerd op eindige-elementenmethoden (FEM), die vertrouwen op lokale roostergeometrie en interpolatiefoutschattingen (constanten $C_h$) die geen direct spectraal tegenhanger hebben. Dit artikel adresseert de kloof door een rigoureus, onvoorwaardelijk kader voor ondergrenzen te vestigen voor conforme spectrale Galerkin-methoden toegepast op Schrödinger-operatoren ($-\Delta + V$).

Methodologie en Kader
Het artikel bouwt voort op een abstract Hilbert-ruimte kader (oorspronkelijk ontwikkeld door Liu [16, 17]) dat een discrete Galerkin-eigenwaarde $\lambda_{k,N}$ relateert aan de exacte eigenwaarde $\lambda_k$ via een projectiefoutconstante $C_N$. De kernongelijkheid is:
$$\lambda_k \ge \frac{\lambda_{k,N}}{1 + C_N^2 \lambda_{k,N}}$$
De methodologie verloopt via vier verschillende mechanismen om $C_N$ voor spectrale ruimtes te bepalen:

1. Exacte eigenbasisruimtes: Voor proefruimtes gespannen door exacte eigenfuncties van de operator, leidt het artikel een gesloten vorm, optimale projectieconstante af: $C_N = \lambda_{M+1}^{-1/2}$, waarbij $\lambda_{M+1}$ de eerste weggelaten eigenwaarde is. Dit wordt vastgesteld via een Parseval-argument, onafhankelijk van de roostergeometrie.
2. Projectiekloof-mechanisme (Standaard matrix): Voor Schrödinger-operatoren met niet-negatieve potentialen ($V \ge 0$) waarbij de proefruimte (bijv. trigonometrische polynomen) geen eigenbasis is, verschilt de standaard Ritz-projectie $P_N$ van de Laplaciaanse modale projectie $P_N^0$. Het artikel introduceert een "projectiekloof"-analyse, waarbij het verschil $e = P_N^0 u - P_N u$ wordt gecontroleerd met behulp van een hulp-discrete operator. Dit levert een constante $C_N^{(V)}$ op die afhankelijk is van een potentiaalvormfactor $\eta_V$. Voor begrensde potentialen hangt deze factor initieel af van $\|V\|_{L^\infty}$, wat leidt tot een constante $C_N^{(V)} = \sqrt{(1 + \min(X^2, X))/\nu^*}$, waarbij $X = \|V\|_{L^\infty}/\lambda_{1,N}$ en $\nu^*$ de spectrale kloof van de Laplaciaan is.
3. Composiet Discretisatie (Verwijderen van $\|V\|_{L^\infty}$-afhankelijkheid): Om de inflatie van de constante veroorzaakt door grote potentialen (veelvoorkomend bij domeintruncatie voor onbegrensde domeinen) te elimineren, stelt het artikel een composiet discretisatie voor. Dit omvat:
   • Het construeren van een bandgelimiteerde onderbenadering $W$ van de potentiaal $V$.
   • Het definiëren van een gereduceerde-graad gewogen projectie $\Pi_m^{\tilde{W}}$ op een subruimte $V_m$.
   • Het gebruik van een "slaving identity" (een exacte aliasing-eigenschap) waarbij de gewogen projectie de interpolatierest van de hogere-graads ruimte annihileert.
     Deze constructie herstelt de scherpe, potentiaal-onafhankelijke constante $C_N = \nu^{*-1/2}$, waardoor de kwaliteit van de grens effectief ontkoppeld wordt van de grootte van de potentiaal.
4. Domeintruncatie: Voor operatoren op $\mathbb{R}^d$ combineert het artikel de spectrale ondergrenzen op een getrunceerde box $\Omega$ met Neumann-randvoorwaarden en a posteriori confinement-controles (het vergelijken van de berekende bovengrens met het minimum van de externe potentiaal $\sigma(\Omega)$) om gecertificeerde twee-zijdige grenzen voor de gehele ruimte te produceren.

Belangrijkste Bijdragen

• Optimale Spectrale Constante: De afleiding van de gesloten vorm constante $C_N = \lambda_{M+1}^{-1/2}$ voor eigenbasis proefruimtes, die als theoretische benchmark dient.
• Projectiekloof-analyse: De uitbreiding van de hulp-operator techniek naar spectrale methoden, waarbij expliciete constanten worden geleverd voor de standaard Galerkin-matrix van Schrödinger-operatoren met begrensde potentialen.
• Composiet Spectrale Methode: De ontwikkeling van een composiet discretisatie die de optimale spectrale-kloofconstante bereikt zonder afhankelijkheid van $\|V\|_{L^\infty}$, waarmee de beperkingen van de standaard matrixaanpak voor confinerende potentialen worden overwonnen.
• Uitbreiding naar Singuliere Potentialen: Een vooruitblik op hoe het hulp-projector mechanisme zich uitbreidt naar singuliere potentialen (bijv. Coulomb-singulariteiten) via een gelokaliseerde Hardy-ongelijkheid, waardoor de vereiste van begrensde $L^\infty$ normen wordt weggenomen.
• Twee-zijdige Enclosures: Het bieden van volledig berekenbare, gecertificeerde twee-zijdige grenzen voor Schrödinger-operatoren op $\mathbb{R}^d$ met behulp van Neumann domeintruncatie.

Numerieke Resultaten
Het artikel valideert de theorie op twee benchmark-potentialen in twee dimensies: een ring-potentiaal $V_1 = (|x|^2-1)^2$ en een anisotrope dubbele put $V_2 = (x_1^2-1)^2 + x_2^2$.

• Efficiëntie: De spectrale methode bereikt gecertificeerde grenzen die vergelijkbaar zijn met of superieur zijn aan die van gecertificeerde Crouzeix–Raviart verrijkte (CECR) eindige-elementenmethoden (uit [18]), met ongeveer twee ordes van grootte minder vrijheidsgraden (bijv. ~2.400 VGD's vs. ~400.000 VGD's).
• Nauwkeurigheid: Voor de ring-potentiaal ligt de spectrale enclosure bij $N=48$ strikt binnen de FEM-enclosure. Voor de dubbele put zijn de spectrale grenzen vergelijkbaar met de FEM-grenzen bij $N=48$.
• Mechanisme Vergelijking: In regimes met gematigde potentialen is de standaard hulp-projector methode (Theorema 4.9) efficiënt en bijna optimaal. Echter, voor grote confinerende potentialen (waar $X \gg 1$) zwelt de standaard kloof aanzienlijk op (bijv. 72% kloof vs. 0,48% voor de composiet methode). De composiet methode behoudt de optimale convergentieorde $O(N^{-2})$ ongeacht de grootte van de potentiaal.
• Tail Refinement: Het artikel demonstreert dat een "tail-refined" vormfactor (Theorema 4.5), berekenbaar vanuit de standaard matrix, de grenzen van de standaard methode aanzienlijk verbetert (de kloof verminderen van 72% naar 7,7% voor $V_1$), hoewel deze inferieur blijft aan de composiet methode voor zeer grote potentialen.

Betekenis en Claims
Het artikel claimt het langlopende probleem op te lossen van het verkrijgen van rigoureuze ondergrenzen direct uit spectrale discretisaties zonder a priori spectrale informatie. Door het projectie-gebaseerde kader uit te breiden naar spectrale methoden, biedt het:

1. Onvoorwaardelijke Grenzen: Grenzen die geen kennis vereisen van naburige eigenwaarden.
2. Optimale Constanten: Een constante die bepaald wordt door de spectrale kloof van de Laplaciaan in plaats van mesh-afhankelijke interpolatieconstanten.
3. Hoge Efficiëntie: Het vermogen om gecertificeerde nauwkeurigheid te bereiken met aanzienlijk minder vrijheidsgraden dan lage-orde FEM, waarbij gebruik wordt gemaakt van de spectrale convergentie van de bovengrenzen.
4. Robuustheid: De composiet methode zorgt ervoor dat de kwaliteit van de certificering niet verslechtert voor grote of confinerende potentialen, een cruciale vereiste voor kwantumchemische toepassingen waarbij domeintruncatie wordt gebruikt.

De auteurs merken op dat terwijl de bovengrenzen spectraal convergeren, de gecertificeerde ondergrenzen een algebraïsche snelheid van $O(N^{-2})$ bereiken. Het overbrengen van exponentiële nauwkeurigheid naar de ondergrenzen blijft een open probleem. De uitbreiding naar singuliere potentialen (bijv. Coulomb-systemen) wordt erkend als een onderwerp voor een begeleidend artikel.