Search

Table-top x-ray spectroscopy measurements at the carbon K-edge complemented by ab initio calculations are used to investigate the influence of the bromine atom on the carbon core-valence transitions in bromobenzene cation (BrBz+). The electronic ground state of the cation is prepared by resonance-enhanced two-photon ionization of neutral bromobenzene (BrBz) and probed by x-rays produced by high-harmonic generation (HHG). Replacing one of the hydrogen atoms in benzene with a bromine atom shifts the transition from the 1sC* orbital of the carbon atom (C*) bonded to bromine by ~1 eV to higher energy in the x-ray spectrum compared to the other carbon atoms (C). The x-ray absorption spectra of neutral BrBz reveals the influence of the bromine atom on the x-ray absorption spectrum due to the 1sC* orbital, while the influence of the bromine atom on the cation spectrum is even more prominent. The x-ray spectrum of the cation is dominated by two relatively intense transitions, the 1sC→ π* and the 1sC*→ σ*(C*-Br), where the second transition is enhanced relative to the neutral BrBz. In addition, a doublet peak shape for these two transitions is observed in the experiment. The 1sC→ π* doublet peak shape arises due to the spin coupling of the unpaired electron in the partially vacant π orbital (from ionization) with the two other unpaired electrons resulting from the transition from the 1sC core orbital to the fully vacant π* orbitals. The 1sC*→ σ* doublet peak shape results from several transitions involving σ* and vibrational C*-Br mode activation following the UV ionization, which demonstrates the impact of the C*-Br bond length on the core-valence transition as well as the relaxation geometry of BrBz+.

We present a new computational framework to describe polaritons, which treats photons and electrons on the same footing using coupled-cluster theory. As a proof of concept, we study the coupling between the first electronically excited state of carbon monoxide and an optical cavity. In particular, we focus on how the interaction with the photonic mode changes the vibrational spectroscopic signature of the electronic state and how this is affected when tuning the cavity frequency and the light–matter coupling strength. For this purpose, we consider different methodologies and investigate the validity of the Born–Oppenheimer approximation in such situations.

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. 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.

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- (CH3I, HX, and X-, X = Cl-At) and xenon-containing (Xe, XeF2) 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 (2DCM), 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 (2DCGM) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.

The recent development of liquid jet and liquid leaf sample delivery systems allows for accurate measurements of soft X-ray absorption spectra in transmission mode of solutes in a liquid environment. As this type of measurement becomes increasingly accessible, there is a strong need for reliable theoretical methods for assisting in the interpretation of the experimental data. Coupled cluster methods have been extensively developed over the past decade to simulate X-ray absorption in the gas phase. Their performance for solvated species, on the contrary, remains largely unexplored. Here, we investigate the current state of the art of coupled cluster modeling of nitrogen K-edge X-ray absorption of aqueous ammonia and ammonium based on quantum mechanics/molecular mechanics, where both the level of coupled cluster calculations and polarizable embedding are scrutinized. The results are compared to existing experimental data as well as simulations based on transition potential density functional theory.

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.

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. Ahttp://dx.doi.org/10.1021/acs.jpca.0c08736]. We show that the 1sC → π 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 1sC → π* 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 1sC → π* 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.

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 as well as the experimental measurement ofthe carbon K-edge spectra of both species using a table-top high-harmonic generation (HHG) source have been described in a companion paper [M. Epshtein et al., J. Am.Chem. Soc., submitted. Available on ChemRxiv]. We show that the 1sC -> pi 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 presenceof the unpaired (spectator) electron in the pi-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 1sC ->pi* 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 1sC -> pi* band of the cation is attributed to the interplay between the coupling of the core -> pi* excitation with the unpaired electronin the pi-subshell and the Jahn-Teller distortion. The calculations attribute most ofthe splitting (~1-1.2 eV) to the spin coupling, which is visible already at the Franck-Condon structure, and estimate the additional splitting due to structural relaxation tobe 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-Tellereffect in benzene cation.

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 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 argon atoms to small molecules containing sulfur and silicon.

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, submitted to J. Phys. Chem. Lett.; ChemRxiv link: doi XXXXX], 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 1sC core orbital to the fully vacant π* orbitals.

Physical review letters 124(23), 236001 (2020). doi:10.1103/PhysRevLett.124.236001, 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., Published by APS, College Park, Md.

DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator theory. At the self-consistent-field level a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding, and frozen density embedding models. DIRAC was one of the earliest codes for relativistic molecular calculations and remains a reference in its field., The following article has been submitted to The Journal of Chemical Physics. After it is published, it will be found at https://aip.scitation.org/toc/jcp/current

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree–Fock, Kohn–Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

Correction for ‘Dyson Orbitals within the fc-CVS-EOM-CCSD framework: theory and application to X-ray photoelectron spectroscopy of ground and excited states’ by M. L. Vidal et al., Phys. Chem. Chem. Phys., 2020, DOI: 10.1039/C9CP03695D

We report on the implementation of Dyson orbitals within the recently introduced frozen-core (fc) core-valence separated (CVS) equation-of-motion (EOM) coupled-cluster singles and doubles (CCSD) method, which enables efficient and reliable characterization of core-level states. The ionization potential (IP) variant of fc-CVS-EOM-CCSD, in which the EOM target states have one electron less than the reference, gives access to core-ionized states thus enabling modeling ofX-ray photoelectron spectra (XPS) and its time-resolved variant (TR-XPS). Dyson orbitals are reduced quantities that can be interpreted as correlated states of the ejected/attached electron; they enter the expressions of various experimentally relevant quantities. In the context of photoelectron spectroscopy, Dyson orbitals can be used to estimate the strengths of photoionization transitions. We illustrate the utility of Dyson orbitals and fc-CVS-EOM-IP-CCSD by calculating XPS of the ground state of adenine and TR-XPS of the excited states of uracil.

We present a novel approach for computing resonant inelastic X-ray scattering (RIXS) cross sections within the equation-of-motion coupled-cluster (EOM-CC) framework. The approach is based on recasting the sum-over-states expressions for RIXS moments into closed-form expressions by using damped response theory. Damped response formalism allows one to circumvent problems of divergent behavior of response equations in the resonant regime. However, the convergence of response equations in the X-ray frequency range is often erratic due to the electronically metastable (i.e., resonant) nature of the virtual core-excited states embedded in the valence ionization continuum. We circumvent this problematic behavior by extending the core–valence separation (CVS) scheme, which decouples the valence-excited and core-excited configurations of the excitation manifold, into the response domain. The accuracy of the CVS-enabled damped response theory, implemented within the EOM-EE-CCSD (EOM-CC for excitation energies with single and double excitations) framework, is assessed by comparison against standard damped EOM-EE-CCSD response calculations. The capabilities of the new approach are illustrated by calculations of RIXS cross sections for benzene and benzene radical cation.

We report on the implementation and illustrative applications of Dyson orbitals within the recently proposed frozen-core (fc) core-valence separated (CVS) equation-of-motion (EOM) coupled-cluster singles and doubles (CCSD) method, which enables efficient and accurate characterization of core-ionized states. Dyson orbitals are reduced quantities that can be interpreted as correlated states of the ejected/attached electron.Dyson orbitals enter the expressions of various experimental observables, such as photoionization cross sections; thus, they are necessary for modeling photoelectron spectra.Here we discuss the simulations of X-ray photoelectron spectra (XPS) and propose an approach to simulate time-resolved (TR-)XPS for probing excited states. As illustrative examples, we present the simulation of the XPS of the ground state of adenine and of TR-XPS of the excited states of uracil.