Prediction of Excited-State Energies and Singlet–Triplet Gaps of Charge-Transfer States Using a Restricted Open-Shell Kohn–Sham Approach

Computed Gyration Tensors of Knotted Chiral and Achiral Topological Stereoisomers of C60 Cyclocarbons

The electronic origins of the computed optical rotations of the simplest chiral and achiral chemical knots with comparatively simple compositions and large, anticipated magnetoelectric polarizabilities are provided. Linear response theory (LRT) is used to calculate the gyration at 1064 nm of two knotted polyyne chains, topological stereoisomers of cyclo[60]carbon. One isomer is analogous to the trefoil knot with approximate D3 symmetry and the other to the figure eight knot with approximate S4 symmetry. The response in each case can be attributed largely to the magnetic dipole term that arises in a near degenerate E-like excited state. An oriented achiral figure eight knot is as optically active in some directions as the chiral knot in any direction, and its absolute eigenvalues are larger.

Orbital-Optimized Versus Time-Dependent Density Functional Calculations of Intramolecular Charge Transfer Excited States

The performance of time-independent, orbital-optimized calculations of excited states is assessed with respect to charge transfer excitations in organic molecules in comparison to the linear-response time-dependent density functional theory (TD-DFT) approach. A direct optimization method to converge on saddle points of the electronic energy surface is used to carry out calculations with the local density approximation (LDA) and the generalized gradient approximation (GGA) functionals PBE and BLYP for a set of 27 excitations in 15 molecules. The time-independent approach is fully variational and provides a relaxed excited state electron density from which the extent of charge transfer is quantified. The TD-DFT calculations are generally found to provide larger charge transfer distances compared to the orbital-optimized calculations, even when including orbital relaxation effects with the Z-vector method. While the error on the excitation energy relative to theoretical best estimates is found to increase with the extent of charge transfer up to ca. -2 eV for TD-DFT, no correlation is observed for the orbital-optimized approach. The orbital-optimized calculations with the LDA and the GGA functionals provide a mean absolute error of ∼0.7 eV, outperforming TD-DFT with both local and global hybrid functionals for excitations with a long-range charge transfer character. Orbital-optimized calculations with the global hybrid functional B3LYP and the range-separated hybrid functional CAM-B3LYP on a selection of states with short- and long-range charge transfer indicate that inclusion of exact exchange has a small effect on the charge transfer distance, while it significantly improves the excitation energy, with the best-performing functional CAM-B3LYP providing an absolute error typically around 0.15 eV.

How does theory compare to experiment for oscillator strengths in electronic spectra? Proposing range-separated hybrids with reliable accountability

As an important quantity in atomic and molecular spectroscopy, oscillator strength should be mentioned. Oscillator strength is linked to the transition dipole moment and consequently to the transition probability between two states, where its magnitude is directly connected to the intensity of the peaks in ultraviolet-visible spectra. However, accurately accounting for oscillator strengths still remains one of the greatest challenges in theory and experiment. Given previous efforts in the context of investigations into oscillator strengths, the related theoretical treatments are relatively limited and have proven to be challenging. In this work, the oscillator strengths in the electronic spectra of organic compounds have thoroughly been investigated with the help of optimally tuned range-separated hybrids (OT-RSHs). In particular, variants of the OT-RSHs combined with the polarizable continuum model (PCM), OT-RSHs-PCM, as well as their screened versions accounting for the screening effects by the electron correlation through the dielectric constant, OT-SRSHs-PCM, are proposed for reliable prediction of the oscillator strengths. The role of the involved ingredients in the proposed methods, namely the underlying density functional approximations, short-range and long-range Hartree-Fock (HF) exchange, as well as the range-separation parameter, has been examined in detail. It is shown that any combination of the parameters in the proposed approximations does not render the reliable oscillator strengths, but a particular compromise among them is needed to describe the experimental data well. Perusing all the results of our developed methods, the best ones are found to be the generalized gradient approximation-based OT-RSHs-PCM, coupled with the linear response theory in the non-equilibrium solvation regime, with the correct asymptotic behavior and incorporating no (low) HF exchange contributions in the short-range part. The best proposed approximations also reveal superior performances not only with respect to their standard counterparts with the default parameters but also as compared to earlier range-separated functionals. Finally, the applicability of the best approximation is also put into broader perspective, where it is used for predicting the oscillator strengths in other sets of compounds not included in the process of developing the approximations. Hopefully, our proposed method can function as an affordable alternative to the expensive wave function-based methods for both theoretical modeling and confirming the experimental observations in the field of electronic spectroscopy.

Orbital-optimized density functional calculations of molecular Rydberg excited states with real space grid representation and self-interaction correction

Density functional calculations of Rydberg excited states up to high energy are carried out for several molecules using an approach where the orbitals are variationally optimized by converging on saddle points on the electronic energy surface within a real space grid representation. Remarkably good agreement with experimental estimates of the excitation energy is obtained using the generalized gradient approximation (GGA) functional of Perdew, Burke, and Ernzerhof (PBE) when Perdew-Zunger self-interaction correction is applied in combination with complex-valued orbitals. Even without the correction, the PBE functional gives quite good results despite the fact that corresponding Rydberg virtual orbitals have positive energy in the ground state calculation. Results obtained using the Tao, Perdew, Staroverov, and Scuseria (TPSS) and r2SCAN meta-GGA functionals are also presented, but they do not provide a systematic improvement over the results from the uncorrected PBE functional. The grid representation combined with the projector augmented-wave approach gives a simpler and better representation of diffuse Rydberg orbitals than a linear combination of atomic orbitals with commonly used basis sets, the latter leading to an overestimation of the excitation energy due to confinement of the excited states.

Triplet states in the reaction center of Photosystem II

In oxygenic photosynthesis sunlight is harvested and funneled as excitation energy into the reaction center (RC) of Photosystem II (PSII), the site of primary charge separation that initiates the photosynthetic electron transfer chain. The chlorophyll ChlD1 pigment of the RC is the primary electron donor, forming a charge-separated radical pair with the vicinal pheophytin PheoD1 (ChlD1+PheoD1-). To avert charge recombination, the electron is further transferred to plastoquinone QA, whereas the hole relaxes to a central pair of chlorophylls (PD1PD2), subsequently driving water oxidation. Spin-triplet states can form within the RC when forward electron transfer is inhibited or back reactions are favored. This can lead to formation of singlet dioxygen, with potential deleterious effects. Here we investigate the nature and properties of triplet states within the PSII RC using a multiscale quantum-mechanics/molecular-mechanics (QM/MM) approach. The low-energy spectrum of excited singlet and triplet states, of both local and charge-transfer nature, is compared using range-separated time-dependent density functional theory (TD-DFT). We further compute electron paramagnetic resonance properties (zero-field splitting parameters and hyperfine coupling constants) of relaxed triplet states and compare them with available experimental data. Moreover, the electrostatic modulation of excited state energetics and redox properties of RC pigments by the semiquinone QA- is described. The results provide a detailed electronic-level understanding of triplet states within the PSII RC and form a refined basis for discussing primary and secondary electron transfer, charge recombination pathways, and possible photoprotection mechanisms in PSII.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Exciton Nature of Plasma Phase Transition in Warm Dense Fluid Hydrogen: ROKS Simulation

The transition of warm dense fluid hydrogen from an insulator to a conducting state at pressures of about 20-400 GPa and temperatures of 500-5000 K has been the subject of active scientific research over the past few decades. However, various experimental and theoretical methods do not provide consistent results. In this work, we have applied the restricted open-shell Kohn-Sham (ROKS) method for first principles molecular dynamics of dense hydrogen after thermal excitation to the first singlet excited state. The Wannier localization method has allowed us to analyze the exciton dynamics in this system. The model shows that a key mechanism of the transition is associated with the dissociation of electron-hole pairs, which allows explaining several stages of the transition of fluid H₂ from molecular state to plasma. This mechanism is able to give a quantitative description of several experimental results as well as to resolve the discrepancies between experimental studies.

Triplets with a Twist: Ultrafast Intersystem Crossing in a Series of Electron Acceptor Materials Driven by Conformational Disorder

Control over the populations of singlet and triplet excitons is key to organic semiconductor technologies. In different contexts, triplets can represent an energy loss pathway that must be managed (i.e., solar cells, light-emitting diodes, and lasers) or provide avenues to improve energy conversion (i.e., photon upconversion and multiplication systems). A key consideration in the interplay of singlet and triplet exciton populations in these systems is the rate of intersystem crossing (ISC). In this work, we design, measure, and model a series of new electron acceptor molecules and analyze them using a combination of ultrafast transient absorption and ultrafast broadband photoluminescence spectroscopies. We demonstrate that intramolecular triplet formation occurs within several hundred picoseconds in solution and is accelerated considerably in the solid state. Importantly, ISC occurs with sufficient rapidity to compete with charge formation in modern organic solar cells, implicating triplets in intrinsic exciton loss channels in addition to charge recombination. Density functional theory calculations reveal that ISC occurs in triplet excited states characterized by local deviations from orbital π-symmetry associated with rotationally flexible thiophene rings. In disordered films, structural distortions, therefore, result in significant increases in spin-orbit coupling, enabling rapid ISC. We demonstrate the generality of this proposal in an oligothiophene model system where ISC is symmetry-forbidden and show that conformational disorder introduced by the formation of a solvent glass accelerates ISC, outweighing the lower temperature and increased viscosity. This proposal sheds light on the factors responsible for facile ISC and provides a simple framework for molecular control over spin states.

Either Accurate Singlet-Triplet Gaps or Excited-State Structures: Testing and Understanding the Performance of TD-DFT for TADF Emitters

The energy gap between the lowest singlet and triplet excited states (ΔE_ST) is a key property of thermally activated delayed fluorescence (TADF) emitters, where these states are dominated by charge-transfer (CT) character. Despite its well-known shortcomings concerning CT states, time-dependent density functional theory (TD-DFT) is widely used to predict this gap and study TADF. Moreover, polar CT states exhibit a strong interaction with their molecular environment, which further complicates their computational description. Addressing these two major challenges, this work studies the performance of Tamm-Dancoff-approximated TD-DFT (TDA-DFT) on the recent STGABS27 benchmark set,¹ exploring different strategies to include orbital and structural relaxation, as well as dielectric embedding. The results show that the best-performing strategy is to calculate ΔE_ST at the ground-state structure using functionals with a surprisingly small amount of Fock exchange of ≈10% and without a (complete) solvent model. However, as this approach heavily relies on error cancellation to mimic dielectric relaxation, it is not robust and exhibits large systematic deviations in excited state energies, state characters, and structures. More rigorous approaches, including state-specific solvation, do not share these systematic deviations, but their predicted ΔE_ST values exhibit larger statistical errors. We thus conclude that for the description of CT states in dielectric environments, none of the tested TDA-DFT methods is competitive with the recently presented ROKS/PCM approach regarding robustness, accuracy, and computational efficiency.

Toward highly efficient hyperfluorescence-based emitters through excited-states alignment using novel optimally tuned range-separated models

Hyperfluorescence has recently been introduced as a promising strategy to achieve organic light-emitting diodes (OLEDs) with high color purity and enhanced stability. In this approach, fluorescent emitters (FEs) with strong and narrow band fluorescence are integrated in thin films containing sensitizers exhibiting thermally activated delayed fluorescence (TADF). Toward highly efficient hyperfluorescence-based emitters, the excited-states ordering of the FEs should be well-aligned. Given some recent endeavors in this context, the related theoretical explorations are relatively limited and have proven to be challenging. In this work, alignments of the corresponding excited-states, crucial for both the fast Förster resonance energy transfer and suppression of the Dexter energy transfer from TADF sensitizers to FEs, have theoretically been investigated using optimally tuned range-separated hybrid functionals (OT-RSHs). We have proposed and validated several variants of the models including OT-RSHs, their coupled versions with the polarizable continuum model, OT-RSHs-PCM, as well as the screened versions accounting for the screening effects by the electron correlation through the scalar dielectric constant, OT-SRSHs, for a reliable description of the excited-states ordering in the FEs of the hyperfluorescence-based materials. Particular attention is paid to the influence of the underlying density functional approximations as well as the short- and long-range Hartree-Fock (HF) exchange contributions and the range-separation parameter. Considering a series of experimentally known hyperfluorescence-based emitters as working models, it is unveiled that any combination of the ingredients in the proposed models does not render the correct order of the excited-states of the FEs, but a particular compromise among the involved parameters is needed to more accurately account for the relevant excited-states alignment. Perusing the results of our developed methods, the best ones are found to be the generalized gradient approximation-based OT-RSHs-PCM with the correct asymptotic behavior and incorporating no (low) HF exchange contribution at the short-range regime. The proposed models show superior performances not only with respect to their standard counterparts with the default parameters but also as compared to other range-separated approximations. Accountability of the best-proposed model is also put into broader perspective, where it has been employed for the computational design of several molecules as promising FE candidates prone to be utilized in hyperfluorescence-based materials. Summing up, the proposed models in this study can be recommended for both the theoretical modeling and confirming the experimental observations in the field of hyperfluorescence-based OLEDs.

Recent Advances in Cartesian-Grid DFT in Atoms and Molecules

In the past several decades, density functional theory (DFT) has evolved as a leading player across a dazzling variety of fields, from organic chemistry to condensed matter physics. The simple conceptual framework and computational elegance are the underlying driver for this. This article reviews some of the recent developments that have taken place in our laboratory in the past 5 years. Efforts are made to validate a viable alternative for DFT calculations for small to medium systems through a Cartesian coordinate grid- (CCG-) based pseudopotential Kohn–Sham (KS) DFT framework using LCAO-MO ansatz. In order to legitimize its suitability and efficacy, at first, electric response properties, such as dipole moment (μ), static dipole polarizability (α), and first hyperpolarizability (β), are calculated. Next, we present a purely numerical approach in CCG for proficient computation of exact exchange density contribution in certain types of orbital-dependent density functionals. A Fourier convolution theorem combined with a range-separated Coulomb interaction kernel is invoked. This takes motivation from a semi-numerical algorithm, where the rate-deciding factor is the evaluation of electrostatic potential. Its success further leads to a systematic self-consistent approach from first principles, which is desirable in the development of optimally tuned range-separated hybrid and hyper functionals. Next, we discuss a simple, alternative time-independent DFT procedure, for computation of single-particle excitation energies, by means of “adiabatic connection theorem” and virial theorem. Optical gaps in organic chromophores, dyes, linear/non-linear PAHs, and charge transfer complexes are faithfully reproduced. In short, CCG-DFT is shown to be a successful route for various practical applications in electronic systems.

Transition-Based Constrained DFT for the Robust and Reliable Treatment of Excitations in Supramolecular Systems

Despite the variety of available computational approaches, state-of-the-art methods for calculating excitation energies, such as time-dependent density functional theory (TDDFT), are computationally demanding and thus limited to moderate system sizes. Here, we introduce a new variation of constrained DFT (CDFT), wherein the constraint corresponds to a particular transition (T), or a combination of transitions, between occupied and virtual orbitals, rather than a region of the simulation space as in traditional CDFT. We compare T-CDFT with TDDFT and ΔSCF results for the low-lying excited states (S1 and T1) of a set of gas-phase acene molecules and OLED emitters and with reference results from the literature. At the PBE level of theory, T-CDFT outperforms ΔSCF for both classes of molecules, while also proving to be more robust. For the local excitations seen in the acenes, T-CDFT and TDDFT perform equally well. For the charge transfer (CT)-like excitations seen in the OLED molecules, T-CDFT also performs well, in contrast to the severe energy underestimation seen with TDDFT. In other words, T-CDFT is equally applicable to both local excitations and CT states, providing more reliable excitation energies at a much lower computational cost than TDDFT cost. T-CDFT is designed for large systems and has been implemented in the linear-scaling BigDFT code. It is therefore ideally suited for exploring the effects of explicit environments on excitation energies, paving the way for future simulations of excited states in complex realistic morphologies, such as those which occur in OLED materials.

Revisiting the Performance of Time-Dependent Density Functional Theory for Electronic Excitations: Assessment of 43 Popular and Recently Developed Functionals from Rungs One to Four

In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between computational efficiency and accuracy. The functionals ωB97X-D and BMK are found to offer the best performance overall with a root-mean square error (RMSE) of around 0.27 eV, better than the computationally more demanding CIS(D) wave function method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for time-dependent density functional theory excitation energies, though hybrid meta generalized-gradient approximations (meta-GGAs) are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.

The ΔSCF method for non-adiabatic dynamics of systems in the liquid phase

Computational studies of ultrafast photoinduced processes give valuable insights into the photochemical mechanisms of a broad range of compounds. In order to accurately reproduce, interpret, and predict experimental results, which are typically obtained in a condensed phase, it is indispensable to include the condensed phase environment in the computational model. However, most studies are still performed in vacuum due to the high computational cost of state-of-the-art non-adiabatic molecular dynamics (NAMD) simulations. The quantum mechanical/molecular mechanical (QM/MM) solvation method has been a popular model to perform photodynamics in the liquid phase. Nevertheless, the currently used QM/MM embedding techniques cannot sufficiently capture all solute-solvent interactions. In this Perspective, we will discuss the efficient ΔSCF electronic structure method and its applications with respect to the NAMD of solvated compounds, with a particular focus on explicit quantum mechanical solvation. As more research is required for this method to reach its full potential, some challenges and possible directions for future research are presented as well.

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.

Charge-Transfer Excitation within a Hybrid-(G)KS Framework through Cartesian Grid DFT

Organic molecules that exhibit charge-transfer (CT) excited states are known to play an important role in processes linked to electron transfer properties and molecular conductance. In this article, we present a simple technique based on "Becke's excitation theorem" that offers an accurate picture of these electronic states. It expresses the correlated energy splitting between triplet and its corresponding singlet states by a two-electron integral, which is numerically evaluated by our recently developed strategy on Cartesian grid. We first examine the consistency of our adopted numerical strategy to evaluate the integral with the original prescribed technique. Then we assess the method on weakly bound CT complexes with three different functionals (BLYP, B3LYP, and LC-BLYP). The accuracy on asymptotic limit of CT excitation is also explored. Finally in order to illustrate the strength and feasibility, it is further extended to a few "challenging" molecules. The method, when employed with hybrid B3LYP functional, turns out to be quite accurate to describe CT excitation energy.

Electronic excitations through the prism of mean-field decomposition techniques

The potential of mean-field decomposition techniques in interpreting electronic transitions in molecules is explored, in particular, the usefulness of these for offering computational signatures of different classes of such excitations. When viewed as a conceptual lens for this purpose, decomposed results are presented for ground- and excited-state energies and dipole moments of selected prototypical organic dyes, and the discrete nature of these properties as well as how they change upon transitioning from one state to another is analyzed without recourse to a discussion based on the involved molecular orbitals. On the basis of results obtained both with and without an account of continuum solvation, our work is further intended to shed new light on practical and pathological differences in between various functional approximations in orbital-optimized Kohn-Sham density functional theory for excited states, equipping practitioners and developers in the field with new probes and possible validation tools.

A two-step MM and QM/MM approach to model AIEE of aryloxy benzothiadiazole derivatives for optoelectronic applications

Aryloxy-benzothiadiazole (ArO-Btz) derivatives show aggregation-induced enhanced emission (AIEE) in the solid-state and are promising candidates for optoelectronic applications. However, understanding the AIEE is a challenging task and is necessary for the rational molecular design of emitters. Therefore, in the present study, electron acceptors (-F, -CN, -NO₂, and -COOH) on the benzothiadiazole ring have been screened for emission in solution and aggregated phases. Herein, we report QM (DFT/TDDFT) and ONIOM (QM/MM) studies on the four ArO-Btz derivatives in comparison with the parent molecule with typical characteristics of AIEE, optoelectronic and non-linear optical properties. Starting from the optimized crystal structure of the parent compound, the structures of the designed clusters have been pre-optimized with MM and then with QM/MM to explore their absorption and emission in the solid phase. The results indicate that in the aggregated phase, the surrounding environment reduces intra-molecular rotations and molecular motion that lead to enhanced emission. Natural bond orbital (NBO) analyses reveal that the ground state structure is stabilized from electron delocalization and operative push-pull effects. Interestingly, nitro-benzothiadiazole exhibits prominent AIEE phenomena, with an emission wavelength beyond 700 nm in solution and in the cluster, reinforced by the magnification of its oscillatory strength by 100 times when aggregated. This dinitro-aryloxy-benzothiadiazole derivative is proposed as a near-infrared emitter for dye-sensitized solar cell, optoelectronic, and non-linear optical applications.

Next generation quantum theory of atoms in molecules for the design of emitters exhibiting thermally activated delayed fluorescence with laser irradiation

The effect of a static electric (E)-field and an unchirped and chirped laser pulse field on the cycl[3.3.3]azine molecule was investigated using next-generation quantum theory of atoms in molecules (NG-QTAIM). Despite the magnitude of the E-field of the laser pulses being an order of magnitude lower than for the static E-field, the variation of the energy gap between the lowest lying singlet (S₁ ) and triplet (T₁ ) excited states was orders of magnitude greater for the laser pulse than for the static E-field. Insights into the response of the electronic structure were captured by NG-QTAIM, where differences in the inverted singlet-triplet gap due to the laser pulses were significant larger compared to those induced by the static E-field. The response of the S₁ and T₁ excited states, as determined by NG-QTAIM, switched discontinuously between weak and strong chemical character for the static E-field. In contrast, the response to the laser pulses, determined by NG-QTAIM, is to induce a continuous range of chemical character, indicating the unique ability of the laser pulses to induce polarization effects in the form of "mixed" bond types. Our analysis demonstrates that NG-QTAIM is a useful tool for understanding the response to laser irradiation of the lowest-lying singlet S₁ and triplet T₁ excited states of emitters exhibiting thermally activated delayed fluorescence. The chirped laser pulse led to more frequent instances of the desired outcome of an inverted singlet-triplet gap than the unchirped pulse, indicating its usefulness as a tool to design more efficient organic light-emitting diode devices.

Simulation of Low-Lying Singlet and Triplet Excited States of Multiple-Resonance-Type Thermally Activated Delayed Fluorescence Emitters by Delta Self-Consistent Field (ΔSCF) Method

The delta self-consistent field (ΔSCF) method was applied to the simulation of low-lying singlet and triplet excited states of multiple-resonance (MR)-type thermally activated delayed fluorescence (TADF) molecules, which form a promising group for organic light-emitting diode (OLED) emitters. A comparison with the experimental values of 13 emitters from the literature showed that ΔSCF gave fairly accurate S₁ and T₁ excitation energies (mean absolute errors (MAEs) of 0.092 and 0.055 eV, respectively) as well as quite accurate ΔE_ST (S₁-T₁ gap, MAE of 0.041 eV), which could not be calculated with sufficient accuracy by the conventional time-dependent density functional theory (TDDFT). ΔSCF also demonstrated its utility for the analysis of photophysical properties through a simulation of the reverse intersystem crossing (RISC) process of an MR-type emitter.

Self-Improving Photosensitizer Discovery System via Bayesian Search with First-Principle Simulations

Artificial intelligence (AI) based self-learning or self-improving material discovery system will enable next-generation material discovery. Herein, we demonstrate how to combine accurate prediction of material performance via first-principle calculations and Bayesian optimization-based active learning to realize a self-improving discovery system for high-performance photosensitizers (PSs). Through self-improving cycles, such a system can improve the model prediction accuracy (best mean absolute error of 0.090 eV for singlet-triplet spitting) and high-performance PS search ability, realizing efficient discovery of PSs. From a molecular space with more than 7 million molecules, 5357 potential high-performance PSs were discovered. Four PSs were further synthesized to show performance comparable with or superior to commercial ones. This work highlights the potential of active learning in first-principle-based materials design, and the discovered structures could boost the development of photosensitization related applications.

PCM-ROKS for the Description of Charge-Transfer States in Solution: Singlet-Triplet Gaps with Chemical Accuracy from Open-Shell Kohn-Sham Reaction-Field Calculations

The adiabatic energy gap between the lowest singlet and triplet excited states ΔE_ST is a central property of thermally activated delayed fluorescence (TADF) emitters. Since these states are dominated by a charge-transfer character, causing strong orbital-relaxation and environmental effects, an accurate prediction of ΔE_ST is very challenging, even with modern quantum-chemical excited-state methods. Addressing this major challenge, we present an approach that combines spin-unrestricted (UKS) and restricted open-shell Kohn-Sham (ROKS) self-consistent field calculations with a polarizable-continuum model and range-separated hybrid functionals. Tests on a new representative benchmark set of 27 TADF emitters with accurately known ΔE_ST values termed STGABS27 reveal a robust and unprecedented performance with a mean absolute deviation of only 0.025 eV (∼0.5 kcal/mol) and few deviations greater than 0.05 eV (∼1 kcal/mol), even in electronically challenging cases. Requiring only two geometry optimizations per molecule at the ROKS/UKS level in a compact double-ζ basis, the approach is computationally efficient and can routinely be applied to molecules with more than 100 atoms.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Orbital Optimized Density Functional Theory for Electronic Excited States

Density functional theory (DFT) based modeling of electronic excited states is of importance for investigation of the photophysical/photochemical properties and spectroscopic characterization of large systems. The widely used linear response time-dependent DFT (TDDFT) approach is, however, not effective at modeling many types of excited states, including (but not limited to) charge-transfer states, doubly excited states, and core-level excitations. In this perspective, we discuss state-specific orbital optimized (OO) DFT approaches as an alterative to TDDFT for electronic excited states. We motivate the use of OO-DFT methods and discuss reasons behind their relatively restricted historical usage (vs TDDFT). We subsequently highlight modern developments that address these factors and allow efficient and reliable OO-DFT computations. Several successful applications of OO-DFT for challenging electronic excitations are also presented, indicating their practical efficacy. OO-DFT approaches are thus increasingly becoming a useful route for computing excited states of large chemical systems. We conclude by discussing the limitations and challenges still facing OO-DFT methods, as well as some potential avenues for addressing them.

Improved Description of Charge-Transfer Potential Energy Surfaces via Spin-Component-Scaled CC2 and ADC(2) Methods

The molecular level understanding of electronic transport properties depends on the reliable theoretical description of charge-transfer (CT)-type electronic states. In this paper, the performance of spin-component-scaled variants of the popular CC2 and ADC(2) methods is evaluated for CT states, following benchmark strategies of earlier studies that revealed a compromised accuracy of the unmodified models. In addition to statistics on the accuracy of vertical excitation energies at equilibrium and infinite separation of bimolecular complexes, potential energy surfaces of the ammonia-fluorine complex are also reported. The results show the capability of spin-component-scaled approaches to reduce the large errors of their regular counterparts to a significant extent, outperforming even the coupled-cluster single and double method in many cases. The cost-effective scaled-opposite-spin variants are found to provide a remarkably good agreement with the CCSDT-3 reference data, thereby being recommended methods of choice in the study of charge-transfer states.

Elucidating the Magnitude of Internal Reorganization Energy of Molecular Excited States from the Perspective of Transition Density

Quantifying vibronic couplings in molecular excited states is crucial for the elucidation of a broad range of photophysical phenomena. In this study, we compare different theoretical approaches for the calculation of reorganization energy, a measure of vibronic coupling strength, and provide a rigorous derivation to show that molecular transition density characterizing electron-hole excitation could be used to quantify the magnitude of reorganization energy. The theory enables a descriptor based on molecular-orbital coefficients and atomic transition densities to quantify the magnitude of reorganization energies in molecular excited states. Applying the approach to low-lying excited states of polyacenes, we demonstrate that transition density distribution explains the difference in the magnitude of the reorganization energy of different excited states. Furthermore, to clarify the applicability of the transition density descriptor in molecular design for small-reorganization energy molecules, we investigate a broad range of molecular chromophores to show the effectiveness of the proposed theory. With this perspective on the relationship between reorganization energy and transition density, we successfully provide a quantitative rule to identify π-conjugated systems with small reorganization energy in the excited state, which should be useful for the development of novel optoelectronic materials.

Assessment of the Density-Fitted Second-Order Quasidegenerate Perturbation Theory for Transition Energies: Accurate Computations of Singlet-Triplet Gaps for Charge-Transfer Compounds

Organic light-emitting diodes (OLEDs) have been of significant interest because of their superior performance and low cost of production. Thermally activated delayed fluorescence (TADF) has attracted significant interest in the OLED technology because it improves the efficiency of OLEDs by harvesting triplet excitons. Therefore, the accurate computation of singlet-triplet transition energies (ΔE_S1-T1) of charge-transfer molecules is very important. However, the accurate computation of the ΔE_S1-T1 values is a challenging problem for single-reference methods because of the multireference character of excited states. In this research, an assessment of density-fitted second-order quasidegenerate perturbation theory (DF-QDPT2) [Bozkaya, U.; J. Chem. Theory Comput. 2019, 15, 4415-4429] for singlet-triplet transition energies (ΔE_S1-T1) of charge-transfer compounds is presented. The performance of the DF-QDPT2 method has been compared with those of several density-functional theory functionals, such as B3LYP, PBE0, M06-2X, ωB97X-D, and MN15; density-fitted state-averaged CASSCF (DF-SA-CASSCF); and single-state single-reference second-order perturbation theory (SS-SR-CASPT2) methods. For the TADF molecules considered, the DF-QDPT2 method provides a mean absolute error (MAE) of 0.13 eV, while the MAE values of DF-SA-CASSCF and SS-SR-CASPT2 are 0.65 and 0.74 eV, respectively. The performances of B3LYP and PBE0 are slightly better than that of DF-QDPT2, while M06-2X and ωB97X-D provide noticeably higher errors compared with DF-QDPT2. Furthermore, the standard CASSCF without state-averaging yields dramatic errors with an MAE value of 3.0 eV. Our results demonstrate that eigenvalues of the DF-QDPT2-effective Hamiltonian can be reliably used for the prediction of singlet-triplet transition energies, while eigenvalues of DF-CASSCF/DF-SA-CASSCF fail to provide accurate predictions. Overall, we conclude that the DF-QDPT2 method emerges as a very useful tool for the computation of excited-state properties.

A New Benchmark Set for Excitation Energy of Charge Transfer States: Systematic Investigation of Coupled Cluster Type Methods

The numerous existing publications on benchmarking quantum chemistry methods for excited states rarely include Charge Transfer (CT) states, although many interesting phenomena in, e.g., biochemistry and material physics involve the transfer of electrons between fragments of the system. Therefore, it is timely to test the accuracy of quantum chemical methods for CT states, as well. In this study we first propose a new benchmark set consisting of dimers having low-energy CT states. On this set, the vertical excitation energy has been calculated with Coupled Cluster methods including triple excitations (CC3, CCSDT-3, CCSD(T)(a)*), as well as with methods including full or approximate doubles (CCSD, STEOM-CCSD, CC2, ADC(2), EOM-CCSD(2)). The results show that the popular CC2 and ADC(2) methods are much less accurate for CT states than for valence states. On the other hand, EOM-CCSD seems to have similar systematic overestimation of the excitation energies for both types of states. Among the triples methods the novel EOM-CCSD(T)(a)* method including noniterative triple excitations is found to stand out with its consistently good performance for all types of states, delivering essentially EOM-CCSDT quality results.

Constructing Spin-Adiabatic States for the Modeling of Spin-Crossing Reactions. I. A Shared-Orbital Implementation

In the modeling of spin-crossing reactions, it has become popular to directly explore the spin-adiabatic surfaces. Specifically, through constructing spin-adiabatic states from a two-state Hamiltonian (with spin-orbit coupling matrix elements) at each geometry, one can readily employ advanced geometry optimization algorithms to acquire a "transition state" structure, where the spin crossing occurs. In this work, we report the implementation of a fully-variational spin-adiabatic approach based on Kohn-Sham density functional theory spin states (sharing the same set of molecular orbitals) and the Breit-Pauli one-electron spin-orbit operator. For three model spin-crossing reactions [predissociation of N₂O, singlet-triplet conversion in CH₂, and CO addition to Fe(CO)₄], the spin-crossing points were obtained. Our results also indicated the Breit-Pauli one-electron spin-orbit coupling can vary significantly along the reaction pathway on the spin-adiabatic energy surface. On the other hand, due to the restriction that low-spin and high-spin states share the same set of molecular orbitals, the acquired spin-adiabatic energy surface shows a cusp (i.e. a first-order discontinuity) at the crossing point, which prevents the use of standard geometry optimization algorithms to pinpoint the crossing point. An extension with this restriction removed is being developed to achieve the smoothness of spin-adiabatic surfaces.

Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.

Highly Accurate Prediction of Core Spectra of Molecules at Density Functional Theory Cost: Attaining Sub-electronvolt Error from a Restricted Open-Shell Kohn-Sham Approach

We present the use of the recently developed square gradient minimization (SGM) algorithm for excited-state orbital optimization to obtain spin-pure restricted open-shell Kohn-Sham (ROKS) energies for core excited states of molecules. The SGM algorithm is robust against variational collapse and offers a reliable route to converging orbitals for target excited states at only 2-3 times the cost of ground-state orbital optimization (per iteration). ROKS/SGM with the modern SCAN/ωB97X-V functionals is found to predict the K-edge of C, N, O, and F to a root mean squared error of ∼0.3 eV. ROKS/SGM is equally effective at predicting L-edge spectra of third period elements, provided a perturbative spin-orbit correction is employed. This high accuracy can be contrasted with traditional time-dependent density functional theory (TDDFT), which typically has greater than 10 eV error and requires translation of computed spectra to align with experiment. ROKS is computationally affordable (having the same scaling as ground-state DFT and a slightly larger prefactor) and can be applied to geometry optimizations/ab initio molecular dynamics of core excited states, as well as condensed phase simulations. ROKS can also model doubly excited/ionized states with one broken electron pair, which are beyond the ability of linear response based methods.

Beyond the Coulson–Fischer point: characterizing single excitation CI and TDDFT for excited states in single bond dissociations

HF/DFT orbitals spin-polarize when single bonds are stretched past the Coulson–Fischer point. We report unphysical features in the excited state potential energy surfaces predicted by CIS/TDDFT in this regime, and characterize their origin.

Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number

Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps, and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C₁₂H₁₄ polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge-based data point could therefore permit fitting/tuning of functionals with lower delocalization error.

Unveiling the Photophysical Properties of Boron Heptaaryldipyrromethene Derivatives

Increased interest has been devoted to the discovery of multifunctional materials with desirable properties, as continuous performance enhancement of various devices mainly depends on high-performance materials. Now, density functional theory has become a powerful tool to design new materials and rationalize experimental observations. In this work, we explored the photophysical properties origin of chiral boron heptaaryldipyrromethene (heptaaryl-BODIPY), which has charming optoelectronic properties. At the same time, we designed the other five compounds on the basis of heptaaryl-BODIPY. The simulated electronic absorption and emission spectra of heptaaryl-BODIPY are in agreement with experimental ones, allowing us to reliably assign its electronic transition property. The designed compound 6 shows remarkably large first hyperpolarizability value up to 82.78×10^-30  esu. For this kind of compounds, their NLO response values associate with not only position but also electronic nature of substituent groups. Moreover, electron reorganization energies of compounds 1-4 are comparable to tris(8-hydroxyquinolinato)aluminium(III) which is a typical electron transport material. Intriguingly, the studied compounds are the excellent fluorescent probe materials from the standpoint of large Stokes shift and high emission efficiency. Our work enables an opportunity for understanding the relationship between microelectronic structure and macroscopic performance of BODIPY derivatives.

Computational Design of Thermally Activated Delayed Fluorescence Materials: The Challenges Ahead

Thermally activated delayed fluorescence (TADF) offers promise for all-organic light-emitting diodes with quantum efficiencies competing with those of transition-metal-based phosphorescent devices. While computational efforts have so far largely focused on gas-phase calculations of singlet and triplet excitation energies, the design of TADF materials requires multiple methodological developments targeting among others a quantitative description of electronic excitation energetics, fully accounting for environmental electrostatics and molecular conformational effects, the accurate assessment of the quantum mechanical interactions that trigger the elementary electronic processes involved in TADF, and a robust picture for the dynamics of these fundamental processes. In this Perspective, we describe some recent progress along those lines and highlight the main challenges ahead for modeling, which we hope will be useful to the whole TADF community.

Accurate Treatment of Charge-Transfer Excitations and Thermally Activated Delayed Fluorescence Using the Particle-Particle Random Phase Approximation

Thermally activated delayed florescence (TADF) is a mechanism that increases the electroluminescence efficiency in organic light-emitting diodes by harnessing both singlet and triplet excitons. TADF is facilitated by a small energy difference between the first singlet (S1) and triplet (T1) excited states (Δ E(ST)), which is minimized by spatial separation of the donor and acceptor moieties. The resultant charge-transfer (CT) excited states are difficult to model using time-dependent density functional theory (TDDFT) because of the delocalization error present in standard density functional approximations to the exchange-correlation energy. In this work we explore the application of the particle-particle random phase approximation (pp-RPA) for the determination of both S1 and T1 excitation energies. We demonstrate that the accuracy of the pp-RPA is functional dependent and that, when combined with the hybrid functional B3LYP, the pp-RPA computed Δ E(ST) have a mean absolute deviation (MAD) of 0.12 eV for the set of examined molecules. A key advantage of the pp-RPA approach is that the S1 and T1 states are characterized as CT states for all of experimentally reported TADF molecules examined here, which allows for an estimate of the singlet-triplet CT excited state energy gap (Δ E(ST) = 1CT - 3CT). For experimentally known TADF molecules with a small (<0.2 eV) Δ E(ST) in this data set, a high accuracy is demonstrated for the prediction of both the S1 (MAD = 0.18 eV) and T1 (MAD = 0.20 eV) excitation energies as well as Δ E(ST) (MAD = 0.05 eV). This result is attributed to the consideration of correct antisymmetry in the particle-particle interaction leading to the use of full exchange kernel in addition to the Coulomb contribution, as well as a consistent treatment of both singlet and triplet excited states. The computational efficiency of this approach is similar to that of TDDFT, and the cost can be reduced significantly by using the active-space approach.

Modeling TADF in organic emitters requires a careful consideration of the environment and going beyond the Franck-Condon approximation

The origin of the thermally-activated delayed fluorescence (TADF) of three organic emitters is investigated, focusing on the nature of the lowest excited states, their transition properties, as well as the role of the environment. For this purpose, the algebraic-diagrammatic construction for the polarization propagator at second order of perturbation theory [ADC(2)], time-dependent density-functional theory in the Tamm-Dancoff approximation (TDA) and unrestricted Kohn-Sham DFT in combination with the maximum-overlap method (ΔDFT) are employed. The influence of the dielectric environment is rigorously included using different variants of the polarizable continuum model. The calculations reveal the lowest excited singlet and triplet states of all studied emitters to be dominated by charge-transfer (CT) character already in the most apolar environment corresponding to cyclo-hexane. The dielectric stabilization entails a drastic reduction of the singlet-triplet gaps, increasing the calculated TADF rates by several orders of magnitude. Another ingredient for accurate TADF rates is hidden in the excited-state potential-energy surfaces along the donor-acceptor twisting angle. A presence of large, shallow plateaus in apolar environments causes the transition properties to be governed by thermal fluctuations rather than the minimum-energy geometries. This leads to a large increase of the oscillator strengths, as well as a breakdown of the Franck-Condon approximation. The last ingredient is a small but significant spin-orbit coupling (SOC) between the singlet and triplet CT states, which is traced back to a delocalization of the excitation hole or excited electron of the triplet CT state. Taking into account all of these effects, a reasonable agreement with experimental TADF and fluorescence lifetimes is obtained. For this, it proves to be sufficient to consider only the lowest lying singlet and triplet excited states.

How accurate are static polarizability predictions from density functional theory? An assessment over 132 species at equilibrium geometry

Static polarizabilities are the first response of the electron density to electric fields, and offer a formally exact measure of the accuracy of excited states. We have developed a benchmark database of polarizabilities and have assessed the performance of 60 popular and recent functionals in predicting them.

Synthesis and Exciton Dynamics of Donor-Orthogonal Acceptor Conjugated Polymers: Reducing the Singlet-Triplet Energy Gap

The presence of energetically low-lying triplet states is a hallmark of organic semiconductors. Even though they present a wealth of interesting photophysical properties, these optically dark states significantly limit optoelectronic device performance. Recent advances in emissive charge-transfer molecules have pioneered routes to reduce the energy gap between triplets and "bright" singlets, allowing thermal population exchange between them and eliminating a significant loss channel in devices. In conjugated polymers, this gap has proved resistant to modification. Here, we introduce a general approach to reduce the singlet-triplet energy gap in fully conjugated polymers, using a donor-orthogonal acceptor motif to spatially separate electron and hole wave functions. This new generation of conjugated polymers allows for a greatly reduced exchange energy, enhancing triplet formation and enabling thermally activated delayed fluorescence. We find that the mechanisms of both processes are driven by excited-state mixing between π-π*and charge-transfer states, affording new insight into reverse intersystem crossing.

Dissecting the accountability of parameterized and parameter-free single-hybrid and double-hybrid functionals for photophysical properties of TADF-based OLEDs

Organic light emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) emitters are an attractive category of materials that have witnessed a booming development in recent years. In the present contribution, we scrutinize the accountability of parameterized and parameter-free single-hybrid (SH) and double-hybrid (DH) functionals through the two formalisms, full time-dependent density functional theory (TD-DFT) and Tamm-Dancoff approximation (TDA), for the estimation of photophysical properties like absorption energy, emission energy, zero-zero transition energy, and singlet-triplet energy splitting of TADF molecules. According to our detailed analyses on the performance of SHs based on TD-DFT and TDA, the TDA-based parameter-free SH functionals, PBE0 and TPSS0, with one-third of exact-like exchange turned out to be the best performers in comparison to other functionals from various rungs to reproduce the experimental data of the benchmarked set. Such affordable SH approximations can thus be employed to predict and design the TADF molecules with low singlet-triplet energy gaps for OLED applications. From another perspective, considering this point that both the nonlocal exchange and correlation are essential for a more reliable description of large charge-transfer excited states, applicability of the functionals incorporating these terms, namely, parameterized and parameter-free DHs, has also been evaluated. Perusing the role of exact-like exchange, perturbative-like correlation, solvent effects, and other related factors, we find that the parameterized functionals B2π-PLYP and B2GP-PLYP and the parameter-free models PBE-CIDH and PBE-QIDH have respectable performance with respect to others. Lastly, besides the recommendation of reliable computational protocols for the purpose, hopefully this study can pave the way toward further developments of other SHs and DHs for theoretical explorations in the field of OLEDs technology.

Benchmark of Bethe-Salpeter for Triplet Excited-States

We have evaluated the accuracy of the Bethe-Salpeter singlet-triplet transition energies as well as singlet-triplet and triplet-triplet splittings for 20 organic molecules, using as reference the CC3 values determined by Thiel and co-workers with both the TZVP and aug-cc-pVTZ atomic basis sets. Our excitation energies are obtained on the basis of GW quasiparticle energy levels that are self-consistently converged with respect to the starting DFT eigenvalues. In its current form, BSE/GW is often unable to provide a balanced description of both singlet and triplet excited-states. While the singlet-singlet and triplet-triplet energy separations are obtained accurately, triplets are located too close in energy from the ground-state, by typically -0.55 eV when using standard functionals to generate the starting eigenstates. Applying the Tamm-Dancoff approximation upshifts the BSE triplet energies and allows reducing this error to ca. -0.40 eV, while using M06-HF eigenstates allows a further increase and hence a reduction of the error for triplet states, but at the cost of larger errors for the singlet excited-states. At this stage, the most accurate TD-DFT estimates therefore remain competitive for computing singlet-triplet transition energies. Indeed, with M06-2X, irrespective of the application or not of the Tamm-Dancoff approximation and of the selected atomic basis set, the deviations obtained with TD-DFT are rather small.

Theoretical tuning of the singlet–triplet energy gap to achieve efficient long-wavelength thermally activated delayed fluorescence emitters: the impact of substituents

Long-wavelength TADF emitters could be achieved with both small ΔE_ST and high k_r through introducing small electron-withdrawing substituents!