Communication: Fingerprints of reaction mechanisms in product distributions: Eley-Rideal-type reactions between D and CD3/Cu(111)

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.

Adsorption of Atomic Hydrogen on Hydrogen Boride Sheets Studied by Photoelectron Spectroscopy

Hydrogen boride (HB) sheets are emerging as a promising two-dimensional (2D) boron material, with potential applications as unique electrodes, substrates, and hydrogen storage materials. The 2D layered structure of HB was successfully synthesized using an ion-exchange method. The chemical bonding and structure of the HB sheets were investigated using Fourier Transform Infrared (FT-IR) spectroscopy and Transmission Electron Microscopy (TEM), respectively. X-ray photoelectron spectroscopy (XPS) was employed to study the chemical states and transformation of the components before and after atomic hydrogen adsorption, thereby elucidating the atomic hydrogen adsorption process on HB sheets. Our results indicate that, upon atomic hydrogen adsorption onto the HB sheets, the B-H-B bonds were broken and converted into B-H bonds. This research highlights and demonstrates the changes in chemical states and component transformations of the boron element on the HB sheets' surface before and after atomic hydrogen adsorption, thus providing a clearer understanding of the unique bonding and structural characteristics of the HB sheets.

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.

Mechanism of Ultraviolet-Induced CO Desorption from CO Ice: Role of Vibrational Relaxation Highlighted

Although UV photon-induced CO ice desorption is clearly observed in many cold regions of the Universe as well as in the laboratory, the fundamental question of the mechanisms involved at the molecular scale remains debated. In particular, the exact nature of the involved energy transfers in the indirect desorption pathway highlighted in previous experiments is not explained. Using ab initio molecular dynamics simulations, we explore a new indirect desorption mechanism in which a highly vibrationally excited CO (v=40) within an aggregate of 50 CO molecules triggers the desorption of molecules at the surface. The desorption originates first from a mutual attraction between the excited molecule and the surrounding molecule(s), followed by a cascade of energy transfers, ultimately resulting in the desorption of vibrationally cold CO (∼95% in v=0). The theoretical vibrational distribution, along with the kinetic energy one, which peaks around 25 meV for CO with low rotational levels (v=0, J<7), is in excellent agreement with the results obtained from VUV laser induced desorption (157 nm) of CO (v=0, 1) probed using REMPI.

Computational approaches to dissociative chemisorption on metals: towards chemical accuracy

We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (E_b, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate E_b can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol^-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for E_b for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH₄. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of E_b for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate E_b for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O₂, H₂O, NH₃, CO₂, and CH₃OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.

High-Dimensional Atomistic Neural Network Potential to Study the Alignment-Resolved O2 Scattering from Highly Oriented Pyrolytic Graphite

A high dimensional and accurate atomistic neural network potential energy surface (ANN-PES) that describes the interaction between one O₂ molecule and a highly oriented pyrolytic graphite (HOPG) surface has been constructed using the open-source package (aenet). The validation of the PES is performed by paying attention to static characteristics as well as by testing its performance in reproducing previous ab initio molecular dynamics simulation results. Subsequently, the ANN-PES is used to perform quasi-classical molecular dynamics calculations of the alignment-dependent scattering of O₂ from HOPG. The results are obtained for 200 meV O₂ molecules with different initial alignments impinging with a polar incidence angle with respect to the surface normal of 22.5° on a thermalized (110 and 300 K) graphite surface. The choice of these initial conditions in our simulations is made to perform comparisons to recent experimental results on this system. Our results show that the scattering of O₂ from the HOPG surface is a rather direct process, that the angular distributions are alignment dependent, and that the final translational energy of end-on molecules is around 20% lower than that of side-on molecules. Upon collision with the surface, the molecules that are initially aligned perpendicular to the surface become highly rotationally excited, whereas a very small change in the rotational state of the scattered molecules is observed for the initial parallel alignments. The latter confirms the energy transfer dependence on the stereodynamics for the present system. The results of our simulations are in overall agreement with the experimental observations regarding the shape of the angular distributions and the alignment dependence of the in-plane reflected molecules.

Insights into adsorption, diffusion, and reactions of atomic nitrogen on a highly oriented pyrolytic graphite surface

To gain insight into the nitrogen-related gas-surface reaction dynamics on carbon-based thermal protection systems of hypersonic vehicles, we have investigated the adsorption, diffusion, and reactions of atomic nitrogen, N(⁴S), on the (0001) face of graphite using periodic density functional theory with a dispersion corrected functional. The atomic nitrogen is found to bind with pristine graphite at a bridge site, with a barrier of 0.88 eV for diffusing to an adjacent bridge site. Its adsorption energy at defect sites is significantly higher, while that between graphene layers is lower. The formation of N₂ via Langmuir-Hinshelwood (LH) and Eley-Rideal (ER) mechanisms was also investigated. In the LH pathway, the recombinative desorption of N₂ proceeds via a transition state with a relatively low barrier (0.53 eV). In addition, there is a metastable surface species, which is capable of trapping the nascent N₂ at low surface temperatures as a result of the large energy disposal into the N-N vibration. The desorbed N₂ is highly excited in both of its translational and vibrational degrees of freedom. The ER reaction is direct and fast, and it also leads to translationally and internally excited N₂. Finally, the formation of CN from a defect site is calculated to be endoergic by 2.75 eV. These results are used to rationalize the results of recent molecular beam experiments.

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.

Bridging the Gap between Direct Dynamics and Globally Accurate Reactive Potential Energy Surfaces Using Neural Networks

Direct dynamics simulations become increasingly popular in studying reaction dynamics for complex systems where analytical potential energy surfaces (PESs) are unavailable. Yet, the number and/or the propagation time of trajectories are often limited by high computational costs, and numerous energies and forces generated on-the-fly become wasted after simulations. We demonstrate here an example of reusing only a very small portion of existing direct dynamics data to reconstruct a 90-dimensional globally accurate reactive PES describing the interaction of CO₂ with a movable Ni(100) surface based on a machine learning approach. In addition to reproducing previous results with much better statistics, we predict scattering probabilities of CO₂ at the state-to-state level, which is extremely demanding for direct dynamics. We propose this unified way to investigate gaseous and gas-surface reactions of medium size, initiating with hundreds of preliminary direct dynamics trajectories, followed by low-cost and high-quality simulations on full-dimensional analytical PESs.