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

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.

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

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

Computational approaches to dissociative chemisorption on metals: towards chemical accuracy

We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (E_b, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate E_b can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol^-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for E_b for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH₄. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of E_b for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate E_b for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O₂, H₂O, NH₃, CO₂, and CH₃OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.

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

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

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.

Scaling Platinum-Catalyzed Hydrogen Dissociation on Corrugated Surfaces

We determine absolute reactivities for dissociation at low coordinated Pt sites. Two curved Pt(111) single-crystal surfaces allow us to probe either straight or highly kinked step edges with molecules impinging at a low impact energy. A model extracts the average reactivity of inner and outer kink atoms, which is compared to the reactivity of straight A- and B-type steps. Local surface coordination numbers do not adequately capture reactivity trends for H2 dissociation. We utilize the increase of reactivity with step density to determine the area over which a step causes increased dissociation. This step-type specific reactive area extends beyond the step edge onto the (111) terrace. It defines the reaction cross-section for H2 dissociation at the step, bypassing assumptions about contributions of individual types of surface atoms. Our results stress the non-local nature of H2 interaction with a surface and provide insight into reactivity differences for nearly identical step sites.

Chiral Surface Characterisation and Reactivity Toward H–D Exchange of a Curved Platinum Crystal

Understanding heterogeneous catalysis at the atomic level requires detailed knowledge of the reactivity of different surface sites toward specific bond breaking and bond making events. We illustrate a new method in such investigations. We use a macroscopically curved Pt single crystal containing a large variation in density of highly kinked steps of two different chiralities. Scanning tunneling microscopy maps the entire range of surface structures present on the 31° section surrounding the Pt(111) apex. Whereas most of the surface shows the expected characteristic arrays of parallel steps, hexagonally-shaped, single-atom deep pits remain after cleaning procedures near the apex. Their orientation is indicative of the different chiralities present on the two sides of the crystal’s apex. These unintended defects locally raise the surface defect concentration, but are of little consequence to subsequent reactivity measurements for $$\text {D}_2$$ D 2 dissociation and H–D exchange as probed by supersonic molecular beam techniques. We quantify absolute elementary dissociation and relative isotopic exchange rates across the surface with high spatial resolution. At low incident energies, elementary dissociation of the homonuclear isotoplogues is dominated by the kinked steps. H–D exchange kinetics depend also mostly linearly on step density. The changing ratio of D2 dissociation to H–D formation, however, suggests that anisotropic diffusion of H(D) atoms is of influence to the measured HD production rate.

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.