Critical analysis of spectral solvent shifts calculated by the contemporary PCM approaches of a representative series of charge-transfer complexes between tetracyanoethylene and methylated benzenes

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.

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.

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.

Modifications of the chromophore of Spinach aptamer based on QM:MM calculations

Spinach aptamer was developed as an RNA analog of the green fluorescent protein. The aptamer interacts with its ligand and modifies its electronic spectrum so that it fluoresces brightly at the wavelength of 501 nm. Song et al. investigated modifications of the ligand in their experimental study and found a molecule emitting at 523 nm upon creating a complex with the Spinach aptamer. The crystal structure of the aptamer in complex with its original ligand has been published, which enabled us to study the system computationally. In this article, we suggest several new modifications of the ligand that shift the emission maximum of the complex to even longer wavelengths. Our results are based on combined quantum mechanical/molecular mechanical calculations with DFT method used for geometry optimization and TD-DFT for calculations of absorption and emission energies.

On the accuracy of the general, state-specific polarizable-continuum model for the description of correlated ground- and excited states in solution

Equilibrium and non-equilibrium formulations of the state-specific PCM are evaluated in combination with correlated ground- and excited-state densities provided by ADC/ISR approach of up to third order of perturbation theory.