Functionalization of Clar's Goblet Diradical with Heteroatoms: Tuning the Excited-State Energies to Promote Triplet-to-Singlet Conversion

The ground-state spin multiplicity as well as the energy difference between the lowest-energy spin-singlet (S1) and spin-triplet (T1) excited states of topologically frustrated organic (diradical) molecules can be tuned by doping with a pair of heteroatoms (N or B atoms). We have thus systematically studied here a set of Clar's Goblet derivatives upon a controlled substitution at different C sites, to alter the electronic structure of the molecules and disclose the positions at which: (i) the ground-state multiplicity becomes a closed-shell singlet and (ii) the energy difference between S1 and T1 is considerably small (i.e., below 0.1-0.2 eV to induce a triplet exciton recovery upon thermal effects). This electronic structure outcome is driven by strong correlation effects; therefore, we have here applied a variety of single-reference [TD-DFT, CIS(D), SCS-CC2] and multireference [CASSCF, NEVPT2, RAS-srDFT] methods. For TD-DFT, we have covered global hybrid (PBE0, M06-2X), range-separated hybrid (ωB97X), and double-hybrid (PBE-QIDH, SOS1-PBE-QIDH, and PBE0-2) functionals to ascertain whether the results were highly dependent on the functional choice. Overall, we found that the heterosubstitution strategy could largely modify the electronic and optical properties of the pristine diradical system, with these organic forms thus constituting a new set of compounds with further optoelectronic applications.

The Golden Age of Thermally Activated Delayed Fluorescence Materials: Design and Exploitation

Since the seminal report by Adachi and co-workers in 2012, there has been a veritable explosion of interest in the design of thermally activated delayed fluorescence (TADF) compounds, particularly as emitters for organic light-emitting diodes (OLEDs). With rapid advancements and innovation in materials design, the efficiencies of TADF OLEDs for each of the primary color points as well as for white devices now rival those of state-of-the-art phosphorescent emitters. Beyond electroluminescent devices, TADF compounds have also found increasing utility and applications in numerous related fields, from photocatalysis, to sensing, to imaging and beyond. Following from our previous review in 2017 ( Adv. Mater. 2017, 1605444), we here comprehensively document subsequent advances made in TADF materials design and their uses from 2017-2022. Correlations highlighted between structure and properties as well as detailed comparisons and analyses should assist future TADF materials development. The necessarily broadened breadth and scope of this review attests to the bustling activity in this field. We note that the rapidly expanding and accelerating research activity in TADF material development is indicative of a field that has reached adolescence, with an exciting maturity still yet to come.

Reduced Scaling Correlated Natural Transition Orbitals for Multilevel Coupled Cluster Calculations

Multilevel coupled cluster theory offers reduced scaling computation of intensive properties in systems that are too large for standard coupled cluster calculations. A significant benefit of the multilevel coupled cluster framework is the possibility of calculating intensive properties that are not tightly localized if an appropriate set of active orbitals is used. Correlated natural transition orbitals (CNTOs) are tailored to describe excitation processes. For multilevel coupled cluster singles and doubles (MLCCSD) and singles and perturbative doubles (MLCC2) calculations, the construction of CNTOs generally becomes the computational bottleneck. Here, we demonstrate how CNTOs can be obtained with O ( N 3 ) operations, eliminating the O ( N 5 ) -scaling steps involved in the original approach. This reduction in scaling moves the bottleneck of MLCC2 and MLCCSD calculations from the active orbital space preparation to the MLCC2 and MLCCSD equations with O ( N 4 ) -scaling.

Performance of Density Functionals for Excited-State Properties of Isolated Chromophores and Exciplexes: Emission Spectra, Solvatochromic Shifts, and Charge-Transfer Character

This study assesses the performance of various meta-generalized gradient approximation (meta-GGA), global hybrid, and range-separated hybrid (RSH) density functionals in capturing the excited-state properties of organic chromophores and their excited-state complexes (exciplexes). Motivated by their uses in solar energy harvesting and photoredox CO2 reduction, we use oligo-(p-phenylenes) and their excited-state complexes with triethylamine as model systems. We focus on the fluorescence properties of these systems, specifically emission energies. We also consider solvatochromic shifts and wave function characteristics. The latter is described by using reduced quantities such as natural transition orbitals (NTOs) and exciton descriptors. The functionals are benchmarked against the experimental fluorescence spectra and the equation-of-motion coupled-cluster method with single and double excitations. Both in isolated chromophores and in exciplexes, meta-GGA functionals drastically underestimate the emission energies and exhibit significant exciton delocalization and anticorrelation between electron and hole motion. The performance of global hybrid functionals is strongly dependent on the percentage of exact exchange. Our study identifies RSH GGAs as the best-performing functionals, with ωPBE demonstrating the best agreement with experimental results. RSH meta-GGAs often overestimate emission energies in exciplexes and yield larger hole NTOs. Their performance can be improved by optimally tuning the range-separation parameter.

Triplet Excited States with Multilevel Coupled Cluster Theory

We extend the multilevel coupled cluster framework with triplet excitation energies at the singles and perturbative doubles (MLCC2) and singles and doubles (MLCCSD) levels of theory. In multilevel coupled cluster theory, we partition the orbitals and restrict the higher-order excitations in the cluster operator to a set of active orbitals. With an appropriate choice of these orbitals, the multilevel approach can give significant computational savings while maintaining the high accuracy of standard coupled cluster theory. In this work, we generated active orbitals from approximate correlated natural transition orbitals (CNTOs). The CNTOs form a compact orbital space specifically tailored to describe the triplet excited states of interest. We compare the performance of MLCCSD and MLCC2, in terms of cost and accuracy, to those of their standard coupled cluster counterparts (CC2 and CCSD) and finally show proof-of-concept calculations of the singlet-triplet gaps of molecules that are of interest for their potential use in organic light-emitting diodes.

Intramolecular locking and coumarin insertion: a stepwise approach for TADF design

Three novel TADF (thermally activated delayed fluorescence) emitters based on the well-studied Qx-Ph-DMAC fluorophore are designed and synthesized. The photophysical properties of these materials are studied from a theoretical and experimental point of view, demonstrating the cumulative effects of multiple small modifications that combine to afford significantly improved TADF performance. First, an extra phenyl ring is added to the acceptor part of Qx-Ph-DMAC to increase the conjugation length, resulting in BQx-Ph-DMAC, which acts as an intermediate molecular structure. Next, an electron-deficient coumarin unit is incorporated to fortify the electron accepting ability, affording ChromPy-Ph-DMAC with red-shifted emission. Finally, the conjugated system is further enlarged by 'locking' the molecular structure, generating DBChromQx-DMAC with further red-shifted emission. The addition of the coumarin unit significantly impacts the charge-transfer excited state energy levels with little effect on the locally excited states, resulting in a decrease of the singlet-triplet energy gap. As a result, the two coumarin-based emitters show considerably improved TADF performance in 1 w/w% zeonex films when compared to the initial Qx-Ph-DMAC structure. 'Locking' the molecular structure further lowers the singlet-triplet energy gap, resulting in more efficient reverse intersystem crossing and increasing the contribution of TADF to the total emission.

Benchmarking DFT Functionals for Excited-State Calculations of Donor-Acceptor TADF Emitters: Insights on the Key Parameters Determining Reverse Inter-System Crossing

The importance of intermediate triplet states and the nature of excited states has gained interest in recent years for the thermally activated delayed fluorescence (TADF) mechanism. It is widely accepted that simple conversion between charge transfer (CT) triplet and singlet excited states is too crude, and a more complex route involving higher-lying locally excited triplet excited states has to be invoked to witness the magnitude of the rate of reverse inter-system crossing (RISC) rates. The increased complexity has challenged the reliability of computational methods to accurately predict the relative energy between excited states as well as their nature. Here, we compare the results of widely used density functional theory (DFT) functionals, CAM-B3LYP, LC-ωPBE, LC-ω*PBE, LC-ω*HPBE, B3LYP, PBE0, and M06-2X, against a wavefunction-based reference method, Spin-Component Scaling second-order approximate Coupled Cluster (SCS-CC2), in 14 known TADF emitters possessing a diversity of chemical structures. Overall, the use of the Tamm-Dancoff Approximation (TDA) together with CAM-B3LYP, M06-2X, and the two ω-tuned range-separated functionals LC-ω*PBE and LC-ω*HPBE demonstrated the best agreement with SCS-CC2 calculations in predicting the absolute energy of the singlet S₁, and triplet T₁ and T₂ excited states and their energy differences. However, consistently across the series and irrespective of the functional or the use of TDA, the nature of T₁ and T₂ is not as accurately captured as compared to S₁. We also investigated the impact of the optimization of S₁ and T₁ excited states on ΔE_ST and the nature of these states for three different functionals (PBE0, CAM-B3LYP, and M06-2X). We observed large changes in ΔE_ST using CAM-B3LYP and PBE0 functionals associated with a large stabilization of T₁ with CAM-B3LYP and a large stabilization of S₁ with PBE0, while ΔE_ST is much less affected considering the M06-2X functional. The nature of the S₁ state barely evolves after geometry optimization essentially because this state is CT by nature for the three functionals tested. However, the prediction of the T₁ nature is more problematic since these functionals for some compounds interpret the nature of T₁ very differently. SCS-CC2 calculations on top of the TDA-DFT optimized geometries lead to a large variation in terms of ΔE_ST and the excited-state nature depending on the chosen functionals, further stressing the large dependence of the excited-state features on the excited-state geometries. The presented work highlights that despite good agreement of energies, the description of the exact nature of the triplet states should be undertaken with caution.

Modeling of Multiresonant Thermally Activated Delayed Fluorescence Emitters─Properly Accounting for Electron Correlation Is Key!

With the surge of interest in multiresonant thermally activated delayed fluorescent (MR-TADF) materials, it is important that there exist computational methods to accurately model their excited states. Here, building on our previous work, we demonstrate how the spin-component scaling second-order approximate coupled-cluster (SCS-CC2), a wavefunction-based method, is robust at predicting the ΔE_ST (i.e., the energy difference between the lowest singlet S₁ and triplet T₁ excited states) of a large number of MR-TADF materials, with a mean average deviation (MAD) of 0.04 eV compared to experimental data. Time-dependent density functional theory calculations with the most common DFT functionals as well as the consideration of the Tamm-Dancoff approximation (TDA) consistently predict a much larger ΔE_ST as a result of a poorer account of Coulomb correlation as compared to SCS-CC2. Very interestingly, the use of a metric to assess the importance of higher order excitations in the SCS-CC2 wavefunctions shows that Coulomb correlation effects are substantially larger in the lowest singlet compared to the corresponding triplet and need to be accounted for a balanced description of the relevant electronic excited states. This is further highlighted with coupled cluster singles-only calculations, which predict very different S₁ energies as compared to SCS-CC2 while T₁ energies remain similar, leading to very large ΔE_ST, in complete disagreement with the experiments. We compared our SCS-CC2/cc-pVDZ with other wavefunction approaches, namely, CC2/cc-pVDZ and SOS-CC2/cc-pVDZ leading to similar performances. Using SCS-CC2, we investigate the excited-state properties of MR-TADF emitters showcasing large ΔE_T2T1 for the majority of emitters, while π-electron extension emerges as the best strategy to minimize ΔE_ST. We also employed SCS-CC2 to evaluate donor-acceptor systems that contain a MR-TADF moiety acting as the acceptor and show that the broad emission observed for some of these compounds arises from the solvent-promoted stabilization of a higher-lying charge-transfer singlet state (S₂). This work highlights the importance of using wavefunction methods in relation to MR-TADF emitter design and associated photophysics.

Benchmarking of Density Functionals for the Description of Optical Properties of Newly Synthesized π-Conjugated TADF Blue Emitters

Computational modeling of the optical characteristics of organic molecules with potential for thermally activated delayed fluorescence (TADF) may assist markedly the development of more efficient emitting materials for organic light-emitting diodes. Recent theoretical studies in this area employ mostly methods from density functional theory (DFT). In order to obtain accurate predictions within this approach, the choice of a proper functional is crucial. In the current study, we focus on testing the performance of a set of DFT functionals for estimation of the excitation and emission energy and the excited singlet-triplet energy gap of three newly synthesized compounds with capacity for TADF. The emitters are designed specifically to enable charge transfer by π-electron conjugation, at the same time possessing high-energy excited triplet states. The functionals chosen for testing are from various groups ranging from gradient-corrected through global hybrids to range-separated ones. The results show that the monitored optical properties are especially sensitive to how the long-range part of the exchange energy is treated within the functional. The accurate functional should also be able to provide well balanced distribution of the π-electrons among the molecular fragments. Global hybrids with moderate (less than 0.4) share of exact exchange (B3LYP, PBE0) and the meta-GGA HSE06 are outlined as the best performing methods for the systems under study. They can predict all important optical parameters correctly, both qualitatively and quantitatively.

Photophysical Properties of Benzophenone-Based TADF Emitters in Relation to Their Molecular Structure

Thermally activated delayed fluorescence (TADF) materials are commonly used in various apparatus, including organic light-emitting device-based displays, as they remarkably improve the internal quantum efficiencies. Although there is a wide range of donor-acceptor-based compounds possessing TADF properties, in this computational study, we investigated TADF and some non-TADF chromophores, containing benzophenone or its structural derivatives as the acceptor core, together with various donor moieties. Following the computational modeling of the emitters, several excited state properties, such as the absorption spectra, singlet-triplet energy gaps (ΔE_ST), natural transition orbitals, and the topological Φ_S indices, have been computed. Along with the donor-acceptor torsion angles and spin-orbit coupling values, these descriptors have been utilized to investigate potential TADF efficiency. Our study has shown that on the one hand, our photophysical/structural descriptors and computational methodologies predict the experimental results quite well, and on the other hand, our extensive benchmark can be useful to pinpoint the most promising functionals and descriptors for the study of benzophenone-based TADF emitters.

Solvent effects on excited-state relaxation dynamics of paddle-wheel BODIPY-Hexaoxatriphenylene conjugates: Insights from non-adiabatic dynamics simulations

Understanding the excited state dynamics of donor-acceptor (D-A) complexes is of fundamental importance both experimentally and theoretically. Herein, we have first explored the photoinduced dynamics of a recently synthesized paddle-wheel BODIPY-hexaoxatriphenylene (BODIPY is the abbreviation for BF2-chelated dipyrromethenes) conjugates D-A complexes with the combination of both electronic structure calculations and non-adiabatic dynamics simulations. On the basis of computational results, we concluded that the BODIPY-hexaoxatriphenylene (BH) conjugates will be promoted to the local excited (LE) states of the BODIPY fragments upon excitation, which is followed by the ultrafast exciton transfer from LE state to charge transfer (CT). Instead of the photoinduced electron transfer process proposed in previous experimental work, such a exciton transfer process is accompanied with the photoinduced hole transfer from BODIPY to hexaoxatriphenylene. Additionally, solvent effects are found to play an important role in the photoinduced dynamics. Specifically, the hole transfer dynamics is accelerated by the acetonitrile solvent, which can be ascribed to significant influences of the solvents on the charge transfer states, i.e. the energy gaps between LE and CT excitons are reduced greatly and the non-adiabatic couplings are increased in the meantime. Our present work not only provides valuable insights into the underlying photoinduced mechanism of BH, but also can be helpful for the future design of novel donor-acceptor conjugates with better optoelectronic performance.

Spin Orbit Coupling in Orthogonal Charge Transfer States: (TD-)DFT of Pyrene-Dimethylaniline

The conformational dependence of the matrix element for spin-orbit coupling and of the electronic coupling for charge separation are determined for an electron donor-acceptor system containing a pyrene acceptor and a dimethylaniline donor. Different kinetic and energetic aspects that play a role in the spin-orbit charge transfer intersystem crossing (SOCT-ISC) mechanism are discussed. This includes parameters related to initial charge separation and the charge recombination pathways using the Classical Marcus Theory of electron transfer. The spin-orbit coupling, which plays a significant role in charge recombination to the triplet state, can be probed by (TD)-DFT, using the latter as a tool to understand and predict the SOCT-ISC mechanism. The matrix elements for spin-orbit coupling for acetone and 4-thio-thymine are used for benchmarking. (Time Dependent-) Density Functional Theory (DFT and TD-DFT) calculations are applied using the quantum chemical program Amsterdam Density Functional (ADF).

Extremely large static and dynamic nonlinear optical response of small superalkali clusters NM3M' (M, M'=Li, Na, K)

Exploring novel nonlinear optical (NLO) materials with excess electron properties is essential for advancing the use of excess electron compounds in optics. The studied superalkali clusters NM₃M' (M, M' = Li, Na, K) are thermodynamically stable and their binding energies range from -27.10 to -53.84 kcal mol^-1. The observed significant values for VIPs suggest their electronic stabilities. Being excess electron candidate these clusters show significant β_o value (3.9 × 10⁷ au), which nicely correlates the hyperpolarizability reported by a two-level model (β_tl). Furthermore, these clusters exhibit a remarkable static second hyperpolarizability (γ_o) value of 1.1 × 10¹⁰ au for the NK₄ superalkali cluster. The hyper Rayleigh scattering (β_HRS) is also computed where the highest value of 2.9 × 10⁷ is recorded for NNa₃K superalkali. The obtained values of β_vec values (projection of hyperpolarizability on dipole moment vector) also signify the excellent nonlinearity of clusters. Besides, the calculated electro-optica pockel's effect β(-ω; ω,0) and second harmonic generation β(-2ω; ω, ω) values are much pronounced at larger dispersion frequency ω = 1064 nm. Moreover, the frequency-dependent second hyperpolarizability γ(ω) with dc-Kerr effect γ(-ω; ω,0,0) and electric field induced second harmonic generation γ(-2ω; ω,ω,0) show larger values at ω = 1064 nm. Thus the highest value of the dc-Kerr constant increases up to 1.0 × 10¹¹ au which also signifies the larger nonlinear refractive index of the studied cluster. We hope this work could open up new possibilities using superalkali clusters as NLO materials for optoelectronics, laser, second harmonic generation and as frequency doubler.

Heavy-Atom-Free Bay-Substituted Perylene Diimide Donor-Acceptor Photosensitizers

Perylene diimide (PDI) dyes are extensively investigated because of their favorable photophysical characteristics for a wide range of organic material applications. Fine-tuning of the optoelectronic properties is readily achieved by functionalization of the electron-deficient PDI scaffold. Here, we present four new donor-acceptor type dyads, wherein the electron donor units - benzo[1,2-b : 4,5-b']dithiophene, 9,9-dimethyl-9,10-dihydroacridine, dithieno[3,2-b : 2',3'-d]pyrrole, and triphenylamine-are attached to the bay-positions of the PDI acceptor. Intersystem crossing occurs for these systems upon photoexcitation, without the aid of heavy atoms, resulting in singlet oxygen quantum yields up to 80 % in toluene solution. Furthermore, this feature is retained when the system is directly irradiated with energy corresponding to the intramolecular charge-transfer absorption band (at 639 nm). Geometrical optimization and (time-dependent) density functional theory calculations afford more insights into the requirements for intersystem crossing such as spin-orbit coupling, dihedral angles, the involvement of charge-transfer states, and energy level alignment.