Semilocal exchange functionals from the exact-exchange condition for the hydrogen atom: Hydrogenic exactness and recovery of Rydberg-like bound states

Resumo Técnico: Funcionais de Troca Semilocais a partir da Condição de Troca Exata para o Átomo de Hidrogênio

Enunciado do Problema
Os funcionais de densidade semilocais padrão (GGAs como PBE e meta-GGAs como SCAN) sofrem de um defeito estrutural fundamental: seus potenciais de troca ($v_x$) decaem exponencialmente nas regiões externas de sistemas finitos, falhando em reproduzir a cauda coulombiana exata de $-1/r$. Essa deficiência leva a uma família conhecida de falhas, incluindo a ausência de estados de Rydberg ligados, íons negativos instáveis, energias de ionização sistematicamente subestimadas e delocalização excessiva de carga devido ao erro de autointeração (SIE). Embora potenciais de modelo (ex: LB94, mBJ) e híbridos de separação de intervalo (range-separated hybrids) abordem isso, os primeiros carecem de consistência variacional (não possuem um funcional de energia pai), e os últimos incorrem em altos custos computacionais. Tentativas semilocais existentes, como o funcional Armiento-Kümmel (AK13), melhoram as assíntotas, mas não são derivadas de condições exatas para o átomo de hidrogênio.

Metodologia
Este trabalho revisita a equação diferencial ordinária derivada por Gill e Pople (GP93) em 1993, que define a condição que um GGA deve satisfazer para reproduzir o potencial de troca exato do átomo de hidrogênio. Os autores resolvem esta equação completamente como um problema inverso para derivar um fator de aumento único, denominado fator GP93 ($F_{GP93}$).

A metodologia procede em três camadas:

1. Exatidão Hidrogênica: O fator GP93 é construído para ser exato na densidade 1s do hidrogênio. Por escalonamento de coordenadas uniforme, ele produz autovalores exatos ($\epsilon_{1s} = -Z^2/2$) para todos os íons hidrogênicos 1s sem parâmetros empíricos.
2. Comportamento Assintótico: No regime de grande gradiente, o fator GP93 exibe um crescimento ilimitado da forma $O[s(\ln s)^{2/3}]$, distinto do crescimento $s \ln s$ do AK13. Esse crescimento específico garante que o potencial de troca recupere a cauda exata $-1/r$ na densidade do hidrogênio.
3. Mecanismo de Alternância (Switching): Para evitar que o comportamento divergente do fator GP93 desestabilize sistemas de muitos elétrons, os autores o inserem em uma forma meta-GGA com "alternância" (switched).
   • Alternância $\alpha$: Utiliza o indicador $\alpha = (\tau - \tau^W)/\tau_{unif}$ para detectar regiões semelhantes a um elétron. O termo GP93 é ativo apenas onde $\alpha \approx 0$ (regiões de um único orbital).
   • Alternância $q$ (Alternância Tripla): Para distinguir caudas atômicas (onde a cauda é necessária) de centros de ligações covalentes (onde a cauda causa instabilidade), é introduzido um indicador de Laplaciano reduzido $q = \nabla^2 n / (4(3\pi^2)^{2/3}n^{5/3})$. Caudas atômicas possuem $q \to +\infty$, enquanto centros de ligação possuem $q \lesssim 0$. Uma função de alternância $D(q)$ suprime o termo GP93 em centros de ligação, preservando-o em caudas atômicas.
   • Independência de Orbital: A densidade de energia cinética $\tau$ é substituída por um funcional de densidade de energia cinética livre de orbitais (KEDF) baseado na construção de Perdew-Constantin (PC07), que impõe o limite $\tau^W \le \tau$. Isso torna todo o funcional um funcional de densidade pura, eliminando a necessidade de um solver de Kohn–Sham generalizado (gKS).

Principais Resultados

• Exatidão Hidrogênica: O funcional reproduz o autovalor exato $\epsilon_{1s} = -Z^2/2$ para íons hidrogênicos ($Z=1$ a $5$) ao nível do milihartree.
• Estados de Rydberg Ligados: Em sistemas com HOMOs de casca $s$ semelhantes a um elétron (H, He, e Li/Na/K com pseudopotenciais de núcleo grande), o funcional recupera uma série de estados virtuais de Kohn–Sham ligados (4–8 níveis). Em contraste, funcionais semilocais padrão (PBE, SCAN, TPSS, M06-L) e híbridos globais (B3LYx, PBE0) não ligam nenhum estado de Rydberg nesses sistemas.
• Energias de Ionização: Para os sistemas de casca $s$ testados, o autovalor mais ocupado ($-\epsilon_{HOMO}$) está significativamente mais próximo das energias de ionização experimentais (Erro Médio Absoluto $\approx 0,80$ eV) comparado a PBE, SCAN ou B3LYP.
• Seletividade: O funcional suprime corretamente o termo GP93 para átomos de casca $p$ (Ne, Ar) onde $\alpha > 0$, resultando em zero estados de Rydberg ligados, consistente com o projeto.
• Estabilidade Molecular: A alternância tripla suprime com sucesso a divergência numérica em centros de ligações covalentes (ex: $H_2$, $N_2$) onde o fator GP93 não alternado falharia, mantendo a cauda assintótica correta para sistemas atômicos.
• Compatibilidade de Correlação: O funcional é pareado com o funcional de correlação LYP, que desaparece para sistemas de um elétron, preservando assim a exatidão hidrogênica do componente de troca.

Significância e Alegações
O artigo afirma ter construído um funcional de troca semilocal não empírico que recupera a cauda exata $-1/r$ e estados de Rydberg ligados para sistemas de casca $s$ semelhantes a um elétron, uma falha qualitativa de todos os funcionais semilocais existentes, incluindo o SCAN.

Os autores delimitam explicitamente o escopo de suas alegações:

• Domínio de Sucesso: O funcional é projetado para propriedades de um elétron (estados de Rydberg, assíntotas, níveis hidrogênicos, SIE) em sistemas com HOMOs de casca $s$ semelhantes a um elétron.
• Limitações:
  • Termoquímica: O funcional não visa energias de atomização. Impor a exatidão hidrogênica desloca a energia absoluta do átomo de H para longe dos valores compensados de erro otimizados no PBE, tornando-o inadequado para termoquímica envolvendo H.
  • Sistemas de casca $p$: A melhoria assintótica não se estende a HOMOs de casca $p$ (ex: Ne, Ar) devido ao mecanismo de alternância $\alpha$.
  • IPs Quantitativas: Embora preciso para He, as energias de ionização para sistemas mais pesados podem conter erros sistemáticos devido à calibração referenciada no hidrogênio.
  • Propriedades de Resposta: Derivadas superiores (ex: polarizabilidades) estão fora do escopo devido à sensibilidade das fronteiras de alternância.
• Implementação: O artigo demonstra o funcional de forma não autoconsistente em densidades fixas. Nota-se que uma implementação totalmente autoconsistente em bases de Gauss é instável devido ao termo do Laplaciano ($v_{lapl}$), mas argumenta-se que o método é naturalmente adequado para implementações de pseudopotencial em espaço real, onde o Laplaciano pode ser tratado de forma estável.

Este trabalho estabelece que uma única entrada não empírica — a condição de exatidão hidrogênica GP93 — pode gerar simultaneamente a cauda assintótica correta e níveis hidrogênicos exatos, oferecendo um caminho para a física de Rydberg ao custo computacional semelocante, sem a necessidade de troca exata não local.