Charge-transfer excitations in low-gap systems under the influence of solvation and conformational disorder: Exploring range-separation tuning

Nonempirical Prediction of the Length-Dependent Ionization Potential in Molecular Chains

The ionization potential of molecular chains is well-known to be a tunable nanoscale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nanosized systems prediction of the evolution of the electronic structure between the two limits is more difficult. Recently, optimal tuning of range-separated hybrid functionals has emerged as a highly accurate method for predicting ionization potentials. This was first achieved for molecules using the ionization potential theorem (IPT) and more recently extended to solid-state systems, based on an ansatz that generalizes the IPT to the removal of charge from a localized Wannier function. Here, we study one-dimensional molecular chains of increasing size, from the monomer limit to the infinite polymer limit using this approach. By comparing our results with other localization-based methods and where available with experiment, we demonstrate that Wannier-localization-based optimal tuning is highly accurate in predicting ionization potentials for any chain length, including the nanoscale regime.

Thiocarbonyl-Bridged N-Heterotriangulenes for Energy Efficient Triplet Photosensitization: A Theoretical Perspective

Structurally-rigid metal-free organic molecules are of high demand for various triplet harvesting applications. However, inefficient intersystem crossing (ISC) due to large singlet-triplet gap (Δ E S - T ${\Delta {E}_{S-T}}$ ) and small spin-orbit coupling (SOC) between lowest excited singlet and triplet often limits their efficiency. Excited electronic states, fluorescence and ISC rates in several thiocarbonyl-bridged N-heterotriangulene ( m ${m}$ S-HTG) with systematically increased thione content (m = ${m=}$ 0-3) are investigated implementing polarization consistent time-dependent optimally-tuned range-separated hybrid. All m ${m}$ S-HTGs are dynamically stable and also thermodynamically feasible to synthesize. Relative energies of several low-lying singlets (S n ${{S}_{n}}$ ) and triplets (T n ${{T}_{n}}$ ), and their excitation nature (i. e.,n π * ${n{\pi }^{^{\ast}}}$ orπ π * ${\pi {\pi }^{^{\ast}}}$ ) and SOC are determined for these m ${m}$ S-HTGs in dichloromethane. Low-energy optical peak displays gradual red-shift with increasing thione content due to relatively smaller electronic gap resulted from greater degree of orbital delocalization. Significantly large SOC due to different orbital-symmetry and heavy-atom effect produces remarkably high ISC rates (k I S C ${{k}_{ISC}}$ ~1012 s-1) for enthalpically favouredS 1 n π * → T 2 ${{S}_{1}\left(n{\pi }^{^{\ast}}\right)\to {T}_{2}}$ (π π * ${\pi {\pi }^{^{\ast}}}$ ) channel in these m ${m}$ S-HTGs, which outcompete radiative fluorescence rates (~108 s-1) even directly from higher lying optically brightπ π * ${\pi {\pi }^{^{\ast}}}$ singlets. Importantly, high energy triplet excitons of ~1.7 eV resulting from such significantly large ISC rates from non-fluorescentS 1 n π * ${{S}_{1}\left(n{\pi }^{^{\ast}}\right)}$ make these thiocarbonylated HTGs ideal candidates for energy efficient triplet harvest including triplet-photosensitization.

Screened optimally tuned range separated hybrid functional for solvated low bandgap molecular systems

The description of electronic properties of low bandgap molecular system is often performed by using density functional theory (DFT) and time dependent (TD) DFT calculations with the optimally tuned range-separated hybrid (OT-RSH) functional, as it contains the necessary ingredients to reliably predict charge transfer excitations. However, the range separating parameter (ω) is system-dependent and its optimization, including the chemical environment, is intricate. Refaely-Abramson et al. demonstrated that the gap renormalization in molecular crystals, a ground state property, can be represented by an OT-RSH functional screened by ɛstatic [Phys. Rev. B 88, 081204(R) (2013)], the zero frequency scalar dielectric constant. In this study, we propose the use of an OT-RSH functional screened by the scalar dielectric constant in the high frequency limit (OT-sRSH), ɛ∞, an appropriate constraint for vertical ionization energies or excitations in a dielectric environment. We have performed calculations for S,N-heteroacene derivatives in tetrahydrofuran and dichloromethane. The "unscreened" OT-RSH functional tends to underestimate experimental ionization potentials (IPs) and optical gaps (Egs) by up to 1.5 and 0.5 eV, respectively. In contrast, OT-sRSH functional calculations underestimate IPs and Egs by only 0.4 and 0.2 eV. We also compared the OT-sRSH results to explicitly solvated OT-RSH functional calculations for oligothiophenes in dioxane, benzene in ammonia, and methylene blue in water. We observe that both the approaches perform similarly for weakly interacting intermolecular systems and deviate for solvent-solute interacting systems, as expected. In conclusion, the OT-sRSH functional can describe molecular systems with environmental polarization effects accurately, a step toward describing realistic molecular systems.

Nature and energetics of low-lying excited singlets/triplets and intersystem crossing rates in selone analogs of perylenediimide: A theoretical perspective

The structural rigidity and chemical diversity of the highly fluorescent perylenediimide (PDI) provide wide opportunities for developing triplet photosensitizers with sufficiently increased energy efficiency. Remarkably high intersystem crossing (ISC) rates with a complete fluorescence turn-off reported recently for several thione analogs of PDI due to substantially large spin-orbit coupling garners huge attention to develop other potential analogs. Here, several selone analogs of PDI, denoted as mSe-PDIs (m = 1-4) with varied Se content and positions, are investigated to provide a comprehensive and comparative picture down the group-16 using density functional theory (DFT) and time-dependent DFT implementing optimally tuned range-separated hybrid in toluene dielectric. All mSe-PDIs are confirmed to be dynamically stable and also thermodynamically feasible to synthesize from their oxygen and thione congeners. The first excited-state singlet (S1) of mSe-PDI with relatively low Se-content (m = 1, 2) is of nπ* character with an expected fluorescence turn-off. Whereas, the ππ* nature of the S1 for 3Se-PDI and 4Se-PDI suggests a possible fluorescence turn-on in the absence of any other active nonradiative deactivation pathways. However, ∼4-6 orders greater ISC rates (∼1012-1014 s-1) than the fluorescence ones (∼108 s-1) for all mSe-PDIs signify highly efficient triplet harvest. Importantly, significantly higher ISC rates for these mSe-PDIs than their thione congeners render them efficient triplet photosensitizers.

Does the Intersystem Crossing Rate of β-Iodinated Phosphorus Corrole Depend on Iodine Numbers and/or Positions?

The heavy-atom effect is known to enhance the intersystem crossing (ISC) in organic molecular systems. Effects of iodine numbers and positions on the ISC rate of a few meso-difluorophenyl substituted β-iodinated phosphorus corroles (PCs) with axially ligated fluorine atoms (mI-FPC; m = 1-4) are studied using a time-dependent optimally tuned range-separated hybrid. Solvent effects are accounted for through a polarizable continuum model with a toluene dielectric. Calculations suggest similar thermodynamic stability for all mI-FPCs and also reproduce the experimentally measured 0-0 energies for some of the freebase phosphorus corrole (FPC) systems studied here. Importantly, our results reveal that all mI-FPCs display 10 times larger ISC rate (∼109 s-1) than the fluorescence rate (∼108 s-1), and the higher ISC rate stems from the improved spin-orbit coupling (SOC) introduced by lighter heteroatoms like central P and biaxial F rather than the I heavy-atom effect. However, an enhanced SOC is found with increasing I content for El-Sayed forbidden ISC channels. Research findings reported in this study unveil the impact of light heteroatoms and heavy atoms in promoting ISC in several iodinated PCs, which help in designing visible-light-driven efficient triplet photosensitizers.

Solvent effects on the photoinduced charge separation dynamics of directly linked zinc phthalocyanine-perylenediimide dyads: a nonadiabatic dynamics simulation with an optimally tuned screened range-separated hybrid functional

Herein, we have employed a combination of the optimally tuned screened range-separated hybrid (OT-SRSH) functional, the polarizable continuum model (PCM), and nonadiabatic dynamics (NAMD) simulations to investigate the photoinduced dynamics of directly linked donor-acceptor dyads formed using zinc phthalocyanine (ZnPc) and perylenediimide (PDI), in which ZnPc is the donor while PDI is the acceptor. Our simulations aim to analyze the behavior of these dyads upon local excitation of the ZnPc moiety in the gas phase and in benzonitrile. Our findings indicate that the presence of a solvent can significantly influence the excited state dynamics of ZnPc-PDI dyads. Specifically, the polar solvent benzonitrile effectively lowers the vertical excitation energies of the charge transfer (CT) state from ZnPc to PDI. As a result, the energetic order of the locally excited (LE) states of ZnPc and the CT states is reversed compared to the gas phase. Consequently, the photoinduced electron transfer (PET) dynamics from ZnPc to PDI, which is absent in the gas phase, takes place in benzonitrile with a time constant of 10.4 ps. Importantly, our present work not only qualitatively agrees with experimental results but also provides in-depth insights into the underlying mechanisms responsible for the photoinduced dynamics of ZnPc-PDI. Moreover, this study emphasizes the importance of appropriately considering solvent effects in NAMD simulation of organic donor-acceptor systems, taking into account the distinct excited state dynamics observed in the gas phase and benzonitrile. Furthermore, the combination of the OT-SRSH functional, the PCM solvent model, and nonadiabatic dynamics simulations shows promise as a strategy for investigating the complex excited state dynamics of organic donor-acceptor systems in solvents. These findings will be valuable for the future design of novel organic donor-acceptor structures with improved performance.

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0-1 and only five systems (0.07%) with an absolute error larger than 0.01 a0-1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Understanding Primary Charge Separation in the Heliobacterial Reaction Center

The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.

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.

Poisoning density functional theory with benchmark sets of difficult systems

Large benchmark sets like GMTKN55 [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184] let us analyse the performance of density functional theory over a diverse range of systems and bonding types. However, assessing over a large and diverse set can miss cases where approaches fail badly, and can give a misleading sense of security. To this end we introduce a series of 'poison' benchmark sets, P30-5, P30-10 and P30-20, comprising systems with up to 5, 10 and 20 atoms, respectively. These sets represent the most difficult-to-model systems in GMTKN55. We expect them to be useful in developing new approximations, identifying weak points in existing ones, and to aid in selecting appropriate DFAs for computational studies involving difficult physics, e.g. catalysis.

Hybrid functionals with local range separation: Accurate atomization energies and reaction barrier heights

Range-separated hybrid approximations to the exchange-correlation density functional mix exact and semi-local exchange in a position-dependent manner. In their conventional form, the range separation is controlled by a constant parameter. Turning this constant into a density functional leads to a locally space-dependent range-separation function and thus a more powerful and flexible range-separation approach. In this work, we explore the self-consistent implementation of a local range-separated hybrid, taking into account a one-electron self-interaction correction and the behavior under uniform density scaling. We discuss different forms of the local range-separation function that depend on the electron density, its gradient, and the kinetic energy density. For test sets of atomization energies, reaction barrier heights, and total energies of atoms, we demonstrate that our best model is a clear improvement over common global range-separated hybrid functionals and can compete with density functionals that contain multiple empirical parameters. Promising results for equilibrium bond lengths, harmonic vibrational frequencies, and vertical ionization potentials further underline the potential and flexibility of our approach.

Linear fractional charge behavior in density functional theory through dielectric tuning of conductor-like polarizable continuum model

A property of exact density functional theory is linear fractional charge behavior as electrons are added or removed from a molecule. Typical density functional approximations (DFAs) exhibit delocalization error, which overstabilizes this fractional charge. Conversely, solvent corrections have been shown to erroneously destabilize this fractional charge. This work will show that an implicit solvent correction with a tuned dielectric can be used as an ad hoc correction to offset the delocalizing character of DFAs and achieve linear fractional charge behavior. While desirable, in principle, we find that this linear charge behavior degrades the vertical ionization energies reported by DFAs. Our results reveal that the localizing character of the solvent correction and the Hartree-Fock (HF) exchange offset each other. This helps explain the decreased ratios of HF exchange to DFA exchange in long-range hybrid tuning studies that use a solvent correction.

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.

Simplified tuning of long-range corrected density functionals for use in symmetry-adapted perturbation theory

Long considered a failure, second-order symmetry-adapted perturbation theory (SAPT) based on Kohn-Sham orbitals, or SAPT0(KS), can be resurrected for semiquantitative purposes using long-range corrected density functionals whose asymptotic behavior is adjusted separately for each monomer. As in other contexts, correct asymptotic behavior can be enforced via "optimal tuning" based on the ionization energy theorem of density functional theory, but the tuning procedure is tedious, expensive for large systems, and comes with a troubling dependence on system size. Here, we show that essentially identical results are obtained using a fast, convenient, and automated tuning procedure based on the size of the exchange hole. In conjunction with "extended" (X)SAPT methods that improve the description of dispersion, this procedure achieves benchmark-quality interaction energies, along with the usual SAPT energy decomposition, without the hassle of system-specific tuning.

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low-energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. Intensity borrowing by these spurious states can affect intensities of the valence excitations, altering electronic bandshapes. At higher excitation energies, it is difficult to distinguish spurious charge-transfer states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces "natural charge-transfer orbitals" that provide a means to isolate orbitals that are most likely to participate in a CTTS excitation. Projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation and be used as an automated subspace selection algorithm for projection-based embedding of a high-level description of the CTTS state in a lower-level description of its environment. We apply this method to an ab initio molecular dynamics trajectory of I^-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with the experiment, and only one-third of the water molecules in the I^-(H₂O)₉₆ simulation cell need to be described with LR-TDDFT to obtain excitation energies that are converged to <0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments.

Investigating Primary Charge Separation in the Reaction Center of Heliobacterium modesticaldum

We compute the primary charge separation step in the homodimeric reaction center (RC) of Heliobacterium modesticaldum from first principles. Using time-dependent density functional theory with the optimally tuned range-separated hybrid functional ωPBE, we calculate the excitations of a system comprising the special pair, the adjacent accessory bacteriochlorophylls, and the most relevant parts of the surrounding protein environment. The structure of the excitation spectrum can be rationalized from coupling of the individual bacteriochlorophyll pigments similar to molecular J- and H-aggregates. We find excited states corresponding to forward-charge transfer along the individual branches of the RC of H. modesticaldum. In the spectrum, these are located at an energy between the coupled Q_y and Q_x transitions. With ab initio Born-Oppenheimer molecular dynamics simulations, we reveal the influence of thermal vibrations on the excited states. The results show that the energy gap between the coupled Q_y and the forward-charge transfer excitations is ∼0.4 eV, which we consider to conflict with the concept of a direct transfer mechanism. Our calculations, however, reveal a certain spectral overlap of the forward-charge transfer and the coupled Q_x excitations. The reliability and robustness of the results are demonstrated by several numerical tests.

First principles theoretical spectroscopy of methylene blue: Between limitations of time-dependent density functional theory approximations and its realistic description in the solvent

Methylene blue [3,7-Bis(di-methylamino) phenothiazin-5-ium chloride] is a phenothiazine dye with applications as a sensitizer for photodynamic therapy, photoantimicrobials, and dye-sensitized solar cells. Time-dependent density functional theory (TDDFT), based on (semi)local and global hybrid exchange-correlation functionals, fails to correctly describe its spectral features due to known limitations for describing optical excitations of π-conjugated systems. Here, we use TDDFT with a non-empirical optimally tuned range-separated hybrid functional to explore the optical excitations of gas phase and solvated methylene blue. We compute solvated configurations using molecular dynamics and an iterative procedure to account for explicit solute polarization. We rationalize and validate that by extrapolating the optimized range separation parameter to an infinite amount of solvating molecules, the optical gap of methylene blue is well described. Moreover, this method allows us to resolve contributions from solvent-solute intermolecular interactions and dielectric screening. We validate our results by comparing them to first-principles calculations based on the GW+Bethe-Salpeter equation approach and experiment. Vibronic calculations using TDDFT and the generating function method account for the spectra's subbands and bring the computed transition energies to within 0.15 eV of the experimental data. This methodology is expected to perform equivalently well for describing solvated spectra of π-conjugated systems.

Beyond Koopmans' theorem: electron binding energies in disordered materials

The topical review focuses on calculating ionization energies (IE), or electronic polarons in quasi-particle terminology, in large disordered systems, e.g. for a solute dissolved in a molecular solvent. The simplest estimate of the ionization energy is provided by one-electron energies in the Hartree-Fock theory, but the calculated quantities are not accurate. Density functional theory as many-body theory provides a principal opportunity for calculating one-electron energies including correlation and relaxation effects, i.e. the true energies of electronic polarons. We argue that such a principal possibility materializes within the concept of optimally tuned range-separated hybrid functionals (OT-RSH). We describe various schemes for optimal tuning. Importantly, the OT-RSH scheme is investigated for systems capped with dielectric continuum, providing a consistent picture on the QM/dielectric boundary. Finally, some limitations and open issues of the OT-RSH approach are addressed.

Dielectric Screening Meets Optimally Tuned Density Functionals

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

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.

Using non-empirically tuned range-separated functionals with simulated emission bands to model fluorescence lifetimes

Fluorescence lifetimes were evaluated using TD-DFT under different approximations for the emitting molecule and various exchange-correlation functionals, such as B3LYP, BMK, CAM-B3LYP, LC-BLYP, M06, M06-2X, M11, PBE0, ωB97, ωB97X, LC-BLYP*, and ωB97X* where the range-separation parameters in the last two functionals were tuned in a non-empirical fashion. Changes in the optimised molecular geometries between the ground and electronically excited states were found to affect the quality of the calculated lifetimes significantly, while the inclusion of vibronic features led to further improvements over the assumption of a vertical electronic transition. The LC-BLYP* functional was found to return the most accurate fluorescence lifetimes with unsigned errors that are mostly within 1.5 ns of experimental values.

Calculation of linear and nonlinear optical properties of azobenzene derivatives with Kohn–Sham and coupled-cluster methods

A non-empirically tuned generalized Kohn–Sham functional allows the prediction of accurate low-energy excitation energies and linear polarizabilities. Second hyperpolarizabilities are not improved when compared to coupled-cluster benchmark data.

Effect of Solid-State Polarization on Charge-Transfer Excitations and Transport Levels at Organic Interfaces from a Screened Range-Separated Hybrid Functional

We develop a robust approach for the description of the energetics of charge-transfer (CT) excitations and transport levels at organic interfaces based on a screened range-separated hybrid (SRSH) functional. We find that SRSH functionals correctly capture the effect of solid-state electronic polarization on transport gap renormalization and on screening of the electrostatic electron-hole interaction. With respect to calculations based on nonscreened optimally tuned RSH (long-range corrected) functionals, the SRSH-based calculations can be performed for both isolated molecular complexes and systems embedded in a dielectric medium with the same range-separation parameter, which allows a clear physical interpretation of the results in terms of solid-state polarization without any perturbation of the molecular electronic structure. By considering weakly interacting donor/acceptor complexes of pentacene with C₆₀ and poly-3-hexylthiophene (P3HT) with PCBM, we show that this new approach provides CT-state energies that compare very well with experimental data.

Description of the Charge Transfer States at the Pentacene/C60 Interface: Combining Range-Separated Hybrid Functionals with the Polarizable Continuum Model

Density functional theory (DFT) approaches based on range-separated hybrid functionals are currently methods of choice for the description of the charge-transfer (CT) states in organic donor/acceptor solar cells. However, these calculations are usually performed on small-size donor/acceptor complexes and as result do not account for electronic polarization effects. Here, using a pentacene/C60 complex as a model system, we discuss the ability of long-range corrected (LCR) hybrid functionals in combination with the polarizable continuum model (PCM) to determine the impact of the solid-state environment on the CT states. The CT energies are found to be insensitive to the interactions with the dielectric medium when a conventional time-dependent DFT/PCM (TDDFT/PCM) approach is used. However, a decrease in the energy of the CT state in the framework of LRC functionals can be obtained by using a smaller range-separated parameter when going from an isolated donor/acceptor complex to the solid-state case.

Rational computing of energy levels for organic electronics: the case of 2-benzylidene-1,3-indandiones

To compute ionization potential and electron affinity (HOMO & LUMO levels by tradition), diffuse functions are redundant, while range-separated hybrid DFT is not always better than B3LYP. Solvent modelling is essential and can be cheap with CPCM.

Tuned range separated hybrid functionals for solvated low bandgap oligomers

The description of charge transfer excitations has long been a challenge to time dependent density functional theory. The recently developed concept of "optimally tuned range separated hybrid (OT-RSH) functionals" has proven to describe charge transfer excitations accurately in many cases. However, describing solvated or embedded systems is yet a challenge. This challenge is not only computational but also conceptual, because the tuning requires identifying a specific orbital, typically the highest occupied one of the molecule under study. For solvated molecules, this orbital may be delocalized over the solvent. We here demonstrate that one way of overcoming this problem is to use a locally projected self-consistent field diagonalization on an absolutely localized molecular orbital expansion. We employ this approach to determine ionization energies and the optical gap of solvated oligothiophenes, i.e., paradigm low gap systems that are of relevance in organic electronics. Dioxane solvent molecules are explicitly represented in our calculations, and the ambiguities of straightforward parameter tuning in solution are elucidated. We show that a consistent estimate of the optimal range separated parameter (ω) at the limit of bulk solvation can be obtained by gradually extending the solvated system. In particular, ω is influenced by the solvent beyond the first coordination sphere. For determining ionization energies, a considerable number of solvent molecules on the first solvation shell must be taken into account. We demonstrate that accurately calculating optical gaps of solvated systems using OT-RSH can be done in three steps: (i) including the chemical environment when determining the range-separation parameter, (ii) taking into account the screening due to the solvent, and (iii) using realistic molecular geometries.