CO2 reforming of methane, also known as dry reforming of methane, offers promising prospects for the catalytic conversion of CO2 and CH4 to produce H2, but low-cost catalysts with high activity and stability in operation under severe conditions are needed. We now report such a catalyst, consisting of copper-nickel nanoparticles on a support. The catalyst operated with high activity for more than 1200 h of continuous operation in a flow reactor with a feed of CH4:CO2:N2 in a molar ratio of 1:1:1. In-situ X-ray absorption spectra show that the catalyst consisted of a dilute alloy of copper in nickel during operation, being resistant to metal sintering and the formation of carbonaceous deposits even at 700 °C. The data are consistent with the suggestion that the catalyst was a single-atom alloy and point to the possibility of a range of applications of such catalysts.

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function approach [Rehr et al., J. Chem. Phys. 152, 174113 (2020)]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz, ensuring the exactness in the expansion limit. In addition, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green's function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green's function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green's function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green's function approach and other Green's function approaches in the numerical exact limit.

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.

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.

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.

Theoretical descriptions of excited states of molecular systems in high-energy regimes are crucial for supporting and driving many experimental efforts at light source facilities. However, capturing their complicated correlation effects requires formalisms that provide a hierarchical infrastructure of approximations. These approximations lead to an increased overhead in classical computing methods, and therefore, decisions regarding the ranking of approximations and the quality of results must be made on purely numerical grounds. The emergence of quantum computing methods has the potential to change this situation. In this study, we demonstrate the efficiency of Quantum Phase Estimator (QPE) in identifying core-level states relevant to x-ray photoelectron spectroscopy. We compare and validate the QPE predictions with exact diagonalization and real-time equation-of-motion coupled cluster formulations, which are some of the most accurate methods for states dominated by collective correlation effects.

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 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 & horbar;linear, planar, and cubic & horbar;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.

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.

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.

Gallium-based heterogeneous catalysts originally developed for commercial use in propane dehydroaromatization have been explored extensively as potential replacements for Pt- and Cr-based catalysts used in propane dehydrogenation for on-demand propylene production. In a large number of experimental and theoretical studies, active sites with a variety of Ga nuclearities, coordination environments, and oxidation states have been proposed. Isolated Ga(I) ions are often invoked, despite the scarcity of well-defined molecular species and their well-documented instability. In this study, we investigate the appearance of a high-intensity, low-energy white line at the Ga K-edge upon reduction of Ga/HZSM-5 by H2 at a temperature of ca. 500 degrees C, accompanied by a dramatic reduction in the Ga/gamma-Al2O3 does not show such behavior. In order to lay a rigorous foundation for characterizing these types of systems and to establish experimental signatures for the elusive Ga(I) oxidation state, we recorded Ga K-edge X-ray absorption spectra [including high-energy-resolution fluorescence detection-X-ray absorption near-edge spectroscopy (HERFD-XANES)] for several well-defined molecular and crystalline Ga(I) compounds. XANES is essential to establishing the presence of Ga(I), despite the overlap in edge energies with organoGa(III) compounds, because Ga(I)-containing oxide materials show very weak EXAFS scattering. Compared to the XANES of trigonal Ga(III)-containing materials, Ga(I) spectra display a significantly more intense white line feature. Theoretical simulations agree well with this experimental observation and reveal that the strong XANES intensity originates from the superposition of transitions to several empty, nearly degenerate p-like states. These signatures provide compelling evidence for assigning the intense white line and dramatic loss of EXAFS intensity in Ga/HZSM-5 to the near-quantitative reduction of Ga(III) to Ga(I), while the weaker white line and conventional EXAFS signal of Ga/gamma-Al2O3 point to, at most, a minor fraction of Ga(I) sites.

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.

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

We discuss the real-space Green's function (RSGF) approach for core-level photoelectron diffraction (PD) and its comparison to x-ray absorption spectra (XAS). We focus especially on the separable Green's function formalism for efficient calculations of high-order multiple scattering. Computational details such as scattering potentials, self-energy effects, Debye–Waller factors, and inelastic losses are discussed. Finally, we consider prospects for improving ab initio PD calculations.

Many interesting properties of functional materials, such as dynamic response and thermodynamic behavior, depend on their excited state properties. These functional properties are often related to excitations in the system, such as phonons and plasmons, which lead to inelastic losses, lifetime, and other dynamic effects. The excitations are pure many-body correlation effects that are missing from independent particle theories. They are revealed in x-ray spectra such as photoemission and absorption, where they show up as satellites beyond the quasi-particle approximation. Our main focus in this work is the use of Green's function methods to describe these effects. In particular, we discuss how the cumulant Green's function provides a unified treatment of such dynamic correlation effects in many contexts. Besides a robust theoretical framework, these methods also yield widely applicable tools for practical calculations of many functional properties of materials. This methodology is illustrated with a number of applications ranging from optical and x-ray spectra to thermodynamic properties, and dynamic response. Some recent extensions for more correlated systems are also briefly discussed.

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.

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.

Optical constants characterize the interaction of materials with light and are important properties in material design. Here we present a Python-based Corvus workflow for simulations of full spectrum optical constants from the UV-VIS to hard x-ray wavelengths based on the real-space Green's function code FEFF10 and structural data from the Materials Project (MP). The Corvus workflow manager and its associated tools provide an interface to FEFF10 and the MP database. The workflow parallelizes the FEFF computations of optical constants over all absorption edges for each material in the MP database specified by a unique MP-ID. The workflow tools determine the distribution of computational resources needed for that case. Similarly, the optical constants for selected sets of materials can be computed in a single-shot. To illustrate the approach, we present results for nearly all elemental solids in the periodic table, as well as a sample compound, and compared with experimental results. As in x-ray absorption spectra, these results are interpreted in terms of an atomic-like background and fine-structure contributions.

X-ray absorption spectroscopy (XAS) [extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES)] is a key technique within the heterogeneous catalysis community to probe the structure and properties of the active site(s) for a diverse range of catalytic materials. However, the interpretation of the raw experimental data to derive an atomistic picture of the catalyst requires modeling and analysis; the EXAFS data are compared to a model, and a goodness of fit parameter is used to judge the best fit. This EXAFS modeling can often be nontrivial and time-consuming; overcoming or improving these limitations remains a central challenge for the community. Considering these limitations, this Perspective highlights how recent developments in analysis software, increased availability of reliable computational models, and application of data science tools can be used to improve the speed, accuracy, and reliability of EXAFS interpretation. In particular, we emphasize the advantages of combining theory and EXAFS as a unified technique that should be treated as a standard (when applicable) to identify catalytic sites and not two separate complementary methods. Building on the recent trends in the computational catalysis community, we also present a community-driven approach to adopt FAIR Guiding Principles for the collection, analysis, dissemination, and storage of XAS data. Written with both the experimental and theory audience in mind, we provide a unified roadmap to foster collaborations between the two communities.