Identification of reaction intermediates in the decomposition of formic acid on Pd.

Uncovering the role of reaction intermediates is crucial to developing an understanding of heterogeneous catalysis because catalytic reactions often involve complex networks of elementary steps. Identifying the reaction intermediates is often difficult because their short lifetimes and low concentrations make it difficult to observe them with surface sensitive spectroscopic techniques. In this paper we report a different approach to identify intermediates for the formic acid decomposition reaction on Pd(111) and Pd(332) based on accurate measurements of isotopologue specific thermal reaction rates. At low surface temperatures (∼400 K) CO2 formation is the major reaction pathway. The CO2 kinetic data show this occurs via two temporally resolved reaction processes. Thus, there must be two parallel pathways which we attribute to the participation of two intermediate species in the reaction. Isotopic substitution reveals large and isotopologue specific kinetic isotope effects that allow us to identify the two key intermediates as bidentate formate and carboxyl. The decomposition of the bidentate formate is substantially slower than that of carboxyl. On Pd(332), at high surface temperatures (643 K to 693 K) we observe both CO and CO2 production. The observation of CO formation reinforces the conclusion of calculations that suggest the carboxyl intermediate plays a major role in the water-gas shift reaction, where carboxyl exhibits temperature dependent branching between CO2 and CO.

Velocity-resolved kinetics of site-specific carbon monoxide oxidation on platinum surfaces.

Catalysts are widely used to increase reaction rates. They function by stabilizing the transition state of the reaction at their active site, where the atomic arrangement ensures favourable interactions 1 . However, mechanistic understanding is often limited when catalysts possess multiple active sites-such as sites associated with either the step edges or the close-packed terraces of inorganic nanoparticles2-4-with distinct activities that cannot be measured simultaneously. An example is the oxidation of carbon monoxide over platinum surfaces, one of the oldest and best studied heterogeneous reactions. In 1824, this reaction was recognized to be crucial for the function of the Davy safety lamp, and today it is used to optimize combustion, hydrogen production and fuel-cell operation5,6. The carbon dioxide products are formed in a bimodal kinetic energy distribution7-13; however, despite extensive study 5 , it remains unclear whether this reflects the involvement of more than one reaction mechanism occurring at multiple active sites12,13. Here we show that the reaction rates at different active sites can be measured simultaneously, using molecular beams to controllably introduce reactants and slice ion imaging14,15 to map the velocity vectors of the product molecules, which reflect the symmetry and the orientation of the active site 16 . We use this velocity-resolved kinetics approach to map the oxidation rates of carbon monoxide at step edges and terrace sites on platinum surfaces, and find that the reaction proceeds through two distinct channels11-13: it is dominated at low temperatures by the more active step sites, and at high temperatures by the more abundant terrace sites. We expect our approach to be applicable to a wide range of heterogeneous reactions and to provide improved mechanistic understanding of the contribution of different active sites, which should be useful in the design of improved catalysts.

Highly Rotationally Excited N2 Reveals Transition-State Character in the Thermal Decomposition of N2O on Pd(110).

We employ time-slice and velocity map ion imaging methods to explore the quantum-state resolved dynamics in thermal N2O decomposition on Pd(110). We observe two reaction channels: a thermal channel that is ascribed to N2 products initially trapped at surface defects and a hyperthermal channel involving a direct release of N2 to the gas phase from N2O adsorbed on bridge sites oriented along the [001] azimuth. The hyperthermal N2 is highly rotationally excited up to J = 52 (v″ = 0) with a large average translational energy of 0.62 eV. Between 35 and 79% of the estimated barrier energy (1.5 eV) released upon dissociation of the transition state (TS) is taken up by the desorbed hyperthermal N2. The observed attributes of the hyperthermal channel are interpreted by post-transition-state classical trajectories on a density functional theory-based high-dimensional potential energy surface. The energy disposal pattern is rationalized by the sudden vector projection model, which attributes to unique features of the TS. Applying detailed balance, we predict that in the reverse Eley-Rideal reaction, both N2 translational and rotational excitation promote N2O formation.

Binding Energy and Diffusion Barrier of Formic Acid on Pd(111).

Velocity-resolved kinetics is used to measure the thermal rate of formic acid desorption from Pd(111) between 228 and 273 K for four isotopologues: HCOOH, HCOOD, DCOOH, DCOOD. Upon molecular adsorption, formic acid undergoes decomposition to CO2 and H2 and thermal desorption. To disentangle the contributions of individual processes, we implement a mass-balance-based calibration procedure from which the branching ratio between desorption and decomposition for formic acid is determined. From experimentally derived elementary desorption rate constants, we obtain the binding energy 639 ± 8 meV and the diffusion barrier 370 ± 130 meV using the detailed balance rate model (DBRM). The DBRM explains the observed kinetic isotope effects.

Quantum state-resolved methane scattering from Ni(111) and NiO(111) by bolometer infrared laser tagging: The effect of surface oxidation.

We describe a novel ultrahigh vacuum state-to-state molecule/surface scattering apparatus with quantum state preparation of the incident molecular beam and angle-resolved quantum state detection of the scattered molecules. State-resolved detection is accomplished using a tunable mid-infrared laser source combined with a cryogenic bolometer detector and is applicable to any molecule with an infrared-active vibrational transition. Results on rotationally inelastic scattering of CH4 methane from a Ni(111) surface and NiO(111)/Ni(111) oxide film, obtained by the new apparatus, are presented. Molecules scattering from the oxidized surface, compared to those scattering from the bare nickel surface, are more highly excited rotationally and scatter into a broader distribution of angles. The internal alignment of molecular rotation is in addition found to be stronger in molecules scattering from the bare surface. Furthermore, the maxima of the state-resolved angular distributions shift toward and away from surface normal with increasing rotational quantum number J for the oxidized and bare surface, respectively. Finally, the rotational state populations produced in scattering from the oxidized surface are well-described by a Boltzmann distribution, while those produced in scattering from the bare surface exhibit large deviations from their best-fit Boltzmann distributions. These results point toward a marked enhancement in molecule-surface collisional energy exchange induced by oxidation of the nickel surface.

Steric Hindrance of NH3 Diffusion on Pt(111) by Co-Adsorbed O-Atoms.

A detailed velocity-resolved kinetics study of NH3 thermal desorption rates from p(2 × 2) O/Pt(111) is presented. We find a large reduction in the NH3 desorption rate due to adsorption of O-atoms on Pt(111). A physical model describing the interactions between adsorbed NH3 and O-atoms explains these observations. By fitting the model to the derived desorption rate constants, we find an NH3 stabilization on p(2 × 2) O/Pt(111) of 0.147-0.014+0.023 eV compared to Pt(111) and a rotational barrier of 0.084-0.022+0.049 eV, which is not present on Pt(111). The model also quantitatively predicts the steric hindrance of NH3 diffusion on Pt(111) due to co-adsorbed O-atoms. The derived diffusion barrier of NH3 on p(2 × 2) O/Pt(111) is 1.10-0.13+0.22 eV, which is 0.39-0.14+0.22 eV higher than that on pristine Pt(111). We find that Perdew Burke Ernzerhof (PBE) and revised Perdew Burke Ernzerhof (RPBE) exchange-correlation functionals are unable to reproduce the experimentally observed NH3-O adsorbate-adsorbate interactions and NH3 binding energies at Pt(111) and p(2 × 2) O/Pt(111), which indicates the importance of dispersion interactions for both systems.

Kinetics of NH3 Desorption and Diffusion on Pt: Implications for the Ostwald Process.

We report accurate time-resolved measurements of NH3 desorption from Pt(111) and Pt(332) and use these results to determine elementary rate constants for desorption from steps, from (111) terrace sites and for diffusion on (111) terraces. Modeling the extracted rate constants with transition state theory, we find that conventional models for partition functions, which rely on uncoupled degrees of freedom (DOFs), are not able to reproduce the experimental observations. The results can be reproduced using a more sophisticated partition function, which couples DOFs that are most sensitive to NH3 translation parallel to the surface; this approach yields accurate values for the NH3 binding energy to Pt(111) (1.13 ± 0.02 eV) and the diffusion barrier (0.71 ± 0.04 eV). In addition, we determine NH3's binding energy preference for steps over terraces on Pt (0.23 ± 0.03 eV). The ratio of the diffusion barrier to desorption energy is ∼0.65, in violation of the so-called 12% rule. Using our derived diffusion/desorption rates, we explain why established rate models of the Ostwald process incorrectly predict low selectivity and yields of NO under typical reactor operating conditions. Our results suggest that mean-field kinetics models have limited applicability for modeling the Ostwald process.

The Barrier for CO2 Functionalization to Formate on Hydrogenated Pt.

Understanding heterogeneous catalysis is based on knowing the energetic stability of adsorbed reactants, intermediates, and products as well as the energetic barriers separating them. We report an experimental determination of the barrier to CO2 functionalization to form bidentate formate on a hydrogenated Pt surface and the corresponding reaction energy. This determination was possible using velocity resolved kinetics, which simultaneously provides information about both the dynamics and rates of surface chemical reactions. In these experiments, a pulse of isotopically labeled formic acid (DCOOH) doses the Pt surface rapidly forming bidentate formate (DCO*O*). We then record the (much slower) rate of decomposition of DCO*O* to form adsorbed D* and gas phase CO2. We establish the reaction mechanism by dosing with O2 to form adsorbed O*, which efficiently converts H* or D* to gas phase water. H2O is formed immediately reflecting rapid loss of the acidic proton associated with formation of formate, while D2O formation proceeds more slowly and on the same time scale as the CO2 production. The temperature dependence of the reaction rate yields an activation energy that reflects the energy of the transition state with respect to DCO*O*. The derived heat of formation for DCO*O* on Pt(111) agrees well with results of microcalorimetry. The maximum release of translational energy of the formed CO2 provides a measure of the energy of the transition state with respect to the products and the barrier to the reverse process, functionalization of CO2. The comparison between the results on Pt(111) and Pt(332) shows that the barrier for CO2 functionalization is reduced by the presence of steps. The approach taken here could provide a method to optimize catalysts for CO2 functionalization.

First principles rates for surface chemistry employing exact transition state theory: application to recombinative desorption of hydrogen from Cu(111).

We present first principles calculations of the reactive flux for thermal recombinative desorption of hydrogen from Cu(111). We follow a theoretical paradigm used successfully for gas phase reactions, where electronic structure theory (DFT-GGA) is combined with transition state theory (TST). Classical ab initio molecular dynamics trajectories initiated from a thermal distribution near the transition state provide dynamical corrections to the desorption rate. We use this to calculate and study the recrossing error of TST and to directly simulate thermal desorption experiments based on a high temperature permeation method. Transition state recrossing is strongly temperature dependent and is even important in a frozen Cu-atom model. It is not influenced by inclusion of electron-hole pair excitation at the level of the local density electronic friction approximation. We also present the kinetic energy resolved flux of desorbing H2 at elevated temperature. This provides a more direct way to compare first principles theory to experiment, with no need to invoke detailed balance.

Associative desorption of hydrogen isotopologues from copper surfaces: Characterization of two reaction mechanisms.

We report quantum-state resolved measurements of angular and velocity distributions of the associative desorption of H2, HD, and D2 from Cu(111) and Cu(211) surfaces. The desorbing molecules have bimodal velocity distributions comprising a "fast" channel and a "slow" channel on both facets. The "fast channel" is promoted by both hydrogen incidence translational and vibrational energy, while the "slow channel" is promoted by vibrational energy but inhibited by translational energy. Using detailed balance, we determine state-specific reaction probabilities for dissociative adsorption and compare these to theoretical calculations. The results for the activation barrier for the "fast channel" on Cu(111) are in agreement with theory within "chemical accuracy" (1 kcal/mole). Results on the Cu(211) facet provide direct information on the effect of increasing step density, which is commonly believed to increase reactivity. Differences in reactivity on the (111) and (211) facets are subtle - quantum state specific reactivity on the (211) surface is characterized by a broader distribution of barrier heights whose average values are higher than for reaction on (111). We fully characterize the "slow channel," which has not been found in theoretical calculations although it makes up a large fraction of the reactivity in these experiments.

The dynamics of molecular interactions and chemical reactions at metal surfaces: testing the foundations of theory.

We review studies of molecular interactions and chemical reactions at metal surfaces, emphasizing progress toward a predictive theory of surface chemistry and catalysis. For chemistry at metal surfaces, a small number of central approximations are typically made: (a) the Born-Oppenheimer approximation of electronic adiabaticity, (b) the use of density functional theory at the generalized gradient approximation level, (c) the classical approximation for nuclear motion, and (d) various reduced-dimensionality approximations. Together, these approximations constitute a provisional model for surface chemical reactivity. We review work on some carefully studied examples of molecules interacting at metal surfaces that probe the validity of various aspects of the provisional model.

Observation of orientation-dependent electron transfer in molecule-surface collisions.

Molecules typically must point in specific relative directions to participate efficiently in energy transfer and reactions. For example, Förster energy transfer favors specific relative directions of each molecule's transition dipole [Förster T (1948) Ann Phys 2(1-2):55-75] and electron transfer between gas-phase molecules often depends on the relative orientation of orbitals [Brooks PR, et al. (2007) J Am Chem Soc 129(50):15572-15580]. Surface chemical reactions can be many orders of magnitude faster than their gas-phase analogs, a fact that underscores the importance of surfaces for catalysis. One reason surface reactions can be so fast is the labile change of oxidation state that commonly takes place upon adsorption, a process involving electron transfer between a solid metal and an approaching molecule. By transferring electrons to or from the adsorbate, the process of bond weakening and/or cleavage is initiated, chemically activating the reactant [Yoon B, et al. (2005) Science 307(5708):403-407]. Here, we show that the vibrational relaxation of NO--an example of electronically nonadiabatic energy transfer that is driven by an electron transfer event [Gadzuk JW (1983) J Chem Phys 79(12):6341-6348]--is dramatically enhanced when the molecule approaches an Au(111) surface with the N atom oriented toward the surface. This represents a rare opportunity to investigate the steric influences on an electron transfer reaction happening at a surface.

CO desorption from a catalytic surface: elucidation of the role of steps by velocity-selected residence time measurements.

Directly measuring the rate of a surface chemical reaction remains a challenging problem. For example, even after more than 30 years of study, there is still no agreement on the kinetic parameters for one of the simplest surface reactions: desorption of CO from Pt(111). We present a new experimental technique for determining rates of surface reactions, the velocity-selected residence time method, and demonstrate it for thermal desorption of CO from Pt(111). We use UV−UV double resonance spectroscopy to record surface residence times at selected final velocities of the desorbing CO subsequent to dosing with a pulsed molecular beam. Velocity selection differentiates trapping-desorption from direct scattering and removes influences on the temporal profile arising from the velocity distribution of the desorbing CO. The kinetic data thus obtained are of such high quality that bi-exponential desorption kinetics of CO from Pt(111) can be clearly seen. We assign the faster of the two rate processes to desorption from (111) terraces, and the slower rate process to sequential diffusion from steps to terraces followed by desorption. The influence of steps, whose density may vary from crystal to crystal, accounts for the diversity of previously reported (single exponential) kinetics results. Using transition-state theory, we derive the binding energy of CO to Pt(111) terraces, D(0)(terr) (Pt−CO) = 34 ± 1 kcal/mol (1.47 ± 0.04 eV) for the low coverage limit (≤0.03 ML) where adsorbate−adsorbate interactions are negligible. This provides a useful benchmark for electronic structure theory of adsorbates on metal surfaces.

Using Ion Imaging to Measure Velocity Distributions in Surface Scattering Experiments.

We present a new implementation of ion imaging for the study of surface scattering processes. The technique uses a combination of spatial ion imaging with laser slicing and delayed pulsed extraction. The scattering velocities of interest are parallel to the imaging plane, allowing speed and angular distributions to be extracted from a single image. The first results of direct scattering of N2 from a clean, single-crystal Au(111) surface are reported, and the speed resolution is shown to be competitive with current state-of-the-art time-of-flight methods for velocity measurements while providing simultaneous measurements of in-plane angular distributions.

An accurate full-dimensional potential energy surface for H-Au(111): Importance of nonadiabatic electronic excitation in energy transfer and adsorption.

We have constructed a potential energy surface (PES) for H-atoms interacting with fcc Au(111) based on fitting the analytic form of the energy from Effective Medium Theory (EMT) to ab initio energy values calculated with density functional theory. The fit used input from configurations of the H-Au system with Au atoms at their lattice positions as well as configurations with the Au atoms displaced from their lattice positions. It reproduces the energy, in full dimension, not only for the configurations used as input but also for a large number of additional configurations derived from ab initio molecular dynamics (AIMD) trajectories at finite temperature. Adiabatic molecular dynamics simulations on this PES reproduce the energy loss behavior of AIMD. EMT also provides expressions for the embedding electron density, which enabled us to develop a self-consistent approach to simulate nonadiabatic electron-hole pair excitation and their effect on the motion of the incident H-atoms. For H atoms with an energy of 2.7 eV colliding with Au, electron-hole pair excitation is by far the most important energy loss pathway, giving an average energy loss ≈3 times that of the adiabatic case. This increased energy loss enhances the probability of the H-atom remaining on or in the Au slab by a factor of 2. The most likely outcome for H-atoms that are not scattered also depends prodigiously on the energy transfer mechanism; for the nonadiabatic case, more than 50% of the H-atoms which do not scatter are adsorbed on the surface, while for the adiabatic case more than 50% pass entirely through the 4 layer simulation slab.

Incidence energy dependent state-to-state time-of-flight measurements of NO(v = 3) collisions with Au(111): the fate of incidence vibrational and translational energy.

We report measurements of translational energy distributions when scattering NO(vi = 3, Ji = 1.5) from a Au(111) surface into vibrational states vf = 1, 2, 3 and rotational states up to Jf = 32.5 for various incidence energies ranging from 0.11 eV to 0.98 eV. We observed that the vibration-to-translation as well as the translation-to-rotation coupling depend on translational incidence energy, EI. The vibration-to-translation coupling, i.e. the additional recoil energy observed for vibrationally inelastic (v = 3 → 2, 1) scattering, is seen to increase with increasing EI. The final translational energy decreases approximately linearly with increasing rotational excitation. At incidence energies EI > 0.5 eV, the slopes of these dependencies are constant and identical for the three vibrational channels. At lower incidence energies, the slopes gradually approach zero for the vibrationally elastic channel while they exhibit more abrupt transitions for the vibrationally inelastic channels. We discuss possible mechanisms for both effects within the context of nonadiabatic electron-hole pair mediated energy transfer and orientation effects.

Dissociation and recombination of D2 on Cu(111): ab initio molecular dynamics calculations and improved analysis of desorption experiments.

Obtaining quantitative agreement between theory and experiment for dissociative adsorption of hydrogen on and associative desorption of hydrogen from Cu(111) remains challenging. Particularly troubling is the fact that theory gives values for the high energy limit to the dissociative adsorption probability that is as much as two times larger than experiment. In the present work we approach this discrepancy in three ways. First, we carry out a new analysis of the raw experimental data for D2 associatively desorbing from Cu(111). We also perform new ab initio molecular dynamics (AIMD) calculations that include effects of surface atom motion. Finally, we simulate time-of-flight (TOF) spectra from the theoretical reaction probability curves and we directly compare them to the raw experimental data. The results show that the use of more flexible functional forms for fitting the raw TOF spectra gives fits that are in slightly better agreement with the raw data and in considerably better agreement with theory, even though the theoretical reaction probabilities still achieve higher values at high energies. The mean absolute error (MAE) for the energy E0 at which the reaction probability equals half the experimental saturation value is now lower than 1 kcal/mol, the limit that defines chemical accuracy, while a MAE of 1.5 kcal/mol was previously obtained. The new AIMD results are only slightly different from the previous static surface results and in slightly better agreement with experiment.

The importance of accurate adiabatic interaction potentials for the correct description of electronically nonadiabatic vibrational energy transfer: a combined experimental and theoretical study of NO(v = 3) collisions with a Au(111) surface.

We present a combined experimental and theoretical study of NO(v = 3 → 3, 2, 1) scattering from a Au(111) surface at incidence translational energies ranging from 0.1 to 1.2 eV. Experimentally, molecular beam-surface scattering is combined with vibrational overtone pumping and quantum-state selective detection of the recoiling molecules. Theoretically, we employ a recently developed first-principles approach, which employs an Independent Electron Surface Hopping (IESH) algorithm to model the nonadiabatic dynamics on a Newns-Anderson Hamiltonian derived from density functional theory. This approach has been successful when compared to previously reported NO/Au scattering data. The experiments presented here show that vibrational relaxation probabilities increase with incidence energy of translation. The theoretical simulations incorrectly predict high relaxation probabilities at low incidence translational energy. We show that this behavior originates from trajectories exhibiting multiple bounces at the surface, associated with deeper penetration and favored (N-down) molecular orientation, resulting in a higher average number of electronic hops and thus stronger vibrational relaxation. The experimentally observed narrow angular distributions suggest that mainly single-bounce collisions are important. Restricting the simulations by selecting only single-bounce trajectories improves agreement with experiment. The multiple bounce artifacts discovered in this work are also present in simulations employing electronic friction and even for electronically adiabatic simulations, meaning they are not a direct result of the IESH algorithm. This work demonstrates how even subtle errors in the adiabatic interaction potential, especially those that influence the interaction time of the molecule with the surface, can lead to an incorrect description of electronically nonadiabatic vibrational energy transfer in molecule-surface collisions.

Dynamical steering in an electron transfer surface reaction: oriented NO(v = 3, 0.08 < Ei < 0.89 eV) relaxation in collisions with a Au(111) surface.

We report measurements of the incidence translational energy dependence of steric effects in collisions of NO(v = 3) molecules with a Au(111) surface using a recently developed technique to orient beams of vibrationally excited NO molecules at incidence energies of translation between 0.08 and 0.89 eV. Incidence orientation dependent vibrational state distributions of scattered molecules are detected by means of resonance enhanced multiphoton ionization spectroscopy. Molecules oriented with the N-end towards the surface exhibit a higher vibrational relaxation probability than those oriented with the O-end towards the surface. This strong orientation dependence arises from the orientation dependence of the underlying electron transfer reaction responsible for the vibrational relaxation. At reduced incidence translational energy, we observe a reduced steric effect. This reflects dynamical steering and re-orientation of the NO molecule upon its approach to the surface.

Ab initio molecular dynamics calculations on scattering of hyperthermal H atoms from Cu(111) and Au(111).

Energy loss from the translational motion of an atom or molecule impinging on a metal surface to the surface may determine whether the incident particle can trap on the surface, and whether it has enough energy left to react with another molecule present at the surface. Although this is relevant to heterogeneous catalysis, the relative extent to which energy loss of hot atoms takes place to phonons or electron-hole pair (ehp) excitation, and its dependence on the system's parameters, remain largely unknown. We address these questions for two systems that present an extreme case of the mass ratio of the incident atom to the surface atom, i.e., H + Cu(111) and H + Au(111), by presenting adiabatic ab initio molecular dynamics (AIMD) predictions of the energy loss and angular distributions for an incidence energy of 5 eV. The results are compared to the results of AIMDEFp calculations modeling energy loss to ehp excitation using an electronic friction ("EF") model applied to the AIMD trajectories, so that the energy loss to the electrons is calculated "post" ("p") the computation of the AIMD trajectory. The AIMD calculations predict average energy losses of 0.38 eV for Cu(111) and 0.13-0.14 eV for Au(111) for H-atoms that scatter from these surfaces without penetrating the surface. These energies closely correspond with energy losses predicted with Baule models, which is suggestive of structure scattering. The predicted adiabatic integral energy loss spectra (integrated over all final scattering angles) all display a lowest energy peak at an energy corresponding to approximately 80% of the average adiabatic energy loss for non-penetrative scattering. In the adiabatic limit, this suggests a way of determining the approximate average energy loss of non-penetratively scattered H-atoms from the integral energy loss spectrum of all scattered H-atoms. The AIMDEFp calculations predict that in each case the lowest energy loss peak should show additional energy loss in the range 0.2-0.3 eV due to ehp excitation, which should be possible to observe. The average non-adiabatic energy losses for non-penetrative scattering exceed the adiabatic losses to phonons by 0.9-1.0 eV. This suggests that for scattering of hyperthermal H-atoms from coinage metals the dominant energy dissipation channel should be to ehp excitation. These predictions can be tested by experiments that combine techniques for generating H-atom beams that are well resolved in translational energy and for detecting the scattered atoms with high energy-resolution.