Mode and Bond Selectivities in Methane Dissociative Chemisorption: Quasi-Classical Trajectory Studies on Twelve-Dimensional Potential Energy Surface

Advances in the Dynamics of Adsorbate Diffusion on Metal Surfaces: Focus on Hydrogen and Oxygen

Adsorbates on metal surfaces are typically formed from the dissociative chemisorption of molecules occurring at gas-solid interfaces. These adsorbed species exhibit unique diffusion behaviors on metal surfaces, which are influenced by their translational energy. They play crucial roles in various fields, including heterogeneous catalysis and corrosion. This review examines recent theoretical advancements in understanding the diffusion dynamics of adsorbates on metal surfaces, with a specific emphasis on hydrogen and oxygen atoms. The diffusion processes of adsorbates on metal surfaces involve two energy transfer mechanisms: surface phonons and electron-hole pair excitations. This review also surveys new theoretical methods, including the characterization of the electron-hole pair excitation within electronic friction models, the acceleration of quantum chemistry calculations through machine learning, and the treatment of atomic nuclear motion from both quantum mechanical and classical perspectives. Furthermore, this review offers valuable insights into how energy transfer, nuclear quantum effects, supercell sizes, and the topography of potential energy surfaces impact the diffusion behavior of hydrogen and oxygen species on metal surfaces. Lastly, some preliminary research proposals are presented.

Mode Selective Chemistry for the Dissociation of Methane on Efficient Ni/Pt-Bimetallic Alloy Catalysts

The mode selectivity of methane dissociation is studied on three different Ni/Pt-bimetallic alloy surfaces using a fully quantum approach based on reaction path Hamiltonian. Dissociative sticking probability depends on the...

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.

Reagent-Triggered Isomerization of Fluxional Cluster Catalyst via Dynamic Coupling

Metallic cluster catalysts have many thermodynamically accessible isomers with diverse active sites and low reaction barriers, and lately a strong hypothesis emerged that the many catalyst states collectively drive the catalysis. However, it remained a hypothesis that catalyst isomerization is actually kinetically feasible under the current reaction conditions. Using high-temperature dynamics simulations and sampling, a range of orientations, and vibrational energy distributions, we probe how thermal effects and molecular events affect cluster catalyst dynamics. We show that even such a delicate affair as the dissociation or scattering of a methane molecule on the heavy and thus slow Pt₁₃ cluster triggers substantial isomerization of the catalyst, far beyond thermal at 700 K. A kinetic coupling between the methane activity and cluster catalyst dynamics is observed. In return, the thermal dynamics of the cluster affects the methane reaction and scattering probabilities. Hence, molecular events at the surfaces of fluxional cluster catalysts should facilitate the population of an ensemble of catalyst states under the current reaction conditions, with implications for available active sites, reaction mechanisms, and apparent rates.

Plasma-enhanced catalysis for the upgrading of methane: a review of modelling and simulation methods

Modelling methods and simulation works on the upgrading of methane via plasma and plasma-enhanced catalysis reviewed.

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.

Mode-specific quantum dynamics and kinetics of the hydrogen abstraction reaction OH + H2O → H2O + OH

The synergistic effect between the reactant stretching and bending modes on promoting the reaction.

Vibrational control of selective bond cleavage in dissociative chemisorption of methanol on Cu(111)

Controlling product branching ratios in a chemical reaction represents a desired but difficult achievement in chemistry. In this work, we demonstrate the first example of altering the branching ratios in a multichannel reaction, i.e., methanol dissociative chemisorption on Cu(111), via selectively exciting specific vibrational modes. To this end, we develop a globally accurate full-dimensional potential energy surface for the CH3OH/Cu(111) system and perform extensive vibrational state-selected molecular dynamics simulations. Our results show that O-H/C-H/C-O stretching vibrational excitations substantially enhance the respective bond scission processes, representing extraordinary bond selectivity. At a given total energy, the branching ratio of C-O/C-H dissociation can increase by as large as 100 times by exciting the C-O stretching mode which possesses an unprecedentedly strong vibrational efficacy on reactivity. This vibrational control can be realized by the well-designed experiment using a linearly polarized laser.

Rotational and steric effects in water dissociative chemisorption on Ni(111)

Powerful laser techniques have recently enabled quantum-state resolved molecular beam experiments for investigating gas-surface reactions, which have unveiled intriguing vibrational, rotational, and also steric effects. For reactions involving polyatomic molecules, e.g., the dissociative chemisorption of methane and water, the rotational and related steric effects are far less understood despite a large body of theoretical work having been able to reproduce the observed vibrational mode specificity and related bond selectivity semi-quantitatively or even within chemical accuracy. Herein, we report a high dimensional quantum dynamics study of water dissociation on Ni(111) on a first-principles potential energy surface, focusing on the reactivities of D2O in various rotational quantum states with different spatial orientations. Through an accurate quantum mechanical description of this asymmetric top, remarkable dependence of the reactivity on the orientation is observed. This dependence is site specific and rotational state specific. These single site rotational and steric effects are partially justified by a sudden model on the basis of the overlap between the rotational wavefunctions and the angular potential near the transition state, but rotational steering also plays a significant role which complicates the dynamics. Although site averaging weakens the influence of initial rotational excitations and leads to minor effects to the reactivity, steric effects are predicted to be observable if the water molecule is selectively excited and aligned by a linearly polarized laser.

Dynamics of carbon monoxide dissociation on Co(112̄0)

The dissociative chemisorption dynamics of CO on rigid Co(112̄0) is investigated using a quasi-classical trajectory method on a new global six-dimensional potential energy surface.

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.

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.

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.

Quantum dynamics of polyatomic dissociative chemisorption on transition metal surfaces: mode specificity and bond selectivity

Recent advances in quantum dynamical characterization of polyatomic dissociative chemisorption on accurate global potential energy surfaces are critically reviewed.

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.

Reaction Rate Constants of CH4(ads) ⇌ CH3(ads) + H(ads) on Ni(111): The Effect of Lattice Motion

Methane dissociation on metal surfaces is of great commercial importance. The dissociation and recombination rate constants of CH4 on Ni(111) are calculated using the quantum instanton approach with the path integral Monte Carlo method. The Ni(111) lattice is treated rigidly, classically, and quantum mechanically to reveal the effects of lattice motion and quantum tunneling. For the dissociation of CH4, the rates have the smallest value on the rigid lattice, while they possess the largest value on the quantum lattice. For instance, at 300 K, the rates on the classical and quantum lattices are 5 and 12 times larger than that on the rigid lattice, respectively. The curve of the Arrhenius plot for the dissociation rates on the rigid lattice demonstrates that the quantum tunneling effect of the ruptured H atom is remarkable, while the nearly invariable dissociation rates at low temperatures on the quantum lattice confirm that the thermally assisted tunneling should be dominant at low temperatures. For the recombination of CH4, the quantum lattice still has rates that are much larger than that of the rigid lattice. For instance, the ratio of the recombination rates on the quantum and rigid lattices is 12 at 300 K. The quantum tunneling effect seems to play a minor role in the recombination rates on the rigid lattice; however, the thermally assisted tunneling is still very significant for the recombination process.

Dynamics of water dissociative chemisorption on Ni(111): effects of impact sites and incident angles

The dissociative chemisorption of water on rigid Ni(111) is investigated using a quasiclassical trajectory method on a nine-dimensional global potential energy surface based on a faithful permutation invariant fit of ∼25 000 density functional theory points. This full-dimensional model not only confirms the validity of our earlier reduced-dimensional model with 6 degrees of freedom, but also allows the examination of the influence of impact sites and incident angles. It is shown that the reactivity depends on the site of impact in a complex fashion controlled by the topography of the potential energy surface rather than the barrier height alone. In addition, the reaction is promoted by momenta both parallel and perpendicular to the surface, as predicted by the recently proposed sudden vector projection model.

Six-dimensional quantum dynamics of dissociative chemisorption of H2 on Co(0001) on an accurate global potential energy surface

Six-dimensional quantum dynamics of hydrogen dissociative chemisorption on Co(0001) is investigated on a DFT based potential energy surface.

The sudden vector projection model for reactivity: mode specificity and bond selectivity made simple

CONSPECTUS: Mode specificity is defined by the differences in reactivity due to excitations in various reactant modes, while bond selectivity refers to selective bond breaking in a reaction. These phenomena not only shed light on reaction dynamics but also open the door for laser control of reactions. The existence of mode specificity and bond selectivity in a reaction indicates that not all forms of energy are equivalent in promoting the reactivity, thus defying a statistical treatment. They also allow the enhancement of reactivity and control product branching ratio. As a result, they are of central importance in chemistry. This Account discusses recent advances in our understanding of these nonstatistical phenomena. In particular, the newly proposed sudden vector projection (SVP) model and its applications are reviewed. The SVP model is based on the premise that the collision in many direct reactions is much faster than intramolecular vibrational energy redistribution in the reactants. In such a sudden limit, the coupling of a reactant mode with the reaction coordinate at the transition state, which dictates its ability to promote the reaction, is approximately quantified by the projection of the former onto the latter. The SVP model can be considered as a generalization of the venerable Polanyi's rules, which are based on the location of the barrier. The SVP model is instead based on properties of the saddle point and as a result capable of treating the translational, rotational, and multiple vibrational modes in reactions involving polyatomic reactants. In case of surface reactions, the involvement of surface atoms can also be examined. Taking advantage of microscopic reversibility, the SVP model has also been used to predict product energy disposal in reactions. This simple yet powerful rule of thumb has been successfully demonstrated in many reactions including uni- and bimolecular reactions in the gas phase and gas-surface reactions. The success of the SVP model underscores the importance of the transition state in controlling mode-specific and bond-selective chemistry.

Methane dissociative chemisorption and detailed balance on Pt(111): dynamical constraints and the modest influence of tunneling

A dynamically biased (d-) precursor mediated microcanonical trapping (PMMT) model of the activated dissociative chemisorption of methane on Pt(111) is applied to a wide range of dissociative sticking experiments, and, by detailed balance, to the methane product state distributions from the thermal associative desorption of adsorbed hydrogen with coadsorbed methyl radicals. Tunneling pathways were incorporated into the d-PMMT model to better replicate the translational energy distribution of the desorbing methane product from the laser induced thermal reaction of coadsorbed hydrogen and methyl radicals occurring near T(s) = 395 K. Although tunneling is predicted to be inconsequential to the thermal dissociative chemisorption of CH4 on Pt(111) at the high temperatures of catalytic interest, once the temperature drops to 395 K the tunneling fraction of the reactive thermal flux reaches 15%, and as temperatures drop below 275 K the tunneling fraction exceeds 50%. The d-PMMT model parameters of {E0 = 58.9 kJ/mol, s = 2, η(v) = 0.40} describe the apparent threshold energy for CH4/Pt(111) dissociative chemisorption, the number of surface oscillators involved in the precursor complex, and the efficacy of molecular vibrational energy to promote reaction, relative to translational energy directed along the surface normal. Molecular translations parallel to the surface and rotations are treated as spectator degrees of freedom. Transition state vibrational frequencies are derived from generalized gradient approximation-density functional theory electronic structure calculations. The d-PMMT model replicates the diverse range of experimental data available with good fidelity, including some new effusive molecular beam and ambient gas dissociative sticking measurements. Nevertheless, there are some indications that closer agreement between theory and experiments could be achieved if a surface efficacy less than one was introduced into the modeling as an additional dynamical constraint.

Vibrationally promoted dissociation of water on Ni(111)

Water dissociation on transition-metal catalysts is an important step in steam reforming and the water-gas shift reaction. To probe the effect of translational and vibrational activation on this important heterogeneous reaction, we performed state-resolved gas/surface reactivity measurements for the dissociative chemisorption of D2O on Ni(111), using molecular beam techniques. The reaction occurs via a direct pathway, because both the translational and vibrational energies promote the dissociation. The experimentally measured initial sticking probabilities were used to calibrate a first-principles potential energy surface based on density functional theory. Quantum dynamical calculations on the scaled potential energy surface reproduced the experimental results semiquantitatively. The larger increase of the dissociation probability by vibrational excitation than by translation per unit of energy is consistent with a late barrier along the O-D stretch reaction coordinate.

Vibrational Activation of Methane Chemisorption: The Role of Symmetry

Quantum state-resolved reactivity measurements probe the role of vibrational symmetry on the vibrational activation of the dissociative chemisorption of CH4 on Pt(111). IR-IR double resonance excitation in a molecular beam is used to prepare CH4 in all three different vibrational symmetry components A1, E, and F2 of the 2ν3 antisymmetric stretch overtone vibration. Methyl dissociation products chemisorbed on the cold Pt(111) surface are detected via reflection absorption infrared spectroscopy (RAIRS). We observe similar reactivity for CH4 prepared in the A1 and F2 sublevels but up to a factor of 2 lower reactivity for excitation of the E sublevel. It is suggested that differences in the localization of the C-H stretch amplitudes for the three states at the transition state leads to the observed difference in reactivity rather than state-specific vibrational energy transfer to electronic excitation of the metal.

Ab Initio Molecular Dynamics Calculations versus Quantum-State-Resolved Experiments on CHD3 + Pt(111): New Insights into a Prototypical Gas-Surface Reaction

The dissociative chemisorption of methane on metal surfaces is of fundamental and practical interest, being a rate-limiting step in the steam reforming process. The reaction is best modeled with quantum dynamics calculations, but these are currently not guaranteed to produce accurate results because they rely on potential energy surfaces based on untested density functionals and on untested dynamical approximations. To help overcome these limitations, here we present for the first time statistically accurate reaction probabilities obtained with ab initio molecular dynamics (AIMD) for a polyatomic gas-phase molecule reacting with a metal surface. Using a general purpose density functional, the AIMD reaction probabilities are in semiquantitative agreement with new quantum-state-resolved experiments on CHD3 + Pt(111). The comparison suggests the use of the sudden approximation for treating the rotations even though CHD3 has large rotational constants and yields an estimated reaction barrier of 0.9 eV for CH4 + Pt(111).

Mode specificity and product energy disposal in unimolecular reactions: insights from the sudden vector projection model

A simple model is proposed to predict mode specificity and product energy disposal in unimolecular dissociation reactions. This so-called Sudden Vector Projection (SVP) model quantifies the coupling of a reactant or product mode with the reaction coordinate at the transition state by projecting the corresponding normal mode vector onto the imaginary frequency mode at the saddle point. Due to the sudden assumption, SVP predictions for mode specificity are expected to be valid only when the reactant molecule has weak intermodal coupling. On the other hand, the sudden limit is generally satisfied for its predictions of product energy disposal in unimolecular reactions with a tight barrier. The SVP model is applied to several prototypical systems and the agreement with available experimental and theoretical results is satisfactory.