Extension of local response dispersion method to excited-state calculation based on time-dependent density functional theory

How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models

Time-dependent double hybrids with spin-component or spin-opposite scaling to their second-order perturbative correlation correction have demonstrated competitive robustness in the computation of electronic excitation energies. Some of the most robust are those recently published by our group (M. Casanova-Páez, L. Goerigk, J. Chem. Theory Comput. 2021, 20, 5165). So far, the implementation of these functionals has not allowed correctly calculating their ground-state total energies. Herein, we define their correct spin-scaled ground-state energy expressions which enables us to test our methods on the noncovalent excited-state interaction energies of four aromatic excimers. A range of 22 double hybrids with and without spin scaling are compared to the reasonably accurate wavefunction reference from our previous work (A. C. Hancock, L. Goerigk, RSC Adv. 2023, 13, 35964). The impact of spin scaling is highly dependent on the underlying functional expression, however, the smallest overall errors belong to spin-scaled functionals with range separation: SCS- and SOS- ω PBEPP86, and SCS-RSX-QIDH. We additionally determine parameters for DFT-D3(BJ)/D4 ground-state dispersion corrections of these functionals, which reduce errors in most cases. We highlight the necessity of dispersion corrections for even the most robust TD-DFT methods but also point out that ground-state based corrections are insufficient to completely capture dispersion effects for excited-state interaction energies.

Excitation configuration analysis for divide-and-conquer excited-state calculation method using dynamical polarizability

The authors previously developed a divide-and-conquer (DC)-based non-local excited-state calculation method for large systems using dynamical polarizability [Nakai and Yoshikawa, J. Chem. Phys. 146, 124123 (2017)]. This method evaluates the excitation energies and oscillator strengths using information on the dynamical polarizability poles. This article proposes a novel analysis of the previously developed method to obtain further configuration information on excited states, including excitation and de-excitation coefficients of each excitation configuration. Numerical applications to simple molecules, such as ethylene, hydrogen molecule, ammonia, and pyridazine, confirmed that the proposed analysis could accurately reproduce the excitation and de-excitation coefficients. The combination with the DC scheme enables both the local and non-local excited states of large systems with an excited nature to be treated.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Efficient Calculation of the Dispersion Energy for Multireference Systems with Cholesky Decomposition: Application to Excited-State Interactions

Accurate and efficient prediction of dispersion interactions in excited-state complexes poses a challenge due to the complex nature of electron correlation effects that need to be simultaneously considered. We propose an algorithm for computing the dispersion energy in nondegenerate ground- or excited-state complexes with arbitrary spin. The algorithm scales with the fifth power of the system size due to employing Cholesky decomposition of Coulomb integrals and a recently developed recursive formula for density response functions of the monomers. As a numerical illustration, we apply the new algorithm in the framework of multiconfigurational symmetry adapted perturbation theory, SAPT(MC), to study interactions in dimers with localized excitons. The SAPT(MC) analysis reveals that the dispersion energy may be the main force stabilizing excited-state dimers.

Chromone-methanol clusters in the electronic ground and lowest triplet state: a delicate interplay of non-covalent interactions

Chromone offers two energetically almost equivalent docking sites for alcohol molecules, in which the hydroxyl group is hydrogen bonded to one of the free electron pairs of the carbonyl O atom. Here, the delicate balance between these two competing arrangements is studied by combining IR/R2PI and UV/IR/UV spectroscopy in a molecular beam supported by quantum-chemical calculations. Most interestingly, chromone undergoes an efficient intersystem crossing into the triplet manifold upon electronic excitation, so that the studies on aromatic molecule-solvent complexes are for the first time extended to such a cluster in a triplet state. As the lowest triplet state (T₁) is of ground state character, powerful energy decomposition approaches such as symmetry-adapted perturbation theory (SAPT) and local energy decomposition using the domain-based local pair natural orbital coupled-cluster method (DLPNO-CCSD(T)/LED) are applied. From the theoretical analysis we infer for the T₁ state a loss of planarity (puckering) of the 4-pyrone ring of the chromone unit, which considerably affects the interplay between different types of non-covalent interactions at the two possible binding sites.

Dispersion Interactions in Exciton-Localized States. Theory and Applications to π-π* and n-π* Excited States

We address the problem of intermolecular interaction energy calculations in molecular complexes with localized excitons. Our focus is on the correct representation of the dispersion energy. We derive an extended Casimir-Polder formula for direct computation of this contribution through second order in the intermolecular interaction operator V̂. An alternative formula, accurate to infinite order in V̂, is derived within the framework of the adiabatic connection (AC) theory. We also propose a new parametrization of the VV10 nonlocal correlation density functional, so that it corrects the CASSCF energy for the dispersion contribution and can be applied to excited-state complexes. A numerical investigation is carried out for benzene, pyridine, and peptide complexes with the local exciton corresponding to the lowest π–π* or n– π* states. The extended Casimir-Polder formula is implemented in the framework of multiconfigurational symmetry-adapted perturbation theory, SAPT(MC). A SAPT(MC) analysis shows that the creation of a localized exciton affects mostly the electrostatic component of the interaction energy of investigated complexes. Nevertheless, the changes in Pauli repulsion and dispersion energies cannot be neglected. We verify the performance of several perturbation- and AC-based methods. Best results are obtained with a range-separated variant of an approximate AC approach employing extended random phase approximation and CASSCF wave functions.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for Time-Dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Symmetry-Adapted Perturbation Theory Based on Multiconfigurational Wave Function Description of Monomers

We present a formulation of the multiconfigurational (MC) wave function symmetry-adapted perturbation theory (SAPT). The method is applicable to noncovalent interactions between monomers which require a multiconfigurational description, in particular when the interacting system is strongly correlated or in an electronically excited state. SAPT(MC) is based on one- and two-particle reduced density matrices of the monomers and assumes the single-exchange approximation for the exchange energy contributions. Second-order terms are expressed through response properties from extended random phase approximation (ERPA). The dispersion components of SAPT(MC) have been introduced in our previous works [HapkaM.J. Chem. Theory Comput.2019, 15, 1016−102730525591; HapkaM.J. Chem. Theory Comput.2019, 15, 6712–672331670950]. SAPT(MC) is applied either with generalized valence bond perfect pairing (GVB) or with complete active space self-consistent field (CASSCF) treatment of the monomers. We discuss two model multireference systems: the H₂ ··· H₂ dimer in out-of-equilibrium geometries and interaction between the argon atom and excited state of ethylene. Using the C₂H₄* ··· Ar complex as an example, we examine second-order terms arising from negative transitions in the linear response function of an excited monomer. We demonstrate that the negative-transition terms must be accounted for to ensure qualitative prediction of induction and dispersion energies and develop a procedure allowing for their computation. Factors limiting the accuracy of SAPT(MC) are discussed in comparison with other second-order SAPT schemes on a data set of small single-reference dimers.

Second-Order Dispersion Energy Based on Multireference Description of Monomers

We propose a method for calculating a second-order dispersion energy for weakly interacting multireference systems in arbitrary electronic states. It is based on response properties obtained from extended random phase approximation equations. The introduced formalism is general and requires only one- and two-particle reduced density matrices of monomers. We combine the new method with either generalized valence bond perfect pairing (GVB) or complete active space (CAS) self-consistent field description of the interacting systems. In addition to a general scheme, three approximations, leading to significant reduction of the computational cost, are developed by exploiting Dyall partitioning of the monomer Hamiltonians. For model multireference systems (H₂···H₂ and Be···Be) the method is accurate, unlike its single-reference-based counterpart. Neither GVB nor CAS description of single-reference monomers improves the dispersion energy with respect to the Hartree-Fock-based results.

The effect of electronic excitation on London dispersion

Atomic and molecular dispersion coefficients can now be calculated routinely using density-functional theory. In this work, we present the first determination of how electronic excitation affects molecular C6 London dispersion coefficients from the exchange-hole dipole moment (XDM) dispersion model. Excited states are typically found to have larger dispersion coefficients than the corresponding ground states, due to their more diffuse electron densities. A particular focus is both intramolecular and intermolecular charge-transfer excitations, which have high absorbance intensities and are important in organic dyes, light-emitting diodes, and photovoltaics. In these classes of molecules, the increase in C6 for the electron-accepting moiety is largely offset by a decrease in C6 for the electron-donating moiety. As a result, the change in dispersion energy for a chromophore interacting with neighbouring molecules in the condensed phase is minimal.

How do London Dispersion Interactions Impact the Photochemical Processes of Molecular Switches?

In the last two decades, linear-response time-dependent density functional theory (LR-TDDFT) has become one of the most widely used approaches for the computation of the excited-state properties of atoms and molecules. Despite its success in describing the photochemistry and the photophysics of a vast majority of molecular systems, its domain of applicability has been limited by several substantial drawbacks. Commonly identified problems of LR-TDDFT include the correct description of Rydberg states, charge-transfer excited states, doubly excited states, and nearly degenerate states. In addition to these widely recognized shortcomings, the approximate functionals used in LR-TDDFT are unable to fully describe London dispersion interactions. In this work, we aim at understanding the impact of van der Waals interactions on the properties of chemical systems beyond their electronic ground state. For this purpose, we compare the results of excited-state energy profiles and dynamic trajectories for the prototypical cis-stilbene molecule with its 3-3',5-5'-tetra-tert-butyl derivative. While the explicit treatment of London dispersion interactions results in negligible changes for the cis-stilbene, we show that these attractive forces have a substantial influence on the energetics and structural evolution of the substituted derivative. In the latter case, intramolecular dispersion interactions impact the outcome of the simulation qualitatively, leading to an increased preference for the photocyclization pathway. The methodological consequences of this work are not uniquely applicable to the illustrative stilbene case. In fact, this molecule is representative of a whole class of chemical situations, where dispersion forces dominate the interactions between the unexcited substituents of a photoexcited chromophore. This is, for instance, a common situation in organic photovoltaics where donor molecules are usually functionalized with long alkyl side chains to improve solubility and assembly.

Dispersion Interactions in QM/EFP

The dispersion energy term between quantum-mechanical (QM) and classical (represented by effective fragment potentials, EFP) subsystems is developed and implemented. A new formulation is based on long-range perturbation theory and uses dynamic polarizability tensors of the effective fragments and electric field integrals and orbital energies of the quantum-mechanical subsystem. No parametrization is involved. The accuracy of the QM-EFP dispersion energy is tested on a number of model systems; the average mean unsigned error is 0.8 kcal/mol or 13% with respect to the symmetry adapted perturbation theory on the S22 data set of noncovalent interactions. The computational cost of the dispersion energy computation is low compared to the self-consistent field calculation of the QM subsystem. The dispersion energy is sensitive to the level of theory employed for the QM part and to the electrostatic interactions in the system. The latter means that the dispersion interactions in the QM/EFP method are not purely two-body but have more complex many-body behavior.

Structure-based modeling of energy transfer in photosynthesis

We provide a minimal model for a structure-based simulation of excitation energy transfer in pigment-protein complexes (PPCs). In our treatment, the PPC is assembled from its building blocks. The latter are defined such that electron exchange occurs only within, but not between these units. The variational principle is applied to investigate how the Coulomb interaction between building blocks changes the character of the electronic states of the PPC. In this way, the standard exciton Hamiltonian is obtained from first principles and a hierarchy of calculation schemes for the parameters of this Hamiltonian arises. Possible extensions of this approach are discussed concerning (i) the inclusion of dispersive site energy shifts and (ii) the inclusion of electron exchange between pigments. First results on electron exchange within the special pair of photosystem II of cyanobacteria and higher plants are presented and compared with earlier results on purple bacteria. In the last part of this mini-review, the coupling of electronic and nuclear degrees of freedom is considered. First, the standard exciton-vibrational Hamiltonian is parameterized with the help of a normal mode analysis of the PPC. Second, dynamical theories are discussed that exploit this Hamiltonian in the study of dissipative exciton motion.

Theoretical study of the electronic excitations of free-base porphyrin-Ar2 van der Waals complexes

The intermolecular interaction of free-base porphine (FBP)-Ar2 and free-base tetraazaporphyrin (FBPz)-Ar2 van der Waals (vdW) complexes was calculated in the ground state and vertical excitations that correspond to the Q- and B-bands using the many-body wavefunction theory of the symmetry-adapted cluster-configuration interaction (SAC-CI) method and time-dependent density functional theory (TDDFT). For the 1(1)B3u state of FBP-Ar2 a blueshift (high-energy shift) of excitation energy was calculated using the SAC-CI method; such a blueshift was not obtained by TDDFT calculations. This calculated blueshift corresponds to the experimentally observed blueshift in the Qx-band of FBP for FBP-Arn complexes. For FBPz-Ar2, blueshifts of the Q-band were not obtained using SAC-CI and TDDFT. These behaviors of the energy shift of the Q-bands could not be explained by the point dipole-point dipole interaction model. Large redshifts (low-energy shift) were obtained for the B-band states (2(1)B3u and 2(1)B2u) of FBP and FBPz. The energy shift showed the inverse sixth-power dependence on the intermolecular distance. The point dipole-point dipole interaction model can describe the redshift of the B-band. For the excited states that exhibit large redshifts, the TDDFT can qualitatively describe the vdW interaction in the excited states by supermolecular calculations. The solvatochromic shifts for FBP and FBPz in an Ar matrix were examined by the linear-response polarizable continuum model and TDDFT. The magnitude of calculated solvatochromic redshifts is proportional to the square of the transition dipole moment.