Dynamical reaction pathways in Eley-Rideal recombination of nitrogen from W(100)

Isotope effects in Eley-Rideal abstraction of hydrogen from tungsten surfaces: the role of dissipation

Molecular dynamics simulations are performed to investigate the influence of isotope substitutions on the Eley-Rideal recombination dynamics of hydrogen isotopes from the (100) and (110) surfaces of tungsten. Dissipation of electrons and phonons is taken into account by, respectively, the local density friction approximation and the general Langevin oscillator, effective models which have been intensively used in recent years. As the coupling to surface phonons and electrons might be altered by the mass combination, the main objective of the paper is to assess the role of dissipation to the surface in the course of abstraction.

Scattering of CO from Vacant-MoSe2 with O Adsorbates: Is CO2 Formed?

Using ab initio molecular dynamics (AIMD) simulations, based on density functional theory that also accounts for van der Waals interactions, we study the oxidation of gas-phase CO on MoSe2 with a Se vacancy and oxygen coverage of 0.125 ML. In the equilibrium configuration, one of the O atoms is adsorbed on the vacancy and the other one atop one Se atom. Recombination of the CO molecule with the second of these O atoms to form CO2 is a highly exothermic reaction, with an energy gain of around 3 eV. The likeliness of the CO oxidation reaction on this surface is next examined by calculating hundreds of AIMD trajectories for incidence energies that suffice to overcome the energy barriers in the entrance channel of the CO oxidative recombination. In spite of this, no CO2 formation event is obtained. In most of the calculated trajectories, the incoming CO molecule is directly reflected, and in some cases, mainly at low energies, the molecules remain trapped at the surface but without reacting. As an important conclusion, our AIMD simulations show that the recombination of CO molecules with adsorbed O atoms is a very unlikely reaction in this system, despite its large exothermicity.

Plasma Catalysis Modeling: How Ideal Is Atomic Hydrogen for Eley-Rideal?

Plasma catalysis is an emerging technology, but a lot of questions about the underlying surface mechanisms remain unanswered. One of these questions is how important Eley-Rideal (ER) reactions are, next to Langmuir-Hinshelwood reactions. Most plasma catalysis kinetic models predict ER reactions to be important and sometimes even vital for the surface chemistry. In this work, we take a critical look at how ER reactions involving H radicals are incorporated in kinetic models describing CO2 hydrogenation and NH3 synthesis. To this end, we construct potential energy surface (PES) intersections, similar to elbow plots constructed for dissociative chemisorption. The results of the PES intersections are in agreement with ab initio molecular dynamics (AIMD) findings in literature while being computationally much cheaper. We find that, for the reactions studied here, adsorption is more probable than a reaction via the hot atom (HA) mechanism, which in turn is more probable than a reaction via the ER mechanism. We also conclude that kinetic models of plasma-catalytic systems tend to overestimate the importance of ER reactions. Furthermore, as opposed to what is often assumed in kinetic models, the choice of catalyst will influence the ER reaction probability. Overall, the description of ER reactions is too much "ideal" in models. Based on our findings, we make a number of recommendations on how to incorporate ER reactions in kinetic models to avoid overestimation of their importance.

Dynamics in reactions on metal surfaces: A theoretical perspective

Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.

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

Kinetic analysis of interaction between N atoms and O-covered Ru(0001)

Eley-Rideal (ER) reactions involving neutral atoms heavier than hydrogen reacting with adsorbed atoms of similar mass were first observed in recent molecular beam experiments by Zaharia et al. [Phys. Rev. Lett. 113, 053201 (2014)]. Through analysis of two types of measurements, they obtained different estimations for the N-O ER reaction cross section, one of which is unexpectedly high. This was qualitatively accounted for by invoking a secondary effect whereby the presence of N adatoms on the surface acted to "shield" O adatoms from prompt recombinative desorption. We apply a rate equation model that includes two ER processes involving different adsorbed species (N-Oad and N-Nad) and an N-adsorption process to the full-beam exposure subset of the experimental data in order to study the reaction kinetics. Values for the individual reaction cross sections are derived. The measured N2 response can be well described by the model, but it is insufficient to completely describe the NO response. Modeling of different exposures is used to evaluate the qualitative picture presented by Zaharia et al.

Eley-rideal reactions with N atoms at Ru(0001): formation of NO and N(2)

Forward-directed NO molecules with large translational energies are formed upon exposure of an O-covered Ru(0001) surface to a nitrogen (N+N_{2}) beam. This is an unequivocal experimental demonstration of the Eley-Rideal reaction for a "heavy" (i.e., nonhydrogenated) neutral system. The time dependence of prompt NO formation exhibits an exceptionally fast decay as a consequence of shifting reaction pathways and probabilities over the course of the exposure. Prompt production shuts down as the O coverage decreases due to competition from more favorable Eley-Rideal production of N_{2}.

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.