Quantifying the Role of Higher-Lying Excited States in Organic Emitters via Multistate Ab Initio Kinetic Modeling

Cet article présente KinLuv, un modèle cinétique ab initio multistat qui quantifie l'influence des états excités supérieurs sur les rendements quantiques et les durées de vie des émetteurs organiques, permettant ainsi de sélectionner des modèles simplifiés pour la conception rationnelle de matériaux photophysiques performants.

L'essentiel

🌟 Le Grand Voyage des Électrons : Pourquoi les modèles simples ne suffisent plus

Imaginez que vous essayez de comprendre comment fonctionne une ville très animée (une molécule organique) quand elle est excitée par la lumière. Traditionnellement, les scientifiques utilisaient une carte très simple, comme un plan de métro à trois stations seulement :

1. La Gare du Départ (État Fondamental S0) : La molécule au repos.
2. La Station Express (État Singulet S1) : La molécule excitée qui brille immédiatement (fluorescence).
3. La Station de Repos (État Triplet T1) : Une version "endormie" de la molécule qui peut se réveiller plus tard pour briller à nouveau (fluorescence retardée).

Pour beaucoup de molécules, ce plan à trois stations fonctionnait bien. Mais les chercheurs ont découvert que pour certaines molécules modernes (utilisées dans les écrans OLED de nouvelle génération), ce plan était incomplet. Il manquait des gares cachées (les états excités supérieurs, S2 et T2) qui jouent un rôle crucial, mais qu'on ignorait.

🚀 L'Innovation : Le Nouveau Logiciel "KinLuv"

L'équipe de chercheurs, dirigée par Daniel Escudero, a créé un nouveau logiciel appelé KinLuv. Au lieu de regarder seulement les trois stations principales, KinLuv dessine une carte complète incluant toutes les gares possibles, même celles situées plus haut dans le ciel (les états S2 et T2).

L'analogie du trafic routier :
Imaginez que les électrons sont des voitures.

• L'ancien modèle (3 états) : Disait que les voitures allaient directement de la gare A à la gare B, ou parfois faisaient un détour par la gare C.
• Le nouveau modèle (KinLuv) : Réalise que parfois, les voitures doivent passer par une autoroute très rapide (l'état T2) ou un pont suspendu (l'état S2) pour arriver à destination. Si on ignore ces routes, on ne comprend pas pourquoi il y a des embouteillages ou pourquoi certaines voitures arrivent plus vite que prévu.

🔍 La Révélation : Quand faut-il regarder les gares cachées ?

En testant ce logiciel sur plusieurs types de molécules, les chercheurs ont découvert une règle d'or pour savoir quand il faut utiliser la carte complète ou se contenter de la carte simple :

1. Le Cas "Tranquille" (Exemple : DOBNA) :

   • La situation : Les voitures (électrons) vont lentement. Elles passent beaucoup de temps à la gare principale (S1) et à la gare de repos (T1).
   • Le résultat : Les gares cachées (S2, T2) sont là, mais personne ne les utilise vraiment.
   • Conclusion : Pour ces molécules, l'ancien modèle à 3 stations suffit. On peut faire des économies de temps de calcul sans perdre en précision.

2. Le Cas "Autoroute" (Exemple : DiKTa et DBT) :

   • La situation : Les voitures vont à toute vitesse ! Elles passent par les gares cachées (T2) comme des éclairs pour changer de voie.
   • Le résultat : Si on utilise l'ancien modèle, on se trompe complètement sur le nombre de voitures qui arrivent à destination (le rendement lumineux) et sur le temps qu'elles mettent.
   • Conclusion : Ici, il est obligatoire d'inclure les gares cachées (au moins T2) pour avoir un résultat juste. Sans cela, on prédit que la molécule brille beaucoup, alors qu'en réalité, elle perd son énergie sur ces routes rapides.

3. Le Cas "Pont Suspendu" (Exemple : PPDs-1) :

   • La situation : Il existe un pont très haut (l'état S2) qui permet aux voitures de sauter directement hors de la ville sans jamais s'arrêter.
   • Le résultat : Si on ignore ce pont, on pense que la ville est pleine de voitures, alors qu'elles ont toutes sauté par-dessus bord !
   • Conclusion : Pour ces molécules, il faut inclure l'état S2 pour comprendre pourquoi elles ne brillent presque pas.

🎨 Le Secret de la Vitesse : La "Danse" des Atomes

Un autre point clé de l'article est pourquoi ces transitions sont si rapides.
Les chercheurs ont montré que ce n'est pas seulement une question d'énergie, mais de vibration.

• Imaginez que pour passer d'une gare à l'autre, la molécule doit faire une petite danse.
• Les chercheurs ont découvert que les mouvements des atomes (les vibrations) aident énormément les électrons à changer de voie, un peu comme si le sol vibrait pour aider les voitures à faire un saut périlleux.
• Le logiciel KinLuv prend en compte cette "danse" (appelée couplage vibronique), ce qui rend les prédictions beaucoup plus réalistes que les anciens modèles qui ignoraient ce mouvement.

🏁 En Résumé : Pourquoi c'est important ?

Cette recherche est comme un manuel de mise à jour pour les ingénieurs qui créent les écrans de demain (téléphones, TV OLED).

• Avant : On utilisait une règle simple pour tout le monde. Ça marchait pour certains, mais ça échouait pour les meilleurs.
• Maintenant : Grâce à KinLuv, les scientifiques peuvent dire : "Pour cette molécule précise, on a besoin de la carte complète avec les gares cachées, sinon on va gaspiller de l'argent en fabriquant un écran qui ne brille pas assez."

C'est une avancée majeure pour concevoir des matériaux plus brillants, plus efficaces et plus durables, en sachant exactement quelles "routes" les électrons empruntent dans leur voyage lumineux.

Résumé technique

1. Problématique et Contexte

La dynamique des états excités des chromophores organiques est traditionnellement modélisée par une approximation à trois états : l'état fondamental ($S_0$), le premier état singulet excité ($S_1$) et le premier état triplet excité ($T_1$). Bien que cette approche soit souvent suffisante pour des systèmes où les états supérieurs ($S_{n>1}, T_{n>1}$) sont énergétiquement bien séparés, elle devient inadéquate dans de nombreux cas modernes, notamment pour les émetteurs à fluorescence retardée activée thermiquement (TADF) et les réactions photochimiques anti-Kasha.

Les limites de ce modèle simplifié sont multiples :

• Influence des états supérieurs : Des études récentes suggèrent que les états excités supérieurs jouent un rôle crucial dans les processus de conversion intersystème (ISC) et de conversion intersystème inverse (rISC), agissant parfois comme des canaux intermédiaires de mélange de spin.
• Déficience des modèles existants : La plupart des analyses cinétiques quantitatives reposent sur des ajustements de données expérimentales (décaissements bi-exponentiels) ou sur des calculs théoriques qui négligent le couplage vibronique de type Herzberg-Teller (HT). Par exemple, les travaux antérieurs de Shizu et al. sous-estiment considérablement les taux de rISC en omettant l'effet HT, conduisant à des écarts de plusieurs ordres de grandeur par rapport à l'expérience.
• Manque de cadre ab initio complet : Il existe un vide dans la compréhension quantitative de l'impact explicite des états supérieurs sur les observables macroscopiques (rendement quantique de photoluminescence - PLQY, durées de vie) au sein d'un cadre de modélisation cinétique entièrement ab initio.

2. Méthodologie

Pour combler ce vide, les auteurs ont développé une approche intégrée combinant la chimie quantique de haut niveau et la modélisation cinétique multistat.

A. Développement du logiciel KinLuv

Les auteurs ont créé KinLuv, un package Python conçu pour résoudre les systèmes d'équations différentielles ordinaires (EDO) décrivant la cinétique des états excités.

• Modèles hiérarchiques : Le logiciel permet de comparer quatre modèles cinétiques d'complexité croissante :
  1. Deux états : $S_0, S_1$ (fluorescence prompte uniquement).
  2. Trois états : $S_0, S_1, T_1$ (modèle TADF standard).
  3. Quatre états : $S_0, S_1, T_1, T_2$ (incluant un triplet supérieur).
  4. Cinq états : $S_0, S_1, S_2, T_1, T_2$ (incluant un singulet et un triplet supérieurs).
• Calcul des observables : KinLuv calcule les durées de vie de la fluorescence prompte et retardée, ainsi que le PLQY, en résolvant les EDOs pour la période de décroissance ou en appliquant l'approximation de l'état stationnaire (SSA).

B. Calcul des constantes de vitesse Ab Initio

Les constantes de vitesse pour tous les processus (fluorescence, phosphorescence, conversion interne IC, ISC, et leurs réverses) sont calculées en utilisant la Règle d'Or de Fermi (FGR).

• Couplage Vibronique : Contrairement aux approches précédentes, cette étude inclut explicitement l'effet Herzberg-Teller (HT). Cela permet de prendre en compte la dépendance des éléments de matrice de couplage (dipôle de transition, couplage spin-orbite, couplage non adiabatique) par rapport aux coordonnées nucléaires, brisant ainsi l'approximation de Condon.
• Niveaux de théorie :
  • Géométries et fréquences : TD(A)-CAM-B3LYP (avec correction de dispersion D3).
  • Énergies d'excitation adiabatiques : Calculs de points uniques SCS-ADC(2) pour une précision accrue des gaps énergétiques ($\Delta E_{ST}$).
  • Couplages : Calculs des éléments de matrice de couplage spin-orbite (SOCME) et non adiabatiques (NACME) avec des corrections appropriées.

3. Résultats Clés

L'approche a été validée sur six chromophores organiques : deux émetteurs MR-TADF prototypes (DOBNA, DiKTa), leurs dérivés (DABNA-1, DQAO), et deux chromophores non-TADF (DBT, PPDs-1).

A. Rôle du couplage vibronique (HT)

• Les calculs montrent que le couplage HT est dominant pour les processus ISC et rISC. Pour la plupart des transitions, la contribution HT dépasse 90-100% de la constante de vitesse totale.
• L'inclusion de l'effet HT corrige les sous-estimations massives observées dans les travaux antérieurs (ex: pour DABNA-1, le taux rISC passe de ~2,5 s⁻¹ sans HT à ~2,4 $\times$ 10⁴ s⁻¹ avec HT, se rapprochant des valeurs expérimentales).

B. Dynamique des émetteurs MR-TADF (DOBNA vs DiKTa)

• DOBNA (Ralentissement des processus) : Dans ce système, les taux ISC/rISC sont lents par rapport à la fluorescence et à la conversion interne. L'inclusion des états $T_2$ et $S_2$ n'altère pas significativement le PLQY ou les durées de vie prédites par le modèle à trois états. Le modèle simplifié reste suffisant pour reproduire les observables expérimentaux.
• DiKTa (Processus rapides) : Ce système présente des taux ISC/rISC très rapides. Ici, le modèle à trois états échoue à prédire le PLQY expérimental (surestimation). L'inclusion explicite de l'état $T_2$ est cruciale car il introduit une voie de décaissement non radiative compétitive ($S_1 \to T_2 \to T_1 \to S_0$) qui réduit le PLQY. Le modèle à quatre ou cinq états reproduit alors avec précision le PLQY et la durée de vie retardée.

C. Émetteurs Non-TADF (DBT et PPDs-1)

• DBT : La formation de triplets via $S_1 \to T_2$ est rapide. L'inclusion de $T_2$ est nécessaire pour reproduire le faible PLQY expérimental (~9%).
• PPDs-1 : Ce système est dominé par une conversion interne rapide via l'état singulet supérieur $S_2$ ($S_1 \to S_2 \to S_0$). Seul le modèle à cinq états (incluant $S_2$) permet de reproduire le PLQY extrêmement faible (0,9%) et la durée de vie ultra-courte observés expérimentalement.

4. Contributions Majeures

1. Développement de KinLuv : Un outil logiciel open-source permettant une modélisation cinétique multistat flexible et automatisée, intégrant des constantes de vitesse calculées ab initio.
2. Validation de l'effet Herzberg-Teller : Démonstration systématique que l'omission de l'effet HT conduit à des erreurs quantitatives majeures dans la prédiction des taux ISC/rISC, rendant les modèles sans HT inapplicables pour une prédiction précise.
3. Critères de sélection de modèles : Établissement de critères clairs pour déterminer la complexité minimale requise d'un modèle cinétique :
   • Si les processus (r)ISC sont lents $\rightarrow$ Modèle à 3 états suffisant.
   • Si les processus (r)ISC sont rapides et peuplent les triplets supérieurs $\rightarrow$ Nécessité d'inclure $T_2$ (modèle à 4 états).
   • Si la conversion interne via $S_2$ est dominante $\rightarrow$ Nécessité d'inclure $S_2$ (modèle à 5 états).
4. Précision Quantitative : Pour la première fois, une modélisation cinétique ab initio complète reproduit simultanément le PLQY et les durées de vie (prompte et retardée) pour une variété de systèmes organiques complexes, y compris ceux où les états supérieurs sont déterminants.

5. Signification et Impact

Ce travail marque une avancée significative dans la conception rationnelle de matériaux émetteurs organiques pour les OLEDs et autres applications optoélectroniques.

• Au-delà de l'intuition qualitative : Il transforme le rôle des états excités supérieurs d'une simple curiosité mécanistique en un paramètre quantifiable essentiel pour la prédiction des performances.
• Optimisation du compromis : Il offre une stratégie pour choisir le modèle cinétique le plus simple possible tout en restant physiquement robuste, évitant ainsi la sur-paramétrisation inutile tout en garantissant la précision numérique.
• Design in silico : En fournissant un cadre fiable pour prédire le PLQY et les durées de vie, cette méthodologie permet d'accélérer le criblage virtuel de nouveaux émetteurs TADF et non-TADF, en identifiant dès la phase de conception si des états supérieurs risquent de dégrader les performances.

En conclusion, l'article démontre que l'inclusion explicite des états excités supérieurs et du couplage vibronique n'est pas seulement un raffinement théorique, mais une nécessité quantitative pour comprendre et prédire le comportement photophysique d'une large classe de matériaux organiques émetteurs.