Search

Theoretical optical and x-ray spectra of model structures of water and ice are calculated using a many-body perturbation theory, Bethe-Salpeter equation (BSE) approach implemented in the valence- and core-excitation codes AI2NBSE and OCEAN. These codes use ab initio density functional theory wave functions from a plane-wave, pseudopotential code, quasiparticle self-energy corrections, and a BSE treatment of particle-hole interactions. This approach improves upon independent-particle methods through the inclusion of a complex, energy-dependent self-energy and screened particle-hole interactions to account for inelastic losses and excitonic effects. These many-body effects are found to be crucial for quantitative calculations of ice and water spectra. [ABSTRACT FROM AUTHOR]

We present results for realistic modeling of the x-ray absorption near edge structure (XANES) of the overdoped high-Tc, superconductor La2-xSrxCuO4 in the hole doping range x = 0.20 - 0.30. Our computations are based on a real-space Green's function approach in which strong-correlation effects are taken into account in terms of a doping-dependent self-energy. The predicted O K-edge XANES is found to be in good accord with the corresponding experimental results in this overdoped regime. We find that the low energy spectra are dominated by the contribution of O atoms in the cuprate planes, with little contribution from apical O atoms. [ABSTRACT FROM AUTHOR]

We present ab initio calculations of the photoemission spectra of bulk sodium using different flavors of the cumulant expansion approximation for the Green's function. In particular, we study the dispersion and intensity of the plasmon satellites. We show that the satellite spectrum is much more sensitive to many details than the quasiparticle spectrum, which suggests that the experimental investigation of satellites could yield additional information beyond the usual studies of the band structure. In particular, a comparison to the homogeneous electron gas shows that the satellites are influenced by the crystal environment, although the crystal potential in sodium is weak. Moreover, the temperature dependence of the lattice constant is reflected in the position of the satellites. Details of the screening also play an important role; in particular, the contribution of transitions from 2 s and 2 p semicore levels influences the satellites, but not the quasiparticle. Moreover, inclusion of contributions to the screening beyond the random-phase approximation has an effect on the satellites. Finally, we elucidate the importance of the coupling of electrons and holes by comparing the results of the time-ordered and the retarded cumulant expansion approximations. Again, we find small but noticeable differences. Since all the small effects add up, our most advanced calculation yields a satellite position which is improved with respect to previous calculations by almost 1 eV. This stresses the fact that the calculation of satellites is much more delicate than the calculation of a quasiparticle band structure. [ABSTRACT FROM AUTHOR]

The electronic structure of stoichiometric tin dioxide (SnO2) is studied by probing its unoccupied states using the fine structure in the electron energy-loss spectra (EELS) at the oxygen-K (O-K) edge. The spectral measurements were performed both at room and at high temperatures (773 K) and compared to ab initio calculations carried out using the real-space multiple-scattering and linearized augmented-plane-wave methods. Important many-body effects are included via quasiparticle corrections calculated within the many-pole GW self-energy approximation. An additional energy-dependent damping is calculated to account for vibrational effects. Results from this paper demonstrated that quantitative agreement between theoretical and experimental spectra can be obtained when nonspherical potentials and quasiparticle self-energy effects are considered and vibrational broadening is included. Modifications of the electronic structure by single oxygen vacancies, both in the bulk and at the (110) surface, also are predicted. Our predictions support the use of O-K EELS as a probe of the defect structures in SnO2 surfaces and nanoparticles. [ABSTRACT FROM AUTHOR]

We present calculated valence and C 1s near-edge excitation spectra of solid C60 and experimental results measured with high-resolution electron energy-loss spectroscopy. The near-edge calculations are carried out using three different methods: solution of the Bethe-Salpeter equation (BSE) as implemented in the ocean suite (Obtaining Core Excitations with Ab Initio methods and the NIST BSE solver), the excited-electron core-hole approach, and the constrained-occupancy method using the Stockholm-Berlin core excitation code, StoBe. The three methods give similar results and are in good agreement with experiment, though the BSE results are the most accurate. The BSE formalism is also used to carry out valence level calculations using the NIST BSE solver. Theoretical results include self-energy corrections to the band gap and bandwidths, lifetime-damping effects, and Debye-Waller effects in the core excitation case. A comparison of spectral features to those observed experimentally illustrates the sensitivity of certain features to computational details, such as self-energy corrections to the band structure and core-hole screening. [ABSTRACT FROM AUTHOR]

High resolution inelastic x-ray scattering experiments are widely used to study the electronic and chemical properties of materials under a range of conditions, from ambient temperature to the warm dense matter regime. We use the real-space multiple scattering (RSMS) Green's function formalism coupled with density functional theory molecular dynamics (DFT-MD) to study thermal effects on the Compton profile (CP) of disordered systems. The RSMS method is advantageous for calculations of highly disordered, aperiodic systems because it places no restriction on symmetry. As a test, we apply our approach to thermally disordered Be, Li, and Si in both liquid and solid phases. We find good agreement with experimental and other theoretical results, showing that the real-space multiple scattering approach coupled with DFT-MD is an efficient and reliable method for calculating the CP of disordered systems at finite temperatures. [ABSTRACT FROM AUTHOR]

The treatment of electronic correlations in open-shell systems is among the most challenging problems of condensed matter theory. Current approximations are only partly successful. Ligand-field multiplet theory has been widely successful in describing intra-atomic correlation effects in x-ray spectra, but typically ignores itinerant states. The cumulant expansion for the one-electron Green's function has been successful in describing shake-up effects but ignores atomic multiplets. More complete methods, such as dynamic mean-field theory can be computationally demanding. Here, we show that separating the dynamic Coulomb interactions into local and longer-range parts with ab initio parameters yields a combined multiplet-plus-cumulant approach that accounts for both local atomic multiplets and satellite excitations. The approach is illustrated in transition metal oxides and explains the multiplet peaks, charge-transfer satellites, and distributed background features observed in XPS experiment.

Many-body excitations in X-ray photoemission spectra have been difficult to simulate from first principles. We have recently developed a cumulant-based one-electron Green's function method using the real-time coupled-cluster-singles equation-of-motion approach (RT-EOM-CCS) that provides a general framework for treating these problems. Here we extend this approach to include double excitations in the ground-state energy and in the coupled cluster amplitudes, which have been implemented using subroutines generated by the Tensor Contraction Engine (TCE). As in the case of the singles approximation, RT-EOM-CCSD yields a nonperturbative cumulant form of the Green's function in terms of the time-dependent cluster amplitudes, adding nonlinear corrections to the traditional cumulant forms. The extended approach is applied to the core-hole spectral function for small molecular systems. We find that, when core-optimized basis sets are used, the doubles contributions reduce the mean absolute errors in the core binding energies of the 10 e systems from 0.8 to 0.3 eV. They also significantly improve the quasiparticle-satellite gap by reducing its overestimation from about 3-5 to about 0-1 eV in CH 4 , NH 3 , and H 2 O, and also improving the overall shape of the satellite features. Finally, we demonstrate the application of the new implementation to the larger, classical XPS ESCA series of molecules and show that the singles approximation can be paired with a modest basis set to study carbon speciation.

Green's function methods provide a robust, general framework within many-body theory for treating electron correlation in both excited states and x-ray spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green's function approach has been introduced, where the cumulant is obtained as the solution to a system of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green's function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function for the core Green's function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then calculated using a convolution of the core spectral function and an effective, one-body XAS. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green's function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. The method is illustrated for the XPS and XAS of NH 3 ., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Vila, Kas, Rehr, Kowalski and Peng.)

Interaction effects can change materials properties in intriguing ways, and they have, in general, a huge impact on electronic spectra. In particular, satellites in photoemission spectra are pure many-body effects, and their study is of increasing interest in both experiment and theory. However, the intrinsic spectral function is only a part of a measured spectrum, and it is notoriously difficult to extract this information, even for simple metals. Our joint experimental and theoretical study of the prototypical simple metal aluminum demonstrates how intrinsic satellite spectra can be extracted from measured data using angular resolution in photoemission. A nondispersing satellite is detected and explained by electron-electron interactions and the thermal motion of the atoms. Additional nondispersing intensity comes from the inelastic scattering of the outgoing photoelectron. The ideal intrinsic spectral function, instead, has satellites that disperse both in energy and in shape. Theory and the information extracted from experiment describe these features with very good agreement., Competing Interests: The authors declare no competing interest.

Green's function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Here, we develop the cumulant form of the one-electron Green's function using a real-time coupled-cluster equation-of-motion approach, in an extension of our previous study (Rehr J.; et al. J. Chem. Phys. 2020, 152, 174113). The approach yields a nonperturbative expression for the cumulant in terms of the solution to a set of coupled first-order, nonlinear differential equations. The method thereby adds nonlinear corrections to traditional cumulant methods, which are linear in the self-energy. The approach is applied to the core-hole Green's function and is illustrated for a number of small molecular systems. For these systems, we find that the nonlinear contributions yield significant improvements, both for quasiparticle properties such as core-level binding energies and for inelastic losses that correspond to satellites observed in photoemission spectra.

Electron-core hole interactions are critical for proper interpretation of core-level spectroscopies commonly used as analytical tools in materials science. Here we utilize resonant Auger-electron spectroscopy to uniquely identify exciton, shake, and charge-transfer processes that result from the sudden creation of the core hole in both x-ray-absorption and photoemission spectra. These effects are captured for the transition-metal compounds SrTiO 3 and MoS 2 by fully ab initio , combined real-time cumulant, and Bethe-Salpeter equation approaches to account for core hole dynamics and screening. Atomic charges and excited-state electron-density fluctuations reflect materials' solid-state electronic structure, loss of translational symmetry around the core hole, and breakdown of the sudden approximation. They also demonstrate competition between long- and short-range screening in a solid.

First-principles, real-time-cumulant, and Bethe-Salpeter-equation calculations fully capture the detailed satellite structure that occurs in response to the sudden creation of the core hole in both photoemission and x-ray absorption spectra of the transition-metal compounds SrTiO 3 and rutile TiO 2 . Analysis of the excited-state, real-space charge-density fluctuations betrays the physical nature of these many electron excitations that are shown to reflect the materials' solid-state electronic structure and chemical bonding. This first-principles development of the cumulant-based core hole spectral function is generally applicable to other systems and should become a standard tool for all similar spectroscopic analysis going beyond the quasiparticle physics of the photoelectric effect.

We present a combined equation of motion coupled-cluster cumulant Green's function approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS) and x-ray photoemission spectra. The method is based on a factorization of the transition amplitude in the time domain, which leads to a convolution of an effective one-body absorption spectrum and the core-hole spectral function. The spectral function characterizes intrinsic losses in terms of shake-up excitations and satellites using a cumulant representation of the core-hole Green's function that simplifies the interpretation. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.

We present calculated valence and C 1s near-edge excitation spectra of solid C 60 and experimental results measured with high-resolution electron energy-loss spectroscopy. The near-edge calculations are carried out using three different methods: solution of the Bethe-Salpeter equation (BSE) as implemented in the OCEAN suite (Obtaining Core Excitations with ab initio methods and the NIST BSE solver), the excited-electron core-hole approach (XCH), and the constrained-occupancy method using the Stockholm-Berlin core-excitation code, StoBe. The three methods give similar results and are in good agreement with experiment, though the BSE results are the most accurate. The BSE formalism is also used to carry out valence level calculations using the NIST Bethe-Salpeter Equation solver (NBSE). Theoretical results include self-energy corrections to the band gap and band widths, lifetime-damping effects, and Debye-Waller effects in the core-excitation case. A comparison of spectral features to those observed experimentally illustrates the sensitivity of certain features to computational details, such as self-energy corrections to the band structure and core-hole screening.

One of the big challenges of theoretical condensed-matter physics is the description, understanding, and prediction of the effects of the Coulomb interaction on materials properties. In electronic spectra, the Coulomb interaction causes a renormalization of energies and change of spectral weight. Most importantly, it can lead to new structures, often called satellites. These can be linked to the coupling of excitations, also termed dynamical effects. State-of-the-art methods in the framework of many-body perturbation theory, in particular, the widely used GW approximation, often fail to describe satellite spectra. Instead, approaches based on a picture of electron-boson coupling such as the cumulant expansion are promising for the description of plasmon satellites. In this work, we give a unified derivation of the GW approximation and the cumulant expansion for the one-body Green's function. Using the example of bulk sodium, we compare the resulting spectral functions both in the valence and in the core region, and we discuss the dispersion of quasi-particles and satellites. We show that self-consistency is crucial to obtain meaningful results, in particular, at large binding energies. Very good agreement with experiment is obtained when the intrinsic spectral function is corrected for extrinsic and interference effects. Finally, we sketch how one can approach the problem in the case of the two-body Green's function, and we discuss the cancellation of various dynamical effects that occur in that case.

We show how the classic approach of Nozières and di Domenicis for treating the edge singularity in x-ray absorption and emission can be generalized to treat the more complex case of the resonant inelastic x-ray scattering (RIXS) process, including effects of the intermediate states involved therein in the presence of the core hole. We solve our lattice model essentially exactly (numerically) to obtain a novel form of edge singularity at the RIXS threshold energy. Our RIXS spectrum naturally includes both the well and poorly screened spectral components and their dispersions and allows its separation into pair and multiple-pair excitations.

We investigate the origin and physical effects of "dynamic structural disorder" (DSD) in supported nano-scale catalysts. DSD refers to the intrinsic fluctuating, inhomogeneous structure of such nano-scale systems. In contrast to bulk materials, nano-scale systems exhibit substantial fluctuations in structure, charge, temperature, and other quantities, as well as large surface effects. The DSD is driven largely by the stochastic librational motion of the center of mass and fluxional bonding at the nanoparticle surface due to thermal coupling with the substrate. Our approach for calculating and understanding DSD is based on a combination of real-time density functional theory/molecular dynamics simulations, transient coupled-oscillator models, and statistical mechanics. This approach treats thermal and dynamic effects over multiple time-scales, and includes bond-stretching and -bending vibrations, and transient tethering to the substrate at longer ps time-scales. Potential effects on the catalytic properties of these clusters are briefly explored. Model calculations of molecule-cluster interactions and molecular dissociation reaction paths are presented in which the reactant molecules are adsorbed on the surface of dynamically sampled clusters. This model suggests that DSD can affect both the prefactors and distribution of energy barriers in reaction rates, and thus can significantly affect catalytic activity at the nano-scale.

We present a hybrid approach for Bethe-Salpeter equation (BSE) calculations of core excitation spectra, including x-ray absorption (XAS), electron energy loss spectra (EELS), and nonresonant inelastic x-ray scattering (NRIXS). The method is based on ab initio wave functions from the plane-wave pseudopotential code ABINIT; atomic core-level states and projector augmented wave (PAW) transition matrix elements; the NIST core-level BSE solver; and a many-pole self-energy model to account for final-state broadening and self-energy shifts. Multiplet effects are also approximately accounted for. The approach is implemented using an interface dubbed OCEAN (Obtaining Core Excitations using ABINIT and NBSE). To demonstrate the utility of the code we present results for the K edges in LiF as probed by XAS and NRIXS, the K edges of KCl as probed by XAS, the Ti L 2,3 edge in SrTiO 3 as probed by XAS, and the Mg L 2,3 edge in MgO as probed by XAS. These results are compared with experiment and with other theoretical approaches.

New theoretical and experimental investigations of the occupied and unoccupied local electronic densities of states (DOS) are reported for alpha-Li(3)N. Band-structure and density-functional theory calculations confirm the absence of covalent bonding character. However, real-space full-multiple-scattering (RSFMS) calculations of the occupied local DOS find less extreme nominal valences than have previously been proposed. Nonresonant inelastic x-ray scattering, RSFMS calculations, and calculations based on the Bethe-Salpeter equation are used to characterize the unoccupied electronic final states local to both the Li and N sites. There is a good agreement between experiment and theory. Throughout the Li 1s near-edge region, both experiment and theory find strong similarities in the s-and p-type components of the unoccupied local final DOS projected onto an orbital angular momentum basis (l-DOS). An unexpected, significant correspondence exists between the near-edge spectra for the Li 1s and N 1s initial states. We argue that both spectra are sampling essentially the same final DOS due to the combination of long core-hole lifetimes, long photoelectron lifetimes, and the fact that orbital angular momentum is the same for all relevant initial states. Such considerations may be generally applicable for low atomic number compounds.

We discuss the application of the ab initio program FEFF8 to calculations of electron energy-loss spectroscopy (EELS), focusing in particular on core-loss spectra. FEFF8 is based on a self-consistent, real space multiple scattering formalism. We focus on issues relevant to practical simulations, including the construction of well-converged potentials, the treatment of inelastic losses and exchange-correlation potentials and the core-hole. We also discuss how to account for experimental conditions, for example, sample orientation and finite temperature effects such as Debye-Waller factors. Finally we discuss the interpretation of the spectra in terms of electronic structure and local projected density of states (LDOS). As an explicit example, we illustrate various features of the code by application to the ionization edges of GaN.

There has been dramatic progress over the past decade both in theory and in ab initio calculations of X-ray absorption fine structure. Significant progress has also been made in understanding X-ray absorption near-edge structure (XANES). This contribution briefly reviews the developments in this field leading up to the current state. One of the key advances has been the development of several ab initio codes such as FEFF, which permit an interpretation of the spectra in terms of geometrical and electronic properties of a material. Despite this progress, XANES calculations have remained challenging both to compute and to interpret. However, recent advances based on parallel Lanczos multiple-scattering algorithms have led to speed increases of typically two orders of magnitude, making fast calculations practicable. Improvements in the interpretation of near-edge structure have also been made. It is suggested that these developments can be advantageous in structural biology, e.g. in post-genomics studies of metalloproteins.

A quantitative analysis of the XAS spectra of the tetracyanonickelate complex [Ni(CN)4]2- has been carried out. The simultaneous study of the EXAFS and XANES regions yielded complementary information regarding the geometric and electronic structures of the complex. XANES spectra were modeled by applying recently developed self-consistent, full multiple-scattering algorithms in the FEFF8 code (version 8 x 34). XANES spectra for clusters of different sizes (from 9 to 125 atoms) were computed and compared with experimental spectra. This region of the spectra was proportional to a broadened Ni p-density of states diagram above the Fermi level. Although the main features of the XANES spectra were reasonably reproduced by computations, the weak dependence of the theoretical spectra on cluster size contrasts with the close similarity between the experimental spectra of the solid and solution systems. Because of the special geometry of the complex, calculations with polarized light parallel and perpendicular to the molecular plane were carried out, yielding a reasonable reproduction of the experimental data from another report for cluster sizes equal to or higher than 45 atoms. The highly symmetric square planar structure of the complex was found to be responsible for the unusual amplitude of the multiple-scattering (MS) contributions to the EXAFS spectra. Spectra in this region were fitted using the FEFFIT EXAFS analysis program, taking into account only the MS paths that simultaneously have both a high amplitude, as calculated with the ab initio code FEFF, and a small Debye-Waller factor, as estimated by the independent-vibration approximation model. Fitting results yielded very similar structures for the Ni2+ complex in the solid state and in solution, though the larger Debye-Waller factors found for the solid suggest higher static disorder in this state.