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

Cet article présente la dérivation, l'implémentation et l'évaluation des éléments de matrice de couplage non adiabatique analytiques du premier ordre pour la méthode X-TDDFT adaptée au spin, démontrant qu'elle réduit considérablement les erreurs par rapport à la U-TDDFT standard et fournit des insights qualitativement corrects sur la photophysique des systèmes à couche ouverte tels que la porphyrine de cuivre(II).

L'essentiel

Imaginez que vous essayez de prédire comment une machine complexe, comme un toupie constituée d'aimants, réagira lorsque vous la secouerez. Dans le monde de la chimie, cette « machine » est une molécule possédant des électrons non appariés (comme les radicaux organiques ou les complexes de métaux de transition), et le « secoussement » est la lumière qui la frappe, provoquant son saut vers un état excité.

Les scientifiques utilisent un outil appelé Théorie de la Fonctionnelle de la Densité Dépendante du Temps (TDDFT) pour simuler ces réactions. Considérez la TDDFT comme une prévision météorologique sophistiquée pour les molécules. Elle prédit comment la molécule se déplace et change d'énergie.

Cependant, il y a un problème. La TDDFT standard (appelons-la U-TDDFT) est comme une prévision météorologique qui suppose que le vent souffle toujours en ligne droite. Elle fonctionne correctement pour les molécules simples, mais pour les molécules complexes avec des électrons « non appariés » (comme notre toupie magnétique), elle se trompe. Elle traite les deux « spins » des électrons (appelons-les Spin A et Spin B) comme s'ils étaient indépendants, ce qui conduit à des erreurs. C'est comme essayer de décrire une danse où deux partenaires se tiennent par la main, alors que la prévision suppose qu'ils dansent seuls.

La Nouvelle Solution : X-TDDFT

Les auteurs de cet article ont développé une mise à niveau appelée X-TDDFT. C'est comme un nouveau modèle météorologique qui comprend que les partenaires se tiennent par la main. Il force les mathématiques à respecter les règles du « spin » de la mécanique quantique. Ils l'avaient déjà utilisée pour mieux prédire l'énergie et la forme de ces molécules, mais il leur manquait une pièce cruciale : les Éléments de Matrice de Couplage Non Adiabatique (EMCNA).

Qu'est-ce qu'un EMCNA ?
Imaginez que la molécule est une voiture roulant sur une route cahoteuse.

• L'Énergie vous indique à quelle vitesse la voiture va.
• Les Gradients vous indiquent dans quelle direction la route s'incline.
• Les EMCNA vous indiquent à quel point il est probable que la voiture change de voie ou crase dans un état différent.

En chimie, ce « changement de voie » s'appelle la Conversion Interne (CI). C'est le processus par lequel une molécule absorbe de l'énergie, devient excitée, puis rejette rapidement cette énergie vers l'état fondamental, libérant souvent de la chaleur au lieu de lumière. Si votre calcul d'EMCNA est erroné, vous pourriez penser que la voiture restera dans sa voie, alors qu'en réalité, elle est sur le point de virer violemment dans un fossé.

Que Ont-ils Fait ?

L'équipe a dérivé les formules mathématiques pour calculer ces probabilités de « changement de voie » (EMCNA) en utilisant leur nouvelle méthode X-TDDFT, sensible au spin. Ils l'ont ensuite testée de deux manières :

1. Le Petit Test (Radical Formaldéhyde) : Ils ont comparé leur nouvelle méthode à un calcul de référence « or » ultra-précis (comme vérifier un nouveau GPS par rapport à une carte satellite). Ils ont constaté que l'ancienne méthode (U-TDDFT) était souvent erronée de manière massive — parfois décalée d'un tiers à deux tiers. La nouvelle méthode (X-TDDFT) a corrigé la plupart de ces erreurs, rendant la prédiction de la vitesse à laquelle la molécule « refroidit » (taux de Conversion Interne) beaucoup plus précise. Dans certains cas, la nouvelle méthode a prédit que la vitesse de refroidissement était 100 fois plus lente que ce que l'ancienne méthode prédisait.

2. Le Grand Test (Porphyrines de Cuivre) : Ils ont examiné des molécules complexes à base de cuivre (similaires à l'hème dans le sang, mais avec du cuivre).

   • L'Ancienne Vue (U-TDDFT) : Prédit que lorsque la molécule est excitée, elle a une chance égale de refroidir directement ou de faire un détour par des états intermédiaires.
   • La Nouvelle Vue (X-TDDFT) : Prédit que la molécule presque jamais ne refroidit directement. Elle prend presque toujours le détour.
   • Le Résultat : Cela a complètement changé l'histoire du comportement de ces molécules. L'ancienne méthode ne s'est pas seulement trompée légèrement sur les chiffres ; elle a raconté la histoire de travers. Elle a également faussé la comparaison entre différentes versions de la molécule (avec différentes décorations chimiques), donnant l'impression qu'une version était plus rapide qu'une autre alors que l'inverse était vrai.

L'Essentiel

L'article conclut que pour les molécules avec des électrons non appariés (comme les radicaux ou les métaux de transition), on ne peut pas faire confiance aux anciennes mathématiques « en ligne droite » (U-TDDFT) pour prédire comment elles changent d'états d'énergie.

Tout comme vous n'utiliseriez pas une carte plate pour naviguer dans une chaîne de montagnes, vous ne devriez pas utiliser l'ancienne méthode TDDFT pour ces molécules complexes. La nouvelle méthode X-TDDFT agit comme une carte topographique 3D, révélant que les « routes » (voies d'énergie) sont très différentes de ce que l'on pensait auparavant. C'est crucial pour les scientifiques essayant de concevoir de meilleures cellules solaires, DEL ou catalyseurs, car si vous ne savez pas dans quelle « voie » la molécule va sauter, vous ne pouvez pas contrôler son comportement.

En bref : Les auteurs ont construit une meilleure règle pour mesurer comment les molécules « sautent » entre les états d'énergie. Ils ont prouvé que l'ancienne règle était si imprécise qu'elle racontait des histoires complètement différentes sur le fonctionnement de ces molécules, en particulier pour celles impliquant le cuivre et d'autres métaux de transition.

Résumé technique

Résumé technique : Éléments de matrice de couplage non adiabatique analytiques d'ordre premier de la théorie de la fonctionnelle de la densité dépendante du temps adaptée aux spins pour les systèmes à couche ouverte

Énoncé du problème
La théorie de la fonctionnelle de la densité dépendante du temps adaptée aux spins (X-TDDFT) a récemment été établie comme une méthode robuste pour le calcul des énergies d'excitation et des gradients analytiques des molécules à couche ouverte (telles que les radicaux et les complexes de métaux de transition), offrant des améliorations significatives par rapport à la TDDFT non contrainte (U-TDDFT) en traitant correctement la contamination de spin et en incluant implicitement des excitations doubles spécifiques. Cependant, l'effet de l'adaptation aux spins sur les éléments de matrice de couplage non adiabatique (NACMEs) — critiques pour la simulation des taux de conversion interne (IC) et de la dynamique non adiabatique — restait inconnu. Bien que des théories rigoureuses de NACMEs U-TDDFT existent, elles souffrent d'erreurs importantes dues à la contamination de spin, en particulier pour les états dominés par des excitations de couche fermée vers couche vacante (CV). Les théories de NACMEs adaptées aux spins existantes sont soit limitées à des hamiltoniens semi-empiriques, soit à des excitations spécifiques de retournement de spin, soit elles manquent du formalisme nécessaire pour les systèmes à couche ouverte généraux. Par conséquent, la simulation précise de la photophysique et de la photochimie dans les systèmes à couche ouverte utilisant des cadres conservant le spin était entravée par l'absence d'une théorie cohérente de NACMEs X-TDDFT.

Méthodologie
Les auteurs dérivent et implémentent des NACMEs analytiques d'ordre premier pour la méthode X-TDDFT, spécifiquement pour les transitions de l'état fondamental vers l'état excité (ge) et d'état excité vers état excité (ee). La dérivation suit une approche de mouvement d'équation (EOM), remplaçant les opérateurs d'excitation et de désexcitation standards de l'approximation de phase aléatoire (RPA) par leurs équivalents adaptés aux spins.

Les composants théoriques clés comprennent :

1. Opérateurs adaptés aux spins : Le formalisme utilise un état de référence Kohn-Sham à couche ouverte restreinte (ROKS). L'opérateur d'excitation $O^\dagger_I$ est construit en utilisant des excitations CV(0) et CV(1) adaptées aux spins, où cette dernière est mélangée avec des excitations doubles spécifiques ($\Psi^{\bar{t}a}_{\bar{i}t}$) pour assurer une symétrie de spin appropriée.
2. Formulation lagrangienne : Pour gérer la dépendance implicite des vecteurs d'excitation et des coefficients d'orbitales moléculaires (OM) par rapport aux coordonnées nucléaires, une approche lagrangienne est employée. Cela permet le calcul de dérivées analytiques sans différencier explicitement les coefficients d'OM, en incorporant les termes de Pulay via des multiplicateurs de Lagrange (vecteurs $Z$ et $W$).
3. Termes de correction : Les NACMEs X-TDDFT sont exprimés comme les NACMEs U-TDDFT standards plus des termes de correction spécifiques découlant de :
   • L'utilisation d'une référence adaptée aux spins (ROKS) au lieu d'une référence non contrainte (UKS).
   • La modification de l'opérateur d'excitation CV(1) pour inclure des excitations doubles.
   • Les changements résultants dans les matrices de densité de transition (TDM) et la matrice de Casida (blocs $A$ et $B$).
4. Stratégie d'implémentation : L'implémentation est conçue pour être construite sur la base de codes NACMEs U-TDDFT existants. Les auteurs démontrent que la surcharge de calcul supplémentaire est acceptable, comparable à celle des calculs de gradients analytiques, car elle implique principalement la résolution d'équations de vecteur Z modifiées et le calcul de matrices de densité intermédiaires.

Contributions clés

• Dérivation théorique : L'article fournit la première dérivation rigoureuse de NACMEs analytiques pour la X-TDDFT conservant le spin, couvrant à la fois les transitions ge et ee.
• Implémentation algorithmique : La méthode est implémentée dans le package logiciel BDF. Elle exploite l'infrastructure existante U-TDDFT, ne nécessitant que l'ajout de termes de correction spécifiques aux matrices de densité de transition et aux multiplicateurs lagrangiens.
• Étalonnage : Les auteurs étalonnent la méthode par rapport à des calculs multiréférence de haut niveau (XMS-CASPT2) pour le cation radical formaldéhyde ($CH_2O^+$).
• Application aux complexes de métaux de transition : La méthode est appliquée aux porphyrines de cuivre(II) (CuP, CuOEP, CuTPP) pour simuler les voies de relaxation d'état excité et les taux de conversion interne.

Résultats

• Précision de l'étalonnage ($CH_2O^+$) : Par rapport aux données de référence XMS-CASPT2, les NACMEs X-TDDFT réduisent l'écart absolu moyen (MAD) des normes des NACMEs U-TDDFT d'environ 32 % à 65 % (selon le fonctionnel). Les améliorations les plus significatives sont observées pour les états ayant un caractère CV(1) important, où l'U-TDDFT échoue souvent en raison de la contamination de spin. Cependant, les auteurs notent que pour certains états mixtes (par exemple, $1^2A_2$), la X-TDDFT peut occasionnellement produire des erreurs plus grandes que l'U-TDDFT, bien que cela soit attribué à une annulation fortuite d'erreurs dans les résultats U-TDDFT plutôt qu'à un échec de la théorie adaptée aux spins.
• Porphyrines de cuivre(II) :
  • Magnitudes des NACMEs : La X-TDDFT prédit des magnitudes de NACMEs significativement plus faibles (jusqu'à un ordre de grandeur) pour la transition $2T_1 \to 2S_0$ par rapport à l'U-TDDFT. Cette réduction est attribuée à l'annulation presque parfaite des contributions $\alpha$ et $\beta$ dans les états CV(1) purs au sein du cadre adapté aux spins, une annulation qui est incomplète dans l'U-TDDFT en raison des différences de forme des orbitales.
  • Taux de conversion interne : La réduction des normes de NACMEs entraîne des réductions drastiques des taux d'IC calculés (jusqu'à deux ordres de grandeur).
  • Voies de relaxation : Pour la CuP non substituée, l'U-TDDFT prédit que la voie de relaxation dominante implique une IC directe vers l'état fondamental et une IC indirecte via des états de champ de ligand approximativement égales. En revanche, la X-TDDFT prédit que la voie indirecte via les états de champ de ligand est dominante.
  • Effets des substituants : Les tendances relatives des taux d'IC pour les porphyrines substituées (CuP, CuOEP, CuTPP) diffèrent entre les deux méthodes. L'U-TDDFT prédit une tendance CuOEP < CuTPP < CuP, tandis que la X-TDDFT prédit CuTPP < CuOEP < CuP. Cela indique que les erreurs de l'U-TDDFT sont non seulement importantes mais aussi non systématiques, échouant à capturer les tendances relatives même pour des molécules étroitement apparentées.

Signification et revendications
L'article revendique que le développement des NACMEs X-TDDFT complète la boîte à outils essentielle pour l'étude des processus photophysiques et photochimiques des molécules à couche ouverte en utilisant la TDDFT adaptée aux spins. Les auteurs soulignent que l'adaptation aux spins n'est pas simplement une correction pour les énergies et les gradients, mais est cruciale pour les NACMEs et, par conséquent, pour la précision de la dynamique non adiabatique simulée et des taux d'IC.

L'étude conclut que le recours aux NACMEs U-TDDFT pour les systèmes à couche ouverte peut conduire à des prédictions qualitativement incorrectes des mécanismes de relaxation et à des erreurs quantitatives dans les constantes de vitesse qui sont à la fois importantes et non systématiques. Par conséquent, les auteurs recommandent l'utilisation de NACMEs X-TDDFT pour les systèmes à couche ouverte en général, même dans les cas où les énergies d'excitation U-TDDFT et X-TDDFT semblent similaires. Ce travail jette les bases de futures applications en dynamique moléculaire non adiabatique et dans l'étude de la photophysique ultrarapide des radicaux et des complexes de métaux de transition.