Formation and dynamics of hot-precursor hydrogen atoms on metal surfaces: Trajectory simulations and stochastic models

Hydrogen abstraction from metal surfaces: when electron-hole pair excitations strongly affect hot-atom recombination

Using molecular dynamics simulations, we predict that the inclusion of nonadiabatic electronic excitations influences the dynamics of preadsorbed hydrogen abstraction from the W(110) surface by hydrogen scattering. The hot-atom recombination, which involves hyperthermal diffusion of the impinging atom on the surface, is significantly affected by the dissipation of energy mediated by electron-hole pair excitations at low coverage and low incidence energy. This issue is of importance as this abstraction mechanism is thought to largely contribute to molecular hydrogen formation from metal surfaces.

Energy dissipation to tungsten surfaces upon hot-atom and Eley–Rideal recombination of H2

Adiabatic and nonadiabatic quasi-classical molecular dynamics simulations are performed to investigate the role of electron–hole pair excitations in hot-atom and Eley–Rideal H₂ recombination mechanisms on H-covered W(100). The influence of the surface structure is analyzed by comparing with previous results for W(110).

Electronic friction dominates hydrogen hot-atom relaxation on Pd(100)

We study the dynamics of transient hot H atoms on Pd(100) that originated from dissociative adsorption of H2. The methodology developed here, denoted AIMDEF, consists of ab initio molecular dynamics simulations that include a friction force to account for the energy transfer to the electronic system. We find that the excitation of electron-hole pairs is the main channel for energy dissipation, which happens at a rate that is five times faster than energy transfer into Pd lattice motion. Our results show that electronic excitations may constitute the dominant dissipation channel in the relaxation of hot atoms on surfaces.

The formation of ice mantles on interstellar grains revisited – the effect of exothermicity

Modelling of grain surface chemistry generally deals with the simulation of rare events. Usually deterministic methods or statistical approaches such as the kinetic Monte Carlo technique are applied for these simulations. All assume that the surface processes are memoryless, the Markov chain assumption, and usually also that their rates are time independent. In this paper we investigate surface reactions for which these assumptions are not valid, and discuss what the effect is on the formation of water on interstellar grains. We will particularly focus on the formation of two OH radicals by the reaction H + HO₂. Two reaction products are formed in this exothermic reaction and the resulting momentum gained causes them to move away from each other. What makes this reaction special is that the two products can undergo a follow-up reaction to form H₂O₂. Experimentally, OH has been observed, which means that the follow-up reaction does not proceed with 100% efficiency, even though the two OH radicals are formed in each other's vicinity in the same reaction. This can be explained by a combined effect of the directionality of the OH radical movement together with energy dissipation. Both effects are constrained by comparison with experiments, and the resulting parametrised mechanism is applied to simulations of the formation of water ice under interstellar conditions.

Surface temperature effects on the dynamics of N2 Eley-Rideal recombination on W(100)

Quasiclassical trajectories simulations are performed to study the influence of surface temperature on the dynamics of a N atom colliding a N-preadsorbed W(100) surface under normal incidence. A generalized Langevin surface oscillator scheme is used to allow energy transfer between the nitrogen atoms and the surface. The influence of the surface temperature on the N(2) formed molecules via Eley-Rideal recombination is analyzed at T = 300, 800, and 1500 K. Ro-vibrational distributions of the N(2) molecules are only slightly affected by the presence of the thermal bath whereas kinetic energy is rather strongly decreased when going from a static surface model to a moving surface one. In terms of reactivity, the moving surface model leads to an increase of atomic trapping cross section yielding to an increase of the so-called hot atoms population and a decrease of the direct Eley-Rideal cross section. The energy exchange between the surface and the nitrogen atoms is semi-quantitatively interpreted by a simple binary collision model.

Relaxation of hot atoms following H2 dissociation on a Pd(111) surface

We study the relaxation of hot H atoms produced by dissociation of H2 molecules on the Pd111 surface. Ab initio density-functional theory calculations and the "corrugation reducing procedure" are used to determine the interaction potential for a H atom in front of a rigid surface as well as its modification under surface-atom vibrations. A slab of 80 Pd atoms is used to model the surface together with "generalized Langevin oscillators" to account for energy dissipation to the bulk. We show that the energy relaxation is fast, about 75% of the available energy being lost by the hot atoms after 0.5 ps. As a consequence, the hot atoms do not travel more than a few angstroms along the surface before being trapped into the potential well located over the hollow site.

Quasiclassical study of Eley–Rideal and hot atom reactions of H atoms with Cl adsorbed on a Au(111) surface

Using quasiclassical methods and a potential energy surface based on total energy calculations, we have found that H atoms react with Cl atoms adsorbed onto a Au(111) surface to produce HCl via Eley-Rideal (ER), hot atom (HA), and Langmuir-Hinschelwood (LH) pathways. We observe two ER mechanisms. At small normal incidence energies reaction results from a more or less direct collision with Cl, leading to a large amount of product vibration (nu=8), and relatively cold rotation and translation. In the second mechanism, more dominant at near-normal incidence and/or large incident energies, the H atom passes near Cl, recoils from the metal, and is pulled into orbit about Cl. This leads to broader product state distributions, and a more even distribution of the 3.0 eV of available energy among the product degrees of freedom, similar to products formed via the HA pathway. Overall, ER processes tend to contribute less than 10% to the reactivity, and most of the HCl is formed via HA processes. There is an increase in HCl formation with surface temperature for both the ER and HA mechanisms, but this increase is relatively weak. We observe typically about 12% H atom sticking, which would lead to HCl formation via a LH process in the experiments, above 140 K. We observe a weak forward scattering due to the direct ER component, as in the experiments. However, unlike the experiments, we observe a dip in our product angular distributions about thetaf=0 degrees, which we ascribe to our quasiclassical approximation. While we tend to see more energy in the hot products than in the experiments, our product translational, rotational, and vibrational distributions are in relatively reasonable agreement with those measured. One major disagreement with experiment is that there is apparently a significant sticking of the H atom at low temperatures, leading to a large LH component. In addition, the ER and HA components increase much more strongly with temperature than in the calculations. It is possible that electon-hole pair excitations in the metal strongly relax both the H atom and the excited HCl molecules formed.

Hot-atom versus Eley–Rideal dynamics in hydrogen recombination on Ni(100). I. The single-adsorbate case

We compare the efficiency of the Eley-Rideal (ER) reaction with the formation of hot-atom (HA) species in the simplest case, i.e., the scattering of a projectile off a single adsorbate, considering the Hydrogen and Hydrogen-on-Ni(100) system. We use classical mechanics and the accurate embedded diatomics-in-molecules potential to study the collision system over a wide range of collision energies (0.10-1.50 eV), both with a rigid and a nonrigid Ni substrate and for impact on the occupied and neighboring empty cells. In the rigid model metastable and truly bound hot-atoms occur and we find that the cross section for the formation of bound hot-atoms is considerably higher than that for the ER reaction over the whole range of collision energies examined. Metastable hot-atoms form because of the inefficient energy transfer to the adsorbate and have lifetimes of the order 0.1-0.7 ps, depending on the collision energy. When considering the effects of lattice vibrations we find, on average, a consistent energy transfer to the substrate, say 0.1-0.2 eV, which forced us to devise a two-step dynamical model to get rid of the problems associated with the use of periodic boundary conditions. Results for long-lived HA formation due to scattering on the occupied cell at a surface temperature of 120 K agree well with those of the rigid model, suggesting that in the above process the substrate plays only a secondary role and further calculations at surface temperatures of 50 and 300 K are in line with these findings. However, considerably high cross sections for formation of long-lived hot-atoms result also from scattering off the neighboring cells where the energy transfer to the lattice cannot be neglected. Metastable hot-atoms are reduced in number and have usually lifetimes shorter than those of the rigid-model, say less than 0.3 ps. In addition, ER cross sections are only slightly affected by the lattice motion and show a little temperature dependence. Finally, we find also that absorption and reflection strongly depend on the correct consideration of lattice vibrations and the occurrence of trapping.

Efficient Eley-Rideal reactions of H atoms with single Cl adsorbates on Au(111)

Density functional theory is used to construct an interaction model for H atoms with Cl over Au(111). Single-adsorbate Eley-Rideal reactions are investigated with quantum and quasiclassical methods. The reaction cross sections, amounting to 2-3 A(2), are much larger than for HD recombinations on metals. This can be traced to the adsorbed Cl being relatively far above the surface, the H-Cl interaction prevailing over the H-substrate attraction for a sizable range of impact parameters.