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

Technische Zusammenfassung: Explizite zweiseitige Eigenwertabschätzungen für Schrödinger-Operatoren mit singulären Potentialen mittels der Finite-Elemente-Methode

Problemstellung
Der Artikel behandelt die Berechnung expliziter, berechenbarer zweiseitiger Abschätzungen für die Eigenwerte $\lambda_k$ des Schrödinger-Operators $H = -\Delta + V$ auf unbeschränkten Gebieten $\mathbb{R}^N$ ($N=2, 3$). Die Haupt Herausforderung besteht im Umgang mit singulären anziehenden Potentialen vom Coulomb-Typ, insbesondere $V(x) = -Z/|x|$ (Modellierung des Wasserstoffatoms) und zweizentrierten Potentialen $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$ (Modellierung des $H_2^+$-Molekülions).

Bestehende numerische Methoden liefern typischerweise obere Schranken über die konforme Galerkin-Methode (Rayleigh–Ritz). Die Gewinnung rigoroser unterer Schranken ist erheblich schwieriger, insbesondere für unbeschränkte Gebiete und singuläre Potentiale, bei denen das Potential nach unten unbeschränkt ist. Klassische Techniken für untere Schranken (z. B. Lehmann–Maehly–Goerisch, Kato–Temple) erfordern oft nicht-konforme Methoden oder spezifische Bedingungen, die versagen, wenn das Potential $V$ negativ und singulär ist. Darüber hinaus stützen sich Standard-Finite-Elemente-Ansätze für untere Schranken häufig auf eine Verschiebung des Potentials um die $L^\infty$-Norm des negativen Teils ($\gamma_h = \|c_h^-\|_{L^\infty}$). Für Coulomb-Singularitäten divergiert $\gamma_h$ bei Verfeinerung des Gitters ($O(h^{-1})$), wodurch die untere Schranke versagt, zu konvergieren.

Methodik
Der vorgeschlagene Algorithmus kombiniert Gebietsabschneidung mit einer Erweiterung der Composite Enriched Crouzeix–Raviart (CECR)-Finite-Elemente-Methode. Der Ansatz ist in drei Hauptkomponenten strukturiert:

1. Gebietsabschneidung: Das unbeschränkte Gebiet $\mathbb{R}^N$ wird auf ein beschränktes Gebiet $D(R)$ abgeschnitten.

   • Obere Schranken: Werden über Dirichlet-Randbedingungen ($u=0$ auf $\partial D(R)$) unter Verwendung standardkonformer $P_1$-Galerkin-Elemente ermittelt.
   • Untere Schranken: Werden über Neumann-Randbedingungen ($\partial u/\partial n = 0$) ermittelt. Ein zentrales Lemma stellt fest, dass, wenn das Minimum des äußeren Potentials $\sigma_{ext}(D(R))$ den $k$-ten Eigenwert übersteigt, der Neumann-Eigenwert auf $D(R)$ eine gültige untere Schranke für den Eigenwert im Vollraum liefert. Diese Bedingung wird a posteriori unter Verwendung der berechneten oberen Schranke überprüft.

2. CECR-Rahmenwerk für Vorzeichen-wechselnde Potentiale:

   • Die Methode verwendet ein Paarraum-Rahmenwerk, bei dem Funktionen zu Paaren $(u, \Pi_{0,h}u)$ gehoben werden, wodurch die kinetische Energie (Gradient) von der Reaktionsenergie (Potential) getrennt wird.
   • Um die Singularität ohne divergierende Verschiebungen zu handhaben, führen die Autoren eine elementweise Reaktionskonstante $\Gamma_h := \max_K c_{h|K}^- \cdot h_K^2 / \pi^2$ ein. Im Gegensatz zur globalen $L^\infty$-Verschiebung $\gamma_h$ geht $\Gamma_h$ bei Verfeinerung des Gitters gegen Null ($O(h^2)$ oder $O(h)$, je nach Vernetzung), selbst für Coulomb-Singularitäten.
   • Eine feste Verschiebung $\sigma > 0$ (unabhängig von der Gittergröße $h$) wird auf das Potential angewendet, um die Koerzitivität der verschobenen Bilinearform sicherzustellen. Dies vermeidet die Notwendigkeit von gitterabhängigen Verschiebungen, die divergieren.

3. Störungsanalyse:

   • Das kontinuierliche Potential $V$ wird durch seinen stückweise konstanten Elementmittelwert $c_h$ approximiert.
   • Ein Störungslemma schätzt den Fehler ab, der durch den Ersatz von $V$ durch $c_h$ entsteht. Das Störungsverhältnis $\epsilon_h$ konvergiert unter spezifischen Annahmen zur Gitterverfeinerung in der Nähe der Singularitäten gegen Null ($O(h^2)$ auf verfeinerten Gittern).

Hauptbeiträge
Der Artikel beansprucht folgende spezifische Beiträge:

• Erste explizite konvergente untere Schranken: Er stellt den ersten numerischen Algorithmus vor, der explizite, konvergente zweiseitige Eigenwertabschätzungen für Schrödinger-Operatoren mit unbeschränkten singulären Potentialen auf unbeschränkten Gebieten liefert.
• Paarraum-abstraktes Theorem: Er entwickelt ein Paarraum-Rahmenwerk, das die CECR-Analyse auf einen zweikomponentigen Raum hebt und einen in sich geschlossenen Max-Min-Beweis für das Theorem der unteren Schranke für vorzeichenindefinite Potentiale ermöglicht.
• Elementweise Reaktionskonstante ($\Gamma_h$): Sie führt $\Gamma_h$ ein, um die divergente $L^\infty$-Verschiebung zu ersetzen. Dies ermöglicht, dass die CECR-Lücke auf verfeinerten Gittern mit einer Rate von $O(h^2)$ schließt und die Notwendigkeit nicht berechenbarer Form-Norm-Konstanten umgeht.
• Zweiseitige Schranke mit fester Verschiebung: Sie stellt eine vollständig berechenbare Einschließung $L_k \le \lambda_k \le U_k$ zusammen, die keine $L^\infty$-Kontrolle des negativen Teils des Potentials erfordert.
• Anwendung auf $H_2^+$: Sie spezialisiert die Theorie auf das zweizentrierte Coulomb-Potential und liefert explizite Formeln für zulässige Verschiebungen und Störungsparameter sowohl in 2D als auch in 3D.

Numerische Ergebnisse
Der Algorithmus wird an vier Testfällen validiert:

1. 2D-Wasserstoffatom ($V = -1/|x|$): Die Methode umschließt erfolgreich den exakten Grundzustand $\lambda_1 = -1$. Das zertifizierte Intervall der unteren Schranke verengt sich bei Gitterverfeinerung erheblich, wobei die Grenzen mit optimaler Verschiebung auf dem feinsten Niveau eine Breite von $0,002$ erreichen.
2. 2D-$H_2^+$-Ion: Ergebnisse für das zweizentrierte Potential zeigen Konvergenz, wobei das zertifizierte Intervall für den Grundzustand sich von $[-4,43, -1,317]$ auf $[-2,36, -1,317]$ verengt, während das Gitter verfeinert wird.
3. 3D-Wasserstoffatom: Auf einem Kugelgebiet umschließt die Methode den exakten Wert $\lambda_1 = -0,25$. Die Breite der zertifizierten Neumann-Einschließung konvergiert mit $O(h^2)$.
4. 3D-$H_2^+$-Ion: Für das Molekülion auf einem Boxgebiet liefert die Methode ein zertifiziertes Intervall $[-1,045, -0,544]$ (Breite $0,501$) für den Grundzustand, das den Referenzwert $\approx -0,5513$ enthält.

In allen Fällen wird beobachtet, dass sich die CECR-Lücke der unteren Schranke $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ auf verfeinerten Gittern mit der vorhergesagten Rate $O(h^2)$ schließt. Der Artikel stellt fest, dass, obwohl die 3D-Lücken aufgrund des langsameren Abfalls der elementweisen Konstante auf groben Tetraedergittern breiter sind als in 2D, das Konvergenzverhalten mit der Theorie übereinstimmt.

Bedeutung und Behauptungen
Der Artikel behauptet, dass diese Arbeit das erste Finite-Elemente-Rahmenwerk darstellt, das in der Lage ist, explizite und konvergente untere Schranken für Eigenwerte im Zusammenhang mit unbeschränkten singulären Potentialen zu liefern. Die Bedeutung liegt in der Überwindung der durch anziehende singuläre Potentiale verursachten „Koerzitivitätslücke", die zuvor dazu führte, dass Standardtechniken für untere Schranken divergierten.

Die Autoren betonen, dass die Methode explizit ist: Alle Konstanten (Interpolationskonstanten, Störungsschranken, Verschiebungsparameter) werden in geschlossener Form aus den Gitter- und Potentialdaten abgeleitet. Obwohl die aktuelle Implementierung Gleitkommazahlen mit doppelter Genauigkeit verwendet, geben die Autoren ausdrücklich an, dass die Umwandlung in vollständig verifizierte (IEEE-754-sichere, Intervallarithmetik) Einschließungen ein separater Verifikationsschritt ist, der zukünftigen Arbeiten vorbehalten bleibt. Der Artikel stellt auch fest, dass, obwohl die aktuellen Schranken projektionsbasiert sind, sie als rigoroser Ausgangspunkt ($\rho$) für schärfere Nachverarbeitungstechniken wie Lehmann–Goerisch-Schranken dienen, die die Intervalle weiter verengen könnten.

Die Arbeit beansprucht nicht, das Mehr-Elektronen-Problem direkt zu lösen, identifiziert aber die Erweiterung auf magnetische Hamilton-Operatoren und Mehr-Elektronen-Systeme mit Antisymmetrie-Bedingungen als natürliche zukünftige Richtungen.