Although important to heterogeneous catalysis, the ability to accurately model reactions of polyatomic molecules with metal surfaces has not kept pace with developments in gas phase dynamics. Partnering the specific reaction parameter (SRP) approach to density functional theory with ab initio molecular dynamics (AIMD) extends our ability to model reactions with metals with quantitative accuracy from only the lightest reactant, H2, to essentially all molecules. This is demonstrated with AIMD calculations on CHD3 + Ni(111) in which the SRP functional is fitted to supersonic beam experiments, and validated by showing that AIMD with the resulting functional reproduces initial-state selected sticking measurements with chemical accuracy (4.2 kJ/mol ≈ 1 kcal/mol). The need for only semilocal exchange makes our scheme computationally tractable for dissociation on transition metals.

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.

The theoretical description of methane dissociating on metal surfaces is a current frontier in the field of gas-surface dynamics. Dynamical models that aim at achieving a highly accurate description of this reaction rely on potential energy surfaces based on density functional theory calculations at the generalized gradient approximation. We focus here on the effect that the exchange-correlation functional has on the reactivity of methane on a metal surface, using CHD 3 + Pt(111) as a test case. We present new ab initio molecular dynamics calculations performed with various density functionals, looking also at functionals that account for the van der Waals (vdW) interaction. While searching for a semi-empirical specific reaction parameter density functional for this system, we find that the use of a weighted average of the PBE and the RPBE exchange functionals * Email: f.nattino@chem.leidenuniv.nl 1 2 together with a vdW-corrected correlation functional leads to an improved agreement with quantum state-resolved experimental data for the sticking probability, compared to previous PBE calculations. With this semi-empirical density functional we have also investigated the surface temperature dependence of the methane dissociation reaction and the influence of the rotational alignment on the reactivity, and compared our results with experiments.

In order to model accurately reactions of polyatomic molecules with metal surfaces important for heterogeneous catalysis in industry, the Specific Reaction Parameter (SRP) approach to density functional theory has been developed. This approach has been shown to describe the dissociation of CHD 3 on Ni(111), Pt(111), and Pt(211) with chemical accuracy. In this work, predictions have been made for the reaction of CHD 3 on Cu(111) and Cu(211) using barriers, elbow plots, and ab initio molecular dynamics. Future experiments could hopefully prove the transferability of the SRP functional to systems in which methane reacts with flat and stepped surfaces of adjacent groups of the periodic table, by comparison with our predictions. Moreover, the effect of a so-called single atom alloy on the reactivity of methane is investigated by making predictions for CHD 3 on Pt-Cu(111) and Pd-Cu(111). It is found that the reactivity is only increased for Pt-Cu(111) near the alloyed atom, which is not only caused by the lowering of the barrier height but also by changes in the dynamical pathway and reduction of energy transfer from methane to the surface.

The simulation of the dissociation of molecules on metal surfaces is a cornerstone for the understanding of heterogeneously catalyzed processes. However, due to high computational demand, the accurate dynamical simulation of the dissociative chemisorption of polyatomic molecules has been limited mostly to flat low-index metal surfaces. The study of surfaces that feature "defected" sites, such as steps, is crucial to improve the understanding of the overall catalytic process due to the high reactivity of under-coordinated sites for this kind of reaction. In this work, we have extensively analyzed more than 10 000 molecular dynamics trajectories where a CHD molecule is impinging either on the flat Pt(111) surface or on the stepped Pt(211) surface for different initial rovibrational states and collision energies. The results have been compared in order to get insight into the effect of the step in the dissociation of methane. We have found that, despite a large difference in the activation barrier and consequently in reactivity, the geometry of the lowest transition states is very similar on the two surfaces and this results in a similar dissociation dynamics. Furthermore, the trapping observed on the Pt(211) surface can be explained with energy transfer to parallel translational motion induced by the geometry of the slab and by a larger energy transfer to phonons for the stepped Pt(211) surface.