Dynamics of H2 Eley-Rideal abstraction from W(110): Sensitivity to the representation of the molecule-surface potential

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.

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.

Random Force in Molecular Dynamics with Electronic Friction

Originally conceived to describe thermal diffusion, the Langevin equation includes both a frictional drag and a random force, the latter representing thermal fluctuations first seen as Brownian motion. The random force is crucial for the diffusion problem as it explains why friction does not simply bring the system to a standstill. When using the Langevin equation to describe ballistic motion, the importance of the random force is less obvious and it is often omitted, for example, in theoretical treatments of hot ions and atoms interacting with metals. Here, friction results from electronic nonadiabaticity (electronic friction), and the random force arises from thermal electron-hole pairs. We show the consequences of omitting the random force in the dynamics of H-atom scattering from metals. We compare molecular dynamics simulations based on the Langevin equation to experimentally derived energy loss distributions. Despite the fact that the incidence energy is much larger than the thermal energy and the scattering time is only about 25 fs, the energy loss distribution fails to reproduce the experiment if the random force is neglected. Neglecting the random force is an even more severe approximation than freezing the positions of the metal atoms or modelling the lattice vibrations as a generalized Langevin oscillator. This behavior can be understood by considering analytic solutions to the Ornstein-Uhlenbeck process, where a ballistic particle experiencing friction decelerates under the influence of thermal fluctuations.

Dissociative Chemisorption of O2 on Al(111): Dynamics on a Correlated Wave-Function-Based Potential Energy Surface

Dissociative chemisorption of O₂ on the Al(111) surface represents an extensively studied prototype for understanding the interaction between O₂ and metal surfaces. It is well known that the experimentally observed activation barrier for O₂ dissociation is not captured by conventional density functional theory. The interpretation of this barrier as a result of spin transitions along the reaction path has been challenged by recent embedded correlated wave function (ECW) calculations that naturally yield an adiabatic barrier. However, the ECW calculations have been limited to a static analysis of the reaction pathways and have not yet been tested by dynamics simulations. We present a global six-dimensional potential energy surface (PES) for this system parametrized with ECW data points. This new PES provides a reasonable description of the site-specific and orientation-dependent activation barriers. Quasi-classical trajectory calculations on this PES semiquantitatively reproduce both the observed translational energy dependence of the sticking probability and steric effects with aligned O₂ molecules.

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).