Micah P. Prange

Parameter-free calculations of X-ray spectra with FEFF9

We briefly review our implementation of the real-space Greens function (RSGF) approach for calculations of X-ray spectra, focusing on recently developed parameter free models for dominant many-body effects. Although the RSGF approach has been widely used both for near edge (XANES) and extended (EXAFS) ranges, previous implementations relied on semi-phenomenological methods, e.g., the plasmon-pole model for the self-energy, the final-state rule for screened core hole effects, and the correlated Debye model for vibrational damping. Here we describe how these approximations can be replaced by efficient ab initio models including a many-pole model of the self-energy, inelastic losses and multiple-electron excitations; a linear response approach for the core hole; and a Lanczos approach for Debye-Waller effects. We also discuss the implementation of these models and software improvements within the FEFF9 code, together with a number of examples.

Scientific Computing in the Cloud

Large, virtualized pools of computational resources raise the possibility of a new, advantageous computing paradigm for scientific research. To help achieve this, new tools make the cloud platform behave virtually like a local homogeneous computer cluster, giving users access to high-performance clusters without requiring them to purchase or maintain sophisticated hardware.

Atomic-Resolution Imaging of Spin-State Superlattices in Nanopockets within Cobaltite Thin Films

Certain cobalt oxides are known to exhibit ordered Co spin states, as determined from macroscopic techniques. Here we report real-space atomic-resolution imaging of Co spin-state ordering in nanopockets of La(0.5)Sr(0.5)CoO(3-δ) thin films. Unlike the bulk material, where no Co spin-state ordering is found, thin films present a strain-induced domain structure due to oxygen vacancy ordering, inside of which some nanometer sized domains show high-spin Co ions in the planes containing O vacancies and low-spin Co ions in the stoichiometric planes. First-principles calculations provide support for this interpretation.

Many-pole model of inelastic losses in x-ray absorption spectra

Inelastic losses are crucial to a quantitative analysis of x-ray absorption spectra. However, current treatments are semi-phenomenological in nature. Here a first-principles, many-pole generalization of the plasmon-pole model is developed for improved calculations of inelastic losses. The method is based on the GW approximation for the self-energy and real space multiple scattering calculations of the dielectric function for a given system. The model retains the efficiency of the plasmonpole model and is applicable both to periodic and aperiodic materials over a wide energy range. The same many-pole model is applied to extended GW calculations of the quasiparticle spectral function. This yields estimates of multi-electron excitation effects, e.g., the many-body amplitude factor S 2 0 due to intrinsic losses. Illustrative calculations are compared with other GW calculations of the self-energy, the inelastic mean free path, and experimental x-ray absorption spectra.

Ab initio calculations of electron inelastic mean free paths and stopping powers

A method is presented for first-principles calculations of inelastic mean free paths and stopping powers in condensed matter over a broad energy range. The method is based on {\it ab initio} calculations of the dielectric function in the long wavelength limit using a real-space Greens function formalism, together with extensions to finite momentum transfer. From these results we obtain the loss function and related quantities such as optical-oscillator strengths and mean excitation energies. From a many-pole representation of the dielectric function we then obtain the electron self-energy and inelastic mean free paths (IMFP). Finally using our calculated dielectric function and the optical-data model of Fernandez-Varea {\it et al}., we obtain collision stopping powers (CSP) and penetration ranges. The results are consistent with semi-empirical approaches and with experiment.

Physically motivated global alignment method for electron tomography

Electron tomography is widely used for nanoscale determination of 3-D structures in many areas of science. Determining the 3-D structure of a sample from electron tomography involves three major steps: acquisition of sequence of 2-D projection images of the sample with the electron microscope, alignment of the images to a common coordinate system, and 3-D reconstruction and segmentation of the sample from the aligned image data. The resolution of the 3-D reconstruction is directly influenced by the accuracy of the alignment, and therefore, it is crucial to have a robust and dependable alignment method. In this paper, we develop a new alignment method which avoids the use of markers and instead traces the computed paths of many identifiable ‘local’ center-of-mass points as the sample is rotated. Compared with traditional correlation schemes, the alignment method presented here is resistant to cumulative error observed from correlation techniques, has very rigorous mathematical justification, and is very robust since many points and paths are used, all of which inevitably improves the quality of the reconstruction and confidence in the scientific results.

Reducing radiation damage in macromolecular crystals at synchrotron sources

A new strategy is presented to reduce primary X-ray damage in macromolecular crystallography. The strategy is based on separating the diffracting and damaged regions as much as feasible. The source of the radiation damage to macromolecular crystals is from two primary mechanisms: the direct excitations of electrons by absorption, and inelastic scattering of the X-rays. The first produces photoelectrons with their accompanying Auger electrons from relaxation of the core hole and the second creates Compton electrons. The properties of these two mechanisms and calculations of primary X-ray damage quantify how to modify the spatial distribution of X-rays to reduce the deleterious effects of radiation damage. By focusing the incident X-rays into vertical stripes, it is estimated that the survival (the time during which quality diffraction data can be obtained with a given X-ray flux) of large crystals can be increased by at least a factor of 1.6, while for very small platelet crystals the survival can be increased by up to a factor of 14.

Monte Carlo simulation of gamma-ray response of BaF2 and CaF2

We have employed a Monte Carlo (MC) method to study intrinsic properties of two alkaline-earth halides, namely, BaF2 and CaF2, relevant to their use as radiation detector materials. The MC method follows the fate of individual electron-hole (e-h) pairs and thus allows for a detailed description of the microscopic structure of ionization tracks created by incident γ-ray radiation. The properties of interest include the mean energy required to create an e-h pair, W, Fano factor, F, the maximum theoretical light yield, and the spatial distribution of e-h pairs resulting from γ-ray excitation. Although W and F vary with incident photon energy at low energies, they tend to constant values at energies higher than 1 keV. W is determined to be 18.9 and 19.8 eV for BaF2 and CaF2, respectively, in agreement with published data. The e-h pair spatial distributions exhibit a linear distribution along the fast electron tracks with high e-h pair densities at the end of the tracks. Most e-h pairs are created by interband...

Real space calculation of optical constants from optical to x-ray frequencies

We present a theory of linear optical constants based on the single-particle density operator and implemented in an extension of the real space multiple scattering code known as 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, which is applicable to arbitrary aperiodic systems. The approach is illustrated with calculations of optical properties and applications for several materials and compared with existing tabulations.

Radiation response of inorganic scintillators: Insights from Monte Carlo simulations

The spatial and temporal scales of hot particle thermalization in inorganic scintillators are critical factors determining the extent of second- and third-order nonlinear quenching in regions with high densities of electron-hole pairs, which, in turn, leads to the light yield nonproportionality observed, to some degree, for all inorganic scintillators. Therefore, kinetic Monte Carlo simulations were performed to calculate the distances traveled by hot electrons and holes as well as the time required for the particles to reach thermal energy following γ-ray irradiation. CsI, a common scintillator from the alkali halide class of materials, was used as a model system. Two models of quasi-particle dispersion were evaluated, namely, the effective mass approximation model and a model that relied on the group velocities of electrons and holes determined from band structure calculations. Both models predicted rapid electron-hole pair recombination over short distances (a few nanometers) as well as a significant extent of charge separation between electrons and holes that did not recombine and reached thermal energy. However, the effective mass approximation model predicted much longer electron thermalization distances and times than the group velocity model. Comparison with limited experimental data suggested that the group velocity model provided more accurate predictions. Nonetheless, both models indicated that hole thermalization is faster than electron thermalization and thus is likely to be an important factor determining the extent of third-order nonlinear quenching in high-density regions. The merits of different models of quasi-particle dispersion are also discussed.