Condensed-phase isomerization through tunnelling gateways.

Quantum mechanical tunnelling describes transmission of matter waves through a barrier with height larger than the energy of the wave1. Tunnelling becomes important when the de Broglie wavelength of the particle exceeds the barrier thickness; because wavelength increases with decreasing mass, lighter particles tunnel more efficiently than heavier ones. However, there exist examples in condensed-phase chemistry where increasing mass leads to increased tunnelling rates2. In contrast to the textbook approach, which considers transitions between continuum states, condensed-phase reactions involve transitions between bound states of reactants and products. Here this conceptual distinction is highlighted by experimental measurements of isotopologue-specific tunnelling rates for CO rotational isomerization at an NaCl surface3,4, showing nonmonotonic mass dependence. A quantum rate theory of isomerization is developed wherein transitions between sub-barrier reactant and product states occur through interaction with the environment. Tunnelling is fastest for specific pairs of states (gateways), the quantum mechanical details of which lead to enhanced cross-barrier coupling; the energies of these gateways arise nonsystematically, giving an erratic mass dependence. Gateways also accelerate ground-state isomerization, acting as leaky holes through the reaction barrier. This simple model provides a way to account for tunnelling in condensed-phase chemistry, and indicates that heavy-atom tunnelling may be more important than typically assumed.

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.

Vibrational energy transfer in collisions of molecules with metal surfaces.

The Born-Oppenheimer approximation (BOA), which serves as the basis for our understanding of chemical bonding, reactivity and dynamics, is routinely violated for vibrationally inelastic scattering of molecules at metal surfaces. The title-field therefore represents a fascinating challenge to our conventional wisdom calling for new concepts that involve explicit electron dynamics occurring in concert with nuclear motion. Here, we review progress made in this field over the last decade, which has witnessed dramatic advances in experimental methods, thereby providing a much more extensive set of diverse observations than has ever before been available. We first review the experimental methods used in this field and then provide a systematic tour of the vast array of observations that are currently available. We show how these observations - taken together and without reference to computational simulations - lead us to a simple and intuitive picture of BOA failure in molecular dynamics at metal surfaces, one where electron transfer between the molecule and the metal plays a preeminent role. We also review recent progress made in the theory of electron transfer mediated BOA failure in molecule-surface interactions, describing the most important methods and their ability to reproduce experimental observation. Finally, we outline future directions for research and important unanswered questions.

Vibrationally Mode-Specific Molecular Energy Transfer to Surface Electrons in Metastable Formaldehyde Scattering from Cesium-Covered Au(111).

Nonadiabatic interaction of adsorbate nuclear motion with the continuum of electronic states is known to affect the dynamics of chemical reactions at metal surfaces. A large body of work has probed the fundamental mechanisms of such interactions for atomic and diatomic molecules at surfaces. In polyatomic molecules, the possibility of mode-specific damping of vibrational motion due to the effects of electronic friction raises the question of whether such interactions could profoundly affect the outcome of chemistry at surfaces by selectively removing energy from a particular intramolecular adsorbate mode. However, to date, there have not been any fundamental experiments demonstrating nonadiabatic electron-vibration coupling in a polyatomic molecule at a surface. In this work, we scatter excited metastable formaldehyde and formaldehyde-d2 from a low work function surface and detect ejected exoelectrons that accompany molecular relaxation. The exoelectron ejection efficiency exhibits a strong dependence on the vibrational mode that is excited: out-of-plane bending excitation (ν4) leads to significantly more exoelectrons than does CO stretching excitation (ν2). The results provide clear evidence for mode-specific energy transfer from vibration to surface electrons.

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 and classical molecular dynamics for H atom scattering from graphene.

This work presents systematic comparisons between classical molecular dynamics (cMD) and quantum dynamics (QD) simulations of 15-dimensional and 75-dimensional models in their description of H atom scattering from graphene. We use an experimentally validated full-dimensional neural network potential energy surface of a hydrogen atom interacting with a large cell of graphene containing 24 carbon atoms. For quantum dynamics simulations, we apply Monte Carlo canonical polyadic decomposition to transform the original potential energy surface (PES) into a sum of products form and use the multi-layer multi-configuration time-dependent Hartree method to simulate the quantum scattering of a hydrogen or deuterium atom with an initial kinetic energy of 1.96 or 0.96 eV and an incident angle of 0°, i.e., perpendicular to the graphene surface. The cMD and QD initial conditions have been carefully chosen in order to be as close as possible. Our results show little differences between cMD and QD simulations when the incident energy of the H atom is equal to 1.96 eV. However, a large difference in sticking probability is observed when the incident energy of the H atom is equal to 0.96 eV, indicating the predominance of quantum effects. To the best of our knowledge, our work provides the first benchmark of quantum against classical simulations for a system of this size with a realistic PES. Additionally, new projectors are implemented in the Heidelberg multi-configuration time-dependent Hartree package for the calculation of the atom scattering energy transfer distribution as a function of outgoing angles.

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.

Effective medium theory for bcc metals: electronically non-adiabatic H atom scattering in full dimensions.

We report a newly derived Effective Medium Theory (EMT) formalism for bcc metals and apply it for the construction of a full-dimensional PES for H atoms interacting with molybdenum (Mo) and tungsten (W). We construct PESs for the (111) and (110) facets of both metals. The EMT-PESs have the advantage that they automatically provide the background electron density on the fly which makes incorporation of ehp excitation within the framework of electronic friction straightforward. Using molecular dynamics with electronic friction (MDEF) with these new PESs, we simulated 2.76 eV H atoms scattering and adsorption. The large energy losses at a surface temperature of 300 K is very similar those seen for H atom scattering from the late fcc metals and is dominated by ehp excitation. We see significant differences in the scattering from different surface facets of the same metal. For the (110) facet, we see strong evidence of sub-surface scattering, which should be observable in experiment and we predict the best conditions for observing this novel type of scattering process. At low temperatures the MD simulations predict that H atom scattering is surface specific due to the reduced influence of the random force.

H atom scattering from W(110): A benchmark for molecular dynamics with electronic friction.

Molecular dynamics with electronic friction (MDEF) at the level of the local density friction approximation (LDFA) has been applied to describe electronically non-adiabatic energy transfer accompanying H atom collisions with many solid metal surfaces. When implemented with full dimensional potential energy and electron density functions, excellent agreement with experiment is found. Here, we compare the performance of a reduced dimensional MDEF approach involving a simplified description of H atom coupling to phonons to that of full dimensional MDEF calculations known to yield accurate results. Both approaches give remarkably similar results for H atom energy loss distributions with a 300 K W(110) surface. At low surface temperature differences are seen; but, quantities like average energy loss are still accurately reproduced. Both models predict similar conditions under which H atoms that have penetrated into the subsurface regions could be observed in scattering experiments.

Stereodynamics of adiabatic and non-adiabatic energy transfer in a molecule surface encounter.

Molecular energy transfer and reactions at solid surfaces depend on the molecular orientation relative to the surface. While such steric effects have been largely understood in electronically adiabatic processes, the orientation-dependent energy transfer in NO scattering from Au(111) was complicated by electron-mediated nonadiabatic effects, thus lacking a clear interpretation and posing a great challenge for theories. Herein, we investigate the stereodynamics of adiabatic and nonadiabatic energy transfer via molecular dynamics simulations of NO(v = 3) scattering from Au(111) using realistic initial orientation distributions based on accurate neural network fitted adiabatic potential energy surface and electronic friction tensor. Our results reproduce the observed stronger vibrational relaxation for N-first orientation and enhanced rotational rainbow for O-first orientation, and demonstrate how adiabatic anisotropic interactions steer molecules into the more attractive N-first orientation to experience more significant energy transfer. Remaining disagreements with experiment suggest the direction for further developments of nonadiabatic theories for gas-surface scattering.

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.

Inelastic Scattering of H Atoms from Surfaces.

We have developed an instrument that uses photolysis of hydrogen halides to produce nearly monoenergetic hydrogen atom beams and Rydberg atom tagging to obtain accurate angle-resolved time-of-flight distributions of atoms scattered from surfaces. The surfaces are prepared under strict ultrahigh vacuum conditions. Data from these experiments can provide excellent benchmarks for theory, from which it is possible to obtain an atomic scale understanding of the underlying dynamical processes governing H atom adsorption. In this way, the mechanism of adsorption on metals is revealed, showing a penetration-resurfacing mechanism that relies on electronic excitation of the metal by the H atom to succeed. Contrasting this, when H atoms collide at graphene surfaces, the dynamics of bond formation involving at least four carbon atoms govern adsorption. Future perspectives of H atom scattering from surfaces are also outlined.

Multibounce and Subsurface Scattering of H Atoms Colliding with a van der Waals Solid.

We report the results of inelastic differential scattering experiments and full-dimensional molecular dynamics trajectory simulations for 2.76 eV H atoms colliding at a surface of solid xenon. The interaction potential is based on an effective medium theory (EMT) fit to density functional theory (DFT) energies. The translational energy-loss distributions derived from experiment and theory are in excellent agreement. By analyzing trajectories, we find that only a minority of the scattering results from simple single-bounce dynamics. The majority comes from multibounce collisions including subsurface scattering where the H atoms penetrate below the first layer of Xe atoms and subsequently re-emerge to the gas phase. This behavior leads to observable energy-losses as large as 0.5 eV, much larger than a prediction of the binary collision model (0.082 eV), which is often used to estimate the highest possible energy-loss in direct inelastic surface scattering. The sticking probability computed with the EMT-PES (0.15) is dramatically reduced (5 × 10-6) if we employ a full-dimensional potential energy surface (PES) based on Lennard-Jones (LJ) pairwise interactions. Although the LJ-PES accurately describes the interactions near the H-Xe and Xe-Xe energy minima, it drastically overestimates the effective size of the Xe atom seen by the colliding H atom at incidence energies above about 0.1 eV.

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.

Unified description of H-atom-induced chemicurrents and inelastic scattering.

The Born-Oppenheimer approximation (BOA) provides the foundation for virtually all computational studies of chemical binding and reactivity, and it is the justification for the widely used "balls and springs" picture of molecules. The BOA assumes that nuclei effectively stand still on the timescale of electronic motion, due to their large masses relative to electrons. This implies electrons never change their energy quantum state. When molecules react, atoms must move, meaning that electrons may become excited in violation of the BOA. Such electronic excitation is clearly seen for: (i) Schottky diodes where H adsorption at Ag surfaces produces electrical "chemicurrent;" (ii) Au-based metal-insulator-metal (MIM) devices, where chemicurrents arise from H-H surface recombination; and (iii) Inelastic energy transfer, where H collisions with Au surfaces show H-atom translation excites the metal's electrons. As part of this work, we report isotopically selective hydrogen/deuterium (H/D) translational inelasticity measurements in collisions with Ag and Au. Together, these experiments provide an opportunity to test new theories that simultaneously describe both nuclear and electronic motion, a standing challenge to the field. Here, we show results of a recently developed first-principles theory that quantitatively explains both inelastic scattering experiments that probe nuclear motion and chemicurrent experiments that probe electronic excitation. The theory explains the magnitude of chemicurrents on Ag Schottky diodes and resolves an apparent paradox--chemicurrents exhibit a much larger isotope effect than does H/D inelastic scattering. It also explains why, unlike Ag-based Schottky diodes, Au-based MIM devices are insensitive to H adsorption.

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.

An experimentally validated neural-network potential energy surface for H-atom on free-standing graphene in full dimensionality.

We present a first principles-quality potential energy surface (PES) describing the inter-atomic forces for hydrogen atoms interacting with free-standing graphene. The PES is a high-dimensional neural network potential that has been parameterized to 75 945 data points computed with density-functional theory employing the PBE-D2 functional. Improving over a previously published PES [Jiang et al., Science, 2019, 364, 379], this neural network exhibits a realistic physisorption well and achieves a 10-fold reduction in the RMS fitting error, which is 0.6 meV per atom. The chemisorption barrier is 172 meV, which is lower than that of the REBO-EMFT PES (260 meV). We used this PES to calculate about 1.5 million classical trajectories with carefully selected initial conditions to allow for direct comparison to results of H- and D-atom scattering experiments performed at incidence translational energy of 1.9 eV and a surface temperature of 300 K. The theoretically predicted scattering angular and energy loss distributions are in good agreement with experiment, despite the fact that the experiments employed graphene grown on Pt(111). Compared to previous calculations, the agreement with experiments is improved. The remaining discrepancies between experiment and theory are likely due to the influence of the Pt substrate only present in the experiment.

Vibrational Relaxation Lifetime of a Physisorbed Molecule at a Metal Surface.

Previous measurements of vibrational relaxation lifetimes for molecules adsorbed at metal surfaces yielded values of 1-3 ps; however, only chemisorbed molecules have been studied. We report the first measurements of the vibrational relaxation lifetime of a molecule physisorbed to a metal surface. For CO(v=1) adsorbed on Au(111) at 35 K the vibrational lifetime of the excited stretching mode is 49±3 ps. The long lifetime seen here is likely to be a general feature of physisorption, which involves weaker electronic coupling between the adsorbate and the solid due to bonding at larger distances.