Boron Switch for Selectivity of Catalytic Dehydrogenation on Size-Selected Pt Clusters on Al2O3

Dehydrogenation of diborane on small Nbn+ clusters

The reactivity of Nbn+ (1 ≤ n ≤ 21) clusters with B2H6 is studied by using a self-developed multiple-ion laminar flow tube reactor combined with a triple quadrupole mass spectrometer (MIFT-TQMS). The Nbn+ clusters were generated by a magnetron sputtering source and reacted with the B2H6 gas under fully thermalized conditions in the downstream flow tube where the reaction time was accurately controlled and adjustable. The complete and partial dehydrogenation products NbnB1-4+ and NbnB1-4H1,2,4+ were detected, indicative of the removal of H2 and likely BHx moieties. Interestingly, these NbnB1-4+ and NbnB1-4H1,2,4+ products are limited to 3 ≤ n ≤ 6, suggesting that the small Nbn+ clusters are relatively more reactive than the larger Nbn>6+ clusters under the same conditions. By varying the B2H6 gas concentrations and the reactant doses introduced into the flow tube, and by changing the reaction time, we performed a detailed analysis of the reaction dynamics in combination with the DFT-calculated thermodynamics. It is demonstrated that the lack of cooperative active sites on the Nb1+ cations accounts for the weakened dehydrogenation efficiency. Nb2+ forms partial dehydrogenation products at a faster rate. In contrast, the Nbn>6+ clusters are subject to more flexible vibrational relaxation which disperse the energy gain of B2H6-adsorption and thus are unable to overcome the energy barriers for subsequent hydrogen atom transfer and H2 release.

Photocatalytic non-oxidative conversion of methane

The direct conversion of methane to hydrogen and high-value hydrocarbons under mild conditions is an ideal, carbon-neutral method for utilizing natural gas resources. Compared with traditional high-temperature thermal catalytic methods, using clean light energy to activate inert C-H bonds in methane can not only significantly reduce the reaction temperature and avoid catalyst deactivation, but also surpass the limitations of thermodynamic equilibrium and provide new reaction pathways. This paper provides a comprehensive review of developments in the field of photocatalytic non-oxidative conversion of methane (PNOCM), while also highlighting our contributions, particularly focusing on catalyst design, product selectivity, and the underlying photophysical and chemical mechanisms. The challenges and potential solutions are also evaluated. The goal of this feature article is to establish a foundational understanding and stimulate further research in this emerging area.

H-Induced Restructuring on Cu(111) Triggers CO Electroreduction in an Acidic Electrolyte

In acidic conditions, the electroreduction of CO or CO2 (noted CO(2)RR) on metal surfaces is conventionally hindered by intense competition with the hydrogen evolution reaction (HER). In this study, we present first-principles calculations of a mechanism wherein the formation of H-induced Cu adatoms on Cu(111) serves as a pivotal trigger for CORR in acidic environments. Through an analysis of the grand canonical surface state population, we elucidate that these newly formed adatoms create an array of active sites essential for both CO adsorption and subsequent reduction. Our ensemble-based kinetic models unveil the role of adatoms, enhancing the HER while simultaneously initiating CORR. Notably, the cumulative activity of the HER and CORR is contingent upon the combination of various surface states, with their individual contributions varying based on the electrode potential and pH. The interplay between surface state dynamics and electrochemical activity sheds new light on the potential-dependent nature of the active site and reaction kinetics governing CORR on Cu(111) in acidic media.

Activity Origin of the Nickel Cluster on TiC Support for Nonoxidative Methane Conversion

Designing an active and selective catalyst for nonoxidative conversion of methane under mild conditions is critical for natural gas utilization as a chemical feedstock. Here, we demonstrate that the origin of the selective nonoxidative conversion of methane by the titanium carbide supported nickel cluster arises from the formation of a nickel carbide site under the reaction conditions, which could stabilize the CH_x intermediate to facilitate the C-C coupling, but further coking is rather limited. The reaction mechanism reveals that the C₂ products can be formed via a key -CH_x-CH₃ intermediate. In addition, we demonstrate that boration of the nickel cluster site can improve the methane conversion toward C₂ products. That higher activity and selectivity from the moderate rise in d orbital energy levels can therefore be considered as a descriptor of the catalyst effectiveness. These findings provide an understanding of the dynamic behavior of the single nickel cluster toward methane conversion to C₂ products and guidance for their future rational design.

Electrocatalytic Hydrogen Evolution at Full Atomic Utilization over ITO-Supported Sub-nano-Ptn Clusters: High, Size-Dependent Activity Controlled by Fluxional Pt Hydride Species

A combination of density functional theory (DFT) and experiments with atomically size-selected Pt_n clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Pt_n size on the electrocatalytic activity of Pt_n (n = 1, 4, 7, and 8) for the hydrogen evolution reaction (HER). Activity is found to be negligible for isolated Pt atoms on ITO, increasing rapidly with Pt_n size such that Pt₇/ITO and Pt₈/ITO have roughly double the activity per Pt atom compared to atoms in the surface layer of polycrystalline Pt. Both the DFT and experiment find that hydrogen under-potential deposition (H_upd) results in Pt_n/ITO (n = 4, 7, and 8) adsorbing ∼2H atoms/Pt atom at the HER threshold potential, equal to ca. double the H_upd observed for Pt bulk or nanoparticles. The cluster catalysts under electrocatalytic conditions are hence best described as a Pt hydride compound, significantly departing from a metallic Pt cluster. The exception is Pt₁/ITO, where H adsorption at the HER threshold potential is energetically unfavorable. The theory combines global optimization with grand canonical approaches for the influence of potential, uncovering the fact that several metastable structures contribute to the HER, changing with the applied potential. It is hence critical to include reactions of the ensemble of energetically accessible Pt_nH_x/ITO structures to correctly predict the activity vs Pt_n size and applied potential. For the small clusters, spillover of H_ads from the clusters to the ITO support is significant, resulting in a competing channel for loss of H_ads, particularly at slow potential scan rates.

Enhancing the Performance of Global Optimization of Platinum Cluster Structures by Transfer Learning in a Deep Neural Network

The global optimization of metal cluster structures is an important research field. The traditional deep neural network (T-DNN) global optimization method is a good way to find out the global minimum (GM) of metal cluster structures, but a large number of samples are required. We developed a new global optimization method which is the combination of the DNN and transfer learning (DNN-TL). The DNN-TL method transfers the DNN parameters of the small-sized cluster to the DNN of the large-sized cluster to greatly reduce the number of samples. For the global optimization of Pt₉ and Pt₁₃ clusters in this research, the T-DNN method requires about 3-10 times more samples than the DNN-TL method, and the DNN-TL method saves about 70-80% of time. We also found that the average amplitude of parameter changes in the T-DNN training is about 2 times larger than that in the DNN-TL training, which rationalizes the effectiveness of transfer learning. The average fitting errors of the DNN trained by the DNN-TL method can be even smaller than those by the T-DNN method because of the reliability of transfer learning. Finally, we successfully obtained the GM structures of Pt_n (n = 8-14) clusters by the DNN-TL method.

Enhancing the Catalytic Activity of Mo(110) Surface via Its Alloying with Submonolayer to Multilayer Boron Films and Oxidation of the Alloy: A Case of (CO + O2) to CO2 Conversion

In-situ formation of boron thin films on the Mo(110) surface, as well as the formation of the molybdenum boride and its oxide and the trends of carbon monoxide catalytic oxidation on the substrates formed, have been studied in an ultra-high vacuum (UHV) by a set of surface-sensitive characterization techniques: Auger and X-ray photoelectron spectroscopy (AES, XPS), low-energy ion scattering (LEIS), reflection-absorption infrared spectroscopy (RAIRS), temperature-programmed desorption (TPD), electron energy loss spectroscopy (EELS) and work function measurements using the Anderson method. The boron deposited at Mo(110) via electron-beam deposition at a substrate temperature of 300 K grows as a 2D layer, at least in submonolayer coverage. Such a film is bound to the Mo(110) via polarized chemisorption bonds, dramatically changing the charge density at the substrate surface manifested by the Mo(110) surface plasmon damping. Upon annealing of the B-Mo(110) system, the boron diffuses into the Mo(110) bulk following a two-mode regime: (1) quite easy dissolution, starting at a temperature of about 450 K with an activation energy of 0.4 eV; and (2) formation of molybdenum boride at a temperature higher than 700 K with M-B interatomic bonding energy of 3.8 eV. The feature of the formed molybdenum boride is that there is quite notable carbon monoxide oxidation activity on its surface. A further dramatic increase of such an activity is achieved when the molybdenum boride is oxidized. The latter is attributed to more activated states of molecular orbitals of coadsorbed carbon monoxide and oxygen due to their enhanced interaction with both boron and oxygen species for Mo_xB_yO_z ternary compound, compared to only boron for the Mo_x'B_y' double alloy.

Ultrasmall Cluster Model for Investigating Single Atom Catalysis: Dehydrogenation of 1-Propanamine by Size-Selected Pt1Zr2O7 Clusters Supported on HOPG

The selective dehydrogenation of hydrocarbons and their functionalized derivatives is a promising pathway in the realization of endothermic fuel systems for powering important technologies such as hypersonic aircraft. The recent surge in interest in single atom catalysts (SACs) over the past decade offers the opportunity to achieve the ultimate levels of selectivity through the subnanoscale design tailoring of novel catalysts. Experimental techniques capable of investigating the fundamental nature of the active sites of novel SACs in well-controlled model studies offer the chance to reveal promising insights. We report here an approach to accomplish this through the soft landing of mass-selected, ultrasmall metal oxide cluster ions, in which a single noble metal atom bound to a metal oxide moiety serves as a model SAC active site. This method allows the preparation of model catalysts in which monodispersed neutral SAC model active sites are decorated across an inert electrically conductive support at submonolayer surface coverage, in this case, Pt₁Zr₂O₇ clusters supported on highly oriented pyrolytic graphite (HOPG). The results contained herein show the characterization of the Pt₁Zr₂O₇/HOPG model catalyst by X-ray photoelectron spectroscopy (XPS), along with an investigation of its reactivity toward the functionalized hydrocarbon molecule, 1-propanamine. Through temperature-programmed desorption/reaction (TPD/R) experiments it was shown that Pt₁Zr₂O₇/HOPG decomposes 1-propanamine exclusively into propionitrile and H₂, which desorb at 425 and 550 K, respectively. Conversely, clusters without the single platinum atom, that is, Zr₂O₇/HOPG, exhibited no reactivity toward 1-propanamine. Hence, the single platinum atom in Pt₁Zr₂O₇/HOPG was found to play a critical role in the observed reactivity.

Effect of nanostructuring on the interaction of CO2 with molybdenum carbide nanoparticles

Transition metal carbides are increasingly used as catalysts for the transformation of CO₂ into useful chemicals. Recently, the effect of nanostructuring of such carbides has started to gain relevance in tailoring their catalytic capabilities. Catalytic materials based on molybdenum carbide nanoparticles (MoC_y) have shown a remarkable ability to bind CO₂ at room temperature and to hydrogenate it into oxygenates or light alkanes. However, the involved chemistry is largely unknown. In the present work, a systematic computational study is presented aiming to elucidate the chemistry behind the bonding of CO₂ with a representative set of MoC_y nanoparticles of increasing size, including stoichiometric and non-stoichiometric cases. The obtained results provide clear trends to tune the catalytic activity of these systems and to move towards more efficient CO₂ transformation processes.

Navigate Flying Molecular Elephants Safely to the Ground: Mass-Selective Soft Landing up to the Mega-Dalton Range by Electrospray Controlled Ion-Beam Deposition

The prototype of a highly versatile and efficient preparative mass spectrometry system used for the deposition of molecules in ultrahigh vacuum (UHV) is presented, along with encouraging performance data obtained using four model species that are thermolabile or not sublimable. The test panel comprises two small organic compounds, a small and very large protein, and a large DNA species covering a 4-log mass range up to 1.7 MDa as part of a broad spectrum of analyte species evaluated to date. Three designs of innovative ion guides, a novel digital mass-selective quadrupole (dQMF), and a standard electrospray ionization (ESI) source are combined to an integrated device, abbreviated electrospray controlled ion-beam deposition (ES-CIBD). Full control is achieved by (i) the square-wave-driven radiofrequency (RF) ion guides with steadily tunable frequencies, including a dQMF allowing for investigation, purification, and deposition of a virtually unlimited m/z range, (ii) the adjustable landing energy of ions down to ∼2 eV/z enabling integrity-preserving soft landing, (iii) the deposition in UHV with high ion beam intensity (up to 3 nA) limiting contaminations and deposition time, and (iv) direct coverage control via the deposited charge. The maximum resolution of R = 650 and overall efficiency up to T_total = 4.4% calculated from the solution to UHV deposition are advantageous, whereby the latter can be further enhanced by optimizing ionization performance. In the setup presented, a scanning tunneling microscope (STM) is attached for in situ UHV investigations of deposited species, demonstrating a selective, structure-preserving process and atomically clean layers.

Theoretical Perspective on Operando Spectroscopy of Fluxional Nanocatalysts

Improvements in operando spectroscopy have enabled the catalysis community to investigate the dynamic nature of catalysts under operating conditions with increasing detail. Still, the highly dynamic nature of some catalysts, such as fluxional supported subnano clusters, presents a formidable challenge even for the most state-of-the-art techniques. The reason is that such fluxional catalytic interfaces contain a variety of thermally accessible states. Operando spectroscopies used in catalysis generally fall into two categories: ensemble-based techniques, which provide spectra containing the signals of the entire ensemble of states of the catalyst and are not necessarily dominated by the most active species, and localized techniques, which provide atomistic-level information about the dynamics of active sites in a very small area, which might not include the most active species. Combining many different kinds of techniques can provide detailed insight; however, we propose that effective utilization of specific computational techniques and approaches within the fluxionality paradigm can fill the gap and enable atomistic characterization of the most relevant catalytic sites.

Size-dependent phase transitions boost catalytic activity of sub-nanometer gold clusters

The characterization and identification of the dynamics of cluster catalysis are crucial to unraveling the origin of catalytic activity. However, the dynamical catalytic effects during the reaction process remain unclear. Herein, we investigate the dynamic coupling effect of elementary reactions with the structural fluctuations of sub-nanometer Au clusters with different sizes using ab initio molecular dynamics and the free energy calculation method. It was found that the adsorption-induced solid-to-liquid phase transitions of the cluster catalysts give rise to abnormal entropy increase, facilitating the proceeding of reaction, and this phase transition catalysis exists in a range of clusters with different sizes. Moreover, clusters with different sizes show different transition temperatures, resulting in a non-trivial size effect. These results unveil the dynamic effect of catalysts and help understand cluster catalysis to design better catalysts rationally.

Addressing Dynamics at Catalytic Heterogeneous Interfaces with DFT-MD: Anomalous Temperature Distributions from Commonly Used Thermostats

Density functional theory-based molecular dynamics (DFT-MD) has been widely used for studying the chemistry of heterogeneous interfacial systems under operational conditions. We report frequently overlooked errors in thermostated or constant-temperature DFT-MD simulations applied to study (electro)catalytic chemistry. Our results demonstrate that commonly used thermostats such as Nosé-Hoover, Berendsen, and simple velocity-rescaling methods fail to provide a reliable temperature description for systems considered. Instead, nonconstant temperatures and large temperature gradients within the different parts of the system are observed. The errors are not a "feature" of any particular code but are present in several ab initio molecular dynamics implementations. This uneven temperature distribution, due to inadequate thermostatting, is well-known in the classical MD community, where it is ascribed to the failure in kinetic energy equipartition among different degrees of freedom in heterogeneous systems (Harvey et al. J. Comput. Chem. 1998, 726-740) and termed the flying ice cube effect. We provide tantamount evidence that interfacial systems are susceptible to substantial flying ice cube effects and demonstrate that the traditional Nosé-Hoover and Berendsen thermostats should be applied with care when simulating, for example, catalytic properties or structures of solvated interfaces and supported clusters. We conclude that the flying ice cube effect in these systems can be conveniently avoided using Langevin dynamics.

Nonprecious Metal Borides: Emerging Electrocatalysts for Hydrogen Production

ConspectusThe development of highly active noble-metal-free catalysts for the hydrogen evolution reaction (HER) is the focus of current fundamental research, aiming for a more efficient and economically affordable water-splitting process. While most HER catalysts are studied only at the nanoscale (small particle size and high surface area), metal borides (MBs) are mostly studied in bulk form. This offers a unique opportunity for designing highly efficient and nonprecious HER MBs electrocatalysts based on structure-activity relationships, especially because of their rich compositional and structural diversity.In this Account, we focus on the importance of boron and its substructures in achieving extraordinary HER performances and the importance of using structure-activity relationships to design next-generation MBs electrocatalysts. Studying the Mo-B system, we found that the HER activity of molybdenum borides increases with increasing boron content: from Mo₂B (no B-B bonds in the structure, least active) to α-MoB and β-MoB (zigzag boron chains, intermediate activity) and MoB₂ (planar graphene-like boron layer, most active). Density functional theory (DFT) calculations have shown that the (001) boron layer in hexagonal MoB₂ (α-MoB₂) is the most active surface and has similar HER activity behavior like the benchmark Pt(111) surface. However, puckering this flat boron layer to the chair-like configuration (phosphorene-like layer) drastically reduces its activity, thereby making the rhombohedral modification of MoB₂ (Mo₂B₄ or β-MoB₂) less active than α-MoB₂. This discovery was then further supported by studies of the Mo-W-B system. In fact, the binary WB₂, which also contains the puckered boron layer, is also less active than α-MoB₂, despite containing the more active transition metal W, which performs better in elemental form than Mo. To design a superior catalyst, the more active W was then stabilized in the hexagonal α-MoB₂ structure through the synthesis of α-Mo_0.7W_0.3B₂ ,which indeed proved to be a better HER electrocatalyst than α-MoB₂. Using the isoelectronic Cr instead of W led to the α-Cr_1-xMo_xB₂ solid solution, the HER activity of which followed unexpected canonic-like (or volcano-like) behavior that perfectly matched that of the c lattice parameter trend, thereby producing the best catalyst α-Cr_0.4Mo_0.6B₂ that outperformed Pt/C at high current density, thus underscoring the effectiveness of the structure-activity concept in designing highly active catalysts. This concept was further applied to the V-B system, leading to the discovery of an unexpected boron chain dependency of the HER activity that ultimately led to the prediction of new V_xB_y catalysts and their crystal structures and overpotentials. Finally, reducing the particle sizes of all of these bulk crystalline catalysts is also possible and offers an even greater potential as demonstrated for nanoscale a-MoB₂ and VB₂. Nevertheless, these crystalline nanomaterials are still highly agglomerated due to the high temperature required for their synthesis, thus the synthesis of highly dispersed MBs is an urgent goal that will enable the fulfillment of their extraordinary potential in the future.

Exploring the materials space in the smallest particle size range: From heterogeneous catalysis to electrocatalysis and photocatalysis

Ultrasmall clusters of subnanometer size can possess unique and even unexpected physical and chemical propensities which make them interesting in various fields of basic science and for potential applications, such...

When Fluxionality Beats Size Selection: Acceleration of Ostwald Ripening of Sub-Nano Clusters

Size selection was demonstrated to suppress Ostwald ripening of supported catalytic nanoparticles. When the supported clusters are subnanometer in size and highly fluxional, such as Pt clusters on the rutile TiO₂ (110) surface, this paradigm breaks down, and the established theory of sintering needs a revision. At temperatures characteristic of catalysis (i.e. 700 K), sub-nano clusters thermally populate many low-energy metastable isomers. As these isomers all have different geometric and electronic structures, and thus, formation and dissociation energies (in lieu of surface energy), Ostwald ripening is not suppressed, despite the size-selection. However, some clusters arise as magic numbers in terms of sintering stability at the ensemble level. Acceleration of sintering by metastable species persists though weakens in polydisperse cluster systems. We propose a competing pathways theory for sintering, which at the atomistic level describes the found size-specific sintering behavior.

Dehydrogenation of light alkanes to mono-olefins

In the past several decades, light alkane dehydrogenation to mono-olefins, especially propane dehydrogenation to propylene has gained widespread attention and much development in the field of research and commercial application. Under suitable conditions, the supported Pt-Sn and CrO_x catalysts widely used in industry exhibit satisfactory dehydrogenation activity and selectivity. However, the high cost of Pt and the potential environmental problems of CrO_x have driven researchers to improve the coking and sintering resistance of Pt catalysts, and to find new non-noble metal and environment-friendly catalysts. As for the development of the reactor, it should be noted that low operation pressure is beneficial for improving the single-pass conversion, decreasing the amount of unconverted alkane recycled back to the reactor, and reducing the energy consumption of the whole process. Therefore, the research direction of reactor improvement is towards reducing the pressure drop. This review is aimed at introducing the characteristics of the dehydrogenation reaction, the progress made in the development of catalysts and reactors, and a new understanding of reaction mechanism as well as its guiding role in the development of catalyst and reactor.

Unravelling the reactivity of metastable molybdenum carbide nanoclusters in the C-H bond activation of methane, ethane and ethylene

C-H bond activation steps in non-oxidative methane dehydroaromatization (MDA), constitute a key functionalization of the reactant and adsorbed species to form aromatics. Previous studies have focused on studying the energetics of these steps at the most stable active sites involving molybdenum carbide species. Herein, a different paradigm is presented via studying the reactivity of a metastable molybdenum carbide (Mo₂C₆) nanocluster for the C-H bond activation of methane, ethane, and ethylene and comparing it with the reactivity of the lowest energy Mo₂C₆ nanocluster. Interestingly, the metastable nanocluster is observed to result in a consistent reduction (by half) in the C-H bond activation barrier of the respective alkane and alkene molecules compared to the global minimum isomer. This specific metastable form of the nanocluster is identified from a cascade genetic algorithm search, which facilitated a rigorous scan of the potential energy surface. We attribute this significant lowering of the C-H bond activation barrier to unique co-planar orbital overlap between the reactant molecule and active centers on the metastable nanocluster. Based on geometrical and orbital analysis of the transition states arising during the C-H bond activation of methane, ethane, and ethylene, a proton-coupled electron transfer mechanism is proposed that facilitated C-H bond cleavage. Motivated by the high reactivity for C-H bond activation observed on the metastable species, a contrasting framework to analyze the elementary-step rate contributions is presented. This is based on the statistical ensemble analysis of nanocluster isomers, where the calculated rates on respective isomers are normalized with respect to the Boltzmann probability distribution. From this framework, the metastable isomer is observed to provide significant contributions to the ensemble average representations of the rate constants calculated for C-H bond activation during the MDA reaction.

Metal-based catalysts for the non-oxidative dehydrogenation of light alkanes to light olefins

This review provides an overview of metal-based catalysts, including Pt-, Pd-, Rh- and Ni-based bimetallic catalysts for non-oxidative dehydrogenation of light alkanes to olefins.

Propane dehydrogenation: catalyst development, new chemistry, and emerging technologies

This review describes recent advances in the propane dehydrogenation process in terms of emerging technologies, catalyst development and new chemistry.

Principles of Operation of a Rotating Wall Mass Analyzer for Preparative Mass Spectrometry

In this contribution, we describe the principles of operation of a rotating wall mass analyzer (RWMA), a mass-dispersive device for preparative mass spectrometry. Ions of different m/z are spatially separated by RWMA and deposited onto ring-shaped areas of distinct radii on a surface. We use a combination of an analytical equation for predicting the radius of the deposition ring and SIMION simulations to understand how to optimize the experimental conditions for the separation of multicomponent mixtures. The results of these simulations are compared with the experimental data. We introduce a universal mass calibration procedure, based on a series of polyacrylamide ions, which is subsequently used to predict the deposition radii of unknown analytes. The calibration is independent of the polarity, kinetic energy, and charge state of the ion as demonstrated by assigning m/z values of different analytes including multiply charged ubiquitin ions. We demonstrate that mass resolution of the RWMA is affected by the width and kinetic energy distribution of the ion beam. The best mass resolution obtained in this study is m/Δm = ∼20. Preparative mass spectrometry using RWMA provides the advantages of simplicity, compactness, and low fabrication cost, which are particularly promising for the development of miniaturized instrumentation. The results presented in this work can be readily adapted to preparative separation of a variety of charged species of interest to the broad scientific community.

Catalytic conversion of ethane to valuable products through non-oxidative dehydrogenation and dehydroaromatization

Chemical utilization of ethane to produce valuable chemicals has become especially attractive since the expanded utilization of shale gas in the United States and associated petroleum gas in the Middle East. Catalytic conversion to ethylene and aromatic hydrocarbons through non-oxidative dehydrogenation and dehydroaromatization of ethane (EDH and EDA) are potentially beneficial technologies because of their high selectivity to products. The former represents an attractive alternative to conventional thermal cracking of ethane. The latter can produce valuable aromatic hydrocarbons from a cheap feedstock. Nevertheless, further progress in catalytic science and technology is indispensable to implement these processes beneficially. This review summarizes progress that has been achieved with non-oxidative EDH and EDA in terms of the nature of active sites and reaction mechanisms. Briefly, platinum-, chromium- and gallium-based catalysts have been introduced mainly for EDH, including effects of carbon dioxide co-feeding. Efforts to use EDA have emphasized zinc-modified MFI zeolite catalysts. Finally, some avenues for development of catalytic science and technology for ethane conversion are summarized.

Reagent-Triggered Isomerization of Fluxional Cluster Catalyst via Dynamic Coupling

Metallic cluster catalysts have many thermodynamically accessible isomers with diverse active sites and low reaction barriers, and lately a strong hypothesis emerged that the many catalyst states collectively drive the catalysis. However, it remained a hypothesis that catalyst isomerization is actually kinetically feasible under the current reaction conditions. Using high-temperature dynamics simulations and sampling, a range of orientations, and vibrational energy distributions, we probe how thermal effects and molecular events affect cluster catalyst dynamics. We show that even such a delicate affair as the dissociation or scattering of a methane molecule on the heavy and thus slow Pt₁₃ cluster triggers substantial isomerization of the catalyst, far beyond thermal at 700 K. A kinetic coupling between the methane activity and cluster catalyst dynamics is observed. In return, the thermal dynamics of the cluster affects the methane reaction and scattering probabilities. Hence, molecular events at the surfaces of fluxional cluster catalysts should facilitate the population of an ensemble of catalyst states under the current reaction conditions, with implications for available active sites, reaction mechanisms, and apparent rates.

CO2 Hydrogenation to Formate and Formic Acid by Bimetallic Palladium-Copper Hydride Clusters

Mass spectrometric analysis of the anionic products of interaction between bimetallic palladium-copper tetrahydride anions, PdCuH₄^-, and carbon dioxide, CO₂, in a reaction cell shows an efficient generation of the PdCuCO₂H₄^- intermediate and formate/formic acid complexes. Multiple structures of PdCuH₄^- and PdCuCO₂H₄^- are identified by a synergy between anion photoelectron spectroscopy and quantum chemical calculations. The higher energy PdCuH₄^- isomer is shown to drive the catalytic hydrogenation of CO₂, emphasizing the importance of accounting for higher energy isomers for cluster catalytic activity. This study represents the first example of CO₂ hydrogenation by bimetallic hydride clusters.

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of the chemical industry. Being able to tune and design efficient catalysts for processes of interest is of the utmost importance, and for this, a molecular-level understanding of heterogeneous catalysts is the first step and indeed a prime focus of modern catalysis research. For a long time, the single most thermodynamically stable structure of the catalytic interface attained under the reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures and pressures and binding and changing reagents. Among particularly dynamic interfaces are such widely used catalysts as crystalline and amorphous surfaced supporting (sub)nanometallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states under catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many structures (rather than one structure). Every member in the ensemble contributes to the properties of the catalyst differently, in proportion to its probability of being populated. This new notion flips the established paradigm and calls for a new theory, new modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported subnanocluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, and CuPd clusters, which are catalytic toward oxidative and nonoxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous "every atom counts" phenomenon. The support and adsorbates affect catalyst structures, and the state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima and pathways sampling as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and prediction of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

An Eight-Atom Iridium-Aluminum Oxide Cluster IrAlO6+ Catalytically Oxidizes Six CO Molecules

Fundamental understanding regarding oxygen storage capacity involving how and why an active site can buffer a large number of oxygen atoms in redox processes is vital to the design of advanced oxygen storage materials, while it is challenging because of the complexity of heterogeneous catalysis. Herein, we identified that an eight-atom iridium-aluminum oxide cluster IrAlO₆⁺ can transfer all the oxygen atoms to catalytically oxidize six CO molecules. This finding represents a breakthrough in cluster catalysis where at most three oxygen atoms from a heteronuclear metal oxide cluster can be catalytically involved in CO oxidation. We found that oxygen prefers to be stored on aluminum to form an O₃^-• radical in the energetically unfavorable IrAlO₆⁺ isomer and generate the low-coordinated iridium that is pivotal to capturing CO and triggering the catalysis. The powerful electron cycling capability of iridium and the cooperative iridium-aluminum interplay are emphasized to drive the oxygen atom-transfer behavior.

Dynamic charge and oxidation state of Pt/CeO2 single-atom catalysts

The catalytic activity of metals supported on oxides depends on their charge and oxidation state. Yet, the determination of the degree of charge transfer at the interface remains elusive. Here, by combining density functional theory and first-principles molecular dynamics on Pt single atoms deposited on the CeO₂ (100) surface, we show that the common representation of a static metal charge is oversimplified. Instead, we identify several well-defined charge states that are dynamically interconnected and thus coexist. The origin of this new class of strong metal-support interactions is the relative position of the Ce(4f) levels with respect to those of the noble metal, allowing electron injection to (or recovery from) the support. This process is phonon-assisted, as the Ce(4f) levels adjust by surface atom displacement, and appears for other metals (Ni) and supports (TiO₂). Our dynamic model explains the unique reactivity found for activated single Pt atoms on ceria able to perform CO oxidation, meeting the Department of Energy 150 °C challenge for emissions.