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

Resumen Técnico: Límites explícitos bilaterales de autovalores para operadores de Schrödinger con potenciales singulares mediante el Método de Elementos Finitos

Enunciado del Problema
El artículo aborda el cálculo de límites bilaterales explícitos y computables para los autovalores $\lambda_k$ del operador de Schrödinger $H = -\Delta + V$ en dominios no acotados $\mathbb{R}^N$ ($N=2, 3$). El desafío principal radica en manejar potenciales atractivos singulares de tipo Coulomb, específicamente $V(x) = -Z/|x|$ (que modela el átomo de hidrógeno) y potenciales de dos centros $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$ (que modelan el ion molecular $H_2^+$).

Los métodos numéricos existentes suelen proporcionar cotas superiores mediante el método de Galerkin conforme (Rayleigh–Ritz). Obtener cotas inferiores rigurosas es significativamente más difícil, particularmente para dominios no acotados y potenciales singulares donde el potencial no está acotado inferiormente. Las técnicas clásicas de cota inferior (por ejemplo, Lehmann–Maehly–Goerisch, Kato–Temple) a menudo requieren métodos no conformes o condiciones específicas que fallan cuando el potencial $V$ es negativo y singular. Además, los enfoques estándar de elementos finitos para cotas inferiores suelen depender de desplazar el potencial por su norma $L^\infty$ de la parte negativa ($\gamma_h = \|c_h^-\|_{L^\infty}$). Para singularidades de Coulomb, $\gamma_h$ diverge a medida que se refina la malla ($O(h^{-1})$), lo que provoca que la cota inferior no converja.

Metodología
El algoritmo propuesto combina la truncación del dominio con una extensión del método de elementos finitos Crouzeix–Raviart Enriquecido Compuesto (CECR). El enfoque se estructura en tres componentes principales:

1. Truncación del Dominio: El dominio no acotado $\mathbb{R}^N$ se trunca a un dominio acotado $D(R)$.

   • Cotas Superiores: Se obtienen mediante condiciones de contorno de Dirichlet ($u=0$ en $\partial D(R)$) utilizando elementos de Galerkin conformes $P_1$ estándar.
   • Cotas Inferiores: Se obtienen mediante condiciones de contorno de Neumann ($\partial u/\partial n = 0$). Un lema clave establece que si el mínimo del potencial exterior $\sigma_{ext}(D(R))$ excede el $k$-ésimo autovalor, el autovalor de Neumann en $D(R)$ proporciona una cota inferior válida para el autovalor en todo el espacio. Esta condición se verifica a posteriori utilizando la cota superior calculada.

2. Marco CECR para Potenciales que Cambian de Signo:

   • El método emplea un marco de espacio de pares donde las funciones se elevan a pares $(u, \Pi_{0,h}u)$, separando la energía cinética (gradiente) de la energía de reacción (potencial).
   • Para manejar la singularidad sin desplazamientos divergentes, los autores introducen una constante de reacción por elemento $\Gamma_h := \max_K c_{h|K}^- \cdot h_K^2 / \pi^2$. A diferencia del desplazamiento global $L^\infty$ $\gamma_h$, $\Gamma_h$ se anula ($O(h^2)$ u $O(h)$ dependiendo de la graduación) a medida que se refina la malla, incluso para singularidades de Coulomb.
   • Se aplica un desplazamiento fijo $\sigma > 0$ (independiente del tamaño de malla $h$) al potencial para asegurar la coercividad de la forma bilineal desplazada. Esto evita la necesidad de desplazamientos dependientes de la malla que divergen.

3. Análisis de Perturbación:

   • El potencial continuo $V$ se aproxima por su promedio elemental constante por partes $c_h$.
   • Un lema de perturbación acota el error introducido al reemplazar $V$ por $c_h$. Se demuestra que la razón de perturbación $\epsilon_h$ converge a cero ($O(h^2)$ en mallas graduadas) bajo suposiciones específicas de graduación de malla cerca de las singularidades.

Contribuciones Clave
El artículo afirma las siguientes contribuciones específicas:

• Primeras Cotas Inferiores Convergentes Explícitas: Presenta el primer algoritmo numérico que proporciona límites bilaterales de autovalores explícitos y convergentes para operadores de Schrödinger con potenciales singulares no acotados en dominios no acotados.
• Teorema Abstracto de Espacio de Pares: Desarrolla un marco de espacio de pares que eleva el análisis CECR a un espacio de dos componentes, permitiendo una demostración autocontenida de tipo max-min del teorema de cota inferior para potenciales indefinidos en signo.
• Constante de Reacción por Elemento ($\Gamma_h$): Introduce $\Gamma_h$ para reemplazar el desplazamiento divergente $L^\infty$. Esto permite que la brecha CECR se cierre a una tasa de $O(h^2)$ en mallas graduadas, evitando la necesidad de constantes de norma de forma no computables.
• Cota Bilateral con Desplazamiento Fijo: Ensambla un recinto completamente computable $L_k \le \lambda_k \le U_k$ que no requiere control $L^\infty$ de la parte negativa del potencial.
• Aplicación a $H_2^+$: Especializa la teoría al potencial de Coulomb de dos centros, proporcionando fórmulas explícitas para desplazamientos admisibles y parámetros de perturbación tanto en 2D como en 3D.

Resultados Numéricos
El algoritmo se valida en cuatro casos de prueba:

1. Átomo de Hidrógeno 2D ($V = -1/|x|$): El método acota con éxito el estado fundamental exacto $\lambda_1 = -1$. El intervalo de cota inferior certificado se estrecha significativamente con el refinamiento de la malla, alcanzando anchuras de $0.002$ en el nivel más fino con los límites de desplazamiento óptimo.
2. Ion $H_2^+$ 2D: Los resultados para el potencial de dos centros muestran convergencia, con el intervalo certificado para el estado fundamental estrechándose de $[-4.43, -1.317]$ a $[-2.36, -1.317]$ a medida que la malla se refina.
3. Átomo de Hidrógeno 3D: En un dominio esférico, el método acota el $\lambda_1 = -0.25$ exacto. La anchura del recinto certificado de Neumann converge a $O(h^2)$.
4. Ion $H_2^+$ 3D: Para el ion molecular en un dominio tipo caja, el método proporciona un intervalo certificado $[-1.045, -0.544]$ (anchura $0.501$) para el estado fundamental, que contiene el valor de referencia $\approx -0.5513$.

En todos los casos, se observa que la brecha de la cota inferior CECR $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ se cierra a la tasa predicha de $O(h^2)$ en mallas graduadas. El artículo señala que, aunque las brechas en 3D son más amplias que en 2D debido a la descomposición más lenta de la constante por elemento en mallas tetraédricas gruesas, el comportamiento de convergencia es consistente con la teoría.

Significado y Afirmaciones
El artículo afirma que este trabajo constituye el primer marco de elementos finitos capaz de proporcionar cotas inferiores explícitas y convergentes para autovalores asociados a potenciales singulares no acotados. El significado radica en superar la "brecha de coercividad" causada por potenciales atractivos singulares, lo que anteriormente hacía que las técnicas estándar de cota inferior divergieran.

Los autores enfatizan que el método es explícito: todas las constantes (constantes de interpolación, cotas de perturbación, parámetros de desplazamiento) se derivan en forma cerrada a partir de los datos de la malla y del potencial. Si bien la implementación actual utiliza aritmética de doble precisión, los autores declaran explícitamente que la conversión a recintos totalmente verificados (seguros IEEE-754, aritmética de intervalos) es un paso de verificación separado dejado para trabajos futuros. El artículo también señala que, aunque las cotas actuales se basan en proyecciones, sirven como un punto de partida riguroso ($\rho$) para técnicas de post-procesamiento más precisas como los límites de Lehmann–Goerisch, que podrían estrechar aún más los intervalos.

El trabajo no afirma resolver directamente el problema de múltiples electrones, pero identifica la extensión a hamiltonianos magnéticos y sistemas de múltiples electrones con restricciones de antisimetría como direcciones naturales futuras.