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.

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.

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.

Adsorption and Absorption Energies of Hydrogen with Palladium.

Thermal recombinative desorption rates of HD on Pd(111) and Pd(332) are reported from transient kinetic experiments performed between 523 and 1023 K. A detailed kinetic model accurately describes the competition between recombination of surface-adsorbed hydrogen and deuterium atoms and their diffusion into the bulk. By fitting the model to observed rates, we derive the dissociative adsorption energies (E 0, ads H2 = 0.98 eV; E 0, ads D2 = 1.00 eV; E 0, ads HD = 0.99 eV) as well as the classical dissociative binding energy ϵads = 1.02 ± 0.03 eV, which provides a benchmark for electronic structure theory. In a similar way, we obtain the classical energy required to move an H or D atom from the surface to the bulk (ϵsb = 0.46 ± 0.01 eV) and the isotope specific energies, E 0, sb H = 0.41 eV and E 0, sb D = 0.43 eV. Detailed insights into the process of transient bulk diffusion are obtained from kinetic Monte Carlo simulations.

Quantum effects in thermal reaction rates at metal surfaces.

There is wide interest in developing accurate theories for predicting rates of chemical reactions that occur at metal surfaces, especially for applications in industrial catalysis. Conventional methods contain many approximations that lack experimental validation. In practice, there are few reactions where sufficiently accurate experimental data exist to even allow meaningful comparisons to theory. Here, we present experimentally derived thermal rate constants for hydrogen atom recombination on platinum single-crystal surfaces, which are accurate enough to test established theoretical approximations. A quantum rate model is also presented, making possible a direct evaluation of the accuracy of commonly used approximations to adsorbate entropy. We find that neglecting the wave nature of adsorbed hydrogen atoms and their electronic spin degeneracy leads to a 10× to 1000× overestimation of the rate constant for temperatures relevant to heterogeneous catalysis. These quantum effects are also found to be important for nanoparticle catalysts.

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.

NO Binding Energies to and Diffusion Barrier on Pd Obtained with Velocity-Resolved Kinetics.

We report nitric oxide (NO) desorption rates from Pd(111) and Pd(332) surfaces measured with velocity-resolved kinetics. The desorption rates at the surface temperatures from 620 to 800 K span more than 3 orders of magnitude, and competing processes, like dissociation, are absent. Applying transition state theory (TST) to model experimental data leads to the NO binding energy E 0 = 1.766 ± 0.024 eV and diffusion barrier D T = 0.29 ± 0.11 eV on the (111) terrace and the stabilization energy for (110)-steps ΔE ST = 0.060-0.030 +0.015 eV. These parameters provide valuable benchmarks for theory.

Adsorbate modification of electronic nonadiabaticity: H atom scattering from p(2 × 2) O on Pt(111).

We report inelastic differential scattering experiments for energetic H and D atoms colliding at a Pt(111) surface with and without adsorbed O atoms. Dramatically, more energy loss is seen for scattering from the Pt(111) surface compared to p(2 × 2) O on Pt(111), indicating that O adsorption reduces the probability of electron-hole pair (EHP) excitation. We produced a new full-dimensional potential energy surface for H interaction with O/Pt that reproduces density functional theory energies accurately. We then attempted to model the EHP excitation in H/D scattering with molecular dynamics simulations employing the electronic density information from the Pt(111) to calculate electronic friction at the level of the local density friction approximation (LDFA). This approach, which assumes that O atoms simply block the Pt atom from the approaching H atom, fails to reproduce experiment due to the fact that the effective collision cross section of the O atom is only 10% of the area of the surface unit cell. An empirical adiabatic sphere model that reduces electronic nonadiabaticity within an O-Pt bonding length scale of 2.8 Å reproduces experiment well, suggesting that the electronic structure changes induced by chemisorption of O atoms nearly remove the H atom's ability to excite EHPs in the Pt. Alternatives to LDFA friction are needed to account for this adsorbate effect.

Following the microscopic pathway to adsorption through chemisorption and physisorption wells.

Adsorption involves molecules colliding at the surface of a solid and losing their incidence energy by traversing a dynamical pathway to equilibrium. The interactions responsible for energy loss generally include both chemical bond formation (chemisorption) and nonbonding interactions (physisorption). In this work, we present experiments that revealed a quantitative energy landscape and the microscopic pathways underlying a molecule's equilibration with a surface in a prototypical system: CO adsorption on Au(111). Although the minimum energy state was physisorbed, initial capture of the gas-phase molecule, dosed with an energetic molecular beam, was into a metastable chemisorption state. Subsequent thermal decay of the chemisorbed state led molecules to the physisorption minimum. We found, through detailed balance, that thermal adsorption into both binding states was important at all temperatures.

Fundamental mechanisms for molecular energy conversion and chemical reactions at surfaces.

The dream of theoretical surface chemistry is to predict the outcome of reactions in order to find the ideal catalyst for a certain application. Having a working ab initio theory in hand would not only enable these predictions but also provide insights into the mechanisms of surface reactions. The development of theoretical models can be assisted by experimental studies providing benchmark data. Though for some reactions a quantitative agreement between experimental observations and theoretical calculations has been achieved, theoretical surface chemistry is in general still far away from gaining predictive power. Here we review recent experimental developments towards the understanding of surface reactions. It is demonstrated how quantum-state resolved scattering experiments on reactive and nonreactive systems can be used to test front-running theoretical approaches. Two challenges for describing dynamics at surfaces are addressed: nonadiabaticity in diatomic molecule surface scattering and the increasing system size when observing and describing the dynamics of polyatomic molecules at surfaces. Finally recent experimental studies on reactive systems are presented. It is shown how elementary steps in a complex surface reaction can be revealed experimentally.

Roaming Dynamics in the Photodissociation of Formic Acid at 230 nm.

Roaming dynamics is observed in the photodissociation of formic acid (HCOOH) at 230 nm by using the slice imaging method. In combination with rotational state selective (2 + 1) resonance-enhanced multiphoton ionization of the CO fragments, the speed distributions of the CO fragments exhibit a low recoil velocity at low rotational levels of J = 9 and 20, while the velocity distributions of CO at high rotational levels of J = 30 and 48 show a relatively large recoil velocity. The experimental results indicate that the roaming of OH radical should be related with the formation of CO + H2O channel at the present photolysis energy. Unlike the roaming pathways occurring in H2CO that can be described by loose flat potential, our CO speed distribution analysis suggests the presence of a "tight" flat potential in the roaming dynamics of HCOOH molecules.

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.

Assessment of UVA-Riboflavin Corneal Cross-Linking Using Small Amplitude Oscillatory Shear Measurements.

PURPOSE: The effect of ultraviolet (UV)-riboflavin cross-linking (CXL) has been measured primarily using the strip extensometry technique. We propose a simple and reliable methodology for the assessment of CXL treatment by using an established rheologic protocol based on small amplitude oscillatory shear (SAOS) measurements. It provides information on the average cross-link density and the elastic modulus of treated cornea samples. METHODS: Three fresh postmortem porcine corneas were used to study the feasibility of the technique, one serving as control and two receiving corneal collagen cross-linking treatment. Subsequently, five pairs of fresh postmortem porcine corneas received corneal collagen cross-linking treatment with riboflavin and UVA-irradiation (370 nm; irradiance of 3 mW/cm2) for 30 minutes (Dresden protocol); the contralateral porcine corneas were used as control samples. After the treatment, the linear viscoelastic moduli of the corneal samples were measured using SAOS measurements and the average cross-linking densities extracted. RESULTS: For all cases investigated, the dynamic moduli of the cross-linked corneas were higher compared to those of the corresponding control samples. The increase of the elastic modulus of the treated samples was between 122% and 1750%. The difference was statistically significant for all tested samples (P = 0.018, 2-tailed t-test). CONCLUSIONS: We report a simple and accurate methodology for quantifying the effects of cross-linking on porcine corneas treated with the Dresden protocol by means of SAOS measurements in the linear regime. The measured dynamic moduli, elastic and viscous modulus, represent the energy storage and energy dissipation, respectively. Hence, they provide a means to assess the changing physical properties of the cross-linked collagen networks after CXL treatment.

New insights into the photodissociation of methyl iodide at 193 nm: stereodynamics and product branching ratios.

The stereodynamics of methyl iodide photodissociation after excitation at 193 nm has been studied using a combination of slice imaging and resonance enhanced multiphoton ionization (REMPI) detection of the methyl and iodine products. A weak anisotropic ring appearing in the image corresponding to vibrationally excited CH3(ν1 = 1) confirms the production of ground state I((2)P3/2) atoms at this excitation wavelength as a signature of the predissociation channel reported previously [M. G. González et al., J. Chem. Phys., 2011, 135, 021102] tentatively assigned to the coupling between the B-band (3)R1 Rydberg state and the A-band (1)Q1 repulsive state. Direct REMPI detection of ground state iodine atoms indicates that most of the I((2)P3/2) species are produced in correlation with highly internally excited methyl radicals, in excellent agreement with the recent results of Xu and Pratt [Xu et al., J. Chem. Phys., 2013, 139, 214310; Xu et al., J. Phys. Chem. A, 2015, 119, 7548]. From the comparison between the CH3(ν) second order Dixon's bipolar moments ß(2)(0)(20), ß(0)(0)(22), ß(2)(0)(02) and ß(2)(0)(22) measured in this work and those reported previously for the B-band origin and the A-band, a general picture of the CH3I photodissociation stereodynamics in terms of different effects, such as the breakdown of the unique recoil direction (URD) approximation, the non-adiabatic curve crossings and the depolarization induced by the parent molecule rotation, is drawn.

Slice imaging of methyl bromide photofragmentation at 193 nm.

The photodissociation of methyl bromide at 193 nm is studied using slice imaging. From the measured photofragment translational energy and angular distributions we were able to extract methyl-state-specific dissociation channel yields and crossing probabilities between (3)Q0 and (1)Q1 surfaces. The angular distributions for the Bromine photofragments show a strong dependence on the total translational energy release. Nonadiabatic transition probabilities from the (3)Q0 to (1)Q1 surface dominate the dynamics in this excitation energy region and for most of the methyl vibrational states probed.

Observation of direct vibrational excitation in gas-surface collisions of CO with Au(111): a new model system for surface dynamics.

We report vibrational excitation of CO from its ground (v = 0) to first excited (v = 1) vibrational state in collision with Au(111) at an incidence energy of translation of E(I) = 0.45 eV. Unlike past work, we can exclude an excitation mechanism involving temporary adsorption on the surface followed by thermalization and desorption. The angular distributions of the scattered CO molecules are narrow, consistent with direct scattering occurring on a sub-ps time scale. The absolute excitation probabilities are about 3% of those expected from thermal accommodation. The surface temperature dependence of excitation, which was measured between 373 and 973 K, is Arrhenius-like with an activation energy equal to the energy required for vibrational excitation. Our measurements are consistent with a vibrational excitation mechanism involving coupling of thermally excited electron-hole pairs of the solid to CO vibration.

Driving photochemistry by clustering: the ICl-Xe case.

We present slice imaging data demonstrating the influence of clustering on the photodissociation dynamics of a diatomic molecule: iodine monochloride (ICl) was dissociated at 235 nm in He and Xe seed gasses, probing both Cl and I photofragment energy and angular distributions. We observe that the kinetic energy releases of both Cl and I fragments change from He to Xe seeding. For Cl fragments, the seeding in Xe increases the kinetic energy release of some Cl fragments with a narrow kinetic energy distribution, and leads to some fragments with rather broad statistical distribution falling off exponentially from near-zero energies up to about 2.5 eV. Iodine fragment distribution changes even more dramatically from He to Xe seeding: sharp features essentially disappear and a broad distribution arises reaching to about 2.5 eV. Both these observations are rationalized by a simple qualitative cluster model assuming ICl dissociation inside larger xenon clusters and "on surface" of smaller Xe species.

Ultraviolet photodissociation of iodine monochloride (ICl) at 235, 250, and 265 nm.

ICl photolysis in the ultraviolet region of the spectrum (235-265 nm) is studied using the Slice Imaging technique. The Cl∗((2)P(1/2))/Cl((2)P(3/2)) and the I∗((2)P(1/2))/I((2)P(3/2)) branching ratio between the I((2)P(3/2)) + Cl((2)P(3/2))∕Cl∗((2)P(1/2)) and I∗((2)P(1/2)) + Cl((2)P(3∕/2))∕Cl∗((2)P(1/2)) channels is extracted from the respective iodine and chlorine photofragment images. We find that ground state chlorine atoms (Cl((2)P(3/2))) are formed nearly exclusively with excited state iodine atoms (I∗((2)P(1/2))), while excited spin-orbit chlorine atoms (Cl∗((2)P(1/2))) are concurrently produced only with ground state iodine atoms (I((2)P(3/2))). We conclude that photolysis of ICl in this UV region is a relatively "clean" source of spin-orbit excited chlorine atoms that can be used in crossed molecular beam experiments.

Branching ratios and anisotropy parameters in ICl photolysis from 400 to 570 nm using slice imaging.

ICl photolysis in the visible region of the spectrum (400-570 nm) is studied using the slice imaging technique. The Cl((2)P(1/2))/Cl((2)P(3/2)) branching ratio between the I((2)P(3/2))+Cl((2)P(3/2)) and I((2)P(3/2))+Cl((2)P(1/2)) channels is extracted from the iodine photofragment images and it is found to range from 0 to 2.5, rising from 570 to 490 nm and dropping at higher photolysis energies. The I+Cl angular distribution exhibits a similar trend, changing from purely perpendicular at 570 nm to isotropic at 545 nm, fairly parallel at 490 nm and again perpendicular at 440 nm. Following previous work, we discuss these changes in light of avoided curve crossing and determine the crossing probability as a function of wavelength. The angular anisotropy parameter beta of the second channel ranges between 0.6 and 1.4.