Communication: The pole structure of the dynamical polarizability tensor in equation-of-motion coupled-cluster theory

Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-Accelerated Computer Architectures

We present the development and implementation of relativistic coupled cluster linear response theory (CC-LR), which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin for the magnetic properties) as well as taking into account the finite lifetime of excited states in the framework of damped response theory. We showcase our implementation, which is capable to offload the computationally intensive tensor contractions characteristic of coupled cluster theory onto graphical processing units, in the calculation of (a) frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a particular emphasis on the calculation of valence absorption cross sections for the I2 molecule; (b) indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H2Se-H2O dimer as a prototypical system containing hydrogen bonds; and (c) optical rotations at the sodium D line for hydrogen peroxide analogues (H2Y2, Y = O, S, Se, Te). Thanks to this implementation, we are able to show the similarities in performance, but often the significant discrepancies, between CC-LR and approximate methods such as density functional theory. Comparing standard CC response theory with the flavor based upon the equation of motion formalism, we find that for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region, such as spin-spin couplings, on the other hand, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table, it may become increasingly difficult to measure pure optical rotation at the sodium D line due to the appearance of absorbing states.

Theory, implementation, and disappointing results for two-photon absorption cross sections within the doubly electron-attached equation-of-motion coupled-cluster framework

The equation-of-motion coupled-cluster singles and doubles method with double electron attachment (EOM-DEA-CCSD) is capable of computing reliable energies, wave functions, and first-order properties of excited states in diradicals and polyenes that have a significant doubly excited character with respect to the ground state, without the need for including the computationally expensive triple excitations. Here, we extend the capabilities of the EOM-DEA-CCSD method to the calculations of a multiphoton property, two-photon absorption (2PA) cross sections. Closed-form expressions for the 2PA cross sections are derived within the expectation-value approach using response wave functions. We analyze the performance of this new implementation by comparing the EOM-DEA-CCSD energies and 2PA cross sections with those computed using the CC3 quadratic response theory approach. As benchmark systems, we consider transitions to the states with doubly excited character in twisted ethene and in polyenes, for which EOM-EE-CCSD (EOM-CCSD for excitation energies) performs poorly. The EOM-DEA-CCSD 2PA cross sections are comparable with the CC3 results for twisted ethene; however, the discrepancies between the two methods are large for hexatriene. The observed trends are explained by configurational analysis of the 2PA channels.

Cherry-Picking Resolvents: Recovering the Valence Contribution in X-ray Two-Photon Absorption within the Core-Valence-Separated Equation-of-Motion Coupled-Cluster Response Theory

Calculations of first-order response wave functions in the X-ray regime often diverge within correlated frameworks such as equation-of-motion coupled-cluster singles and doubles (EOM-CCSD), a consequence of the coupling with the valence ionization continuum. Here, we extend our strategy of introducing a hierarchy of approximations to the EOM-EE-CCSD resolvent (or, inversely, the model Hamiltonian) involved in the response equations for the calculation of X-ray two-photon absorption (X2PA) cross sections. We exploit the frozen-core core-valence separation (fc-CVS) scheme to first decouple the core and valence Fock spaces, followed by a separate approximate treatment of the valence resolvent. We demonstrate the robust convergence of X-ray response calculations within this framework and compare X2PA spectra of small benchmark molecules with the previously reported density functional theory results.

Probing Molecular Chirality of Ground and Electronically Excited States in the UV-vis and X-ray Regimes: An EOM-CCSD Study

We present several strategies for computing electronic circular dichroism (CD) spectra across different frequency ranges at the equation-of-motion coupled-cluster singles and doubles level of theory. CD spectra of both ground and electronically excited states are discussed. For selected cases, the approach is compared with coupled-cluster linear response results as well as time-dependent density functional theory. The extension of the theory to include the effect of spin-orbit coupling is presented and illustrated by calculations of X-ray CD spectra at the L-edge.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Short Iterative Lanczos Integration in Time-Dependent Equation-of-Motion Coupled-Cluster Theory

A time-dependent (TD) formulation of equation-of-motion coupled-cluster (EOM-CC) theory can provide excited-state information over an arbitrarily wide energy window with a reduced memory footprint relative to conventional, frequency-domain EOM-CC theory. However, the floating-point costs of the time-integration required by TD-EOM-CC are generally far larger than those of the frequency-domain form of the approach. This work considers the potential of the short iterative Lanczos (SIL) integration scheme [J. Chem. Phys. 1986, 85, 5870-5876] to reduce the floating-point costs of TD-EOM-CC simulations. Low-energy and K-edge absorption features for small molecules are evaluated using TD-EOM-CC with single and double excitations, with the time-integrations carried out via SIL and fourth-order Runge-Kutta (RK4) schemes. Spectra derived from SIL- and RK4-driven simulations are nearly indistinguishable, and with an appropriately chosen subspace dimension, the SIL requires far fewer floating-point operations than are required by RK4. For K-edge spectra, SIL is the more efficient scheme by an average factor of 7.2.

The orbital picture of the first dipole hyperpolarizability from many-body response theory

We present an approach for obtaining a molecular orbital picture of the first dipole hyperpolarizability (β) from correlated many-body electronic structure methods. Ab initio calculations of β rely on quadratic response theory, which recasts the sum-over-all-states expression of β into a closed-form expression by calculating a handful of first- and second-order response states; for resonantly enhanced β, damped response theory is used. These response states are then used to construct second-order response reduced one-particle density matrices (1PDMs), which, upon visualization in terms of natural orbitals (NOs), facilitate a rigorous and black-box mapping of the underlying electronic structure with β. We explain the interpretation of different components of the response 1PDMs and the corresponding NOs within both the undamped and damped response theory framework. We illustrate the utility of this new tool by deconstructing β for cis-difluoroethene, para-nitroaniline, and hemibonded OH^· + H₂O complex, computed within the framework of coupled-cluster singles and doubles response theory, in terms of the underlying response 1PDMs and NOs for a range of frequencies.

Interpretation of Coupled-Cluster Many-Electron Dynamics in Terms of Stationary States

We demonstrate theoretically and numerically that laser-driven many-electron dynamics, as described by bivariational time-dependent coupled-cluster (CC) theory, may be analyzed in terms of stationary-state populations. Projectors heuristically defined from linear response theory and equation-of-motion CC theory are proposed for the calculation of stationary-state populations during interaction with laser pulses or other external forces, and conservation laws of the populations are discussed. Numerical tests of the proposed projectors, involving both linear and nonlinear optical processes for He and Be atoms and for LiH, CH+, and LiF molecules show that the laser-driven evolution of the stationary-state populations at the coupled-cluster singles-and-doubles (CCSD) level is very close to that obtained by full configuration interaction (FCI) theory, provided that all stationary states actively participating in the dynamics are sufficiently well approximated. When double-excited states are important for the dynamics, the quality of the CCSD results deteriorates. Observing that populations computed from the linear response projector may show spurious small-amplitude, high-frequency oscillations, the equation-of-motion projector emerges as the most promising approach to stationary-state populations.

A simple molecular orbital picture of RIXS distilled from many-body damped response theory

Ab initio calculations of resonant inelastic x-ray scattering (RIXS) often rely on damped response theory, which prevents the divergence of response solutions in the resonant regime. Within the damped response theory formalism, RIXS moments are expressed as the sum over all electronic states of the system [sum-over-states (SOS) expressions]. By invoking resonance arguments, this expression can be reduced to a few terms, an approximation commonly exploited for the interpretation of computed cross sections. We present an alternative approach: a rigorous formalism for deriving a simple molecular orbital picture of the RIXS process from many-body calculations using the damped response theory. In practical implementations, the SOS expressions of RIXS moments are recast in terms of matrix elements between the zero-order wave functions and first-order frequency-dependent response wave functions of the initial and final states such that the RIXS moments can be evaluated using complex response one-particle transition density matrices (1PTDMs). Visualization of these 1PTDMs connects the RIXS process with the changes in electronic density. We demonstrate that the real and imaginary components of the response 1PTDMs can be interpreted as contributions of the undamped off-resonance and damped near-resonance SOS terms, respectively. By analyzing these 1PTDMs in terms of natural transition orbitals, we derive a rigorous, black-box mapping of the RIXS process into a molecular orbital picture. We illustrate the utility of the new tool by analyzing RIXS transitions in the OH radical, benzene, para-nitroaniline, and 4-amino-4'-nitrostilbene. These examples highlight the significance of both the near-resonance and off-resonance channels.

How to stay out of trouble in RIXS calculations within equation-of-motion coupled-cluster damped response theory? Safe hitchhiking in the excitation manifold by means of core–valence separation

We present a novel approach with robust convergence of the response equations for computing resonant inelastic X-ray scattering (RIXS) cross sections within the equation-of-motion coupled-cluster (EOM-CC) framework.

New and Efficient Equation-of-Motion Coupled-Cluster Framework for Core-Excited and Core-Ionized States

We present a fully analytical implementation of the core-valence separation (CVS) scheme for the equation-of-motion (EOM) coupled-cluster singles and doubles (CCSD) method for calculations of core-level states. Inspired by the CVS idea as originally formulated by Cederbaum, Domcke, and Schirmer, pure valence excitations are excluded from the EOM target space and the frozen-core approximation is imposed on the reference-state amplitudes and multipliers. This yields an efficient, robust, practical, and numerically balanced EOM-CCSD framework for calculations of excitation and ionization energies as well as state and transition properties (e.g., spectral intensities, natural transition, and Dyson orbitals) from both the ground and excited states. The errors in absolute excitation/ionization energies relative to the experimental reference data are on the order of 0.2-3.0 eV, depending on the K-edge considered and on the basis set used, and the shifts are systematic for each edge. Compared to a previously proposed CVS scheme where CVS was applied as a posteriori projection only during the solution of the EOM eigenvalue equations, the new scheme is computationally cheaper. It also achieves better cancellation of errors, yielding similar spectral profiles but with absolute core excitation and ionization energies that are systematically closer to the corresponding experimental data. Among the presented results are calculations of transient-state X-ray absorption spectra, relevant for interpretation of UV-pump/X-ray probe experiments.