Source author record

J. J. Rehr

J. J. Rehr appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Ab initio multiplet plus cumulant approach for correlation effects in x-ray photoelectron spectroscopy

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 (LFMT) 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 Greens 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 problematic. 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.

preprint2016arXiv

Efficient implementation of core-excitation Bethe Salpeter equation calculations

We present an efficient implementation of the Bethe-Salpeter equation (BSE) method for obtaining core-level spectra including x-ray absorption (XAS), x-ray emission (XES), and both resonant and non-resonant inelastic x-ray scattering spectra (N/RIXS). Calculations are based on density functional theory (DFT) electronic structures generated either by abinit or Quantumespresso, both plane-wave basis, pseudopotential codes. This electronic structure is improved through the inclusion of a GW self energy. The projector augmented wave technique is used to evaluate transition matrix elements between core-level and band states. Final two-particle scattering states are obtained with the NIST core-level BSE solver (NBSE). We have previously reported this implementation, which we refer to as ocean (Obtaining Core Excitations from Ab initio electronic structure and NBSE) [Phys. Rev. B 83, 115106 (2011)]. Here, we present additional efficiencies that enable us to evaluate spectra for systems ten times larger than previously possible; containing up to a few thousand electrons. These improvements include the implementation of optimal basis functions that reduce the cost of the initial DFT calculations, more complete parallelization of the screening calculation and of the action of the BSE Hamiltonian, and various memory reductions. Scaling is demonstrated on supercells of SrTiO_3 and example spectra for the organic light emitting molecule Tris-(8-hydroxyquinoline)aluminum (Alq_3 ) are presented. The ability to perform large-scale spectral calculations is particularly advantageous for investigating dilute or non-periodic systems such as doped materials, amorphous systems, or complex nano-structures.

preprint2016arXiv

Particle-hole cumulant approach for inelastic losses in x-ray spectra

Inelastic losses in core level x-ray spectra arise from many-body excitations, leading to broadening and damping as well as satellite peaks in x-ray photoemission (XPS) and x-ray absorption (XAS) spectra. Here we present a practical approach for calculating these losses based on a cumulant representation of the particle-hole Green's function, a quasi-boson approximation, and a partition of the cumulant into extrinsic, intrinsic and interference terms. The intrinsic losses are calculated using real-time, time-dependent density functional theory while the extrinsic losses are obtained from the GW approximation of the photo-electron self-energy and the interference terms are approximated. These effects are included in the spectra using a convolution with an energy dependent particle-hole spectral function. The approach elucidates the nature of the spectral functions in XPS and XAS and explains the significant cancellation between extrinsic and intrinsic losses. Edge-singularity effects in metals are also accounted for. Illustrative results are presented for the XPS and XAS for both weakly and more correlated systems.

preprint2014arXiv

Cumulant expansion for phonon contributions to the electron spectral function

We describe an approach for calculations of phonon contributions to the electron spectral function, including both quasiparticle properties and satellites. The method is based on a cumulant expansion for the retarded one-electron Green's function and a many-pole model for the electron self-energy. The electron-phonon couplings are calculated from the Eliashberg functions, and the phonon density of states is obtained from a Lanczos representation of the phonon Green's function. Our calculations incorporate ab initio dynamical matrices and electron-phonon couplings from the density functional theory code ABINIT. Illustrative results are presented for several elemental metals and for Einstein and Debye models with a range of coupling constants. These are compared with experiment and other theoretical models. Estimates of corrections to Migdal's theorem are obtained by comparing with leading order contributions to the self-energy, and are found to be significant only for large electron-phonon couplings at low temperatures.

preprint2014arXiv

Generalized particle/hole cumulant approximation for the electron Green's function

The cumulant expansion is a powerful approach for including correlation effects in electronic structure calculations beyond the GW approximation. However, current implementations are incomplete since they ignore terms that lead to partial occupation numbers and satellites both above and below the Fermi energy. These limitations are corrected here with a generalized cumulant approximation that includes both particle and hole contributions within a retarded Green's function formalism. The computational effort is still comparable to GW, and the method can be extended easily to finite temperature. The approach is illustrated with calculations for the homogeneous electron gas and comparisons to experiment and other methods.

preprint2014arXiv

Real time cumulant approach for charge transfer satellites in x-ray photoemission spectra

X-ray photoemission spectra generally exhibit satellite features in addition to the quasi-particle peaks due to many-body excitations, which have been of considerable theoretical and experimental interest. However, the satellites attributed to charge-transfer (CT) excitations in correlated materials have proved difficult to calculate from first principles. Here we report a real-time, real-space approach for such calculations based on a cumulant representation of the core-hole Green's function and time-dependent density functional theory. This approach also yields an interpretation of CT satellites in terms of a complex oscillatory, transient response to a suddenly created core hole. Illustrative results for TiO$_2$ and NiO are in good agreement with experiment.

preprint2013arXiv

Accuracy of generalized gradient approximation functionals for density functional perturbation theory calculations

We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We investigate a wide variety of materials including a semiconductor (silicon), a metal (copper), and various insulators (SiO$_2$ $α$-quartz and stishovite, ZrSiO$_4$ zircon, and MgO periclase). For the structural properties, we find that PBEsol and WC are the closest to the experiments and AM05 performs only slightly worse. All three functionals actually improve over LDA and PBE in contrast with HTBS, which is shown to fail dramatically for $α$-quartz. For the vibrational and thermodynamical properties, LDA performs surprisingly very good. In the majority of the test cases, it outperforms PBE significantly and also the WC, PBEsol and AM05 functionals though by a smaller margin (and to the detriment of structural parameters). On the other hand, HTBS performs also poorly for vibrational quantities. For the dielectric properties, none of the functionals can be put forward. They all (i) fail to reproduce the electronic dielectric constant due to the well-known band gap problem and (ii) tend to overestimate the oscillator strengths (and hence the static dielectric constant).

preprint2013arXiv

Charge transfer satellites in x-ray spectra of transition metal oxides

Strongly correlated materials such as transition metal oxides (TMOs) often exhibit large satellites in their x-ray photoemission (XPS) and x-ray absorption spectra (XAS). These satellites arise from localized charge-transfer (CT) excitations that accompany the sudden creation of a core hole. Here we use a two-step approach to treat such excitations in a localized system embedded in a condensed system and coupled to a photoelectron. The total XAS is then given by a convolution of a spectral function representing the localized excitations and the XAS of the extended system. The local system is modeled roughly in terms of a simple three-level model, leading to a double-pole approximation for the spectral function that represents dynamically weighted contributions from the dominant neutral and charge-transfer excitations. This method is implemented using a resolvent approach, with potentials, radial wave-functions and matrix elements from the real-space Green's function code FEFF, and parameters fitted to XPS experiments. Representative calculations for several TMOs are found to be in reasonable agreement with experiment.

preprint2013arXiv

Shake-Rattle-and-Roll: A Model of Dynamic Structural Disorder in Supported Nanoscale Catalysts

We investigate the effects of "dynamic structural disorder" (DSD) on the behavior of 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 energy, charge, and other extensive quantities as well as large surface effects. The DSD is driven 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 DSD is based on a combination of statistical mechanics, transient coupled-oscillator models, and real-time DFT/MD simulations. This approach treats thermal and dynamic effects over multiple time-scales, including bond-stretching and -bending vibrations, DSD, and transient tethering to the substrate at longer ps time-scales. 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 affects both the prefactors and distribution of energy barriers in reaction rates, and thus can strongly affect catalytic activity at the nano-scale.

preprint2011arXiv

GW+U real-space Green's function calculations of x-ray spectra

The Hubbard model is implemented in real-space Green's function calculations of x-ray spectra using an effective self-energy adapted from the LSDA+U method of Anisimov et al. This self-energy consists of an energy-dependent many-pole approximation to the GW self-energy with an additive correction due to on-site Coulomb repulsion among the partially filled localized-electron states. This leads to a GW+U approach which provides an efficient procedure to account for local correlation effects on x-ray spectra. Results are presented for the spin and angular momentum projected density of states of MnO, NiO, and La_(2-x)Sr_xCuO_4 (LSCO), for the K-edge x-ray spectra of O atoms in MnO and NiO, and the unoccupied electronic states and O K-edge spectra of undoped LSCO. The method is found to yield reasonable agreement with experiment.

preprint2011arXiv

Real Space Green's Function Approach to RIXS

We present an ab initio theory of core- and valence resonant inelastic x-ray scattering (RIXS) based on a real-space multiple scattering Green's function formalism and a quasi-boson model Hamiltonian. Simplifying assumptions are made which lead to an approximation of the RIXS spectrum in terms of a convolution of an effective x-ray absorption signal with the x-ray emission signal. Additional many body corrections are incorporated in terms of an effective energy dependent spectral function. Example calculations of RIXS are found to give qualitative agreement with experimental data. Our approach also yields simulations of lifetime-broadening suppressed XAS, as observed in high energy resolutionfluorescence detection experiment (HERFD). Finally possible improvements to our approach are briefly discussed.

preprint2011arXiv

Theoretical optical and x-ray spectra of liquid and solid H_2O

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, quasi-particle self energy corrections, and a BSE treatment of particle-hole interactions. The 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.

preprint2011arXiv

X-ray absorption Debye-Waller factors from ab initio molecular dynamics

An ab initio equation of motion method is introduced to calculate the temperature-dependent mean square vibrational amplitudes which appear in the Debye-Waller factors in x-ray absorption, x-ray scattering, and related spectra. The approach avoids explicit calculations of phonon-modes, and is based instead on calculations of the displacement-displacement time correlation function from ab initio density functional theory molecular dynamics simulations. The method also yields the vibrational density of states and thermal quantities such as the lattice free energy. Illustrations of the method are presented for a number of systems and compared with other methods and experiment.

preprint2011arXiv

X-ray absorption near-edge spectra of overdoped La_2-xSr_xCuO_4 high-T_c superconductors

We present results for realistic modeling of the x-ray absorption near edge structure (XANES) of the overdoped high-T_c superconductor La_2-xSr_xCuO_4 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.

preprint2010arXiv

Ab initio analysis of the x-ray absorption spectrum of the myoglobin-carbon monoxide complex: Structure and vibrations

We present a comparison between Fe K-edge x-ray absorption spectra of carbonmonoxy-myoglobin and its simulation based on density-functional theory determination of the structure and vibrations and spectral simulation with multiple-scattering theory. An excellent comparison is obtained for the main part of the molecular structure without any structural fitting parameters. The geometry of the CO ligand is reliably determined using a synergic approach to data analysis. The methodology underlying this approach is expected to be especially useful in similar situations in which high-resolution data for structure and vibrations are available.

preprint2010arXiv

Basis set effects on the hyperpolarizability of CHCl_3: Gaussian-type orbitals, numerical basis sets and real-space grids

Calculations of the hyperpolarizability are typically much more difficult to converge with basis set size than the linear polarizability. In order to understand these convergence issues and hence obtain accurate ab initio values, we compare calculations of the static hyperpolarizability of the gas-phase chloroform molecule (CHCl_3) using three different kinds of basis sets: Gaussian-type orbitals, numerical basis sets, and real-space grids. Although all of these methods can yield similar results, surprisingly large, diffuse basis sets are needed to achieve convergence to comparable values. These results are interpreted in terms of local polarizability and hyperpolarizability densities. We find that the hyperpolarizability is very sensitive to the molecular structure, and we also assess the significance of vibrational contributions and frequency dispersion.

preprint2010arXiv

Bethe-Salpeter Equation Calculations of Core Excitation Spectra

We present a hybrid approach for GW/Bethe-Salpeter Equation (BSE) calculations of core excitation spectra, including x-ray absorption (XAS), electron energy loss spectra (EELS), and non-resonant inelastic x-ray scattering (NRIXS). The method is based on {\it ab initio} wavefunctions 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 GW self-energy model to account for final-state broadening and self-energy shifts. Multiplet effects are also 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. We compare the results to experiments and results obtained using other theoretical approaches.

preprint2010arXiv

Phase Space Wannier Functions in Electronic Structure Calculations

We consider the applicability of phase space Wannier functions" to electronic structure calculations. These generalized Wannier functions are analogous to localized plane waves and constitute a complete, orthonormal set which is exponentially localized both in position and momentum space. Their properties are described and an illustrative application to bound states in one dimension is presented. Criteria for choosing basis set size and lattice constant are discussed based on semi-classical, phase space considerations. Convergence of the ground state energy with respect to basis size is evaluated. Comparison with plane-waves basis sets indicates that the number of phase space Wannier functions needed for convergence can be signicantly smaller in three dimensions. PACS: 71.10.+x, 71.50.+t

preprint2008arXiv

Density matrix calculation of optical constants from optical to x-ray frequencies

We present a theory of linear optical constants based on a single-particle density matrix and implemented in an extension of the real-space multiple scattering code FEFF. This approach avoids the need to compute wave-functions explicitly, and yields efficient calculations for frequencies ranging from the IR to hard x-rays, and applicable to arbitrary aperiodic systems. Our approach is illustrated with calculations of optical properties and applications for several materials.

preprint2006arXiv

Background Proportional Enhancement of the Extended Fine Structure in Nonresonant Inelastic X-ray Scattering

We report new measurements and calculations of the non-resonant inelastic x-ray scattering (NRIXS) from Mg and Al for a wide range of momentum transfers, q. Extended oscillations in the dynamic structure factor S(q,w) due to scattering from the 2p and 2s orbitals (i.e. L-edges) are observed out to more than 150 eV past the binding energy. These results are discussed in context of the recently proposed representation of S(q,w) for core shells as an atomic background modulated by interference between different photoelectron scattering paths, in analogy to the standard treatment of extended x-ray absorption fine structure. In agreement with this representation, we find a strong increase in the atomic background with increasing q with a concomitant enhancement in the amplitude of the extended fine structure. This effect should be generic and hence may enable improved measurement of the extended fine structure in a wide range of materials containing low-Z elements.

preprint2004arXiv

Local environment of Nitrogen in GaN{y}As{1-y} epilayers on GaAs (001) studied using X-ray absorption near edge spectroscopy

X-ray absorption near-edge spectroscopy (XANES) is used to study the N environment in bulk GaN and in GaN{y}As{1-y} epilayers on GaAs (001), for y \~5%. Density-functional optimized structures were used to predict XANES via multiple-scattering theory. We obtain striking agreement for pure GaN. An alloy model with nitrogen pairs on Ga accurately predicts the threshold energy, the width of the XANES ``white line'', and features above threshold, for the given X-ray polarization. The presence of N-pairs may point to a role for molecular N_2 in epitaxial growth kinetics.