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

Sintesi Tecnica: Limiti Espliciti Bifacciali per gli Autovalori degli Operatori di Schrödinger con Potenziali Singolari tramite il Metodo degli Elementi Finiti

Enunciato del Problema
Il lavoro affronta il calcolo di limiti bifacciali espliciti e calcolabili per gli autovalori $\lambda_k$ dell'operatore di Schrödinger $H = -\Delta + V$ su domini illimitati $\mathbb{R}^N$ ($N=2, 3$). La sfida principale risiede nella gestione di potenziali attrattivi singolari di tipo Coulombiano, specificamente $V(x) = -Z/|x|$ (che modella l'atomo di idrogeno) e potenziali a due centri $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$ (che modellano lo ione molecolare $H_2^+$).

I metodi numerici esistenti forniscono tipicamente limiti superiori tramite il metodo di Galerkin conforme (Rayleigh–Ritz). Ottenere limiti inferiori rigorosi è significativamente più difficile, in particolare per domini illimitati e potenziali singolari in cui il potenziale è illimitato dal basso. Le tecniche classiche per i limiti inferiori (ad es. Lehmann–Maehly–Goerisch, Kato–Temple) richiedono spesso metodi non conformi o condizioni specifiche che falliscono quando il potenziale $V$ è negativo e singolare. Inoltre, gli approcci standard agli elementi finiti per i limiti inferiori si basano spesso sullo spostamento del potenziale mediante la norma $L^\infty$ della sua parte negativa ($\gamma_h = \|c_h^-\|_{L^\infty}$). Per le singolarità Coulombiane, $\gamma_h$ diverge al raffinarsi della mesh ($O(h^{-1})$), causando il fallimento della convergenza del limite inferiore.

Metodologia
L'algoritmo proposto combina la troncatura del dominio con un'estensione del metodo agli elementi finiti Crouzeix–Raviart Composito Arricchito (CECR). L'approccio è strutturato in tre componenti principali:

1. Troncatura del Dominio: Il dominio illimitato $\mathbb{R}^N$ viene troncato in un dominio limitato $D(R)$.

   • Limiti Superiori: Ottenuti tramite condizioni al contorno di Dirichlet ($u=0$ su $\partial D(R)$) utilizzando elementi conformi standard $P_1$ di Galerkin.
   • Limiti Inferiori: Ottenuti tramite condizioni al contorno di Neumann ($\partial u/\partial n = 0$). Un lemma chiave stabilisce che se il minimo del potenziale esterno $\sigma_{ext}(D(R))$ supera il $k$-esimo autovalore, l'autovalore di Neumann su $D(R)$ fornisce un limite inferiore valido per l'autovalore dello spazio completo. Questa condizione è verificata a posteriori utilizzando il limite superiore calcolato.

2. Quadro CECR per Potenziali a Segno Variabile:

   • Il metodo impiega un quadro a spazio di coppie in cui le funzioni sono sollevate in coppie $(u, \Pi_{0,h}u)$, separando l'energia cinetica (gradiente) dall'energia di reazione (potenziale).
   • Per gestire la singolarità senza spostamenti divergenti, gli autori introducono una costante di reazione elementare $\Gamma_h := \max_K c_{h|K}^- \cdot h_K^2 / \pi^2$. A differenza dello spostamento globale $L^\infty$ $\gamma_h$, $\Gamma_h$ si annulla ($O(h^2)$ o $O(h)$ a seconda della griglia) al raffinarsi della mesh, anche per singolarità Coulombiane.
   • Viene applicato uno spostamento fisso $\sigma > 0$ (indipendente dalla dimensione della mesh $h$) al potenziale per garantire la coercività della forma bilineare spostata. Ciò evita la necessità di spostamenti dipendenti dalla mesh che divergono.

3. Analisi della Perturbazione:

   • Il potenziale continuo $V$ è approssimato dalla sua media elementare costante a tratti $c_h$.
   • Un lemma di perturbazione delimita l'errore introdotto sostituendo $V$ con $c_h$. Il rapporto di perturbazione $\epsilon_h$ è mostrato convergere a zero ($O(h^2)$ su mesh grigliate) sotto specifiche ipotesi di grigliatura vicino alle singolarità.

Contributi Chiave
Il lavoro rivendica i seguenti contributi specifici:

• Primi Limiti Inferiori Convergenti Espliciti: Presenta il primo algoritmo numerico che fornisce limiti bifacciali espliciti e convergenti per gli autovalori degli operatori di Schrödinger con potenziali singolari illimitati su domini illimitati.
• Teorema Astratto a Spazio di Coppie: Sviluppa un quadro a spazio di coppie che solleva l'analisi CECR a uno spazio a due componenti, permettendo una prova autocontenuta max-min del teorema del limite inferiore per potenziali a segno indefinito.
• Costante di Reazione Elementare ($\Gamma_h$): Introduce $\Gamma_h$ per sostituire lo spostamento divergente $L^\infty$. Ciò permette al gap CECR di chiudersi a un tasso di $O(h^2)$ su mesh grigliate, aggirando la necessità di costanti di norma della forma non calcolabili.
• Limite Bifacciale a Spostamento Fisso: Assembla un'inclusione completamente calcolabile $L_k \le \lambda_k \le U_k$ che non richiede il controllo $L^\infty$ della parte negativa del potenziale.
• Applicazione a $H_2^+$: Specializza la teoria per il potenziale Coulombiano a due centri, fornendo formule esplicite per gli spostamenti ammissibili e i parametri di perturbazione sia in 2D che in 3D.

Risultati Numerici
L'algoritmo è validato su quattro casi di prova:

1. Atomo di Idrogeno 2D ($V = -1/|x|$): Il metodo inquadra con successo lo stato fondamentale esatto $\lambda_1 = -1$. L'intervallo inferiore certificato si restringe significativamente con il raffinarsi della mesh, con i limiti a spostamento ottimale che raggiungono larghezze di $0.002$ al livello più fine.
2. Ione $H_2^+$ 2D: I risultati per il potenziale a due centri mostrano convergenza, con l'intervallo certificato per lo stato fondamentale che si restringe da $[-4.43, -1.317]$ a $[-2.36, -1.317]$ mentre la mesh si raffina.
3. Atomo di Idrogeno 3D: Su un dominio a sfera, il metodo inquadra l'esatto $\lambda_1 = -0.25$. La larghezza dell'inclusione di Neumann certificata converge a $O(h^2)$.
4. Ione $H_2^+$ 3D: Per lo ione molecolare su un dominio a scatola, il metodo fornisce un intervallo certificato $[-1.045, -0.544]$ (larghezza $0.501$) per lo stato fondamentale, che contiene il valore di riferimento $\approx -0.5513$.

In tutti i casi, il gap del limite inferiore CECR $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ è osservato chiudersi al tasso previsto $O(h^2)$ su mesh grigliate. Il lavoro nota che, sebbene i gap 3D siano più ampi rispetto al 2D a causa del decadimento più lento della costante elementare su mesh tetraedriche grezze, il comportamento di convergenza è coerente con la teoria.

Significato e Rivendicazioni
Il lavoro afferma che questo studio costituisce il primo quadro agli elementi finiti capace di fornire limiti inferiori espliciti e convergenti per gli autovalori associati a potenziali singolari illimitati. Il significato risiede nel superare il "gap di coercività" causato dai potenziali attrattivi singolari, che in precedenza causavano la divergenza delle tecniche standard per i limiti inferiori.

Gli autori sottolineano che il metodo è esplicito: tutte le costanti (costanti di interpolazione, limiti di perturbazione, parametri di spostamento) sono derivate in forma chiusa dai dati della mesh e del potenziale. Sebbene l'implementazione attuale utilizzi l'aritmetica in doppia precisione, gli autori dichiarano esplicitamente che la conversione in inclusioni completamente verificate (sicure IEEE-754, ad aritmetica intervallare) è un passo di verifica separato lasciato per lavori futuri. Il lavoro nota inoltre che, sebbene i limiti attuali siano basati sulla proiezione, servono come punto di partenza rigoroso ($\rho$) per tecniche di post-processing più precise come i limiti di Lehmann–Goerisch, che potrebbero ulteriormente stringere gli intervalli.

Il lavoro non rivendica di risolvere direttamente il problema multi-elettronico, ma identifica l'estensione agli Hamiltoniani magnetici e ai sistemi multi-elettronici con vincoli di antisimmetria come direzioni future naturali.