Calculating X-ray Absorption Spectra of Open-Shell Molecules with the Unrestricted Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator

GUGA-based MRCI approach with core-valence separation approximation (CVS) for the calculation of the core-excited states of molecules

We develop and demonstrate how to use the Graphical Unitary Group Approach (GUGA)-based MRCISD with Core-Valence Separation (CVS) approximation to compute the core-excited states. First, perform a normal Self-Consistent-Field (SCF) or valence MCSCF calculation to optimize the molecular orbitals. Second, rotate the optimized target core orbitals and append to the active space, form an extended CVS active space, and perform a CVS-MCSCF calculation for core-excited states. Finally, construct the CVS-MRCISD expansion space and perform a CVS-MRCISD calculation to optimize the CI coefficients based on the variational method. The CVS approximation with GUGA-based methods can be implemented by flexible truncation of the Distinct Row Table. Eliminating the valence-excited configurations from the CVS-MRCISD expansion space can prevent variational collapse in the Davidson iteration diagonalization. The accuracy of the CVS-MRCISD scheme was investigated for excitation energies and compared with that of the CVS-MCSCF and CVS-CASPT2 methods using the same active space. The results show that CVS-MRCISD is capable of reproducing well-matched vertical core excitation energies that are consistent with experiments by combining large basis sets and a rational reference space. The calculation results also highlight the fact that the dynamic correlation between electrons makes an undeniable contribution in core-excited states.

Algebraic Diagrammatic Construction Theory for Simulating Charged Excited States and Photoelectron Spectra

Charged excitations are electronic transitions that involve a change in the total charge of a molecule or material. Understanding the properties and reactivity of charged species requires insights from theoretical calculations that can accurately describe orbital relaxation and electron correlation effects in open-shell electronic states. In this Review, we describe the current state of algebraic diagrammatic construction (ADC) theory for simulating charged excitations and its recent developments. We start with a short overview of ADC formalism for the one-particle Green's function, including its single- and multireference formulations and extension to periodic systems. Next, we focus on the capabilities of ADC methods and discuss recent findings about their accuracy for calculating a wide range of excited-state properties. We conclude our Review by outlining possible directions for future developments of this theoretical approach.

A study of core-excited states of organic molecules computed with the generalized active space driven similarity renormalization group

This work examines the accuracy and precision of x-ray absorption spectra computed with a multireference approach that combines generalized active space (GAS) references with the driven similarity renormalization group (DSRG). We employ the x-ray absorption benchmark of organic molecule (XABOOM) set, consisting of 116 transitions from mostly organic molecules [Fransson et al., J. Chem. Theory Comput. 17, 1618 (2021)]. Several approximations to a full-valence active space are examined and benchmarked. Absolute excitation energies and intensities computed with the GAS-DSRG truncated to second-order in perturbation theory are found to systematically underestimate experimental and reference theoretical values. Third-order perturbative corrections significantly improve the accuracy of GAS-DSRG absolute excitation energies, bringing the mean absolute deviation from experimental values down to 0.32 eV. The ozone molecule and glyoxylic acid are particularly challenging for second-order perturbation theory and are examined in detail to assess the importance of active space truncation and intruder states.

Double-Hybrid Density Functional Theory for Core Excitations: Theory and Benchmark Calculations

The double-hybrid (DH) time-dependent density functional theory is extended to core excitations. Two different DH formalisms are presented utilizing the core-valence separation (CVS) approximation. First, a CVS-DH variant is introduced relying on the genuine perturbative second-order correction, while an iterative analogue is also presented using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH ansatz. The performance of the new approaches is tested for the most popular DH functionals using the recently proposed XABOOM [J. Chem. Theory Comput.2021, 17, 1618] benchmark set. In order to make a careful comparison, the accuracy and precision of the methods are also inspected. Our results show that the genuine approaches are highly competitive with the more advanced CVS-ADC(2)-based methods if only excitation energies are required. In contrast, as expected, significant differences are observed in oscillator strengths; however, the precision is acceptable for the genuine functionals as well. Concerning the performance of the CVS-DH approaches, the PBE0-2/CVS-ADC(2) functional is superior, while its spin-opposite-scaled variant is also recommended as a cost-effective alternative. For these approaches, significant improvements are realized in the error measures compared with the popular CVS-ADC(2) method.

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.

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.

Vibronic Coupling Effect on the Vibrationally Resolved Electronic Spectra and Intersystem Crossing Rates of a TADF Emitter: 7-PhQAD

Assessing and improving the performance of organic light-emitting diode (OLED) materials require quantitative prediction of rate coefficients for the intersystem crossing (ISC) and reverse ISC (RISC) processes, which are determined not only by the energy gap and the direct spin-orbit coupling (SOC) between the first singlet and triplet excited-states at a thermal equilibrium position of the initial electronic state but also by the non-Condon effects such as the Herzberg-Teller-like vibronic coupling (HTVC) and the spin-vibronic coupling (SVC). Here we apply the time-dependent correlation function approaches to quantitatively calculate the vibrationally resolved absorption and fluorescence spectra and ISC/RISC rates of a newly synthesized multiple-resonance-type (MR-type) thermally activated delayed fluorescence (TADF) emitter, 7-phenylquinolino[3,2,1-de]acridine-5,9-dione (7-PhQAD), with the inclusion of the Franck-Condon (FC), HTVC, and Duschinsky rotation (DR) effects. The SVC effect on the rates has also been approximately evaluated. We find that the experimentally measured ISC rates of 7-PhQAD originate predominantly from the vibronic coupling, consistent with the previous reports on other MR-type TADF emitters. The SVC effect on ISC rates is about 10 times larger than the HTVC effect, and the latter increases the ISC rates by more than 1 order of magnitude while it slightly affects the vibrationally resolved absorption and fluorescence spectra. The discrepancy between the theoretical and experimental results is attributed to inaccurately describing excited-states calculated by the time-dependent density functional theory as well as to not fully accounting for the complex experimental conditions. This work provides a demonstration of what proportion of ISC and RISC rate coefficients of a MR-type TADF emitter can be covered by the HTVC effect, and it opens design routes that go beyond the FC approximation for the future development of high-performance OLED devices.

Intruder Peak-Free Transient Inner-Shell Spectra Using Real-Time Simulations

Real-time methods are convenient for simulating core-level absorption spectra but suffer from nonphysical intruder peaks when using atom-centered basis sets. In transient absorption spectra, these peaks exhibit highly nonphysical time-dependent modulations in their energies and oscillator strengths. In this paper, we address the origins of these intruder peaks and propose a straightforward and effective solution based on a filtered dipole operator. In combination with real-time time-dependent density functional theory (RT-TDDFT), we demonstrate how to compute intruder-free attosecond transient X-ray absorption spectra for the aminophenol (C₆H₇NO) oxygen and nitrogen K-edges and the α-quartz (SiO₂) silicon L-edge. Without filtering, the computed spectra are qualitatively wrong. This procedure is suitable for both static and transient inner-shell spectroscopy studies and can easily be implemented in a range of real-time methodologies.

Theoretical Calculation of Core-Excited States along Dissociative Pathways beyond Second-Order Perturbation Theory

We extend the multireference driven similarity renormalization (MR-DSRG) method to compute core-excited states by combining it with a GASSCF treatment of orbital relaxation and static electron correlation effects. We consider MR-DSRG treatments of dynamical correlation truncated at the level of perturbation theory (DSRG-MRPT2/3) and iterative linearized approximations with one- and two-body operators [MR-LDSRG(2)] in combination with a spin-free exact-two-component (X2C) one-electron treatment of scalar relativistic effects. This approach is calibrated and tested on a series of 16 core-excited states of five closed- and open-shell diatomic molecules containing first-row elements (C, N, and O). All GASSCF-MR-DSRG theories show excellent agreement with experimental adiabatic transitions energies, with mean absolute errors ranging between 0.17 and 0.35 eV, even for the challenging partially doubly excited states of the N₂⁺ molecule. The vibrational structure of all these transitions, obtained from using a full potential energy scan, shows a mean absolute error as low as 25 meV for DSRG-MRPT2 and 12/13 meV for DSRG-MRPT3 and MR-LDSRG(2). We generally find that a treatment of dynamical correlation that goes beyond the second-order level in perturbation theory improves the accuracy of the potential energy surface, especially in the bond-dissociation region.

Equation of motion coupled-cluster study of core excitation spectra II: Beyond the dipole approximation

We present the time-independent (TI) and time-dependent (TD) equation of motion coupled-cluster (EOM-CC) oscillator strengths not limited to those obtained by the dipole approximation. For the conventional TI-EOM-CC, we implement all the terms in the multipole expansion through second order that contributes to the oscillator strength. These include contributions such as magnetic dipole, electric quadrupole, electric octupole, and magnetic quadrupole. In TD-EOM-CC, we only include the quadrupole moment contributions. This augments our previous work [Y. C. Park, A. Perera, and R. J. Bartlett, J. Chem. Phys. 151, 164117 (2019)]. The inclusion of the quadrupole contributions (and all the other contributions through second order in the case of TI-EOM-CCSD) enables us to obtain the intensities for the pre-edge transitions in the metal K-edge spectra, which are dipole inactive. The TI-EOM-CCSD and TD-EOM-CCSD spectra of Ti⁴⁺ atoms are used to showcase the implementation of the second-order oscillator strengths. The origin of 1s → e and 1s → t₂ in core spectra from iron tetrachloride and titanium tetrachloride is discussed and compared with the experiment.

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.

Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules

Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the M_S = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.

Analytical gradients for core-excited states in the algebraic diagrammatic construction (ADC) framework

Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core-valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn-Teller distorted geometry of the 1s → π* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states.

Unitary coupled-cluster approach for the calculation of core-excited states and x-ray absorption spectra

In this work, we present the core-valence separation (CVS) approximation applied to unitary coupled-cluster (UCC) theory for the calculation of core-excited states and the simulation of x-ray absorption spectroscopy (XAS). Excitation energies and oscillator strengths of small- to medium-sized organic molecules have been computed using the second-order and extended second-order UCC schemes (CVS-UCC2 and CVS-UCC2-x) as well as the third-order scheme (CVS-UCC3). All results are compared to the corresponding algebraic-diagrammatic construction methods and experimental data. The agreement between CVS-UCC and experimental data demonstrates its potential as a new approach for the calculation of XAS.

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.

Quantifying the error of the core-valence separation approximation

For the calculation of core-excited states probed through X-ray absorption spectroscopy, the core-valence separation (CVS) scheme has become a vital tool. This approach allows us to target such states with high specificity, albeit introducing an error. We report the implementation of a post-processing step for CVS excitations obtained within the algebraic-diagrammatic construction scheme for the polarization propagator, which removes this error. Based on this, we provide a detailed analysis of the CVS scheme, identifying its accuracy to be dominated by an error balance between two neglected couplings, one between core and valence single excitations and the other between single and double core excitations. The selection of the basis set is shown to be vital for a proper description of both couplings, with tight polarizing functions being necessary for a good balance of errors. The CVS error is confirmed to be stable across multiple systems, with an element-specific spread for K-edge spectrum calculations of only about ±0.02 eV. A systematic lowering of the CVS error by 0.02 eV-0.03 eV is noted when considering excitations to extremely diffuse states, emulating ionization.

Simulated field-modulated x-ray absorption in titania

In this paper, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin-orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO₂) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t_2g peaks in the Ti L-edge was observed to red shift and the e_g peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.

Third-Order Unitary Coupled Cluster (UCC3) for Excited Electronic States: Efficient Implementation and Benchmarking

The efficient implementation of the third-order unitary coupled-cluster scheme (UCC3) for the calculation of excited electronic states is reported. The UCC3 scheme and its second-order UCC2 variant have been benchmarked and compared to Jacquemin's recently introduced, as well as Thiel's well-established, benchmark sets for excitation energies and oscillator strengths. For the latter, the calculation of 134 excited singlet and 71 excited triplet states of 28 small- to medium-sized organic molecules has revealed that UCC2 exhibits a mean error and standard deviation of 0.36 ± 0.41 eV for singlet states and 0.22 ± 0.21 eV for triplet states, whereas UCC3 revealed an accuracy of 0.06 ± 0.27 eV for singlet and -0.22 ± 0.15 eV for triplet states. In addition, the oscillator strengths obtained with effective transition moments correct through second order in perturbation theory are in very good agreement with literature data.

PSIXAS: A Psi4 plugin for efficient simulations of X-ray absorption spectra based on the transition-potential and Δ-Kohn-Sham method

Near edge X-ray absorption fine structure (NEXAFS) spectra and their pump-probe extension (PP-NEXAFS) offer insights into valence- and core-excited states. We present PSIXAS, a recent implementation for simulating NEXAFS and PP-NEXAFS spectra by means of the transition-potential and the Δ-Kohn-Sham method. The approach is implemented in form of a software plugin for the Psi4 code, which provides access to a wide selection of basis sets as well as density functionals. We briefly outline the theoretical foundation and the key aspects of the plugin. Then, we use the plugin to simulate PP-NEXAFS spectra of thymine, a system already investigated by others and us. It is found that larger, extended basis sets are needed to obtain more accurate absolute resonance positions. We further demonstrate that, in contrast to ordinary NEXAFS simulations, where the choice of the density functional plays a minor role for the shape of the spectrum, for PP-NEXAFS simulations the choice of the density functional is important. Especially hybrid functionals (which could not be used straightforwardly before to simulate PP-NEXAFS spectra) and their amount of "Hartree-Fock like" exact exchange affects relative resonance positions in the spectrum.

Generalized single excitation configuration interaction: an investigation into the impact of the inclusion of non-orthogonality on the calculation of core-excited states

In this paper, we investigate different non-orthogonal generalizations of the configuration interaction with single substitutions (CIS) method and their impact on the calculation of core-excited states.

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.

Equation of motion coupled-cluster for core excitation spectra: Two complementary approaches

This paper presents core excitation spectra from coupled-cluster (CC) theory obtained from both a time-independent and a new time-dependent formalism. The conventional time-independent CC formulation for excited states is the equation-of-motion (EOM-CC) method whose eigenvalues and eigenvectors describe the core excited states. An alternative computational procedure is offered by a time-dependent CC description. In that case, the dipole transition operator is expressed in the CC effective Hamiltonian form and propagated with respect to time. The absorption spectrum is obtained from the CC dipole autocorrelation function via a Fourier transformation. Comparisons are made among the time-dependent results obtained from second-order perturbation theory, to coupled cluster doubles and their linearized forms (CCD and LCCD), to CC singles and doubles (CCSD) and the linearized form (LCCSD). In the time-independent case, considerations of triples (EOM-CCSDT) and quadruples (EOM-CCSDTQ) are used to approach sub-electron volt accuracy. A particular target is the allyl radical, as an example of an open-shell molecule. As the results have to ultimately be the same, the two procedures offer a complementary approach toward analyzing experimental results.

Approximate Green's Function Coupled Cluster Method Employing Effective Dimension Reduction

The Green's function coupled cluster (GFCC) method, originally proposed in the early 1990s, is a powerful many-body tool for computing and analyzing the electronic structure of molecular and periodic systems, especially when electrons of the system are strongly correlated. However, in order for the GFCC to become a method that may be routinely used in the electronic structure calculations, robust numerical techniques and approximations must be employed to reduce its extremely high computational overhead. In our recent studies, it has been demonstrated that the GFCC equations can be solved directly in the frequency domain using iterative linear solvers, which can be easily distributed in a massively parallel environment. In the present work, we demonstrate a successful application of model-order-reduction (MOR) techniques in the GFCC framework. Briefly speaking, for a frequency regime of interest that requires high-resolution descriptions of spectral function, instead of solving the GFCC linear equation of full dimension for every single frequency point of interest, an efficiently solvable linear system model of a reduced dimension may be built upon projecting the original GFCC linear system onto a subspace. From this reduced order model is obtained a reasonable approximation to the full dimensional GFCC linear equations in both interpolative and extrapolative spectral regions. Here, we show that the subspace can be properly constructed in an iterative manner from the auxiliary vectors of the GFCC linear equations at some selected frequencies within the spectral region of interest. During the iterations, the quality of the subspace, as well as the linear system model, can be systematically improved. The method is tested in this work in terms of the efficiency and accuracy of computing spectral functions for some typical molecular systems such as carbon monoxide, 1,3-butadiene, benzene, and adenine. To reach the same level of accuracy as that of the original GFCC method, the application of MOR in the GFCC method is able to significantly lower the original computational cost for the aforementioned molecules in designated frequency regimes. As a byproduct, the reduced order model obtained by this method is found to provide a high-quality initial guess, which improves the convergence rate for the existing iterative linear solver.

Algebraic diagrammatic construction for the polarisation propagator in combination with effective fragment potentials

The effective fragment potential (EFP) method for the efficient inclusion of solvation effects is combined with the algebraic diagrammatic construction (ADC) scheme for the second- and third-order polarisation propagator.

Modeling core-level excitations with variationally optimized reduced-density matrices and the extended random phase approximation

The information contained within ground-state one- and two-electron reduced-density matrices (RDMs) can be used to compute wave functions and energies for electronically excited states through the extended random phase approximation (ERPA). The ERPA is an appealing framework for describing excitations out of states obtained via the variational optimization of the two-electron RDM (2-RDM), as the variational 2-RDM (v2RDM) approach itself can only be used to describe the lowest-energy state of a given spin symmetry. The utility of the ERPA for predicting near-edge features relevant to x-ray absorption spectroscopy is assessed for the case that the 2-RDM is obtained from a ground-state v2RDM-driven complete active space self-consistent field (CASSCF) computation. A class of killer conditions for the CASSCF-specific ERPA excitation operator is derived, and it is demonstrated that a reliable description of core-level excitations requires an excitation operator that fulfills these conditions; the core-valence separation (CVS) scheme yields such an operator. Absolute excitation energies evaluated within the CASSCF/CVS-ERPA framework are slightly more accurate than those obtained from the usual random phase approximation (RPA), but the CVS-ERPA is not more accurate than RPA for predicting the relative positions of near-edge features. Nonetheless, CVS-ERPA is established as a reasonable starting point for the treatment of core-level excitations using variationally optimized 2-RDMs.

Green's function coupled cluster formulations utilizing extended inner excitations

In this paper, we analyze new approximations of the Green's function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i(n, m) method, where the excitation level of the inner auxiliary operators (m) used to describe the ionization potential and electron affinity effects in the N - 1 and N + 1 particle spaces is higher than the excitation level (n) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called "n + 1" rule in this category [or the GFCC-i(n, n + 1) method], which states that in order to maintain size-extensivity of the Green's function matrix elements, the excitation level of inner auxiliary operators X _p (ω) and Y _q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N₂ and CO molecules in the inner and outer valence regimes. In comparison with the Green's function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.

The All-Seeing Eye of Resonant Auger Electron Spectroscopy: A Study on Aqueous Solution Using Tender X-rays

X-ray absorption and Auger electron spectroscopies are demonstrated to be powerful tools to unravel the electronic structure of solvated ions. In this work for the first time, we use a combination of these methods in the tender X-ray regime. This allowed us to address electronic transitions from deep core levels, to probe environmental effects, specifically in the bulk of the solution since the created energetic Auger electrons possess large mean free paths, and moreover, to obtain dynamical information about the ultrafast delocalization of the core-excited electron. In the considered exemplary aqueous KCl solution, the solvated isoelectronic K⁺ and Cl^- ions exhibit notably different Auger electron spectra as a function of the photon energy. Differences appear due to dipole-forbidden transitions in aqueous K⁺ whose occurrence, according to the performed ab initio calculations, becomes possible only in the presence of solvent water molecules.

Births in the United States, 2017

This report presents selected highlights from 2017 final birth data on key demographic, health care utilization, and infant health indicators. General fertility rates (the number of births per 1,000 females aged 15-44 years) and teen birth rates are presented by race and Hispanic origin. The use of Medicaid as the source of payment for the delivery and preterm birth rates are presented by the age of the mother. Data for 2017 are compared with 2016 for each indicator.

Model Order Reduction Algorithm for Estimating the Absorption Spectrum

The ab initio description of the spectral interior of the absorption spectrum poses both a theoretical and computational challenge for modern electronic structure theory. Due to the often spectrally dense character of this domain in the quantum propagator's eigenspectrum for medium-to-large sized systems, traditional approaches based on the partial diagonalization of the propagator often encounter oscillatory and stagnating convergence. Electronic structure methods which solve the molecular response problem through the solution of spectrally shifted linear systems, such as the complex polarization propagator, offer an alternative approach which is agnostic to the underlying spectral density or domain location. This generality comes at a seemingly high computational cost associated with solving a large linear system for each spectral shift in some discretization of the spectral domain of interest. In this work, we present a novel, adaptive solution to this high computational overhead based on model order reduction techniques via interpolation. Model order reduction reduces the computational complexity of mathematical models and is ubiquitous in the simulation of dynamical systems and control theory. The efficiency and effectiveness of the proposed algorithm in the ab initio prediction of X-ray absorption spectra is demonstrated using a test set of challenging water clusters which are spectrally dense in the neighborhood of the oxygen K-edge. On the basis of a single, user defined tolerance we automatically determine the order of the reduced models and approximate the absorption spectrum up to the given tolerance. We also illustrate that, for the systems studied, the automatically determined model order increases logarithmically with the problem dimension, compared to a linear increase of the number of eigenvalues within the energy window. Furthermore, we observed that the computational cost of the proposed algorithm only scales quadratically with respect to the problem dimension.