Search

The d electronic states of Pt3M intermetallic alloys, probed by Resonant Inelastic X-ray scattering and Density Functional Theory, reveal new descriptors for molecular adsorption.

The interplay between covalent and non-covalent interactions at the solid-liquid interface strongly influences electrocatalytic reactions. Although methods to determine the former interactions have been rigorously developed, the latter are often described with static bilayer models or similarly approximate methods. In this study, we account for disorder and dynamics at complex electrochemical interfaces by proposing a simple theory to estimate the enthalpy of solvation for adsorbed intermediates. In a strategy reminiscent of Born-Haber cycles, the enthalpy of solvation is expressed in terms of two sub-processes: vacancy creation by water reorganization, and adsorbate interaction with solvent molecules. The magnitude of the solvation enthalpy is then determined as a mean value from statistical sampling of hydrated adsorbate-catalyst configurations obtained from simulated annealing with ab initio molecular dynamic (AIMD) simulations. This theory is generalizable for many combinations of surfaces and adsorbates. Its application improves treatment of energetics at electrified solid-liquid interfaces, as well as corresponding structure-activity-stability predictions, as demonstrated for electrochemical oxygen reduction on Pt(111) and on Fe-N-C catalysts and for ethanol electrooxidation on Pt(111).

Efficient and direct conversion of methane to value-added products has been a long-term challenge in shale gas applications. Here, we show that atomically thin nanolayers of Pt with a single or double atomic layer thickness, supported on a two-dimensional molybdenum titanium carbide (MXene), catalyse non-oxidative coupling of methane to ethane/ethylene (C2). Kinetic and theoretical studies, combined with in-situ spectroscopic and microscopic characterizations, demonstrate that Pt nanolayers anchored at the hexagonal close-packed sites of the MXene support can activate the first C–H bond of methane to form methyl radicals that favour desorption over further dehydrogenation and thus suppress coke deposition. At 750 °C and 7% methane conversion, the catalyst runs for 72 hours of continuous operation without deactivation and exhibits >98% selectivity towards C2 products, with a turnover frequency of 0.2–0.6 s−1. Our findings provide insights into the design of highly active and stable catalysts for methane activation and create a platform for developing atomically thin supported metal catalysts. The challenge in non-oxidative coupling of methane lies in the activation of the first C–H bond while avoiding further dehydrogenations, which lead to the formation of coke. Here, atomically thin platinum nanolayers on two-dimensional molybdenum titanium carbides are reported as a superior catalyst for this reaction owing to reduced coke formation.

Layered double hydroxides (LDHs) are among the most active and studied catalysts for the oxygen evolution reaction (OER) in alkaline electrolytes. However, previous studies have generally either focused on a small number of LDHs, applied synthetic routes with limited structural control, or used non‐intrinsic activity metrics, thus hampering the construction of consistent structure–activity‐relations. Herein, by employing new individually developed synthesis strategies with atomic structural control, we obtained a broad series of crystalline α‐MA(II)MB(III) LDH and β‐MA(OH)2 electrocatalysts (MA=Ni, Co, and MB=Co, Fe, Mn). We further derived their intrinsic activity through electrochemical active surface area normalization, yielding the trend NiFe LDH > CoFe LDH > Fe‐free Co‐containing catalysts > Fe‐Co‐free Ni‐based catalysts. Our theoretical reactivity analysis revealed that these intrinsic activity trends originate from the dual‐metal‐site nature of the reaction centers, which lead to composition‐dependent synergies and diverse scaling relationships that may be used to design catalysts with improved performance., Catalytic activities for oxygen evolution on crystalline 3d transition metal layered double hydroxides are derived using electrochemical surface area based normalization. Density functional calculations reveal a dual‐metal‐site feature of the reaction centers that provides opportunities to design new catalysts with improved performance.

The structure and mobility of gold nanoclusters entrained in the titanosilicate zeolite TS-1 were investigated using ab initio molecular dynamics to characterize the stability and diffusion of cata...

Alkali metal–O2 batteries (i.e., Li/Na–O2) with high specific energies are promising alternatives to state-of-the-art metal-ion batteries. However, they are plagued by challenges arising from the u...

The interaction between catalyst and support is well known to influence the reactivity and stability of heterogeneous catalysts, and electrochemical hydrogen evolution catalysts based on amorphous ...

A series of cobalt-promoted Pt catalysts supported on multi-walled carbon nanotubes was synthesized, and their performance was evaluated for the water-gas shift (WGS) reaction. Compared to the monometallic Pt catalyst, the WGS turnover rate (TOR) at 300 °C was promoted by a factor of 10 at a Pt:Co molar ratio of 1:3. X-ray absorption spectroscopy and XRD showed the presence of a Pt3Co alloy along with a partially oxidized cobalt and a free cobalt metal phase after reduction pretreatment. In order to determine the dominant active site over the Co-promoted catalysts, selective leaching of partially oxidized Co (designated as CoOxHy) and Co metal was performed with a 5% wt. acetic acid solution, while preserving the Pt-rich phases. The WGS TOR at 300 °C for the Co-promoted catalysts after leaching was observed to be even lower than that of the monometallic Pt catalyst. Thus, the alloy formation was determined to be inconsequential towards promotion in the WGS TOR, while the dominant active site was determined to be a PtCo alloy in intimate contact with the CoOxHy phase. Combined Density Functional Theory (DFT) calculations and ab-initio thermodynamic phase diagrams point to a monolayer of CoOH being the most stable Co phase on Pt(1 1 1) under WGS conditions. Calculations of OH binding energies on Pt(1 1 1), Pt3Co(1 1 1), and at the interface between CoOH overlayers and Pt(1 1 1) show that trends in the WGS activity of these catalysts are linked to the strength of OH binding, with the strongest OH binding found at the interface between CoOH and Pt, supporting the conclusion that a similar interface is the source of enhanced WGS activity in the PtCo bimetallic system.

Active centers in porous solid catalysts are multifaceted in structure, comprised of primary sites that bind intermediates, secondary environments that confine intermediates and transition states, ...

The impact of metal/oxide interfaces on the catalytic properties of oxide-supported metal nanoparticles is a topic of longstanding interest in the field of heterogeneous catalysis. The significance of the metal/oxide interaction has been shown to vary according to both the inherent reactivity of the metal nanoparticle and the properties of the oxide support, with effects such as the metal d-band center, the nanoparticle shape, and the reducibility of the oxide believed to contribute to the overall system reactivity. In recent years, the water gas shift (WGS) reaction, wherein carbon monoxide and water are converted to carbon dioxide and hydrogen, has emerged as a model chemistry to probe the molecular-level details of how catalysis can be promoted in such environments, and this reaction is the focus of the present contribution. Using a combination of periodic Density Functional Theory calculations and microkinetic modeling, we present a comprehensive analysis of the WGS mechanism at the interface between a quasi-one dimensional platinum nanowire and an irreducible MgO support. The nanowire is lattice matched to the MgO support to remove spurious strain at the metal/oxide interface, and reactions both on the nanowire and at the three-phase boundary itself are considered in the mechanistic analysis. Additionally, to elucidate the consequences of adsorbate–adsorbate interactions on the WGS chemistry, an ab-initio thermodynamic analysis of CO coverage is performed, and the impact of the higher coverage CO states on the reaction chemistry is explicitly evaluated. These results are combined with detailed calculations of adsorbate entropies and dual-site microkinetic modeling to determine the kinetically significant features of the WGS reaction network which are subsequently, validated through experimental measurements of apparent reaction orders and activation barrier. The analysis demonstrates the important role that the metal/oxide interface plays in the reaction, with the water dissociation step being facile at the interface compared to the pure metal or oxide surfaces. Further, explicit consideration of CO interactions with other adsorbates at the metal/oxide interface is found to be central to correctly determining reaction mechanisms, rate determining steps, reaction orders, and effective activation barriers. These results are captured in a closed-form Langmuir–Hinshelwood model, derived from a simplified version of the complete microkinetic analysis, which reveals, among other results, that the celebrated carboxyl mechanism of Mavrikakis and coworkers is the governing pathway when accounting for reaction-relevant CO coverages.

Transition metal dichalcogenides (TMDCs) have garnered much attention recently due to their remarkable performance for different electrochemical systems. In this study, we report on the synthesis a...

Protective coatings have been shown to effectively suppress Mn ion dissolution from the spinel LiMn2O4 lithium-ion battery cathode by stabilizing the surface against undesired side reactions with t...

High-entropy nanoparticles have become a rapidly growing area of research in recent years. Because of their multielemental compositions and unique high-entropy mixing states (i.e., solid-solution) that can lead to tunable activity and enhanced stability, these nanoparticles have received notable attention for catalyst design and exploration. However, this strong potential is also accompanied by grand challenges originating from their vast compositional space and complex atomic structure, which hinder comprehensive exploration and fundamental understanding. Through a multidisciplinary view of synthesis, characterization, catalytic applications, high-throughput screening, and data-driven materials discovery, this review is dedicated to discussing the important progress of high-entropy nanoparticles and unveiling the critical needs for their future development for catalysis, energy, and sustainability applications.

Intermixing of atomic species at the electrode-electrolyte boundaries can impact the properties of the interfaces in solid-state batteries. Herein, this work uses first-principles statistical mechanics along with experimental characterization to understand intermixing at the electrode-electrolyte interface. For the model presented in this work, lithium manganese oxide (LiMn

Heterogeneous catalytic reactions are influenced by a subtle interplay of atomic-scale factors, ranging from the catalysts’ local morphology to the presence of high adsorbate coverages. Describing such phenomena via computational models requires generation and analysis of a large space of atomic configurations. To address this challenge, we present Adsorbate Chemical Environment-based Graph Convolution Neural Network (ACE-GCN), a screening workflow that accounts for atomistic configurations comprising diverse adsorbates, binding locations, coordination environments, and substrate morphologies. Using this workflow, we develop catalyst surface models for two illustrative systems: (i) NO adsorbed on a Pt3Sn(111) alloy surface, of interest for nitrate electroreduction processes, where high adsorbate coverages combined with low symmetry of the alloy substrate produce a large configurational space, and (ii) OH* adsorbed on a stepped Pt(221) facet, of relevance to the Oxygen Reduction Reaction, where configurational complexity results from the presence of irregular crystal surfaces, high adsorbate coverages, and directionally-dependent adsorbate-adsorbate interactions. In both cases, the ACE-GCN model, trained on a fraction (~10%) of the total DFT-relaxed configurations, successfully describes trends in the relative stabilities of unrelaxed atomic configurations sampled from a large configurational space. This approach is expected to accelerate development of rigorous descriptions of catalyst surfaces under in-situ conditions.

We demonstrate the synthesis of a Pt3V alloy and Pt/Pt3V core/shell catalysts, which are highly selective for propane dehydrogenation. The selectivity is a result of the Pt3V intermetallic phase, w...

Recent work has identified enhanced charge-storage capacity in the spinel lithium manganese oxide (LiMn2O4; LMO) lithium-ion battery cathode upon a single atomic layer deposition (ALD) cycle compos...

Tuning catalysis at solid–solid interfaces is critical for the development of next-generation energy storage devices such as Li–O2 batteries, where solid lithium–oxygen species are formed and dissociated on a solid catalyst. Herein, atomically controlled synthesis is combined with theoretical calculations, electrochemical studies, and detailed characterization measurements to show that the interface between an oxide catalyst and the solid products is key to selectively control discharge product distribution, consequently affecting charge overpotentials. A surface structure-dependent electrochemical performance for nonprecious metal-containing, nanostructured lanthanum nickelate oxide (La2NiO4+δ, LNO) electrocatalysts is demonstrated. LNO nanostructures with (001) NiO-terminated surfaces exhibit lower charge overpotentials, as opposed to irregularly terminated polyhedral-shaped oxides of the same composition. It is found that these LNO nanostructures, with controlled surface structure, enhance the performa...

Measurements of turnover rates of gas-phase bimolecular ethanol dehydration to diethyl ether (404–438 K) on a suite of hydrophobic and hydrophilic Sn-zeolites (Sn-Beta, Sn-BEC, Sn-MFI) of varying S...

Oxide conversion reactions are known to have substantially higher specific capacities than intercalation materials used in Li-ion batteries, but universally suffer from large overpotentials associated with the formation of interfaces between the resulting nanoscale metal and Li2O products. Here we use the interfacial sensitivity of operando X-ray reflectivity to visualize the structural evolution of ultrathin NiO electrodes and their interfaces during conversion. We observe two additional reactions prior to the well-known bulk, three-dimensional conversion occurring at 0.6 V: an accumulation of lithium at the buried metal/oxide interface (at 2.2 V) followed by interfacial lithiation of the buried NiO/Ni interface at the theoretical potential for conversion (at 1.9 V). To understand the mechanisms for bulk and interfacial lithiation, we calculate interfacial energies using density functional theory to build a potential-dependent nucleation model for conversion. These calculations show that the additional space charge layer of lithium is a crucial component for reducing energy barriers for conversion in NiO.

The discovery of more efficient, economical, and selective catalysts for oxidative dehydrogenation is of immense economic importance. However, the temperatures required for this reaction are typically high, often exceeding 400 °C. Herein, we report the discovery of subnanometer sized cobalt oxide clusters for oxidative dehydrogenation of cyclohexane that are active at lower temperatures than reported catalysts, while they can also eliminate the combustion channel. These results found for the two cluster sizes suggest other subnanometer size (CoO)x clusters will also be active at low temperatures. The high activity of the cobalt clusters can be understood on the basis of density functional studies that reveal highly active under-coordinated cobalt atoms in the clusters and show that the oxidized nature of the clusters substantially decreases the binding energy of the cyclohexene species which desorb from the cluster at low temperature., Current oxidative dehydrogenation processes are based on petroleum cracking that is indirect, environmentally unfriendly, and energy intensive. Here, the authors discover that subnanometer sized cobalt oxide clusters are active for oxidative dehydrogenation of cyclohexane at lower temperatures than reported catalysts.

The NiFe layered double hydroxides (LDHs) are among the most active electrocatalysts for the oxygen evolution reaction (OER) in alkaline electrolytes.1 The incorporation of Fe dramatically increases the catalytic activity of Ni hydroxides.2 However, even the most active NiFe LDH catalyst still shows a considerable overpotential for the OER. Understanding the nature of the catalytic active sites and the catalytic mechanism are key challenges to develop better OER electrocatalysts. In this contribution, atomic-scale details of the catalytic active phase will be presented, showing that NiFe LDHs are oxidized under applied anodic potentials from as-prepared α-phases to activated γ-phases.3 The interlayer and in-plane lattice parameters of the OER-active γ-phase were obtained by performing wide angle X-ray scattering (WAXS) measurements on NiFe LDH nanoplatelets during operating electrochemical conditions and were characterized by a contraction of about 8%. Operando WAXS was then performed for other selected catalysts belonging to the transition metal LDH family of materials. Structural similarities of the catalytically active phases will be highlighted. Finally, in order to derive activity-structure relationships, an approach is presented to calculate intrinsic catalytic activities, which are challenging to evaluate for this class of materials. The presented method is based on electrochemical active surface area normalization obtained by impedance spectroscopy measurements under OER.4 References F. Dionigi and P. Strasser, Advanced Energy Materials, 2016, 6 1600621. L. Trotochaud, S. L. Young, J. K. Ranney and S. W. Boettcher, Journal of the American Chemical Society, 2014, 136: 6744-6753. F. Dionigi, Z. Zeng, I. Sinev, T. Merzdorf, S. Deshpande, M. B. Lopez, S. Kunze, I. Zegkinoglou, H. Sarodnik, D. Fan, A. Bergmann, J. Drnec, J. Ferreira de Araujo, M. Gliech, D. Teschner, J. Zhu, W.-X. Li, J. Greeley, B. Roldan Cuenya and P. Strasser, Nat Commun, 2020, 11 2522. F. Dionigi, J. Zhu, Z. Zeng, T. Merzdorf, H. Sarodnik, M. Gliech, L. Pan, W.-X. , Li, J. Greeley, and P. Strasser, Angew Chem Int Edit, 60, 14446 – 14457 (2021).

Structural and chemical transformations of ultrathin oxide films on transition metals lie at the heart of many complex phenomena in heterogeneous catalysis, such as the strong metal-support interaction (SMSI). However, there is limited atomic-scale understanding of these transformations, especially for irreducible oxides such as ZnO. Here, by combining density functional theory calculations and surface science techniques, including scanning tunneling microscopy, X-ray photoelectron spectroscopy, high-resolution electron energy loss spectroscopy, and low-energy electron diffraction, we investigated the interfacial interaction of well-defined ultrathin ZnOxHy films on Pd(111) under varying gas-phase conditions [ultrahigh vacuum (UHV), 5 × 10-7 mbar of O2, and a D2/O2 mixture] to shed light on the SMSI effect of irreducible oxides. Sequential treatment of submonolayer zinc oxide films in a D2/O2 mixture (1:4) at 550 K evoked reversible structural transformations from a bilayer to a monolayer and further to a Pd-Zn near-surface alloy, demonstrating that zinc oxide, as an irreducible oxide, can spread on metal surfaces and show an SMSI-like behavior in the presence of hydrogen. A mixed canonical-grand canonical phase diagram was developed to bridge the gap between UHV conditions and true SMSI environments, revealing that, in addition to surface alloy formation, certain ZnOxHy films with stoichiometries that do not exist in bulk are stabilized by Pd in the presence of hydrogen. Based on the combined theoretical and experimental observations, we propose that SMSI metal nanoparticle encapsulation for irreducible oxide supports such as ZnO involves both surface (hydroxy)oxide and surface alloy formation, depending on the environmental conditions.

In heterogeneous catalysis, olefin oligomerization is typically performed on immobilized transition metal ions, such as Ni2+ and Cr3+. Here we report that silica-supported, single site catalysts containing immobilized, main group Zn2+ and Ga3+ ion sites catalyze ethylene and propylene oligomerization to an equilibrium distribution of linear olefins with rates similar to that of Ni2+. The molecular weight distribution of products formed on Zn2+ is similar to Ni2+, while Ga3+ forms higher molecular weight olefins. In situ spectroscopic and computational studies suggest that oligomerization unexpectedly occurs by the Cossee-Arlman mechanism via metal hydride and metal alkyl intermediates formed during olefin insertion and β-hydride elimination elementary steps. Initiation of the catalytic cycle is proposed to occur by heterolytic C-H dissociation of ethylene, which occurs at about 250 °C where oligomerization is catalytically relevant. This work illuminates new chemistry for main group metal catalysts with potential for development of new oligomerization processes. Silica-supported, single site, main group Zn(II) and Ga(III) ions catalyze ethylene and propylene oligomerization. Here, experimental and theoretical evidence suggests a Cossee-Arlman reaction mechanism similar to that for transition metal catalysts.

The development of industrial electrocatalysts with excellent performance in fuel cells strongly relies on the fundamental understanding of governing principles of catalysis from fundamental research on well-defined model systems. One particularly successful strategy has been to tune surface strain and thereby modulate the surface electronic structure and energetics of reaction intermediates, leading to a fundamental understanding and identification of catalysts with improved activity. For the well-known example of the oxygen reduction reaction (ORR) on platinum in fuel cells, previous studies indicated that a 1% compressive strain could improve the activity by more than 300%. While the concept of tailoring the reactivity of metal catalysts through surface strain has been discussed over two decades, fundamental challenges remain regarding both strain generation and strain tuning. In the present contribution, We will share our recent progress on these topics, including [1] tunable intrinsic strain in two-dimensional transition metal electrocatalysts, [2] tunable intrinsic strain in stepped transition metal electrocatalysts, and [3] strain tuning of supported transition metal overlayers. We will demonstrate [1] how these fundamental understanding can be used to understand terrace-width dependent ORR activity of stepped Pt surfaces that have puzzled the community over a decade; and [2] how they pave the way toward developing both model electrocatalysts and industrial electrocatalysts with significantly improved performance toward the oxygen reduction reaction. References Wang, Z. Zeng*†, W. Gao, T. Maxson, D. Raciti, M. Giroux, X. Pan, C. Wang, J. Greeley,. Science 2019, 363, 870-874. Wang, W. Gao, Z. Liu, Z. Zeng, Y. Liu, M. Giroux, M. Chi, G. Wang, J. Greeley, X. Pan, C. Wang, ACS Catal. 2018, 8, 35-42. J. Steinbach, C. Duru, A. T. Haug, A. E. Hester, M. Kuznia, K. A. Lewinski, S. M. Luopa, J. T. Petrin, G. M. Thoma, A. J. Kropf, D. J. Myers, D. Yang, D. A. Cullen, J. Greeley, Z. Zeng, ECS Transactions 2017, 80, 659-676. Zeng, J. Greeley, Nano Energy 2016, 29, 369.

Electrochemical water splitting using renewable energy is a key route to generate green H2 fuel. In this process, however, O2 generation at the anode through the oxygen evolution reaction (OER) is inherently slower by over four orders of magnitude than H2 generation at the cathode. Thus, improving OER efficiency has been a majority effort in electrolysis. Ni-based and Co-based layered double hydroxides (LDHs) are among the most active and studied catalysts for the oxygen evolution reaction (OER) in alkaline electrolytes. Because it happens under extremely oxidative aqueous conditions, however, in-situ crystal structures of the OER active phase are still largely unknown, significantly hindering the establishment of structure-property relationships. In this talk, we provide the first direct atomic-scale evidence that, under applied anodic potentials, NiFe and CoFe LDHs oxidize from as-prepared α-phases to activated γ-phases. The OER-active γ-phases are characterized by about 8% contraction of the lattice spacing and switching of the intercalated ions from carbonate to potassium. The calculated surface phase diagrams indicate that surface O sites are saturated with H by forming bridge OH, and coordinatively unsaturated metal sites are poisoned by OH adsorption under OER conditions. These structures, and the associated reaction free energies, suggest that the OER proceeds via a Mars van Krevelen mechanism, starting with the oxidation of bridge OH at the reaction centers with dual metal sites, i.e., M1-OH-M2. Our study suggests that the compound-dependent activity originates from the dual-metal site feature of the reaction centers. While this feature does not influence the OH-OOH scaling relationship, it leads to diverse OH-O scaling relationships, including those with near-zero slopes and negative slopes. Breaking OH-OOH scaling relationships were frequently discussed in the literature, as it determines the minimum overpotential. However, our study showed that, to approach the minimum overpotential dictated by a specific OH-OOH scaling relationship, the key is to break the OH-O scaling relationship. A possible route is to form binary metal oxyhydroxides with dual metal sites at the reaction centers or introduce a third element into NiFe LDH or CoFe LDH. References: Dionigi, Z. Zeng, I. Sinev, T. Merzdorf, S. Deshpande, M. B. Lopez, S. Kunze, I. Zegkinoglou, H. Sarodnik, D. Fan, A. Bergmann, J. Drnec, J. F. d. Araujo, M. Gliech, D. Teschner, J. Zhu, W.-X. Li, J. Greeley, B. R. Cuenya, P. Strasser, Nature Communications 2020, 11, 2522. Dionigi, J. Zhu, Z. Zeng, T. Merzdorf, H. Sarodnik, M. Gliech, L. Pan, W.-X. Li, J. Greeley, P. Strasser, Angew. Chem., Int. Ed. 2021, 60, 14446-14457.

In heterogeneous catalysis, olefin oligomerization is typically performed on immobilized transition metal ions, such as Ni2+ and Cr3+. Here we report that silica-supported, single site catalysts containing immobilized, main group Zn2+ and Ga3+ ion sites catalyze ethylene and propylene oligomerization to an equilibrium distribution of linear olefins with rates similar to that of Ni2+. The molecular weight distribution of products formed on Zn2+ is similar to Ni2+, while Ga3+ forms higher molecular weight olefins. In situ spectroscopic and computational studies suggest that oligomerization unexpectedly occurs by the Cossee-Arlman mechanism via metal hydride and metal alkyl intermediates formed during olefin insertion and β-hydride elimination elementary steps. Initiation of the catalytic cycle is proposed to occur by heterolytic C-H dissociation of ethylene, which occurs at about 250 °C where oligomerization is catalytically relevant. This work illuminates new chemistry for main group metal catalysts with potential for development of new oligomerization processes., Silica-supported, single site, main group Zn(II) and Ga(III) ions catalyze ethylene and propylene oligomerization. Here, experimental and theoretical evidence suggests a Cossee-Arlman reaction mechanism similar to that for transition metal catalysts.

The insertion and removal of Li

Heterogeneous catalysts constitute a crucial component of many industrial processes, and to gain an understanding of the atomic-scale features of such catalysts, ab initio density functional theory is widely employed. Recently, growing computational power has permitted the extension of such studies to complex reaction networks involving either high adsorbate coverages or multidentate adsorbates, which bind to the surface through multiple atoms. Describing all possible adsorbate configurations for such systems, however, is often not possible based on chemical intuition alone. To systematically treat such complexities, we present a generalized Python-based graph theory approach to convert atomic scale models into undirected graph representations. These representations, when combined with workflows such as evolutionary algorithms, can systematically generate high coverage adsorbate models and classify unique minimum energy multidentate adsorbate configurations for surfaces of low symmetry, including multi-elemental alloy surfaces, steps, and kinks. Two case studies are presented which demonstrate these capabilities; first, an analysis of a coverage-dependent phase diagram of absorbate NO on the Pt3Sn(111) terrace surface, and second, an investigation of adsorption energies, together with identifying unique minimum energy configurations, for the reaction intermediate propyne (CHCCH3*) adsorbed on a PdIn(021) step surface. The evolutionary algorithm approach reproduces high coverage configurations of NO on Pt3Sn(111) using only 15% of the number of simulations required for a brute force approach. Furthermore, the screening of potentially hundreds of multidentate adsorbates is shown to be possible without human intervention. The strategy presented is quite general and can be applied to a spectrum of complex atomic systems.

Water networks confined within zeolites solvate clustered reactive intermediates and must rearrange to accommodate transition states that differ in size and polarity, with thermodynamic penalties that depend on the shape of the confining environment., Aqueous-phase reactions within microporous Brønsted acids occur at active centers comprised of water-reactant-clustered hydronium ions, solvated within extended hydrogen-bonded water networks that tend to stabilize reactive intermediates and transition states differently. The effects of these diverse clustered and networked structures were disentangled here by measuring turnover rates of gas-phase ethanol dehydration to diethyl ether (DEE) on H-form zeolites as water pressure was increased to the point of intrapore condensation, causing protons to become solvated in larger clusters that subsequently become solvated by extended hydrogen-bonded water networks, according to in situ IR spectra. Measured first-order rate constants in ethanol quantify the stability of SN2 transition states that eliminate DEE relative to (C2H5OH)(H+)(H2O)n clusters of increasing molecularity, whose structures were respectively determined using metadynamics and ab initio molecular dynamics simulations. At low water pressures (2–10 kPa H2O), rate inhibition by water (–1 reaction order) reflects the need to displace one water by ethanol in the cluster en route to the DEE-formation transition state, which resides at the periphery of water–ethanol clusters. At higher water pressures (10–75 kPa H2O), water–ethanol clusters reach their maximum stable size ((C2H5OH)(H+)(H2O)4–5), and water begins to form extended hydrogen-bonded networks; concomitantly, rate inhibition by water (up to –3 reaction order) becomes stronger than expected from the molecularity of the reaction, reflecting the more extensive disruption of hydrogen bonds at DEE-formation transition states that contain an additional solvated non-polar ethyl group compared to the relevant reactant cluster, as described by non-ideal thermodynamic formalisms of reaction rates. Microporous voids of different hydrophilic binding site density (Beta; varying H+ and Si–OH density) and different size and shape (Beta, MFI, TON, CHA, AEI, FAU), influence the relative extents to which intermediates and transition states disrupt their confined water networks, which manifest as different kinetic orders of inhibition at high water pressures. The confinement of water within sub-nanometer spaces influences the structures and dynamics of the complexes and extended networks formed, and in turn their ability to accommodate the evolution in polarity and hydrogen-bonding capacity as reactive intermediates become transition states in Brønsted acid-catalyzed reactions.

Kinetic Monte Carlo simulations of electrochemical oxidation and reduction are presented that match many features of the experimentally observed electrochemical and morphological response of Pt(111). Included in the simulation are all relevant microscopic transitions, including the formation of Pt-OH and Pt-O from Pt, surface diffusion of all three species, as well as an effective place exchange diffusion at high potential. A detailed description of this approach to modeling such a complex surface is also presented. Overall, it is found that many features of the Pt(111) CV, including hydroxylation, hysteresis, and surface roughening, can be correlated to events associated with n-coordinated surface species, such as the hydroxylation wave corresponding to a one-electron oxidation of 9-coordinated terrace sites. Oxidation to Pt-O species at potentials above 1.0 V are shown to correlate to the presence of growing surface roughness, and the simulations suggest the onset of Pt-O formation in steady-state cyclic voltammetry is dominated by the oxidation of 8-coordinated step edges rather than terrace sites. Implications for the stability of Pt(111) catalysts after thousands of voltammetric cycles are discussed.

Artificial photosynthesis can be used to store solar energy and reduce CO2 into fuels to potentially alleviate global warming and the energy crisis. Compared to the generation of gaseous products, it remains a great challenge to tune the product distribution of artificial photosynthesis to liquid fuels, such as CH3 OH, which are suitable for storage and transport. Herein, we describe the introduction of metallic Cu nanoparticles (NPs) on Cu2 O films to change the product distribution from gaseous products on bare Cu2 O to predominantly CH3 OH by CO2 reduction in aqueous solutions. The specifically designed Cu/Cu2 O interfaces balance the binding strengths of H* and CO* intermediates, which play critical roles in CH3 OH production. With a TiO2 model photoanode to construct a photoelectrochemical cell, a Cu/Cu2 O dark cathode exhibited a Faradaic efficiency of up to 53.6 % for CH3 OH production. This work demonstrates the feasibility and mechanism of interface engineering to enhance the CH3 OH production from CO2 reduction in aqueous electrolytes.

Linear scaling relationships (SRs), which relate binding energies of adsorbates across a space of catalyst surfaces, have been extensively explored for metal and oxide surfaces, but little is known about their properties at interfaces between metal nanoparticles and oxide supports, which are ubiquitous in heterogeneous catalysis. Using periodic DFT calculations, scaling principles are extended to bifunctional Au/oxide interfaces. Adopting a Au nanorod on doped MgO (100) as a model, SRs for species participating in water gas shift, methanol synthesis, and oxidation reactions are reported. SR slopes are not constrained by the bond order conservation rule postulated for metals, oxides, and zeolites, potentially permitting greater flexibility in catalyst design strategies. The deviation from bond counting, along with the physical origin of scaling behavior at interfaces, are explored using a conceptual framework involving electrostatic interactions at the Au/oxide interface.

Supported multimetallic nanoparticles (NPs) are widely used in industrial catalytic processes, where the relation between surface structure and function is well-known. However, the effect of subsurface layers on such catalysts remains mostly unstudied. Here, we demonstrate a clear subsurface effect on supported 2 nm core-shell NPs with atomically precise and high temperature stable Pt3Mn intermetallic surface measured by in situ synchrotron X-ray Diffraction, difference X-ray Absorption Spectroscopy, and Energy Dispersive X-ray Spectroscopy. The NPs with a Pt3Mn subsurface have 98% selectivity to C-H over C-C bond activation during propane dehydrogenation at 550 °C compared with 82% for core-shell NPs with a Pt subsurface. The difference is correlated with significant reduction in the heats of reactant adsorption due to the Pt3Mn intermetallic subsurface as discerned by theory as well as experiment. The findings of this work highlight the importance of subsurface for supported NP catalysts, which can be tuned via controlled intermetallic formation. Such approach is generally applicable to modifying multimetallic NPs, adding another dimension to the tunability of their catalytic performance.