A search for vibrational contributions to the activated adsorption of H2 on copper

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.

Machine Learning Interatomic Potentials for Reactive Hydrogen Dynamics at Metal Surfaces Based on Iterative Refinement of Reaction Probabilities

The reactive chemistry of molecular hydrogen at surfaces, notably dissociative sticking and hydrogen evolution, plays a crucial role in energy storage and fuel cells. Theoretical studies can help to decipher underlying mechanisms and reaction design, but studying dynamics at surfaces is computationally challenging due to the complex electronic structure at interfaces and the high sensitivity of dynamics to reaction barriers. In addition, ab initio molecular dynamics, based on density functional theory, is too computationally demanding to accurately predict reactive sticking or desorption probabilities, as it requires averaging over tens of thousands of initial conditions. High-dimensional machine learning-based interatomic potentials are starting to be more commonly used in gas-surface dynamics, yet robust approaches to generate reliable training data and assess how model uncertainty affects the prediction of dynamic observables are not well established. Here, we employ ensemble learning to adaptively generate training data while assessing model performance with full uncertainty quantification (UQ) for reaction probabilities of hydrogen scattering on different copper facets. We use this approach to investigate the performance of two message-passing neural networks, SchNet and PaiNN. Ensemble-based UQ and iterative refinement allow us to expose the shortcomings of the invariant pairwise-distance-based feature representation in the SchNet model for gas-surface dynamics.

Simulating Highly Activated Sticking of H2 on Al(110): Quantum versus Quasi-Classical Dynamics

We evaluate the importance of quantum effects on the sticking of H₂ on Al(110) for conditions that are close to those of molecular beam experiments that have been done on this system. Calculations with the quasi-classical trajectory (QCT) method and with quantum dynamics (QD) are performed using a model in which only motion in the six molecular degrees of freedom is allowed. The potential energy surface used has a minimum barrier height close to the value recently obtained with the quantum Monte Carlo method. Monte Carlo averaging over the initial rovibrational states allowed the QD calculations to be done with an order of magnitude smaller computational expense. The sticking probability curve computed with QD is shifted to lower energies relative to the QCT curve by 0.21 to 0.05 kcal/mol, with the highest shift obtained for the lowest incidence energy. Quantum effects are therefore expected to play a small role in calculations that would evaluate the accuracy of electronic structure methods for determining the minimum barrier height to dissociative chemisorption for H₂ + Al(110) on the basis of the standard procedure for comparing results of theory with molecular beam experiments.

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.

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.

Performance of Made Simple Meta-GGA Functionals with rVV10 Nonlocal Correlation for H2 + Cu(111), D2 + Ag(111), H2 + Au(111), and D2 + Pt(111)

Accurately modeling heterogeneous catalysis requires accurate descriptions of rate-controlling elementary reactions of molecules on metal surfaces, but standard density functionals (DFs) are not accurate enough for this. The problem can be solved with the specific reaction parameter approach to density functional theory (SRP-DFT), but the transferability of SRP DFs among chemically related systems is limited. We combine the MS-PBEl, MS-B86bl, and MS-RPBEl semilocal made simple (MS) meta-generalized gradient approximation (GGA) (mGGA) DFs with rVV10 nonlocal correlation, and we evaluate their performance for the hydrogen (H2) + Cu(111), deuterium (D2) + Ag(111), H2 + Au(111), and D2 + Pt(111) gas-surface systems. The three MS mGGA DFs that have been combined with rVV10 nonlocal correlation were not fitted to reproduce particular experiments, nor has the b parameter present in rVV10 been reoptimized. Of the three DFs obtained the MS-PBEl-rVV10 DF yields an excellent description of van der Waals well geometries. The three original MS mGGA DFs gave a highly accurate description of the metals, which was comparable in quality to that obtained with the PBEsol DF. Here, we find that combining the three original MS mGGA DFs with rVV10 nonlocal correlation comes at the cost of a slightly less accurate description of the metal. However, the description of the metal obtained in this way is still better than the descriptions obtained with SRP DFs specifically optimized for individual systems. Using the Born-Oppenheimer static surface (BOSS) model, simulations of molecular beam dissociative chemisorption experiments yield chemical accuracy for the D2 + Ag(111) and D2 + Pt(111) systems. A comparison between calculated and measured E 1/2(ν, J) parameters describing associative desorption suggests chemical accuracy for the associative desorption of H2 from Au(111) as well. Our results suggest that ascending Jacob's ladder to the mGGA rung yields increasingly more accurate results for gas-surface reactions of H2 (D2) interacting with late transition metals.

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.

Designing new SRP density functionals including non-local vdW-DF2 correlation for H2 + Cu(111) and their transferability to H2 + Ag(111), Au(111) and Pt(111)

Specific reaction parameter density functionals (SRP-DFs) that can describe dissociative chemisorption molecular beam experiments of hydrogen (H₂) on cold transition metal surfaces with chemical accuracy have so far been shown to be only transferable among different facets of the same metal, but not among different metals. We design new SRP-DFs that include non-local vdW-DF2 correlation for the H₂ + Cu(111) system, and evaluate their transferability to the highly activated H₂ + Ag(111) and H₂ + Au(111) systems and the non-activated H₂ + Pt(111) system. We design our functionals for the H₂ + Cu(111) system since it is the best studied system both theoretically and experimentally. Here we demonstrate that a SRP-DF fitted to reproduce molecular beam sticking experiments for H₂ + Cu(111) with chemical accuracy can also describe such experiments for H₂ + Pt(111) with chemical accuracy, and vice versa. Chemically accurate functionals have been obtained that perform very well with respect to reported van der Waals well geometries, and which improve the description of the metal over current generalized gradient approximation (GGA) based SRP-DFs. From a systematic comparison of our new SRP-DFs that include non-local correlation to previously developed SRP-DFs, for both activated and non-activated systems, we identify non-local correlation as a key ingredient in the construction of transferable SRP-DFs for H₂ interacting with transition metals. Our results are in excellent agreement with experiment when accurately measured observables are available. It is however clear from our analysis that, except for the H₂ + Cu(111) system, there is a need for more, more varied, and more accurately described experiments in order to further improve the design of SRP-DFs. Additionally, we confirm that, when including non-local correlation, the sticking of H₂ on Cu(111) is still well described quasi-classically.

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.

Correction to "Quantum Dynamics of Dissociative Chemisorption of H2 on the Stepped Cu(211) Surface″

[This corrects the article DOI: 10.1021/acs.jpcc.9b06539.].

Quantum Dynamics of Dissociative Chemisorption of H2 on the Stepped Cu(211) Surface

Reactions on stepped surfaces are relevant to heterogeneous catalysis, in which a reaction often takes place at the edges of nanoparticles where the edges resemble steps on single-crystal stepped surfaces. Previous results on H₂ + Cu(211) showed that, in this system, steps do not enhance the reactivity and raised the question of whether this effect could be, in any way, related to the neglect of quantum dynamical effects in the theory. To investigate this, we present full quantum dynamical molecular beam simulations of sticking of H₂ on Cu(211), in which all important rovibrational states populated in a molecular beam experiment are taken into account. We find that the reaction of H₂ with Cu(211) is very well described with quasi-classical dynamics when simulating molecular beam sticking experiments, in which averaging takes place over a large number of rovibrational states and over translational energy distributions. Our results show that the stepped Cu(211) surface is distinct from its component Cu(111) terraces and Cu(100) steps and cannot be described as a combination of its component parts with respect to the reaction dynamics when considering the orientational dependence. Specifically, we present evidence that, at translational energies close to the reaction threshold, vibrationally excited molecules show a negative rotational quadrupole alignment parameter on Cu(211), which is not found on Cu(111) and Cu(100). The effect arises because these molecules react with a site-specific reaction mechanism at the step, that is, inelastic rotational enhancement, which is only effective for molecules with a small absolute value of the magnetic rotation quantum number. From a comparison to recent associative desorption experiments as well as Born-Oppenheimer molecular dynamics calculations, it follows that the effects of surface atom motion and electron-hole pair excitation on the reactivity fall within chemical accuracy, that is, modeling these effect shifts extracted reaction probability curves by less than 1 kcal/mol translational energy. We found no evidence in our fully state-resolved calculations for the "slow" reaction channel that was recently reported for associative desorption of H₂ from Cu(111) and Cu(211), but our results for the fast channel are in good agreement with the experiments on H₂ + Cu(211).

Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H2 + Cu(111) and H2 + Ag(111)

Specific reaction parameter density functionals (SRP-DFs), which can describe dissociative chemisorption reactions on metals to within chemical accuracy, have so far been based on exchange functionals within the generalized gradient approximation (GGA) and on GGA correlation functionals or van der Waals correlation functionals. These functionals are capable of describing the molecule-metal surface interaction accurately, but they suffer from the general GGA problem that this can be done only at the cost of a rather poor description of the metal. Here, we show that it is possible also to construct SRP-DFs for H2 dissociation on Cu(111) based on meta-GGA functionals, introducing three new functionals based on the "made-simple" (MS) concept. The exchange parts of the three functionals (MS-PBEl, MS-B86bl, and MS-RPBEl) are based on the expressions for the PBE, B86b, and RPBE exchange functionals. Quasi-classical trajectory (QCT) calculations performed with potential energy surfaces (PESs) obtained with the three MS functionals reproduce molecular beam experiments on H2, D2 + Cu(111) with chemical accuracy. Therefore, these three non-empirical functionals themselves are also capable of describing H2 dissociation on Cu(111) with chemical accuracy. Similarly, QCT calculations performed on the MS-PBEl and MS-B86bl PESs reproduced molecular beam and associative desorption experiments on D2, H2 + Ag(111) more accurately than was possible with the SRP48 density functional for H2 + Cu(111). Also, the three new MS functionals describe the Cu, Ag, Au, and Pt metals more accurately than the all-purpose Perdew-Burke-Ernzerhof (PBE) functional. The only disadvantage we noted of the new MS functionals is that, as found for the example of H2 + Cu(111), the reaction barrier height obtained by taking weighted averages of the MS-PBEl and MS-RPBEl functionals is tunable over a smaller range (9 kJ/mol) than possible with the standard GGA PBE and RPBE functionals (33 kJ/mol). As a result of this restricted tunability, it is not possible to construct an SRP-DF for H2 + Ag(111) on the basis of the three examined MS meta-GGA functionals.

Hydrogen adsorption and desorption from Cu(111) and Cu(211)

We present a combined experimental-theoretical study on structural and coverage dependences of the adsorption and desorption of molecular hydrogen on atomically flat Cu(111) and highly stepped Cu(211) surfaces. For molecules with identical incident energy from supersonic molecular beams, we find a reduced dissociative sticking probability for the stepped surface compared to Cu(111). DFT calculations of activation barriers to dissociation for the clean and partially precovered surfaces, as well as quantitative analysis of TPD spectra, support that the A-type step of the (211) surface causes an upward shift in activation barriers to dissociation and lowering of the desorption barrier. The new data allow us to determine low sticking probabilities at conditions where King and Wells measurements fail to determine the reactivity. They are also fully consistent with the unexpected observation that monoatomic steps on a surface lower the reactivity toward the dissociation of a diatomic molecule.

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.

Toward a Database of Chemically Accurate Barrier Heights for Reactions of Molecules with Metal Surfaces

Being able to calculate reaction barrier heights to within chemical accuracy (errors < 1 kcal/mol) is crucial to the accurate modeling of chemical reactions. Although accurate databases exist that can help theorists with benchmarking new electronic structure theories on gas-phase chemical reactions, no such databases exist for reactions of molecules with metal surfaces. Nonetheless, most chemicals are made in heterogeneously catalyzed processes, of which many take place over metal particles. Presently, barrier heights for molecule-metal surface reactions have been determined with chemical accuracy for only two systems, that is, H2 + Cu(111) and H2 + Cu(100). This has been done with semiempirically determined density functionals, which were fitted through comparisons of dynamics results with molecular beam sticking probabilities. The prospects of extending the database with chemically accurate data for other molecule-metal reactions, either with the use of semiempirical density functional theory or with first-principles theory, are discussed.

Significant quantum effects in hydrogen activation

Dissociation of molecular hydrogen is an important step in a wide variety of chemical, biological, and physical processes. Due to the light mass of hydrogen, it is recognized that quantum effects are often important to its reactivity. However, understanding how quantum effects impact the reactivity of hydrogen is still in its infancy. Here, we examine this issue using a well-defined Pd/Cu(111) alloy that allows the activation of hydrogen and deuterium molecules to be examined at individual Pd atom surface sites over a wide range of temperatures. Experiments comparing the uptake of hydrogen and deuterium as a function of temperature reveal completely different behavior of the two species. The rate of hydrogen activation increases at lower sample temperature, whereas deuterium activation slows as the temperature is lowered. Density functional theory simulations in which quantum nuclear effects are accounted for reveal that tunneling through the dissociation barrier is prevalent for H2 up to ∼190 K and for D2 up to ∼140 K. Kinetic Monte Carlo simulations indicate that the effective barrier to H2 dissociation is so low that hydrogen uptake on the surface is limited merely by thermodynamics, whereas the D2 dissociation process is controlled by kinetics. These data illustrate the complexity and inherent quantum nature of this ubiquitous and seemingly simple chemical process. Examining these effects in other systems with a similar range of approaches may uncover temperature regimes where quantum effects can be harnessed, yielding greater control of bond-breaking processes at surfaces and uncovering useful chemistries such as selective bond activation or isotope separation.

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.

Interaction of Hydrogen with Bimetallic Surfaces

The interaction of hydrogen with bimetallic surfaces is assessed and discussed using structurally well defined PdAu and PtRu single crystal surfaces as model systems. The description is based on a correlation between the hydrogen adsorption behavior and the presence and abundance of specific bimetallic nanostructures. The latter are determined quantitatively by high resolution scanning tunneling microscopy. Following the changes in the hydrogen adsorption behavior on surfaces with systematically varied composition and thus varying concentrations of specific structural elements and adsorption ensembles allows a direct correlation between structural and chemical properties and thus to identify the adsorption characteristics of individual adsorption ensembles. This provides a structural basis for parallel theoretical studies, which are essential for a systematic understanding of the complex changes in the adsorption properties.Based on results from experimental studies and calculations, the role of electronic strain and ligand effects, including lateral and vertical ligands, as well the influence of geometric ensemble effects on the hydrogen adsorption characteristics are discussed. The distinct differences between surfaces modified by a monolayer or monolayer islands of a second metal and mixed surfaces such as surface alloys or alloy surfaces are discussed in these terms.

Dissociative chemisorption of H2 on the Cu(110) surface: a quantum and quasiclassical dynamical study

Six-dimensional quantum dynamical and quasiclassical trajectory (QCT) calculations are reported for the reaction and vibrationally inelastic scattering of (v = 0,1,j = 0) H(2) scattering from Cu(110), and for the reaction and rovibrationally elastic and inelastic scattering of (v = 1,j = 1) H(2) scattering from Cu(110). The dynamics results were obtained using a potential energy surface obtained with density functional theory using the PW91 functional. The reaction probabilities computed with quantum dynamics for (v = 0,1,j = 0) were in excellent agreement with the QCT results obtained earlier for these states, thereby validating the QCT approach to sticking of hydrogen on Cu(110). The vibrational de-excitation probability P(v=1,j = 0 --> v = 0) computed with the QCT method is in remarkably good agreement with the quantum dynamical results for normal incidence energies E(n) between 0.2 and 0.6 eV. The QCT result for the vibrational excitation probability P(v = 0,j = 0 --> v = 1) is likewise accurate for E(n) between 0.8 and 1 eV, but the QCT method overestimates vibrational excitation for lower E(n). The QCT method gives probabilities for rovibrationally (in)elastic scattering, P(v = 1,j = 1 --> v('),j(')), which are in remarkably good agreement with quantum dynamical results. The rotationally averaged, initial vibrational state-selective reaction probability obtained with QCT agrees well with the initial vibrational state-selective reaction probability extracted from molecular beam experiments for v = 1, for the range of collision energies for which the v=1 contribution to the measured total sticking probability dominates. The quantum dynamical probabilities for rovibrationally elastic scattering of (v = 1,j = 1) H(2) from Cu(110) are in good agreement with experiment for E(n) between 0.08 and 0.25 eV.

Microcanonical Transition State Theory for Activated Gas−Surface Reaction Dynamics: Application to H2/Cu(111) with Rotation as a Spectator

A microcanonical unimolecular rate theory (MURT) model incorporating quantized surface vibrations and Rice-Ramsperger-Kassel-Marcus rate constants is applied to a benchmark system for gas-surface reaction dynamics, the activated dissociative chemisorption and associative desorption of hydrogen on Cu(111). Both molecular translation parallel to the surface and rotation are treated as spectator degrees of freedom. MURT analysis of diverse experiments indicates that one surface oscillator participates in the dissociative transition state and that the threshold energy for H2 dissociation on Cu(111) is E0 = 62 kJ/mol. The spectator approximation for rotation holds well at thermally accessible rotational energies (i.e., for Er less than approximately 40 kJ/mol). Over the temperature range from 300 to 1000 K, the calculated thermal dissociative sticking coefficient is ST = S0 exp(-Ea/kBT) where S0 = 1.57 and Ea = 62.9 kJ/mol. The sigmoid shape of rovibrational eigenstate-resolved dissociative sticking coefficients as a function of normal translational energy is shown to derive from an averaging of the microcanonical sticking coefficient, with threshold energy E0, over the thermal surface oscillator distribution of the gas-surface collision complexes. Given that H2/Cu(111) is one of the most dynamically biased of gas-surface reactive systems, the simple statistical MURT model simulates and broadly rationalizes the H2/Cu(111) reactive behavior with remarkable fidelity.

Seven-dimensional microcanonical treatment of hydrogen dissociation dynamics on Cu(111): Clarifying the essential role of surface phonons

A simple picture of the hydrogen dissociation/associative desorption dynamics on Cu(111) emerges from a two-parameter, full dimensionality microcanonical unimolecular rate theory (MURT) model of the gas-surface reactivity. Vibrational frequencies for the reactive transition state were taken from density functional theory calculations of a six-dimensional potential energy surface [Hammer et al., Phys. Rev. Lett. 73, 1400 (1994)]. The two remaining parameters required by the MURT were fixed by simulation of experiments. These parameters are the dissociation threshold energy, E(0)=79 kJmol, and the number of surface oscillators involved in the localized H(2)Cu(111) collision complex, s=1. The two-parameter MURT quantitatively predicts much of the varied behavior observed for the H(2) and D(2)Cu(111) reactive systems, including the temperature-dependent associative desorption angular distributions, mean translational energies of the associatively desorbing hydrogen as a function of rovibrational eigenstate, etc. The divergence of the statistical theory's predictions from experimental results at low rotational quantum numbers, J < or approximately 5, suggests that either (i) rotational steering is important to the dissociation dynamics at low J, an effect that washes out at high J, or (ii) molecular rotation is approximately a spectator degree of freedom to the dissociation dynamics for these low J states, the states that dominate the thermal reactivity. Surface vibrations are predicted to provide approximately 30% of the energy required to surmount the activation barrier to H(2) dissociation under thermal equilibrium conditions. The MURT with s=1 is used to analytically confirm the experimental finding that partial differential "E(a)(T(s))" partial differential E(t)= -1 for eigenstate-resolved dissociative sticking at translational energies E(t)<E(0)-E(v)-E(r). Explicit treatment of the surface motion (i.e., surface not frozen at T(s)=0 K) is a relatively novel aspect of the MURT theoretical approach.

Theoretical study of hydrogen dissociative adsorption on the Cu(110) surface

We have calculated the six-dimensional (6D) potential energy surface for H2 in front of a frozen Cu(110) surface using density functional theory for 22 H2-surface configurations and the corrugation reducing procedure to interpolate between them. We carry out classical trajectory calculations on the dissociative adsorption process and find excellent agreement with measurements. We find that it is of prominent importance to account for the rovibrational state distribution in the incident H2 beam. A straightforward analysis leads to the conclusion that the motion along the surface does not play an appreciable role in the dynamics whereas the dynamical role of molecular rotation is crucial. The latter fact precludes any interpretation of dissociation in terms of a static concept such as "barrier distributions."

Six-dimensional quantum dynamics of dissociative chemisorption of H2 on Ru(0001)

Six-dimensional quantum dynamics calculations on dissociative chemisorption of H(2) on Ru(0001) are performed. The six-dimensional potential energy surface is generated using density functional theory. Two different generalized gradient approximations are used, i.e., RPBE and PW91, to allow the results to be compared. The dissociation probability for normally incident H(2) on a clean Ru(0001) surface is calculated. Large differences between the reaction probabilities calculated using the RPBE and PW91 are seen, with the PW91 results showing a much narrower reaction probability curve and a much higher reactivity. Using the reaction probabilities and assuming normal energy scaling reaction rates are generated for temperatures between 300 and 800 K. The rate generated using the PW91 results is higher by about a factor 5 than the rate based on the RPBE results in the range of temperatures relevant to ammonia production.