Search

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.

The transformative impact of modern computational paradigms and technologies, such as high-performance computing (HPC), quantum computing, and cloud computing, has opened up profound new opportunities for scientific simulations. Scalable computational chemistry is one beneficiary of this technological progress. The main focus of this paper is on the performance of various quantum chemical formulations, ranging from low-order methods to high-accuracy approaches, implemented in different computational chemistry packages and libraries, such as NWChem, NWChemEx, Scalable Predictive Methods for Excitations and Correlated Phenomena, ExaChem, and Fermi-Löwdin orbital self-interaction correction on Azure Quantum Elements, Microsoft's cloud services platform for scientific discovery. We pay particular attention to the intricate workflows for performing complex chemistry simulations, associated data curation, and mechanisms for accuracy assessment, which is demonstrated with the Arrows automated workflow for high throughput simulations. Finally, we provide a perspective on the role of cloud computing in supporting the mission of leadership computational facilities., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Supported noble metal catalysts, ubiquitous in chemical technology, often undergo dynamic transformations between reduced and oxidized states-which influence the metal nuclearities, oxidation states, and catalytic properties. In this investigation, we report the results of in situ X-ray absorption spectroscopy, scanning transmission electron microscopy, and other physical characterization techniques, bolstered by density functional theory, to elucidate the structural transformations of a set of MgO-supported palladium catalysts under oxidative treatment conditions. As the calcination temperature increased, the as-synthesized supported metallic palladium nanoparticles underwent oxidation to form palladium oxides (at approximately 400 °C), which, at approximately 500 °C, were oxidatively fragmented to form mixtures of atomically dispersed palladium cations. The data indicate two distinct types of atomically dispersed species: palladium cations located at MgO steps and those embedded in the first subsurface layer of MgO. The former exhibit significantly higher (>500 times) catalytic activity for ethylene hydrogenation than the latter. The results pave the way for designing highly active and stable supported palladium hydrogenation catalysts with optimized metal utilization., Competing Interests: The authors declare no conflict of interests., (This journal is © The Royal Society of Chemistry.)

Photoelectron spectroscopy (PES) is a standard experimental method for material characterization, but its interpretation can be hampered by its reliance on standard materials. To facilitate the study of unknown systems, theoretical methods are desirable. Here, we present a real-time equation-of-motion coupled cluster (RT-EOM-CC) approach for valence PES, extending our core-level development. We demonstrate that RT-EOM-CC yields ionization energies and spectral functions in good agreement with experimental and CI-based results, even for some more correlated cases.

Soft-oxidant-assisted methane coupling has emerged as a promising pathway to upgrade methane from natural gas sources to high-value commodity chemicals, such as ethylene, at selectivities higher than those associated with oxidative (O 2 ) methane coupling (OCM). To date, few studies have reported investigations into the electronic structure and the microscopic physical structure of catalytic active sites present in the binary metal oxide catalyst systems that are known to be effective for this reaction. Correlating the catalyst activity to specific active site structures and electronic properties is an essential aspect of catalyst design. Here, we used X-ray absorption spectroscopy at the Ca K-edge to ascertain the most probable local environment of Ca in the ZnO-supported Ca oxide catalysts. These catalysts are shown here to be active for N 2 O-assisted methane coupling (N 2 O-OCM) and have previously been reported to be active for CO 2 -assisted methane coupling (CO 2 -OCM). X-ray absorption near edge structure features at multiple Ca loadings are interpreted through simulated spectra derived from ab initio full multiple scattering calculations. These simulations included consideration of CaO structures organized in multiple spatial arrangements-linear, planar, and cubic-with separate analyses of Ca atoms in the surfaces and bulk of the three-dimensional structures. The morphology of the oxide clusters was found to influence the various regions of the X-ray absorption spectrum differently. Experiment and theory show that for low-Ca-loading catalysts (≤1 mol %), which contain sites particularly active for methane coupling, Ca primarily exists in an oxidized state that is consistent with the coordination environment of Ca ions in one- and two-dimensional clusters. In addition to their unique nanoscale structures, the spectra also indicate that these clusters have varying degrees of undercoordinated surface Ca atoms that could further influence their catalytic activities. The local Ca structure was correlated to methane coupling activity from N 2 O-OCM and previously reported CO 2 -OCM reactor studies. This study provides a unique perspective on the relationship between the catalyst physical and electronic structure and active sites for soft-oxidant-assisted methane coupling, which can be used to inform future catalyst development., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.

Electrolyte cations can have significant effects on the kinetics and selectivity of electrocatalytic reactions. We show an atypical mechanism through which electrolyte cations can impact electrocatalyst performance─direct incorporation of the cation into the oxide electrocatalyst lattice. We investigate the transformations of copper electrodes in alkaline electrochemistry through operando X-ray absorption spectroscopy in KOH and Ba(OH) 2 electrolytes. In KOH electrolytes, both the near-edge structure and extended fine-structure agree with previous studies; however, the X-ray absorption spectra vary greatly in Ba(OH) 2 electrolytes. Through a combination of electronic structure modeling, near-edge simulation, and postreaction characterization, we propose that Ba 2+ cations are directly incorporated into the lattice and form an ordered BaCuO 2 phase at potentials more oxidizing than 200 mV vs the normal hydrogen electrode (NHE). BaCuO 2 formation is followed by further oxidation to a bulk Cu 3+ -like Ba x Cu y O z phase at 900 mV vs NHE. Additionally, during reduction in Ba(OH) 2 electrolyte, we find both Cu-O bonds and Cu-Ba scattering persist at potentials as low as -400 mV vs NHE. To our knowledge, this is the first evidence for direct oxidative incorporation of an electrolyte cation into the bulk lattice to form a mixed oxide electrode. The oxidative incorporation of electrolyte cations to form mixed oxides could open a new route for the in situ formation of active and selective oxidation electrocatalysts.

Sustainable production of rare earth elements (REEs) is critical for technologies needed for climate change mitigation, including wind turbines and electric vehicles. However, separation technologies currently used in REE production have large environmental footprints, necessitating more sustainable strategies. Aqueous, affinity-based separations are examples of such strategies. To make these technologies feasible, it is imperative to connect aqueous ligand structure to ligand selectivity for individual REEs. As a step toward this goal, we analyzed the extended X-ray absorption fine structure (EXAFS) of four lanthanides (La, Ce, Pr, and Nd) complexed by a common REE chelator, ethylenediaminetetraacetic acid (EDTA) to determine the aqueous-phase structure. Reference structures from density functional theory (DFT) were used to help fit the EXAFS spectra. We found that all four Ln-EDTA coordination complexes formed 9-coordinate structures with 6 coordinating atoms from EDTA (4 carboxyl oxygen atoms and 2 nitrogen atoms) and 3 oxygen atoms from water molecules. All EXAFS fits were of high quality ( R -factor < 0.02) and showed decreasing average first-shell coordination distance across the series (2.62-2.57 Å from La-Nd), in agreement with DFT (2.65-2.56 Å from La-Nd). The insights determined herein will be useful in the development of ligands for sustainable rare earth elements (REE) separation technologies.

We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [ J. Chem. Phys. 2020, 152, 174113] within the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a massively parallel heterogeneous tensor library designed for utilizing forthcoming exascale computing resources. The two-body electron repulsion matrix elements are Cholesky-decomposed, and we imposed spin-explicit forms of the various operators when evaluating the tensor contractions. Unlike our previous real algebra Tensor Contraction Engine (TCE) implementation, the TAMM implementation supports fully complex algebra. The RT-EOM-CC singles (S) and doubles (D) time-dependent amplitudes are propagated using a first-order Adams-Moulton method. This new implementation shows excellent scalability tested up to 500 GPUs using the Zn-porphyrin molecule with 655 basis functions, with parallel efficiencies above 90% up to 400 GPUs. The TAMM RT-EOM-CCSD was used to study core photoemission spectra in the formaldehyde and ethyl trifluoroacetate (ESCA) molecules. Simulations of the latter involve as many as 71 occupied and 649 virtual orbitals. The relative quasiparticle ionization energies and overall spectral functions agree well with available experimental results.

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N - 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

X-Ray and related spectroscopies are powerful probes of atomic, vibrational, and electronic structure. In order to unlock the full potential of such experimental techniques, accurate and efficient theoretical and computational approaches are essential. Here we review the status of a variety of first-principles and nearly first principles techniques for X-ray spectroscopies such as X-ray absorption, X-ray emission, and X-ray photoemission, with a focus on Green's function based methods. In particular, we describe the current state of multiple scattering Green's function techniques available in the FEFF10 code and cumulant Green's function techniques for including the effects of many-body electronic excitations. Illustrative examples are shown for a variety of materials and compared with other theoretical and experimental results.

Atomically dispersed metals on metal oxide supports are a rapidly growing class of catalysts. Developing an understanding of where and how the metals are bonded to the supports is challenging because support surfaces are heterogeneous, and most reports lack a detailed consideration of these points. Herein, we report two atomically dispersed CO oxidation catalysts having markedly different metal-support interactions: platinum in the first layer of crystalline MgO powder and platinum in the second layer of this support. Structural models have been determined on the basis of data and computations, including those determined by extended X-ray absorption fine structure and X-ray absorption near edge structure spectroscopies, infrared spectroscopy of adsorbed CO, and scanning transmission electron microscopy. The data demonstrate the transformation of surface to subsurface platinum as the temperature of sample calcination increased. Catalyst performance data demonstrate the lower activity but greater stability of the subsurface platinum than of the surface platinum.

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.

Atomically dispersed supported metal catalysts offer new properties and the benefits of maximized metal accessibility and utilization. The characterization of these materials, however, remains challenging. Using atomically dispersed platinum supported on crystalline MgO (chosen for its well-defined bonding sites) as a prototypical example, we demonstrate how systematic density functional theory calculations for assessing all the potentially stable platinum sites, combined with automated analysis of extended X-ray absorption fine structure (EXAFS) spectra, leads to unbiased identification of isolated, surface-enveloped platinum cations as the catalytic species for CO oxidation. The catalyst has been characterized by atomic-resolution imaging and EXAFS and high-energy resolution fluorescence detection X-ray absorption near edge spectroscopy. The proposed platinum sites are in agreement with experiment. This theory-guided workflow leads to rigorously determined structural models and provides a more detailed picture of the structure of the catalytically active site than what is currently possible with conventional EXAFS analyses. As this approach is efficient and agnostic to the metal, support, and catalytic reaction, we posit that it will be of broad interest to the materials characterization and catalysis communities.

The real-space Green's function code FEFF has been extensively developed and used for calculations of X-ray and related spectra, including X-ray absorption (XAS), X-ray emission (XES), inelastic X-ray scattering, and electron energy-loss spectra. The code is particularly useful for the analysis and interpretation of the XAS fine-structure (EXAFS) and the near-edge structure (XANES) in materials throughout the periodic table. Nevertheless, many applications, such as non-equilibrium systems, and the analysis of ultra-fast pump-probe experiments, require extensions of the code including finite-temperature and auxiliary calculations of structure and vibrational properties. To enable these extensions, we have developed in tandem a new version FEFF10 and new FEFF-based workflows for the Corvus workflow manager, which allow users to easily augment the capabilities of FEFF10 via auxiliary codes. This coupling facilitates simplified input and automated calculations of spectra based on advanced theoretical techniques. The approach is illustrated with examples of high-temperature behavior, vibrational properties, many-body excitations in XAS, super-heavy materials, and fits of calculated spectra to experiment.

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.)

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.

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.

Rh-based catalysts have shown promise for the direct conversion of syngas to higher oxygenates. Although improvements in higher oxygenate yield have been achieved by combining Rh with metal oxide promoters, details of the structure of the promoted catalyst and the role of the promoter in enhancing catalytic performance are not well understood. In this work, we show that MoO 3 -promoted Rh nanoparticles form a novel catalyst structure in which Mo substitutes into the Rh surface, leading to both a 66-fold increase in turnover frequency and an enhancement in oxygenate yield. By applying a combination of atomically controlled synthesis, in situ characterization, and theoretical calculations, we gain an understanding of the promoter-Rh interactions that govern catalytic performance for MoO 3 -promoted Rh. We use atomic layer deposition to modify Rh nanoparticles with monolayer-precise amounts of MoO 3 , with a high degree of control over the structure of the catalyst. Through in situ X-ray absorption spectroscopy, we find that the atomic structure of the catalytic surface under reaction conditions consists of Mo-OH species substituted into the surface of the Rh nanoparticles. Using density functional theory calculations, we identify two roles of MoO 3 : first, the presence of Mo-OH in the catalyst surface enhances CO dissociation and also stabilizes a methanol synthesis pathway not present in the unpromoted catalyst; and second, hydrogen spillover from Mo-OH sites to adsorbed species on the Rh surface enhances hydrogenation rates of reaction intermediates.

We demonstrate that the possibility of monitoring relative photoionization cross sections over a large photon energy range allows us to study and disentangle shake processes and intramolecular inelastic scattering effects. In this gas-phase study, relative intensities of the carbon 1s photoelectron lines from chemically inequivalent carbon atoms in the same molecule have been measured as a function of the incident photon energy in the range of 300-6000 eV. We present relative cross sections for the chemically shifted carbon 1s lines in the photoelectron spectra of ethyl trifluoroacetate (the "ESCA" molecule). The results are compared with those of methyl trifluoroacetate and S -ethyl trifluorothioacetate as well as a series of chloro-substituted ethanes and 2-butyne. In the soft X-ray energy range, the cross sections show an extended X-ray absorption fine structure type of wiggles, as was previously observed for a series of chloroethanes. The oscillations are damped in the hard X-ray energy range, but deviations of cross-section ratios from stoichiometry persist, even at high energies. The current findings are supported by theoretical calculations based on a multiple scattering model. The use of soft and tender X-rays provides a more complete picture of the dominant processes accompanying photoionization. Such processes reduce the main photoelectron line intensities by 20-60%. Using both energy ranges enabled us to discern the process of intramolecular inelastic scattering of the outgoing electron, whose significance is otherwise difficult to assess for isolated molecules. This effect relates to the notion of the inelastic mean free path commonly used in photoemission studies of clusters and condensed matter.

Corvus, a Python-based package designed for managing workflows of physical simulations that utilize multiple scientific software packages, is presented. Corvus can be run as an executable script with an input file and automatically generated or custom workflows, or interactively, in order to build custom workflows with a set of Corvus-specific tools. Several prototypical examples are presented that link density functional, vibrational and X-ray spectroscopy software packages and are of interest to the synchrotron community. These examples highlight the simplification of complex spectroscopy calculations that were previously limited to expert users, and demonstrate the flexibility of the Corvus infrastructure to tackle more general problems in other research areas.

Extended x-ray absorption fine structure (EXAFS) is well-suited for investigations of structure and disorder of complex materials. Recently, experimental measurements and analysis of EXAFS have been carried out to elucidate the mechanisms responsible for the negative thermal expansion (NTE) in zirconium tungstate (ZrW 2 O 8 ). In contrast to previous work suggesting that transverse O-displacements are largely responsible, the EXAFS analysis suggested that correlated rotations and translations of octahedra and tetrahedra within the structure are a major source. In an effort to resolve this controversy, we have carried out ab initio calculations of the structure, lattice vibrations, and EXAFS of ZrW 2 O 8 based on real-space multiple-scattering calculations using the FEFF9 code and auxiliary calculations of structure and Debye-Waller factors. We find that the theoretical simulations are consistent with observed EXAFS, and show that both of the above mechanisms contribute to the dynamical structure of ZrW 2 O 8 .

We present a study of the origin of the negative thermal expansion (NTE) on ZrW 2 O 8 by combining an efficient approach for computing the dynamical matrix with the Lanczos algorithm for generating the phonon density of states in the quasi-harmonic approximation. The simulations show that the NTE arises primarily from the motion of the O-sublattice, and in particular, from the transverse motion of the O atoms in the W-O and W-O-Zr bonds. In the low frequency range these combine to keep the WO 4 tetrahedra rigid and induce internal distortions in the ZrO 6 octahedra. The force constants associated with these distortions become stronger with expansion, resulting in negative Grüneisen parameters and NTE from the low frequency modes that dominate the positive contributions from the high frequency modes. This leads us to propose an anharmonic, two-frequency Einstein model that quantitatively captures the NTE behavior.

Supported Pt nanocatalysts generally exhibit anomalous behavior, including negative thermal expansion and large structural disorder. Finite temperature DFT/MD simulations reproduce these properties, showing that they are largely explained by a combination of thermal vibrations and low-frequency disorder. We show here that a full interpretation is more complex and that the DFT/MD mean-square relative displacements (MSRD) can be further separated into vibrational disorder, "dynamic structural disorder" (DSD), and long-time equilibrium fluctuations of the structure dubbed "anomalous structural disorder" (ASD). We find that the vibrational and DSD components behave normally, increasing linearly with temperature while the ASD decreases, reflecting the evolution of mean nanoparticle geometry. As a consequence the usual procedure of fitting the MSRD to normal vibrations plus temperature-independent static disorder results in unphysical bond strengths and Grüneisen parameters.

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.

First principle calculations of extended X-ray absorption fine structure (EXAFS) data have seen widespread use in bioinorganic chemistry, perhaps most notably for modeling the Mn4Ca site in the oxygen evolving complex (OEC) of photosystem II (PSII). The logic implied by the calculations rests on the assumption that it is possible to a priori predict an accurate EXAFS spectrum provided that the underlying geometric structure is correct. The present study investigates the extent to which this is possible using state of the art EXAFS theory. The FEFF program is used to evaluate the ability of a multiple scattering-based approach to directly calculate the EXAFS spectrum of crystallographically defined model complexes. The results of these parameter free predictions are compared with the more traditional approach of fitting FEFF calculated spectra to experimental data. A series of seven crystallographically characterized Mn monomers and dimers is used as a test set. The largest deviations between the FEFF calculated EXAFS spectra and the experimental EXAFS spectra arise from the amplitudes. The amplitude errors result from a combination of errors in calculated S0(2) and Debye-Waller values as well as uncertainties in background subtraction. Additional errors may be attributed to structural parameters, particularly in cases where reliable high-resolution crystal structures are not available. Based on these investigations, the strengths and weaknesses of using first-principle EXAFS calculations as a predictive tool are discussed. We demonstrate that a range of DFT optimized structures of the OEC may all be considered consistent with experimental EXAFS data and that caution must be exercised when using EXAFS data to obtain topological arrangements of complex clusters.

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.

A transferable potential energy function for describing the interaction between water molecules is presented. The electrostatic interaction is described rigorously using a multipole expansion. Only one expansion center is used per molecule to avoid the introduction of monopoles. This single center approach turns out to converge and give close agreement with ab initio calculations when carried out up to and including the hexadecapole. Both dipole and quadrupole polarizability are included. All parameters in the electrostatic interaction as well as the dispersion interaction are taken from ab initio calculations or experimental measurements of a single water molecule. The repulsive part of the interaction is parametrized to fit ab initio calculations of small water clusters and experimental measurements of ice Ih. The parametrized potential function was then used to simulate liquid water and the results agree well with experiment, even better than simulations using some of the point charge potentials fitted to liquid water. The evaluation of the new interaction potential for condensed phases is fast because point charges are not present and the interaction can, to a good approximation, be truncated at a finite range.

Carbon 1s photoelectron spectra for 2-butyne (CH3C≡CCH3) measured in the photon energy range from threshold to 150 eV above threshold show oscillations in the intensity ratio C2,3/C1,4. Similar oscillations have been seen in chloroethanes, where the effect has been attributed to EXAFS-type scattering from the substituent chlorine atoms. In 2-butyne, however, there is no high-Z atom to provide a scattering center and, hence, oscillations of the magnitude observed are surprising. The results have been analyzed in terms of two different theoretical models: a density-functional model with B-spline atom-centered functions to represent the continuum electrons and a multiple-scattering model using muffin-tin potentials to represent the scattering centers. Both methods give a reasonable description of the energy dependence of the intensity ratios.

Recent studies of the crystallization of cyclotrimethylene-trinitramine (RDX) have shown that the presence of the α- and β-phases of the compound is sensitive to the substrate when using drop cast crystallization methods. The specific phase has potential consequences for measurements of the nitrogen K X-ray emission spectrum (XES) that were recently reported for this compound using samples crystallized on In metal substrates. We have determined that the crystallization of RDX on a clean In metal substrate starts out completely as the β-phase but progressively incorporates the α-phase as the film thickens. In addition, we have carried out additional molecular orbital calculations of the N 1s X-ray fluorescence from the valence band, comparing the results expected from the α-and β- phases. The differences due to the presence of the β-phase instead of, or in addition to, the α-phase appear to be minimal.

X-ray photoemission spectroscopy is used in a great variety of research fields; one observable is the sample's stoichiometry. The stoichiometry can be deduced based on the expectation that the ionization cross sections for innershell orbitals are independent of the molecular composition. Here we used chlorine-substituted ethanes in the gas phase to investigate the apparent carbon stoichiometry. We observe a nonstoichiometric ratio for a wide range of photon energies, the ratio exhibits x-ray-absorption fine structure spectroscopy (EXAFS)-like oscillations and hundreds of eV above the C1s ionization approaches a value far from 1. These effects can be accounted for by considering the scattering of the outgoing photoelectron, which we model by multiple-scattering EXAFS calculations, and by considering the effects of losses due to monopole shakeup and shakeoff and to intramolecular inelastic scattering processes.

Nonresonant X-ray emission spectroscopy was used to compare the nitrogen-rich compounds ammonium nitrate, trinitrotoluene, and cyclotrimethylene-trinitramine. They are representative of crystalline and molecular structures of special importance in industrial and military applications. The spectral signature of each substance was analyzed and correlated with features in the electronic structure of the systems. This analysis was accomplished by means of theoretical simulations of the emission spectra and a detailed examination of the molecular orbitals and densities of states. We find that the two theoretical methods used (frozen-orbital density functional theory and real-space Green's function simulations) account semiquantitatively for the observed spectra and are able to predict features arising from distinct chemical complexes. A comparison of the calculations and the data provides insight into the electronic contributions of specific molecular orbitals, as well as the features due to bandlike behavior. With some additional refinements, these methods could be used as an alternative to reference compounds.

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.

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.

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 report nonresonant inelastic X-ray scattering (NRIXS) measurement of core-shell excitations from both B 1s and C 1s initial states in all three isomers of the dicarba-closo-dodecarboranes C2B10H12. First, these data yield an experimental determination of the angular-momentum-projected final local density of states (l-DOS). We find low-energy resonances with distinctive local s- or p-type character, providing a more complete experimental characterization of bond hybridization than is available from dipole-transition limited techniques, such as X-ray absorption spectroscopies. This analysis is supported by independent density functional theory and real-space full multiple scattering calculation of the l-DOS which yield a clear distinction between tangential and radial contributions. Second, we investigate the isomer sensitivity of the NRIXS signal and compare and contrast these results with prior electron energy loss spectroscopy measurements. This work establishes NRIXS as a valuable tool for borane chemistry, not only for the unique spectroscopic capabilities of the technique but also through its compatibility with future studies in solution or in high-pressure environments. In addition, this work also establishes the real-space full multiple scattering approach as a useful alternative to traditional approaches for excited states calculations of aromatic polyhedral boranes and related systems.

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.

The prediction of nonlinear electro-optic (EO) behavior of molecules with quantum methods is the first step in the development of organic-based electro-optic devices. Typical EO molecules may require calculations with several hundred electrons, which prevents all but the fastest methods (semiempirical and density functional theory (DFT)) from being used for EO estimation. To test the reliability of these methods, we compare dipole moments, polarizabilities, and first-order hyperpolarizabilities for a wide range of structures of experimental interest with Hartree-Fock (HF), intermediate neglect of differential overlap (INDO), and DFT methods. The relative merits of molecules are consistently predictable with every method.