Excited State Dipole Moments in Solution: Comparison between State-Specific and Linear-Response TD-DFT Values

Benchmarking Charge-Transfer Excited States in TADF Emitters: ΔDFT Outperforms TD-DFT for Emission Energies

Charge-transfer (CT) excited states are crucial to organic light-emitting diodes (OLEDs), particularly to those based on thermally activated delayed fluorescence (TADF). However, accurately modeling CT states remains challenging, even with modern implementations of (time-dependent) density functional theory [(TD-)DFT], especially in a dielectric environment. To identify shortcomings and improve the methodology, we previously established the STGABS27 benchmark set with highly accurate experimental references for the adiabatic energy gap between the lowest singlet and triplet excited states (ΔEST). Here, we diversify this set to the STGABS27-EMS benchmark by including experimental emission energies (Eem) and use this new set to (re)-evaluate various DFT-based approaches. Surprisingly, these tests demonstrate that a state-specific (un)restricted open-shell Kohn-Sham (U/ROKS) DFT coupled with a polarizable continuum model for perturbative state-specific nonequilibrium solvation (ptSS-PCM) provides exceptional accuracy for predicting Eem over a wide range of density functionals. In contrast, the main workhorse of the field, Tamm-Dancoff-approximated TD-DFT (TDA-DFT) paired with the same ptSS-PCM, is distinctly less accurate and strongly functional-dependent. More importantly, while TDA-DFT requires the choice of two very different density functionals for good performance on either ΔEST or Eem, the time-independent U/ROKS/PCM approaches deliver excellent accuracy for both quantities with a wide variety of functionals.

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.

Benchmarking Density Functional Approximations for Excited-State Properties of Fluorescent Dyes

This study presents an extensive analysis of the predictive power of time-dependent density functional theory in determining the excited-state properties of two groups of important fluorescent dyes, difluoroboranes and hydroxyphenylimidazo[1,2-a]pyridine derivatives. To ensure statistically meaningful results, the data set is comprised of 85 molecules manifesting diverse photophysical properties. The vertical excitation energies and dipole moments (in the electronic ground and excited states) of the aforementioned dyes were determined using the RI-CC2 method (reference) and with 18 density functional approximations (DFA). The set encompasses DFAs with varying amounts of exact exchange energy (EEX): from 0% (e.g., SVWN, BLYP), through a medium (e.g., TPSSh, B3LYP), up to a major contribution of EEX (e.g., BMK, MN15). It also includes range-separated hybrids (CAM-B3LYP, LC-BLYP). Similar error profiles of vertical energy were obtained for both dye groups, although the errors related to hydroxyphenylimidazopiridines are significantly larger. Overall, functionals including 40-55% of EEX (SOGGA11-X, BMK, M06-2X) ensure satisfactory agreement with the reference vertical excitation energies obtained using the RI-CC2 method; however, MN15 significantly outperforms them, providing a mean absolute error of merely 0.04 eV together with a very high correlation coefficient (R2 = 0.98). Within the investigated set of functionals, there is no single functional that would equally accurately determine ground- and excited-state dipole moments of difluoroboranes and hydroxyphenylimidazopiridine derivatives. Depending on the chosen set of dyes, the most accurate μGS predictions were delivered by MN15 incorporating a major EEX contribution (difluoroboranes) and by PBE0 containing a minor EEX fraction (hydroxyphenylimidazopiridines). Reverse trends are observed for μES, i.e., for difluoroboranes the best results were obtained with functionals including a minor fraction of EEX, specifically PBE0, while in the case of hydroxyphenylimidazopiridines, much more accurate predictions were provided by functionals incorporating a major EEX contribution (BMK, MN15).

Simple Protocol for Capturing Both Linear-Response and State-Specific Effects in Excited-State Calculations with Continuum Solvation Models

We present an effective computational protocol (cLR²) to describe both solvatochromism and fluorosolvatochromism. This protocol, which couples the polarizable continuum model to time-dependent density functional theory, simultaneously accounts for both linear-response and state-specific solvation effects. A series of test cases, including solvatochromic and fluorosolvatochromic compounds and excited-state intramolecular proton transfers, are used to highlight that cLR² is especially beneficial for modeling bright excitations possessing a significant charge-transfer character, as well as cases in which an accurate balance between states of various polarities should be restored.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Solvent effect on the photophysical properties of thermally activated delayed fluorescence molecules

As the third-generation organic electroluminescent materials, thermally activated delayed fluorescence (TADF) molecules have become the research focus recently. Significant solvent effect on TADF molecules were found experimentally, while theoretical investigations are quite limited. In this work, the solvent effect on photophysical properties of DCBPy and DTCBPy are investigated with first-principles calculations. The solvent polarity has slight influence on the molecular geometries and orbitals, while it can decrease the energy gap between the first singlet excited state (S1) and first triplet excited state (T1) significantly. Both the oscillator strength and the radiation rates of S1 increase with larger solvent polarity. The large energy gap between S1 and T1 induce negligible intersystem crossing (ISC) and reverse ISC rates between them, which also indicates higher triplet excited states are involved in the up-conversion process. Our results provide valuable information about solvent influence on the light-emitting properties of TADF molecules, which could help one better understand the light-emitting mechanism of them and favor the design of TADF molecules.

Electronic transitions for a fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles: Linear and corrected linear response regimes

The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFμ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.

The role of solvation models on the computed absorption and emission spectra: the case of fireflies oxyluciferin

Absorption and emission energies calculation covering both implicit and explicit solvation models using oxyluciferin as the case of study.

First-principles investigation of the double ESIPT process in a thiophene-based dye

Explanation of the experimental triple emission with theoretical tools requires advanced solvent models.

The electronic transitions of analogs of red wine pyranoanthocyanin pigments

High level theoretical methods have been used in order to understand structural effects on the absorption spectra of pyranoanthocyanins.

The Bethe-Salpeter formalism with polarisable continuum embedding: reconciling linear-response and state-specific features

The Bethe-Salpeter equation (BSE) formalism has been recently shown to be a valuable alternative to time-dependent density functional theory (TD-DFT) with the same computing time scaling with system size. In particular, problematic transitions for TD-DFT such as charge-transfer, Rydberg and cyanine-like excitations were shown to be accurately described with BSE. We demonstrate here that combining the BSE formalism with the polarisable continuum model (PCM) allows us to include simultaneously linear-response and state-specific contributions to solvatochromism. This is confirmed by exploring transitions of various natures (local, charge-transfer, etc.) in a series of solvated molecules (acrolein, indigo, p-nitro-aniline, donor-acceptor complexes, etc.) for which we compare BSE solvatochromic shifts to those obtained by linear-response and state-specific TD-DFT implementations. Such a remarkable and unique feature is particularly valuable for the study of solvent effects on excitations presenting a hybrid localised/charge-transfer character.

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.