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

Resumo Técnico: Limites Explícitos Bilaterais para Autovalores de Operadores de Schrödinger com Potenciais Singulares via Método de Elementos Finitos

Enunciação do Problema
O artigo aborda o cálculo de limites bilaterais explícitos e computáveis para os autovalores $\lambda_k$ do operador de Schrödinger $H = -\Delta + V$ em domínios ilimitados $\mathbb{R}^N$ ($N=2, 3$). O principal desafio reside no manuseio de potenciais atrativos singulares do tipo Coulomb, especificamente $V(x) = -Z/|x|$ (modelando o átomo de hidrogênio) e potenciais de dois centros $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$ (modelando o íon molecular $H_2^+$).

Métodos numéricos existentes tipicamente fornecem limites superiores via o método de Galerkin conformante (Rayleigh–Ritz). Obter limites inferiores rigorosos é significativamente mais difícil, particularmente para domínios ilimitados e potenciais singulares onde o potencial é ilimitado inferiormente. Técnicas clássicas de limite inferior (por exemplo, Lehmann–Maehly–Goerisch, Kato–Temple) frequentemente exigem métodos não conformes ou condições específicas que falham quando o potencial $V$ é negativo e singular. Além disso, abordagens padrão de elementos finitos para limites inferiores frequentemente dependem do deslocamento do potencial por sua norma $L^\infty$ da parte negativa ($\gamma_h = \|c_h^-\|_{L^\infty}$). Para singularidades de Coulomb, $\gamma_h$ diverge conforme a malha é refinada ($O(h^{-1})$), fazendo com que o limite inferior falhe em convergir.

Metodologia
O algoritmo proposto combina truncamento de domínio com uma extensão do método de elementos finitos Crouzeix–Raviart Enriquecido Composto (CECR). A abordagem é estruturada em três componentes principais:

1. Truncamento de Domínio: O domínio ilimitado $\mathbb{R}^N$ é truncado para um domínio limitado $D(R)$.

   • Limites Superiores: Obtidos via condições de contorno de Dirichlet ($u=0$ em $\partial D(R)$) usando elementos de Galerkin conformes $P_1$ padrão.
   • Limites Inferiores: Obtidos via condições de contorno de Neumann ($\partial u/\partial n = 0$). Um lema chave estabelece que, se o mínimo do potencial exterior $\sigma_{ext}(D(R))$ exceder o $k$-ésimo autovalor, o autovalor de Neumann em $D(R)$ fornece um limite inferior válido para o autovalor do espaço completo. Esta condição é verificada a posteriori usando o limite superior computado.

2. Estrutura CECR para Potenciais que Mudam de Sinal:

   • O método emprega uma estrutura de espaço de pares onde as funções são elevadas a pares $(u, \Pi_{0,h}u)$, separando a energia cinética (gradiente) da energia de reação (potencial).
   • Para lidar com a singularidade sem deslocamentos divergentes, os autores introduzem uma constante de reação por elemento $\Gamma_h := \max_K c_{h|K}^- \cdot h_K^2 / \pi^2$. Diferentemente do deslocamento global $L^\infty$ $\gamma_h$, $\Gamma_h$ anula-se ($O(h^2)$ ou $O(h)$ dependendo do gradiente) conforme a malha é refinada, mesmo para singularidades de Coulomb.
   • Um deslocamento fixo $\sigma > 0$ (independente do tamanho da malha $h$) é aplicado ao potencial para garantir a coercividade da forma bilinear deslocada. Isso evita a necessidade de deslocamentos dependentes da malha que divergem.

3. Análise de Perturbação:

   • O potencial contínuo $V$ é aproximado por sua média elementar constante por partes $c_h$.
   • Um lema de perturbação limita o erro introduzido pela substituição de $V$ por $c_h$. A razão de perturbação $\epsilon_h$ é mostrada convergir para zero ($O(h^2)$ em malhas graduadas) sob suposições específicas de gradiente de malha próximas às singularidades.

Principais Contribuições
O artigo reivindica as seguintes contribuições específicas:

• Primeiros Limites Inferiores Convergentes Explícitos: Apresenta o primeiro algoritmo numérico que fornece limites bilaterais explícitos e convergentes para autovalores de operadores de Schrödinger com potenciais singulares ilimitados em domínios ilimitados.
• Teorema Abstrato de Espaço de Pares: Desenvolve uma estrutura de espaço de pares que eleva a análise CECR a um espaço de dois componentes, permitindo uma prova autocontida de máximo-mínimo do teorema do limite inferior para potenciais de sinal indefinido.
• Constante de Reação por Elemento ($\Gamma_h$): Introduz $\Gamma_h$ para substituir o deslocamento divergente $L^\infty$. Isso permite que a lacuna CECR feche a uma taxa de $O(h^2)$ em malhas graduadas, contornando a necessidade de constantes de norma de forma não computáveis.
• Limite Bilateral com Deslocamento Fixo: Monta um envoltório totalmente computável $L_k \le \lambda_k \le U_k$ que não requer controle $L^\infty$ da parte negativa do potencial.
• Aplicação a $H_2^+$: Especializa a teoria para o potencial de Coulomb de dois centros, fornecendo fórmulas explícitas para deslocamentos admissíveis e parâmetros de perturbação tanto em 2D quanto em 3D.

Resultados Numéricos
O algoritmo é validado em quatro casos de teste:

1. Átomo de Hidrogênio 2D ($V = -1/|x|$): O método enquadra com sucesso o estado fundamental exato $\lambda_1 = -1$. O intervalo inferior certificado estreita-se significativamente com o refinamento da malha, com os limites de deslocamento ótimo atingindo larguras de $0,002$ no nível mais fino.
2. Íon $H_2^+$ 2D: Resultados para o potencial de dois centros mostram convergência, com o intervalo certificado para o estado fundamental estreitando-se de $[-4,43, -1,317]$ para $[-2,36, -1,317]$ à medida que a malha é refinada.
3. Átomo de Hidrogênio 3D: Em um domínio esférico, o método enquadra o $\lambda_1 = -0,25$ exato. A largura do envoltório de Neumann certificado converge a $O(h^2)$.
4. Íon $H_2^+$ 3D: Para o íon molecular em um domínio cúbico, o método fornece um intervalo certificado $[-1,045, -0,544]$ (largura $0,501$) para o estado fundamental, que contém o valor de referência $\approx -0,5513$.

Em todos os casos, a lacuna do limite inferior CECR $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ é observada fechando à taxa prevista de $O(h^2)$ em malhas graduadas. O artigo nota que, embora as lacunas 3D sejam mais amplas do que em 2D devido ao decaimento mais lento da constante por elemento em malhas tetraédricas grosseiras, o comportamento de convergência é consistente com a teoria.

Significado e Reivindicações
O artigo afirma que este trabalho constitui o primeiro quadro de elementos finitos capaz de fornecer limites inferiores explícitos e convergentes para autovalores associados a potenciais singulares ilimitados. O significado reside em superar a "lacuna de coercividade" causada por potenciais atrativos singulares, o que anteriormente fazia com que técnicas padrão de limite inferior divergissem.

Os autores enfatizam que o método é explícito: todas as constantes (constantes de interpolação, limites de perturbação, parâmetros de deslocamento) são derivadas em forma fechada a partir dos dados da malha e do potencial. Embora a implementação atual utilize aritmética de dupla precisão, os autores afirmam explicitamente que a conversão para envoltórios totalmente verificados (seguros IEEE-754, aritmética intervalar) é uma etapa de verificação separada deixada para trabalhos futuros. O artigo também observa que, embora os limites atuais sejam baseados em projeção, eles servem como um ponto de partida rigoroso ($\rho$) para técnicas de pós-processamento mais afiadas, como limites de Lehmann–Goerisch, que poderiam estreitar ainda mais os intervalos.

O trabalho não reivindica resolver diretamente o problema de múltiplos elétrons, mas identifica a extensão para hamiltonianos magnéticos e sistemas de múltiplos elétrons com restrições de antissimetria como direções futuras naturais.