Assessment of the ΔSCF density functional theory approach for electronic excitations in organic dyes

State-Specific Coupled-Cluster Methods for Excited States

We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.

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.

Aufbau Suppressed Coupled Cluster Theory for Electronically Excited States

We introduce an approach to improve single-reference coupled cluster theory in settings where the Aufbau determinant is absent from or plays only a small role in the true wave function. Using a de-excitation operator that can be efficiently hidden within a similarity transform, we create a coupled cluster wave function in which de-excitations work to suppress the Aufbau determinant and produce wave functions dominated by other determinants. Thanks to an invertible and fully exponential form, the approach is systematically improvable, size consistent, size extensive, and, interestingly, size intensive in a granular way that should make the adoption of some ground state techniques, such as local correlation, relatively straightforward. In this initial study, we apply the general formalism to create a state-specific method for orbital-relaxed, singly excited states. We find that this approach matches the accuracy of similar-cost equation-of-motion methods in valence excitations while offering improved accuracy for charge transfer states. We also find the approach to be more accurate than excited-state-specific perturbation theory in both types of states.

Seniority and Hierarchy Configuration Interaction for Radicals and Excited States

Hierarchy configuration interaction (hCI) has recently been introduced as an alternative configuration interaction (CI) route combining excitation degree and seniority number and has been shown to efficiently recover both dynamic and static correlations for closed-shell molecular systems [ J. Phys. Chem. Lett. 2022, 13, 4342]. Here we generalize hCI for an arbitrary reference determinant, allowing calculations for radicals and excited states in a state-specific way. We gauge this route against excitation-based CI (eCI) and seniority-based CI (sCI) by evaluating how different ground-state properties of radicals converge to the full CI limit. We find that hCI outperforms or matches eCI, whereas sCI is far less accurate, in line with previous observations for closed-shell molecules. Employing second-order Epstein-Nesbet (EN2) perturbation theory as a correction significantly accelerates the convergence of hCI and eCI. We further explore various hCI and sCI models to calculate the excitation energies of closed- and open-shell systems. Our results underline that the choice of both the reference determinant and the set of orbitals drives the fine balance between correlation of ground and excited states. State-specific hCI2 and higher-order models perform similarly to their eCI counterparts, whereas lower orders of hCI deliver poor results unless supplemented by the EN2 correction, which substantially improves their accuracy. In turn, sCI1 produces decent excitation energies for radicals, encouraging the development of related seniority-based coupled-cluster methods.

Photophysics of Solvated Molecules: Computational Protocol Combining Nuclear Ensemble and Nonequilibrium State-Specific Solvation Methods

Efficient calculations of the photophysical properties of molecules are essential for both understanding experimental results and accelerating materials discovery through computational simulations. However, to achieve accurate results, the effects of the surrounding medium must be taken into account. Here, we present a computational protocol that combines the nuclear ensemble method with a nonequilibrium state-specific polarizable continuum model to simulate absorption, fluorescence, phosphorescence, and intersystem crossing processes. Additionally, we introduced an extrapolation strategy that enables predictions for multiple solvents without incurring additional computational costs. We demonstrate the method's effectiveness by modeling the photophysical properties of a molecule that exhibits thermally activated delayed fluorescence, showcasing how these properties vary with solvent polarity. We also provide insight into the relationship between the solvent and photophysics by using ensemble analysis to rationalize simulation results. Furthermore, we introduce a metric for the intensity of the charge transfer character of electronic states and demonstrate how vibrations can significantly mix the electronic character of excited states. Overall, this work presents a computational approach that offers new insights into the photophysics of molecules and has the potential to advance materials discovery.

Excited-State-Specific Pseudoprojected Coupled-Cluster Theory

We present an excited-state-specific coupled-cluster approach in which both the molecular orbitals and cluster amplitudes are optimized for an individual excited state. The theory is formulated via a pseudoprojection of the traditional coupled-cluster wavefunction that allows correlation effects to be introduced atop an excited-state mean field starting point. The approach shares much in common with ground-state CCSD, including size extensivity and an N6 cost scaling. Preliminary numerical tests show that, when augmented with N5 cost perturbative corrections for key terms, the method can improve over excited-state-specific second-order perturbation theory in valence, charge transfer, and Rydberg states.

The application of the SCAN density functional to color centers in diamond

Detailed characterization of deep-level color centers requires understanding their electronic and atomic structure, which is most commonly investigated utilizing the Kohn-Sham density functional theory. Standard semilocal functionals based on the generalized gradient approximation (GGA) are inclined toward an imprecise quantitative description of defects' electronic structure. Hybrid functionals provide an improved prediction of electronic properties, albeit at a much higher computational cost. In this work, we test the newly developed Strongly Constrained and Appropriately Normed (SCAN) family of meta-GGA density functionals for selected color centers in diamond. In particular, we study nitrogen-, silicon-, germanium-, and tin-vacancy centers that have been recently investigated for their use in quantum technological applications. We show that SCAN and its derivatives, the rSCAN and r2SCAN functionals, significantly improve the calculated energies of optical transitions within the delta-self-consistent-field approach, almost reaching the accuracy of the hybrid Heyd-Scuseria-Ernzerhof (HSE) functional. In the case of the NV- center, we also show that the SCAN family of functionals improves the description of the adiabatic potential energy surfaces compared to both GGA and hybrid functionals, improving calculated luminescence lineshapes. As a result of these findings, we recommend using the SCAN family of functionals as a promising alternative for studying color centers in solids.

Reverse Designing the Wavelength-Specific Thermally Activation Delayed Fluorescent Molecules Using a Genetic Algorithm Coupled with Cheap QM Methods

Genetic algorithm (GA) optimization coupled with the semiempirical intermediate neglect of differential overlap (INDO)/CIS method is presented to inversely design the red thermally activation delayed fluorescent (TADF) molecules. According to the predefined donor-acceptor (DA) library to build an AD_n-type TADF candidate, we utilized the chemical notation language SMILES code to generate a TADF molecule and apply the RDKit program to produce the initial 3D molecular structure. A combined fitness function is proposed to evaluate the performance of the functional-lead TADF molecule. The fitness function includes three key parameters, i.e., the emission wavelength, the energy gap (ΔE_ST) between the lowest singlet (S₁)- and triplet (T₁)-excited states, and the oscillator strengths for electron transition from S₀ and S₁. A cheap QM method, i.e., INDO/CIS, on the basis of an xTB-optimized molecular geometry is applied to quickly calculate the fitness function. Finally, the GA approach is utilized to globally search for the wavelength-specific TADF molecules under our predefined DA library, and the optimum 630 nm red and 660 nm deep red TADF molecules are inversely designed according to the evolution of molecular fitness functions.

State-Specific Configuration Interaction for Excited States

We introduce and benchmark a systematically improvable route for excited-state calculations, labeled state-specific configuration interaction (ΔCI), which is a particular realization of multiconfigurational self-consistent field and multireference configuration interaction. Starting with a reference built from optimized configuration state functions, separate CI calculations are performed for each targeted state (hence, state-specific orbitals and determinants). Accounting for single and double excitations produces the ΔCISD model, which can be improved with second-order Epstein-Nesbet perturbation theory (ΔCISD+EN2) or a posteriori Davidson corrections (ΔCISD+Q). These models were gauged against a vast and diverse set of 294 reference excitation energies. We have found that ΔCI is significantly more accurate than standard ground-state-based CI, whereas close performances were found between ΔCISD and EOM-CC2 and between ΔCISD+EN2 and EOM-CCSD. For larger systems, ΔCISD+Q delivers more accurate results than EOM-CC2 and EOM-CCSD. The ΔCI route can handle challenging multireference problems, singly and doubly excited states, from closed- and open-shell species, with overall comparable accuracy and thus represents a promising alternative to more established methodologies. In its current form, however, it is reliable only for relatively low-lying excited states.

Gold(i)-containing light-emitting molecules with an inverted singlet–triplet gap†

Delayed fluorescence from molecules with an inverted singlet–triplet gap (DFIST) is the consequence of the unusual reverse order of the lowest excited singlet (S₁) and triplet (T₁) states of thermally activated delayed fluorescence (TADF) emitters. Heptazine (1,3,4,6,7,9,9b-heptaazaphenalene) derivatives have an inverted singlet–triplet gap thanks to the combination of multiple resonance (MR) effects and a significant double excitation character. Here, we study computationally the effect of gold(i) metalation and coordination on the optical properties of heptazine (molecule 4) and the phosphine-functionalized 2,5,8-tris(dimethylphosphino)heptazine derivatives (molecules 1–3). Ab initio calculations at the approximate second-order coupled cluster (CC2) and extended multiconfigurational quasi degenerate perturbation theory at the second order (XMC-QDPT2) levels show that molecules 1–4 have an inverted singlet–triplet gap due to the alternating spatial localization of the electron and hole of the exciton in the heptazine core. A non-vanishing one-electron spin–orbit coupling operator matrix element between T₁ and and a fast S₁ ← T₁ intersystem crossing rate constant (k_ISC) calculated at the XMC-QDPT2(12,12) level of theory for molecule 4 suggest that this new family of complexes may be the first organometallic DFIST emitters reported.

Substitution with gold(i)-containing moieties results in non-vanishing oscillator strengths and spin–orbit coupling leading to fast intersystem crossing in light-emitting heptazine derivates with an inverted singlet–triplet gap.

Combining Time-Dependent Density Functional Theory and the ΔSCF Approach for Accurate Core-Electron Spectra

Spectroscopies that probe electronic excitations from core levels into unoccupied orbitals, such as X-ray absorption spectroscopy and electron energy loss spectroscopy, are widely used to gain insight into the electronic and chemical structure of materials. To support the interpretation of experimental spectra, we assess the performance of a first-principles approach that combines linear-response time-dependent density (TDDFT) functional theory with the Δ self-consistent field (ΔSCF) approach. In particular, we first use TDDFT to calculate the core-level spectrum and then shift the spectrum such that the lowest excitation energy from TDDFT agrees with that from ΔSCF. We apply this method to several small molecules and find encouraging agreement between calculated and measured spectra.

Diabatic Decomposition Perspective on the Role of Charge Transfer and Local Excitations in Thermally Activated Delayed Fluorescence

Thermally activated delayed fluorescence (TADF) is a phenomenon that relies on the upconversion of triplet excitons to singlet excitons by means of reverse intersystem crossing (rISC). It has been shown both experimentally and theoretically that the TADF mechanism depends on the interplay between charge transfer and local excitations. However, the difference between the diabatic and adiabatic character of the involved excited states is rarely discussed in the literature. Here we develop a diabatization procedure to implement a four-state model Hamiltonian to a set of TADF molecules. We provide physical interpretations of the Hamiltonian elements and show their dependence on the electronic state of the equilibrium geometry. We also demonstrate how vibrations affect the TADF efficiency by modifying the diabatic decomposition of the molecule. Finally, we provide a simple model that connects the diabatic Hamiltonian to the electronic properties relevant to TADF and show how this relationship translates into different optimization strategies for rISC, fluorescence, and overall TADF performance.

On the origin of the inverted singlet-triplet gap of the 5th generation light-emitting molecules

Excitation energies of the lowest singlet and triplet state of molecules whose first excited singlet state lies energetically below the first triplet state have been studied computationally at (time-dependent) density functional theory, coupled-cluster, and second-order multiconfiguration perturbation theory levels. The calculations at the ab initio levels show that the singlet-triplet gap is inverted as compared to the one expected from Hund's rule, whereas all density functionals yield the triplet state as the lowest excited state. Double excitations responsible for the inverted singlet-triplet gap have been identified. Employing the spin-flip and ΔSCF methods, singlet-triplet inversion was obtained at the density functional theory level for some of the studied molecules.

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.

Quantum simulations of thermally activated delayed fluorescence in an all-organic emitter

We investigate the prototypical NAI-DMAC thermally activated delayed fluorescence (TADF) emitter in the gas phase- and high-packing fraction limits at finite temperature, by combining first principles molecular dynamics with a quantum thermostat to account for nuclear quantum effects (NQE). We find a weak dependence of the singlet-triplet energy gap (ΔE_ST) on temperature in both the solid and the molecule, and a substantial effect of packing. While the ΔE_ST vanishes in the perfect crystal, it is of the order of ∼0.3 eV in the molecule, with fluctuations ranging from 0.1 to 0.4 eV at 300 K. The transition probability between the HOMOs and LUMOs has a stronger dependence on temperature than the singlet-triplet gap, with a desirable effect for thermally activated fluorescence; such temperature effect is weaker in the condensed phase than in the molecule. Our results on ΔE_ST and oscillator strengths, together with our estimates of direct and reverse intersystem crossing rates, show that optimization of packing and geometrical conformation is critical to increase the efficiency of TADF compounds. Our findings highlight the importance of considering thermal fluctuations and NQE to obtain robust predictions of the electronic properties of NAI-DMAC.

Robust ΔSCF calculations with direct energy functional minimization methods and STEP for molecules and materials

The direct energy functional minimization method using the orbital transformation (OT) scheme in the program package CP2K has been employed for Δ self-consistent field (ΔSCF) calculations. The OT method for non-uniform molecular orbitals occupations allows us to apply the ΔSCF method for various kinds of molecules and periodic systems. Vertical excitation energies of heteroaromatic molecules and condensed phase systems, such as solvated ethylene and solvated uracil obeying periodic boundary conditions, are reported using the ΔSCF method. In addition, a Re-phosphate molecule attached to the surface of anatase (TiO₂) has been investigated. Additionally, we have implemented a recently proposed state-targeted energy projection ΔSCF algorithm [K. Carter-Fenk and J. M. Herbert, J. Chem. Theory Comput. 16(8), 5067-5082 (2020)] for diagonalization based SCF in CP2K. It is found that the OT scheme provides a smooth and robust SCF convergence for all investigated excitation energies and (non-)periodic systems.

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.

Origin of the Failure of Density Functional Theories in Predicting Inverted Singlet-Triplet Gaps

Recent experimental and theoretical studies have shown several new organic molecules that violate Hund’s rule and have the first singlet excited state lower in energy than the first triplet excited state. While many correlated single reference wave function methods have successfully predicted excited-state energetics of these low-lying states, conventional linear-response time-dependent density functional theory (TDDFT) fails to predict the correct excited-state energy ordering. In this article, we have explored the performance of combined DFT and wave function methods like doubles-corrected TDDFT and multiconfiguration pair-density functional theory for the calculation of inverted singlet–triplet gaps. We have also tested the performance of the excited-state DFT (eDFT) method for this problem. Our results have shown that it is possible to obtain inverted singlet–triplet gaps both by using doubles-corrected TDDFT with a proper choice of double-hybrid functionals or by using eDFT.

The photodissociation of solvated cyclopropanone and its hydrate explored via non-adiabatic molecular dynamics using ΔSCF

The decay of cyclopropanone is a typical example of a photodecomposition process. Ethylene and carbon monoxide are formed following the excitation to the first singlet excited state through a symmetrical or asymmetrical pathway. The results obtained with non-adiabatic molecular dynamics (NAMD) using the delta self-consistent field (ΔSCF) method correspond well to previous experimental and multireference theoretical studies carried out in the gas phase. Moreover, this efficient methodology allows NAMD simulations of cyclopropanone in aqueous solution to be performed, which reveal analogue deactivation mechanisms, but a shorter lifetime and reduced photodissociation as compared to the gas-phase. The excited state dynamics of cyclopropanone hydrate, an enzyme inhibitor, in an aqueous environment are reported as well. Cyclopropanone hydrate strongly interacts with the surrounding solvent via the formation of hydrogen bonds. Excitation to the first singlet excited state shows an asymmetric pathway with cyclopropanone hydrate and propionic acid as the main photoproducts.

The lifetime and photodissociation of cyclopropanone are reduced in aqueous solution, while the excitation of solvated cyclopropanone hydrate yields a range of photoproducts.

Taking Advantage of a Systematic Energy Non-linearity Error in Density Functional Theory for the Calculation of Electronic Energy Levels

We present an approximate approach for the calculation of ionization potential (IP) and electron affinity (EA) by exploiting the complementary energy non-linearity errors for a species M and its one-electron-ionized counterpart (M⁺). Reasonable IPs and EAs are thus obtained by averaging the orbital energies of M and M⁺, even with a low-level method such as BLYP/6-31G(d). By combining the corrected IPs and EAs, we can further obtain reasonable excitation energies. The errors in uncorrected valence IPs and uncorrected virtual-orbital energies show systematic trends. These characteristics provide a convenient and computationally efficient avenue for qualitative estimation of these properties with single corrections for multiple IPs and excitation energies.

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.

Ensemble Density Functional Theory of Neutral and Charged Excitations : Exact Formulations, Standard Approximations, and Open Questions

Recent progress in the field of (time-independent) ensemble density-functional theory (DFT) for excited states are reviewed. Both Gross-Oliveira-Kohn (GOK) and N-centered ensemble formalisms, which are mathematically very similar and allow for an in-principle-exact description of neutral and charged electronic excitations, respectively, are discussed. Key exact results, for example, the equivalence between the infamous derivative discontinuity problem and the description of weight dependencies in the ensemble exchange-correlation density functional, are highlighted. The variational evaluation of orbital-dependent ensemble Hartree-exchange (Hx) energies is discussed in detail. We show in passing that state-averaging individual exact Hx energies can lead to severe (although solvable) v-representability issues. Finally, we explore the possibility of using the concept of density-driven correlation, which has been introduced recently and does not exist in regular ground-state DFT, for improving state-of-the-art correlation density-functional approximations for ensembles. The present review reflects the efforts of a growing community to turn ensemble DFT into a rigorous and reliable low-cost computational method for excited states. We hope that, in the near future, this contribution will stimulate new formal and practical developments in the field.

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.

Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules

Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the M_S = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.

A Simple Effective Δ SCF Method for Computing Optical Gaps in Organic Chromophores

Photoluminescence effects in organic chromophores are of significant importance and requires precise description of low lying excited states. In this communication, we put forward an alternative time-independent DFT scheme for computing lowest single-particle excitation energy, especially for singlet excited state. This adopts a recently developed "virial"-theorem based model of singlet-triplet splitting which requires a DFT calculation on closed shell ground state and a restricted open-shell triplet excited state, followed by a simple 2 e - integral evaluation. This produces vertical excitation energies in small molecules, linear and non-linear polycyclic aromatic hydrocarbon and organic dyes in comparable accuracy to the TDDFT. We also explore the functional dependency of present method with three different functionals (B3LYP, wB97X and CAM-B3LYP) for polyenes and linear acenes. A systematic comparison with literature value illustrates the validity and usefulness of the present scheme in determining optical gap with fair computational cost.

Method for Calculating Excited Electronic States Using Density Functionals and Direct Orbital Optimization with Real Space Grid or Plane-Wave Basis Set

A direct orbital optimization method is presented for density functional calculations of excited electronic states using either a real space grid or a plane-wave basis set. The method is variational, provides atomic forces in the excited states, and can be applied to Kohn-Sham (KS) functionals as well as orbital-density-dependent (ODD) functionals including explicit self-interaction correction. The implementation for KS functionals involves two nested loops: (1) An inner loop for finding a stationary point in a subspace spanned by the occupied and a few virtual orbitals corresponding to the excited state; (2) an outer loop for minimizing the energy in a tangential direction in the space of the orbitals. For ODD functionals, a third loop is used to find the unitary transformation that minimizes the energy functional among occupied orbitals only. Combined with the maximum overlap method, the algorithm converges in challenging cases where conventional self-consistent field algorithms tend to fail. The benchmark tests presented include two charge-transfer excitations in nitrobenzene and an excitation of CO to degenerate π* orbitals where the importance of complex orbitals is illustrated. The application of this method to several metal-to-ligand charge-transfer and metal-centered excited states of an Fe^II photosensitizer complex is described, and the results are compared to reported experimental estimates. This method is also used to study the effect of the Perdew-Zunger self-interaction correction on valence and Rydberg excited states of several molecules, both singlet and triplet states, and the performance compared to semilocal and hybrid functionals.

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.

Perturbation theory made efficient and effective for predictions of ionization potential and electron affinity

Ionization potential and electron affinity are essential molecular properties. The most straightforward method is to calculate them by taking the total energy differences of the initial and final states according to the definition. However, it often suffers from a serious convergence problem due to the requirement of the self-consistent field (SCF) calculations for the ionic states with non-Aufbau choices of occupations. In the present work, we have constructed a theoretical framework in view of perturbation theory to bypass the SCF calculations of the ionic states. To address the imbalance issue that arises from the precisely treated neutral ground state followed by the truncated perturbative treatment of the ionic states, an accurate yet effective method has been developed here, which adds back some terms from the higher order perturbations into the lower order to cancel out the most computationally cost terms in the truncated expansion, thus reaching a better convergence with less computation. The validity of the present methodology has been tested out by applying it to the Hartree-Fock (HF) method in combination with the correlation effect described at the second-order Møller-Plesset level in a frozen-orbital approximation. All the derivations in this work are given in a general framework, which are applicable not only to HF but also to a wide range of density functional theory methods from semi-local functionals to hybrid and doubly hybrid functionals.

Donor acceptor fluorophores: synthesis, optical properties, TD-DFT and cytotoxicity studies

Donor-π-acceptor (D-π-A) fluorophores consisting of a donor unit, a π linker, and an acceptor moiety have attracted attention in the last decade. In this study, we report the synthesis, characterization, optical properties, TD-DFT, and cytotoxicity studies of 17 near infrared (NIR) D-π-A analogs which have not been reported so far to the best of our knowledge. These fluorophores have chloroacrylic acid as the acceptor unit and various donor units such as indole, benzothiazole, benzo[e]indole, and quinoline. The fluorophores showed strong absorption in the NIR (700-970 nm) region due to their enhanced intramolecular charge transfer (ICT) between chloroacrylic acid and the donor moieties connected with the Vilsmeier-Haack linker. The emission wavelength maxima of the fluorophores were in between 798 and 870 nm. Compound 20 with a 4-quinoline donor moiety showed an emission wavelength above 1000 nm in the NIR II window. The synthesized fluorophores were characterized by 1H NMR and 13C NMR, and their optical properties were studied. Time dependent density functional theory (TD-DFT) calculations showed that the charge transfer occurs from the donor groups (indole, benzothiazole, benzo[e]indole, and quinoline) to the acceptor chloroacrylic acid moiety. Fluorophores with [HOMO] to [LUMO+1] transitions were shown to possess a charge separation character. The cytotoxicity of selected fluorophores, 4, 7, 10 and 12 was investigated against breast cancer cell lines and they showed better activity than the anti-cancer agent docetaxel.

Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer

Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time-dependent density functional theory, and exorbitant computational cost when combined with high-level wavefunction-based methods. Here, we introduce a density functional theory (DFT)-based diabatization scheme that directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs), which we denote as Δ-ALMO(MSDFT2). We demonstrate that our method, which combines ALMO calculations with the ΔSCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.

An ab initio exciton model for singlet fission

We present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet-multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

Third-Order Møller-Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error [Bertels, L. W.; J. Phys. Chem. Lett. 2019, 10, 4170 4176]. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N⁶) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

Neutral excitation density-functional theory: an efficient and variational first-principles method for simulating neutral excitations in molecules

We introduce neutral excitation density-functional theory (XDFT), a computationally light, generally applicable, first-principles technique for calculating neutral electronic excitations. The concept is to generalise constrained density functional theory to free it from any assumptions about the spatial confinement of electrons and holes, but to maintain all the advantages of a variational method. The task of calculating the lowest excited state of a given symmetry is thereby simplified to one of performing a simple, low-cost sequence of coupled DFT calculations. We demonstrate the efficacy of the method by calculating the lowest single-particle singlet and triplet excitation energies in the well-known Thiel molecular test set, with results which are in good agreement with linear-response time-dependent density functional theory (LR-TDDFT). Furthermore, we show that XDFT can successfully capture two-electron excitations, in principle, offering a flexible approach to target specific effects beyond state-of-the-art adiabatic-kernel LR-TDDFT. Overall the method makes optical gaps and electron-hole binding energies readily accessible at a computational cost and scaling comparable to that of standard density functional theory. Owing to its multiple qualities beneficial to high-throughput studies where the optical gap is of particular interest; namely broad applicability, low computational demand, and ease of implementation and automation, XDFT presents as a viable candidate for research within materials discovery and informatics frameworks.

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

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.