We briefly review our implementation of the real-space Green's 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.

No abstract

We present ab initio calculations of frequency-dependent linear and nonlinear optical responses based on real-time time-dependent density functional theory for arbitrary photonic molecules. This approach is based on an extension of an approach previously implemented for a linear response using the electronic structure program SIESTA. Instead of calculating excited quantum states, which can be a bottleneck in frequency-space calculations, the response of large molecular systems to time-varying electric fields is calculated in real time. This method is based on the finite field approach generalized to the dynamic case. To speed the nonlinear calculations, our approach uses Gaussian enveloped quasimonochromatic external fields. We thereby obtain the frequency-dependent second harmonic generation beta(-2omega;omega,omega), the dc nonlinear rectification beta(0;-omega,omega), and the electro-optic effect beta(-omega;omega,0). The method is applied to nanoscale photonic nonlinear optical molecules, including p-nitroaniline and the FTC chromophore, i.e., 2-[3-Cyano-4-(2-{5-[2-(4-diethylamino-phenyl)-vinyl]-thiophen-2-yl}-vinyl)-5,5-dimethyl-5H-furan-2-ylidene]-malononitrile, and yields results in good agreement with experiment.

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.

An approach is presented for theoretical calculations of the Debye-Waller factors in x-ray absorption spectra. These factors are represented in terms of the cumulant expansion up to third order. They account respectively for the net thermal expansion σ(1) (T ), the mean-square relative displacements σ 2 (T ), and the asymmetry of the pair distribution function σ (3) (T ). Similarly, we obtain Debye-Waller factors for x-ray and neutron scattering in terms of the mean-square vibrational amplitudes u 2 (T ). Our method is based on density functional theory calculations of the dynamical matrix, together with an efficient Lanczos algorithm for projected phonon spectra within the quasiharmonic approximation. Due to anharmonicity in the interatomic forces, the results are highly sensitive to variations in the equilibrium lattice constants, and hence to the choice of exchangecorrelation potential. In order to treat this sensitivity, we introduce two prescriptions: one based on the local density approximation, and a second based on a modified generalized gradient approximation. Illustrative results for the leading cumulants are presented for several materials and compared with experiment and with correlated Einstein and Debye models. We also obtain Born-von Karman parameters and corrections due to perpendicular vibrations.

The nature of local atomic and electronic structure at the nanoscale is of both fundamental and technological importance. For example, supported metal nanoclusters exhibit a number of unusual phenomena including large structural disorder and bond-length contraction with increasing temperature. We investigate this behavior for a prototypical ten atom Pt cluster supported on ␥ alumina using temperature-dependent, real-time simulations based on density functional theory/molecular-dynamics and x-ray spectroscopy theory. The simulations reveal a complex dynamical structure on multiple-time scales including librational motion of the center of mass and fluctuating bonding characteristics, which explain many of the unusual properties.

We report autodetachment spectra of the mass-selected, anionic water clusters, (H2O)n−, n=2, 3, 5–9, 11 in the OH stretching region (3000–4000 cm−1), and interpret the spectra with the aid of ab initio calculations. For n⩾5, the spectra are structured and are generally dominated by an intense doublet, split by about 100 cm−1, which gradually shifts toward lower energy with increasing cluster size. This behavior indicates that the n=5–11 clusters share a common structural motif. The strong bands appear in the frequency region usually associated with single-donor vibrations of water molecules embedded in extended networks, and theoretical calculations indicate that the observed spectra are consistent with linear “chainlike” (H2O)n− species. We test this assignment by recording the spectral pattern of the cooled (argon solvated) HDO⋅(D2O)5− isotopomer over the entire OH stretching frequency range.

X-ray absorption spectroscopy (XAS) is a widely used materials characterization technique to determine oxidation states, coordination environment, and other local atomic structure information. Analysis of XAS relies on comparison of measured spectra to reliable reference spectra. However, existing databases of XAS spectra are highly limited both in terms of the number of reference spectra available as well as the breadth of chemistry coverage. In this work, we report the development of XASdb, a large database of computed reference XAS, and an Ensemble-Learned Spectra IdEntification (ELSIE) algorithm for the matching of spectra. XASdb currently hosts more than 800,000 K-edge X-ray absorption near-edge spectra (XANES) for over 40,000 materials from the open-science Materials Project database. We discuss a high-throughput automation framework for FEFF calculations, built on robust, rigorously benchmarked parameters. FEFF is a computer program uses a real-space Green's function approach to calculate X-ray absorption spectra. We will demonstrate that the ELSIE algorithm, which combines 33 weak "learners" comprising a set of preprocessing steps and a similarity metric, can achieve up to 84.2% accuracy in identifying the correct oxidation state and coordination environment of a test set of 19 K-edge XANES spectra encompassing a diverse range of chemistries and crystal structures. The XASdb with the ELSIE algorithm has been integrated into a web application in the Materials Project, providing an important new public resource for the analysis of XAS to all materials researchers. Finally, the ELSIE algorithm itself has been made available as part of veidt, an open source machine-learning library for materials science.npj Computational Materials (2018) 4:12 ; doi:10.1038/s41524-018-0067-x INTRODUCTION X-ray absorption spectroscopy (XAS) is a widely used technique in the study of the properties, physical states, and local environments of materials.1-3 When incident X-ray photons with energy greater than the binding energy are absorbed by an atom, a corelevel electron is removed from its quantum level. In XAS, the absorption coefficient, μ(E) is measured as a function of X-ray energy E. Detailed descriptions of X-ray absorption theory and equation have been included in many excellent books and review papers.