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

Technische Samenvatting: Semilokale Uitwisselingsfunctionalen gebaseerd op de Exacte-Uitwisselingsconditie voor het Waterstofatoom

Probleemstelling
Standaard semilocale dichtheidsfunctionalen (GGA's zoals PBE en meta-GGA's zoals SCAN) lijden aan een fundamenteel structureel defect: hun uitwisselingspotentialen ($v_x$) vervallen exponentieel in de buitenste regio's van eindige systemen, waardoor ze er niet in slagen de exacte Coulombische $-1/r$ staart te reproduceren. Dit tekortschieten leidt tot een bekende familie van fouten, waaronder de afwezigheid van gebonden Rydberg-toestanden, instabiele negatieve ionen, systematisch onderschatte ionisatie-energieën en excessieve ladingsdelokalisatie door zelfinteractiefout (SIE). Hoewel modelpotentialen (bijv. LB94, mBJ) en range-separated hybriden dit adresseren, missen de eerste een variationele consistentie (geen ouder energiefunctional) en de laatste hoge computationele kosten. Bestaande semilocale pogingen, zoals de Armiento-Kümmel (AK13) functional, verbeteren de asymptotica maar zijn niet afgeleid van exacte condities voor het waterstofatoom.

Methodologie
Dit werk herbezoekt de gewone differentiaalvergelijking afgeleid door Gill en Pople (GP93) in 1993, die de conditie definieert die een GGA moet vervullen om het exacte uitwisselingspotentiaal van het waterstofatoom te reproduceren. De auteurs lossen deze vergelijking volledig op als een invers probleem om een unieke verbeterfactor af te leiden, aangeduid als de GP93-factor ($F_{GP93}$).

De methodologie verloopt in drie lagen:

1. Hydrogene Exactheid: De GP93-factor wordt geconstrueerd om exact te zijn op de waterstof 1s-dichtheid. Door uniforme coördinaatschaling levert het exacte eigenwaarden ($\epsilon_{1s} = -Z^2/2$) voor alle waterstofachtige 1s-ionen zonder empirische parameters.
2. Asymptotisch Gedrag: In het regime met grote gradiënten vertoont de GP93-factor een onbegrensde groei van de vorm $O[s(\ln s)^{2/3}]$, die verschilt van de $s \ln s$ groei van AK13. Deze specifieke groei zorgt ervoor dat het uitwisselingspotentiaal de exacte $-1/r$ staart herstelt op de waterstofdichtheid.
3. Schakelmechanisme: Om te voorkomen dat het divergente gedrag van de GP93-factor veel-elektronen-systemen destabiliseert, integreren de auteurs deze in een "geschakelde" meta-GGA vorm.
   • $\alpha$-Switch: Gebruikt de indicator $\alpha = (\tau - \tau^W)/\tau_{unif}$ om één-elektron-achtige regio's te detecteren. De GP93-term is alleen actief waar $\alpha \approx 0$ (één-orbitaal regio's).
   • $q$-Switch (Triple Switch): Om atomaire staarten (waar de staart nodig is) te onderscheiden van covalente bindingscentra (waar de staart instabiliteit veroorzaakt), wordt een gereduceerde Laplaciaan-indicator $q = \nabla^2 n / (4(3\pi^2)^{2/3}n^{5/3})$ geïntroduceerd. Atomaire staarten hebben $q \to +\infty$, terwijl bindingscentra $q \lesssim 0$ hebben. Een schakelfunctie $D(q)$ onderdrukt de GP93-term in bindingscentra terwijl deze behouden blijft in atomaire staarten.
   • Orbitaalonafhankelijkheid: De kinetische energiedichtheid $\tau$ wordt vervangen door een orbital-free kinetische-energie-dichtheidsfunctional (KEDF) gebaseerd op de Perdew-Constantin (PC07) constructie, die de begrenzing $\tau^W \le \tau$ afdwingt. Dit maakt de gehele functional een zuivere dichtheidsfunctional, waardoor de noodzaak voor een generalized Kohn–Sham (gKS) solver vervalt.

Belangrijkste Resultaten

• Hydrogene Exactheid: De functional reproduceert de exacte eigenwaarde $\epsilon_{1s} = -Z^2/2$ voor waterstofachtige ionen ($Z=1$ tot $5$) op millihartree-niveau.
• Rydberg-achtige Gebonden Toestanden: In systemen met één-elektron-achtige $s$-schil HOMOs (H, He, en Li/Na/K met grote-kern pseudopotentialen), herstelt de functional een reeks gebonden virtuele Kohn–Sham toestanden (4–8 niveaus). In contrast hiermee binden standaardfunctionalen (PBE, SCAN, TPSS, M06-L) en globale hybriden (B3LYP, PBE0) nul Rydberg-toestanden in deze systemen.
• Ionisatie-energieën: Voor de geteste $s$-schil systemen ligt de hoogste bezette eigenwaarde ($-\epsilon_{HOMO}$) aanzienlijk dichter bij de experimentele ionisatiepotentialen (Gemiddelde Absolute Fout $\approx 0,80$ eV) vergeleken met PBE, SCAN of B3LYP.
• Selectiviteit: De functional onderdrukt correct de GP93-term voor $p$-schil atomen (Ne, Ar) waar $\alpha > 0$, wat resulteert in nul gebonden Rydberg-toestanden, consistent met het ontwerp.
• Moleculaire Stabiliteit: De triple switch onderdrukt succesvol numerieke divergentie in covalente bindingscentra (bijv. $H_2$, $N_2$) waar de ongeschakelde GP93-factor zou falen, terwijl de correcte asymptotische staart voor atomaire systemen behouden blijft.
• Correlatie Compatibiliteit: De functional wordt gekoppeld aan de LYP correlatiefunctional, die verdwijnt voor één-elektronsystemen, waardoor de hydrogene exactheid van de uitwisselingscomponent behouden blijft.

Betekenis en Claims
Het artikel claimt dat zij een niet-empirische, semilocale uitwisselingsfunctional hebben geconstrueerd die de exacte $-1/r$ staart en gebonden Rydberg-toestanden herstelt voor één-elektron-achtige $s$-schil systemen, een kwalitatieve fout van alle bestaande semilocale functionalen inclusief SCAN.

De auteurs beperken de reikwijdte van hun claims expliciet:

• Succesdomein: De functional is ontworpen voor één-elektronen eigenschappen (Rydberg-toestanden, asymptotica, hydrogene niveaus, SIE) in systemen met één-elektron-achtige $s$-schil HOMOs.
• Beperkingen:
  • Thermochemie: De functional richt zich niet op atomisatie-energieën. Het afdwingen van hydrogene exactheid verschuift de absolute energie van het H-atoom weg van de fout-gecompenseerde waarden die geoptimaliseerd zijn in PBE, waardoor het ongeschikt is voor thermochemie waarbij H betrokken is.
  • $p$-schil Systemen: De asymptotische verbetering strekt zich niet uit tot $p$-schil HOMOs (bijv. Ne, Ar) vanwege de $\alpha$-switch mechanisme.
  • Kwantitatieve IP's: Hoewel accuraat voor He, kunnen de ionisatie-energieën voor zwaardere systemen systematische fouten bevatten door de op waterstof gebaseerde kalibratie.
  • Respons-eigenschappen: Hogere derivaten (bijv. polariseerbaarheid) vallen buiten de scope vanwege de gevoeligheid van de schakelgrenzen.
• Implementatie: Het artikel demonstreert de functional niet-zelfconsistent op vaste dichtheden. Er wordt opgemerkt dat een volledig zelfconsistente implementatie in Gaussische basissen instabiel is vanwege de Laplaciaan-term ($v_{lapl}$), maar zij argumenteren dat de methode van nature geschikt is voor real-space pseudopotentiaal implementaties waar de Laplaciaan stabiel kan worden afgehandeld.

Dit werk stelt vast dat een enkele niet-empirische input—de GP93 hydrogene exactheidconditie—gelijktijdig de correcte asymptotische staart en exacte hydrogene niveaus kan genereren, wat een pad biedt naar Rydberg-fysica tegen semilocale computationele kosten zonder de noodzaak van niet-lokale exacte uitwisseling.