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

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.