Dissociative Chemisorption of Methane on Ni and Pt Surfaces: Mode-Specific Chemistry and the Effects of Lattice Motion

On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K

We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.

Accurate Reaction Probabilities for Translational Energies on Both Sides of the Barrier of Dissociative Chemisorption on Metal Surfaces

Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule-metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.

SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces

Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.

Effect of temperature on CO oxidation over Pt(111) in two-dimensional confinement

Confined catalysis between a two-dimensional (2D) cover and metal surfaces has provided a unique environment with enhanced activity compared to uncovered metal surfaces. Within this 2D confinement, weakened adsorption and lowered activation energies were observed using surface science experiments and density functional theory (DFT) calculations. Computationally, the role of electronic and mechanical factors responsible for the improved activity was deduced only from static DFT calculations. This demands a detailed investigation on the dynamics of reactions under 2D confinement, including temperature effects. In this work, we study CO oxidation on a 2D graphene covered Pt(111) surface at 90 and 593 K using DFT-based ab initio molecular dynamics simulations starting from the transition state configuration. We show that CO oxidation in the presence of a graphene cover is substantially enhanced (2.3 times) at 90 K. Our findings suggest that 2D confined spaces can be used to enhance the activity of chemical reactions, especially at low temperatures.

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.

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...

Effects of alloying for steam or dry reforming of methane: a review of recent studies

A survey on the catalytic nature of Ni-based alloy catalysts in recent years provides a direction for future catalyst development.

The trapping of methane on Ir(111): A first-principles quantum study

We implement a fully quantum mechanical study of methane trapping on Ir(111), where the phonons, the molecule-surface interaction, and the molecule-phonon coupling are all computed from first-principles. We find that both the surface corrugation and the phonon coupling vary strongly with molecular orientation and that there is a "chemical" aspect to this due to the catalytic nature of the metal. For example, molecules with reactive orientations can approach close to surface sites with low barriers to dissociation. Moreover, lattice motion can lower the barrier to dissociation, leading to unusual behavior for the phonon coupling. We find good agreement with experiment and two recent classical studies if we average our potential energy surface over several orientations of the molecule. We also find reasonable agreement with a recent study of methane diffraction, although we show that diffraction does not play a major role in trapping on the smooth Ir(111) surface and that trapping obeys normal energy scaling, consistent with experiment. We show that the trapping probability can be sensitive to the temperature at both high and low incidence energies. Relaxation and desorption of trapped particles are examined.

Mechanical Vibrational Relaxation of NO Scattering from Metal and Insulator Surfaces: When and Why They Are Different

NO scattering from metallic and insulating surfaces represents contrasting benchmark systems for understanding energy transfer at gas-surface interface. Strikingly different behaviors of highly vibrationally excited NO scattered from Au(111) and LiF(001) were observed and attributed to disparate electronic structures between metals and insulators. Here, we reveal an alternative mechanical origin of this discrepancy by comparative molecular dynamics simulations with globally accurate adiabatic neural network potentials of both systems. We find that highly vibrating NO can reach for the high-dissociation barrier on Au(111), by which vibrational energy can largely transfer to translation or rotation and further dissipate into substrate phonons. This mechanical energy transfer channel is forbidden in the purely repulsive NO/LiF(001) system or for low-vibrating NO on Au(111), where molecular vibration is barely coupled to other degrees of freedom. Our results emphasize that the initial state and potential energy landscape concurrently influence the mechanical energy transfer dynamics of gas-surface scattering.

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.

Time-of-flight measurements of the low-energy scattering of CH4 from Ir(111)

We have measured high-resolution time-of-flight (TOF) spectra of methane scattered from an Ir(111) surface at an incident energy of 81 meV. The angular distributions of scattered CH₄ reveal the presence of a sharp and intense specular peak in addition to sharp features corresponding to rotationally inelastic diffraction (RID) peaks along the two main symmetry directions of Ir(111). TOF spectra have been recorded at several RID positions for the two high-symmetry directions. The data show that the scattering dynamics of CH₄ is more complex than the one reported for H₂/D₂, where energy losses in TOF correspond to the expected excitation/deexcitation RID energy transitions. For CH₄, this is the case only for RID peaks showing up far from the specular peak, whereas those appearing close to the specular position present different behaviors, depending on the incident direction. The results are compared with Ne scattering TOF data, which allows to assess the relevance of multiphonon scattering in the energy-exchange process. Finally, we report experimental evidence of selective adsorption resonances detected with CH₄ beams. This will allow characterizing the CH₄-metal surface physisorption well by measuring angular distributions with CH₄ beams.

Towards bridging the structure gap in heterogeneous catalysis: the impact of defects in dissociative chemisorption of methane on Ir surfaces

The negatively activated region in CH₄ dissociation is attributed to a precursor-mediated mechanism involving surface defects.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Theoretical research on a coke-resistant catalyst for the partial oxidation of methane: Pt/Cu single-atom alloys

Cu can prevent carbon deposition on a surface due to weak adsorption, but it exhibits a high energy barrier to C–H bond activation, which means that it is not practical.

Fundamental mechanisms for molecular energy conversion and chemical reactions at surfaces

The dream of theoretical surface chemistry is to predict the outcome of reactions in order to find the ideal catalyst for a certain application. Having a working ab initio theory in hand would not only enable these predictions but also provide insights into the mechanisms of surface reactions. The development of theoretical models can be assisted by experimental studies providing benchmark data. Though for some reactions a quantitative agreement between experimental observations and theoretical calculations has been achieved, theoretical surface chemistry is in general still far away from gaining predictive power. Here we review recent experimental developments towards the understanding of surface reactions. It is demonstrated how quantum-state resolved scattering experiments on reactive and nonreactive systems can be used to test front-running theoretical approaches. Two challenges for describing dynamics at surfaces are addressed: nonadiabaticity in diatomic molecule surface scattering and the increasing system size when observing and describing the dynamics of polyatomic molecules at surfaces. Finally recent experimental studies on reactive systems are presented. It is shown how elementary steps in a complex surface reaction can be revealed experimentally.

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.

Hydrogen Accumulation and Distribution in Pipeline Steel in Intensified Corrosion Conditions

Hydrogen accumulation and distribution in pipeline steel under conditions of enhanced corrosion has been studied. The XRD analysis, optical spectrometry and uniaxial tension tests reveal that the corrosion environment affects the parameters of the inner and outer surface of the steel pipeline as well as the steel pipeline bulk. The steel surface becomes saturated with hydrogen released as a reaction product during insignificant methane dissociation. Measurements of the adsorbed hydrogen concentration throughout the steel pipe bulk were carried out. The pendulum impact testing of Charpy specimens was performed at room temperature in compliance with national standards. The mechanical properties of the steel specimens were found to be considerably lower, and analogous to the properties values caused by hydrogen embrittlement.

Diffraction of CH4 from a Metal Surface

Diffraction with matter waves has been reported since the beginning of quantum mechanics. In free space, diffraction effects have been observed even with objects as large as C₆₀ molecules. However, in scattering from a solid surface, pure elastic diffraction features have never been observed with molecules larger than D₂. Here we report the observation of pure molecular diffraction for CH₄ scattered off of an Ir(111) surface. These results prove that quantum coherence is preserved, despite the small separation between rotational levels and the interaction with surface phonons. Density functional theory calculations of the potential energy surface provide some clues to understand the larger corrugation sampled by CH₄ molecules in comparison to Ne atoms. Our results show that isotope separation of polyatomic molecules may be possible using gas-surface diffraction.

Test of the Transferability of the Specific Reaction Parameter Functional for H2 + Cu(111) to D2 + Ag(111)

The accurate description of the dissociative chemisorption of a molecule on a metal surface requires a chemically accurate description of the molecule-surface interaction. Previously, it was shown that the specific reaction parameter approach to density functional theory (SRP-DFT) enables accurate descriptions of the reaction of dihydrogen with metal surfaces in, for instance, H₂ + Pt(111), H₂ + Cu(111), and H₂ + Cu(100). SRP-DFT likewise allowed a chemically accurate description of dissociation of methane on Ni(111) and Pt(111), and the SRP functional for CH₄ + Ni(111) was transferable to CH₄ + Pt(111), where Ni and Pt belong to the same group. Here, we investigate whether the SRP density functional derived for H₂ + Cu(111) also gives chemically accurate results for H₂ + Ag(111), where Ag belongs to the same group as Cu. To do this, we have performed quasi-classical trajectory calculations using the six-dimensional potential energy surface of H₂ + Ag(111) within the Born-Oppenheimer static surface approximation. The computed reaction probabilities are compared with both state-resolved associative desorption and molecular beam sticking experiments. Our results do not yet show transferability, as the computed sticking probabilities and initial-state selected reaction probabilities are shifted relative to experiment to higher energies by about 2-3 kcal/mol. The lack of transferability may be due to the different character of the SRP functionals for H₂ + Cu and CH₄ + group 10 metals, the latter containing a van der Waals correlation functional and the former not.

Incident Angle Dependence of CHD3 Dissociation on the Stepped Pt(211) Surface

The dissociation of methane on transition metal surfaces is not only of fundamental interest but also of industrial importance as it represents a rate-controlling step in the steam-reforming reaction used commercially to produce hydrogen. Recently, a specific reaction parameter functional (SRP32-vdW) has been developed, which describes the dissociative chemisorption of CHD₃ at normal incidence on Ni(111), Pt(111), and Pt(211) within chemical accuracy (4.2 kJ/mol). Here, we further test the validity of this functional by comparing the initial sticking coefficients (S₀), obtained from ab-initio molecular dynamics calculations run using this functional, with those measured with the King and Wells method at different angles of incidence for CHD₃ dissociation on Pt(211). The two sets of data are in good agreement, demonstrating that the SRP32-vdW functional also accurately describes CHD₃ dissociation at off-normal angles of incidence. When the direction of incidence is perpendicular to the step edges, an asymmetry is seen in the reactivity with respect to the surface normal, with S₀ being higher when the molecule is directed toward the (100) step rather than the (111) terrace. Although there is a small shadowing effect, the trends in S₀ can be attributed to different activation barriers for different surface sites, which in turn is related to the generalized co-ordination numbers of the surface atom to which the dissociating molecule is adsorbed in the transition state. Consequently, most reactivity is seen on the least co-ordinated step atoms at all angles of incidence.

How Au Outperforms Pt in the Catalytic Reduction of Methane towards Ethane and Molecular Hydrogen

Within the context of a "hydrogen economy", it is paramount to guarantee a stable supply of molecular hydrogen to devices such as fuel cells. At the same time, catalytic conversion of the environmentally harmful methane into ethane, with a significantly lower Global Warming Potential, turns into a highly desirable challenge. Herein we propose a first-step novel proof-of-concept mechanism to accomplish both tasks simultaneously. For that purpose we provide transition-state barriers and reaction Helmholtz free energies obtained from first-principles Density Functional Theory by taking account vibrations for 2CH₄(g) → C₂H₆(g) + H₂(g) to show that molecular hydrogen can be produced by subnanometer Pt₃₈ and Au₃₈ nanoparticles from natural gas. Interestingly, the active sites for the reaction are located on different planes on the two nanoparticles, effectively differentiating the working principle of the two metals. The analysis shows that the complete cycle to reduce CH₄ can be performed on Au and Pt with similar efficiencies, but Au requires only half the working temperature of Pt. This substantial decrease of temperature can be traced back to several intermediate steps, but most crucially to the final one where the catalyst must be cleaned from H(⋆) to be able to restart the catalytic cycle. This simple study case provides useful guidelines to capitalize on finite-size effects in small nanoparticles for the design of new and more efficient catalysts. Interestingly, present results obtained for the intermediate steps of the catalytic cycle show an excellent agreement with previous experimental evidence. Finally, we stress the importance of including the final cleaning steps to start a new fresh catalytic cycle.

Vibrational Energy Redistribution in a Gas-Surface Encounter: State-to-State Scattering of CH_{4} from Ni(111)

The fate of vibrational energy in the collision of methane (CH_{4}) in its antisymmetric C-H stretch vibration (ν_{3}) with a Ni(111) surface has been studied in a state-to-state scattering experiment. Laser excitation in the incident molecular beam prepared the J=1 rotational state of ν_{3}, and a bolometer in combination with selective laser excitation detected the scattered methane. The rovibrationally resolved scattering distributions reveal very efficient vibrational energy redistribution from ν_{3} to the symmetric C-H stretch vibration (ν_{1}). The branching ratio ν_{1}/ν_{3} is near 0.4 and insensitive to changes in incident kinetic energy in the range from 100 to 370 meV. State-resolved angular distributions and measurements of the residual Doppler linewidths prove that the scattering is direct. The observed vibrationally inelastic scattering provides direct experimental evidence for surface-induced vibrational energy redistribution.

Surface Reaction Barriometry: Methane Dissociation on Flat and Stepped Transition-Metal Surfaces

Accurately simulating heterogeneously catalyzed reactions requires reliable barriers for molecules reacting at defects on metal surfaces, such as steps. However, first-principles methods capable of computing these barriers to chemical accuracy have yet to be demonstrated. We show that state-resolved molecular beam experiments combined with ab initio molecular dynamics using specific reaction parameter density functional theory (SRP-DFT) can determine the molecule-metal surface interaction with the required reliability. Crucially, SRP-DFT exhibits transferability: the functional devised for methane reacting on a flat (111) face of Pt (and Ni) also describes its reaction on stepped Pt(211) with chemical accuracy. Our approach can help bridge the materials gap between fundamental surface science studies on regular surfaces and heterogeneous catalysis in which defected surfaces are important.

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.

Coherent quantum scattering of CH4 from Ni(111)

We have measured high-resolution angular distributions of methane scattered from a Ni(111) surface at incident energies between 68.9 meV and 108.6 meV. A sharp and intense specular peak has been observed, in addition to sharp features corresponding to rotationally inelastic diffraction (RID) peaks along the two main symmetry directions of Ni(111). The intensity of the most intense RID peaks is ca. 50% of the specular one. The observation of sharp, coherent elastic peaks at such low incident energies suggests that single scattering dominates over trapping at these energies, and that the depth of the van der Waals well should be lower than 60 meV. In contrast, a broad angular distribution shifted from the specular position is observed from a graphene-covered Ni(111) surface under identical incident conditions. These results open up the possibility of studying the physisorption well between CH₄ and a transition metal surface using high-resolution molecular beams.

The dissociation and recombination rates of CH4 through the Ni(111) surface: The effect of lattice motion

Methane dissociation is a prototypical system for the study of surface reaction dynamics. The dissociation and recombination rates of CH₄ through the Ni(111) surface are calculated by using the quantum instanton method with an analytical potential energy surface. The Ni(111) lattice is treated rigidly, classically, and quantum mechanically so as to reveal the effect of lattice motion. The results demonstrate that it is the lateral displacements rather than the upward and downward movements of the surface nickel atoms that affect the rates a lot. Compared with the rigid lattice, the classical relaxation of the lattice can increase the rates by lowering the free energy barriers. For instance, at 300 K, the dissociation and recombination rates with the classical lattice exceed the ones with the rigid lattice by 6 and 10 orders of magnitude, respectively. Compared with the classical lattice, the quantum delocalization rather than the zero-point energy of the Ni atoms further enhances the rates by widening the reaction path. For instance, the dissociation rate with the quantum lattice is about 10 times larger than that with the classical lattice at 300 K. On the rigid lattice, due to the zero-point energy difference between CH₄ and CD₄, the kinetic isotope effects are larger than 1 for the dissociation process, while they are smaller than 1 for the recombination process. The increasing kinetic isotope effect with decreasing temperature demonstrates that the quantum tunneling effect is remarkable for the dissociation process.