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.

Efficient Protocol for Computing MCD Spectra in a Broad Frequency Range Combining Resonant and Damped CC2 Quadratic Response Theory

Coupled cluster response theory offers a path to high-accuracy calculations of spectroscopic properties, such as magnetic circular dichroism (MCD). However, divergence or slow convergence issues are often encountered for electronic transitions in high-energy regions with a high density of states. This is here addressed for MCD by an implementation of damped quadratic response theory for resolution-of-identity coupled cluster singles-and-approximate-doubles (RI-CC2), along with an implementation of the MCD A term from resonant response theory. Combined, damped and resonant response theory calculations provide an efficient strategy to obtain MCD spectra over a broad frequency range and for systems that include highly symmetric molecules with degenerate excited states. The protocol is illustrated by application to zinc tetrabenzoporphyrin in the energy region of 2-8 eV and comparison to experimental data. Timings are reported for the resonant and damped approaches, showing that a greater part of the calculation time is consumed by the construction of the building blocks for the final MCD ellipticity. A recommendation on how to use the procedure is outlined.

Ab Initio Investigation of the Auger Spectra of Methane, Ethane, Ethylene, and Acetylene

We present an ab initio computational study of the Auger spectra of methane, ethane, ethylene, and acetylene. Auger spectroscopy is an established technique to probe the electronic structure of molecules and exploits the Auger-Meitner effect that core-ionized states undergo. We compute partial decay widths using coupled-cluster theory with single and double substitutions (CCSD) and equation-of-motion CCSD theory combined with complex-scaled basis functions and Feshbach-Fano projection. We generate Auger spectra from these partial widths and draw conclusions about the strength of particular decay channels and trends among the four molecules. A connection to experimental results about fragmentation pathways of the electronic states produced by Auger decay is also made.

Solving response expressions in the ADC/ISR framework

We present an implementation for the calculation of molecular response properties using the algebraic-diagrammatic construction (ADC)/intermediate state representation approach. For the second-order ADC model [ADC(2)], a memory-efficient ansatz avoiding the storage of double excitation amplitudes is investigated. We compare the performance of different numerical algorithms for the solution of the underlying response equations for ADC(2) and show that our approach also strongly improves the convergence behavior for the investigated algorithms compared with the standard implementation. All routines are implemented in an open-source Python library.

Theoretical Investigation of the X-ray Stark Effect in Small Molecules

We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.

New Implementation of an Equation-of-Motion Coupled-Cluster Damped-Response Framework with Illustrative Applications to Resonant Inelastic X-ray Scattering

We present an implementation of a damped response framework for calculating resonant inelastic X-ray scattering (RIXS) at the equation-of-motion coupled-cluster singles and doubles (CCSD) and second-order approximate coupled-cluster singles and doubles (CC2) levels of theory in the open-source program e^T. This framework lays the foundation for future extension to higher excitation methods (notably, the coupled-cluster singles and doubles with perturbative triples, CC3) and to multilevel approaches. Our implementation adopts a fully relaxed ground state and different variants of the core-valence separation projection technique to address convergence issues. Illustrative results are compared with those obtained within the frozen-core core-valence separated approach, available in Q-Chem, as well as with experiment. The performance of the CC2 method is evaluated in comparison with that of CCSD. It is found that, while the CC2 method is noticeably inferior to CCSD for X-ray absorption spectra, the quality of the CC2 RIXS spectra is often comparable to that of the CCSD level of theory, when the same valence excited states are probed. Finally, we present preliminary RIXS results for a solvated molecule in aqueous solution.

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.

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.

Molecular Chirality and Its Monitoring by Ultrafast X-ray Pulses

Major advances in X-ray sources including the development of circularly polarized and orbital angular momentum pulses make it possible to probe matter chirality at unprecedented energy regimes and with Ångström and femtosecond spatiotemporal resolutions. We survey the theory of stationary and time-resolved nonlinear chiral measurements that can be carried out in the X-ray regime using tabletop X-ray sources or large scale (XFEL, synchrotron) facilities. A variety of possible signals and their information content are discussed.

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.

Pair Natural Orbital Equation-of-Motion Coupled-Cluster Method for Core Binding Energies: Theory, Implementation, and Benchmark

We present the theory and implementation of a lower scaling core-valence separated equation-of-motion coupled-cluster approach based on domain-based local pair natural orbitals for core binding energies. The accuracy of the new method has been compared with that of the standard equation-of-motion coupled-cluster method and experimentally measured results. The use of pair natural orbitals significantly reduces the computation cost and can be applied to large molecules.

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.

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

Molecular Auger Decay Rates from Complex-Variable Coupled-Cluster Theory

The emission of an Auger electron is the predominant relaxation mechanism of core-vacant states in molecules composed of light nuclei. In this non-radiative decay process, one valence electron fills the core vacancy while a second valence electron is emitted into the ionization continuum. Because of this coupling to the continuum, core-vacant states represent electronic resonances that can be tackled with standard quantum-chemical methods only if they are approximated as bound states, meaning that Auger decay is neglected. Here, we present an approach to compute Auger decay rates of core-vacant states from coupled-cluster and equation-of-motion coupled-cluster wave functions combined with complex scaling of the Hamiltonian or, alternatively, complex-scaled basis functions. Through energy decomposition analysis, we illustrate how complex-scaled methods are capable of describing the coupling to the ionization continuum without the need to model the wave function of the Auger electron explicitly. In addition, we introduce in this work several approaches for the determination of partial decay widths and Auger branching ratios from complex-scaled coupled-cluster wave functions. We demonstrate the capabilities of our new approach by computations on core-ionized states of neon, water, dinitrogen, and benzene. Coupled-cluster and equation-of-motion coupled-cluster theory in the singles and doubles approximation both deliver excellent results for total decay widths, whereas we find partial widths more straightforward to evaluate with the former method. We also observe that the requirements towards the basis set are less arduous for Auger decay than for other types of resonances so that extensions to larger molecules are readily possible.

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.

All-Electron many-body approach to resonant inelastic x-ray scattering

We present a formalism for the resonant inelastic x-ray scattering (RIXS) cross section. The resulting compact expression in terms of polarizability matrix elements, particularly lends itself to the implementation in...

A Core-Valence Separated Similarity Transformed EOM-CCSD Method for Core-Excitation Spectra

We present the theory and implementation of a core-valence separated similarity transformed EOM-CCSD (STEOM-CCSD) method for K-edge core excitation spectra. The method can select an appropriate active space using CIS natural orbitals and near "black box" to use. The second similarity transformed Hamiltonian is diagonalized in the space of single excitation. Therefore, the final diagonalization step is free from the convergence problem arising due to the coupling of the core-excited states with the continuum of doubly excited states. Convergence trouble can appear for the preceding core-ionized state calculation in STEOM-CCSD. A core-valence separation (CVS) scheme compatible with the natural orbital based active space selection (CVS-STEOM-CCSD-NO) is implemented to overcome the problem. The CVS-STEOM-CCSD-NO has a similar accuracy to that of the standard CVS-EOM-CCSD method but comes with a lower computational cost. The modification required in the CVS scheme to make use of the CIS natural orbital is highlighted. The suitability of the CVS-STEOM-CCSD-NO method for chemical application is demonstrated by simulating the K-edge spectra of glycine and thymine.

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.

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.

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.

Probing Basis Set Requirements for Calculating Core Ionization and Core Excitation Spectra Using Correlated Wave Function Methods

We investigate the basis set requirements for the accurate calculation of core excitations and core ionizations using correlated wave functions of coupled cluster type and linear response methods for describing the excitation. When a core excitation is described as an energy difference calculated using density functional theory, the basis set can be tailored to provide a balanced description of the reference- and excited-hole states. When the core excitation process is described by coupled cluster linear response methods, however, the basis set requirements are somewhat different. A systematic study of the sensitivity of the result to the basis set parameters suggests that a relatively large set of s- and p-type basis functions in combination with a careful selection of valence and core polarization functions is required. Based on these results, we propose a hierarchical sequence of basis sets, denoted ccX-nZ (n = D, T, Q, 5) for the atoms B-Ne, which are suitable for the calculation of core excitations by the correlated wave function linear response and equation-of-motion methods. The ccX-nZ series provides lower basis set errors for a given cardinal number or number of basis functions than other existing basis sets. For large systems, the ccX-nZ basis sets can be combined with the standard basis sets by placing the ccX-nZ only on the atoms where core excitations are of interest, but the accuracy of such mixed basis sets appears to be system-dependent.

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.

Damped (linear) response theory within the resolution-of-identity coupled cluster singles and approximate doubles (RI-CC2) method

An implementation of a complex solver for the solution of the linear equations required to compute the complex response functions of damped response theory is presented for the resolution-of-identity (RI) coupled cluster singles and approximate doubles (CC2) method. The implementation uses a partitioned formulation that avoids the storage of double excitation amplitudes to make it applicable to large molecules. The solver is the keystone element for the development of the damped coupled cluster response formalism for linear and nonlinear effects in resonant frequency regions at the RI-CC2 level of theory. Illustrative results are reported for the one-photon absorption cross section of C₆₀, the electronic circular dichroism of n-helicenes (n = 5, 6, 7), and the C₆ dispersion coefficients of a set of selected organic molecules and fullerenes.

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.

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.

Table-Top X-ray Spectroscopy of Benzene Radical Cation

Ultrafast table-top X-ray spectroscopy at the carbon K-edge is used to measure the X-ray spectral features of benzene radical cations (Bz⁺). The ground state of the cation is prepared selectively by two-photon ionization of neutral benzene, and the X-ray spectra are probed at early times after the ionization by transient absorption using X-rays produced by high harmonic generation (HHG). Bz⁺ is well-known to undergo Jahn-Teller distortion, leading to a lower symmetry and splitting of the π orbitals. Comparison of the X-ray absorption spectra of the neutral and the cation reveals a splitting of the two degenerate π* orbitals as well as an appearance of a new peak due to excitation to the partially occupied π-subshell. The π* orbital splitting of the cation, elucidated on the basis of high-level calculations in a companion theoretical paper [Vidal et al. J. Phys. Chem. A. http://dx.doi.org/10.1021/acs.jpca.0c08732], is discovered to be due to both the symmetry distortion and even more dominant spin coupling of the unpaired electron in the partially vacant π orbital (from ionization) with the unpaired electrons resulting from the transition from the 1s_C core orbital to the fully vacant π* orbitals.

Interplay of Open-Shell Spin-Coupling and Jahn-Teller Distortion in Benzene Radical Cation Probed by X-ray Spectroscopy

We report a theoretical investigation and elucidation of the X-ray absorption spectra of neutral benzene and of the benzene cation. The generation of the cation by multiphoton ultraviolet (UV) ionization and the measurement of the carbon K-edge spectra of both species using a table-top high-harmonic generation source are described in the companion experimental paper [Epshtein, M.; et al. J. Phys. Chem. A http://dx.doi.org/10.1021/acs.jpca.0c08736]. We show that the 1s_C → π transition serves as a sensitive signature of the transient cation formation, as it occurs outside of the spectral window of the parent neutral species. Moreover, the presence of the unpaired (spectator) electron in the π-subshell of the cation and the high symmetry of the system result in significant differences relative to neutral benzene in the spectral features associated with the 1s_C → π* transitions. High-level calculations using equation-of-motion coupled-cluster theory provide the interpretation of the experimental spectra and insight into the electronic structure of benzene and its cation. The prominent split structure of the 1s_C → π* band of the cation is attributed to the interplay between the coupling of the core → π* excitation with the unpaired electron in the π-subshell and the Jahn-Teller distortion. The calculations attribute most of the splitting (∼1-1.2 eV) to the spin coupling, which is visible already at the Franck-Condon structure, and we estimate the additional splitting due to structural relaxation to be around ∼0.1-0.2 eV. These results suggest that X-ray absorption with increased resolution might be able to disentangle electronic and structural aspects of the Jahn-Teller effect in the benzene cation.

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.

Density Functional Theory Based Methods for the Calculation of X-ray Spectroscopy

The availability of new light sources combined with the realization of the unique capabilities of spectroscopy in the X-ray region has driven tremendous advances in the field of X-ray spectroscopy. Currently, these techniques are emerging as powerful analytical tools for the study of a wide range of problems encompassing liquids, materials, and biological systems. Time-resolved measurements add a further dimension to X-ray spectroscopy, opening up the potential to resolve ultrafast chemical processes at an atomic level. X-ray spectroscopy encompasses a range of techniques which provide complementary information, and these include X-ray photoelectron spectroscopy (XPS), X-ray absorption spectroscopy (XAS), X-ray emission spectroscopy (XES), and resonant inelastic X-ray scattering (RIXS). In many studies, the interpretation of the experimental data relies upon calculations to enable the nature of the underlying molecular structure, electronic structure, and bonding to be revealed. Density functional theory (DFT) based methods are some of the most widely used methods for the simulation of X-ray spectra. In this Account, we focus on our recent contributions to the simulation of a range of X-ray spectroscopic techniques using DFT and linear-response time-dependent density functional theory (TDDFT) and show how these methods can provide a computational toolkit for the simulation of X-ray spectroscopy. The importance of the exchange-correlation functional for the calculation of XAS is discussed, and the introduction of short-range corrected functionals is described. The application of these calculations to study large systems through the use of efficient implementations of TDDFT will be highlighted, with the use of these methods illustrated through studies of ionic liquids and transition metal complexes. The extension of TDDFT to calculate XES through the use of a reference determinant for the core-ionized state will be described, and the factors that affect the accuracy of the computed spectra discussed. The application of these approaches will be illustrated through the study of a range of organic molecules and transition metal complexes, which also show how going beyond the dipole approximation in determining the transition intensities can be critical. The application of these approaches to the simulation of the RIXS spectrum of water will also be described, highlighting how ultrafast dynamics on the femtoscale time scale are evident in the measured spectra. In these calculations, the description of the core-ionized and core-excited states becomes increasingly important, and the role of the basis set in accurately describing these states will be explored.

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.

Resonant Inelastic X-Ray Scattering Reveals Hidden Local Transitions of the Aqueous OH Radical

Resonant inelastic x-ray scattering (RIXS) provides remarkable opportunities to interrogate ultrafast dynamics in liquids. Here we use RIXS to study the fundamentally and practically important hydroxyl radical in liquid water, OH(aq). Impulsive ionization of pure liquid water produced a short-lived population of OH(aq), which was probed using femtosecond x-rays from an x-ray free-electron laser. We find that RIXS reveals localized electronic transitions that are masked in the ultraviolet absorption spectrum by strong charge-transfer transitions-thus providing a means to investigate the evolving electronic structure and reactivity of the hydroxyl radical in aqueous and heterogeneous environments. First-principles calculations provide interpretation of the main spectral features.

Dyson orbitals within the fc-CVS-EOM-CCSD framework: theory and application to X-ray photoelectron spectroscopy of ground and excited states

Ionization energies and Dyson orbitals within frozen-core core–valence separated equation-of-motion coupled cluster singles and doubles (fc-CVS-EOM-CCSD) enable efficient and reliable calculations of standard XPS and of UV-pump/XPS probe spectra.