Excited-state Properties Beyond the Excitation Energy from Orbital-Optimized Density Functional Calculations II: Absorption Spectra

Deze studie breidt Löwdins formalisme voor niet-orthogonale determinanten uit naar orbitaal-geoptimaliseerde dichtheidsfunctionaaltheoretische berekeningen binnen het projector augmented-wave kader, waarbij wordt aangetoond dat hoewel de methode absorptiespectra en piekintensiteiten voor single-determinant aangeslagen toestanden kwalitatief reproduceert, het moeite heeft met multi-configurational toestanden en geen systematische verbetering vertoont door exacte uitwisseling of zelfinteractiecorrecties.

De kern

Het Grote Plaatje: Een Snapshot van een Dansend Molecuul

Stel je een molecuul voor als een kleine, complexe dansgroep. Wanneer je licht op het schijnt, springen de dansers (elektronen) naar hogere energieniveaus, waardoor hun danspassen veranderen. Wetenschappers willen twee dingen voorspellen over deze sprong:

1. De Energie: Hoeveel energie is er nodig om de sprong te maken? (Zoals de hoogte van een sprong).
2. De Helderheid: Hoe helder is de lichtflits wanneer ze springen? (Dit wordt de "oscillatorsterkte" of "absorptie-intensiteit" genoemd).

Lange tijd waren computerprogramma's goed in het voorspellen van de hoogte van de sprong, maar kregen ze de helderheid vaak fout. Dit paper introduceert een nieuwe manier om deze sprongen te berekenen, genaamd Orbital-Optimized (OO) Density Functional Theory.

Het Probleen: De "Niet-Orthogonale" Dansvloer

In standaard computermodellen gaan wetenschappers er meestal van uit dat de "grondtoestand" (de rustende dans) en de "geëxciteerde toestand" (de springende dans) volledig gescheiden zijn en elkaar niet overlappen, alsof twee mensen in verschillende kamers staan. Dit maakt de wiskunde eenvoudig.

Deze nieuwe "Orbital-Optimized" methode is echter realistischer. Het laat de elektronen zichzelf specifiek herordenen voor de sprong. Het probleem is dat deze herordening betekent dat de "rustende" toestand en de "springende" toestand niet langer in aparte kamers zijn; ze bevinden zich nu in dezelfde kamer, met een lichte overlap.

De Analogie: Stel je voor dat je de afstand probeert te meten tussen twee mensen die elkaar omhelzen. Als je ervan uitgaat dat ze apart van elkaar staan, is je meting fout. Omdat deze toestanden "elkaar omhelzen" (overlap hebben), wordt het berekenen van hoe helder de lichtflits zal zijn, erg lastig. Eerdere studies controleerden vooral of de spronghoogte klopte, maar ze controleerden niet of de berekening van de helderheid accuraat was wanneer de toestanden overlapten.

Wat Ze Deden: De Nieuwe Methode Testen

De onderzoekers namen deze nieuwe methode en testten deze op een kleine groep moleculen (water, ammoniak, formaldehyde, methanol en ethyleen). Ze wilden zien:

1. Voorspelt deze methode de helderheid van de lichtflits correct?
2. Maakt het uit welke "gereedschappen" (wiskundige basissets) ze gebruiken voor de berekening?
3. Fixen het veranderen van de "spelregels" (de wiskundige formules die het gedrag van elektronen beschrijven) de fouten in de helderheid?

De Resultaten: Een Gemengd Resultaat

Dit is wat ze vonden, simpel uitgelegd:

1. De "Gereedschapskist" Maakt Uit (Basissets)
Beschouw de "basisset" als de resolutie van een camera.

• Lage Resolutie (Simpele tools): Als je een simpele camera gebruikt, mis je de fijne details van de "diffuse" elektronen (die ver weg van het molecuul zweven, zoals een Rydberg-toestand). De berekening van de helderheid is dan een beetje af.
• Hoge Resolutie (Complexe tools): Wanneer ze zeer gedetailleerde tools gebruikten (zoals Plane Waves of double-augmented sets), werden de resultaten voor de helderheid veel consistenter.
• Takeaway: Om de helderheid juist te krijgen, heb je een camera met een hoge resolutie nodig, vooral voor die elektronen die ver weg zweven.

2. De "Regels" Fixen Niet Alles (Exchange-Correlation Functionals)
Wetenschappers probeerden de "regels" van de simulatie te veranderen (door verschillende wiskundige formules zoals PBE, PBE0 of correcties voor zelfinteractie te gebruiken).

• Het Resultaat: Het veranderen van de regels hielp in sommige gevallen bij het corrigeren van de hoogte van de sprong (energie), maar het heeft de helderheid niet consistent verbeterd.
• Analogie: Stel je voor dat je probeert een wazige foto te repareren. Je hebt geprobeerd de lens, de belichting en het filter te veranderen. Soms werd de foto scherper, maar vaak bleef de helderheid nog steeds fout. Er was geen enkele "magische regel" die de helderheid voor elk molecuul oploste.

3. De Echte Boosdoener: De "Solo" versus de "Groep" (Single vs. Multi-Configurational)
Dit is de belangrijkste bevinding. De methode werkte geweldig wanneer de geëxciteerde toestand als een solist was (een enkele, duidelijke configuratie).

• Solisten: Voor eenvoudige sprongen (zoals in ammoniak) voorspelde de methode de helderheid perfect.
• Groepsdansers: Voor complexe sprongen waarbij de elektronendans een mix is van veel verschillende mogelijkheden tegelijkertijd (multi-configurational), faalde de methode in het correct voorspellen van de helderheid.
• De Specifieke Fouten: Ze vonden grote fouten bij water, formaldehyde en ethyleen. In deze gevallen is de geëxciteerde toestand een rommelige mix van verschillende danspassen. Omdat het computermodel de toestand dwingt om eruit te zien als een enkele, schone "solo-pas", krijgt het de helderheid fout.
• Het Overlap-Probleem: Ze controleerden of de "omhelzing" (overlap) tussen de grondtoestand en de geëxciteerde toestand de fout veroorzaakte. Ze ontdekten dat zelfs toen ze de overlap veranderden, de fout in de helderheid hetzelfde bleef. De overlap was dus niet het hoofdzakelijke probleem; de aard van de dans ("solo" versus "groep") was dat wel.

Vergelijking met de Oude Manier (LR-TDDFT)

Ze vergeleken hun nieuwe methode met de standaardmethode (LR-TDDFT):

• Standaardmethode: Goed in het voorspellen van de helderheid van complexe "groepsdansen" (zoals de heldere flitsen in ethyleen), maar slecht in het voorspellen van de energie van de "solisten" (Rydberg-toestanden) die ver weg zweven.
• Nieuwe Methode (OO): Geweldig in het voorspellen van de energie van de "solisten" (Rydberg-toestanden), maar heeft moeite met de helderheid van de complexe "groepsdansen".

De Kern van het Verhaal

Dit paper laat zien dat de nieuwe "Orbital-Optimized" methode een krachtig hulpmiddel is voor het voorspellen van hoe moleculen licht absorberen, maar met een kanttekening:

• Het werkt zeer goed voor eenvoudige, enkelvoudige danspassen (Rydberg-toestanden).
• Het heeft moeite met situaties waarin de dans een complexe mix is van vele passen (multi-configurational states).
• Het simpelweg veranderen van de wiskundige formules of de "cameraresolutie" lost de fouten voor de complexe dansen niet op. Om die te fixen, zouden we een methode nodig hebben die complexe groepsdansen kan aan, en niet alleen solo's.

Kortom: de methode is een grote stap voorwaarts voor het begrijpen van de "hoogte" van elektronensprongen en de helderheid van eenvoudige sprongen, maar het heeft nog steeds hulp nodig om de helderheid van de meest complexe, chaotische elektronendansen nauwkeurig te voorspellen.

Technische samenvatting

Technische Samenvatting: Eigenschappen van de aangeslagen toestand voorbij de excitatie-energie uit baan-geoptimaliseerde dichtheidsfunctionaalberekeningen II: Absorptiespectra

Probleemstelling
Baan-geoptimaliseerde (OO) dichtheidsfunctionaalberekeningen bieden een staat-specifieke benadering voor het beschrijven van aangeslagen toestanden door gebruik te maken van variationeel geoptimaliseerde orbitalen voor elke toestand, in plaats van deze te behandelen als lineaire responsen van de grondtoestand. Hoewel deze benadering heeft aangetoond de excitatie-energieën te verbeteren, met name voor Rydberg-toestanden waar standaard lineaire-respons tijd-afhankelijke dichtheidsfunctionaaltheorie (LR-TDDFT) vaak faalt vanwege de incorrecte asymptotische potentiaalgedraging, introduceert het aanzienlijke complicaties voor het berekenen van transitie-eigenschappen. Omdat elke elektronische toestand wordt beschreven door een distinct pakket moleculaire orbitalen, zijn OO-toestanden over het algemeen niet-orthogonaal. Deze niet-orthogonaliteit bemoeilijkt de evaluatie van transitiedipoolmomenten (TDM's) en oscillatorsterktes. Eerdere benchmarks hebben zich grotendeels gericht op excitatie-energieën of beperkte oscillatorsterkte-berekeningen tot de eerste aangeslagen toestand (HOMO–LUMO-transitie). Bijgevolg blijft de prestatie van OO-methoden voor oscillatorsterktes over een breder spectrum van valentie- en Rydberg-excitaties, en de impact van niet-orthogonaliteit op deze eigenschappen, onvoldoende begrepen.

Methodologie
Deze studie breidt de toepassing van OO-dichtheidsfunctionaalberekeningen uit naar de computation van absorptiespectra voor een set kleine moleculen (H₂O, CH₂O, NH₃, CH₃OH en C₂H₄), waarbij zowel valentie- als Rydberg-toestanden tot 10 eV worden behandeld.

• Theoretisch kader: De auteurs passen een directe baan-optimalisatieschema toe om aangeslagen toestanden te lokaliseren als stationaire punten (zadelpunten) van de energiefunctie, om variatieve instorting te vermijden. Om de niet-orthogonaliteit van de resulterende toestanden te behandelen, wordt Löwdin's formalisme voor niet-orthogonale determinanten uitgebreid naar het projector augmented-wave (PAW) kader. Dit maakt de rigoureuze berekening van matrixelementen tussen niet-orthogonale Slater-determinanten mogelijk.
• Implementatie: De methodologie is geïmplementeerd in de GPAW-software. Transitie-matrixelementen worden afgeleid binnen het PAW-formalisme, waarbij rekening wordt gehouden met de transformatie tussen gladde pseudo-golffuncties en alle-elektronen golffuncties. Het formalisme bevat correcties voor de overlapmatrix en één-elektronen-operatoren (specifiek de dipooloperator) binnen de augmentatie-sferen.
• Computationele details: Berekeningen werden uitgevoerd met verschillende uitwisselings-correlatie (xc) functionalen: de gegeneraliseerde gradiëntbenadering (PBE), de hybride functionaal (PBE0), en PBE met expliciete Perdew–Zunger zelf-interactie correctie (SIC) met globale schalingfactoren ($\alpha = 1$ en $\alpha = 1/2$). Zowel plane-wave (PW) als linear combination of atomic orbitals (LCAO) basissets (aug-cc-pVDZ en d-aug-cc-pVDZ) werden gebruikt binnen het PAW-kader. Resultaten werden vergeleken met hoogwaardige multireferentie referentiedata (MS-CASPT2, CCSD, CCSD(T) en CCSDT) en LR-TDDFT-berekeningen.
• Spin-zuivering: Voor open-shell singlet-toestanden werden energieën en TDM's verkregen met behulp van spin-zuiveringsformules afgeleid van mixed-spin en triplet-berekeningen.

Belangrijkste Bijdragen

1. Uitbreiding van Löwdin's regels naar PAW: Het artikel leidt Löwdin's regels voor niet-orthogonale determinanten af en implementeert deze specifiek binnen het PAW-formalisme, wat de berekening van transitie-eigenschappen voor OO-toestanden met plane-wave en grid-gebaseerde representaties mogelijk maakt.
2. Systematische Benchmarking: De studie biedt een uitgebreide beoordeling van OO-oscillatorsterktes voor meerdere aangeslagen toestanden (voorbij de HOMO–LUMO-transitie) over een diverse set moleculen, waarbij verschillende basissets en xc-functionalen worden vergeleken.
3. Analyse van Niet-Orthogonaliteit: Het werk onderzoekt expliciet of de niet-nul overlap tussen de grondtoestand en aangeslagen toestanden in OO-berekeningen systematische fouten in oscillatorsterktes introduceert, een punt van zorg die in eerdere literatuur is geuit maar niet volledig is opgelost voor volledige absorptiespectra.

Resultaten

• Basisset-afhankelijkheid: Oscillatorsterktes voor diffuse Rydberg-toestanden zijn gevoelig voor de basisset-representatie. Singel-geaugmenteerde LCAO-basissets (aug) overschatten vaak excitatie-energieën en leveren onnauwkeurige oscillatorsterktes voor zeer diffuse toestanden vergeleken met dubbel-geaugmenteerde (d-aug) en plane-wave (PW) representaties. Bright valence-transities zijn echter grotendeels ongevoelig voor de basisset.
• Functionaal-afhankelijkheid: Het wijzigen van de xc-functionaal (PBE, PBE0, PBE-SIC/2, PBE-SIC) verandert de excitatie-energieën en specifieke oscillatorsterktes, maar produceert geen systematische hiërarchie van nauwkeurigheid voor intensiteiten over de gehele dataset. Hoewel SIC de asymptotische potentiaal verbetert, corrigeert het de fouten in oscillatorsterkte niet consistent.
• Single- vs. Multi-Configurational Karakter: Er ontstaat een duidelijk onderscheid op basis van het elektronische karakter van de aangeslagen toestanden.
  • Single-Configurational Toestanden: Voor toestanden die worden gedomineerd door een enkel determinant (bijv. veel Rydberg-toestanden in ammoniak), leveren OO-berekeningen een uitstekende overeenkomst met referentie-oscillatorsterktes en piekposities.
  • Multi-Configurational Toestanden: Voor toestanden met een sterk multi-configurational karakter of significante configuratie-mixing (bijv. H₂O S3/S4, CH₂O S5, C₂H₄ S1), vertonen OO-berekeningen substantiële discrepanties. Specifiek faalt de single-determinant beschrijving om de correcte intensiteitsverdeling te vatten, wat leidt tot overschatting of onderschatting van de oscillatorsterktes die standhoudt over alle geteste functionalen.
• Niet-Orthogonaliteit en Overlap: Analyse laat zien dat er geen duidelijke correlatie bestaat tussen de grootte van de grondtoestand/aangeslagen toestand overlap en de fout in de oscillatorsterktes. Zo leidde het verminderen van de overlap via hybride functionalen of SIC niet systematisch tot een verbetering van de oscillatorsterkte voor problematische toestanden zoals H₂O S3. Daarom wordt niet-orthogonaliteit niet geïdentificeerd als de primaire bron van de fout bij de grootste afwijkingen.
• Vergelijking met LR-TDDFT: OO-methoden presteren over het algemeen beter dan LR-TDDFT bij het beschrijven van de energieën en intensiteiten van Rydberg-toestanden. Daarentegen biedt LR-TDDFT vaak een betere beschrijving van bright valence-transities met sterke configuratie-mixing, waar de single-determinant OO-benadering moeite mee heeft.

Betekenis en Conclusies
Het artikel concludeert dat hoewel OO-dichtheidsfunctionaalberekeningen een robuust instrument zijn voor het voorspellen van excitatie-energieën en absorptiespectra van single-configurational Rydberg-toestanden, hun kwantitatieve betrouwbaarheid voor oscillatorsterktes beperkt wordt door het single-determinant karakter van de Kohn–Sham beschrijving. De methode presteert goed wanneer de aangeslagen toestand adequaat wordt beschreven door één dominante configuratie, maar faalt voor toestanden met een sterk multi-configurational karakter, ongeacht de keuze van de xc-functionaal of basisset.

De auteurs stellen dat de grote fouten waargenomen in oscillatorsterktes voor specifieke transities (bijv. water S3/S4, formaldehyde S5, ethyleen S1) inherent zijn aan de single-determinant benadering en geen artefacten zijn van niet-orthogonaliteit of basisset-onvolledigheid. Bijgevolg moeten OO-KS spectra voor systemen waar multi-configurational karakter significant is, kwalitatief worden geïnterpreteerd (bijv. voor piekposities) in plaats van kwantitatief met betrekking tot intensiteiten. De studie suggereert dat toekomstige verbeteringen een uitbreiding van het OO-kader vereisen om expliciet multi-configurational aangeslagen toestanden te behandelen, in plaats van enkel te vertrouwen op functionaal- of basisset-modificaties.