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

Resumen Técnico: Funcionales de Intercambio Semilocales a partir de la Condición de Intercambio Exacto para el Átomo de Hidrógeno

Planteamiento del Problema
Los funcionales de densidad semilocales estándar (GGA como PBE y meta-GGA como SCAN) sufren un defecto estructural fundamental: sus potenciales de intercambio ($v_x$) decaen exponencialmente en las regiones externas de sistemas finitos, fallando al reproducir la cola de Coulomb exacta de $-1/r$. Esta deficiencia conduce a una familia conocida de fallos, incluyendo la ausencia de estados de Rydberg ligados, iones negativos inestables, energías de ionización sistemáticamente subestimadas y una delocalización de carga excesiva debido al error de autointeracción (SIE). Si bien los potenciales de modelo (p. ej., LB94, mBJ) y los híbridos de rango separado abordan esto, los primeros carecen de consistencia variacional (no poseen un funcional de energía padre), y los segundos incurren en altos costos computacionales. Los intentos semilocales existentes, como el funcional de Armiento-Kümmel (AK13), mejoran las asíntotas pero no se derivan de condiciones exactas para el átomo de hidrógeno.

Metodología
Este trabajo revisita la ecuación diferencial ordinaria derivada por Gill y Pople (GP93) en 1993, la cual define la condición que un GGA debe satisfacer para reproducir el potencial de intercambio exacto del átomo de hidrógeno. Los autores resuelven esta ecuación completamente como un problema inverso para derivar un factor de mejora único, denominado factor GP93 ($F_{GP93}$).

La metodología procede en tres capas:

1. Exactitud Hidrogénica: El factor GP93 se construye para ser exacto en la densidad del hidrógeno 1s. Mediante el escalamiento de coordenadas uniforme, produce autovalores exactos ($\epsilon_{1s} = -Z^2/2$) para todos los iones hidrogénicos 1s sin parámetros empíricos.
2. Comportamiento Asintótico: En el régimen de alto gradiente, el factor GP93 exhibe un crecimiento no acotado de la forma $O[s(\ln s)^{2/3}]$, distinto del crecimiento $s \ln s$ de AK13. Este crecimiento específico asegura que el potencial de intercambio recupere la cola exacta de $-1/r$ en la densidad de hidrógeno.
3. Mecanismo de Conmutación (Switching): Para evitar que el comportamiento divergente del factor GP93 desestabilice sistemas de muchos electrones, los autores lo integran en una forma meta-GGA "conmutada".
   • $\alpha$-Switch: Utiliza el indicador $\alpha = (\tau - \tau^W)/\tau_{unif}$ para detectar regiones de tipo un-electrón. El término GP93 está activo solo donde $\alpha \approx 0$ (regiones de un solo orbital).
   • $q$-Switch (Conmutación Triple): Para distinguir las colas atómicas (donde la cola es necesaria) de los centros de enlaces covalentes (donde la cola causa inestabilidad), se introduce un indicador de Laplaciano reducido $q = \nabla^2 n / (4(3\pi^2)^{2/3}n^{5/3})$. Las colas atómicas tienen $q \to +\infty$, mientras que los centros de enlace tienen $q \lesssim 0$. Una función de conmutación $D(q)$ suprime el término GP93 en los centros de enlace mientras lo preserva en las colas atómicas.
   • Independencia de Orbitales: La densidad de energía cinética $\tau$ se reemplaza por un funcional de densidad de energía cinética orbital-libre (KEDF) basado en la construcción de Perdew-Constantin (PC07), que impone el límite $\tau^W \le \tau$. Esto convierte a todo el funcional en un funcional de densidad puro, eliminando la necesidad de un resolvedor de Kohn–Sham generalizado (gKS).

Resultos Clave

• Exactitud Hidrogénica: El funcional reproduce el autovalor exacto $\epsilon_{1s} = -Z^2/2$ para iones hidrogénicos ($Z=1$ a $5$) al nivel de milihartree.
• Estados de Rydberg Ligados: En sistemas con HOMOs de tipo capa-$s$ de un solo electrón (H, He, y Li/Na/K con pseudopotenciales de núcleo grande), el funcional recupera una serie de estados virtuales de Kohn–Sham ligados (4–8 niveles). En contraste, los funcionales semilocales estándar (PBE, SCAN, TPSS, M06-L) e híbridos globales (B3LY de B3LYP, PBE0) no vinculan cero estados de Rydberg en estos sistemas.
• Energías de Ionización: Para los sistemas de capa-$s$ probados, el autovalor más alto ocupado ($-\epsilon_{HOMO}$) está significativamente más cerca de los potenciales de ionización experimentales (Error Absoluto Medio $\approx 0.80$ eV) en comparación con PBE, SCAN o B3LYP.
• Selectividad: El funcional suprime correctamente el término GP93 para átomos de capa-$p$ (Ne, Ar) donde $\alpha > 0$, resultando en cero estados de Rydberg ligados, consistente con el diseño.
• Estabilidad Molecular: La conmutación triple logra suprimir la divergencia numérica en los centros de enlaces covalentes (p. ej., $H_2$, $N_2$) donde el factor GP93 sin conmutar fallaría, manteniendo al mismo tiempo la cola asintótica correcta para sistemas atómicos.
• Compatibilidad de Correlación: El funcional se empareja con el funcional de correlación LYP, que se anula para sistemas de un solo electrón, preservando así la exactitud hidrogénica del componente de intercambio.

Significado y Reivindicaciones
El artículo afirma haber construido un funcional de intercambio semilocal no empírico que recupera la cola exacta de $-1/r$ y los estados de Rydberg ligados para sistemas de capa-$s$ de un solo electrón, un fallo cualitativo de todos los funcionales semilocales existentes, incluido SCAN.

Los autores delimitan explícitamente el alcance de sus afirmaciones:

• Dominio de Éxito: El funcional está diseñado para propiedades de un solo electrón (estados de Rydberg, asíntotas, niveles hidrogénicos, SIE) en sistemas con HOMOs de capa-$s$ de un solo electrón.
• Limitaciones:
  • Termoquímica: El funcional no apunta a las energías de atomización. Imponer la exactitud hidrogénica desplaza la energía absoluta del átomo de H lejos de los valores compensados de error optimizados en PBE, lo que lo hace inadecuado para la termoquímica que involucra a H.
  • Sistemas de capa-$p$: La mejora asintótica no se extiende a los HOMOs de capa-$p$ (p. ej., Ne, Ar) debido al mecanismo de $\alpha$-switch.
  • IP Cuantitativas: Aunque es preciso para He, las energías de ionización para sistemas más pesados pueden contener errores sistemáticos debido a la calibración referenciada al hidrógeno.
  • Propiedades de Respuesta: Las derivadas superiores (p. ej., polarizabilidades) están fuera del alcance debido a la sensibilidad de los límites de conmutación.
• Implementación: El artículo demuestra el funcional de forma no autoconsistente sobre densidades fijas. Señala que una implementación totalmente autoconsistente en bases Gaussianas es inestable debido al término del Laplaciano ($v_{lapl}$), pero argumenta que el método es naturalmente adecuado para implementaciones de pseudopotencial en espacio real donde el Laplaciano puede manejarse de forma estable.

El trabajo establece que un único dato de entrada no empírico —la condición de exactitud hidrogénica de GP93— puede generar simultáneamente la cola asintótica correcta y los niveles hidrogénicos exactos, ofreciendo una vía hacia la física de Rydberg con un costo computacional semilocal sin la necesidad de intercambio exacto no local.