Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Method for accurate experimental determination of singlet and triplet exciton diffusion between thermally activated delayed fluorescence molecules

Marius Jakoby a, Shahriar Heidrich b, Lorenz Graf von Reventlow c, Carl Degitz d, Subeesh Madayanad Suresh d, Eli Zysman-Colman d, Wolfgang Wenzel b, Bryce S. Richards ac and Ian A. Howard *ac
aInstitute of Microstructure Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany. E-mail: ian.howard@kit.edu
bInstitute of Nanotechnology Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany
cLight Technology Institute, Karlsruhe Institute of Technology, Engesserstrasse 13, D-76131 Karlsruhe, Germany
dOrganic Semiconductor Centre, EaStCHEM School of Chemistry, University of St Andrews, St Andrews KY16 9ST, UK

Received 18th September 2020 , Accepted 13th November 2020

First published on 16th November 2020


Abstract

Understanding triplet exciton diffusion between organic thermally activated delayed fluorescence (TADF) molecules is a challenge due to the unique cycling between singlet and triplet states in these molecules. Although prompt emission quenching allows the singlet exciton diffusion properties to be determined, analogous analysis of the delayed emission quenching does not yield accurate estimations of the triplet diffusion length (because the diffusion of singlet excitons regenerated after reverse-intersystem crossing needs to be accounted for). Herein, we demonstrate how singlet and triplet diffusion lengths can be accurately determined from accessible experimental data, namely the integral prompt and delayed fluorescence. In the benchmark materials 4CzIPN and 4TCzBN, we show that the singlet diffusion lengths are (9.1 ± 0.2) and (12.8 ± 0.3) nm, whereas the triplet diffusion lengths are negligible, and certainly less than 1.0 and 1.2 nm, respectively. Theory confirms that the lack of overlap between the shielded lowest unoccupied molecular orbitals (LUMOs) hinders triplet motion between TADF chromophores in such molecular architectures. Although this cause for the suppression of triplet motion does not occur in molecular architectures that rely on electron resonance effects (e.g. DiKTa), we find that triplet diffusion is still negligible when such molecules are dispersed in a matrix material at a concentration sufficiently low to suppress aggregation. The novel and accurate method of understanding triplet diffusion in TADF molecules will allow accurate physical modeling of OLED emitter layers (especially those based on TADF donors and fluorescent acceptors).


Introduction

Thermally activated delayed fluorescence (TADF) molecules are typically characterized by a small energy gap (ΔEST) between the lowest-excited singlet and triplet excitonic states. Consequently, triplets can convert to singlets by reverse inter-system crossing (RISC).1 Such harvesting of the triplet states for emission significantly enhances the efficiency of organic light-emitting diodes (OLEDs) with red, green and blue organic TADF-based OLEDs having been realized with external quantum efficiencies exceeding 25%.2–4 In order to provide accurate models and simulations of device physics (especially in devices that utilize TADF dopants and fluorescent acceptors), the diffusion constants for singlets and triplets must be known. However, despite pioneering studies,5,6 unambiguous access to triplet diffusion lengths in these materials is lacking. The challenge in determining the diffusion lengths arises in these molecules from the excited-state population cycling multiple times through singlet and triplet states. Menke et al. have proposed a method to split the total excited-state quenching (by a quenching layer in a bilayer structure) into contributions from singlet and triplet diffusion based on a theoretical calculation of the Förster energy transfer (FRET) rate to establish the singlet contribution.5 In the context of our present results, we conclude that this FRET contribution was underestimated, leading to an overestimation of the triplet diffusion. Yurash et al. estimated the change in the rate of delayed fluorescence to be solely due to triplet quenching;6 such an estimate (see Section S4) would also lead to overestimation of triplet diffusion length.

Results and discussion

Herein, we present a novel method for accurately determining singlet and triplet diffusion lengths and constants based on easily accessible experimental data. This approach overcomes previous limitations on exciton diffusion in organic TADF molecules (2,4,5,6-tetra(9H-carbazol-9-yl)isophthalonitrile (4CzIPN)). The current work provides an accurate and experimentally accessible method to measure triplet diffusion between TADF molecules, demonstrated here using 4CzIPN, 2,3,5,6-tetrakis(3,6-di-tert-butyl-9H-carbazol-9-yl)benzonitrile (4TCzBN) and quinolino[3,2,1-de]acridine-5,9-dione (DiKTa). Our results on these systems, supported by theoretical calculations, provide the novel insight that triplet diffusion in TADF molecules can be almost negligible. Such limited triplet diffusion is highly favorable for device architectures utilizing fluorescent acceptors alongside TADF donors, and is a likely explanation for the recently observed high efficiencies in such systems.7–10

To begin, we derive the key result of this communication, eqn (9) and (11), that allows for the accurate determination of the triplet diffusion length based on time-resolved photoluminescence (PL) measurements of a series of samples with varying quencher concentrations. We start with the well-known coupled rate equations from whose solution the singlet, cs, and triplet, ct, excited-state concentrations as a function of time can be found:

 
image file: d0sc05190j-t1.tif(1)
Here, ks and kt, the singlet and triplet decay rates, are given by:
 
ks/t([Q]) = krs/t + knrs/t + kQs/t[Q] + kISC/RISC.(2)

These expressions differ from the standard formulation only in the additional term kQs/t[Q], which, as in the standard Stern–Volmer (SV) analysis, represents a quenching rate that scales linearly with the concentration of a quenching molecule distributed randomly through the film, [Q]. The prefactor kQs/t is related to the singlet/triplet diffusion constant and length. The diffusion lengths can be extracted through a SV analysis, wherein the SV constant Ks/t[Q] is defined by:

 
image file: d0sc05190j-t2.tif(3)
where k0s/t is ks/t([Q]) when Q = 0 (the unquenched singlet or triplet decay rate). Employing the Smoluchowski equation, the following relation between the SV constant Ks/t and the diffusion length for the singlet or triplet excitons ls/t can be derived:11
 
image file: d0sc05190j-t3.tif(4)
where R is the interaction radius that was chosen to be 1 nm for the used quencher TADF molecule combinations.6

The commonly used approximation for organic TADF materials that kskP (where kP is the measured rate of the prompt fluorescence decay) holds for most materials in this class and we will use it here (as kRISCks is almost always true).12 Then, the SV constant for singlets (and therefore the singlet diffusion length during a single pass through the singlet excitonic state) can be easily experimentally found by measuring the quenching of the prompt luminescence as a function of quencher concentration via the relationship:

 
image file: d0sc05190j-t4.tif(5)

The SV constant for the triplets, Kt[Q], is given by:

 
image file: d0sc05190j-t5.tif(6)

In this form, it is not immediately obvious how Kt[Q] should be easily and accurately found from experimental data. The triplet decay rate kt([Q]) (sum of the non-radiative, RISC and quenching rate) is not directly accessible by any optical experiment. However, we show in the following how Kt[Q] can be re-expressed in terms of easily observable experimental quantities.

Using the reverse and inter-system crossing efficiencies ϕRISC and ϕISC, respectively, the total PL quantum yield (PLQY) of a TADF molecule, ηtotal, can be expressed as:

 
image file: d0sc05190j-t6.tif(7)
with the prompt PLQY ηP and the average number of passes through the singlet state 〈ns〉. Solving for kt (using ϕISC/RISC = kISC/RISC/ks/t) leads to:
 
image file: d0sc05190j-t7.tif(8)
with the delayed PLQY ηD = ηtotalηP. Inserting eqn (8) in eqn (6) and again using kskP yields:
 
image file: d0sc05190j-t8.tif(9)
Here, χD = ηD/ηtotal is the fraction of delayed over the overall emission. All quantities of eqn (9) can now be measured easily using time-resolved spectroscopy methods and the triplet diffusion length in a single pass through the triplet state can be established.

Fig. 1 illustrates how time-resolved PL data as a function of quencher concentration should be analyzed to determine both the singlet and triplet diffusion lengths, based on eqn (3) and (9). In Fig. 1, the PL intensity curves are simulated with eqn (1) for kQs = 5k0s[Q] and kQt = k0t[Q], i.e. Ks/Kt = 5 (a complete set of the simulated data and other ratios Ks/Kt are shown in Fig. S3 and S4). Firstly, the data is normalized. Given the monomolecular decays, this allows the delayed and prompt lifetimes to also be found by the integral areas under the respective portions of the curve (labeled A and B). In this case we can rewrite eqn (3) as:

 
image file: d0sc05190j-t9.tif(10)
and eqn (9) as:
 
image file: d0sc05190j-t10.tif(11)


image file: d0sc05190j-f1.tif
Fig. 1 Methodology of determining singlet and triplet diffusion length between TADF molecules. (A) Simulated PL kinetics based on eqn (1) for four different quencher concentrations. (B) SV plot for singlet and triplet excitons based on eqn (10) and (11) using the indicated areas shown in (A).

Eqn (10) and (11) allow the singlet and triplet diffusion lengths (during a single pass through the respective state) to be easily determined through two integrals. The total diffusion length (considering all passes through the respective state) and the diffusion constant for the state can then also be easily found by multiplying the single pass diffusion length by the square root of the number of passes through the state, as discussed below in the experimental demonstration. As a side note, in case the prompt PL transient is not monoexponential, e.g. due to inhomogeneities in the film, an average decay rate has to be determined and eqn (9) has to be considered.

We demonstrate this method using three benchmark TADF molecules. As PL quenchers we use the electron accepting molecules indene-C60 bisadduct (ICBA) or [6,6]-phenyl-C61-butyric acid methyl ester (PCBM) that act as an electron acceptor and therefore will quench both singlet and triplet excitons with similar transfer rates. We have used an intensified charge-coupled device (ICCD) and/or a streak camera to measure the PL kinetics of the three sample sets. The organic TADF molecules were dispersed in a host matrix of 1,3-bis(N-carbazolyl)benzene (mCP) to prevent aggregation of the molecules. For 4CzIPN and 4TCzBN we used a guest concentration of 20 wt%, which is on the high side for the active layer of OLEDs involving those molecules,13 in order to determine an upper bound of the triplet diffusion length in these devices. However, due to the tendency of DiKTa to aggregate at mCP concentration of 80% we have lowered the guest concentration to 1 wt% for the DiKTa thin films.

Fig. S5 and S6 show the raw data for the time-resolved PL measurements of the TADF molecules with various quencher concentrations. The data were analyzed according to eqn (10) and (11). Since kPkD for the studied molecules, we have approximated area A and B by integration over the prompt and delayed fraction of the PL kinetic, respectively. In Fig. 2 the extracted SV plots are shown and a summary of the resulting diffusion lengths lD and diffusion constants is given in Table 1. It should be noted that the above analysis leads to the extraction of the diffusion length during one cycle through the respective exciton. That is, the diffusion lengths within the lifetimes 1/k0s and 1/k0t for singlet and triplet excitons, respectively. To determine the total length that the excited-state moves during all cycles through the singlet/triplet state, the diffusion length for a single pass in a state has to be multiplied by the square root of the average number of passes through the state. The number of passes through a state (given the excited-state enters that state once) is given by 〈n0s〉 = (A0 + B0)/A0.


image file: d0sc05190j-f2.tif
Fig. 2 SV plots of the TADF molecules 4CzIPN, 4TCzBN and DiKTa. Molecular structures are displayed in Fig. S1.
Table 1 Singlet and triplet diffusion length and constant for a single cycle and the total number of cycles through the respective state within the lifetime of the exciton. The data is extracted from the linear fits shown in Fig. 2. The determined diffusion parameters correspond to films in which the TADF molecules are dispersed in a mCP host matrix (20 wt% for 4CzIPN and 4TCzBN and 1 wt% for DiKTa)
Singlet excitons Triplet excitons
Single Total D s (cm2 s−1) Single Total D t (cm2 s−1)
l D (nm) l D (nm) l D (nm) l D (nm)
4CzIPN 8.6 ± 0.2 9.1 ± 0.2 (1.0 ± 0.1) × 10−5 <1.0 <1.2 <6 × 10−10
4TCzBN 7.3 ± 0.2 12.8 ± 0.3 (2.9 ± 0.2) × 10−5 <0.7 <1.0 <8 × 10−10
DiKTa 4.3 ± 0.2 8.3 ± 0.3 (7.3 ± 0.5) × 10−6 <0.7 <1.3 <2 × 10−11


The determined 4CzIPN (in mCP) singlet diffusion length of (9.1 ± 0.2) nm agrees well with the singlet diffusion length determined by Yurash et al. in a neat film of 4CzIPN, but there is a substantial difference for the triplet diffusion length.6 While Yurash and coworkers stated a triplet diffusion length of 2.8 nm we determined a triplet diffusion length below 1.2 nm (limited by our measurement accuracy). The smaller distance between 4CzIPN molecules in the neat film could increase the exchange coupling and the triplet diffusion length, explaining the difference with our results (which are of more relevance to a device situation wherein the TADF molecule is dispersed in a host). Another contribution may be a slight overestimate by Yurash et al. in, to our understanding, basing their triplet quenching efficiencies on the raw change in delayed emission lifetimes.6 For all studied molecules a triplet diffusion length of 0 nm is well within the error bounds of the experimental data. The upper limits in Table 1 are determined by using the upper bound of the deviation from the mean value of Kt. However, as will be shown in the following section a triplet diffusion length of essentially 0 nm is the result of our density functional theory (DFT) simulation. As a side note, the slightly negative mean of Kt for 4TCzBN could be the result of relaxation processes leading to a small overestimation of the singlet exciton quenching within the delayed lifetime. The result of a vanishing triplet diffusion length for 4CzIPN are also supported by a low temperature triplet–triplet annihilation (TTA) study by Niwa et al.14 Based on the TTA coefficient they determined a triplet diffusion constant of only about 1 × 10−13 cm2 s−1 for 10 wt% 4CzIPN in mCP (at 6 K).14

In order to compare the experimental results with theory, we employed DFT calculations to extract Dexter transfer rates based on Marcus theory. Details on these simulations can be found in the ESI (Section S8). Here, we have restricted ourselves to the two molecules 4TCzBN and DiKTa due to the structural similarity of 4TCzBN and 4CzIPN. Fig. S8 depicts the determined Dexter transfer rates as a function of distance between the interacting molecules for a film that includes only TADF molecules. The highest occurring rate in the simulation for 4TCzBN was 267 s−1 with a center-of-mass distance between the interacting molecules of 1.4 nm. Assuming that between each neighboring molecule this transfer rate would apply, an upper bound of the triplet diffusion constant can be determined to ∼1 × 10−12 cm2 s−1. For the chosen 20 wt% concentration that is also used in devices, the average distance between molecules increases to about 2 nm. The calculated charge-transfer rates reveal that electron transfer is 3 to 4 orders of magnitude lower than the transfer of the hole (compare Fig. S10(B)), suggesting the restricted LUMO–LUMO overlap is responsible for the negligible triplet transport in these molecules. Whereas the triplet transport for DiKTa can be significantly faster at short center of mass spacing, at spacings similar to those found in devices the triplet transport is equally low. The mechanisms for the limited triplet transport in these molecules are schematically illustrated in Fig. 3.


image file: d0sc05190j-f3.tif
Fig. 3 Schematic of exciton motion in TADF based devices. (A) Suppressed triplet transfer as a result of the molecular architecture restricting the LUMO–LUMO overlap for two donor–acceptor type molecules like 4CzIPN. (B) Suppressed triplet transfer due to a highly dispersion (to prevent aggregation) for electron resonance based TADF molecules.

The negligible triplet diffusion lengths for donor–acceptor TADF molecules like 4CzIPN and 4TCzBN are supported by recent results of Franco et al. where attaching a phenylene ethynylene oligomer to a 4CzIPN core did not lead to quenching of the TADF charge-transfer triplet state by the lower-lying triplet state on the oligomer.15 Along this line, a recent comparison of experimental device characteristics and kinetic Monte Carlo simulations came to the conclusion that in a TADF donor-fluorescent acceptor device based on 4CzIPN-Me triplet transport to fluorescent dopants is a negligible loss channel.16 Furthermore, Narushima et al. recently showed in the context of conjugated molecular crystals that even though the stacked molecules have a good LUMO–LUMO overlap the very small overlap between adjacent highest occupied molecular orbitals (HOMOs) leads to a very small diffusion constant for triplet excitons (3 × 10−9 cm2 s−1).17

Conclusions

In conclusion, we provide an easy and accurate method of measuring triplet diffusion in TADF molecules. This method overcomes the unique challenge posed in TADF systems of the delayed fluorescence lifetime being quenched both by the motion of triplets, and regenerated singlets. We establish that the triplet motion can be negligible compared to the singlet motion in TADF molecules, a situation that is highly favorable for device applications requiring transfer of singlets to a fluorescence acceptor.18 Given this method relies only on standard experimental data, we anticipate that it will be of broad use for increasing the fundamental understanding of TADF materials, and obtaining accurate parameters to allow physical modeling of device performance.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

The authors thank the German Research Foundation (DFG) for financial support through SFB 1176 (Project A4) and Helmholtz Energy Materials Foundry HEMF funding of the Helmholtz Association. M. J. acknowledges support from the Karlsruhe School of Optics & Photonics (KSOP) graduate school. S. M. S. thanks the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement no 838885 (NarrowbandSSL). The St Andrews team thanks the EPSRC for financial support (EP/P010482/1).

Notes and references

  1. F. B. Dias, T. J. Penfold and A. P. Monkman, Methods Appl. Fluoresc., 2017, 5, 012001 CrossRef .
  2. W. Zeng, H.-Y. Lai, W.-K. Lee, M. Jiao, Y.-J. Shiu, C. Zhong, S. Gong, T. Zhou, G. Xie, M. Sarma, K.-T. Wong, C.-C. Wu and C. Yang, Adv. Mater., 2018, 30, 1704961 CrossRef .
  3. B. S. Kim and J. Y. Lee, Adv. Funct. Mater., 2014, 24, 3970–3977 CrossRef CAS .
  4. T.-A. Lin, T. Chatterjee, W.-L. Tsai, W.-K. Lee, M.-J. Wu, M. Jiao, K.-C. Pan, C.-L. Yi, C.-L. Chung, K.-T. Wong and C.-C. Wu, Adv. Mater., 2016, 28, 6976–6983 CrossRef CAS .
  5. S. M. Menke and R. J. Holmes, J. Phys. Chem. C, 2016, 120, 8502–8508 CrossRef CAS .
  6. B. Yurash, H. Nakanotani, Y. Olivier, D. Beljonne, C. Adachi and T.-Q. Nguyen, Adv. Mater., 2019, 31, 1804490 CrossRef .
  7. H. Nakanotani, T. Higuchi, T. Furukawa, K. Masui, K. Morimoto, M. Numata, H. Tanaka, Y. Sagara, T. Yasuda and C. Adachi, Nat. Commun., 2014, 5, 4016 CrossRef CAS .
  8. S. H. Han, J. H. Jeong, J. W. Yoo and J. Y. Lee, J. Mater. Chem. C, 2019, 7, 3082–3089 RSC .
  9. S. Wang, Y. Zhang, W. Chen, J. Wei, Y. Liu and Y. Wang, Chem. Commun., 2015, 51, 11972–11975 RSC .
  10. X.-K. Liu, Z. Chen, C.-J. Zheng, M. Chen, W. Liu, X.-H. Zhang and C.-S. Lee, Adv. Mater., 2015, 27, 2025–2030 CrossRef CAS .
  11. J. R. Lakowicz, Principles of Fluorescence Spectroscopy, Springer, 2011 Search PubMed .
  12. C. Baleizão and M. N. Berberan-Santos, J. Chem. Phys., 2007, 126, 204510 CrossRef .
  13. B. Wex and B. R. Kaafarani, J. Mater. Chem. C, 2017, 5, 8622–8653 RSC .
  14. A. Niwa, S. Haseyama, T. Kobayashi, T. Nagase, K. Goushi, C. Adachi and H. Naito, Appl. Phys. Lett., 2018, 113, 083301 CrossRef .
  15. O. Franco, M. Jakoby, R. V. Schneider, F. Hundemer, D. Hahn, B. S. Richards, S. Bräse, M. A. R. Meier, U. Lemmer and I. A. Howard, Front. Chem., 2020, 8, 126 CrossRef CAS .
  16. S. Gottardi, M. Barbry, R. Coehoorn and H. van Eersel, Appl. Phys. Lett., 2019, 114, 073301 CrossRef .
  17. K. Narushima, Y. Kiyota, T. Mori, S. Hirata and M. Vacha, Adv. Mater., 2019, 31, 1807268 CrossRef .
  18. M. Jakoby, B. S. Richards, U. Lemmer and I. A. Howard, Phys. Rev. B, 2019, 100, 045303 CrossRef CAS .

Footnotes

The research data supporting this publication can be accessed at https://doi.org/10.17630/883f41b6-0a9f-4867-852b-db0940b6220e.
Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc05190j

This journal is © The Royal Society of Chemistry 2021
Click here to see how this site uses Cookies. View our privacy policy here.