Analytic first order nonadiabatic coupling matrix elements of spin-adapted open-shell time-dependent density functional theory

Dit artikel presenteert de afleiding, implementatie en benchmarking van analytische eerste-orde niet-adiabatische koppelingsmatrixelementen voor de spin-aangepaste X-TDDFT-methode, en toont aan dat deze methode de fouten aanzienlijk vermindert in vergelijking met standaard U-TDDFT en kwalitatief correcte inzichten biedt in de fotofysica van open-schil-systemen zoals koper(II)porfyrine.

De kern

Stel je voor dat je probeert te voorspellen hoe een complexe machine, zoals een tol van magneten die draait, zal reageren wanneer je hem schudt. In de wereld van de chemie is deze "machine" een molecuul met ongepaarde elektronen (zoals organische radicalen of overgangsmetaalcomplexen), en de "schok" is licht dat erop valt, waardoor het in een aangeslagen toestand terechtkomt.

Wetenschappers gebruiken een hulpmiddel genaamd Tijd-afhankelijke Dichtheidsfunctionaaltheorie (TDDFT) om deze reacties te simuleren. Denk aan TDDFT als een geavanceerde weersvoorspelling voor moleculen. Het voorspelt hoe het molecuul beweegt en energie verandert.

Er is echter een probleem. Standaard TDDFT (laten we het U-TDDFT noemen) is als een weersvoorspelling die ervan uitgaat dat de wind altijd in een rechte lijn waait. Het werkt redelijk voor eenvoudige moleculen, maar voor complexe moleculen met "ongepaarde" elektronen (zoals onze tol van magneten) raakt het in de war. Het behandelt de twee "spins" van de elektronen (laten we ze Spin A en Spin B noemen) alsof ze onafhankelijk zijn, wat leidt tot fouten. Het is alsof je probeert een dans te beschrijven waarbij twee partners hand in hand houden, maar de voorspelling ervan uitgaat dat ze alleen dansen.

De Nieuwe Oplossing: X-TDDFT

De auteurs van dit artikel hebben een upgrade ontwikkeld genaamd X-TDDFT. Dit is als een nieuw weermodel dat begrijpt dat de partners hand in hand houden. Het dwingt de wiskunde om de "spin"-regels van de kwantummechanica te respecteren. Ze hadden dit al gebruikt om de energie en vorm van deze moleculen beter te voorspellen, maar ze misten één cruciaal stukje: Niet-adiabatische Koppelingsmatrix-elementen (NACME's).

Wat is een NACME?
Stel je voor dat het molecuul een auto is die over een hobbelige weg rijdt.

• Energie vertelt je hoe snel de auto gaat.
• Gradienten vertellen je welke kant de weg afloopt.
• NACME's vertellen je hoe waarschijnlijk het is dat de auto van baan wisselt of in een andere toestand crasht.

In de chemie heet dit "balken wisselen" Interne Conversie (IC). Het is het proces waarbij een molecuul energie absorbeert, aangeslagen raakt en die energie vervolgens snel weer terug naar de grondtoestand afvoert, vaak met warmte in plaats van licht. Als je NACME-berekening verkeerd is, denk je misschien dat de auto in zijn baan blijft, terwijl hij in werkelijkheid op het punt staat heftig de greppel in te sturen.

Wat hebben ze gedaan?

Het team heeft de wiskundige formules afgeleid om deze "baanwissel"-kansen (NACME's) te berekenen met hun nieuwe, spin-bewuste X-TDDFT-methode. Vervolgens hebben ze het op twee manieren getest:

1. De Kleine Test (Formaldehyde-radicaal): Ze hebben hun nieuwe methode vergeleken met een "gouden standaard" van super-accurate berekening (alsof je een nieuwe GPS vergelijkt met een satellietkaart). Ze ontdekten dat de oude methode (U-TDDFT) vaak met een enorme marge verkeerd zat – soms wel een derde tot twee derde afwijkend. De nieuwe methode (X-TDDFT) corrigeerde de meeste van deze fouten, waardoor de voorspelling van hoe snel het molecuul "afkoelt" (Interne Conversie-snelheid) veel nauwkeuriger werd. In sommige gevallen voorspelde de nieuwe methode dat de afkoelsnelheid 100 keer langzamer was dan de oude methode voorspelde.

2. De Grote Test (Koperporfyrines): Ze keken naar complexe koperbevattende moleculen (vergelijkbaar met het heme in bloed, maar dan met koper).

   • Het Oude Beeld (U-TDDFT): Voorspelde dat wanneer het molecuul aangeslagen raakt, het evenveel kans heeft om direct af te koelen of een omweg te nemen via tussenliggende toestanden.
   • Het Nieuwe Beeld (X-TDDFT): Voorspelde dat het molecuul bijna nooit direct afkoelt. Het neemt bijna altijd de omweg.
   • Het Resultaat: Dit veranderde het verhaal over hoe deze moleculen zich gedragen volledig. De oude methode zat niet alleen een beetje fout in de cijfers; ze had het verhaal verkeerd. Het verpestte ook de vergelijking tussen verschillende versies van het molecuul (met verschillende chemische decoraties), waardoor het leek alsof één versie sneller was dan een andere, terwijl het omgekeerde waar was.

De Conclusie

Het artikel concludeert dat voor moleculen met ongepaarde elektronen (zoals radicalen of overgangsmetalen) je de oude "rechte-lijn"-wiskunde (U-TDDFT) niet kunt vertrouwen voor het voorspellen van hoe ze van energietoestand wisselen.

Zoals je geen platte kaart zou gebruiken om een bergketen te navigeren, moet je de oude TDDFT-methode niet gebruiken voor deze complexe moleculen. De nieuwe X-TDDFT-methode fungeert als een 3D-topografische kaart, die onthult dat de "wegen" (energieroutes) heel anders zijn dan eerder gedacht. Dit is cruciaal voor wetenschappers die proberen betere zonnecellen, LED's of katalysatoren te ontwerpen, want als je niet weet in welke "baan" het molecuul zal springen, kun je zijn gedrag niet controleren.

Kortom: De auteurs bouwden een betere liniaal om te meten hoe moleculen "springen" tussen energietoestanden. Ze bewezen dat de oude liniaal zo onnauwkeurig was dat hij compleet verschillende verhalen vertelde over hoe deze moleculen werken, vooral voor diegene die koper en andere overgangsmetalen bevatten.

Technische samenvatting

Technische Samenvatting: Analytische Eerste-orde Nonadiabatische Koppelingsmatrixelementen van Spin-geadapteerde Open-schil Tijd-afhankelijke Dichtheidsfunctionaaltheorie

Probleemstelling
Spin-geadapteerde tijd-afhankelijke dichtheidsfunctionaaltheorie (X-TDDFT) is recentelijk gevestigd als een robuuste methode voor het berekenen van excitatie-energieën en analytische gradiënten van open-schil moleculen (zoals radicalen en overgangsmetaalcomplexen). Deze methode biedt aanzienlijke verbeteringen ten opzichte van onbeperkte TDDFT (U-TDDFT) door spinverontreiniging correct te behandelen en specifiek dubbele excitaties impliciet op te nemen. Het effect van spin-adaptatie op nonadiabatische koppelingsmatrixelementen (NACMEs) – die cruciaal zijn voor het simuleren van interne conversie (IC) snelheden en nonadiabatische dynamica – bleef echter onbekend. Hoewel er rigoureuze U-TDDFT NACME-theorieën bestaan, lijden deze onder grote fouten door spinverontreiniging, met name voor toestanden die worden gedomineerd door excitaties van een gesloten schil naar een lege schil (CV). Bestaande spin-geadapteerde NACME-theorieën zijn beperkt tot semi-empirische Hamiltonianen, specifieke spin-flip excitaties, of missen de nodige formalisme voor algemene open-schil systemen. Bijgevolg werd de accurate simulatie van fotofysica en fotochemie in open-schil systemen met behulp van spinbehoudende kaders gehinderd door het ontbreken van een consistente X-TDDFT NACME-theorie.

Methodologie
De auteurs leiden analytische eerste-orde NACMEs af en implementeren deze voor de X-TDDFT-methode, specifiek voor overgangen van de grondtoestand naar een aangeslagen toestand (ge) en van aangeslagen toestand naar aangeslagen toestand (ee). De afleiding volgt een equation-of-motion (EOM) aanpak, waarbij de standaard random phase approximation (RPA) excitatie- en de-excitatie-operatoren worden vervangen door hun spin-geadapteerde tegenhangers.

Belangrijke theoretische componenten omvatten:

1. Spin-geadapteerde Operatoren: Het formalisme maakt gebruik van een beperkte open-schil Kohn-Sham (ROKS) referentietoestand. De excitatie-operator $O^\dagger_I$ wordt geconstrueerd met behulp van spin-geadapteerde CV(0) en CV(1) excitaties, waarbij de laatste wordt gemengd met specifieke dubbele excitaties ($\Psi^{\bar{t}a}_{\bar{i}t}$) om de juiste spinsymmetrie te waarborgen.
2. Lagrangiaan Formulering: Om de impliciete afhankelijkheid van excitatievectoren en moleculaire orbitaal (MO) coëfficiënten van nucleaire coördinaten te behandelen, wordt een Lagrangiaan-aanpak toegepast. Dit maakt de berekening van analytische afgeleiden mogelijk zonder expliciet de MO-coëfficiënten te differentiëren, waarbij Pulay-termen worden opgenomen via Lagrange-multiplicatoren ($Z$- en $W$-vectoren).
3. Correctietermen: De X-TDDFT NACMEs worden uitgedrukt als de standaard U-TDDFT NACMEs plus specifieke correctietermen die voortvloeien uit:
   • Het gebruik van een spin-geadapteerde referentie (ROKS) in plaats van een onbeperkte (UKS).
   • De modificatie van de CV(1) excitatie-operator om dubbele excitaties op te nemen.
   • De resulterende veranderingen in de overgangsdichtheidsmatrices (TDMs) en de Casida-matrix (de $A$- en $B$-blokken).
4. Implementatiestrategie: De implementatie is ontworpen om te worden gebouwd op bestaande U-TDDFT NACME-codes. De auteurs tonen aan dat de extra rekenkosten aanvaardbaar zijn en vergelijkbaar met die van analytische gradiëntberekeningen, aangezien dit voornamelijk het oplossen van gewijzigde Z-vector-vergelijkingen en het berekenen van intermediaire dichtheidsmatrices omvat.

Belangrijkste Bijdragen

• Theoretische Afleiding: Het artikel biedt de eerste rigoureuze afleiding van analytische NACMEs voor spinbehoudende X-TDDFT, die zowel ge- als ee-overgangen bestrijken.
• Algoritmische Implementatie: De methode is geïmplementeerd in de BDF-softwarepakket. Deze maakt gebruik van bestaande U-TDDFT-infrastructuur, waarbij alleen de toevoeging van specifieke correctietermen aan de overgangsdichtheidsmatrices en de Lagrange-multiplicatoren vereist is.
• Benchmarking: De auteurs benchmarken de methode tegenover hoogwaardige multireference-berekeningen (XMS-CASPT2) voor het formaldehyde-radicaal-kation ($CH_2O^+$).
• Toepassing op Overgangsmetaalcomplexen: De methode wordt toegepast op koper(II) porfyrines (CuP, CuOEP, CuTPP) om relaxatiepaden in aangeslagen toestanden en interne conversiesnelheden te simuleren.

Resultaten

• Benchmarknauwkeurigheid ($CH_2O^+$): In vergelijking met XMS-CASPT2-referentiedata verminderen X-TDDFT NACMEs de gemiddelde absolute afwijking (MAD) van de normen van U-TDDFT NACMEs met ongeveer 32% tot 65% (afhankelijk van de functionaal). De meest significante verbeteringen worden waargenomen voor toestanden met een aanzienlijk CV(1)-karakter, waarbij U-TDDFT vaak faalt door spinverontreiniging. De auteurs merken echter op dat X-TDDFT voor bepaalde gemengde toestanden (bijv. $1^2A_2$) af en toe grotere fouten kan opleveren dan U-TDDFT, hoewel dit wordt toegeschreven aan toevallige foutcompensatie in de U-TDDFT-resultaten en niet aan een falen van de spin-geadapteerde theorie.
• Koper(II) Porfyrines:
  • NACME-magnitudes: X-TDDFT voorspelt aanzienlijk kleinere NACME-magnitudes (tot een orde van grootte) voor de $2T_1 \to 2S_0$ overgang in vergelijking met U-TDDFT. Deze reductie wordt toegeschreven aan de bijna perfecte compensatie van $\alpha$- en $\beta$-bijdragen in pure CV(1)-toestanden binnen het spin-geadapteerde kader, een compensatie die in U-TDDFT onvolledig is door verschillen in orbitaalvorm.
  • Interne Conversiesnelheden: De reductie in NACME-normen leidt tot drastische reducties in de berekende IC-snelheden (tot twee ordes van grootte).
  • Relaxatiepaden: Voor ongesubstitueerd CuP voorspelt U-TDDFT dat het dominante relaxatiepad ongeveer gelijke directe IC naar de grondtoestand en indirecte IC via ligandveldtoestanden omvat. In tegenstelling hiermee voorspelt X-TDDFT dat het indirecte pad via ligandveldtoestanden dominant is.
  • Substituent-effecten: De relatieve trends van IC-snelheden voor gesubstitueerde porfyrines (CuP, CuOEP, CuTPP) verschillen tussen de twee methoden. U-TDDFT voorspelt een trend van CuOEP < CuTPP < CuP, terwijl X-TDDFT CuTPP < CuOEP < CuP voorspelt. Dit geeft aan dat U-TDDFT-fouten niet alleen groot zijn, maar ook onsystematisch, waardoor ze zelfs voor nauw verwante moleculen geen relatieve trends kunnen vastleggen.

Betekenis en Claims
Het artikel claimt dat de ontwikkeling van X-TDDFT NACMEs de essentiële toolkit vervolledigt voor het bestuderen van fotofysische en fotochemische processen van open-schil moleculen met behulp van spin-geadapteerde TDDFT. De auteurs benadrukken dat spin-adaptatie niet slechts een correctie is voor energieën en gradiënten, maar cruciaal is voor NACMEs en bijgevolg voor de nauwkeurigheid van gesimuleerde nonadiabatische dynamica en IC-snelheden.

De studie concludeert dat het vertrouwen op U-TDDFT NACMEs voor open-schil systemen kan leiden tot kwalitatief onjuiste voorspellingen van relaxatiemechanismen en kwantitatieve fouten in snelheidsconstanten die zowel groot als onsystematisch zijn. Daarom raden de auteurs het gebruik van X-TDDFT NACMEs voor open-schil systemen in het algemeen aan, zelfs in gevallen waar U-TDDFT- en X-TDDFT-excitatie-energieën vergelijkbaar lijken. Het werk legt de basis voor toekomstige toepassingen in nonadiabatische moleculaire dynamica en het bestuderen van ultrafast fotofysica in radicalen en overgangsmetaalcomplexen.