Dense CO Adlayers as Enablers of CO Hydrogenation Turnovers on Ru Surfaces

Dynamics of Elementary Steps on Metal Surfaces at High Coverages: The Prevalence and Kinetic Competence of Contiguous Bare-Atom Ensembles

The rate of elementary steps on densely-covered surfaces depends sensitively on repulsive interactions within dense adlayers, situations ubiquitous in practice and with kinetic consequences seldom captured by Langmuirian treatments of surface catalysis. This study develops an ensemble-based method that assesses how such repulsion influences the prevalence and kinetic competence of bare-atom ensembles of different size. Chemisorbed CO (CO*) is used as an example because it forms dense adlayers on metal nanoparticles during CO2 hydrogenation (CO2-H2) and other reactions, leading to significant repulsion that weakens the binding of CO* and kinetically-relevant transition states (TS). This approach is enabled by density functional theory and probability formalisms and describes the prevalence of ensembles of contiguous bare atoms from their formation energy (via CO* desorption); it then determines their competence in stabilizing the TS and mediating the reaction rates. The specific conclusions reflect the extent to which a given TS and CO* desorbed to form bare ensembles "sense" repulsion and the contribution of each ensemble size to each reaction channel mediated by distinct TS structures. These formalisms are illustrated by assessing the relative contributions, kinetic relevance, and ensemble size requirements for two CO2-H2 routes (direct and H-assisted CO2 activation to CO and H2O) on Ru nanoparticles, but they are not restricted to specific bound species or reaction channels. This method is essential to assess the kinetic relevance of elementary steps in a given catalytic sequence and to determine the contributions from parallel reaction channels at the crowded surfaces that prevail in the practice of surface catalysis.

Recycling Valuable Alkylbenzenes from Polystyrene through Methanol-Assisted Depolymerization

The vast bulk of polystyrene (PS), a major type of plastic polymers, ends up in landfills, which takes up to thousands of years to decompose in nature. Chemical recycling promises to enable lower-energy pathways and minimal environmental impacts compared with traditional incineration and mechanical recycling. Herein, we demonstrated that methanol as a hydrogen supplier assisted the depolymerization of PS (denoted as PS-MAD) into alkylbenzenes over a heterogeneous catalyst composed of Ru nanoparticles on SiO2. PS-MAD achieved a high yield of liquid products which accounted for 93.2 wt % of virgin PS at 280 °C for 6 h with the production rate of 118.1 mmolcarbon gcatal. -1 h-1. The major components were valuable alkylbenzenes (monocyclic aromatics and diphenyl alkanes), the sum of which occupied 84.3 wt % of liquid products. According to mechanistic studies, methanol decomposition dominates the hydrogen supply during PS-MAD, thereby restraining PS aromatization which generates by-products of fused polycyclic arenes and polyphenylenes.

Mechanistic Connections between CO2 and CO Hydrogenation on Dispersed Ruthenium Nanoparticles

Catalytic routes for upgrading CO₂ to CO and hydrocarbons have been studied for decades, and yet the mechanistic details and structure-function relationships that control catalytic performance have remained unresolved. This study elucidates the elementary steps that mediate these reactions and examines them within the context of the established mechanism for CO hydrogenation to resolve the persistent discrepancies and to demonstrate inextricable links between CO₂ and CO hydrogenation on dispersed Ru nanoparticles (6-12 nm mean diameter, 573 K). The formation of CH₄ from both CO₂-H₂ and CO-H₂ reactants requires the cleavage of strong C≡O bonds in chemisorbed CO, formed as an intermediate in both reactions, via hydrogen-assisted activation pathways. The C═O bonds in CO₂ are cleaved via direct interactions with exposed Ru atoms in elementary steps that are shown to be facile by fast isotopic scrambling of C¹⁶O₂-C¹⁸O₂-H₂ mixtures. Such CO₂ activation steps form bound CO molecules and O atoms; the latter are removed via H-addition steps to form H₂O. The kinetic hurdles in forming CH₄ from CO₂ do not reflect the inertness of C═O bonds in CO₂ but instead reflect the intermediate formation of CO molecules, which contain stronger C≡O bonds than CO₂ and are present at near-saturation coverages during CO₂ and CO hydrogenation catalysis. The conclusions presented herein are informed by a combination of spectroscopic, isotopic, and kinetic measurements coupled with the use of analysis methods that account for strong rate inhibition by chemisorbed CO. Such methods enable the assessment of intrinsic reaction rates and are essential to accurately determine the effects of nanoparticle structure and composition on reactivity and selectivity for CO₂-H₂ reactions.

A combined theoretical/experimental study highlighting the formation of carbides on Ru nanoparticles during CO hydrogenation

Formation of stable carbides during CO bond dissociation on small ruthenium nanoparticles (RuNPs) is demonstrated, both by means of DFT calculations and by solid state ¹³C NMR techniques. Theoretical calculations of chemical shifts in several model clusters are employed in order to secure experimental spectroscopic assignations for surface ruthenium carbides. Mechanistic DFT investigations, carried out on a realistic Ru₅₅ nanoparticle model (∼1 nm) in terms of size, structure and surface composition, reveal that ruthenium carbides are obtained during CO hydrogenation. Calculations also indicate that carbide formation via hydrogen-assisted hydroxymethylidyne (COH) pathways is exothermic and occurs at reasonable kinetic cost on standard sites of the RuNPs, such as 4-fold ones on flat terraces, and not only in steps as previously suggested. Another novel outcome of the DFT mechanistic study consists of the possible formation of μ₆ ruthenium carbides in the tip-B₅ site, similar examples being known only for molecular ruthenium clusters. Moreover, based on DFT energies, the possible rearrangement of the surface metal atoms around the same tip-site results in a μ-Ru atom coordinated to the remaining RuNP moiety, reminiscent of a pseudo-octahedral metal center on the NP surface.

Kinetic effects of molecular clustering and solvation by extended networks in zeolite acid catalysis

Reactions catalyzed within porous inorganic and organic materials and at electrochemical interfaces commonly occur at high coverage and in condensed media, causing turnover rates to depend strongly on interfacial structure and composition, collectively referred to as "solvent effects". Transition state theory treatments define how solvation phenomena enter kinetic rate expressions, and identify two distinct types of solvent effects that originate from molecular clustering and from the solvation of such clusters by extended solvent networks. We review examples from the recent literature that investigate reactions within microporous zeolite catalysts to illustrate these concepts, and provide a critical appraisal of open questions in the field where future research can aid in developing new chemistry and catalyst design principles.

Structure and solvation of confined water and water-ethanol clusters within microporous Brønsted acids and their effects on ethanol dehydration catalysis

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 S_N2 transition states that eliminate DEE relative to (C₂H₅OH)(H⁺)(H₂O)_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 H₂O), 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 H₂O), water–ethanol clusters reach their maximum stable size ((C₂H₅OH)(H⁺)(H₂O)_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.

Water Poisons H2 Activation at the Au-TiO2 Interface by Slowing Proton and Electron Transfer between Au and Titania

Understanding the dynamic changes at the active site during catalysis is a fundamental challenge that promises to improve catalytic properties. While performing Arrhenius studies during H₂ oxidation over Au/TiO₂ catalysts, we found different apparent activation energies (E_app) depending on the feedwater pressure. This is partially attributed to changing numbers of metal-support interface (MSI) sites as water coverage changes with temperature. Constant water coverage studies showed two kinetic regimes: fast heterolytic H₂ activation directly at the MSI (E_app ∼ 25 kJ/mol) and significantly slower heterolytic H₂ activation mediated by water (E_app ∼ 45 kJ/mol). The two regimes had significantly different kinetics, suggesting a complicated mechanism of water poisoning. Density functional theory (DFT) showed water has minor effects on the reaction thermodynamics, primarily attributable to intrinsic differences in surface reactivity of different Au sites in the DFT model. The DFT model suggested significant surface restructuring of the TiO₂ support during heterolytic H₂ adsorption; evidence for this phenomenon was observed during in situ infrared spectroscopy experiments. A monolayer of water on the hydroxylated TiO₂ surface increased the H₂ dissociation activation barrier by ∼0.2 eV, in good agreement the difference in experimentally measured values. DFT calculations suggested H₂ activation goes through a proton-coupled electron-transfer-like mechanism. During proton transfer to a basic support hydroxyl group, electron density is distributed through the gold nanorod and partially localized on the protonated support hydroxyl group. Water slows H₂ activation by slowing this H⁺ transfer, forcing negative charge buildup on the Au and increasing the transition state energy.

How coverage influences thermodynamic and kinetic isotope effects for H2/D2 dissociative adsorption on transition metals

Hydrogen isotope effects are influenced by adsorbate coverage: at high coverages, isotope effects are lower than at low coverages. This helps to rationalize observed isotope effects, allowing more precise elucidation of reaction mechanisms.

Similarities and differences for atomic and diatomic molecule adsorption on the B-5 type sites of the HCP(101̅6) surfaces of Co, Os, and Ru from DFT calculations

The differences in relative adsorption energies for mono-atomic and diatomic prototype species (C,N,O,S,H,CO,NO,SO,CH,NH,H₂,O₂) relevant to catalytic processes such as Fischer-Tropsch and Ammonia Synthesis chemistry are investigated on the previously un-studied (101¯6) surface(s) of Co, Os, and Ru. Recent work in the literature has confirmed that catalytically relevant nanoparticles of HCP elements such as Co, Os, and Ru typically possess highly active ‘B5’ sites; unfortunately many early and extant theory and model-ing treatments of "stepped HCP surfaces" use ad-hoc created steps via manual deletion of atoms from an ideal HCP(0001) slab model. To date the differences in adsorption energies at various B5 step edge types, and any possible trends across the same type of B5 sites on various HCP catalyst species has not been thoroughly characterized. Our work in this manuscript uses the low energy (101¯6) Miller Index surface of Co, Os, and Ru which exposes 2 distinct and strongly adsorbing step edge sites, the B5B and B5A step edge which have been reported as relevant in the literature for Cobalt nanoparticle catalysis applications. Results from this study should be used to help further understand atomistic processes on the stepped surfaces of catalytically active HCP elements.

The consequences of surface heterogeneity of cobalt nanoparticles on the kinetics of CO methanation

The CO hydrogenation reaction was studied under methanation conditions (H₂/CO >3, 250–300 °C) on Co/SiO₂ catalysts with different mean Co nanoparticle size (d_p = 4 nm, 13 nm and 33 nm).

Microkinetic Analysis and Scaling Relations for Catalyst Design

Microkinetic analysis plays an important role in catalyst design because it provides insight into the fundamental surface chemistry that controls catalyst performance. In this review, we summarize the development of microkinetic models and the inclusion of scaling relationships in these models. We discuss the importance of achieving stoichiometric and thermodynamic consistency in developing microkinetic models. We also outline how analysis of the maximum rates of elementary steps can be used to determine which transition states and adsorbed intermediates are kinetically significant, allowing the derivation of general reaction kinetics rate expressions in terms of changes in binding energies of the relevant transition states and intermediates. Through these analyses, we present how to predict optimal surface coverages and binding energies of adsorbed species, as well as the extent of potential rate improvement for a catalytic system. For systems in which the extent of potential rate improvement is small because of limitations imposed by scaling relations, different approaches, including the addition of promoters and formation of catalysts containing multiple functionalities, can be used to break the scaling relations and obtain further rate enhancement.