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

Sintesi Tecnica: Funzionali di Scambio Semilocali dalla Condizione di Scambio Esatto per l'Atomo di Idrogeno

Problematica
I comuni funzionali di densità semilocali (GGA come PBE e meta-GGA come SCAN) soffrono di un difetto strutturale fondamentale: i loro potenziali di scambio ($v_x$) decadono esponenzialmente nelle regioni esterne dei sistemi finiti, fallendo nel riprodurre la coda coulombiana esatta $-1/r$. Questa carenza porta a una nota famiglia di fallimenti, tra cui l'assenza di stati di Rydberg legati, ioni negativi instabili, energie di ionizzazione sistematicamente sottostimate ed eccessiva delocalizzazione della carica dovuta all'errore di auto-interazione (SIE). Sebbene i potenziali modello (es. LB94, mBJ) e gli ibridi a separazione di raggio affrontino il problema, i primi mancano di consistenza variazionale (non possiedono un funzionale di energia genitore), mentre i secondi comportano costi computazionali elevati. I tentativi semilocali esistenti, come il funzionale Armiento-Kümmel (AK13), migliorano l'asintotica ma non sono derivati dalle condizioni esatte per l'atomo di idrogeno.

Metodologia
Questo lavoro rivisita l'equazione differenziale ordinaria derivata da Gill e Pople (GP93) nel 1993, che definisce la condizione che un GGA deve soddisfare per riprodurre il potenziale di scambio esatto dell'atomo di idrogeno. Gli autori risolvono completamente questa equazione come un problema inverso per derivare un fattore di potenziamento unico, denominato fattore GP93 ($F_{GP93}$).

La metodologia procede in tre strati:

1. Esattezza Idrogenica: Il fattore GP93 è costruito per essere esatto sulla densità dell'idrogeno 1s. Tramite la scalatura uniforme delle coordinate, esso fornisce autovalori esatti ($\epsilon_{1s} = -Z^2/2$) per tutti gli ioni idrogenoidi 1s senza parametri empirici.
2. Comportamento Asintotico: Nel regime di grande gradiente, il fattore GP93 presenta una crescita illimitata della forma $O[s(\ln s)^{2/3}]$, distinta dalla crescita $s \ln s$ di AK13. Questa specifica crescita assicura che il potenziale di scambio recuperi la coda esatta $-1/r$ sulla densità dell'idrogeno.
3. Meccanismo di Commutazione (Switching): Per evitare che il comportamento divergente del fattore GP93 destabilizzi i sistemi multielettronici, gli autori lo inseriscono in una forma meta-GGA "commutata":
   • Switch $\alpha$: Utilizza l'indicatore $\alpha = (\tau - \tau^W)/\tau_{unif}$ per rilevare regioni di tipo mono-orbitale. Il termine GP93 è attivo solo dove $\alpha \approx 0$ (regioni a singolo orbitale).
   • Switch $q$ (Triple Switch): Per distinguere le code atomiche (dove la coda è necessaria) dai centri di legame covalente (dove la coda causa instabilità), viene introdotto un indicatore del Laplaciano ridotto $q = \nabla^2 n / (4(3\pi^2)^{2/3}n^{5/3})$. Le code atomiche hanno $q \to +\infty$, mentre i centri di legame hanno $q \lesssim 0$. Una funzione di commutazione $D(q)$ sopprime il termine GP93 nei centri di legame preservandolo però nelle code atomiche.
   • Indipendenza Orbitale: La densità di energia cinetica $\tau$ è sostituita da un funzionale della densità di energia cinetica orbitale-free (KEDF) basato sulla costruzione di Perdew-Constantin (PC07), che impone il vincolo $\tau^W \le \tau$. Ciò rende l'intero funzionale un puro funzionale di densità, eliminando la necessità di un risolutore di Kohn–Sham generalizzato (gKS).

Risultati Chiave

• Esattezza Idrogenica: Il funzionale riproduce l'autovalore esatto $\epsilon_{1s} = -Z^2/2$ per gli ioni idrogenoidi ($Z=1$ a $5$) a livello di millihartree.
• Stati di Rydberg Legati: Nei sistemi con HOMO a carattere $s$-shell a un solo elettrone (H, He, e Li/Na/K con pseudopotenziali a nucleo grande), il funzionale recupera una serie di stati virtuali di Kohn–Sham legati (4–8 livelli). Al contrario, i funzionali standard (PBE, SCAN, TPSS, M06-L) e gli ibridi globali (B3LYP, PBE0) non legano alcuno stato di Rydberg in questi sistemi.
• Energie di Ionizzazione: Per i sistemi $s$-shell testati, l'autovalore più alto occupato ($-\epsilon_{HOMO}$) è significativamente più vicino alle energie di ionizzazione sperimentali (Errore Assoluto Medio $\approx 0.80$ eV) rispetto a PBE, SCAN o B3LYP.
• Selettività: Il funzionale sopprime correttamente il termine GP93 per gli atomi a $p$-shell (Ne, Ar) dove $\alpha > 0$, risultando in zero stati di Rydberg legati, in linea con il design.
• Stabilità Molecolare: Il "triple switch" riesce a sopprimere la divergenza numerica nei centri di legame covalente (es. $H_2$, $N_2$) dove il fattore GP93 non commutato fallirebbe, mantenendo al contempo la corretta coda asintotica per i sistemi atomici.
• Compatibilità con la Correlazione: Il funzionale è accoppiato con il funzionale di correlazione LYP, che svanisce per i sistemi a un solo elettrone, preservando così l'esattezza idrogenica della componente di scambio.

Significato e Rivendicazioni
Il lavoro afferma di aver costruito un funzionale di scambio semilocale non empirico che recupera la coda asintotica esatta $-1/r$ e gli stati di Rydberg legati per i sistemi $s$-shell a un solo elettrone, un fallimento qualitativo di tutti i funzionali semilocali esistenti, incluso SCAN.

Gli autori delimitano esplicitamente l'ambito delle loro rivendicazioni:

• Dominio di Successo: Il funzionale è progettato per le proprietà a un solo elettrone (stati di Rydberg, asintotica, livelli idrogenoidi, SIE) in sistemi con HOMO $s$-shell a un solo elettrone.
• Limitazioni:
  • Termochimica: Il funzionale non mira alle energie di atomizzazione. Imporre l'esattezza idrogenica sposta l'energia assoluta dell'atomo di H lontano dai valori ottimizzati per la compensazione dell'errore in PBE, rendendolo inadatto alla termochimica che coinvolge l'idrogeno.
  • Sistemi a $p$-shell: Il miglioramento asintotico non si estende agli HOMO a $p$-shell (es. Ne, Ar) a causa del meccanismo di switch $\alpha$.
  • IP Quantitative: Sebbene accurato per l'elio, le energie di ionizzazione per sistemi più pesanti possono contenere errori sistematici dovuti alla calibrazione riferita all'idrogeno.
  • Proprietà di Risposta: Le derivate superiori (es. polarizzabilità) sono fuori dallo scopo a causa della sensibilità dei confini dello switch.
• Implementazione: Il documento presenta il funzionale in modo non auto-coerente su densità fisse. Si nota che un'implementazione completamente auto-coerente in basi Gaussiane è instabile a causa del termine Laplaciano ($v_{lapl}$), ma si argomenta che il metodo è naturalmente adatto a implementazioni in spazio reale con pseudopotenziali, dove il Laplaciano può essere gestito in modo stabile.

Il lavoro stabilisce che un singolo input non empirico — la condizione di esattezza idrogenica GP93 — può generare simultaneamente la corretta coda asintotica e i livelli idrogenoidi esatti, offrendo una via verso la fisica di Rydberg al costo computazionale semilocale senza la necessità di uno scambio esatto non locale.