Dynamic structural evolution of MgO-supported palladium catalysts: from metal to metal oxide nanoparticles to surface then subsurface atomically dispersed cations.

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.

RT-EOM-CCSD Calculations of Inner and Outer Valence Ionization Energies and Spectral Functions.

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.

Impact of Local Structure in Supported CaO Catalysts for Soft-Oxidant-Assisted Methane Coupling Assessed through Ca K-Edge X-ray Absorption Spectroscopy.

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 (O2) 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 N2O-assisted methane coupling (N2O-OCM) and have previously been reported to be active for CO2-assisted methane coupling (CO2-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 N2O-OCM and previously reported CO2-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.

Cation Incorporation into Copper Oxide Lattice at Highly Oxidizing Potentials.

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 Ba2+ cations are directly incorporated into the lattice and form an ordered BaCuO2 phase at potentials more oxidizing than 200 mV vs the normal hydrogen electrode (NHE). BaCuO2 formation is followed by further oxidation to a bulk Cu3+-like BaxCuyOz 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.

Aqueous Structure of Lanthanide-EDTA Coordination Complexes Determined by a Combined DFT/EXAFS Approach.

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.

NWChem: Recent and Ongoing Developments.

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.

Real-Time Equation-of-Motion Coupled-Cluster Cumulant Green's Function Method: Heterogeneous Parallel Implementation Based on the Tensor Algebra for Many-Body Methods Infrastructure.

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.

Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer.

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.

Ab initio calculation of X-ray and related core-level spectroscopies: Green's function approaches.

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 Platinum in Surface and Subsurface Sites on MgO Have Contrasting Catalytic Properties for CO Oxidation.

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.

A Theory-Guided X-ray Absorption Spectroscopy Approach for Identifying Active Sites in Atomically Dispersed Transition-Metal Catalysts.

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.

Understanding Structure-Property Relationships of MoO3-Promoted Rh Catalysts for Syngas Conversion to Alcohols.

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 MoO3-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 MoO3-promoted Rh. We use atomic layer deposition to modify Rh nanoparticles with monolayer-precise amounts of MoO3, 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 MoO3: 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.

Energy-Dependent Relative Cross Sections in Carbon 1s Photoionization: Separation of Direct Shake and Inelastic Scattering Effects in Single Molecules.

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.

Extended X-Ray Absorption Fine Structure of ZrW2O8: Theory vs. Experiment.

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

Efficient Calculation of the Negative Thermal Expansion in ZrW2O8.

We present a study of the origin of the negative thermal expansion (NTE) on ZrW2O8 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 WO4 tetrahedra rigid and induce internal distortions in the ZrO6 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.

Anomalous Structural Disorder in Supported Pt Nanoparticles.

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.

How Accurately Can Extended X-ray Absorption Spectra Be Predicted from First Principles? Implications for Modeling the Oxygen-Evolving Complex in Photosystem II.

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.

Intensity oscillations in the carbon 1s ionization cross sections of 2-butyne.

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.

Surface effects on the crystallization of cyclo-1,3,5-trimethylene-2,4,6-trinitramine (RDX) and the consequences for its N K X-ray emission spectrum.

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 emission spectroscopy of nitrogen-rich compounds.

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.