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

Dieser Beitrag stellt die Herleitung, Implementierung und Benchmarking analytischer nichtadiabatischer Kopplungsmatrixelemente erster Ordnung für die spinangepasste X-TDDFT-Methode vor und zeigt, dass sie die Fehler im Vergleich zur Standard-U-TDDFT-Methode signifikant reduziert und qualitativ korrekte Einblicke in die Photophysik offenschaliger Systeme wie Kupfer(II)-Porphyrin liefert.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen vorherzusagen, wie eine komplexe Maschine, wie ein aus Magneten gefertigtes Kreisel, reagiert, wenn Sie sie schütteln. In der Welt der Chemie ist diese „Maschine" ein Molekül mit ungepaarten Elektronen (wie organische Radikale oder Übergangsmetallkomplexe), und das „Schütteln" ist Licht, das auf es trifft und es in einen angeregten Zustand versetzt.

Wissenschaftler verwenden ein Werkzeug namens zeitabhängige Dichtefunktionaltheorie (TDDFT), um diese Reaktionen zu simulieren. Betrachten Sie TDDFT als eine hochentwickelte Wettervorhersage für Moleküle. Sie sagt voraus, wie sich das Molekül bewegt und wie sich seine Energie verändert.

Es gibt jedoch ein Problem. Die Standard-TDDFT (nennen wir sie U-TDDFT) ist wie eine Wettervorhersage, die davon ausgeht, dass der Wind immer in einer geraden Linie weht. Sie funktioniert für einfache Moleküle einigermaßen, verwirrt sich aber bei komplexen Molekülen mit „ungepaarten" Elektronen (wie unserem magnetischen Kreisel). Sie behandelt die beiden „Spins" der Elektronen (nennen wir sie Spin A und Spin B) so, als wären sie unabhängig voneinander, was zu Fehlern führt. Es ist, als würde man versuchen, einen Tanz zu beschreiben, bei dem zwei Partner Händchen halten, die Vorhersage jedoch davon ausgeht, dass sie allein tanzen.

Die neue Lösung: X-TDDFT

Die Autoren dieses Papiers entwickelten ein Upgrade namens X-TDDFT. Dies ist wie ein neues Wettermodell, das versteht, dass die Partner Händchen halten. Es zwingt die Mathematik, die „Spin"-Regeln der Quantenmechanik zu respektieren. Sie hatten dies bereits verwendet, um die Energie und die Form dieser Moleküle besser vorherzusagen, aber ihnen fehlte ein entscheidendes Puzzleteil: Nichtadiabatische Kopplungsmatrixelemente (NACMEs).

Was ist ein NACME?
Stellen Sie sich vor, das Molekül ist ein Auto, das auf einer holprigen Straße fährt.

• Die Energie sagt Ihnen, wie schnell das Auto fährt.
• Die Gradienten sagen Ihnen, in welche Richtung die Straße abfällt.
• Die NACMEs sagen Ihnen, wie wahrscheinlich es ist, dass das Auto die Spur wechselt oder in einen anderen Zustand „crasht".

In der Chemie wird dieses „Spurwechseln" als Interne Konversion (IC) bezeichnet. Es ist der Prozess, bei dem ein Molekül Energie absorbiert, angeregt wird und diese Energie dann schnell wieder in den Grundzustand abgibt, oft unter Freisetzung von Wärme statt von Licht. Wenn Ihre NACME-Berechnung falsch ist, könnten Sie denken, das Auto bleibt in seiner Spur, während es in Wirklichkeit kurz davor ist, gewaltsam in einen Graben zu schleudern.

Was haben sie getan?

Das Team leitete die mathematischen Formeln her, um diese „Spurwechsel"-Wahrscheinlichkeiten (NACMEs) mit ihrer neuen, spinbewussten X-TDDFT-Methode zu berechnen. Anschließend testeten sie dies auf zwei Arten:

1. Der kleine Test (Formaldehyd-Radikal): Sie verglichen ihre neue Methode mit einem „Goldstandard" einer supersicheren Berechnung (wie den Abgleich eines neuen GPS mit einer Satellitenkarte). Sie stellten fest, dass die alte Methode (U-TDDFT) oft um einen riesigen Betrag falsch lag – manchmal um ein Drittel bis zwei Drittel. Die neue Methode (X-TDDFT) korrigierte die meisten dieser Fehler und machte die Vorhersage, wie schnell das Molekül „abkühlt" (Interne-Konversions-Rate), viel genauer. In einigen Fällen sagte die neue Methode voraus, dass die Abkühlgeschwindigkeit 100-mal langsamer war als von der alten Methode vorhergesagt.

2. Der große Test (Kupfer-Porphyrine): Sie untersuchten komplexe kupferbasierte Moleküle (ähnlich dem Häm im Blut, aber mit Kupfer).

   • Die alte Sichtweise (U-TDDFT): Sagte voraus, dass das Molekül, wenn es angeregt wird, eine gleich große Chance hat, direkt abzukühlen oder einen Umweg durch Zwischenzustände zu nehmen.
   • Die neue Sichtweise (X-TDDFT): Sagte voraus, dass das Molekül fast nie direkt abkühlt. Es nimmt fast immer den Umweg.
   • Das Ergebnis: Dies veränderte die Geschichte des Verhaltens dieser Moleküle vollständig. Die alte Methode hatte nicht nur die Zahlen leicht falsch berechnet; sie hatte die Geschichte falsch erzählt. Sie verzerrte auch den Vergleich zwischen verschiedenen Versionen des Moleküls (mit unterschiedlichen chemischen Verzierungen) und ließ eine Version schneller erscheinen als eine andere, während das Gegenteil der Fall war.

Das Fazit

Das Papier kommt zu dem Schluss, dass man für Moleküle mit ungepaarten Elektronen (wie Radikale oder Übergangsmetalle) die alte „geradlinige" Mathematik (U-TDDFT) nicht vertrauen kann, um vorherzusagen, wie sie ihre Energiezustände wechseln.

Genauso wie man keine flache Karte zur Navigation in einem Gebirge verwenden würde, sollte man die alte TDDFT-Methode nicht für diese komplexen Moleküle verwenden. Die neue X-TDDFT-Methode wirkt wie eine 3D-topografische Karte und zeigt, dass die „Straßen" (Energiepfade) sehr anders sind als bisher angenommen. Dies ist entscheidend für Wissenschaftler, die bessere Solarzellen, LEDs oder Katalysatoren entwickeln wollen, denn wenn man nicht weiß, in welche „Spur" das Molekül springen wird, kann man sein Verhalten nicht steuern.

Kurz gesagt: Die Autoren haben einen besseren Lineal gebaut, um zu messen, wie Moleküle zwischen Energiezuständen „springen". Sie bewiesen, dass das alte Lineal so ungenau war, dass es völlig unterschiedliche Geschichten darüber erzählte, wie diese Moleküle funktionieren, insbesondere bei denen, die Kupfer und andere Übergangsmetalle beinhalten.

Technische Zusammenfassung

Technische Zusammenfassung: Analytische Matrixelemente der nichtadiabatischen Kopplung erster Ordnung für spinangepasste Open-Shell-Zeitabhängige Dichtefunktionaltheorie

Problemstellung
Die spinangepasste zeitabhängige Dichtefunktionaltheorie (X-TDDFT) hat sich kürzlich als robuste Methode zur Berechnung von Anregungsenergien und analytischen Gradienten für Open-Shell-Moleküle (wie Radikale und Übergangsmetallkomplexe) etabliert. Sie bietet signifikante Verbesserungen gegenüber der unrestricted TDDFT (U-TDDFT), indem sie Spin-Kontamination korrekt behandelt und spezifische Doppelanregungen implizit einbezieht. Der Einfluss der Spin-Anpassung auf Matrixelemente der nichtadiabatischen Kopplung (NACMEs) – die für die Simulation von Raten der internen Konversion (IC) und nichtadiabatischer Dynamik entscheidend sind – blieb jedoch unbekannt. Zwar existieren rigorose U-TDDFT-NACME-Theorien, doch leiden diese unter großen Fehlern infolge von Spin-Kontamination, insbesondere für Zustände, die durch Anregungen vom besetzten zum unbesetzten Orbital (CV) dominiert werden. Bestehende spinangepasste NACME-Theorien beschränken sich entweder auf semiempirische Hamilton-Operatoren, spezifische Spin-Flip-Anregungen oder fehlen der notwendigen Formalismus für allgemeine Open-Shell-Systeme. Folglich wurde die genaue Simulation von Photophysik und Photochemie in Open-Shell-Systemen unter Verwendung spin-erhaltender Rahmenwerke durch das Fehlen einer konsistenten X-TDDFT-NACME-Theorie behindert.

Methodik
Die Autoren leiten und implementieren analytische NACMEs erster Ordnung für die X-TDDFT-Methode ab, speziell für Übergänge vom Grundzustand zum angeregten Zustand (ge) und zwischen angeregten Zuständen (ee). Die Herleitung folgt einem Ansatz der Bewegungsgleichung (EOM), bei dem die Standard-Anregungs- und De-Anregungsoperatoren der Random-Phase-Näherung (RPA) durch ihre spinangepassten Gegenstücke ersetzt werden.

Zu den wichtigsten theoretischen Komponenten gehören:

1. Spin-Angepasste Operatoren: Der Formalismus nutzt einen restricted open-shell Kohn-Sham (ROKS)-Referenzzustand. Der Anregungsoperator $O^\dagger_I$ wird unter Verwendung von spinangepassten CV(0)- und CV(1)-Anregungen konstruiert, wobei letztere mit spezifischen Doppelanregungen ($\Psi^{\bar{t}a}_{\bar{i}t}$) gemischt werden, um eine korrekte Spinsymmetrie sicherzustellen.
2. Lagrange-Formulierung: Um die implizite Abhängigkeit der Anregungsvektoren und der Molekülorbital-Koeffizienten (MO) von den Kernkoordinaten zu behandeln, wird ein Lagrange-Ansatz verwendet. Dies ermöglicht die Berechnung analytischer Ableitungen ohne explizite Differentiation der MO-Koeffizienten, wobei Pulay-Terme über Lagrange-Multiplikatoren (Z- und W-Vektoren) einbezogen werden.
3. Korrekturterme: Die X-TDDFT-NACMEs werden als die Standard-U-TDDFT-NACMEs plus spezifische Korrekturterme ausgedrückt, die sich ergeben aus:
   • Der Verwendung eines spinangepassten Referenzsystems (ROKS) anstelle eines unrestricted (UKS).
   • Der Modifikation des CV(1)-Anregungsoperators zur Einbeziehung von Doppelanregungen.
   • Den daraus resultierenden Änderungen in den Übergangsdichtematrizen (TDMs) und der Casida-Matrix (A- und B-Blöcke).
4. Implementierungsstrategie: Die Implementierung ist so konzipiert, dass sie auf bestehenden U-TDDFT-NACME-Codes aufbaut. Die Autoren zeigen, dass der zusätzliche Rechenaufwand akzeptabel ist und mit dem analytischer Gradientenberechnungen vergleichbar ist, da er hauptsächlich das Lösen modifizierter Z-Vektor-Gleichungen und die Berechnung intermediärer Dichtematrizen umfasst.

Hauptbeiträge

• Theoretische Herleitung: Der Artikel liefert die erste rigorose Herleitung analytischer NACMEs für spin-erhaltende X-TDDFT, die sowohl ge- als auch ee-Übergänge abdeckt.
• Algorithmische Implementierung: Die Methode wurde im BDF-Softwarepaket implementiert. Sie nutzt die bestehende U-TDDFT-Infrastruktur und erfordert lediglich das Hinzufügen spezifischer Korrekturterme zu den Übergangsdichtematrizen und den Lagrange-Multiplikatoren.
• Benchmarking: Die Autoren testen die Methode an hochrangigen Multireferenz-Rechnungen (XMS-CASPT2) für das Formaldehyd-Radikalkation ($CH_2O^+$).
• Anwendung auf Übergangsmetallkomplexe: Die Methode wird auf Kupfer(II)-Porphyrine (CuP, CuOEP, CuTPP) angewendet, um Relaxationspfade angeregter Zustände und Raten der internen Konversion zu simulieren.

Ergebnisse

• Benchmark-Genauigkeit ($CH_2O^+$): Im Vergleich zu XMS-CASPT2-Referenzdaten reduzieren X-TDDFT-NACMEs die mittlere absolute Abweichung (MAD) der Normen von U-TDDFT-NACMEs um etwa 32 % bis 65 % (abhängig vom Funktional). Die signifikantesten Verbesserungen werden für Zustände mit ausgeprägtem CV(1)-Charakter beobachtet, bei denen U-TDDFT aufgrund von Spin-Kontamination oft versagt. Die Autoren weisen jedoch darauf hin, dass X-TDDFT für bestimmte Mischzustände (z. B. $1^2A_2$) gelegentlich größere Fehler als U-TDDFT liefern kann; dies wird jedoch auf zufällige Fehlerkompensation in den U-TDDFT-Ergebnissen zurückgeführt und nicht auf ein Versagen der spinangepassten Theorie.
• Kupfer(II)-Porphyrine:
  • NACME-Magnituden: X-TDDFT sagt für den Übergang $2T_1 \to 2S_0$ signifikant kleinere NACME-Magnituden (bis zu einer Größenordnung) voraus als U-TDDFT. Diese Reduktion wird auf die nahezu perfekte Auslöschung von $\alpha$- und $\beta$-Beiträgen in reinen CV(1)-Zuständen innerhalb des spinangepassten Rahmens zurückgeführt, eine Auslöschung, die in U-TDDFT aufgrund von Unterschieden in der Orbitalform unvollständig ist.
  • Raten der internen Konversion: Die Reduktion der NACME-Normen führt zu drastischen Verringerungen der berechneten IC-Raten (bis zu zwei Größenordnungen).
  • Relaxationspfade: Für unsubstituiertes CuP sagt U-TDDFT voraus, dass der dominante Relaxationspfad etwa gleich große direkte IC zum Grundzustand und indirekte IC über Ligandenfeldzustände umfasst. Im Gegensatz dazu sagt X-TDDFT voraus, dass der indirekte Pfad über Ligandenfeldzustände dominiert.
  • Substituenteneffekte: Die relativen Trends der IC-Raten für substituierte Porphyrine (CuP, CuOEP, CuTPP) unterscheiden sich zwischen den beiden Methoden. U-TDDFT sagt einen Trend von CuOEP < CuTPP < CuP voraus, während X-TDDFT CuTPP < CuOEP < CuP vorhersagt. Dies zeigt, dass U-TDDFT-Fehler nicht nur groß, sondern auch unsystematisch sind und selbst für eng verwandte Moleküle relative Trends nicht erfassen.

Bedeutung und Behauptungen
Der Artikel behauptet, dass die Entwicklung von X-TDDFT-NACMEs das wesentliche Werkzeugset zur Untersuchung photophysikalischer und photochemischer Prozesse von Open-Shell-Molekülen unter Verwendung spinangepasster TDDFT vervollständigt. Die Autoren betonen, dass Spin-Anpassung nicht nur eine Korrektur für Energien und Gradienten ist, sondern für NACMEs und folglich für die Genauigkeit simulierter nichtadiabatischer Dynamik und IC-Raten entscheidend ist.

Die Studie kommt zu dem Schluss, dass die reliance auf U-TDDFT-NACMEs für Open-Shell-Systeme zu qualitativ falschen Vorhersagen von Relaxationsmechanismen und zu quantitativen Fehlern in Geschwindigkeitskonstanten führen kann, die sowohl groß als auch unsystematisch sind. Daher empfehlen die Autoren die Verwendung von X-TDDFT-NACMEs für Open-Shell-Systeme im Allgemeinen, selbst in Fällen, in denen U-TDDFT- und X-TDDFT-Anregungsenergien ähnlich erscheinen. Die Arbeit legt den Grundstein für zukünftige Anwendungen in der nichtadiabatischen Molekulardynamik und der Untersuchung ultraschneller Photophysik in Radikalen und Übergangsmetallkomplexen.