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

Molecular-Orbital Framework of Two-Electron Processes: Application to Auger and Intermolecular Coulomb Decay

States with core- or inner-shell vacancies, which are commonly created by absorption of high-energy photons, can decay by a two-electron process in which one electron fills the core hole and the second one is ejected. These processes accompany many X-ray spectroscopies. Depending on the nature of the initial core- or inner-shell-hole state and the decay valence-hole state, these processes are called Auger decay, intermolecular Coulomb decay, or electron-transfer-mediated decay. To connect many-body wave functions of the initial and final states with the molecular orbital picture of the decay, we introduce the concept of natural Auger orbitals (NAOs). NAOs are obtained by a two-step singular value decomposition of the two-body Dyson orbitals, reduced quantities that enter the expression of the decay rate in the Feshbach-Fano treatment. NAOs afford chemical insight and interpretation of the high-level ab initio calculations of Auger decay and related two-electron relaxation processes.

Natural orbitals and two-particle correlators as tools for the analysis of effective exchange couplings in solids

Using generalizations of spin-averaged natural orbitals and two-particle charge correlators for solids, we investigate the electronic structure of antiferromagnetic transition-metal oxides with a fully self-consistent, imaginary-time GW method. Our findings disagree with the Goodenough-Kanamori (GK) rules that are commonly used for the qualitative interpretation of such solids. First, we found a strong dependence of the natural orbital occupancies on momenta, contradicting GK assumptions. Second, along the momentum path, the character of natural orbitals changes. In particular, the contributions of oxygen 2s orbitals are important, which has not been considered in the GK rules. To analyze the influence of the electronic correlation on the values of effective exchange coupling constants, we use both natural orbitals and two-particle correlators and show that electronic screening modulates the degree of superexchange by stabilizing the charge-transfer contributions, which greatly affects these coupling constants. Finally, we give a set of predictions and recommendations regarding the use of density functional, Green's function, and wave-function methods for evaluating effective magnetic couplings in molecules and solids.

The Auger spectrum of benzene

We present an ab initio computational study of the Auger electron spectrum of benzene. Auger electron spectroscopy exploits the Auger-Meitner effect, and although it is established as an analytic technique, the theoretical modeling of molecular Auger spectra from first principles remains challenging. Here, we use coupled-cluster theory and equation-of-motion coupled-cluster theory combined with two approaches to describe the decaying nature of core-ionized states: (i) Feshbach-Fano resonance theory and (ii) the method of complex basis functions. The spectra computed with these two approaches are in excellent agreement with each other and also agree well with experimental Auger spectra of benzene. The Auger spectrum of benzene features two well-resolved peaks at Auger electron energies above 260 eV, which correspond to final states with two electrons removed from the 1e_1g and 3e_2g highest occupied molecular orbitals. At lower Auger electron energies, the spectrum is less well resolved, and the peaks comprise multiple final states of the benzene dication. In line with theoretical considerations, singlet decay channels contribute more to the total Auger intensity than the corresponding triplet decay 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.

Computational approaches for XANES, VtC-XES, and RIXS using linear-response time-dependent density functional theory based methods

The emergence of state-of-the-art X-ray light sources has paved the way for novel spectroscopies that take advantage of their atomic specificity to shed light on fundamental physical, chemical, and biological processes both in the static and time domains. The success of these experiments hinges on the ability to interpret and predict core-level spectra, which has opened avenues for theory to play a key role. Over the last two decades, linear-response time-dependent density functional theory (LR-TDDFT), despite various theoretical challenges, has become a computationally attractive and versatile framework to study excited-state spectra including X-ray spectroscopies. In this context, we focus our discussion on LR-TDDFT approaches for the computation of X-ray Near-Edge Structure (XANES), Valence-to-Core X-ray Emission (VtC-XES), and Resonant Inelastic X-ray Scattering (RIXS) spectroscopies in molecular systems with an emphasis on Gaussian basis set implementations. We illustrate these approaches with applications and provide a brief outlook of possible new directions.

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.

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.

Deciphering excited state properties utilizing algebraic diagrammatic construction schemes of decreasing order

Excited state properties are difficult to trace back to the common molecular orbital picture when the excited state wavefunction is a linear combination of two or more Slater determinants. Here, a theoretical methodology is introduced based on the algebraic diagrammatic construction scheme for the polarization propagator (ADC(n)) that allows to make this connection and to eventually derive structure-function relationships. The usefulness of this approach is demonstrated by an analysis of the transition dipole moments of the low-lying 1B_3u and 2B_3u states of anthracene and (1,4,5,8)-tetraazaanthracene.

Resonant Inelastic X-ray Scattering Calculations of Transition Metal Complexes Within a Simplified Time-Dependent Density Functional Theory Framework

We present a time-dependent density functional theory (TDDFT) approach to compute the light-matter couplings between two different manifolds of excited states relative to a common ground state in the context of 4d transition metal systems. These quantities are the necessary ingredients to solve the Kramers-Heisenberg (KH) equation for resonant inelastic X-ray scattering (RIXS) and several other types of two-photon spectroscopies. The procedure is based on the pseudo-wavefunction approach, where the solutions of a TDDFT calculation can be used to construct excited-state wavefunctions, and on the restricted energy window approach, where a manifold of excited states can be rigorously defined based on the energies of the occupied molecular orbitals involved in the excitation process. Thus, the present approach bypasses the need to solve the costly TDDFT quadratic-response equations. We illustrate the applicability of the method to 4d transition metal molecular complexes by calculating the 2p4d RIXS maps of three representative ruthenium complexes and comparing them to experimental results. The method can capture all the experimental features in all three complexes to allow the assignment of the experimental peaks, with relative energies correct to within ∼0.6 eV at the cost of two independent TDDFT calculations.

Feshbach-Fano approach for calculation of Auger decay rates using equation-of-motion coupled-cluster wave functions. I. Theory and implementation

X-ray absorption creates electron vacancies in the core shell. These highly excited states often relax by Auger decay-an autoionization process in which one valence electron fills the core hole and another valence electron is ejected into the ionization continuum. Despite the important role of Auger processes in many experimental settings, their first-principles modeling is challenging, even for small systems. The difficulty stems from the need to describe many-electron continuum (unbound) states, which cannot be tackled with standard quantum-chemistry methods. We present a novel approach to calculate Auger decay rates by combining Feshbach-Fano resonance theory with the equation-of-motion coupled-cluster single double (EOM-CCSD) framework. We use the core-valence separation scheme to define projectors into the bound (square-integrable) and unbound (continuum) subspaces of the full function space. The continuum many-body decay states are represented by products of an appropriate EOM-CCSD state and a free-electron state, described by a continuum orbital. The Auger rates are expressed in terms of reduced quantities, two-body Dyson amplitudes (objects analogous to the two-particle transition density matrix), contracted with two-electron bound-continuum integrals. Here, we consider two approximate treatments of the free electron: a plane wave and a Coulomb wave with an effective charge, which allow us to evaluate all requisite integrals analytically; however, the theory can be extended to incorporate a more sophisticated description of the continuum orbital.

Equation-of-Motion Coupled-Cluster Theory to Model L-Edge X-ray Absorption and Photoelectron Spectra

We present an extension of the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) theory for computing X-ray L-edge spectra, both in the absorption (XAS) and in the photoelectron (XPS) regimes. The approach is based on the perturbative evaluation of spin-orbit couplings using the Breit-Pauli Hamiltonian and nonrelativistic wave functions described by the fc-CVS-EOM-CCSD ansatz (EOM-CCSD within the frozen-core core-valence separated (fc-CVS) scheme). The formalism is based on spinless one-particle density matrices. The approach is illustrated by modeling XAS and XPS of several model systems ranging from Ar to small molecules containing sulfur and silicon.