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

Deze studie introduceert KinLuv, een nieuw ab initio kinetisch model dat hogere aangeslagen toestanden en vibronische koppeling expliciet meeneemt om de fotofysische eigenschappen van organische emitters nauwkeurig te kwantificeren en criteria te leveren voor het selecteren van vereenvoudigde modellen.

De kern

Het Verborgen Verkeer in de Moleculaire Stad

Stel je voor dat een organisch lichtgevend molecuul (zoals die in je OLED-scherm of een gloeilampje) een drukke stad is. In deze stad wonen verschillende soorten "inwoners":

• S0: De rustige, slapende inwoners (de grondtoestand).
• S1: De energieke, springerige inwoners (de eerste opgewekte toestand).
• T1: De wat trage, dromerige inwoners (de eerste triplettoestand).

Vroeger dachten wetenschappers dat het leven in deze stad heel simpel was. Ze dachten: "Als we een lichtje op de stad schijnen, springen de inwoners naar S1. Soms vallen ze direct terug naar S0 (licht uit). Soms veranderen ze van karakter naar T1 en komen ze later weer terug naar S1 om licht te geven."

Dit was het drie-staten model. Het was als een simpele kaart van de stad met alleen de hoofdstraten.

Het Probleem: De Verkeersdrukte op de Hoge Wegen

De auteurs van dit artikel (Yue He en Daniel Escudero) zeggen: "Wacht even, die kaart is niet compleet!"

In werkelijkheid zijn er ook hogere verdiepingen in de stad:

• S2: Een nog energiekere inwoner.
• T2: Een nog trager, maar belangrijkere inwoner.

In sommige steden (moleculen) zijn deze hogere verdiepingen zo belangrijk dat ze het verkeer volledig veranderen. Als je ze negeert, begrijp je niet waarom het licht soms uitvalt of waarom het trager oplicht dan verwacht.

De analogie:
Stel je voor dat je een auto (een energiepakket) in de stad rijdt.

• In een simpele stad (zoals het molecuul DOBNA) kun je prima navigeren met alleen de hoofdstraten. De extra wegen (S2, T2) zijn er wel, maar ze zijn zo ver weg of zo slecht begaanbaar dat je ze nooit echt nodig hebt om je bestemming te bereiken.
• In een complexe stad (zoals het molecuul DiKTa of DBT) zijn de hogere wegen echter de snelste routes. Als je die negeert, beland je in een file. Je auto (de energie) raakt vast in de verkeerde straten en komt nooit meer aan bij het licht (de uitstraling).

De Oplossing: De Nieuwe Navigatie (KinLuv)

De onderzoekers hebben een nieuwe software ontwikkeld genaamd KinLuv.

• Wat doet het? Het is als een superkrachtige GPS die niet alleen de hoofdstraten kent, maar ook elke steeg, elke brug en elke tunnel (alle mogelijke routes tussen S0, S1, S2, T1, T2).
• Hoe werkt het? Het berekent precies hoe snel een auto van de ene weg naar de andere kan springen. Ze gebruiken hiervoor een speciale techniek (de "Herzberg-Teller" koppeling), die je kunt vergelijken met het begrijpen van hoe de trillingen van de auto (de vibraties van de atomen) helpen om over drempels te springen die normaal gesproken te hoog zijn. Zonder deze trillingen zou de auto vastlopen; met ze kan hij vliegen.

Wat hebben ze ontdekt?

Ze hebben dit systeem getest op verschillende "steden" (moleculen) en kwamen tot drie belangrijke conclusies:

1. Soms is de simpele kaart genoeg:
   Voor sommige moleculen (zoals DOBNA) werkt het oude, simpele model prima. De hogere wegen spelen geen grote rol. Je kunt de verkeersdrukte voorspellen zonder de hele stad te hoeven kennen.

2. Soms is de simpele kaart gevaarlijk:
   Voor andere moleculen (zoals DiKTa en DBT) is het cruciaal om de hogere wegen mee te nemen. Als je dat niet doet, denk je dat je 80% van je auto's (energie) veilig aankomt, terwijl in werkelijkheid slechts 26% aankomt. De hogere wegen fungeren hier als "afvalbakken" waar energie verloren gaat, of juist als "snelwegen" die het licht sneller laten ontstaan.

3. De balans vinden:
   Het grootste inzicht is niet dat we altijd alles moeten weten. Het gaat om het vinden van de juiste balans.

   • Als je te simpel kijkt, mis je de waarheid.
   • Als je te complex kijkt (en elke mogelijke weg meetelt), kun je de fouten in je berekeningen zo groot maken dat de uitkomst weer onnauwkeurig wordt.
   • KinLuv helpt wetenschappers te beslissen: "Voor dit specifieke molecuul heb ik de T2-straat nodig, maar de S2-straat kan ik negeren."

Waarom is dit belangrijk?

Dit onderzoek helpt bij het ontwerpen van beter licht.
Stel je voor dat je een nieuwe OLED-televisie of een efficiënter zonnepaneel wilt maken. Je wilt dat het licht zo helder en zuinig mogelijk is.
Met deze nieuwe "GPS" (KinLuv) kunnen ingenieurs in de computer simuleren hoe een nieuw molecuul zich zal gedragen voordat ze het in het lab bouwen. Ze kunnen zien: "Oh, als we deze atoomgroep veranderen, opent zich een nieuwe weg naar T2 die onze energie verslindt. Laten we dat niet doen."

Kort samengevat:
De onderzoekers hebben een slimme rekenmethode bedacht om te kijken of we in de complexe wereld van lichtgevende moleculen moeten kijken naar de "verdiepingen boven de begane grond". Soms zijn ze irrelevant, maar vaak zijn ze de sleutel tot waarom het licht wel of niet werkt. Hierdoor kunnen we in de toekomst betere, helderdere en zuiniger schermen en lampen ontwerpen.

Technische samenvatting

Probleemstelling

Traditionele modellen voor de excitatietoestandsdynamica van organische chromoforen (zoals TADF-materialen) baseren zich vaak op een vereenvoudigd driedelig model: de grondtoestand ($S_0$), de laagste singlet-toestand ($S_1$) en de laagste triplet-toestand ($T_1$). Hoewel deze benadering vaak werkt, falen ze in systemen waar hoger gelegen excitatietoestanden ($S_n$ met $n>1$ en $T_n$ met $n>1$) een cruciale rol spelen.

• De lacune: Hoewel er kwalitatief bekend is dat deze hogere toestanden betrokken zijn bij processen zoals anti-Kasha-emissie, inter-systeem kruising (ISC) via $T_n$, en interne conversie (IC), ontbreekt er een systematische, kwantitatieve beoordeling van hun impact op meetbare grootheden zoals de fotoluminescentie-kwantumopbrengst (PLQY) en levensduur.
• Uitdaging: Experimenteel is het moeilijk om de bijdrage van deze toestanden uniek te bepalen zonder complexe tijdsopgeloste spectroscopie. Bestaande theoretische modellen gebruiken vaak vereenvoudigingen (zoals het negeren van Herzberg-Teller vibronische koppeling), wat leidt tot grote afwijkingen in de berekende snelheidsconstanten (bijvoorbeeld een onderschatting van de $rISC$-snelheid met meerdere ordes van grootte).

Methodologie

De auteurs hebben een nieuw, volledig ab initio raamwerk ontwikkeld om deze kloof te dichten:

1. KinLuv Software: Ze hebben een Python-gebaseerde kinetische modelleringstool genaamd KinLuv ontwikkeld. Deze tool lost gekoppelde differentiaalvergelijkingen (ODE's) op voor multistate-systemen om populatie-evolutie, PLQY en levensduur te voorspellen.
2. Multistate Modellen: In plaats van alleen $S_0, S_1, T_1$, worden modellen geanalyseerd die variëren in complexiteit:
   • 2-toestand: $S_0, S_1$
   • 3-toestand: $S_0, S_1, T_1$
   • 4-toestand: $S_0, S_1, T_1, T_2$
   • 5-toestand: $S_0, S_1, S_2, T_1, T_2$
3. Berekening van Snelheidsconstanten: Alle snelheidsconstanten voor stralende en niet-stralende overgangen (fluorescentie, fosforescentie, ISC, rISC, IC, rIC) worden berekend met de Fermi's Gouden Regel.
   • Cruciaal: De berekeningen omvatten expliciet de Herzberg-Teller (HT) vibronische koppeling. Dit betekent dat de afhankelijkheid van de overgangsmomenten en spin-baan-koppelingsmatrixelementen (SOCME) van de nucleaire coördinaten wordt meegenomen, in plaats van alleen de Franck-Condon (FC) benadering.
   • Niveau van theorie: Geometrieën worden geoptimaliseerd met TD(A)-CAM-B3LYP. Adiabatische excitatie-energieën worden verfijnd met SCS-ADC(2) voor hogere nauwkeurigheid.
4. Testsystemen: Het raamwerk werd toegepast op een hiërarchische reeks moleculen:
   • Prototypische MR-TADF-emitters: DOBNA en DiKTa.
   • Derivaten: DABNA-1 en DQAO.
   • Niet-TADF organische chromoforen: Dibenzothiophene (DBT) en PPDs-1.

Belangrijkste Bijdragen en Resultaten

1. Dominantie van Vibronische Koppeling
De analyse toont aan dat snelle (reverse) inter-systeem kruising ((r)ISC) in deze systemen grotendeels wordt gedreven door Herzberg-Teller vibronische koppeling. De FC-benadering alleen zou de snelheidsconstanten ernstig onderschatten. Bijvoorbeeld, voor DABNA-1 leverde de HT-benadering een $rISC$-snelheid op van $\sim 2.4 \times 10^4 s^{-1}$, wat veel dichter bij experimentele waarden ligt dan eerdere theoretische schattingen die HT negeerden.

2. Kwantitatieve Validatie van Kinematische Modellen
De studie vergelijkt de voorspellingen van de verschillende modellen (2 tot 5 toestanden) met experimentele data:

• DOBNA en DABNA-1 (Snel (r)ISC is traag): Voor deze systemen zijn de (r)ISC-processen traag in vergelijking met fluorescentie en IC. De populatie van hogere toestanden blijft klein.
  • Resultaat: Het eenvoudige 3-toestand model is voldoende om PLQY en levensduur nauwkeurig te voorspellen. Het toevoegen van $T_2$ of $S_2$ verbetert het mechanistische inzicht, maar verandert de kwantitatieve uitkomsten nauwelijks.
• DiKTa en DQAO (Snel (r)ISC): Hier zijn de (r)ISC-snelheden hoog, wat leidt tot aanzienlijke populaties van triplettoestanden.
  • Resultaat: Het 3-toestand model faalt in het voorspellen van de PLQY (voorspelt te hoog). Het 4-toestand model (inclusief $T_2$) is essentieel om de experimentele PLQY en levensduur nauwkeurig te reproduceren. $T_2$ fungeert hier als een concurrerend niet-stralend kanaal of een tussenstap.
• DBT en PPDs-1 (Niet-TADF):
  • Voor DBT is de invoering van $T_2$ (4-toestand model) noodzakelijk omdat $S_1 \to T_2$ een snelle ISC-pijp is die de PLQY verlaagt.
  • Voor PPDs-1 is zelfs het 5-toestand model (inclusief $S_2$) nodig. Hier fungeert $S_2$ als een zeer efficiënt niet-stralend kanaal ($S_1 \to S_2 \to S_0$) dat de PLQY drastisch verlaagt. Zonder $S_2$ zou het model een hoge PLQY voorspellen, wat in strijd is met het experiment.

3. Criteria voor Modelselectie
De auteurs stellen duidelijke criteria op voor wanneer een vereenvoudigd model volstaat en wanneer complexere modellen nodig zijn:

• Een 3-toestand model is voldoende als (r)ISC traag is ten opzichte van fluorescentie/IC.
• Een 4-toestand model (met $T_2$) is vereist als snelle ISC via hogere triplets leidt tot significante populaties die de PLQY beïnvloeden.
• Een 5-toestand model (met $S_2$) is vereist als snelle reverse interne conversie (rIC) of IC via hogere singlettoestanden een dominante niet-stralende route vormt.

Significantie en Conclusie

Dit onderzoek biedt een principieel raamwerk voor het rationaliseren van het ontwerp van organische emitters:

• Balans tussen Inzicht en Robuustheid: Het helpt onderzoekers de juiste balans te vinden tussen mechanistische volledigheid (alle fysiek toegankelijke toestanden opnemen) en numerieke robuustheid (vermijden van foutopbouw door te veel onzekere parameters).
• In silico Ontwerp: Door te identificeren wanneer hogere toestanden kwantitatief relevant zijn, kunnen materialenwetenschappers gerichter complexe kinetische modellen selecteren voor het voorspellen van de prestaties van nieuwe TADF- en niet-TADF-materialen.
• Methodologische Vooruitgang: De integratie van HT-vibronische koppeling in een volledig ab initio kinetisch model lost het probleem op van het onderschatten van (r)ISC-snelheden en biedt een nauwkeurigere beschrijving van de excitatietoestandsdynamica dan eerdere benaderingen.

Samenvattend bewijst deze studie dat hoewel een driedelig model soms toereikend is, het systematisch negeren van hoger gelegen excitatietoestanden kan leiden tot fundamenteel onjuiste voorspellingen van de efficiëntie en kinetiek van organische lichtemitters. Het ontwikkelde KinLuv-raamwerk is een krachtig hulpmiddel om deze complexiteit te kwantificeren.