Angeregte Zustände in Materialien, die Licht emittieren
1. Einleitung
Bei LEDs und OLEDs erfolgt die Lichtemission nicht direkt durch freie Ladungsträger, sondern durch angeregte Zustände (Exzitonen), die bei der Rekombination von Elektronen und Löchern entstehen. Diese angeregten Zustände existieren als Singuletts und Tripletts, die jeweils unterschiedliche Abklingmechanismen, Wechselwirkungsmechanismen und Strahlungseffizienzen aufweisen. Das Gleichgewicht zwischen ihrer Entstehung, Umwandlung und ihrem Verlust bestimmt letztendlich die Helligkeit, Effizienz und das Abfallverhalten des Bauelements. Folglich muss ein realistisches Modell die Populationen dieser angeregten Zustände explizit verfolgen. Da Exzitonen durch die Rekombination von Ladungsträgern erzeugt werden, erfordert dies eine Kopplung an eine Drift-Diffusions-Beschreibung des Ladungstransports, die die lokalen Ladungsträgerdichten und Rekombinationsraten liefert, die als Anregungsterme für das System der angeregten Zustände dienen.
In der Praxis bedeutet dies, dass die Drift-Diffusions-Gleichungen zusammen mit einer Reihe von gekoppelten Geschwindigkeitsgleichungen gelöst werden, die die Konzentrationen von Singulett- und Triplett-Zuständen beschreiben, und gegebenenfalls auch die Dynamik von Photonen oder optischen Moden. Die genaue Kombination der Gleichungen hängt von dem untersuchten Problem ab: So kann beispielsweise die Modellierung der Effizienz von OLEDs möglicherweise nur die Verwendung von Exzitonen-Geschwindigkeitsgleichungen erfordern, während die Modellierung organischer Laser zusätzlich die Berücksichtigung der Kopplung an Photonen-Konzentrationen und Kavitätseffekte erfordert. Der Detaillierungsgrad kann daher je nach Bedarf angepasst werden, d. h. je nachdem, ob eine minimale oder eine umfassendere physikalische Beschreibung des Bauelements erforderlich ist.
💡 Hinweis: Für die meisten OLED-Simulationen ist es nicht erforderlich, angeregte Zustände zu berücksichtigen. Dies ist eine fortgeschrittene Funktion, die hauptsächlich dann von Bedeutung ist, wenn die Effizienz bei hohen Stromdichten modelliert werden soll, da in diesem Fall die Wechselwirkungen zwischen Exzitonen und zwischen Exzitonen und Polaronen eine wichtige Rolle spielen.
2. Überblick
OghmaNano kann bis zu vier gekoppelte Singlet- und Triplett-Raten-Gleichungen zusätzlich zu den Drift-Diffusions-Gleichungen in 1D lösen, wie in ?? schematisch dargestellt. Dies sind die Singlet-Gleichung im Wirtsmaterial, die Triplett-Gleichung im Wirtsmaterial, die Singlet-Gleichung im Dotierungsmaterial und die Triplett-Gleichung im Dotierungsmaterial. Das Modell kann aktiviert werden, indem zunächst die Schaltfläche „OLED-Kinetik“ deaktiviert und anschließend über das Dropdown-Menü ausgewählt wird, welche Gleichungen benötigt werden, wie in ?? gezeigt. Für die meisten OLED-Simulationen sind normalerweise nur die ersten beiden Gleichungen erforderlich, die den Singlet- und Triplett-Konzentrationen im Wirtsmaterial entsprechen. Für Benutzer, die den Laserbetrieb simulieren möchten, steht außerdem eine Photonenraten-Gleichung zur Verfügung. Diese Gleichungen werden im Folgenden ausführlich erläutert.
Dieses Modell wurde in den folgenden Veröffentlichungen vorgestellt, in denen Sie mehr über seine Anwendung auf experimentelle Daten erfahren können: A. Bickerdike et al. Aufklärung der räumlich-zeitlichen Exzitonen-Dynamik in elektrisch gepumpten organischen Laserdioden und A. Bickerdike et al. Inverser Intersystem-Übergang und Exzitonen-Ernte: Ein Weg zu einer niedrigen Schwellenspannung in mehrfach resonierenden TADF-organischen Halbleiterlasern, die elektrisch gepumpt werden.
3. Drift-Diffusions-Gleichungen
Bevor die Ratengleichungen für den angeregten Zustand gelöst werden, ist es zunächst erforderlich, zu bestimmen, wie Ladungsträger in dem Bauelement injiziert, transportiert und rekombiniert werden. Dies geschieht mithilfe der Drift-Diffusions-Gleichungen, die die Elektron- und Lochdichten sowie die Rekombinationsraten liefern, welche als Anregungsterme für das System im angeregten Zustand dienen. Dies sind die Standard-Drift-Diffusions-Gleichungen, die zur Beschreibung des Ladungstransports in Halbleitern verwendet werden, einschließlich der Drift von Ladungsträgern im elektrischen Feld, der Diffusion aufgrund von Konzentrationsgradienten und der bimolekularen Rekombination. Eine vollständige Beschreibung dieser Gleichungen finden Sie hier: Drift-Diffusions-Theorie, und eine detaillierte Ableitung hier: Ableitung der Drift-Diffusions-Gleichungen.
In diesen Gleichungen sind \(n\) und \(p\) die Elektron- bzw. Lochdichten, \(\mathbf{J}_n\) und \(\mathbf{J}_p\) die entsprechenden Stromdichten, \(\mu_n\) und \(\mu_p\) die Beweglichkeiten, \(D_n\) und \(D_p\) die Diffusionskoeffizienten, \(\phi\) das elektrostatische Potenzial, \(G_n\) und \(G_p\) die Erzeugungsraten, \(R_{\mathrm{free}}\) die bimolekulare Rekombinationsrate mit dem Koeffizienten \(k_r\) und \(R_{\mathrm{trap}}\) die durch Fallen unterstützte Rekombinationsrate. Die Größen \(n_0\) und \(p_0\) bezeichnen die Gleichgewichtsträgerdichten.
Elektronische Kontinuität
\[ \frac{\partial n}{\partial t} = \frac{1}{q}\nabla \cdot \mathbf{J}_n \;-\; R_{\mathrm{free}} \;-\; R_{\mathrm{trap}} \;+\; G_n, \qquad \mathbf{J}_n = -\,q\mu_n n \nabla\phi \;+\; q D_n \nabla n \]
Lochkontinuität
\[ \frac{\partial p}{\partial t} = -\,\frac{1}{q}\nabla \cdot \mathbf{J}_p \;-\; R_{\mathrm{free}} \;-\; R_{\mathrm{trap}} \;+\; G_p, \qquad \mathbf{J}_p = \;q\mu_p p \nabla\phi \;-\; q D_p \nabla p \]
Anzahl der Einzelteilchen (Wirt)
\[ \frac{dN_S}{dt} = \frac{1}{4}\big(\gamma_{\mathrm{free}} R_{\mathrm{free}} + \gamma_{\mathrm{trap}} R_{\mathrm{trap}}\big) + \frac{1}{4}\kappa_{TT}N_T^{2} - (\kappa_{\mathrm{FRET}}P_{OD} + \kappa_S + \kappa_{ISC})N_S - \Big(\tfrac{7}{4}\kappa_{SS}N_S + \kappa_{SP}(n_{\mathrm{free}}+p_{\mathrm{free}}) + \kappa_{ST}N_T\Big)N_S \]
Populationsgröße der Dreifach-Wirtsart
\[ \frac{dN_T}{dt} = \frac{3}{4}\big(\gamma_{\mathrm{free}} R_{\mathrm{free}} + \gamma_{\mathrm{trap}} R_{\mathrm{trap}}\big) + \kappa_{ISC}N_S + \frac{3}{4}\kappa_{SS}N_S^{2} - (\kappa_{DEXT}P_{OD} + \kappa_T + \kappa_{TP}(n_{\mathrm{free}}+p_{\mathrm{free}}))N_T - \frac{5}{4}\kappa_{TT}N_T^{2} \]
Anzahl der einzelnen Dotieratome
\[ \frac{dN_{SD}}{dt} = \kappa_{\mathrm{FRET}}P_{OD}N_S + \frac{1}{4}\kappa_{TTD}N_{TD}^{2} - (\kappa_{SD} + \kappa_{ISCD})N_{SD} - \Big(\tfrac{7}{4}\kappa_{SSD}N_{SD} + \kappa_{SPD}(n_{\mathrm{free}}+p_{\mathrm{free}}) + \kappa_{STD}N_{TD}\Big)N_{SD} - \xi P_{HO}\big(N_{SD} - WN_{OD}\big) \]
Triplett-Zustand (Dotierungsmittel)
\[ \frac{dN_{TD}}{dt} = \kappa_{DEXT}P_{OD}N_T + \kappa_{ISCD}N_{SD} + \frac{3}{4}\kappa_{SSD}N_{SD}^{2} - \kappa_{TD}N_{TD} - \frac{5}{4}\kappa_{TTD}N_{TD}^{2} - \kappa_{TPD}N_{TD}(n_{\mathrm{free}}+p_{\mathrm{free}}) \]
Photonenanzahl (optischer Modus)
\[ \frac{dP_{HO}}{dt} = \beta_{sp}\kappa_{SD}N_{SD} + \big(\Gamma \xi (N_{SD} - WN_{OD}) - \kappa_{CAV}\big) P_{HO} \]
Beschränkung der Dotierstoffkonzentration
\[ N_{OD} = N_{DOP} - N_{SD} - N_{TD} \]
💡 Hinweis: Die Gleichungen für das Wirts-Singulett und -Triplett können auch Diffusionsglieder enthalten, die über das Dropdown-Menü oben im elektrischen Editor aktiviert werden können. In der Regel bleiben diese jedoch deaktiviert, da sie die numerische Lösung des Systems erschweren, die Berechnungszeit erhöhen und die Stabilität des Solvers verringern. In den meisten Fällen ist dies nicht erforderlich, da Exzitonen im Allgemeinen kurze Diffusionslängen aufweisen, sodass die Aktivierung der Diffusion nur dann einen begrenzten Vorteil bietet, wenn dies ausdrücklich erforderlich ist.
4. Lichtemission
Die Lichtemission im Modell entsteht durch die Rekombination von Ladungsträgern und/oder durch Exzitonenrelaxationsprozesse. Dieses emittierte Licht wird direkt in die berechneten Ausgaben für die Quanteneffizienz eingespeist, einschließlich der externen Quanteneffizienz (EQE), der Licht-Spannungs-Kennlinien (L–V), der Licht-Strom-Kennlinien (L–I) und des Emissionsspektrums des Materials. Der Punkt, an dem das Licht aus dem Modell extrahiert wird, hängt davon ab, welche physikalischen Prozesse aktiviert sind. Wenn das Exzitonenmodell deaktiviert ist, wird die Emission direkt aus den Rekombinationstermen (\(R_{\mathrm{free}}\) und \(R_{\mathrm{trap}}\)) entnommen. Wenn nach und nach mehr Gleichungen aktiviert werden, wird die Emission stattdessen aus den entsprechenden exzitonischen oder optischen Populationen entnommen. Die spezifischen Emissionswege, die bei verschiedenen Modellkonfigurationen verwendet werden, sind unten zusammengefasst.
Je nachdem, welche Gleichungen aktiviert sind, kann das Modell Licht auf vier verschiedene Arten aussenden.
- Wenn das Singulett- und Triplett-Modell deaktiviert wird, wird Licht direkt durch die Rekombination freier Ladungsträger emittiert; dies umfasst auch die Rekombination über Fallezustände, die als \(R_{\mathrm{trap}}\) dargestellt werden können; dieser Beitrag kann unabhängig aktiviert oder deaktiviert werden.
- Wenn die Gleichungen für den Singulett- und Triplett-Zustand aktiviert sind, wird die Emission anhand des Zerfalls des Singulett-Zustands berechnet.
- Wenn die Gleichungen für Wirt-Dotierstoff-Wechselwirkungen aktiviert sind, wird die Emission aus der Dotierstoff-Singulett-Population berechnet.
- Wenn die Gleichung für die Photonenrate aktiviert ist, wird die Emission anhand des Verlustterms für die Kavität berechnet.
| Modellkonfiguration | Emissionsquelle | Photonenerzeugungsrate |
|---|---|---|
| Kein Exzitonenmodell. | Freie Trägerrekombination | \(G_\gamma = \eta_{\mathrm{ph}}\left(R_{\mathrm{free}} + R_{\mathrm{trap}}\right)\) |
| Singulett-/Triplett-Zustände möglich | Zerfall eines Host-Singlets | \(G_\gamma = \eta_{\mathrm{ph}}\,\kappa_S N_S\) |
| Singulett-/Triplett-Zustand + Dotierung aktiviert | Zerfall eines Dotierungssinguletts | \(G_\gamma = \eta_{\mathrm{ph}}\,\kappa_{SD} N_{SD}\) |
| Photonenzählraten-Gleichung aktiviert. | Verlust durch Hohlräume | \(G_\gamma = \kappa_{\mathrm{CAV}} P_{HO}\) |
5. Ergebnisse/Ausgaben
Das Modell für angeregte Zustände erzeugt eine Reihe zusätzlicher Ausgaben, die zur Analyse des Verhaltens von OLED- und Laserstrukturen verwendet werden können. Räumliche Momentaufnahmen der Singulett- und Triplett-Populationen können im Momentaufnahmen-Verzeichnis eingesehen werden, wie in ?? dargestellt. In diesem Beispiel werden die Dateien EX_Ns.csv und EX_Nt.csv als Funktion der Position dargestellt, was den jeweiligen Singulett- und Triplett-Dichten entspricht. Es können auch zeit- oder spannungsabhängige Ausgaben generiert werden, wie in ?? veranschaulicht, wobei die Triplett-Dichte als Funktion der angelegten Spannung dargestellt wird. Die Wirkung des Modells für angeregte Zustände spiegelt sich auch in den berechneten Leistungskennwerten des Bauelements wider. Beispielsweise sinkt die in ?? gezeigte EQE-Kurve bei höherer Spannung aufgrund der Singulett-Triplett-Annihilation. Ebenso berücksichtigt die in ?? gezeigte Licht-Spannungs-Kurve die Wirkung der Singulett-Triplett-Gleichungen, indem sie darstellt, wie das Licht aus dem Modell herausgekoppelt wird.
EX_Ns.csv und EX_Nt.csv, die den Singulett- bzw. Triplett-Zustand darstellen und als Funktion der Position aufgetragen sind.
6. Variablen und Konstanten
Zustandsvariablen
| Symbol | Beschreibung |
|---|---|
| \(N_S\) | Anzahl der Einzelteilchen (Wirt). |
| \(N_T\) | Populationsgröße der Dreifach-Mutation (Wirt). |
| \(N_{SD}\) | Anzahl der einzelnen Dotieratome. |
| \(N_{TD}\) | Triplett-Population (Dotierungsmittel). |
| \(P_{HO}\) | Photonenanzahl (optischer Modus). |
| \(N_{OD}\) | Anzahl der Dotierstoffatome im Grundzustand. |
| \(N_P\) | Anzahl der Polarone. |
Kinetische Koeffizienten (Wirt)
| Symbol | Beschreibung | JSON-Token |
|---|---|---|
| \(\kappa_S\) | Abklingrate der Singulett-Population (des Wirts). | singlet_k_s |
| \(\kappa_{ISC}\) | Übergangsrate zwischen den Systemen von einem Singulett-Zustand (Gastmolekül) zu einem Triplett-Zustand (Gastmolekül). | singlet_k_isc |
| \(\kappa_{SS}\) | Annihilationsrate von Singulett-Singulett-Paaren für die Singulett-Population (Wirt). | singlet_k_ss |
| \(\kappa_{SP}\) | Annihilationsrate von Singulett-Polaronen für die Singulett-Population (Wirtsmaterial). | singlet_k_sp |
| \(\kappa_{ST}\) | Annihilationsrate zwischen Singulett- und Triplett-Zuständen, wobei sowohl der Singulett- als auch der Triplett-Zustand dem Wirtsmolekül zuzuordnen sind. | singlet_k_st |
| \(\kappa_T\) | Abklingrate der Triplett-Population (Wirt). | singlet_k_t |
| \(\kappa_{TP}\) | Annihilationsrate von Triplett- und Polaronzuständen für die Triplett-Population (Wirtsmaterial). | singlet_k_tp |
| \(\kappa_{TT}\) | Annihilationsrate für Triplett-Triplett-Übergänge bei einer Triplett-Population (Wirtsmolekül). | singlet_k_tt |
Transferkoeffizienten (Wirtsmaterial → Dotierstoff)
| Symbol | Beschreibung | JSON-Token |
|---|---|---|
| \(\kappa_{\mathrm{FRET}}\) | Förster-Energieübertragungsrate von der Singulett-Population (Wirt) zur Singulett-Population (Dotierstoff). | singlet_k_fret |
| \(\kappa_{DEXT}\) | Übertragungsrate von Dexter von der Triplett-Population (Wirt) zur Triplett-Population (Dotierstoff). | singlet_k_dext |
Kinetische Koeffizienten (Dotierstoff)
| Symbol | Beschreibung | JSON-Token |
|---|---|---|
| \(\kappa_{SD}\) | Abklingrate der Singulett-Population (Dotierstoff). | singlet_k_sd |
| \(\kappa_{ISCD}\) | Übergangsrate zwischen den Systemen von einem Singulett- (Dotierstoff-) in einen Triplett-Zustand (Dotierstoff). | singlet_k_iscd |
| \(\kappa_{SPD}\) | Annihilationsrate von Singulett-Polaronen für die Singulett-Konzentration (Dotierstoff). | singlet_k_spd |
| \(\kappa_{STD}\) | Annihilationsrate zwischen Singulett- und Triplett-Zuständen, wobei sowohl der Singulett- als auch der Triplett-Zustand durch ein Dotierungsmittel erzeugt werden. | singlet_k_std |
| \(\kappa_{SSD}\) | Annihilationsrate von Singulett-Singulett-Paaren für die Singulett-Population (Dotierstoff). | singlet_k_ssd |
| \(\kappa_{TD}\) | Abklingrate der Triplett-Zustände (Dotierstoff). | singlet_k_td |
| \(\kappa_{TTD}\) | Annihilationsrate von Triplett-Zuständen für die Triplett-Population (Dotierstoff). | singlet_k_ttd |
| \(\kappa_{TPD}\) | Annihilationsrate von Triplett-Polaronen für die Triplett-Konzentration (Dotierstoff). | singlet_k_tpd |
Optische und Hohlraumkoeffizienten
| Symbol | Beschreibung | JSON-Token |
|---|---|---|
| \(\Gamma\) | Der Einschlussfaktor beeinflusst die Kopplung zwischen der Anzahl der Singulett-Zustände (Dotierungsmittel) und der Anzahl der Photonen. | singlet_gamma |
| \(W\) | Spektraler Überlappungsfaktor zwischen der Photonendichte und der Dichte der Dotierstoffe im Grundzustand. | singlet_W |
| \(\xi\) | Gekoppelte stimulierte Emission zwischen der Population von Singulett-Zuständen (Dotierstoff) und der Population von Photonen. | — |
| \(\kappa_{CAV}\) | Verlustrate der Photonen im optischen Modus. | — |
| \(\beta_{sp}\) | Spontane Emission und ihre Wechselwirkung mit der Photonenpopulation. | — |
| \(\gamma\) | Rekombinationskoeffizient von Polaronen, der zur Bildung von Singulett- und Triplett-Zuständen führt (im Wirtsmaterial). | — |
Material- und Zusammensetzungsparameter
| JSON-Token | Symbol | Beschreibung |
|---|---|---|
singlet_C |
\(C\) | Die Dotierstoffkonzentrationsgröße wird verwendet, um die Menge des in dem emittierenden Medium vorhandenen Dotierstoffs zu definieren. |
| — | \(N_{DOP}\) | Gesamtkonzentration der Dotierstoffe. |
| — | \(P_{OD}\) | Verfügbare Dotierstoffkonzentration oder Belegungsgrad, die in die Terme für den Energietransfer vom Wirtsmaterial zum Dotierstoff einfließen. |
| — | \(n_0,\;p_0\) | Die im Ausdruck für die freie Rekombination verwendeten Gleichgewichtskonzentrationen von Elektronen und Löchern. |
Grafische Benutzeroberfläche und Steuerungstoken
| JSON-Token | Symbol | Beschreibung |
|---|---|---|
singlet_enabled |
— | Ein boolescher Schalter, der das Aktivieren oder Deaktivieren des angeregten Zustandsmodells in der grafischen Benutzeroberfläche und im Solver ermöglicht. |
text_singlet_ |
— | GUI-Beschriftungs-Token, das verwendet wird, um die Überschrift „Angeregte Zustände“ einzufügen. |