Virtual Screening of TADF Emitters for Single-Layer OLEDs

Abstract

Thermally-activated delayed fluorescence (TADF) is a concept which helps to harvest triplet excitations, boosting the efficiency of an organic light-emitting diode. TADF can be observed in molecules with spatially separated donor and acceptor groups with a reduced triplet-singlet energy level splitting. TADF materials with balanced electron and hole transport are attractive for realizing efficient single-layer organic light emitting diodes, greatly simplifying their manufacturing and improving their stability. Our goal here is to computationally screen such materials and provide a comprehensive database of compounds with a range of emission wavelengths, ionization energies, and electron affinities.

Introduction

For obtaining efficient organic light-emitting diodes (OLEDs), it is convenient to tune individual processes, such as charge injection, balanced hole and electron transport, and triplet and singlet exciton harvesting, by using dedicated layers. Every new material adds a degree of freedom and hence flexibility to the OLED design. For instance, doped charge transport layers ensure Ohmic injection, an appropriate host material balances transport inside the emitting layer, and the phosphorescent emitter ensures triplet harvesting. However, every new emitter requires optimization of the surrounding layers, with respect to energy levels, triplet energies, and charge-transport properties, complicating the OLED design.

Recently, it was demonstrated that a complex multilayer design can be substituted by a simple single-layer architecture (Kotadiya et al., 2019a) without compromising the balanced and trap-free electron and hole transport. The ohmic charge injection and the absence of heterojunctions resulted in extremely low operating voltages and thus power efficiency in a single-layer OLED utilizing thermally activated delayed fluorescence, which helps to convert triplet into singlet excitons (Uoyama et al., 2012; Godumala et al., 2019)^. An external quantum efficiency of 19% was achieved. Owing to the broad recombination zone and low operating voltages, one of the key features of the single-layer device is the improved device stability, which can be used to design a stable blue OLED, a grand challenge in OLED research (Heimel et al., 2018; Paterson et al., 2019, 2020). In view of this, it would be useful to understand if the single-layer design can be employed for blue OLEDs: the issue here is the trap-free transport for both holes and electrons, which sets limits on the transport gap. In this paper, we first formulate the chemical design rules for TADF emitters with ambipolar transport. Using these rules, we then computationally pre-screen a set of molecules comprised of acceptor, donor, and bridge blocks and grade them according to the predicted emission wavelength.

Design Criteria

Singlet-Triplet Energy Splitting

The important task of a TADF emitter is to convert triplet into singlet excitations. To do this, the reverse intersystem crossing rate, , should be high, which is only possible if the energy difference between the first singlet and the first triplet excited state is small, . A typical example of a TADF emitter is CzDBA (Wu et al., 2018), shown in Figure 1. CzDBA has a D--A--D architecture: two carbazole (Cz) fragments, two m-xylene bridges and a central 5,10-dihydroboranthrene (DBA) core. The methyl groups on the m-xylene bridge ensure that the core unit is nearly orthogonal to the bridge, leading to a small overlap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) and hence nearly zero .

FIGURE 1

Ambipolar Trap-free Transport

To ensure a broad recombination zone within the emission layer, the thin film of the TADF emitter should provide balanced and trap-free transport of holes and electrons. To realize this, one needs to select compounds with an ionization energy (IE) and electron affinity (EA) lying within the trap-free energy window (Kotadiya et al., 2019b), i.e., with ionization energy (IE) < 6.5 eV and electron affinity (EA) > 2.5 eV. These criteria ensure that contaminants such as oxygen or water do not serve as energetic traps for holes and electrons.

Small Energetic Disorder

From a dipolar glass model, the energetic disorder present in a disordered molecular solid is proportional to the dipole moment of the composing molecule. Therefore, thin organic films with molecules with a small dipole moment (D) normally have a narrower density of states (Novikov and Vannikov, 2009; Lin et al., 2019; Mondal et al., 2021; Stankevych et al., 2021). This design criteria can be enforced by selecting centrosymmetric molecules only of the D--A--D or A--D--A type, similar to CzDBA. This molecular architecture ensures a small dipole moment and hence narrow density of states (Liu et al., 2021).

Building Blocks

With these design rules in mind, and in view of the successful example of CzDBA, we build and characterize a database of emitters that fulfill the aforementioned criteria. To construct the emitters, we start with 97 potential donor and acceptor building blocks, all shown in the Supplementary Note S1. All of them are (quasi-)linear, composed of three (fused) rings and are reported in literature (synthesizable). These building blocks are further pre-screened to ensure the desired donor-acceptor architecture in an emitter. The pre-screening proceeds as follows: knowing that the IE and EA of CzDBA is already quite close to the boundary of the trap-free window (Kotadiya et al., 2019a; Liu et al., 2021) we take the IE_Cz and EA_DBA as the pre-screening criteria for donors and acceptors, respectively. Only the fragments possessing IE < IE_Cz + 0.2 eV (EA > EA_DBA - 0.2 eV) will be chosen as “trap-free” donors (acceptors) and enter the next round, see Supplementary Figure S2. The molecular structures of donors and acceptors that pass the prescreening step are summarized in Figure 2. To build the emitter molecules, only the building blocks with the inversion symmetry are used as core fragments. These are shown in dark colors in Figure 2. This choice helps to fulfill the centrosymmetric requirement for the entire molecule.

FIGURE 2

Computational Workflow

Using the selected building blocks, we constructed the database of D--A--D and A--D--A. The simplified molecular-input-line-entry system (SMILES) strings of compounds were created through combination of the SMILES strings of the composing donor, bridge and acceptor. The initial geometry of each compound was first optimized using a semi-empirical method and then by density functional theory (DFT). Details are given in the Supplementary Note S2.

To obtain reliable predictions of solid-state IE, EA and excited-state energy, we followed the cost-effective -tuning protocol (Sun et al., 2016, 2017). In addition to the , the difference in the characters of the singlet and triplet excited states are crucial to the rISC rate (El-Sayed, 1963). For this reason, the excited-state characters were evaluated using a fragment-based method (Plasser, 2020).

For compounds that pass the screening criteria, the density-of-states distributions for holes and electrons were computed via multi-scale simulations, that include morphology generation using molecular-dynamics simulations, followed by polarizable force-field evaluation of the solid-state contributions to the gas-phase energy levels (Rühle et al., 2011; Poelking and Andrienko, 2016; Andrienko, 2018; Mondal et al., 2021). The entire workflow is illustrated in Figure 3.

FIGURE 3

Results and Discussion

Compounds With Small Singlet-Triplet Splitting

The combination of the core and the arm fragments gives in total 441 A--D--A and 504 D--A--D compounds. Due to convergence problems in geometry optimization, especially in the anionic state with implicit solvent, the final database contained 433 A--D--A and 481 D--A--D compounds.

The IE and EA of all compounds either lies within the “trap-free window” or close to the borderline of the window, showing that the effectiveness of prescreening of the building blocks. Therefore, we put our emphasis on the small criterion. The distributions of the E_S1 and ΔE_ST are shown in Figure 4 (A--D--A) and Supplementary Figure S2 (D--A--D). Around 50% of the compounds (206 out of 433 for A--D--A and 268 out of 481 for D--A--D) have very small singlet-triplet energy level splitting, < 0.1 eV, which illustrates the efficiency of the design strategy, that is the use of the m-xylene bridge. Moreover, the computed energy and of CzDBA is 2.487 and 0.016 eV, which is in excellent agreement with the experimental values of 2.48 and 0.033 eV (Wu et al., 2018; Kotadiya et al., 2019a).

FIGURE 4

Among these small- compounds, we observed a broad distribution in the energy, ranging from 0.2 to 2.9 eV. This indicates the opportunity to design single-layer emitting OLEDs of different colors, including the infrared region. The two branches in Figure 4 represent the rest (50%) of the emitters with > 0.1 eV, where a similar branch is also observed for D--A--D (Supplementary Figure S3). This is counterintuitive as the should be small if the HOMO and the LUMO are separated via the m-xylene bridges.

Analysis of the Excited-State Character

To better understand the origin of the large , we calculated the charge transfer (CT) number ranging from 0 to 1, using the fragment-based analysis (see Supplementary Note S2). We define the core as one fragment (f_C) and two bridge + arm pairs as the other fragment (f_A). If the hole is 100% located at one fragment and the electron is 100% located at the other one, the charge transfer number is 1, representing a 100% CT character. In contrast, if the hole and the electron are both localized on the same fragment, the CT number is 0, featuring a local-excitation (LE) character. In most cases, the CT number is a fraction between 0 and 1 since most adiabatic excited states exhibit a mixture of CT and LE characters. The larger the CT number of the excited state is, the higher the CT character it has.

Figure 5 depicts the 2D histogram based on the CT numbers of T₁ and S₁ states for the 227 A--D--A compounds with > 0.1 eV. Most of the scatter points are located at the upper left corner, meaning that these emitters possess a charge-transfer S₁ and locally-excited T₁ states. A similar result was observed in the D--A--D case (see Supplementary Figure S4).

FIGURE 5

The emergence of the LE states can be explained utilizing the frontier molecular orbital (FMO) energies of the constituent building blocks (Blaskovits et al., 2020), which is illustrated in Figure 6A. The competition between the CT excitation and LE excitation depends on the relative ordering of the FMOs. In this context, we can define two descriptors, , , where are the LUMO/HOMO energies of the acceptor/donor. If the or the is much larger than 1, the CT excitation is more favorable than the LE for the low-lying excited states and vice versa. Figure 6B demonstrates that this simple approximation works quite well for our A--D--A database: For or smaller than ∼1.2, the CT number of T₁ becomes close to 0. The same behavior was also found in the D--A--D database, as shown in Supplementary Figure S5. This indicates that a prescreening step based on the individual building blocks saves the computational cost, similar to pre-screening of singlet fission donor-acceptor copolymers (Blaskovits et al., 2020).

FIGURE 6

Compounds With T_n States Lying Close to S₁

The S₁ and T₁ states of most molecules that pass the first screening step ( < 0.1 eV) exhibit CT character. According to the El-Sayed rule, the is zero between two states having the same excited-state character, which implies that the rISC may not occur for these pre-screened compounds. However, the conformational disorder present in the solid state leads to a distribution of dihedral angles between the constituent donor and acceptor (Weissenseel et al., 2019). This disorder gives rise to different excited-state characters, that is different mixing of CT and LE diabatic states of the S₁ and T₁ states, (de Silva et al., 2019), resulting in non-zero . This explains why TADF could still be observed in the thin film of CzDBA, where CT_S1 and CT_T1 are both close to 1 in the gas phase (Kotadiya et al., 2019a).

In addition, higher triplet states (T_n with ) with different excited-state character from that of S₁, can also assist in the rISC process via a two-step mechanism (Gibson et al., 2016). A large second order coupling can be achieved when the energies of S₁, T₁ and T_n are close to each other. Compounds with close-lying S₁ and T₁ that already show different excited-state characters would possess large first-order coupling and hence high k_rISC. Therefore, we applied additional screening criteria to the as-screened ∼500 molecules: 1) there should be at least one triplet state T_n that is close to S₁ (| E_S1− E_Tn | < 0.1 eV); 2) for the triplet states that are energetically close to S₁, the difference between the CT numbers of S₁ and T_n should be larger than 0.5 (CT_S1−CT_Tn > 0.5) to give reasonable spin-orbit coupling.

Overall, around 100 molecules pass the criteria (49 A--D--A and 46 D--A--D), where the molecular structures are summarized in Figure 7. All of these compounds, except for A--D--A molecules with non-centrosymmetric core boant4 ( Debye), possess nearly zero molecular dipole moment. Therefore, they are considered promising candidates for single-layer OLED emitters. The position of the electroluminescence (EL) spectrum maximum of each compound, as shown in Figure 7, was estimated by subtracting the computed S₁ energy by a value δ, which is defined as , where E_S1 and λ_EL,max is the experimental optical gap and the wavelength of the EL spectrum maximum of CzDBA (Kotadiya et al., 2019a). These values are listed in Supplementary Tables S1,S2. We obtained a series of potential TADF emitters with various EL spectrum maximum, ranging from infrared (0.716 eV) to blue color (2.660 eV), which paves the way for future development of single-layer OLED devices.

FIGURE 7

Charge Carrier Density of States

More sophisticated solid-state simulations can then be performed for the much smaller molecular dataset, which is now only ∼10% of the initial number of compounds. As a proof of concept, we computed the charge carrier density of states for the blue A--D--A emitter, 37bdt1-ant2 (as shown in Figure 8). The amorphous simulated morphology was generated using molecular dynamics, where the details can be found in Supplementary Note S4. The energetic disorder for electrons (0.11 eV) and holes (0.12 eV) is relatively small, which indicates a good hole/electron mobility. This also demonstrates the success of our design strategy regarding small molecular dipole moment. Since the simulated IE and EA of 37bdt1-ant2 lie at the border of the trap-free window, further experimental measurements are necessary to verify if it is really free from universal traps.

FIGURE 8

Conclusion

Using the suggested design rules, we have proposed a set of TADF emitters with a broad range of emission wavelengths, from infrared to sky-blue. We hope that the suggested structures can serve as a clear guide towards further development of efficient and stable single-layer OLEDs.

References

• 1

  Andrienko D. (2018). “Multiscale Concepts in Simulations of Organic Semiconductors,” in Handbook of Materials Modeling : Methods: Theory and Modeling. Editors Andreoni,W.YipS. (Cham: Springer International Publishing), 1–12. 10.1007/978-3-319-42913-7_39-1

  • CrossRef
  • Google Scholar

• 2

  Blaskovits J. T. Fumanal M. Vela S. Corminboeuf C. (2020). Designing Singlet Fission Candidates from Donor-Acceptor Copolymers. Chem. Mater.32, 6515–6524. 10.1021/acs.chemmater.0c01784

  • CrossRef
  • Google Scholar

• 3

  de Silva P. Kim C. A. Zhu T. Van Voorhis T. (2019). Extracting Design Principles for Efficient Thermally Activated Delayed Fluorescence (TADF) from a Simple Four-State Model. Chem. Mater.31, 6995–7006. 10.1021/acs.chemmater.9b01601

  • CrossRef
  • Google Scholar

• 4

  El‐Sayed M. A. (1963). Spin-Orbit Coupling and the Radiationless Processes in Nitrogen Heterocyclics. J. Chem. Phys.38, 2834–2838. 10.1063/1.1733610

  • CrossRef
  • Google Scholar

• 5

  Gibson J. Monkman A. P. Penfold T. J. (2016). The Importance of Vibronic Coupling for Efficient Reverse Intersystem Crossing in Thermally Activated Delayed Fluorescence Molecules. ChemPhysChem17, 2956–2961. 10.1002/cphc.201600662

  • CrossRef
  • Google Scholar

• 6

  Godumala M. Choi S. Cho M. J. Choi D. H. (2019). Recent Breakthroughs in Thermally Activated Delayed Fluorescence Organic Light Emitting Diodes Containing Non-Doped Emitting Layers. J. Mater. Chem. C7, 2172–2198. 10.1039/C8TC06293E

  • CrossRef
  • Google Scholar

• 7

  Heimel P. Mondal A. May F. Kowalsky W. Lennartz C. Andrienko D. et al (2018). Unicolored Phosphor-Sensitized Fluorescence for Efficient and Stable Blue OLEDs. Nat. Commun.9, 4990. 10.1038/s41467-018-07432-2

  • CrossRef
  • Google Scholar

• 8

  Kotadiya N. B. Blom P. W. M. Wetzelaer G.-J. A. H. (2019a). Efficient and Stable Single-Layer Organic Light-Emitting Diodes Based on Thermally Activated Delayed Fluorescence. Nat. Photon.13, 765–769. 10.1038/s41566-019-0488-1

  • CrossRef
  • Google Scholar

• 9

  Kotadiya N. B. Mondal A. Blom P. W. M. Andrienko D. Wetzelaer G.-J. A. H. (2019b). A Window to Trap-Free Charge Transport in Organic Semiconducting Thin Films. Nat. Mater.18, 1182–1186. 10.1038/s41563-019-0473-6

  • CrossRef
  • Google Scholar

• 10

  Lin K.-H. Prlj A. Yao L. Drigo N. Cho H.-H. Nazeeruddin M. K. et al (2019). Multiarm and Substituent Effects on Charge Transport of Organic Hole Transport Materials. Chem. Mater.31, 6605–6614. 10.1021/acs.chemmater.9b00438

  • CrossRef
  • Google Scholar

• 11

  Liu W. Kotadiya N. B. Blom P. W. M. Wetzelaer G. J. A. H. Andrienko D. (2021). Molecular Origin of Balanced Bipolar Transport in Neat Layers of the Emitter CzDBA. Adv. Mater. Technol.6, 2000120. 10.1002/admt.202000120

  • CrossRef
  • Google Scholar

• 12

  Mondal A. Paterson L. Cho J. Lin K.-H. van der Zee B. Wetzelaer G.-J. A. H. et al (2021). Molecular Library of OLED Host Materials-Evaluating the Multiscale Simulation Workflow. Chem. Phys. Rev.2, 031304. 10.1063/5.0049513

  • CrossRef
  • Google Scholar

• 13

  Novikov S. V. Vannikov A. V. (2009). Hopping Charge Transport in Disordered Organic Materials: Where Is the Disorder?J. Phys. Chem. C113, 2532–2540. 10.1021/jp808578b

  • CrossRef
  • Google Scholar

• 14

  Paterson L. Mondal A. Heimel P. Lovrincic R. May F. Lennartz C. et al (2019). Perspectives of Unicolored Phosphor‐Sensitized Fluorescence. Adv. Electron. Mater.5, 1900646. 10.1002/aelm.201900646

  • CrossRef
  • Google Scholar

• 15

  Paterson L. May F. Andrienko D. (2020). Computer Aided Design of Stable and Efficient OLEDs. J. Appl. Phys.128, 160901. 10.1063/5.0022870

  • CrossRef
  • Google Scholar

• 16

  Plasser F. (2020). TheoDORE: A Toolbox for a Detailed and Automated Analysis of Electronic Excited State Computations. J. Chem. Phys.152, 084108. 10.1063/1.5143076

  • CrossRef
  • Google Scholar

• 17

  Poelking C. Andrienko D. (2016). Long-Range Embedding of Molecular Ions and Excitations in a Polarizable Molecular Environment. J. Chem. Theor. Comput.12, 4516–4523. 10.1021/acs.jctc.6b00599

  • CrossRef
  • Google Scholar

• 18

  Rühle V. Lukyanov A. May F. Schrader M. Vehoff T. Kirkpatrick J. et al (2011). Microscopic Simulations of Charge Transport in Disordered Organic Semiconductors. J. Chem. Theor. Comput.7, 3335–3345. 10.1021/ct200388s

  • CrossRef
  • Google Scholar

• 19

  Stankevych A. Vakhnin A. Andrienko D. Paterson L. Genoe J. Fishchuk I. et al (2021). Density of States of OLED Host Materials from Thermally Stimulated Luminescence. Phys. Rev. Appl.15, 044050. 10.1103/PhysRevApplied.15.044050

  • CrossRef
  • Google Scholar

• 20

  Sun H. Ryno S. Zhong C. Ravva M. K. Sun Z. Körzdörfer T. et al (2016). Ionization Energies, Electron Affinities, and Polarization Energies of Organic Molecular Crystals: Quantitative Estimations from a Polarizable Continuum Model (PCM)-Tuned Range-Separated Density Functional Approach. J. Chem. Theor. Comput.12, 2906–2916. 10.1021/acs.jctc.6b00225

  • CrossRef
  • Google Scholar

• 21

  Sun H. Hu Z. Zhong C. Chen X. Sun Z. Brédas J.-L. (2017). Impact of Dielectric Constant on the Singlet-Triplet Gap in Thermally Activated Delayed Fluorescence Materials. J. Phys. Chem. Lett.8, 2393–2398. 10.1021/acs.jpclett.7b00688

  • CrossRef
  • Google Scholar

• 22

  Uoyama H. Goushi K. Shizu K. Nomura H. Adachi C. (2012). Highly Efficient Organic Light-Emitting Diodes from Delayed Fluorescence. Nature492, 234–238. 10.1038/nature11687

  • CrossRef
  • Google Scholar

• 23

  Weissenseel S. Drigo N. A. Kudriashova L. G. Schmid M. Morgenstern T. Lin K.-H. et al (2019). Getting the Right Twist: Influence of Donor-Acceptor Dihedral Angle on Exciton Kinetics and Singlet-Triplet Gap in Deep Blue Thermally Activated Delayed Fluorescence Emitter. J. Phys. Chem. C123, 27778–27784. 10.1021/acs.jpcc.9b08269

  • CrossRef
  • Google Scholar

• 24

  Wu T.-L. Huang M.-J. Lin C.-C. Huang P.-Y. Chou T.-Y. Chen-Cheng R.-W. et al (2018). Diboron Compound-Based Organic Light-Emitting Diodes with High Efficiency and Reduced Efficiency Roll-Off. Nat. Photon12, 235–240. 10.1038/s41566-018-0112-9

  • CrossRef
  • Google Scholar