skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: 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 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 nonlocal core orbitals, and the consequences of breaking symmetry are discussed.  more » « less
Award ID(s):
1848580 2019144
PAR ID:
10420590
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
157
Issue:
8
ISSN:
0021-9606
Page Range / eLocation ID:
084115
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Modeling L-edge spectra at X-ray wavelengths requires consideration of spin–orbit splitting of the 2p orbitals. We introduce a low-cost tool to compute core-level spectra that combines a spin–orbit mean-field description of the Breit–Pauli Hamiltonian with nonrelativistic excited states computed using the semi-empirical density-functional theory configuration-interaction singles (DFT/CIS) method, within the state-interaction approach. Our version of DFT/CIS was introduced recently for K-edge spectra and includes a semi-empirical correction to the core orbital energies, significantly reducing ad hoc shifts that are typically required when time-dependent (TD-)DFT is applied to core-level excitations. In combination with the core/valence separation approximation and spin–orbit couplings, the DFT/CIS method affords semiquantitative L-edge spectra at CIS cost. Spin–orbit coupling has a qualitative effect on the spectra, as demonstrated for a variety of 3d transition metal systems and main-group compounds. The use of different active orbital spaces helps to facilitate spectral assignments. We find that spin–orbit splitting has a negligible effect on M-edge spectra for 3d transition metal species. 
    more » « less
  2. A practical ab initio composite method for modeling x-ray absorption and non-resonant x-ray emission is presented. Vertical K-edge excitation and emission energies are obtained from core-electron binding energies calculated with spin-projected ΔHF/ΔMP and outer-core ionization potentials/electron affinities calculated with electron propagator theory. An assessment of the combined methodologies against experiment is performed for a set of small molecules containing second-row elements. 
    more » « less
  3. Reliable computational methodologies and basis sets for modeling x-ray spectra are essential for extracting and interpreting electronic and structural information from experimental x-ray spectra. In particular, the trade-off between numerical accuracy and computational cost due to the size of the basis set is a major challenge, since molecular orbitals undergo extreme relaxation in the core-hole state. To gain clarity on the changes in electronic structure induced by the formation of a core-hole, the use of sufficiently flexible basis for expanding the orbitals, particularly for the core region, has been shown to be essential. This work focuses on the refinement of core-hole ionized state calculations using the equation-of-motion coupled cluster family of methods through an extensive analysis on the effectiveness of “hybrid” and mixed basis sets. In this investigation, we utilize the CVS-EOMIP-CCSD method in combination and construct hybrid basis sets piecewise from readily available Dunning’s correlation consistent basis sets in order to calculate x-ray ionization energies (IEs) for a set of small gas phase molecules. Our results provide insights into the impact of basis sets on the CVS-EOMIP-CCSD calculations of K-edge IEs of first-row p-block elements. These insights enable us to understand more about the basis set dependence of the core IEs computed and allow us to establish a protocol for deriving reliable and cost-effective theoretical estimates for computing IEs of small molecules containing such elements. 
    more » « less
  4. Strong field ionization of molecules by intense laser pulses can be simulated by time-dependent configuration interaction (TD-CI) with a complex absorbing potential (CAP). Standard molecular basis sets need to be augmented with several sets of diffuse functions for effective interaction with the CAP. This dramatically increases the number of configurations and the cost of the TD-CI simulations as the size of the molecules increases. The cost can be reduced by making use of spin symmetry and by employing an orbital energy cut-off to limit the number of virtual orbitals used to construct the excited configurations. Greater reductions in the number of virtual orbitals can be obtained by examining their interaction with the absorbing potential during simulations and their contributions to the strong field ionization rate. This can be determined from the matrix elements of the absorbing potential and the TD-CI coefficients from test simulations. Compared to a simple 3 hartree cut-off in the orbital energies, these approaches reduce the number of virtual orbitals by 20% - 35% for neutral molecules and 5%-10% for cations. As a result, the cost of simulations is reduced by 35% - 60% for neutral molecules and 5% - 10% for cations. The number of virtual orbitals needed can also be estimated by second-order perturbation theory without the need for test simulations. The number of virtual orbitals can be reduced further by adapting orbitals to the laser field using natural orbitals derived from test simulations. This is particularly effective for cations, yielding reductions of more than 20%. 
    more » « less
  5. 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 of X-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. 
    more » « less