Six-dimensional quantum dynamics study for the dissociative adsorption of HCl on Au(111) surface

Six-dimensional quantum dynamics study for the dissociative chemisorption of H2 on pure and alloyed AgAu surfaces

The 6D time-dependent wave packet calculations were performed to explore H2 dissociation on Ag, Au, and two AgAu alloy surfaces, using four newly fitted potential energy surfaces based on the neural network fitting to density functional theory energy points. The ligand effect resulting from the Ag-Au interaction causes a reduction in the barrier height for H2+Ag/Au(111) compared to H2+Ag(111). However, the scenario is reversed for H2+Au/Ag(111) and H2+Au(111). The 6D dissociation probabilities of H2 on Ag/Au(111) surfaces are significantly higher than those on the pure Ag(111) surface, but the corresponding results for H2 on Au/Ag(111) surfaces are substantially lower than those on the pure Au(111) surface. The reactivity of H2 on Au(111) is larger than that on Ag(111), despite Ag(111) having a slightly lower static barrier height. This can be attributed to the exceptionally small dissociation probabilities at the hcp and fcc regions, which are at least 100 times smaller compared to those at the bridge or top site for H2+Ag(111). Due to the late barrier being more pronounced, the vibrational excitation of H2 on Ag(111) is more effective in promoting the reaction than on Au(111). Moreover, a high degree of alignment dependence is detected for the four reactions, where the H2 dissociation has the highest probability at the helicopter alignment, as opposed to the cartwheel alignment.

Quantum dynamics reveal different ligand effects by vibrational excitation in the dissociative chemisorption of HCl on the Au/Ag(111) surface

The reactivity and selectivity of bimetallic surfaces are of fundamental importance in industrial applications. Here, we report the first six-dimensional (6D) quantum dynamics study for the role of surface strain and ligand effects on the reactivity of HCl on a strained pseudomorphic monolayer of Au deposited onto a Ag(111) substrate, with the aid of accurate machine learning-based potential energy surfaces. The substitute of Au into Ag changes the location of the transition state; however, the static barrier height remains roughly the same as pure Au(111). The 6D quantum dynamics calculations reveal that the surface strain due to lattice expansion slightly enhances the reactivity. The ligand effect due to electronic structure interactions between Au and Ag substantially suppresses the reactivity of HCl in the ground vibrational state but promotes the reactivity via vibrational excitation at high kinetic energies. This finding can be attributed to more close interaction with Ag atoms at the transition state close to the fcc site, as well as the tight transition-state region, making the vibrational excitation highly efficient in enhancing the reactivity. Our study quantitatively unravels the dynamical origin of reactivity control by two metals, which will ultimately provide valuable insight into the selectivity of the catalyst.

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.

Effect of surface temperature on quantum dynamics of H2 on Cu(111) using a chemically accurate potential energy surface

The effect of surface atom vibrations on H₂ scattering from a Cu(111) surface at different temperatures is being investigated for hydrogen molecules in their rovibrational ground state (v = 0, j = 0). We assume weakly correlated interactions between molecular degrees of freedom and surface modes through a Hartree product type wavefunction. While constructing the six-dimensional effective Hamiltonian, we employ (a) a chemically accurate potential energy surface according to the static corrugation model [M. Wijzenbroek and M. F. Somers, J. Chem. Phys. 137, 054703 (2012)]; (b) normal mode frequencies and displacement vectors calculated with different surface atom interaction potentials within a cluster approximation; and (c) initial state distributions for the vibrational modes according to Bose-Einstein probability factors. We carry out 6D quantum dynamics with the so-constructed effective Hamiltonian and analyze sticking and state-to-state scattering probabilities. The surface atom vibrations affect the chemisorption dynamics. The results show physically meaningful trends for both reaction and scattering probabilities compared to experimental and other theoretical results.

Large-Scale Atomic Simulation via Machine Learning Potentials Constructed by Global Potential Energy Surface Exploration

Atomic simulations based on quantum mechanics (QM) calculations have entered into the tool box of chemists over the past few decades, facilitating an understanding of a wide range of chemistry problems, from structure characterization to reactivity determination. Due to the poor scaling and high computational cost intrinsic to QM calculations, one has to either sacrifice accuracy or time when performing large-scale atomic simulations. The battle to find a better compromise between accuracy and speed has been central to the development of new theoretical methods.The recent advances of machine-learning (ML)-based large-scale atomic simulations has shown great promise to the benefit of many branches of chemistry. Instead of solving the Schrödinger equation directly, ML-based simulations rely on a large data set of accurate potential energy surfaces (PESs) and complex numerical models to predict the total energy. These simulations feature both a high speed and a high accuracy for computing large systems. Due to the lack of a physical foundation in numerical models, ML models are often frustrated in their predictivity and robustness, which are key to applications. Focusing on these concerns, here we overview the recent advances in ML methodologies for atomic simulations on three key aspects. Namely, the generation of a representative data set, the extensity of ML models, and the continuity of data representation. While global optimization methods are the natural choice for building a representative data set, the stochastic surface walking method is shown to provide the desired PES sampling for both minima and transition regions on the PES. The current ML models generally utilize local geometrical descriptors as an input and consider the total energy as the sum of atomic energies. There are many flavors of data descriptors and ML models, but the applications for material and reaction predictions are still limited, not least because of the difficulty to train the associated vast global data sets. We show that our recently designed power-type structure descriptors together with a feed-forward neural network (NN) model are compatible with highly complex global PES data, which has led to a large family of global NN (G-NN) potentials.Two recent applications of G-NN potentials in material and reaction simulations are selected to illustrate how ML-based atomic simulations can help the discovery of new materials and reactions.

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.

Stochastic wave packet approach to nonadiabatic scattering of diatomic molecules from metals

In this contribution, we present a quantum dynamical approach to study inelastic scattering of diatomic molecules from metal surfaces at normal incidence. The dissipative dynamics obeys a stochastic Schrödinger equation describing the time-evolution of the system as a piecewise deterministic process. Energy exchange between the molecular vibrational degrees of freedom and the metal electrons is represented using operators in tensor product form, which are coupled via anharmonic transition rates calculated from first-order perturbation theory. Full dimensional observables are obtained by averaging over simulations in 4D-including the internal stretch, the distance to the surface, and the orientation angles-at different surface sites. The method is applied to the state-resolved scattering of vibrationally excited NO from Au(111), revealing important channels for quantized energy relaxation.

Reactive and Nonreactive Scattering of HCl from Au(111): An Ab Initio Molecular Dynamics Study

The HCl + Au(111) system has recently become a benchmark for highly activated dissociative chemisorption, which presumably is strongly affected by electron-hole pair excitation. Previous dynamics calculations, which were based on density functional theory at the generalized gradient approximation level (GGA-DFT) for the molecule-surface interaction, have all overestimated measured reaction probabilities by at least an order of magnitude. Here, we perform ab initio molecular dynamics (AIMD) and AIMD with electronic friction (AIMDEF) calculations employing a density functional that includes the attractive van der Waals interaction. Our calculations model the simultaneous and possibly synergistic effects of surface temperature, surface atom motion, electron-hole pair excitation, the molecular beam conditions of the experiments, and the van der Waals interaction on the reactivity. We find that reaction probabilities computed with AIMDEF and the SRP32-vdW functional still overestimate the measured reaction probabilities, by a factor 18 for the highest incidence energy at which measurements were performed (≈2.5 eV). Even granting that the experiment could have underestimated the sticking probability by about a factor three, this still translates into a considerable overestimation of the reactivity by the current theory. Likewise, scaled transition probabilities for vibrational excitation from ν = 1, j = 1 to ν = 2 are overestimated by the AIMDEF theory, by factors 3-8 depending on the initial conditions modeled. Energy losses to the surface and translational energy losses are, however, in good agreement with experimental values.

Water dissociating on rigid Ni(100): A quantum dynamics study on a full-dimensional potential energy surface

We constructed a nine-dimensional (9D) potential energy surface (PES) for the dissociative chemisorption of H₂O on a rigid Ni(100) surface using the neural network method based on roughly 110 000 energies obtained from extensive density functional theory (DFT) calculations. The resulting PES is accurate and smooth, based on the small fitting errors and the good agreement between the fitted PES and the direct DFT calculations. Time dependent wave packet calculations also showed that the PES is very well converged with respect to the fitting procedure. The dissociation probabilities of H₂O initially in the ground rovibrational state from 9D quantum dynamics calculations are quite different from the site-specific results from the seven-dimensional (7D) calculations, indicating the importance of full-dimensional quantum dynamics to quantitatively characterize this gas-surface reaction. It is found that the validity of the site-averaging approximation with exact potential holds well, where the site-averaging dissociation probability over 15 fixed impact sites obtained from 7D quantum dynamics calculations can accurately approximate the 9D dissociation probability for H₂O in the ground rovibrational state.

High-Dimensional Atomistic Neural Network Potentials for Molecule-Surface Interactions: HCl Scattering from Au(111)

Ab initio molecular dynamics (AIMD) simulations of molecule-surface scattering allow first-principles characterization of the dynamics. However, the large number of density functional theory calculations along the trajectories is very costly, limiting simulations of long-time events and giving rise to poor statistics. To avoid this computational bottleneck, we report here the development of a high-dimensional molecule-surface interaction potential energy surface (PES) with movable surface atoms, using a machine learning approach. With 60 degrees of freedom, this PES allows energy transfer between the energetic impinging molecule and thermal surface atoms. Classical trajectory calculations for the scattering of DCl from Au(111) on this PES are found to agree well with AIMD simulations, with ∼10⁵-fold acceleration. Scattering of HCl from Au(111) is further investigated and compared with available experimental results.

Dissociative chemisorption of methane on Ni(111) using a chemically accurate fifteen dimensional potential energy surface

A new chemically accurate potential energy surface for the dissociative chemisorption of methane on the rigid Ni(111) surface.

An approximate full-dimensional quantum dynamics study of the mode specificity in the dissociative chemisorption of D2O on rigid Cu(111)

The approximate 9D dissociation probabilities for D₂O/Cu(111) are obtained to investigate the influence of mode specificity on reaction dynamics.

Effects of Lattice Motion on Dissociative Chemisorption: Toward a Rigorous Comparison of Theory with Molecular Beam Experiments

The dissociative chemisorption of small molecules such as methane and water on metal surfaces is a key step in many important catalyzed reactions. However, it has only very recently become possible to directly compare theory with molecular beam studies of these reactions. For most experimental conditions, such a comparison requires accurate methods for introducing the effects of lattice motion into quantum reactive scattering calculations. We examine these methods and their recent application to methane and water dissociative chemisorption. New results are presented for CO₂ chemisorption and methane dissociation at step edges. The type of molecule-lattice coupling that leads to a strong variation in the dissociative sticking of methane with temperature is shown to occur for many polyatomic-metal systems. Improvements to these models are discussed. The ability to accurately compare theory with molecular beam experiments should lead to improved density functionals and consequently more accurate thermal rate constants for these important reactions.

First-principles quantum dynamical theory for the dissociative chemisorption of H2O on rigid Cu(111)

Despite significant progress made in the past decades, it remains extremely challenging to investigate the dissociative chemisorption dynamics of molecular species on surfaces at a full-dimensional quantum mechanical level, in particular for polyatomic-surface reactions. Here we report, to the best of our knowledge, the first full-dimensional quantum dynamics study for the dissociative chemisorption of H2O on rigid Cu(111) with all the nine molecular degrees of freedom fully coupled, based on an accurate full-dimensional potential energy surface. The full-dimensional quantum mechanical reactivity provides the dynamics features with the highest accuracy, revealing that the excitations in vibrational modes of H2O are more efficacious than increasing the translational energy in promoting the reaction. The enhancement of the excitation in asymmetric stretch is the largest, but that of symmetric stretch becomes comparable at very low energies. The full-dimensional characterization also allows the investigation of the validity of previous reduced-dimensional and approximate dynamical models.

Activated Dissociation of HCl on Au(111)

We report zero-coverage reaction probabilities (S0) for HCl dissociative adsorption on Au(111) obtained by the seeded molecular beam hot-nozzle method. For measurements at normal incidence with mean translational energies ranging from 0.94 to 2.56 eV (nozzle temperatures 296 to 1060 K), S0 increased from 6 × 10(-6) to 2 × 10(-2). S0 also increased with increasing nozzle temperature for fixed incidence energy associated with the motion normal to the surface. Accounting for the influence of the vibrational state population and translational energy distributions in the incident beam, we are able to compare the experimental results to recent theoretical predictions. These calculations, performed employing 6-D quantum dynamics on an electronically adiabatic potential energy surface obtained using density functional theory at the level of the generalized gradient approximation and the static surface approximation, severely overestimate the reaction probabilities when compared with our experimental results. We discuss some possible reasons for this large disagreement.

A seven-dimensional quantum dynamics study of the dissociative chemisorption of H2O on Cu(111): effects of azimuthal angles and azimuthal angle-averaging

We report the first seven-dimensional quantum dynamics study for the dissociative chemisorption of H2O on Cu(111) using the time-dependent wave-packet approach, based on an accurate nine-dimensional potential energy surface (PES), which is newly developed by neural network fitting to ∼80 000 density functional theory points. This seven-dimensional quantum model allows the examination of the influence of azimuthal angles and also the investigation of the quantitative relationship between the seven-dimensional (7D) dissociation probabilities and those results calculated by the six-dimensional (6D) model with the flat surface approximation. The reactivity strongly depends on the azimuthal rotations due to different barrier heights. Very large differences are seen between the 7D dissociation probabilities and the 6D results with fixed azimuthal angles, at different fixed sites of impact, indicating that the 6D model by neglecting the azimuthal rotation can introduce substantial errors in calculating dissociation probabilities and the 7D quantum dynamics is essential to investigate the dissociation process. A new azimuthal angle-averaging approach is proposed that the 7D dissociation probability can be well reproduced by averaging 6D results over 18 azimuthal angles, in particular in low kinetic energy regions.

Methane dissociation on Ni(111): A fifteen-dimensional potential energy surface using neural network method

In the present work, we develop a highly accurate, fifteen-dimensional potential energy surface (PES) of CH4 interacting on a rigid flat Ni(111) surface with the methodology of neural network (NN) fit to a database consisted of about 194 208 ab initio density functional theory (DFT) energy points. Some careful tests of the accuracy of the fitting PES are given through the descriptions of the fitting quality, vibrational spectrum of CH4 in vacuum, transition state (TS) geometries as well as the activation barriers. Using a 25-60-60-1 NN structure, we obtain one of the best PESs with the least root mean square errors: 10.11 meV for the entrance region and 17.00 meV for the interaction and product regions. Our PES can reproduce the DFT results very well in particular for the important TS structures. Furthermore, we present the sticking probability S0 of ground state CH4 at the experimental surface temperature using some sudden approximations by Jackson's group. An in-depth explanation is given for the underestimated sticking probability.

Mode specificity for the dissociative chemisorption of H2O on Cu(111): a quantum dynamics study on an accurately fitted potential energy surface

The mode-specific dynamics for the dissociative chemisorption of H₂O on Cu(111) is first investigated by seven-dimensional quantum dynamics calculations, based on an accurately fitted potential energy surface (PES) recently developed by neural network fitting to DFT energy points.

Mode specificity of the dissociative chemisorption of HOD on rigid Cu(111): an approximate full-dimensional quantum dynamics study

The approximate 9D dissociation probabilities for HOD/Cu(111) are obtained to investigate the influence of the mode specificity on reaction dynamics.

Electronically non-adiabatic influences in surface chemistry and dynamics

Electronically nonadiabatic interactions between molecules and metal surfaces are now well known. But evidence that such interactions influence reaction rates is still scarce. This paper reviews research related to this topic and proposes pathways forward.

Towards an accurate specific reaction parameter density functional for water dissociation on Ni(111): RPBE versus PW91

Approximated nine dimensional quantum dynamics on a new potential energy surface for water dissociation on Ni(111) computed using the RPBE functional.

The direct and precursor mediated dissociation rates of H2 on a Ni(111) surface

The dissociation and recombination rates of physisorbed H2, and the direct and steady state dissociation (i.e., the precursor mediated dissociation) rates of gas phase H2 on Ni(111), as well as the corresponding kinetic isotope effects, are calculated using the quantum instanton method, together with path integral Monte Carlo and adaptive umbrella sampling techniques. All these rates except the recombination one first decrease and then increase with the increasing temperature, and their minimum values appear at about 250, 300 and 250 K, respectively. These non-monotonic behaviors reveal that the quantum effect of H2 should be very remarkable at low temperatures. The steady state rates are smaller than the direct rates at low temperatures, however, they become larger than the direct ones at high temperatures, these two kinds of rates become equal at about 400 and 300 K on the rigid and quantum lattices, respectively. The quantum motion of the lattice can enhance the direct and steady state rates, and it increases the steady state rate much more than the direct one, for instance, the direct and steady state rates on the quantum lattice are 1.30 and 2.08 times larger than that on the rigid one at 300 K. The calculated kinetic isotope effects are much larger than 1, which reveals that H2 always has a larger rate than that of D2, and the direct process predicts much larger kinetic isotope effects than the steady state process at low temperatures. In addition, the kinetic isotope effects are not affected by the lattice motion.

Six-dimensional quantum dynamics for dissociative chemisorption of H2 and D2 on Ag(111) on a permutation invariant potential energy surface

Quantum dynamics on a permutation invariant potential energy surface for H₂ dissociation on Ag(111) yield satisfactory agreement with experiment.