Visualizing the Promoting Role of Interfacial Water in the Deprotonation of Formic Acid on Cu(111)

Water plays a crucial role in various heterogeneous catalytic reactions, but the atomic-scale characterization of how water participates in these chemical processes remains a significant challenge. Here we directly visualize the promoting role of interfacial water in the deprotonation of formic acid (FA) on a metal surface, using combined scanning tunneling microscopy and qPlus-based noncontact atomic force microscopy. We find the dissociation of FA when coadsorbed with water on the Cu(111) surface, resulting in the formation of hydronium and formate ions. Interestingly, most of the hydrated proton and formate ions exhibit a phase-separated behavior on Cu(111), in which Eigen and Zundel cations assemble into a monolayer hexagonal hydrogen-bonding (H-bonding) network, and bidentate formate ions are solvated with water and aggregate into one-dimensional chains or two-dimensional H-bonding networks. This phase-separated behavior is essential for preventing the proton transfer back from hydronium to formate and the reformation of FA. Density functional theory calculations reveal that the participation of water significantly reduces the deprotonation barrier of FA on Cu(111), in which water catalyzes the decomposition of FA through the Grotthuss proton transfer mechanism. In addition, the separate solvation of hydronium and bidentate formate ions is energetically preferred due to the enhanced interaction with the copper substrate. The promoting role of water in the deprotonation of FA is further confirmed by the temperature-programmed desorption experiment, which shows that the intensity of the H2 desorption peak significantly increases and the desorption of FA declines when water and FA coadsorbed on the Cu(111) surface.

Dehydrogenation and dehydration of formic acid over orthorhombic molybdenum carbide

•

Formic acid (HCOOH) adsorption on β-Mo₂C is exothermic and favours a configuration parallel to the surface.

•

Once adsorbed, thermodynamics favour cleavage of the H—COOH bond to form CO.

•

CO bonds strongly to the surface, potentially poisoning the catalyst.

•

Therefore, kinetics favour dehydrogenation mechanism with CO₂ continuously formed.

The dehydrogenation and dehydration of formic acid is investigated on the β-Mo₂C (100) catalyst surface using time independent density functional theory. The energetics of the two mechanisms are calculated, and the thermochemistry and kinetics are discussed using the transition state theory. Subsequently, microkinetic modelling of the system is conducted, considering the batch reactor model. The potential energy landscape of the reaction shows a thermodynamically favourable cleavage of H—COOH to form CO; however, the kinetics show that the dehydrogenation mechanism is faster and CO₂ is continuously formed. The effect of HCOOH adsorption on the surface is also analysed, in a temperature-programmed desorption, with the conversion proceeding at under 350 K and desorption of CO₂ is observed with a selectivity of about 100 %, in line with the experimental reports.

Controlled synthesis of Bi- and tri-nuclear Cu-oxo nanoclusters on metal-organic frameworks and the structure-reactivity correlations

Precisely tuning the nuclearity of supported metal nanoclusters is pivotal for designing more superior catalytic systems, but it remains practically challenging. By utilising the chemical and molecular specificity of UiO-66-NH2 (a Zr-based metal-organic framework), we report the controlled synthesis of supported bi- and trinuclear Cu-oxo nanoclusters on the Zr6O4 nodal centres of UiO-66-NH2. We revealed the interplay between the surface structures of the active sites, adsorption configurations, catalytic reactivities and associated reaction energetics of structurally related Cu-based 'single atoms' and bi- and trinuclear species over our model photocatalytic formic acid reforming reaction. This work will offer practical insight that fills the critical knowledge gap in the design and engineering of new-generation atomic and nanocluster catalysts. The precise control of the structure and surface sensitivities is important as it can effectively lead to more reactive and selective catalytic systems. The supported bi- and trinuclear Cu-oxo nanoclusters exhibit notably different catalytic properties compared with the mononuclear 'Cu1' analogue, which provides critical insight for the engineering of more superior catalytic systems.

Strong Metal-Support Interaction for 2D Materials: Application in Noble Metal/TiB2 Heterointerfaces and their Enhanced Catalytic Performance for Formic Acid Dehydrogenation

Strong metal-support interaction (SMSI) is a phenomenon commonly observed on heterogeneous catalysts. Here, direct evidence of SMSI between noble metal and 2D TiB₂ supports is reported. The temperature-induced TiB₂ overlayers encapsulate the metal nanoparticles, resulting in core-shell nanostructures that are sintering-resistant with metal loadings as high as 12.0 wt%. The TiO_x -terminated TiB₂ surfaces are the active sites catalyzing the dehydrogenation of formic acid at room temperature. In contrast to the trade-off between stability and activity in conventional SMSI, TiB₂ -based SMSI promotes catalytic activity and stability simultaneously. By optimizing the thickness and coverage of the overlayer, the Pt/TiB₂ catalyst displays an outstanding hydrogen productivity of 13.8 mmol g^-1 _cat h^-1 in 10.0 m aqueous solution without any additive or pH adjustment, with >99.9% selectivity toward CO₂ and H₂ . Theoretical studies suggest that the TiB₂ overlayers are stabilized on different transition metals through an interplay between covalent and electrostatic interactions. Furthermore, the computationally determined trends in metal-TiB₂ interactions are fully consistent with the experimental observations regarding the extent of SMSI on different transition metals. The present research introduces a new means to create thermally stable and catalytically active metal/support interfaces for scalable chemical and energy applications.

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate-determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted-Evans-Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves are discussed. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. We explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. We finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

Combining Computational Modeling with Reaction Kinetics Experiments for Elucidating the In Situ Nature of the Active Site in Catalysis

Microkinetic modeling based on density functional theory (DFT) derived energetics is important for addressing fundamental questions in catalysis. The quantitative fidelity of microkinetic models (MKMs), however, is often insufficient to conclusively infer the mechanistic details of a specific catalytic system. This can be attributed to a number of factors such as an incorrect model of the active site for which DFT calculations are performed, deficiencies in the hypothesized reaction mechanism, inadequate consideration of the surface environment under reaction conditions, and intrinsic errors in the DFT exchange-correlation functional. Despite these limitations, we aim at developing a rigorous understanding of the reaction mechanism and of the nature of the active site for heterogeneous catalytic chemistries under reaction conditions. By achieving parity between experimental and modeling outcomes through robust parameter estimation and by ensuring coverage-consistency between DFT calculations and MKM predictions, it is possible to systematically refine the mechanistic model and, thereby, our understanding of the catalytic active site in situ.Our general approach consists of developing ab initio informed MKM for a given active site and then re-estimating the energies of the transition and intermediate states so that the model predictions match quantities measured in reaction kinetics experiments. If (i) model-experiment parity is high, (ii) the adjustments to the DFT-derived energetics for a given model of the active site are rationalized within the errors of standard DFT exchange-correlation functionals, and (iii) the resultant MKM predicts surface coverages that are consistent with those assumed in the DFT calculations used to initialize the MKM, we conclude that we have correctly identified the active site and the reaction mechanism. If one or more of these requirements are not met, we iteratively refine our model by updating our hypothesis for the structure of the active site and/or by incorporating coverage effects, until we obtain a high-fidelity coverage-self-consistent MKM whose final kinetic and thermodynamic parameters are within error of the values derived from DFT.Using the catalytic reaction of formic acid (FA, HCOOH) decomposition over transition-metal catalysts as an example, here we provide an account of how we applied this algorithm to study this chemistry on powder Au/SiC and Pt/C catalysts. For the case of Au catalysts, on which the FA decomposition occurred exclusively through the dehydrogenation reaction (HCOOH → CO₂+H₂), our approach was used to iteratively refine the model starting from the (111) facet until we found that specific ensembles of Au atoms present in sub-nanometer clusters can describe the active site for this catalysis. For the case of Pt catalysts, wherein both dehydrogenation (HCOOH → CO₂ + H₂) and dehydration (HCOOH → CO + H₂O) reactions were active, our approach identified that a partially CO*-covered (111) surface serves as the active site and that CO*-assisted steps contributed substantially to the overall FA decomposition activity. Finally, we suggest that once the active site and the mechanism are conclusively identified, the model can subsequently serve as a high-quality basis for designing specific goal-oriented experiments and improved catalysts.

The Role of Support in Formic Acid Decomposition on Gold Catalysts

Formic acid (FA) can easily be decomposed, affording molecular hydrogen through a controllable catalytic process, thus attaining great importance as a convenient hydrogen carrier for hydrogen energetics. Supported gold nanoparticles are considered to be among the most promising catalysts for such applications. However, questions remain regarding the influence of the catalyst support on the reaction selectivity. In this study, we have examined the catalytic activity of typical gold catalysts, such as Au/TiO2, Au/SiO2, and Au/Al2O3 in decomposition of FA, and then compared it with the catalytic activity of corresponding supports. The performance of each catalyst and support was evaluated using a gas-flow packed-bed reactor. It is shown that the target reaction, FA → H2 + CO2, is provided by the presence of gold nanoparticles, whereas the concurrent, undesirable pathway, such as FA → H2O + CO, results exclusively from the acid-base behavior of supports.

Hydrogen Production from Formic Acid Attained by Bimetallic Heterogeneous PdAg Catalytic Systems

The production of H2 from the so-called Liquid Organic Hydrogen Carriers (LOHC) has recently received great focus as an auspicious option to conventional hydrogen storage technologies. Among them, formic acid, the simplest carboxylic acid, has recently emerged as one of the most promising candidates. Catalysts based on Pd nanoparticles are the most fruitfully investigated, and, more specifically, excellent results have been achieved with bimetallic PdAg-based catalytic systems. The enhancement displayed by PdAg catalysts as compared to the monometallic counterpart is ascribed to several effects, such as the formation of electron-rich Pd species or the increased resistance against CO-poisoning. Aside from the features of the metal active phases, the properties of the selected support also play an important role in determining the final catalytic performance. Among them, the use of carbon materials has resulted in great interest by virtue of their outstanding properties and versatility. In the present review, some of the most representative investigations dealing with the design of high-performance PdAg bimetallic heterogeneous catalysts are summarised, paying attention to the impact of the features of the support in the final ability of the catalysts towards the production of H2 from formic acid.

Catalytic activity of palladium-doped silver dilute nanoalloys for formate oxidation from a theoretical perspective

The large-scale practical application of formate oxidation reaction (FOR) catalysts is hindered by their low activity and high cost. Herein, for the first time, a series of Pd-doped Ag dilute nanoalloys is demonstrated to have high catalytic activity in FOR with reduced consumption of Pd metals through density functional theory calculations, where the effects of potential, solvent and spin on catalytic performance are discussed. The Pd1Ag(111) single-atom alloy (SAA) exhibits higher FOR catalytic activity as reflected by the low limiting potential of 0.026 eV for the direct association path and a value of 0.084 eV for the direct dissociation path, and the lowest activation energy of 0.774 eV for the rate-determining-step in the direct dissociation path compared with Pd2Ag(111) and Pd3Ag(111) dilute alloys. Pd1Ag(111) SAA exhibits an extremely narrow sharp peak in the partial density of states from -0.75 to -2.0 eV, which is due to the free-atom-like electronic structure of the single Pd atom. The isolated Pd single atom is more stable by -0.041 and -0.097 eV, respectively, than the aggregated Pd2 and Pd3 atom clusters on the Ag(111) surface, which verifies the potential application of Pd1Ag(111) SAA in experiments. Overall, this work further elucidates the theoretical profile of FOR and provides a new strategy for designing the catalytic reaction at the atomic level.

Convolutional Neural Network of Atomic Surface Structures To Predict Binding Energies for High-Throughput Screening of Catalysts

High-throughput screening of catalysts can be performed using density functional theory calculations to predict catalytic properties, often correlated with adsorbate binding energies. However, more complete investigations would require an order of 2 more calculations compared to the current approach, making the computational cost a bottleneck. Recently developed machine-learning methods have been demonstrated to predict these properties from hand-crafted features but have struggled to scale to large composition spaces or complex active sites. Here, we present an application of a deep-learning convolutional neural network of atomic surface structures using atomic and Voronoi polyhedra-based neighbor information. The model effectively learns the most important surface features to predict binding energies. Our method predicts CO and H binding energies after training with 12 000 data for each adsorbate with a mean absolute error of 0.15 eV for a diverse chemical space. Our method is also capable of creating saliency maps that determine atomic contributions to binding energies.

Hydrogen Production from Formic Acid over Au Catalysts Supported on Carbon: Comparison with Au Catalysts Supported on SiO2 and Al2O3

Characteristics and catalytic activity in hydrogen production from formic acid of Au catalysts supported on porous N-free (Au/C) and N-doped carbon (Au/N-C) have been compared with those of Au/SiO2 and Au/Al2O3 catalysts. Among the catalysts examined, the Au/N-C catalyst showed the highest Au mass-based catalytic activity. The following trend was found at 448 K: Au/N-C > Au/SiO2 > Au/Al2O3, Au/C. The trend for the selectivity in hydrogen production was different: Au/C (99.5%) > Au/Al2O3 (98.0%) > Au/N-C (96.3%) > Au/SiO2 (83.0%). According to XPS data the Au was present in metallic state in all catalysts after the reaction. TEM analysis revealed that the use of the N-C support allowed obtaining highly dispersed Au nanoparticles with a mean size of about 2 nm, which was close to those for the Au catalysts on the oxide supports. However, it was by a factor of 5 smaller than that for the Au/C catalyst. The difference in dispersion could explain the difference in the catalytic activity for the carbon-based catalysts. Additionally, the high activity of the Au/N-C catalyst could be related to the presence of pyridinic type nitrogen on the N-doped carbon surface, which activates the formic acid molecule forming pyridinium formate species further interacting with Au. This was confirmed by density functional theory (DFT) calculations. The results of this study may assist the development of novel Au catalysts for different catalytic reactions.

Efficient dehydrogenation of a formic acid-ammonium formate mixture over Au3Pd1 catalyst

A series of AuPd/C catalysts were prepared and tested for the first time for active and stable dehydrogenation of a formic acid-ammonium formate (FA-AF) mixture. The catalysts with different Au-to-Pd molar ratios were prepared using a facile simultaneous reduction method and characterized using transmission electron microscopy (TEM), high-resolution TEM, energy dispersive X-ray spectroscopy, X-ray diffraction, and X-ray photoelectron spectroscopy. It was found that the catalytic activity and stability of the Au3Pd1/C catalyst was the best. The initial turnover frequency for the dehydrogenation of the FA-AF mixture over the Au3Pd1/C catalyst can reach 407.5 h-1 at 365 K. The reaction order with respect to FA and AF is 0.25 and 0.55, respectively. The apparent activation energy of dehydrogenation is 23.3 ± 1.3 kJ mol-1. The catalytic activity of the Au3Pd1/C catalyst remains ca. 88.0% after 4 runs, which is much better than the single Pd/C catalyst. The mechanism for the dehydrogenation is also discussed.

On the nature of active sites for formic acid decomposition on gold catalysts

Atomic scale size-sensitivity of the catalytic properties of sub-nanometer gold clusters for HCOOH decomposition.

Connection between macroscopic kinetic measurables and the degree of rate control

Macroscopic kinetic measurables are linked to elementary reaction steps by the degree of rate control.

Ab initio coverage-dependent microkinetic modeling of benzene hydrogenation on Pd(111)

Coverage-dependent calculations are required for an accurate DFT-based prediction of the activity and a correct mechanistic understanding of catalytic hydrogenation.

Active sites and mechanisms for H₂O₂ decomposition over Pd catalysts

A combination of periodic, self-consistent density functional theory (DFT-GGA-PW91) calculations, reaction kinetics experiments on a SiO2-supported Pd catalyst, and mean-field microkinetic modeling are used to probe key aspects of H2O2 decomposition on Pd in the absence of cofeeding H2 We conclude that both Pd(111) and OH-partially covered Pd(100) surfaces represent the nature of the active site for H2O2 decomposition on the supported Pd catalyst reasonably well. Furthermore, all reaction flux in the closed catalytic cycle is predicted to flow through an O-O bond scission step in either H2O2 or OOH, followed by rapid H-transfer steps to produce the H2O and O2 products. The barrier for O-O bond scission is sensitive to Pd surface structure and is concluded to be the central parameter governing H2O2 decomposition activity.

Support effect for nanosized Au catalysts in hydrogen production from formic acid decomposition

Au catalysts with the same particle sizes demonstrate the following order of activity in formic acid decomposition: Au/Al₂O₃ > Au/ZrO₂ ∼ Au/CeO₂ > Au/La₂O₃ > Au/MgO.