A theoretical and experimental benchmark study of core-excited states in nitrogen

The x-ray absorption spectrum of the tert-butyl radical: An experimental and computational investigation

We report the x-ray absorption spectrum (XAS) of the tert-butyl radical, C4H9. The radical was generated pyrolytically from azo-tert-butane, and the XAS of the pure radical was obtained by subtraction of spectra recorded at different temperatures. The bands in the XAS were assigned by ab initio calculations that are in very good agreement with the experimental data. The lowest energy signal in the XAS is assigned to the C1s electron transition from the central carbon atom to the singly occupied molecular orbital (SOMO), while higher transitions correspond to C1s excitations from terminal carbon atoms. Furthermore, we investigated the fragmentation of the radical following resonant C1s excitation by electron-ion-coincidence spectroscopy. Several fragmentation channels were identified. The C1s excitation of the terminal carbons is associated with a stronger fragmentation tendency compared to the lowest C1s excitation of the central carbon into the SOMO. For this core excited state, we still observe an intact parent ion, C4H9+, and a comparatively higher tendency to dissociate into CH3+ + C3H6+.

Core-ionization spectrum of liquid water

We present state-of-the-art calculations of the core-ionization spectrum of water. Despite significant progress in procedures developed to mitigate various experimental complications and uncertainties, the experimental determination of ionization energies of solvated species involves several non-trivial steps such as assessing the effect of the surface potential, electrolytes, and finite escape depths of photoelectrons. This provides a motivation to obtain robust theoretical values of the intrinsic bulk ionization energy and the corresponding solvent-induced shift. Here we develop theoretical protocols based on coupled-cluster theory and electrostatic embedding. Our value of the intrinsic solvent-induced shift of the 1sO ionization energy of water is -1.79 eV. The computed absolute position and the width of the 1sO peak in photoelectron spectrum of water are 538.47 eV and 1.44 eV, respectively, agreeing well with the best experimental values.

Electronic Characterization of Glycolaldehyde: Experimental and Theoretical Insights from the Core- and Valence-Level Spectroscopy

The core-level electron excitation and ionization spectra of glycolaldehyde have been investigated by photoabsorption and photoemission spectroscopy at both carbon and oxygen K-edges; the valence ionization spectra were also recorded by photoelectron spectroscopy in the UV-vis region. The spectra are interpreted by means of ab initio calculations based on the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) and coupled cluster singles, doubles, and perturbative are in good agreement with the experimental results, and many of the observed features are assigned. The photoabsorption spectra are not only dominated by transitions from core-level orbitals to unoccupied π and σ orbitals but also show structures due to Rydberg transitions.

Signatures of diradicals in x-ray absorption spectroscopy

Theoretical simulations are critical to analyze and interpret the x-ray absorption spectrum of transient open-shell species. In this work, we propose a model of the many-body core-excited states of symmetric diradicals. We apply this model to analyze the carbon K-edge transitions of o-, m-, and p-benzyne, three organic diradicals with diverse and unusual electronic structures. The predictions of our model are compared with high-level multireference computations of the K-edge spectrum of the benzynes obtained with the driven similarity renormalization group truncated to third order. Our model shows the importance of a many-body treatment of the core-excited states of the benzynes and provides a theoretical framework to understand which properties of the ground state of these diradicals can be extracted from their x-ray spectrum.

Oscillator Strengths in the Framework of Equation of Motion Multilevel CC3

We present an efficient implementation of the equation of motion oscillator strengths for the closed-shell multilevel coupled cluster singles and doubles with perturbative triples method (MLCC3) in the electronic structure program e^T. The orbital space is split into an active part treated with CC3 and an inactive part computed at the coupled cluster singles and doubles (CCSD) level of theory. Asymptotically, the CC3 contribution scales as floating-point operations, where n_V is the total number of virtual orbitals while n_v and n_o are the number of active virtual and occupied orbitals, respectively. The CC3 contribution, thus, only scales linearly with the full system size and can become negligible compared to the cost of CCSD. We demonstrate the capabilities of our implementation by calculating the ultraviolet–visible spectrum of azobenzene and a core excited state of betaine 30 with more than 1000 molecular orbitals.

Assessing the performance of ΔSCF and the diagonal second-order self-energy approximation for calculating vertical core excitation energies

Vertical core excitation energies are obtained using a combination of the ΔSCF method and the diagonal second-order (D2) self-energy approximation. These methods are applied to a set of neutral molecules and their anionic forms. An assessment of the results with the inclusion of relativistic effects is presented. For core excitations involving delocalized symmetry orbitals, the applied composite method improves upon the overestimation of ΔSCF by providing approximate values close to experimental K-shell transition energies. The importance of both correlation and relaxation contributions to the vertical core-excited state energies, the concept of local and non-local core orbitals, and the consequences of breaking symmetry are discussed.

Computing x-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn-Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ∼0.6 eV rms error for the K-edges of C-F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

Predictions of Pre-edge Features in Time-Resolved Near-Edge X-ray Absorption Fine Structure Spectroscopy from Hole-Hole Tamm-Dancoff-Approximated Density Functional Theory

Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S₂ (ππ*) and S₁ (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N²) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.

Sulfur Molecules in Space by X-rays: A Computational Study

X-ray astronomy lacks high resolution spectra of interstellar dust analogues and molecules, severely hampering interstellar medium studies based on upcoming X-ray missions. Various theoretical approaches may be used to address this problem, but they must first be shown to reproduce reliable spectra compared to the experiment. In this work, we calculate the sulfur K edge X-ray absorption spectra of H₂S, SO₂, and OCS, whose spectra are already known from X-ray experiments and predict the X-ray spectrum of CS, which as far as we are aware has not been measured, thereby hampering its detection by X-ray telescopes. We chose these four molecules as the astrochemistry of sulfur is an unsolved problem and as the four molecules are already known to exist in space. We consider three types of methods for modeling the X-ray spectra: more accurate calculations with the algebraic-diagrammatic construction (ADC) and the CC2, CCSD, and CC3 coupled cluster (CC) approaches as well as more affordable ones with transition potential density functional theory (TP-DFT). A comparison of our computational results to previously reported experimental spectra shows that the core-valence separation (CVS) approaches CVS-ADC(2)-x and CVS-CC3 generally yield a good qualitative level of agreement with the experiment, suggesting that they can be used for interpreting measured spectra, while the TP-DFT method is not reliable for these molecules. However, quantitative agreement with the experiment is still outside the reach of the computational methods studied in this work.

XABOOM: An X-ray Absorption Benchmark of Organic Molecules Based on Carbon, Nitrogen, and Oxygen 1s → π* Transitions

The performance of several standard and popular approaches for calculating X-ray absorption spectra at the carbon, nitrogen, and oxygen K-edges of 40 primarily organic molecules up to the size of guanine has been evaluated, focusing on the low-energy and intense 1s → π* transitions. Using results obtained with CVS-ADC(2)-x and fc-CVS-EOM-CCSD as benchmark references, we investigate the performance of CC2, ADC(2), ADC(3/2), and commonly adopted density functional theory (DFT)-based approaches. Here, focus is on precision rather than on accuracy of transition energies and intensities-in other words, we target relative energies and intensities and the spread thereof, rather than absolute values. The use of exchange-correlation functionals tailored for time-dependent DFT calculations of core excitations leads to error spreads similar to those seen for more standard functionals, despite yielding superior absolute energies. Long-range corrected functionals are shown to perform particularly well compared to our reference data, showing error spreads in energy and intensity of 0.2-0.3 eV and ∼10%, respectively, as compared to 0.3-0.6 eV and ∼20% for a typical pure hybrid. In comparing intensities, state mixing can complicate matters, and techniques to avoid this issue are discussed. Furthermore, the influence of basis sets in high-level ab initio calculations is investigated, showing that reasonably accurate results are obtained with the use of 6-311++G**. We name this benchmark suite as XABOOM (X-ray absorption benchmark of organic molecules) and provide molecular structures and ground-state self-consistent field energies and spectroscopic data. We believe that it provides a good assessment of electronic structure theory methods for calculating X-ray absorption spectra and will become useful for future developments in this field.

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.

New and Efficient Implementation of CC3

We present a new and efficient implementation of the closed shell coupled cluster singles and doubles with perturbative triples method (CC3) in the electronic structure program eT. Asymptotically, a ground state calculation has an iterative cost of 4nV4nO3 floating point operations (FLOP), where nV and nO are the number of virtual and occupied orbitals, respectively. The Jacobian and transpose Jacobian transformations, required to iteratively solve for excitation energies and transition moments, both require 8nV4nO3 FLOP. We have also implemented equation of motion (EOM) transition moments for CC3. The EOM transition densities require recalculation of triples amplitudes, as nV3nO3 tensors are not stored in memory. This results in a noniterative computational cost of 10nV4nO3 FLOP for the ground state density and 26nV4nO3 FLOP per state for the transition densities. The code is compared to the CC3 implementations in CFOUR, DALTON, and PSI4. We demonstrate the capabilities of our implementation by calculating valence and core excited states of l-proline.

Equation-of-Motion MLCCSD and CCSD-in-HF Oscillator Strengths and Their Application to Core Excitations

We present an implementation of equation-of-motion oscillator strengths for the multilevel CCSD (MLCCSD) model where CCS is used as the lower level method (CCS/CCSD). In this model, the double excitations of the cluster operator are restricted to an active orbital space, whereas the single excitations are unrestricted. Calculated nitrogen K-edge spectra of adenosine, adenosine triphosphate (ATP), and an ATP-water system are used to demonstrate the performance of the model. Projected atomic orbitals (PAOs) are used to partition the virtual space into active and inactive orbital sets. Cholesky decomposition of the Hartree-Fock density is used to partition the occupied orbitals. This Cholesky-PAO partitioning is cheap, scaling as O(N3), and is suitable for the calculation of core excitations, which are localized in character. By restricting the single excitations of the cluster operator to the active space, as well as the double excitations, the CCSD-in-HF model is obtained. A comparison of the two models-MLCCSD and CCSD-in-HF-is presented for the core excitation spectra of the adenosine and ATP systems.

Core excitations with excited state mean field and perturbation theory

We test the efficacy of excited state mean field theory and its excited-state-specific perturbation theory on the prediction of K-edge positions and x-ray peak separations. We find that the mean field theory is surprisingly accurate, even though it contains no accounting of differential electron correlation effects. In the perturbation theory, we test multiple core-valence separation schemes and find that, with the mean field theory already so accurate, electron-counting biases in one popular separation scheme become a dominant error when predicting K-edges. Happily, these appear to be relatively easy to correct for, leading to a perturbation theory for K-edge positions that is lower scaling and more accurate than coupled cluster theory and competitive in accuracy with recent high-accuracy results from restricted open-shell Kohn-Sham theory. For peak separations, our preliminary data show excited state mean field theory to be exceptionally accurate, but more extensive testing will be needed to see how it and its perturbation theory compare to coupled cluster peak separations more broadly.

Highly Accurate Prediction of Core Spectra of Molecules at Density Functional Theory Cost: Attaining Sub-electronvolt Error from a Restricted Open-Shell Kohn-Sham Approach

We present the use of the recently developed square gradient minimization (SGM) algorithm for excited-state orbital optimization to obtain spin-pure restricted open-shell Kohn-Sham (ROKS) energies for core excited states of molecules. The SGM algorithm is robust against variational collapse and offers a reliable route to converging orbitals for target excited states at only 2-3 times the cost of ground-state orbital optimization (per iteration). ROKS/SGM with the modern SCAN/ωB97X-V functionals is found to predict the K-edge of C, N, O, and F to a root mean squared error of ∼0.3 eV. ROKS/SGM is equally effective at predicting L-edge spectra of third period elements, provided a perturbative spin-orbit correction is employed. This high accuracy can be contrasted with traditional time-dependent density functional theory (TDDFT), which typically has greater than 10 eV error and requires translation of computed spectra to align with experiment. ROKS is computationally affordable (having the same scaling as ground-state DFT and a slightly larger prefactor) and can be applied to geometry optimizations/ab initio molecular dynamics of core excited states, as well as condensed phase simulations. ROKS can also model doubly excited/ionized states with one broken electron pair, which are beyond the ability of linear response based methods.

X-ray and UV Spectra of Glycine within Coupled Cluster Linear Response Theory

The coupled cluster models CCSD and CC3 are used to investigate the (core) excited states and ionization energies of glycine in the gas phase. Excited states and ionization energies in the UV spectral range are calculated using a standard coupled cluster linear response, while core-level excited states and ionization potentials are calculated using the core-valence separation approximation. The temperature dependence from different conformers is also assessed.

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.