Dataset Information

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

ABSTRACT: 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 H₂ 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 H₂, D₂ + Cu(111) with chemical accuracy. Therefore, these three non-empirical functionals themselves are also capable of describing H₂ 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 D₂, H₂ + Ag(111) more accurately than was possible with the SRP48 density functional for H₂ + 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 H₂ + 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 H₂ + Ag(111) on the basis of the three examined MS meta-GGA functionals.

SUBMITTER: Smeets EWF

PROVIDER: S-EPMC6600505 | biostudies-literature | 2019 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H<sub>2</sub> + Cu(111) and H<sub>2</sub> + Ag(111).

Smeets Egidius W F EWF Voss Johannes J Kroes Geert-Jan GJ

The journal of physical chemistry. A 20190613 25

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 ...[more]

PMID: 31149824

Similar Datasets

Project description:It is important that theory is able to accurately describe dissociative chemisorption reactions on metal surfaces, as such reactions are often rate-controlling in heterogeneously catalyzed processes. Chemically accurate theoretical descriptions have recently been obtained on the basis of the specific reaction parameter (SRP) approach to density functional (DF) theory (DFT), allowing reaction barriers to be obtained with chemical accuracy. However, being semiempirical, this approach suffers from two basic problems. The first is that sticking probabilities (to which SRP density functionals (DFs) are usually fitted) might show differences across experiments, of which the origins are not always clear. The second is that it has proven hard to use experiments on diffractive scattering of H2 from metals for validation purposes, as dynamics calculations using a SRP-DF may yield a rather poor description of the measured data, especially if the potential used contains a van der Waals well. We address the first problem by performing dynamics calculations on three sets of molecular beam experiments on D2 + Pt(111), using four sets of molecular beam parameters to obtain sticking probabilities, and the SRP-DF recently fitted to one set of experiments on D2 + Pt(111). It is possible to reproduce all three sets of experiments with chemical accuracy with the aid of two sets of molecular beam parameters. The theoretical simulations with the four different sets of beam parameters allow one to determine for which range of incidence conditions the experiments should agree well and for which conditions they should show specific differences. This allows one to arrive at conclusions about the quality of the experiments and about problems that might affect the experiments. Our calculations on diffraction of H2 scattering from Pt(111) show both quantitative and qualitative differences with previously measured diffraction probabilities, which were Debye-Waller (DW)-extrapolated to 0 K. We suggest that DW extrapolation, which is appropriate for direct scattering, might fail if the scattering is affected by the presence of a van der Waals well and that theory should attempt to model surface atom motion for reproducing diffraction experiments performed for surface temperatures of 500 K and higher.

Project description: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.

Dataset Information

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

Publications

Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H<sub>2</sub> + Cu(111) and H<sub>2</sub> + Ag(111).

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets

Dataset Information

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

Publications

Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H<sub>2</sub> + Cu(111) and H<sub>2</sub> + Ag(111).

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets

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