Reference Energies for Double Excitations

Automated Multireference Vertical Excitations for Transition-Metal Compounds

Excited states of transition metal complexes are generally strongly correlated due to the near-degeneracy of the metal d orbitals. Consequently, electronic structure calculations of such species often necessitate multireference approaches. However, widespread use of multireference methods is hindered due to the active space selection problem, which has historically required system-specific chemical knowledge and a trial-and-error approach. Here, we address this issue with an automated method combining the approximate pair coefficient (APC) scheme for estimating orbital entropies with the discrete variational selection (DVS) approach for evaluating active space quality. We apply DVS-APC to the calculation of 67 vertical excitations in transition metal diatomics as well as to two larger complexes. We show DVS-APC generated active spaces yield NEVPT2 mean absolute errors of 0.18 eV, in line with previous accuracies obtained for organic systems, but larger than errors achieved with hand-selected active spaces (0.14 eV). If instead of using DVS we identify the best results from our trial wave functions, we find improved performance (mean absolute error of 0.1 eV) over the manually selected results. We highlight this deviation between DVS and hand selected active spaces as a possible measure of bias introduced when hand selecting active spaces. However, we find that multiconfiguration pair-density functional theory (MC-PDFT) using the tPBE and tPBE0 functionals is roughly 0.15 eV less accurate than NEVPT2 across this class of diatomic systems, potentially accounting for the decreased performance of DVS-APC, which uses MC-PDFT energies to select between active spaces. We also showcase an ability to "down-sample" the DVS-APC wave functions using natural orbital occupancies to achieve smaller minimal active spaces which retain the accuracy of the larger starting active spaces. Finally, DVS-APC and tPBE0 are proven to be effective when applied to modeling excited states in two larger transition metal complexes, suggesting that the transition metal diatomics may be a particular outstanding challenge for DVS-APC and MC-PDFT approaches.

Accurate and efficient prediction of double excitation energies using the particle-particle random phase approximation

Double excitations are crucial to understanding numerous chemical, physical, and biological processes, but accurately predicting them remains a challenge. In this work, we explore the particle-particle random phase approximation (ppRPA) as an efficient and accurate approach for computing double excitation energies. We benchmark ppRPA using various exchange-correlation functionals for 21 molecular systems and two point defect systems. Our results show that ppRPA with functionals containing appropriate amounts of exact exchange provides accuracy comparable to high-level wave function methods such as CCSDT and CASPT2, with significantly reduced computational cost. Furthermore, we demonstrate the use of ppRPA starting from an excited (N - 2)-electron state calculated by ΔSCF for the first time, as well as its application to double excitations in bulk periodic systems. These findings suggest that ppRPA is a promising tool for the efficient calculation of double and partial double excitation energies in both molecular and bulk systems.

QUEST#4X: An Extension of QUEST#4 for Benchmarking Multireference Wave Function Methods

Given a number of data sets for evaluating the performance of single reference methods for the low-lying excited states of closed-shell molecules, a comprehensive data set for assessing the performance of multireference methods for the low-lying excited states of open-shell systems is still lacking. For this reason, we propose an extension (QUEST#4X) of the radical subset of QUEST#4 (J. Chem. Theory Comput. 2020, 16, 3720) to cover 110 doublet and 39 quartet excited states. Near-exact results obtained by iterative configuration interaction with selection and second-order perturbation correction (iCIPT2) are taken as benchmark to calibrate static-dynamic-static configuration interaction (SDSCI) and static-dynamic-static second-order perturbation theory (SDSPT2), which are minimal MRCI and CI-like perturbation theory, respectively. It is found that SDSCI is very close in accuracy to internally contracted multireference configuration interaction with singles and doubles (ic-MRCISD), although its computational cost is just that of one iteration of the latter. Unlike most variants of MRPT2, SDSPT2 treats single and multiple states in the same way and performs similarly to multistate n-electron valence second-order perturbation theory (MS-NEVPT2). These findings put SDSCI and SDSPT2 on a firm basis.

Time-resolved vibronic spectra with nuclear-electronic orbital time-dependent configuration interaction

Time-resolved spectroscopy is an important tool for probing photochemically induced nonequilibrium dynamics and energy transfer. Herein, a method is developed for the ab initio simulation of vibronic spectra and dynamical processes. This framework utilizes the recently developed nuclear-electronic orbital time-dependent configuration interaction (NEO-TDCI) approach, which treats all electrons and specified nuclei quantum mechanically on the same footing. A strategy is presented for calculating time-resolved vibrational and electronic absorption spectra from any initial condition. Although this strategy is general for any TDCI implementation, utilizing the NEO framework allows for the explicit inclusion of quantized nuclei, as illustrated through the calculation of vibrationally hot spectra. Time-resolved spectra produced by either vibrational or electronic Rabi oscillations capture ground-state absorption, stimulated emission, and excited-state absorption between vibronic states. This methodology provides the foundation for fully ab initio simulations of multidimensional spectroscopic experiments.

Capturing the Elusive Curve-Crossing in Low-Lying States of Butadiene with Dressed TDDFT

A striking example of the need to accurately capture states of double-excitation character in molecules is seen in predicting photoinduced dynamics in small polyenes. Due to the coupling of electronic and nuclear motions, the dark 21Ag state, known to have double-excitation character, can be reached after an initial photoexcitation to the bright 11Bu state via crossings of their potential energy surfaces. However, the shapes of the surfaces are so poorly captured by most electronic structure methods, that the crossing is missed or substantially mis-located. We demonstrate that the frequency-dependent kernel of dressed TDDFT beyond Tamm-Dancoff successfully captures the curve-crossing, providing an energy surface close to the highly accurate but more expensive δ-CR-EOMCC(2,3) benchmark reference. This, along with its accurate prediction of the excitation character of the state makes dressed TDDFT a practical and accurate route to electronic structure quantities needed in modeling ultrafast dynamics in molecules.

A Cost-Effective Computational Strategy for the Electronic Layout Characterization of a Second Generation Light-Driven Molecular Rotary Motor in Solution

Light-driven molecular rotary motors are nanometric machines able to convert light into unidirectional motions. Several types of molecular motors have been developed to better respond to light stimuli, opening new avenues for developing smart materials ranging from nanomedicine to robotics. They have great importance in the scientific research across various disciplines, but a detailed comprehension of the underlying ultrafast photophysics immediately after photo-excitation, that is, Franck-Condon region characterization, is not fully achieved yet. For this aim, it is first required to rely on an accurate description at ab initio level of the system in this potential energy region before performing any further step, that is, dynamics. Thus, we present an extensive investigation aimed at accurately describing the electronic structure of low-lying electronic states (electronic layout) of a molecular rotor in the Franck-Condon region, belonging to the class of overcrowded alkenes: 9-(2-methyl-2,3-dihydro-1H-cyclopenta[a]naphthalen-1-ylidene)-9H-fluorene. This system was chosen since its photophysics is very interesting for a more general understanding of similar compounds used as molecular rotors, where low-lying electronic states can be found (whose energetic interplay is crucial in the dynamics) and where the presence of different substituents can tune the HOMO-LUMO gap. For this scope, we employed different theory levels within the time-dependent density functional theory framework, presenting also a careful comparison adopting very accurate post Hartree-Fock methods and characterizing also the different conformations involved in the photocycle. Effects on the electronic layout of different functionals, basis sets, environment descriptions, and the role of the dispersion correction were all analyzed in detail. In particular, a careful treatment of the solvent effects was here considered in depth, showing how the implicit solvent description can be accurate for excited states in the Franck-Condon region by testing both linear-response and state-specific formalisms. As main results, we chose two cost-effective (accurate but relatively cheap) theory levels for the ground and excited state descriptions, and we also verified how choosing these different levels of theory can influence the curvature of the potential via a frequency analysis of the normal modes of vibrations active in the Raman spectrum. This theoretical survey is a crucial step towards a feasible characterization of the early stage of excited states in solution during photoisomerization processes wherein multiple electronic states might be populated upon the light radiation, leading to a future molecular-level interpretation of time-resolved spectroscopies.

Benchmark computations of nearly degenerate singlet and triplet states of N-heterocyclic chromophores. I. Wavefunction-based methods

In this paper, we investigate the role of electron correlation in predicting the S1-S0 and T1-S0 excitation energies and, hence, the singlet-triplet gap (ΔEST) in a set of cyclazines, which act as templates for potential candidates for fifth generation organic light emitting diode materials. This issue has recently garnered much interest with the focus being on the inversion of the ΔEST, although experiments have indicated near degenerate levels with both positive and negative being within the experimental error bar [J. Am. Chem. Soc. 102, 6068 (1980), J. Am. Chem. Soc. 108, 17(1986)]. We have carried out a systematic and exhaustive study of various excited state electronic structure methodologies and identified the strengths and shortcomings of the various approaches and approximations in view of this challenging case. We have found that near degeneracy can be achieved either with a proper balance of static and dynamic correlation in multireference theories or with state-specific orbital corrections, including its coupling with correlation. The role of spin contamination is also discussed. Eventually, this paper seeks to produce benchmark numbers for establishing cost-effective theories, which can then be used for screening derivatives of these templates with desirable optical and structural properties. Additionally, we would like to point out that the use of domain-based local pair natural orbital-similarity transformed EOM-coupled cluster singles and doubles as the benchmark for ΔEST [as used in J. Phys. Chem. A 126(8), 1378 (2022), Chem. Phys. Lett. 779, 138827 (2021)] is not a suitable benchmark for these classes of molecules.

Two-Photon Absorption Strengths of Small Molecules: Reference CC3 Values and Benchmarks

We present a large dataset of highly accurate two-photon transition strengths (δTPA) determined for standard small molecules. Our reference values have been calculated using the quadratic response implementation of the third-order coupled cluster method including iterative triples (Q-CC3). The aug-cc-pVTZ atomic basis set is used for molecules with up to five non-hydrogen atoms, while larger molecules are assessed with aug-cc-pVDZ; the differences due to the basis sets are discussed. This dataset, encompassing 82 singlet transitions of various characters (Rydberg, valence, and double excitations), enables a comprehensive benchmark of smaller basis sets and alternative wavefunction methods when Q-CC3 calculations become beyond reach as well as time-dependent density functional theory (TD-DFT) approaches. The evaluated wavefunction methods include quadratic response and equation-of-motion CCSD approximations, Q-CC2, and second-order algebraic diagrammatic construction in its intermediate state representation (I-ADC2). In the TD-DFT framework, a set of five commonly used exchange-correlation functionals are evaluted. This extensive analysis provides a quantitative assessment of these methods, revealing how different system sizes, response intensities, and types of transitions affect their performances.

Rank-Reduced Equation-of-Motion Coupled Cluster Triples: an Accurate and Affordable Way of Calculating Electronic Excitation Energies

In the present work, we report an implementation of the rank-reduced equation-of-motion coupled cluster method with approximate triple excitations (RR-EOM-CC3). The proposed variant relies on tensor decomposition techniques in order to alleviate the high cost of computing and manipulating the triply excited amplitudes. In the RR-EOM-CC3 method, both ground-state and excited-state triple-excitation amplitudes are compressed according to the Tucker-3 format. This enables factorization of the working equations such that the formal scaling of the method is reduced to N6, where N is the system size. An additional advantage of our method is the fact that the accuracy can be strictly controlled by proper choice of two parameters defining sizes of triple-excitation subspaces in the Tucker decomposition for the ground and excited states. Optimal strategies of selecting these parameters are discussed. The developed method has been tested in a series of calculations of electronic excitation energies and compared to its canonical EOM-CC3 counterpart. Errors several times smaller than the inherent error of the canonical EOM-CC3 method (in comparison to FCI) are straightforward to achieve. This conclusion holds both for valence states dominated by single excitations and for states with pronounced doubly excited character. Taking advantage of the decreased scaling, we demonstrate substantial computational costs reductions (in comparison with the canonical EOM-CC3) in the case of two large molecules - l-proline and heptazine. This illustrates the usefulness of the RR-EOM-CC3 method for accurate determination of excitation energies of large molecules.

Exploring Ground and Excited States Via Single Reference Coupled-Cluster Theory and Algebraic Geometry

The exploration of the root structure of coupled cluster (CC) equations holds both foundational and practical significance for computational quantum chemistry. This study provides insight into the intricate root structures of these nonlinear equations at both the CCD and CCSD level of theory. We utilize computational techniques from algebraic geometry, specifically the monodromy and parametric homotopy continuation methods, to calculate the full solution set. We compare the computed CC roots against various established theoretical upper bounds, shedding light on the accuracy and efficiency of these bounds. We hereby focus on the dissociation processes of four-electron systems such as (H2)2 in both D2h and D∞h configurations, H4 symmetrically distorted on a circle, and lithium hydride. We moreover investigate the ability of single-reference CC solutions to approximate excited state energies. We find that multiple CC roots describe energies of excited states with high accuracy. Our investigations reveal that for systems like lithium hydride, CC not only provides high-accuracy approximations to several excited state energies but also to the states themselves.

Reference CC3 Excitation Energies for Organic Chromophores: Benchmarking TD-DFT, BSE/GW, and Wave Function Methods

To expand the QUEST database of highly accurate vertical transition energies, we consider a series of large organic chromogens ubiquitous in dye chemistry, such as anthraquinone, azobenzene, BODIPY, and naphthalimide. We compute, at the CC3 level of theory, the singlet and triplet vertical transition energies associated with the low-lying excited states. This leads to a collection of more than 120 new highly accurate excitation energies. For several singlet transitions, we have been able to determine CCSDT transition energies with a compact basis set, finding minimal deviations from the CC3 values for most states. Subsequently, we employ these reference values to benchmark a series of lower-order wave function approaches, including the popular ADC(2) and CC2 schemes, as well as time-dependent density-functional theory (TD-DFT), both with and without applying the Tamm-Dancoff approximation (TDA). At the TD-DFT level, we evaluate a large panel of global, range-separated, local, and double hybrid functionals. Additionally, we assess the performance of the Bethe-Salpeter equation (BSE) formalism relying on both G0W0 and evGW quasiparticle energies evaluated from various starting points. It turns out that CC2 and ADC(2.5) are the most accurate models among those with respective O ( N 5 ) and O ( N 6 ) scalings with system size. In contrast, CCSD does not outperform CC2. The best performing exchange-correlation functionals include BMK, M06-2X, M06-SX, CAM-B3LYP, ωB97X-D, and LH20t, with average deviations of approximately 0.20 eV or slightly below. Errors on vertical excitation energies can be further reduced by considering double hybrids. Both SOS-ωB88PP86 and SOS-ωPBEPP86 exhibit particularly attractive performances with overall quality on par with CC2, whereas PBE0-DH and PBE-QIDH are only slightly less efficient. BSE/evGW calculations based on Kohn-Sham starting points have been found to be particularly effective for singlet transitions, but much less for their triplet counterparts.

Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction

Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. We show that the sCI framework enables accurate prediction of band gaps in LiCoO2, closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t2g and eg orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.

Analytic Gradients for Equation-of-Motion Coupled Cluster with Single, Double, and Perturbative Triple Excitations

Understanding the process of molecular photoexcitation is crucial in various fields, including drug development, materials science, photovoltaics, and more. The electronic vertical excitation energy is a critical property, for example in determining the singlet-triplet gap of chromophores. However, a full understanding of excited-state processes requires additional explorations of the excited-state potential energy surface and electronic properties, which is greatly aided by the availability of analytic energy gradients. Owing to its robust high accuracy over a wide range of chemical problems, equation-of-motion coupled cluster with single and double excitations (EOM-CCSD) is a powerful method for predicting excited-state properties, and the implementation of analytic gradients of many EOM-CCSD variants (excitation energies, ionization potentials, electron attachment energies, etc.) along with numerous successful applications highlights the flexibility of the method. In specific cases where a higher level of accuracy is needed or in more complex electronic structures, the inclusion of triple excitations becomes essential, for example, in the EOM-CCSD* approach of Saeh and Stanton. In this work, we derive and implement for the first time the analytic gradients of EOMEE-CCSD*, which also provides a template for analytic gradients of related excited-state methods with perturbative triple excitations. The capabilities of analytic EOMEE-CCSD* gradients are illustrated by several representative examples.

Compactification of determinant expansions via transcorrelation

Although selected configuration interaction (SCI) algorithms can tackle much larger Hilbert spaces than the conventional full CI method, the scaling of their computational cost with respect to the system size remains inherently exponential. In addition, inaccuracies in describing the correlation hole at small interelectronic distances lead to the slow convergence of the electronic energy relative to the size of the one-electron basis set. To alleviate these effects, we show that the non-Hermitian, transcorrelated (TC) version of SCI significantly compactifies the determinant space, allowing us to reach a given accuracy with a much smaller number of determinants. Furthermore, we note a significant acceleration in the convergence of the TC-SCI energy as the basis set size increases. The extent of this compression and the energy convergence rate are closely linked to the accuracy of the correlation factor used for the similarity transformation of the Coulombic Hamiltonian. Our systematic investigation of small molecular systems in increasingly large basis sets illustrates the magnitude of these effects.

Accurate computation of quantum excited states with neural networks

We present an algorithm to estimate the excited states of a quantum system by variational Monte Carlo, which has no free parameters and requires no orthogonalization of the states, instead transforming the problem into that of finding the ground state of an expanded system. Arbitrary observables can be calculated, including off-diagonal expectations, such as the transition dipole moment. The method works particularly well with neural network ansätze, and by combining this method with the FermiNet and Psiformer ansätze, we can accurately recover excitation energies and oscillator strengths on a range of molecules. We achieve accurate vertical excitation energies on benzene-scale molecules, including challenging double excitations. Beyond the examples presented in this work, we expect that this technique will be of interest for atomic, nuclear, and condensed matter physics.

Extended N-centered ensemble density functional theory of double electronic excitations

A recent work (arXiv:2401.04685) has merged N-centered ensembles of neutral and charged electronic ground states with ensembles of neutral ground and excited states, thus providing a general and in-principle exact (so-called extended N-centered) ensemble density functional theory of neutral and charged electronic excitations. This formalism made it possible to revisit the concept of density-functional derivative discontinuity, in the particular case of single excitations from the highest occupied Kohn-Sham (KS) molecular orbital, without invoking the usual "asymptotic behavior of the density" argument. In this work, we address a broader class of excitations, with a particular focus on double excitations. An exact implementation of the theory is presented for the two-electron Hubbard dimer model. A thorough comparison of the true physical ground- and excited-state electronic structures with that of the fictitious ensemble density-functional KS system is also presented. Depending on the choice of the density-functional ensemble as well as the asymmetry of the dimer and the correlation strength, an inversion of states can be observed. In some other cases, the strong mixture of KS states within the true physical system makes the assignment "single excitation" or "double excitation" irrelevant.

Combined Multireference-Multiscale Approach to the Description of Photosynthetic Reaction Centers

A first-principles description of the primary photochemical processes that drive photosynthesis and sustain life on our planet remains one of the grand challenges of modern science. Recent research established that explicit incorporation of protein electrostatics in excited-state calculations of photosynthetic pigments, achieved for example with quantum-mechanics/molecular-mechanics (QM/MM) approaches, is essential for a meaningful description of the properties and function of pigment-protein complexes. Although time-dependent density functional theory has been used productively so far in QM/MM approaches for the study of such systems, this methodology has limitations. Here we pursue for the first time a QM/MM description of the reaction center in the principal enzyme of oxygenic photosynthesis, Photosystem II, using multireference wave function theory for the high-level QM region. We identify best practices and establish guidelines regarding the rational choice of active space and appropriate state-averaging for the efficient and reliable use of complete active space self-consistent field (CASSCF) and the N-electron valence state perturbation theory (NEVPT2) in the prediction of low-lying excited states of chlorophyll and pheophytin pigments. Given that the Gouterman orbitals are inadequate as a minimal active space, we define specific minimal and extended active spaces for the NEVPT2 description of electronic states that fall within the Q and B bands. Subsequently, we apply our multireference-QM/MM protocol to the description of all pigments in the reaction center of Photosystem II. The calculations reproduce the electrochromic shifts induced by the protein matrix and the ordering of site energies consistent with the identity of the primary donor (ChlD1) and the experimentally known asymmetric and directional electron transfer. The optimized protocol sets the stage for future multireference treatments of multiple pigments, and hence for multireference studies of charge separation, while it is transferable to the study of any photoactive embedded tetrapyrrole system.

Reference Energies for Double Excitations: Improvement and Extension

In the realm of photochemistry, the significance of double excitations (also known as doubly excited states), where two electrons are concurrently elevated to higher energy levels, lies in their involvement in key electronic transitions essential in light-induced chemical reactions as well as their challenging nature from the computational theoretical chemistry point of view. Based on state-of-the-art electronic structure methods (such as high-order coupled-cluster, selected configuration interaction, and multiconfigurational methods), we improve and expand our prior set of accurate reference excitation energies for electronic states exhibiting a substantial amount of double excitations [Loos et al. J. Chem. Theory Comput. 2019, 15, 1939]. This extended collection encompasses 47 electronic transitions across 26 molecular systems that we separate into two distinct subsets: (i) 28 "genuine" doubly excited states where the transitions almost exclusively involve doubly excited configurations and (ii) 19 "partial" doubly excited states which exhibit a more balanced character between singly and doubly excited configurations. For each subset, we assess the performance of high-order coupled-cluster (CC3, CCSDT, CC4, and CCSDTQ) and multiconfigurational methods (CASPT2, CASPT3, PC-NEVPT2, and SC-NEVPT2). Using as a probe the percentage of single excitations involved in a given transition (%T1) computed at the CC3 level, we also propose a simple correction that reduces the errors of CC3 by a factor of 3, for both sets of excitations. We hope that this more complete and diverse compilation of double excitations will help future developments of electronic excited-state methodologies.

Excited State-Specific CASSCF Theory for the Torsion of Ethylene

State-specific complete active space self-consistent field (SS-CASSCF) theory has emerged as a promising route to accurately predict electronically excited energy surfaces away from molecular equilibria. However, its accuracy and practicality for chemical systems of photochemical interest have yet to be fully determined. We investigate the performance of the SS-CASSCF theory for the low-lying ground and excited states in the double bond rotation of ethylene. We show that state-specific approximations with a minimal (2e,2o) active space provide comparable accuracy to state-averaged calculations with much larger active spaces, while optimizing the orbitals for each excited state significantly improves the spatial diffusivity of the wave function. However, the incorrect ordering of state-specific solutions causes excited state solutions to coalesce and disappear, creating unphysical discontinuities in the potential energy surface. Our findings highlight the theoretical challenges that must be overcome to realize practical applications of state-specific electronic structure theory for computational photochemistry.

A DFT/MRCI Hamiltonian parameterized using only ab initio data: I. valence excited states

A new combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian parameterized solely using the benchmark ab initio vertical excitation energies obtained from the QUEST databases is presented. This new formulation differs from all previous versions of the method in that the choice of the underlying exchange-correlation (XC) functional employed to construct the one-particle (orbital) basis is considered, and a new XC functional, QTP17, is chosen for its ability to generate a balanced description of core and valence vertical excitation energies. The ability of the new DFT/MRCI Hamiltonian, termed QE8, to furnish accurate excitation energies is confirmed using benchmark quantum chemistry computations, and a mean absolute error of 0.16 eV is determined for the wide range of electronic excitations included in the validation dataset. In particular, the QE8 Hamiltonian dramatically improves the performance of DFT/MRCI for doubly excited states. The performance of fast approximate DFT/MRCI methods, p-DFT/MRCI and DFT/MRCI(2), is also evaluated using the QE8 Hamiltonian, and they are found to yield excitation energies in quantitative agreement with the parent DFT/MRCI method, with the two methods exhibiting a mean difference of 0.01 eV with respect to DFT/MRCI over the entire benchmark set.

Reference Energies for Valence Ionizations and Satellite Transitions

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.

Frozen Natural Orbitals for the State-Averaged Driven Similarity Renormalization Group

We present a reduced-cost implementation of the state-averaged driven similarity renormalization group (SA-DSRG) based on the frozen natural orbital (FNO) approach. The natural orbitals (NOs) are obtained by diagonalizing the one-body reduced density matrix from SA-DSRG second-order perturbation theory (SA-DSRG-PT2). We consider three criteria to truncate the virtual NOs for the subsequent electron correlation treatment beyond SA-DSRG-PT2. An additive second-order correction is applied to the SA-DSRG Hamiltonian to reintroduce correlation effects from the discarded orbitals. The FNO SA-DSRG method is benchmarked on 35 small organic molecules in the QUEST database. When keeping 98-99% of the cumulative occupation numbers, the mean absolute error in the vertical transition energies due to FNO is less than 0.01 eV. Using the same FNO threshold, we observe a speedup of 9 times compared to the conventional SA-DSRG implementation for nickel carbonyl with a quadruple-ζ basis set. The FNO approach enables nonperturbative SA-DSRG computations on chloroiron corrole [FeCl(C19H11N4)] with more than 1000 basis functions, surpassing the current limit of a conventional implementation.

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.

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear-electronic coupling in addition to accurately predicting Franck-Condon vertical excitations. In this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

Rationale for the extrapolation procedure in selected configuration interaction

Selected configuration interaction (SCI) methods have emerged as state-of-the-art methodologies for achieving high accuracy and generating benchmark reference data for ground and excited states in small molecular systems. However, their precision relies heavily on extrapolation procedures to produce a final estimate of the exact result. Using the structure of the exact electronic energy landscape, we provide a rationale for the common linear extrapolation of the variational energy as a function of the second-order perturbative correction. In particular, we demonstrate that the energy gap and the coupling between the so-called internal and external spaces are the key factors determining the rate at which the linear regime is reached. Starting from the first principles, we also derive a new non-linear extrapolation formula that improves the post-processing of data generated from SCI methods and can be applied to both ground- and excited-state energies.

Perspective on Coupled-cluster Theory. The evolution toward simplicity in quantum chemistry

Coupled-cluster theory has revolutionized quantum chemistry. It has provided the framework to effectively solve the problem of electron correlation, the main focus of the field for over 60 years. This has enabled ab initio quantum chemistry to provide predictive quality results for most quantities of interest that are obtainable from first-principle calculations. The best that one can do in a basis is the 'full CI,' the exact solution of the non-relativistic Schrödinger equation or, if need be, the relativistic Dirac equation. With due regard to converging the basis set and adequate consideration of higher clusters and relativity in a calculation, virtually predictive results can be obtained. But in addition to its numerical performance, coupled-cluster theory also offers a conceptually new, many-body foundation for the theory that should be appreciated by all practitioners. The latter is emphasized in this perspective, leading to the 'evolution toward simplicity' in the title. The ultimate theory will benefit from the several features that are uniquely exact in coupled-cluster theory and its equation-of-motion (EOM-CC) extensions.

Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone

We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.

Performance of point charge embedding schemes for excited states in molecular organic crystals

Modeling excited state processes in molecular crystals is relevant for several applications. A popular approach for studying excited state molecular crystals is to use cluster models embedded in point charges. In this paper, we compare the performance of several embedding models in predicting excited states and S1-S0 optical gaps for a set of crystals from the X23 molecular crystal database. The performance of atomic charges based on ground or excited states was examined for cluster models, Ewald embedding, and self-consistent approaches. We investigated the impact of various factors, such as the level of theory, basis sets, embedding models, and the level of localization of the excitation. We consider different levels of theory, including time-dependent density functional theory and Tamm-Dancoff approximation (TDA) (DFT functionals: ωB97X-D and PBE0), CC2, complete active space self-consistent field, and CASPT2. We also explore the impact of selection of the QM region, charge leakage, and level of theory for the description of different kinds of excited states. We implemented three schemes based on distance thresholds to overcome overpolarization and charge leakage in molecular crystals. Our findings are compared against experimental data, G0W0-BSE, periodic TDA, and optimally tuned screened range-separated functionals.

Cluster-Based Approach Utilizing Optimally Tuned TD-DFT to Calculate Absorption Spectra of Organic Semiconductor Thin Films

The photophysics of organic semiconductor (OSC) thin films or crystals has garnered significant attention in recent years since a comprehensive theoretical understanding of the various processes occurring upon photoexcitation is crucial for assessing the efficiency of OSC materials. To date, research in this area has relied on methods using Frenkel-Holstein Hamiltonians, calculations of the GW-Bethe-Salpeter equation with periodic boundaries, or cluster-based approaches using quantum chemical methods, with each of the three approaches having distinct advantages and disadvantages. In this work, we introduce an optimally tuned, range-separated time-dependent density functional theory approach to accurately reproduce the total and polarization-resolved absorption spectra of pentacene, tetracene, and perylene thin films, all representative OSC materials. Our approach achieves excellent agreement with experimental data (mostly ≤0.1 eV) when combined with the utilization of clusters comprising multiple monomers and a standard polarizable continuum model to simulate the thin-film environment. Our protocol therefore addresses a major drawback of cluster-based approaches and makes them attractive tools for OSC investigations. Its key advantages include its independence from external, system-specific fitting parameters and its straightforward application with well-known quantum chemical program codes. It demonstrates how chemical intuition can help to reduce computational cost and still arrive at chemically meaningful and almost quantitative results.

Performance Tests of the Second-Order Approximate Internally Contracted Multireference Coupled-Cluster Singles and Doubles Method icMRCC2

Benchmark results are presented for the second-order approximation of the internally contracted multireference coupled-cluster method with single and double excitations, icMRCC2 [Köhn, Bargholz, J. Chem. Phys. 2019, 151, 041106], which was designed as a multireference analogue of the single-reference second-order approximate coupled-cluster method CC2 [Christiansen, Koch, Jørgensen, Chem. Phys. Lett. 1995, 243, 409-418]. Vertical excitation energies of various small to medium-sized organic molecules are investigated based on established test sets from the literature. Additionally, the spectroscopic constants of ground and excited states of diatomics and the geometric parameters of excited triatomic molecules were determined and compared to the experimental data. The results show that the method clearly extends the applicability of single-reference CC2, including doubly excited states, and also artifacts of CC2 like too low Rydberg excitations and too weak multiple bonds are eliminated. The method is computationally more demanding than standard multireference second-order perturbation theories but improves significantly in accuracy, as shown by the benchmark results. In addition, it is demonstrated that small active spaces are often sufficient to obtain accurate energies with icMRCC2. Example applications like the automerization of cyclobutadiene, the deactivation pathway of ethylene, and the excited states of an iron complex with a noninnocent nitrosyl ligand demonstrate the potential of icMRCC2 in cases with strong multireference character.

Oscillator strengths and excited-state couplings for double excitations in time-dependent density functional theory

Although useful to extract excitation energies of states of double-excitation character in time-dependent density functional theory that are missing in the adiabatic approximation, the frequency-dependent kernel derived earlier [Maitra et al., J. Chem. Phys. 120, 5932 (2004)] was not designed to yield oscillator strengths. These are required to fully determine linear absorption spectra, and they also impact excited-to-excited-state couplings that appear in dynamics simulations and other quadratic response properties. Here, we derive a modified non-adiabatic kernel that yields both accurate excitation energies and oscillator strengths for these states. We demonstrate its performance on a model two-electron system, the Be atom, and on excited-state transition dipoles in the LiH molecule at stretched bond-lengths, in all cases producing significant improvements over the traditional approximations.

Variational Active Space Selection with Multiconfiguration Pair-Density Functional Theory

The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.

Linearized Pair-Density Functional Theory for Vertical Excitation Energies

Multiconfiguration pair-density functional theory (MC-PDFT) is a computationally efficient method that computes the energies of electronic states in a state specific or state average framework via an on-top functional. However, MC-PDFT does not include state interaction among these states since the final energies do not come from the diagonalization of an effective model-space Hamiltonian. Recently, multistate extensions such as linearized PDFT (L-PDFT) have been developed to accurately model the potentials near conical intersections and avoided crossings. However, there has not been any systematic study evaluating their performance for predicting vertical excitations at the equilibrium geometry of a molecule, when the excited states are generally well separated. In this paper, we report the performance of L-PDFT on the extensive QUESTDB data set of vertical excitations using a database of automatically selected active spaces. We show that L-PDFT performs well on all these excitations and successfully reproduces the performance of MC-PDFT. These results further demonstrate the potential of L-PDFT, as its scaling is constant with the number of states included in the state-average manifold, whereas MC-PDFT scales linearly in this regard.

Reference Vertical Excitation Energies for Transition Metal Compounds

To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.

Earth Mover's Distance as a Metric to Evaluate the Extent of Charge Transfer in Excitations Using Discretized Real-Space Densities

This paper presents a novel theoretical measure, μEMD, based on the earth mover's distance (EMD), for quantifying the density shift caused by electronic excitations in molecules. As input, the EMD metric uses only the discretized ground- and excited-state electron densities in real space, rendering it compatible with almost all electronic structure methods used to calculate excited states. The EMD metric is compared against other popular theoretical metrics for describing the extent of electron-hole separation in a wide range of excited states (valence, Rydberg, charge transfer, etc.). The results showcase the EMD metric's effectiveness across all excitation types and suggest that it is useful as an additional tool to characterize electronic excitations. The study also reveals that μEMD can function as a promising diagnostic tool for predicting the failure of pure exchange-correlation functionals. Specifically, we show statistical relationships among the functional-driven errors, the exact exchange content within the functional, and the magnitude of μEMD values.

Generalization of the Tensor Product Selected CI Method for Molecular Excited States

In a recent paper [JCTC, 2020, 16, 6098], we introduced a new approach for accurately approximating full CI ground states in large electronic active-spaces called Tensor Product Selected CI (TPSCI). In TPSCI, a large orbital active space is first partitioned into disjoint sets (clusters) for which the exact, local many-body eigenstates are obtained. Tensor products of these locally correlated many-body states are taken as the basis for the full, global Hilbert space. By folding correlation into the basis states themselves, the low-energy eigenstates become increasingly sparse, creating a more compact selected CI expansion. While we demonstrated that this approach can improve accuracy for a variety of systems, there is even greater potential for applications to excited states, particularly those which have some excited-state character. In this paper, we report on the accuracy of TPSCI for excited states, including a far more efficient implementation in the Julia programming language. In traditional SCI methods that use a Slater determinant basis, accurate excitation energies are obtained only after a linear extrapolation and at a large computational cost. We find that TPSCI with perturbative corrections provides accurate excitation energies for several excited states of various polycyclic aromatic hydrocarbons with respect to the extrapolated result (i.e., near exact result). Further, we use TPSCI to report highly accurate estimates of the lowest 31 eigenstates for a tetracene tetramer system with an active space of 40 electrons in 40 orbitals, giving direct access to the initial bright states and the resulting 18 doubly excited (biexcitonic) states.

Excited states with pair coupled cluster doubles tailored coupled cluster theory

It is known that some non-dynamic effects of electron correlation can be included in coupled cluster theory using a tailoring technique that separates the effects of non-dynamic and dynamic correlations. Recently, the simple pCCD (pair coupled cluster doubles) wavefunction was shown to provide good results for some non-dynamic correlation problems, such as bond-breaking, in a spin-adapted way with no active space selection. In this paper, we report a study of excited states using "tailored coupled cluster singles and doubles," to attempt to use pCCD as a kernel for more complete coupled-cluster singles and doubles (CCSD) results for excited states. Several excited states are explored from those primarily due to single excitations to those dominated by doubly excited states and from singlet-triplet splittings for some diradical states. For the first two situations, tailored pCCD-TCCSD offers no improvement over equation of motion-CCSD. However, when we explore the singlet-triplet gap of diradical molecules that are manifestly multi-reference, a pCCD kernel provides improved results, particularly with generalized valence bond orbitals.

Deciphering Methylation Effects on S2(ππ*) Internal Conversion in the Simplest Linear α,β-Unsaturated Carbonyl

Chemical substituents can influence photodynamics by altering the location of critical points and the topography of the potential energy surfaces (electronic effect) and by selectively modifying the inertia of specific nuclear modes (inertial effects). Using nonadiabatic dynamics simulations, we investigate the impact of methylation on S₂(ππ*) internal conversion in acrolein, the simplest linear α,β-unsaturated carbonyl. Consistent with time constants reported in a previous time-resolved photoelectron spectroscopy study, S₂ → S₁ deactivation occurs on an ultrafast time scale (∼50 fs). However, our simulations do not corroborate the sequential decay model used to fit the experiment. Instead, upon reaching the S₁ state, the wavepacket bifurcates: a portion undergoes ballistic S₁ → S₀ deactivation (∼90 fs) mediated by fast bond-length alternation motion, while the remaining decays on the picosecond time scale. Our analysis reveals that methyl substitution, generally assumed to mainly exert inertial influence, is also manifested in important electronic effects due to its weak electron-donating ability. While methylation at the β C atom gives rise to effects principally of an inertial nature, such as retarding the twisting motion of the terminal -CHCH₃ group and increasing its coupling with pyramidalization, methylation at the α or carbonyl C atom modifies the potential energy surfaces in a way that also contributes to altering the late S₁-decay behavior. Specifically, our results suggest that the observed slowing of the picosecond component upon α-methylation is a consequence of a tighter surface and reduced amplitude along the central pyramidalization, effectively restricting the access to the S₁/S₀-intersection seam. Our work offers new insight into the S₂(ππ*) internal conversion mechanisms in acrolein and its methylated derivatives and highlights site-selective methylation as a tuning knob to manipulate photochemical reactions.

Classification of doubly excited molecular electronic states

Electronic states with partial or complete doubly excited character play a crucial role in many areas, such as singlet fission and non-linear optical spectroscopy. Although doubly excited states have been studied in polyenes and related systems for many years, the assignment as singly vs. doubly excited, even in the simplest case of butadiene, has sparked controversies. So far, no well-defined framework for classifying doubly excited states has been developed, and even more, there is not even a well-accepted definition of doubly excited character as such. Here, we present a solution: a physically motivated definition of doubly excited character based on operator expectation values and density matrices, which works independently of the underlying orbital representation, avoiding ambiguities that have plagued earlier studies. Furthermore, we propose a classification scheme to differentiate three cases: (i) two single excitations occurring within two independent pairs of orbitals leaving four open shells (DOS), (ii) the promotion of both electrons to the same orbital, producing a closed-shell determinant (DCS), and (iii) a mixture of singly and doubly excited configurations not aligning with either one of the previous cases (Dmix). We highlight their differences in underlying energy terms and explain their signatures in practical computations. The three cases are illustrated through various high-level computational methods using dimers for DOS, polyenes for Dmix, and cyclobutane and tetrazine for DCS. The conversion between DOS and DCS is investigated using a well-known photochemical reaction, the photodimerization of ethylene. This work provides a deeper understanding of doubly excited states and may guide more rigorous discussions toward improving their computational description while also giving insight into their fundamental photophysics.

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.

Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States

A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.

Impact of Di- and Poly-Radical Characters on the Relative Energy of the Doubly Excited and La States of Linear Acenes and Cyclacenes

Linear and cyclic acenes are polycyclic aromatic hydrocarbons that can be viewed as building blocks of graphene nanoribbons and carbon nanotubes, respectively. While short linear acenes demonstrated remarkable efficiency in several optoelectronic applications, the longer members are unstable and difficult to synthesize as their cyclic counterparts. Recent progress in on-surface synthesis, a powerful tool to prepare highly reactive species, opens promising perspectives and motivates the computational investigations of these potentially functional molecules. Owing to their di- and poly-radical character, low-lying excited states dominated by doubly excited configurations are expected to become more important for longer members of both linear and cyclic molecules. In this work, we investigate the lowest-lying La and the doubly excited (DE) state of linear acenes and cyclacenes, with different computational approaches, to assess the influence of the di-/poly-radical characters (increasing with the molecular dimensions) on their relative order. We show that DFT/MRCI calculations correctly reproduce the crossing of the two states for longer linear acenes, while TDUDFT calculations fail to predict the correct excitation energy trend of the DE state. The study suggests a similarity in the excited electronic state pattern of long linear and cyclic acenes leading ultimately to a lowest lying dark DE state for both.

R2022: A DFT/MRCI Ansatz with Improved Performance for Double Excitations

A reformulation of the combined density functional theory and multireference configuration interaction method (DFT/MRCI) is presented. Expressions for ab initio matrix elements are used to derive correction terms for a new effective Hamiltonian. On the example of diatomic carbon, the correction terms are derived, focusing on the doubly excited ¹Δ_g state, which was problematic in previous formulations of the method, as were double excitations in general. The derivation shows that a splitting of the parameters for intra- and interorbital interactions is necessary for a concise description of the underlying physics. Results for ¹L_a and ¹L_b states in polyacenes and ¹A_u and ¹A_g states in mini-β-carotenoids suggest that the presented formulation is superior to former effective Hamiltonians. Furthermore, statistical analysis reveals that all the benefits of the previous DFT/MRCI Hamiltonians are retained. Consequently, the here presented formulation should be considered as the new standard for DFT/MRCI calculations.

Excited states of ortho-nitrobenzaldehyde as a challenging case for single- and multi-reference electronic structure theory

We present a large set of vertical excitation calculations for the ortho-nitrobenzaldehyde (oNBA) molecule, which exhibits a very challenging excited-state electronic structure like other nitroaromatic compounds. The single-reference methods produce mostly consistent results up to about 5.5 eV. By contrast, the CAS second-order perturbation theory (CASPT2) results depend sensitively on the employed parameters. At the CAS self-consistent field level, the energies of the bright ππ * $$ {\pi \pi}^{\ast } $$ states are strongly overestimated while doubly excited states appear too low and mix with these ππ * $$ {\pi \pi}^{\ast } $$ states. This mixing hampers the CASPT2 step, leading to inconsistent results. Only by increasing the number of states in the state-averaging step to about 40-to cover all bright ππ * $$ {\pi \pi}^{\ast } $$ states embedded in the double excitations-and employing extended multistate CASPT2 could CASPT2 results consistent with experiment be obtained. We assign the four bands in the molecule's spectrum: The weakest band at 3.7 eV arises from the n NO 2 π * $$ {n}_{\mathrm{NO}2}{\pi}^{\ast } $$ states, the second one at 4.4 eV from the ππ * $$ {\pi \pi}^{\ast } $$ ( L b $$ {L}_b $$ ) state, the shoulder at 5.2 eV from the ππ * $$ {\pi \pi}^{\ast } $$ ( L a $$ {L}_a $$ ) state, and the maximum at 5.7 eV from the ππ * $$ {\pi \pi}^{\ast } $$ ( B a / B b $$ {B}_a/{B}_b $$ ) states. We also highlight the importance of modern wave function analysis techniques in elucidating the absorption spectrum of challenging molecules.

Excitation energies of polycylic aromatic hydrocarbons by double-hybrid functionals: Assessing the PBE0-DH and PBE-QIDH models and their range-separated versions

A family of non-empirical double-hybrid (DH) density functionals, such as Perdew-Burke-Ernzerhof (PBE)0-DH, PBE-QIDH, and their range-separated exchange (RSX) versions RSX-0DH and RSX-QIDH, all using Perdew-Burke-Ernzerhof(PBE) exchange and correlationfunctionals, is applied here to calculate the excitation energies for increasingly longer linear and cyclic acenes as part of their intense benchmarking for excited states of all types. The energies for the two lowest-lying singlet ¹L_a and ¹L_b states of linear oligoacenes as well as the triplet ³L_a and ³L_b states, are calculated and compared with experimental results. These functionals clearly outperform the results obtained from hybrid functionals and favorably compare with other double-hybrid expressions also tested here, such as B2-PLYP, B2GP-PLYP, ωB2-PLYP, and ωB2GP-PLYP. The study is complemented by the computation of adiabatic S₀-T₁ singlet-triplet energy difference for linear acenes as well as the extension of the study to strained cyclic oligomers, showing how the family of non-empirical expressions robustly leads to competitive results.

Dark Doubly Excited States with Modified Coupled Cluster Models: A Reliable Compromise between Cost and Accuracy?

To treat doubly excited states, the treatment of triple excitations is considered necessary in the framework of equation-of-motion coupled cluster (EOM-CC) methods. We investigate models without explicit triples and seek quantitative measure for the performance of EOM based on CC with singles and doubles (CCSD) or modified CCSD (Distinguishable Cluster Approximation) approaches for states with predominant double excitation character. We also test the efficacy of including triples in perturbative manner through EOM-CCSD(T) and in an iterative way through EOM-CCSDT-3 method. Extended similarity transformed EOM-CCSD(EXT-STEOM-CCSD) method is also tested and provides superior quality results at comparatively low cost. We use the QUEST2 benchmark set of double excitations proposed by Loos et al. [ J. Chem. Theory Comput.2019, 15, 1939] to investigate the performance of methods such as EOM-CCSD, EOM-DCSD, EXT-STEOM-CCSD, ΔCCSD, and ΔDCSD. We also test a tailored CC approach, ΔpairCCD-TCCSD.

Assessment of State-Averaged Driven Similarity Renormalization Group on Vertical Excitation Energies: Optimal Flow Parameters and Applications to Nucleobases

We present a comprehensive excited-state benchmark for the state-averaged (SA) driven similarity renormalization group (DSRG) [Li, C.; Evangelista, F. A. J. Chem. Phys. 2018, 148, 124106]. Following the QUEST database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; Loos, P.-F. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2021, 11, e1517], 280 vertical transition energies of 35 medium-sized molecules are computed using the SA-DSRG derived second- and third-order perturbation theories (PT2/PT3) along with a nonperturbative approach [sq-LDSRG(2)]. Comparing to the theoretical best estimates, the optimal flow parameter is found to be 0.35 and 2.0 Eh-2 for SA-DSRG-PT2 and SA-DSRG-PT3, respectively. For SA-sq-LDSRG(2), a flow parameter of 1.5 Eh-2 provides converged equations without compromising the accuracy. We then assess the accuracy of the SA-DSRG hierarchy using these parameters. The SA-DSRG-PT2 scheme outperforms the level-shifted CASPT2 by 0.10 eV in mean absolute error (MAE), yet this accuracy is slightly inferior than that of CASPT2 with the ionization-potential-electron-affinity shift. Both SA-DSRG-PT3 and SA-sq-LDSRG(2) yield a MAE of 0.10 eV, which is comparable to that of CASPT3 (0.09 eV). Finally, we compute vertical excitation energies of several low-lying singlet states of nucleobases. The SA-sq-LDSRG(2) approach provides highly accurate results for π → π* excitations, while n → π* transitions are better described by SA-DSRG-PT3.

Ground- and Excited-State Dipole Moments and Oscillator Strengths of Full Configuration Interaction Quality

We report ground- and excited-state dipole moments and oscillator strengths (computed in different "gauges" or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as Configuration Interaction using a Perturbative Selection made Iteratively (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CC with singles and doubles (CCSD) to CC with singles, doubles, and triples (CCSDT) or from CCSDT to CC with singles, doubles, triples, and quadruples (CCSDTQ)) reduces the average error with respect to the near-FCI reference values by approximately 1 order of magnitude.