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... Show moreAccurate 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. Show less
Fundamental understanding of molecular reactions on metal surfaces is important for improving heterogeneous catalysis. Therefore, the reaction of small molecules on well-defined metal surfaces is... Show moreFundamental understanding of molecular reactions on metal surfaces is important for improving heterogeneous catalysis. Therefore, the reaction of small molecules on well-defined metal surfaces is investigated with state-of-the-art DFT calculations. Efforts are made to improve the agreement between experiment and theory by employing density functionals belonging to a higher level of theory than typically used. Furthermore, molecular dynamics are performed both with ab initio calculations and precomputed potential energy surfaces to investigate reaction mechanisms. This way dynamical aspects of reaction mechanisms can be investigated, e.g., the effect of rovibrational excitation of a molecule on the reaction probability and mechanism. Show less
The specific reaction parameter (SRP) approach to density functional theory has been shown to model reactions of polyatomic molecules with metal surfaces important for heterogeneous catalysis in... Show moreThe specific reaction parameter (SRP) approach to density functional theory has been shown to model reactions of polyatomic molecules with metal surfaces important for heterogeneous catalysis in the industry with chemical accuracy. However, transferability of the SRP functional among systems in which methane interacts with group 10 metals remains unclear for methane + Pd(111). Therefore, in this work, predictions have been made for the reaction of CHD3 on Pd(111) using Born–Oppenheimer molecular dynamics while also performing a rough comparison with experimental data for CH4 + Pd(111) obtained for lower incidence energies. Hopefully, future experiments can test the transferability of the SRP functional among group 10 metals also for Pd(111). We found that the reactivity of CHD3 on Pd(111) is intermediate between and similar to either Pt(111) or Ni(111), depending on the incidence energy and the initial vibrational state distribution. This is surprising because the barrier height and experiments performed at lower incidence energies than investigated here suggest that the reactivity of Pd(111) should be similar to that of Pt(111) only. The relative decrease in the reactivity of Pd(111) at high incidence energies is attributed to site specificity of the reaction and to dynamical effects such as the bobsled effect and energy transfer from methane to the surface. Show less
An accurate description of reactive scattering of molecules on metal surfaces often requires the modeling of energy transfer between the molecule and the surface phonons. Although ab initio... Show moreAn accurate description of reactive scattering of molecules on metal surfaces often requires the modeling of energy transfer between the molecule and the surface phonons. Although ab initio molecular dynamics (AIMD) can describe this energy transfer, AIMD is at present untractable for reactions with reaction probabilities smaller than 1%. Here, we show that it is possible to use a neural network potential to describe a polyatomic molecule reacting on a mobile metal surface with considerably reduced computational effort compared to AIMD. The highly activated reaction of CHD3 on Cu(111) is used as a test case for this method. It is observed that the reaction probability is influenced considerably by dynamical effects such as the bobsled effect and surface recoil. A special dynamical effect for CHD3 + Cu(111) is that a higher vibrational efficacy is obtained for two quanta in the CH stretch mode than for a single quantum. Show less
An important industrial process is methanol steam reforming, which is typically used in conjunction with copper catalysts. However, little agreement exists on the reaction mechanisms involved on a... Show moreAn important industrial process is methanol steam reforming, which is typically used in conjunction with copper catalysts. However, little agreement exists on the reaction mechanisms involved on a copper catalyst. Therefore, we have performed research yielding additional insight into the reaction mechanism for dissociative chemisorption of methanol on Cu(111) using ab initio molecular dynamics, supported by static calculations of the molecule-surface interaction with density functional theory. Our work predicts that after the initial dissociation, formaldehyde is formed through three different mechanisms. Additionally, it is observed that at high energy, CH cleavage is the dominant pathway instead of the formerly presumed OH cleavage pathway. Finally, in order to describe the interaction of methanol with the metal surface, the SRP32-vdW functional is used, which has been previously developed and tested for CHD3 on Ni(111), Pt(111), and Pt(211) using the Specific Reaction Parameter (SRP) approach. In this work, the SRP32-vdW functional is applied to methanol on Cu(111) as well, in the hope that future experiments can validate the transferability of the SRP32-vdW functional to chemically related molecule-metal surface systems. Show less