Excited states via coupled cluster theory without equation-of-motion methods: Seeking higher roots with application to doubly excited states and double core hole states

Quantum Equation-of-Motion Method with Single, Double, and Triple Excitations

The quantum equation-of-motion (qEOM) method with single and double excitations (qEOM-SD) has been proposed to study electronically excited states, but it fails to handle states dominated by double excitations. In this work, we reformulate the qEOM method within the effective Hamiltonian framework that satisfies the killer condition, and then present an efficient implementation incorporating single, double, and triple excitations. To reduce computational complexity, we employ point-group symmetry and perturbation theory to screen triple excitations, effectively reducing the scaling from No6Nv6 to No5Nv5, where No and Nv are the numbers of occupied and virtual spin orbitals, respectively. Furthermore, we account for the effect of neglected triple excitations by introducing a perturbative correction to the excitation energy. We apply this method to challenging cases where the qEOM-SD method exhibits significant errors, such as the 2 1Δ state of CH+ and the 2 1Σ state of HF. Our new method achieves energy errors below 0.18 eV while incorporating less than 8.2% of triple excitations. Additionally, we extend the operator screening technique to the quantum subspace expansion method for the efficient inclusion of selected triple excitations.

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.

Singlet-Triplet Inversion in Triangular Boron Carbon Nitrides

The discovery of singlet-triplet (ST) inversion in some π-conjugated triangle-shaped boron carbon nitrides is a remarkable breakthrough that defies Hund's first rule. Deeply rooted in strong electron-electron interactions, ST inversion has garnered significant interest due to its potential to revolutionize triplet harvesting in organic LEDs. Using the well-established Pariser-Parr-Pople model for correlated electrons in π-conjugated systems, we employ a combination of CISDT and restricted active space configuration interaction calculations to investigate the photophysics of several triangular boron carbon nitrides. Our findings reveal that ST inversion in these systems is primarily driven by a network of alternating electron-donor and electron-acceptor groups in the molecular rim, rather than by the triangular molecular structure itself.

Elucidating Anomalous Intensity Ratios in Chlorine L-Edge X-ray Absorption Spectroscopy: Multiplet Effects and Core Rydberg Transitions

A relativistic core-valence-separated equation-of-motion coupled cluster (CVS-EOM-CC) study of chlorine L2,3-edge X-ray absorption near-edge structure (XANES) spectra using CH3Cl and CH2ICl as representative molecules is reported. The nearly identical intensity for the main features in the L2- and L3-edge XANES spectra is attributed to multiplet effects and the overlap between core-valence and core Rydberg transitions. The multiplet effects originating from the interaction between the core hole and the C-Cl σ* orbitals account for around half of the deviation of the L3 and L2 intensity ratio from the 2:1 ratio of the numbers of 2p3/2 and 2p1/2 electrons. The 2p3/2 → 4s core Rydberg transitions are shown to overlap with the 2p1/2 → σ* transitions and contribute to the other half of the intensity anomaly. We demonstrate that triple excitations in CVS-EOM-CC calculations play important roles in accurate simulation of the overlap between the 2p1/2 → σ* and 2p3/2 → 4s transitions.

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core-electron binding energy (CEBE) analysis. We extend the development of real-space pseudopotential methods to calculating 1s, 2s, and 2p3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn-Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

Toward Accurate Calculation of Excitation Energies on Quantum Computers with ΔADAPT-VQE: A Case Study of BODIPY Derivatives

Quantum chemistry simulations offer a cost-effective way to computationally design BODIPY photosensitizers. However, accurate predictions of excitation energies pose a challenge for time-dependent density functional theory and equation-of-motion coupled-cluster singles and doubles methods. By contrast, reliable predictions can be achieved by multireference quantum chemistry methods; unfortunately, their computational cost increases exponentially with the number of electrons. Alternatively, quantum computing holds potential for an exact simulation of the photophysical properties in a computationally more efficient way. Herein, we introduce the state-specific ΔUCCSD-VQE (unitary coupled-cluster singles and doubles-variational quantum eigensolver) and ΔADAPT-VQE methods in which the electronically excited state is calculated via a non-Aufbau configuration. We show for six BODIPY derivatives that the proposed methods predict accurate excitation energies that are in good agreement with those from experiments. Due to its performance and simplicity, we believe that ΔADAPT will become a useful approach for the simulation of BODIPY photosensitizers on near-term quantum devices.

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.

Electronic Couplings for Triplet-Triplet Annihilation Upconversion in Crystal Rubrene

Triplet-triplet annihilation photon upconversion (TTA-UC) is a process able to repackage two low-frequency photons into light of higher energy. This transformation is typically orchestrated by the electronic degrees of freedom within organic compounds possessing suitable singlet and triplet energies and electronic couplings. In this work, we propose a computational protocol for the assessment of electronic couplings crucial to TTA-UC in molecular materials and apply it to the study of crystal rubrene. Our methodology integrates sophisticated yet computationally affordable approaches to quantify couplings in singlet and triplet energy transfer, the binding of triplet pairs, and the fusion to the singlet exciton. Of particular significance is the role played by charge-transfer states along the b-axis of rubrene crystal, acting as both partial quenchers of singlet energy transfer and mediators of triplet fusion. Our calculations identify the π-stacking direction as holding notable triplet energy transfer couplings, consistent with the experimentally observed anisotropic exciton diffusion. Finally, we have characterized the impact of thermally induced structural distortions, revealing their key role in the viability of triplet fusion and singlet fission. We posit that our approaches are transferable to a broad spectrum of organic molecular materials, offering a feasible means to quantify electronic couplings.

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.

Aufbau Suppressed Coupled Cluster Theory for Electronically Excited States

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

Two determinant distinguishable cluster

A two reference determinant version of the distinguishable cluster with singles and doubles (DCSD) has been developed. We have implemented the two determinant distinguishable cluster (2D-DCSD) and the corresponding traditional 2D-CCSD method in a new open-source package written in Julia called ElemCo.jl. The methods were benchmarked on singlet and triplet excited states of valence and Rydberg character, as well as for singlet-triplet gaps of diradicals. It is demonstrated that the distinguishable cluster approximation improves the accuracy of 2D-CCSD.

Concurrence distribution in excited states of the one-dimensional spin-1/2 transverse-field XY model: Two different regions

We investigate the variation of concurrence in a spin-1/2 transverse field XY chain system in an excited state. Initially, we precisely solve the eigenvalue problem of the system Hamiltonian using the fermionization technique. Subsequently, we calculate the concurrence between nearest-neighbor pairs of spins in all excited states with higher energy than the ground state. Below the factorized field, denoted as h_{f}=sqrt[J^{2}-(Jδ)^{2}], we find no pairwise entanglement between nearest neighbors in excited states. At the factorized field, corresponding to a factorized state, we observe weak concurrence in very low energy states. Beyond h_{f}, the concurrence strengthens, entangling all excited states. The density of entangled states peaks at the center of the excited spectrum. Additionally, the distribution of concurrence reveals that the midpoint of the nonzero concurrence range harbors the most entangled excited states.

Seniority and Hierarchy Configuration Interaction for Radicals and Excited States

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

Unraveling the variational breakdown of core valence separation calculations: Diagnostic and cure to the over relaxation error of double core hole states

The core valence separation (CVS) approximation is the most employed strategy to prevent the variational collapse of standard wave function optimization when attempting to compute electronic states bearing one or more electronic vacancies in core orbitals. Here, we explore the spurious consequences of this approximation on the properties of the computed core hole states. We especially focus on the less studied case of double core hole (DCH) states, whose spectroscopic interest has recently been rapidly growing. We show that the CVS error leads to a systematic underestimation of DCH energies, a property in stark contrast with the case of single core hole states. We highlight that the CVS error can then be interpreted as an over relaxation effect and design a new correction strategy adapted to these specificities.

Characterizing the origin band spectrum of isoquinoline with resonance enhanced multiphoton ionization and electronic structure calculations

We report the experimental resonance enhanced multiphoton ionization spectrum of isoquinoline between 315 and 310 nm, along with correlated electronic structure calculations on the ground and excited states of this species. This spectral region spans the origin transitions to a π-π* excited state, which previous work has suggested to be vibronically coupled with a lower lying singlet n-π* state. Our computational results corroborate previous density functional theory calculations that predict the vertical excitation energy for the n-π* state to be higher than the π-π* state; however, we find an increase in the C-N-C angle brings the n-π* state below the energy of the π-π* state. The calculations find two out-of-plane vibrational modes of the n-π* state, which may be brought into near resonance with the π-π* state as the C-N-C bond angle increases. Therefore, the C-N-C bond angle may be important in activating vibronic coupling between the states. We fit the experimental rotational contour with a genetic algorithm to determine the excited state rotational constants and orientation of the transition dipole moment. The fits show a mostly in-plane polarized transition, and the projection of the transition dipole moment in the a-b plane is about 84° away from the a axis. These results are consistent with the prediction of our electronic structure calculations for the transition dipole moment of the π-π* excited state.

Excited-State-Specific Pseudoprojected Coupled-Cluster Theory

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

The inverted singlet-triplet gap: a vanishing myth?

Molecules with an inverted singlet-triplet gap (STG) between the first excited singlet and triplet states, for example, heptazine, have recently been reported and gained substantial attention since they violate the famous Hund's rule. Utilizing state-of-the-art high-level ab initio methods, the singlet-triplet gap vanishes and approaches zero from below whatever is improved in the theoretical description of the molecules: the basis set or the level of electron correlation. Seemingly, the phenomenon of inverted singlet-triplet gaps tends to vanish the closer we observe.

Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level

The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.

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.

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.

Connections between many-body perturbation and coupled-cluster theories

Here, we build on the works of Scuseria et al. [J. Chem. Phys. 129, 231101 (2008)] and Berkelbach [J. Chem. Phys. 149, 041103 (2018)] to show connections between the Bethe-Salpeter equation (BSE) formalism combined with the GW approximation from many-body perturbation theory and coupled-cluster (CC) theory at the ground- and excited-state levels. In particular, we show how to recast the GW and Bethe-Salpeter equations as non-linear CC-like equations. Similitudes between BSE@GW and the similarity-transformed equation-of-motion CC method are also put forward. The present work allows us to easily transfer key developments and the general knowledge gathered in CC theory to many-body perturbation theory. In particular, it may provide a path for the computation of ground- and excited-state properties (such as nuclear gradients) within the GW and BSE frameworks.

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core-valence separated (FC-CVS) EOM-CCSD approach - a successful variant of EOM-CCSD tailored towards core excitations.

Simulating core electron binding energies of halogenated species adsorbed on ice surfaces and in solution via relativistic quantum embedding calculations

In this work, we investigate the effects of the environment on the X-ray photoelectron spectra of hydrogen chloride and chloride ions adsorbed on ice surfaces, as well as of chloride ions in water droplets. In our approach, we combine a density functional theory (DFT) description of the ice surface with that of halogen species using the recently developed relativistic core-valence separation equation of motion coupled cluster (CVS-EOM-IP-CCSD) via the frozen density embedding formalism (FDE), to determine the K and L_1,2,3 edges of chlorine. Our calculations, which incorporate temperature effects through snapshots from classical molecular dynamics simulations, are shown to reproduce experimental trends in the change of the core binding energies of Cl^- upon moving from a liquid (water droplets) to an interfacial (ice quasi-liquid layer) environment. Our simulations yield water valence band binding energies in good agreement with experiment, which vary little between the droplets and the ice surface. For halide core binding energies there is an overall trend for overestimating experimental values, though good agreement between theory and experiment is found for Cl^- in water droplets and on ice. For HCl on the other hand there are significant discrepancies between experimental and calculated core binding energies when we consider structural models that maintain the H-Cl bond more or less intact. An analysis of models that allow for pre-dissociated and dissociated structures suggests that experimentally observed chemical shifts in binding energies between Cl^- and HCl would require that H⁺ (in the form of H₃O⁺) and Cl^- are separated by roughly 4-6 Å.

The intermediate state approach for doubly excited dark states in EOM-coupled-cluster theory

Solution of dark, doubly excited states using equation-of-motion coupled-cluster (EOM-CC) usually equires at least triple excitations or even quadruples beyond the standard singles and doubles (EOM-CCSD) for an appropriate treatment. A new route to obtain these doubly excited states using EOM-CCSD is demonstrated. Traditionally, EOM-CC is performed on a closed shell reference state that has a well-described single reference CC wavefunction. In this Communication, we attempt to use low spin open-shell states such as the MS = 0 triplet and open-shell singlet as a reference state. Using this intermediate excited state as a reference state provides us with the benefit of obtaining a doubly excited state, as a single excitation at the cost of EOM-CCSD.

Benchmark relativistic delta-coupled-cluster calculations of K-edge core-ionization energies of third-row elements

A benchmark computational study of K-edge core-ionization energies of third-row elements using relativistic delta-coupled-cluster (ΔCC) methods and a revised core-valence separation (CVS) scheme is reported. High-level relativistic (HLR) corrections beyond the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), including the contributions from two-electron picture-change effects, spin-orbit coupling, the Breit term, and quantum electrodynamics effects, have been taken into account and demonstrated to play an important role. Relativistic ΔCC calculations are shown to provide accurate results for core-ionization energies of third-row elements. The SFX2C-1e-CVS-ΔCC results augmented with HLR corrections show a maximum deviation of less than 0.5 eV with respect to experimental values.

Core-hole delocalization for modeling X-ray spectroscopies: A cautionary tale

The influence of core-hole delocalization for X-ray photoelectron, X-ray absorption, and X-ray emission spectrum calculations is investigated in detail, using approaches including response theory, transition-potential methods, and ground state schemes. The question of a localized/delocalized vacancy is relevant for systems with symmetrically equivalent atoms, as well as near-degeneracies which can distribute the core-orbitals over several atoms. We show that issues relating to core-hole delocalization are present for calculations considering explicit core-hole states, e.g. when using a core-excited or core-ionized reference state, or for fractional occupation numbers. Including electron correlation eventually alleviates the issues, but even using CCSD(T) there is a noticable discrepancy between core-ionization energies obtained with a localized and delocalized core-hole (0.5 eV for the carbon K-edge). Within density functional theory, the discrepancy correlates to the exchange interaction involving the core orbitals of the same spin symmetry as the delocalized core-hole. The use of a localized core-hole allows for a reasonably good inclusion of relaxation at lower level of theory, whereas the proper symmetry solution involving a delocalized core-hole requires higher levels of theory to account for the correlation effects involved in orbital relaxation. For linear response methods, we further show that if X-ray absorption spectra are modelled by considering symmetry-unique sets of atoms, care has to be taken such that there are no delocalizations of the core orbitals, which would otherwise introduce shifts in absolute energies and relative features.

Using core-hole reference states for calculating X-ray photoelectron and emission spectra

For the calculation of core-ionization energies (IEs), X-ray photoelectron spectra (XPS), and X-ray emission spectra (XES), a commonly applied approach is to use non-Aufbau reference states with a core-hole as either final (IE and XPS) or initial (XES) state. However, such reference states can introduce numerical instabilities in post-HF methods, relating to the denominator of the energy corrections involved. This may become arbitrarily close to zero if a negative virtual orbital is present, e.g. a core-hole, leading to near-singularities. The resulting instabilities lead to severe convergence issues of the calculation schemes and, in addition, can strongly affect both energies and intensities, with oscillator strengths seen to reach values up to 4 × 10⁷. For the K-edge we propose freezing the highest-energy virtual orbitals which contribute to any denominator below a threshold of 0.1 Hartree. Stable and reliable spectra are then produced, with minimal influence due to freezing energetically high-lying virtual orbitals (typically removing <5% of the total number of MOs). The developed protocol is here tested for Møller-Plesset perturbation theory and for the algebraic diagrammatic construction scheme for the polarization propagator, and it is also relevant for coupled cluster theory and other related methods.

Energy Landscape of State-Specific Electronic Structure Theory

State-specific approximations can provide a more accurate representation of challenging electronic excitations by enabling relaxation of the electron density. While state-specific wave functions are known to be local minima or saddle points of the approximate energy, the global structure of the exact electronic energy remains largely unexplored. In this contribution, a geometric perspective on the exact electronic energy landscape is introduced. On the exact energy landscape, ground and excited states form stationary points constrained to the surface of a hypersphere, and the corresponding Hessian index increases at each excitation level. The connectivity between exact stationary points is investigated, and the square-magnitude of the exact energy gradient is shown to be directly proportional to the Hamiltonian variance. The minimal basis Hartree-Fock and excited-state mean-field representations of singlet H2 (STO-3G) are then used to explore how the exact energy landscape controls the existence and properties of state-specific approximations. In particular, approximate excited states correspond to constrained stationary points on the exact energy landscape, and their Hessian index also increases for higher energies. Finally, the properties of the exact energy are used to derive the structure of the variance optimization landscape and elucidate the challenges faced by variance optimization algorithms, including the presence of unphysical saddle points or maxima of the variance.

Real-Time Equation-of-Motion CCSD Cumulant Green's Function

Many-body excitations in X-ray photoemission spectra have been difficult to simulate from first principles. We have recently developed a cumulant-based one-electron Green's function method using the real-time coupled-cluster-singles equation-of-motion approach (RT-EOM-CCS) that provides a general framework for treating these problems. Here we extend this approach to include double excitations in the ground-state energy and in the coupled cluster amplitudes, which have been implemented using subroutines generated by the Tensor Contraction Engine (TCE). As in the case of the singles approximation, RT-EOM-CCSD yields a nonperturbative cumulant form of the Green's function in terms of the time-dependent cluster amplitudes, adding nonlinear corrections to the traditional cumulant forms. The extended approach is applied to the core-hole spectral function for small molecular systems. We find that, when core-optimized basis sets are used, the doubles contributions reduce the mean absolute errors in the core binding energies of the 10e systems from 0.8 to 0.3 eV. They also significantly improve the quasiparticle-satellite gap by reducing its overestimation from about 3-5 to about 0-1 eV in CH₄, NH₃, and H₂O, and also improving the overall shape of the satellite features. Finally, we demonstrate the application of the new implementation to the larger, classical XPS ESCA series of molecules and show that the singles approximation can be paired with a modest basis set to study carbon speciation.

Specific chemical bond relaxation unraveled by analysis of shake-up satellites in the oxygen single site double core hole spectrum of CO2

We recently developed [A. Ferté, et al., J. Phys. Chem. Lett., 2020, 11, 4359] a method to compute single site double core hole (ssDCH or K^-2) spectra. We refer to that method as NOTA+CIPSI. In the present paper this method is applied to the O K^-2 spectrum of the CO₂ molecule, and we use this as an example to discuss in detail its convergence properties. Using this approach, theoretical spectra in excellent agreement with the experimental one are obtained. Thanks to a thorough interpretation of the shake-up states responsible for the main satellite peaks and through comparison with the O K^-2 spectrum of CO, we can highlight the clear signature of the two non-equivalent carbon oxygen bonds in the oxygen ssDCH CO₂ dication.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Bicyclic Systems

Pursuing our efforts to define highly accurate estimates of the relative energies of excited states in organic molecules, we investigate, with coupled-cluster methods including iterative triples (CC3 and CCSDT), the vertical excitation energies of 10 bicyclic molecules (azulene, benzoxadiazole, benzothiadiazole, diketopyrrolopyrrole, furofuran, phthalazine, pyrrolopyrrole, quinoxaline, tetrathiafulvalene, and thienothiophene). In total, we provide aug-cc-pVTZ reference vertical excitation energies for 91 excited states of these relatively large systems. We use these reference values to benchmark various wave function methods, i.e., CIS(D), EOM-MP2, CC2, CCSD, STEOM-CCSD, CCSD(T)(a)*, CCSDR(3), CCSDT-3, ADC(2), ADC(2.5), and ADC(3), as well as some spin-scaled variants of both CC2 and ADC(2). These results are compared to those obtained previously on smaller molecules. It turns out that while the accuracy of some methods is almost unaffected by system size, e.g., CIS(D) and CC3, the performance of others can significantly deteriorate as the systems grow, e.g., EOM-MP2 and CCSD, whereas others, e.g., ADC(2) and CC2, become more accurate for larger derivatives.

Multi-reference approach to the computation of double core-hole spectra

Double Core-Hole (DCH) states of small molecules are assessed with the restricted active space self-consistent field and multi-state restricted active space perturbation theory of second order approximations. To ensure an unbiased description of the relaxation and correlation effects on the DCH states, the neutral ground-state and DCH wave functions are optimized separately, whereas the spectral intensities are computed with a biorthonormalized set of molecular orbitals within the state-interaction approximation. Accurate shake-up satellite binding energies and intensities of double-core-ionized states (K^-2) are obtained for H₂O, N₂, CO, and C₂H_2n (n = 1-3). The results are analyzed in detail and show excellent agreement with recent theoretical and experimental data. The K^-2 shake-up spectra of H₂O and C₂H_2n molecules are here completely characterized for the first time.

Variational coupled cluster for ground and excited states

In single-reference coupled-cluster (CC) methods, one has to solve a set of non-linear polynomial equations in order to determine the so-called amplitudes that are then used to compute the energy and other properties. Although it is of common practice to converge to the (lowest-energy) ground-state solution, it is also possible, thanks to tailored algorithms, to access higher-energy roots of these equations that may or may not correspond to genuine excited states. Here, we explore the structure of the energy landscape of variational CC and we compare it with its (projected) traditional version in the case where the excitation operator is restricted to paired double excitations (pCCD). By investigating two model systems (the symmetric stretching of the linear H₄ molecule and the continuous deformation of the square H₄ molecule into a rectangular arrangement) in the presence of weak and strong correlations, the performance of variational pCCD (VpCCD) and traditional pCCD is gauged against their configuration interaction (CI) equivalent, known as doubly occupied CI, for reference Slater determinants made of ground- or excited-state Hartree-Fock orbitals or state-specific orbitals optimized directly at the VpCCD level. The influence of spatial symmetry breaking is also investigated.

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: a Comparison between Energy- and Variance-Based Approaches

Self-consistent field methods for excited states offer an attractive low-cost route to study not only excitation energies but also properties of excited states. Here, we present the generalization of two self-consistent field methods, the maximum overlap method (MOM) and the σ-SCF method, to calculate excited states in strong magnetic fields and investigate their stability and accuracy in this context. These methods use different strategies to overcome the well-known variational collapse of energy-based optimizations to the lowest solution of a given symmetry. The MOM tackles this problem in the definition of the orbital occupations to constrain the self-consistent field procedure to converge on excited states, while the σ-SCF method is based on the minimization of the variance instead of the energy. To overcome the high computational cost of the variance minimization, we present a new implementation of the σ-SCF method with the resolution of identity approximation, allowing the use of large basis sets, which is an important requirement for calculations in strong magnetic fields. The accuracy of these methods is assessed by comparison with the benchmark literature data for He, H₂, and CH⁺. The results reveal severe limitations of the variance-based scheme, which become more acute in large basis sets. In particular, many states are not accessible using variance optimization. Detailed analysis shows that this is a general feature of variance optimization approaches due to the masking of local minima in the optimization. In contrast, the MOM shows promising performance for computing excited states under these conditions, yielding results consistent with available benchmark data for a diverse range of electronic states.

Accuracy of the distinguishable cluster approximation for triple excitations for open-shell molecules and excited states

The distinguishable cluster approximation for triple excitations has been applied to calculate thermochemical properties and excited states involving closed-shell and open-shell species, such as small molecules, 3d transition metal atoms, ozone, and an iron-porphyrin model. Excitation energies have been computed using the ΔCC approach by directly optimizing the excited states. A fixed-reference technique has been introduced to target selected spin-states for open-shell molecular systems. The distinguishable cluster approximation consistently improves coupled cluster with singles doubles and triples results for absolute and relative energies. For excited states dominated by a single configuration state function, the fixed-reference approach combined with high-level coupled-cluster methods has a comparable accuracy to the corresponding equation-of-motion coupled-cluster methods with a negligible amount of spin contamination.

Excited States from State-Specific Orbital-Optimized Pair Coupled Cluster

The pair coupled cluster doubles (pCCD) method (where the excitation manifold is restricted to electron pairs) has a series of interesting features. Among others, it provides ground-state energies very close to what is obtained with doubly occupied configuration interaction (DOCI), but with a polynomial cost (compared with the exponential cost of the latter). Here, we address whether this similarity holds for excited states by exploring the symmetric dissociation of the linear H4 molecule. When ground-state Hartree-Fock (HF) orbitals are employed, pCCD and DOCI excited-state energies do not match, a feature that is assigned to the poor HF reference. In contrast, by optimizing the orbitals at the pCCD level (oo-pCCD) specifically for each excited state, the discrepancies between pCCD and DOCI decrease by 1 or 2 orders of magnitude. Therefore, the pCCD and DOCI methodologies still provide comparable energies for excited states, but only if suitable, state-specific orbitals are adopted. We also assessed whether a pCCD approach could be used to directly target doubly excited states, without having to resort to the equation-of-motion (EOM) formalism. In our Δoo-pCCD model, excitation energies are extracted from the energy difference between separate oo-pCCD calculations for the ground state and the targeted excited state. For a set comprising the doubly excited states of CH+, BH, nitroxyl, nitrosomethane, and formaldehyde, we found that Δoo-pCCD provides quite accurate excitation energies, with root-mean-square deviations (with respect to full configuration interaction results) lower than those of CC3 and comparable to those of EOM-CCSDT, two methods with a much higher computational cost.

Relativistic EOM-CCSD for Core-Excited and Core-Ionized State Energies Based on the Four-Component Dirac-Coulomb(-Gaunt) Hamiltonian

We report an implementation of the core-valence separation approach to the four-component relativistic Hamiltonian-based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD) for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants and the use of different Hamiltonians based on the exact two-component (X2C) framework on the energies of different core-ionized and -excited states in halogen- (CH₃I, HX, and X^-, X = Cl-At) and xenon-containing (Xe, XeF₂) species. Our results show that the X2C molecular mean-field approach [Sikkema, J.; J. Chem. Phys. 2009, 131, 124116], based on four-component Dirac-Coulomb mean-field calculations (²DC^M), is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference four-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets lead to non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations (²DCG^M) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.

Perturbation theory in the complex plane: exceptional points and where to find them

We explore the non-Hermitian extension of quantum chemistry in the complex plane and its link with perturbation theory. We observe that the physics of a quantum system is intimately connected to the position of complex-valued energy singularities, known as exceptional points. After presenting the fundamental concepts of non-Hermitian quantum chemistry in the complex plane, including the mean-field Hartree-Fock approximation and Rayleigh-Schrödinger perturbation theory, we provide a historical overview of the various research activities that have been performed on the physics of singularities. In particular, we highlight seminal work on the convergence behaviour of perturbative series obtained within Møller-Plesset perturbation theory, and its links with quantum phase transitions. We also discuss several resummation techniques (such as Padé and quadratic approximants) that can improve the overall accuracy of the Møller-Plesset perturbative series in both convergent and divergent cases. Each of these points is illustrated using the Hubbard dimer at half filling, which proves to be a versatile model for understanding the subtlety of analytically-continued perturbation theory in the complex plane.

Accurate Electronic Excitation Energies in Full-Valence Active Space via Bootstrap Embedding

Fragment embedding has been widely used to circumvent the high computational scaling of using accurate electron correlation methods to describe the electronic ground states of molecules and materials. However, similar applications that utilize fragment embedding to treat electronic excited states are comparably less reported in the literature. The challenge here is twofold. First, most fragment embedding methods are most effective when the property of interest is local, but the change of the wave function upon excitation is nonlocal in general. Second, even for local excitations, an accurate estimate of, for example, the excitation energy can still be challenging owing to the need for a balanced treatment of both the ground and the excited states. In this work, we show that bootstrap embedding (BE), a fragment embedding method developed recently by our group, is promising toward describing general electronic excitations. Numerical simulations show that the excitation energies in full-valence active space (FVAS) can be well-estimated by BE to an error of ∼0.05 eV using relatively small fragments, for both local excitations and the excitations of some large dye molecules that exhibit strong charge-transfer characters. We hence anticipate BE to be a promising solution to accurately describing the excited states of large chemical systems.

An assessment of different electronic structure approaches for modeling time-resolved x-ray absorption spectroscopy

We assess the performance of different protocols for simulating excited-state x-ray absorption spectra. We consider three different protocols based on equation-of-motion coupled-cluster singles and doubles, two of them combined with the maximum overlap method. The three protocols differ in the choice of a reference configuration used to compute target states. Maximum-overlap-method time-dependent density functional theory is also considered. The performance of the different approaches is illustrated using uracil, thymine, and acetylacetone as benchmark systems. The results provide guidance for selecting an electronic structure method for modeling time-resolved x-ray absorption spectroscopy.

Toward Quantum Computing for High-Energy Excited States in Molecular Systems: Quantum Phase Estimations of Core-Level States

This paper explores the utility of the quantum phase estimation (QPE) algorithm in calculating high-energy excited states characterized by the promotion of electrons occupying core-level shells. These states have been intensively studied over the last few decades, especially in supporting the experimental effort at light sources. Results obtained with QPE are compared with various high-accuracy many-body techniques developed to describe core-level states. The feasibility of the quantum phase estimator in identifying classes of challenging shake-up states characterized by the presence of higher-order excitation effects is discussed. We also demonstrate the utility of the QPE algorithm in targeting excitations from specific centers in a molecule. Lastly, we discuss how the lowest-order Trotter formula can be applied to reducing the complexity of the ansatz without affecting the error.

Advanced Computation Method for Double Core Hole Spectra: Insight into the Nature of Intense Shake-up Satellites

Double core hole spectroscopy is an ideal framework for investigating photoionization shake-up satellites. Their important intensity in a single site double core hole (ssDCH) spectrum allows the exploration of the subtle mix of relaxation and correlation effects associated with the inherent multielectronic character of the shake-up process. We present a high-accuracy computation method for single photon double core-shell photoelectron spectra that combines a selected configuration interaction procedure with the use of non-orthogonal molecular orbitals to obtain unbiased binding energy and intensity. This strategy leads to the oxygen ssDCH spectrum of the CO molecule that is in excellent agreement with the experimental result. Through a combined wave function and density analysis, we highlight that the most intense shake-up satellites are characterized by an electronic reorganization that opposes the core hole-induced relaxation.

Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.