Reactive force fields for surface chemical reactions: A case study with hydrogen dissociation on Pd surfaces

The quantum dynamics of H2 on Cu(111) at a surface temperature of 925 K: Comparing state-of-the-art theory to state-of-the-art experiments 2

State-of-the-art 6D quantum dynamics simulations for the dissociative chemisorption of H₂ on a thermally distorted Cu(111) surface, using the static corrugation model, were analyzed to produce several (experimentally available) observables. The expected error, especially important for lower reaction probabilities, was quantified using wavepackets on several different grids as well as two different analysis approaches to obtain more accurate results in the region where a slow reaction channel was experimentally shown to be dominant. The lowest reaction barrier sites for different thermally distorted surface slabs are shown to not just be energetically, but also geometrically, different between surface configurations, which can be used to explain several dynamical effects found when including surface temperature effects. Direct comparison of simulated time-of-flight spectra to those obtained from state-of-the-art desorption experiments showed much improved agreement compared to the perfect lattice BOSS approach. Agreement with experimental rotational and vibrational efficacies also somewhat improved when thermally excited surfaces were included in the theoretical model. Finally, we present clear quantum effects in the rotational quadrupole alignment parameters found for the lower rotationally excited states, which underlines the importance of careful quantum dynamical analyses of this system.

Accurate Description of the Quantum Dynamical Surface Temperature Effects on the Dissociative Chemisorption of H2 from Cu(111)

Accurately describing surface temperature effects for the dissociative scattering of H2 on a metal surface on a quantum dynamical level is currently one of the open challenges for theoretical surface scientists. We present the first quantum dynamical (QD) simulations of hydrogen dissociating on a Cu(111) surface which accurately describe all relevant surface temperature effects, using the static corrugation model (SCM). The reaction probabilities we obtain show very good agreement with those found using quasi-classical dynamics (QCD), both for individual surface slabs and for an averaged, thus Monte-Carlo sampled, set of thermally distorted surface configurations. Rovibrationally elastic scattering probabilities show a much clearer difference between the QCD and QD results, which appears to be traceable back towards thermally distorted surface configurations with very low dissociation probabilities and underlines the importance of investigating more observables than just dissociation. By reducing the number of distorted surface atoms included in the dynamical model, we also show that only including one, or even three, surface atoms is generally not enough to accurately describe the effects of surface temperature on dissociation and elastic scattering. These results are a major step forward in accurately describing hydrogen scattering from a thermally excited Cu(111) surface, and open up a pathway to better describe reaction and scattering from other relevant crystal facets, such as stepped surfaces, at moderately elevated surface temperatures where quantum effects are expected to play a more important role in the dissociation of H2 on Cu.

First-Principles Insights into Adiabatic and Nonadiabatic Vibrational Energy-Transfer Dynamics during Molecular Scattering from Metal Surfaces: The Importance of Surface Reactivity

Energy transfer is ubiquitous during molecular collisions and reactions at gas-surface interfaces. Of particular importance is vibrational energy transfer because of its relevance to bond forming and breaking. In this Perspective, we review recent first-principles studies on vibrational energy-transfer dynamics during molecular scattering from metal surfaces at the state-to-state level. Taking several representative systems as examples, we highlight the intrinsic correlation between vibrational energy transfer in nonreactive scattering and surface reactivity and how it operates in both electronically adiabatic and nonadiabatic pathways. Adiabatically, the presence of a dissociation barrier softens a bond in the impinging molecule and increases its couplings with other molecular modes and surface phonons. In the meantime, the stronger interaction between the molecule and the surface also changes the electronic structure at the barrier, resulting in an increase of nonadiabatic effects. We further discuss future prospects toward a more quantitative understanding of this important surface dynamical process.

Revealing the thermal decomposition mechanism of RDX crystals by a neural network potential

A neural network potential (NNP) is developed to investigate the complex reaction dynamics of 1,3,5-trinitro-1,3,5-triazine (RDX) thermal decomposition.

Accurate Simulations of the Reaction of H2 on a Curved Pt Crystal through Machine Learning

Theoretical studies on molecule-metal surface reactions have so far been limited to small surface unit cells due to computational costs. Here, for the first time molecular dynamics simulations on very large surface unit cells at the level of density functional theory are performed, allowing a direct comparison to experiments performed on a curved crystal. Specifically, the reaction of D2 on a curved Pt crystal is investigated with a neural network potential (NNP). The developed NNP is also accurate for surface unit cells considerably larger than those that have been included in the training data, allowing dynamical simulations on very large surface unit cells that otherwise would have been intractable. Important and complex aspects of the reaction mechanism are discovered such as diffusion and a shadow effect of the step. Furthermore, conclusions from simulations on smaller surface unit cells cannot always be transfered to larger surface unit cells, limiting the applicability of theoretical studies of smaller surface unit cells to heterogeneous catalysts with small defect densities.

Highly Efficient Activation of HCl Dissociation on Au(111) via Rotational Preexcitation

The probability for dissociation of molecules on metal surfaces, which often controls the rate of industrially important catalytic processes, can depend strongly on how energy is partitioned in the incident molecule. There are many example systems where the addition of vibrational energy promotes reaction more effectively than the addition of translational energy, but for rotational pre-excitation similar examples have not yet been discovered. Here, we make an experimentally testable theoretical prediction that adding energy to the rotation of HCl can promote its dissociation on Au(111) 20 times more effectively than increasing its translational energy. In the underlying mechanism, the molecule's initial rotational motion allows it to pass through a critical region of the reaction path, where this path shows a strong and nonmonotonic dependence on the molecular orientation.

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.

Beyond the static corrugation model: Dynamic surfaces with the embedded atom method

The D₂ on Cu(111) system has for many years been one of the major benchmark systems for surface scientists. Generating surface configurations using the embedded atom method (EAM), we investigate the quality of the chemically accurate static corrugation model (SCM) for including surface temperature effects, with a focus on the random displacement approach to its distorted surface generation. With this EAM potential, we also treat the Cu(111) surface of our system fully dynamically and shine a further light on not only the quality of the SCM sudden approach but also the limited effect of energy exchange with the surface. Reaction and (in)elastic scattering probability curves, as well as simulated time-of-flight spectra, show good agreement with both earlier works and experimental results, with surface reactions showing a preference for surface atoms displaced away from the incoming molecule. The good agreement with the non-static surface model also further establishes the limited effect of energy exchange on not only the reaction but also on the elastic and inelastic scattering probabilities, even though some molecular translational energy is deposited into the surface.

Nonthermalized Precursor-Mediated Dissociative Chemisorption at High Catalysis Temperatures

Quasiclassical trajectory calculations and vibrational-state-selected beam-surface measurements of CH₄ chemisorption on Ir(111) reveal a nonthermal, hot-molecule mechanism for C-H bond activation. Low-energy vibrationally excited molecules become trapped in the physisorption well and react before vibrational and translational energies accommodate the surface. The reaction probability is strongly surface-temperature-dependent and arises from the pivotal role of Ir atom thermal motion. In reactive trajectories, the mean outward Ir atom displacement largely exceeds that of the transition-state geometry obtained through a full geometry optimization. The study also highlights a new way for (temporary) surface defects to impact high-temperature heterogeneous catalytic reactivity. Instead of reactants diffusing to and competing for geometrically localized lower barrier sites, transient, thermally activated surface atom displacements deliver low-barrier surface reaction geometries to the physisorbed reactants.

Unified and transferable description of dynamics of H2 dissociative adsorption on multiple copper surfaces via machine learning

Schematic of the developed neural network potential energy surface enabling a unified and transferable description of dynamics of H₂ dissociative adsorption on multiple copper surfaces.

Simple Molecular Reactive Force Field for Metal-Organic Synthesis

For colloidal quantum dots to transition from research laboratories to deployment as optical and electronic products, there will be a need to scale-up their production to large-scale manufacturing processes. This demand increases the need to understand their formation via a molecular representation of the nucleation of lead sulfide (PbS) quantum dot systems passivated by lead oleate complexes. We demonstrate the effectiveness of a new type of reactive potential, custom-made for this system, that is drawn from simple Morse, Lennard-Jones, and Coulombic components, which can reproduce reactions across a broad range of PbS quantum dot sizes with good accuracy. We validate the capability of this model to capture reactive systems by comparison to ab initio calculations for a reaction between two dots.

Quantum and classical dynamics of reactive scattering of H2 from metal surfaces

State-of-the-art theoretical models allow nowadays an accurate description of H₂/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.

Diffusion Rates for Hydrogen on Pd(111) from Molecular Quantum Dynamics Calculations

The van Hove formula for the dynamical structure factor (DSF) related to particle scattering at mobile adsorbates is extended to include the relaxation of the adsorbates' vibrational states. The total rate obtained from the DSF is assumed to be the sum of a diffusion and a relaxation rate. A simple kinetic model to support this assumption is presented. To illustrate its potential applicability, the formula is evaluated using wave functions, energies, and lifetimes of vibrational states obtained for H/Pd(111) from first-principle calculations. Results show that quantum effects can be expected to be important even at room temperature.

Towards bond selective chemistry from first principles: methane on metal surfaces

Controlling bond-selective chemical reactivity is of great importance and has a broad range of applications. Here, we present a molecular dynamics study of bond selective reactivity of methane and its deuterated isotopologues (i.e., CH(4-x)D(x), x=0,1,2,3,4) on Ni(111) and Pt(111) from first principles calculations. Our simulations allow for reproducing the full C-H bond selectivity recently achieved experimentally via mode-specific vibrational excitation and explain its origin. Moreover, we also predict the hitherto unexplored influence of the molecular translational energy on such a selectivity as well as the conditions under which the full selectivity can be realized for the a priori less active C-D bond.

Towards chemically accurate simulation of molecule-surface reactions

This perspective addresses four challenges facing theorists whose aim is to make quantitatively accurate predictions for reactions of molecules on metal surfaces, and suggests ways of meeting these challenges, focusing on dissociative chemisorption reactions of H(2), N(2), and CH(4). Addressing these challenges is ultimately of practical importance to a more accurate description of overall heterogeneously catalysed reactions, which play a role in the production of more than 90% of man-made chemicals. One challenge is to describe the interaction of a molecule with a metal surface with chemical accuracy, i.e., with errors in reaction barrier heights less than 1 kcal mol(-1). In this framework, the potential of a new implementation of specific reaction parameter density functional theory (SRP-DFT) will be discussed, with emphasis on applications to reaction of H(2) with metal surfaces. Two additional challenges are to come up with improved descriptions of the effects of phonons and electron-hole pairs on reaction of molecules like N(2) on metal surfaces. Phonons can be tackled using sudden approximations in quantum dynamics, and through Ab Initio Molecular Dynamics (AIMD) calculations using classical dynamics. To additionally achieve an accurate description of the effect of electron-hole pair excitation on dissociative chemisorption within a classical dynamics framework, it may be possible to combine AIMD with electronic friction. The fourth challenge we will consider is how to achieve an accurate quantum mechanical description of the dissociative chemisorption of a polyatomic molecule, like methane, on a metal surface. A method of potential interest is the Multi-Configuration Time-Dependent Hartree (MCTDH) method.