1
|
Sah MK, Naskar K, Adhikari S, Smits B, Meyer J, Somers MF. On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K. J Chem Phys 2024; 161:014306. [PMID: 38953445 DOI: 10.1063/5.0217639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2024] [Accepted: 06/13/2024] [Indexed: 07/04/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Mantu Kumar Sah
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Bauke Smits
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
2
|
Powell A, Gerrits N, Tchakoua T, Somers MF, Busnengo HF, Meyer J, Kroes GJ, Doblhoff-Dier K. Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals. J Phys Chem Lett 2024; 15:307-315. [PMID: 38169287 PMCID: PMC10788952 DOI: 10.1021/acs.jpclett.3c02972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 12/21/2023] [Accepted: 12/28/2023] [Indexed: 01/05/2024]
Abstract
Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.
Collapse
Affiliation(s)
- Andrew
D. Powell
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Theophile Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Mark F. Somers
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Heriberto F. Busnengo
- Instituto
de Física Rosario (IFIR), CONICET-UNR, 2000 Rosario, Argentina
- Facultad
de Ciencias Exatas, Ingeniería y
Agrimensura, UNR, 2000 Rosario, Argentina
| | - Jörg Meyer
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| |
Collapse
|
3
|
Stark W, Westermayr J, Douglas-Gallardo OA, Gardner J, Habershon S, Maurer RJ. Machine Learning Interatomic Potentials for Reactive Hydrogen Dynamics at Metal Surfaces Based on Iterative Refinement of Reaction Probabilities. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:24168-24182. [PMID: 38148847 PMCID: PMC10749455 DOI: 10.1021/acs.jpcc.3c06648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 11/12/2023] [Accepted: 11/15/2023] [Indexed: 12/28/2023]
Abstract
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.
Collapse
Affiliation(s)
- Wojciech
G. Stark
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Julia Westermayr
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | | | - James Gardner
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Scott Habershon
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| | - Reinhard J. Maurer
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
- Department
of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K.
| |
Collapse
|
4
|
Tchakoua T, Powell AD, Gerrits N, Somers MF, Doblhoff-Dier K, Busnengo HF, Kroes GJ. Simulating Highly Activated Sticking of H 2 on Al(110): Quantum versus Quasi-Classical Dynamics. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:5395-5407. [PMID: 36998253 PMCID: PMC10041643 DOI: 10.1021/acs.jpcc.3c00426] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 02/28/2023] [Indexed: 06/19/2023]
Abstract
We evaluate the importance of quantum effects on the sticking of H2 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 H2 + Al(110) on the basis of the standard procedure for comparing results of theory with molecular beam experiments.
Collapse
Affiliation(s)
- Theophile Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Andrew D. Powell
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F. Somers
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Heriberto F. Busnengo
- Instituto
de Física Rosario (IFIR), CONICET-UNR, Bv. 27 de Febrero 210 bis, 2000 Rosario, Argentina
- Facultad
de Ciencias Exactas, Ingeniería y
Agrimensura, UNR, Av.
Pellegrini 250, 2000 Rosario, Argentina
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
5
|
Smits B, Somers MF. The quantum dynamics of H 2 on Cu(111) at a surface temperature of 925 K: Comparing state-of-the-art theory to state-of-the-art experiments 2. J Chem Phys 2023; 158:014704. [PMID: 36610948 DOI: 10.1063/5.0134817] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
State-of-the-art 6D quantum dynamics simulations for the dissociative chemisorption of H2 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.
Collapse
Affiliation(s)
- B Smits
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Building, 2300 RA Leiden, The Netherlands
| | - M F Somers
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Building, 2300 RA Leiden, The Netherlands
| |
Collapse
|
6
|
Tchakoua T, Gerrits N, Smeets EWF, Kroes GJ. SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces. J Chem Theory Comput 2022; 19:245-270. [PMID: 36529979 PMCID: PMC9835835 DOI: 10.1021/acs.jctc.2c00824] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
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.
Collapse
Affiliation(s)
- T. Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands
| | - N. Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,PLASMANT,
Department of Chemistry, University of Antwerp, BE-2610Antwerp, Belgium
| | - E. W. F. Smeets
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,ALTEN
Nederland, Technology, Fascinatio Boulevard 582, 2909 VACapelle a/d IJssel, The Netherlands
| | - G.-J. Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,E-mail: . Phone: +31 71 527 4396
| |
Collapse
|
7
|
Smits B, Litjens LG, Somers MF. Accurate Description of the Quantum Dynamical Surface Temperature Effects on the Dissociative Chemisorption of H 2 from Cu(111). J Chem Phys 2022; 156:214706. [DOI: 10.1063/5.0094985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Bauke Smits
- Theoretical Chemistry, Leiden University Institute of Chemistry, Netherlands
| | | | - Mark F Somers
- Leiden Institute of Chemistry, Leiden University, Netherlands
| |
Collapse
|
8
|
Formulation of temperature dependent effective Hartree potential incorporating quadratic over linear molecular DOFs-surface modes couplings and its effect on quantum dynamics of D2 (v = 0, j = 0)/D2 (v = 0, j = 2) on Cu(111) metal surface. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
9
|
Roy S, Tiwari A. Mode Selective Chemistry for the Dissociation of Methane on Efficient Ni/Pt-Bimetallic Alloy Catalysts. Phys Chem Chem Phys 2022; 24:16596-16610. [DOI: 10.1039/d2cp02030k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The mode selectivity of methane dissociation is studied on three different Ni/Pt-bimetallic alloy surfaces using a fully quantum approach based on reaction path Hamiltonian. Dissociative sticking probability depends on the...
Collapse
|
10
|
Wei F, Smeets EWF, Voss J, Kroes GJ, Lin S, Guo H. Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Fenfei Wei
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Egidius W. F. Smeets
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Johannes Voss
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park CA 94025, USA
| | - Geert-Jan Kroes
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Sen Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
11
|
Smeets EF, Kroes GJ. Performance of Made Simple Meta-GGA Functionals with rVV10 Nonlocal Correlation for H 2 + Cu(111), D 2 + Ag(111), H 2 + Au(111), and D 2 + Pt(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2021; 125:8993-9010. [PMID: 34084265 PMCID: PMC8162760 DOI: 10.1021/acs.jpcc.0c11034] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Revised: 04/08/2021] [Indexed: 06/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Egidius
W. F. Smeets
- Gorlaeus Laboratories, Leiden
Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Gorlaeus Laboratories, Leiden
Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
12
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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 (Eb, 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 Eb 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 Eb 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 CH4. 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 Eb 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 Eb 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. O2, H2O, NH3, CO2, and CH3OH. 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.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
13
|
Dutta J, Mandal S, Adhikari S, Spiering P, Meyer J, Somers MF. Effect of surface temperature on quantum dynamics of H 2 on Cu(111) using a chemically accurate potential energy surface. J Chem Phys 2021; 154:104103. [PMID: 33722025 DOI: 10.1063/5.0035830] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The effect of surface atom vibrations on H2 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.
Collapse
Affiliation(s)
- Joy Dutta
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Souvik Mandal
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Paul Spiering
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
14
|
Smits B, Somers MF. Beyond the static corrugation model: Dynamic surfaces with the embedded atom method. J Chem Phys 2021; 154:074710. [PMID: 33607896 DOI: 10.1063/5.0036611] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The D2 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.
Collapse
Affiliation(s)
- B Smits
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, RA Leiden 2300, The Netherlands
| | - M F Somers
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, RA Leiden 2300, The Netherlands
| |
Collapse
|
15
|
Powell AD, Kroes GJ, Doblhoff-Dier K. Quantum Monte Carlo calculations on dissociative chemisorption of H2 + Al(110): Minimum barrier heights and their comparison to DFT values. J Chem Phys 2020; 153:224701. [DOI: 10.1063/5.0022919] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrew D. Powell
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Katharina Doblhoff-Dier
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| |
Collapse
|
16
|
Smeets EWF, Kroes GJ. Designing new SRP density functionals including non-local vdW-DF2 correlation for H 2 + Cu(111) and their transferability to H 2 + Ag(111), Au(111) and Pt(111). Phys Chem Chem Phys 2020; 23:7875-7901. [PMID: 33291129 DOI: 10.1039/d0cp05173j] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Specific reaction parameter density functionals (SRP-DFs) that can describe dissociative chemisorption molecular beam experiments of hydrogen (H2) 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 H2 + Cu(111) system, and evaluate their transferability to the highly activated H2 + Ag(111) and H2 + Au(111) systems and the non-activated H2 + Pt(111) system. We design our functionals for the H2 + 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 H2 + Cu(111) with chemical accuracy can also describe such experiments for H2 + 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 H2 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 H2 + 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 H2 on Cu(111) is still well described quasi-classically.
Collapse
Affiliation(s)
- Egidius W F Smeets
- Univeristeit Leiden, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands.
| | | |
Collapse
|
17
|
Galparsoro O, Kaufmann S, Auerbach DJ, Kandratsenka A, Wodtke AM. First principles rates for surface chemistry employing exact transition state theory: application to recombinative desorption of hydrogen from Cu(111). Phys Chem Chem Phys 2020; 22:17532-17539. [DOI: 10.1039/d0cp02858d] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present first principles calculations of the reactive flux for thermal recombinative desorption of hydrogen from Cu(111).
Collapse
Affiliation(s)
- Oihana Galparsoro
- Institute for Physical Chemistry
- Georg-August University of Göttingen
- 37077 Göttingen
- Germany
- Department of Dynamics at Surfaces
| | - Sven Kaufmann
- Department of Dynamics at Surfaces
- Max-Planck-Institute for Biophysical Chemistry
- 37077 Göttingen
- Germany
| | - Daniel J. Auerbach
- Department of Dynamics at Surfaces
- Max-Planck-Institute for Biophysical Chemistry
- 37077 Göttingen
- Germany
| | - Alexander Kandratsenka
- Department of Dynamics at Surfaces
- Max-Planck-Institute for Biophysical Chemistry
- 37077 Göttingen
- Germany
| | - Alec M. Wodtke
- Institute for Physical Chemistry
- Georg-August University of Göttingen
- 37077 Göttingen
- Germany
- Department of Dynamics at Surfaces
| |
Collapse
|
18
|
Zhu L, Zhang Y, Zhang L, Zhou X, Jiang B. Unified and transferable description of dynamics of H2 dissociative adsorption on multiple copper surfaces via machine learning. Phys Chem Chem Phys 2020; 22:13958-13964. [DOI: 10.1039/d0cp02291h] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Schematic of the developed neural network potential energy surface enabling a unified and transferable description of dynamics of H2 dissociative adsorption on multiple copper surfaces.
Collapse
Affiliation(s)
- Lingjun Zhu
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Liang Zhang
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale
- Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes
- Department of Chemical Physics, University of Science and Technology of China
- Hefei
- China
| |
Collapse
|
19
|
Gerrits N, Chadwick H, Kroes GJ. Dynamical Study of the Dissociative Chemisorption of CHD 3 on Pd(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:24013-24023. [PMID: 31602282 PMCID: PMC6778984 DOI: 10.1021/acs.jpcc.9b05757] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Revised: 08/20/2019] [Indexed: 06/10/2023]
Abstract
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 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.
Collapse
Affiliation(s)
- Nick Gerrits
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, the Netherlands
| | - Helen Chadwick
- Department
of Chemistry, Swansea University, Singleton Park, Swansea SA2 8PP, U.K.
| | - Geert-Jan Kroes
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, the Netherlands
| |
Collapse
|
20
|
Smeets EF, Voss J, Kroes GJ. Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H 2 + Cu(111) and H 2 + Ag(111). J Phys Chem A 2019; 123:5395-5406. [PMID: 31149824 PMCID: PMC6600505 DOI: 10.1021/acs.jpca.9b02914] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2019] [Revised: 05/09/2019] [Indexed: 01/08/2023]
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 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.
Collapse
Affiliation(s)
- Egidius
W. F. Smeets
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Johannes Voss
- SLAC
National Accelerator Laboratory, SUNCAT Center Interface Science &
Catalysis, 2575 Sand
Hill Rd, Menlo Park, California 94025, United States
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
21
|
Gerrits N, Shakouri K, Behler J, Kroes GJ. Accurate Probabilities for Highly Activated Reaction of Polyatomic Molecules on Surfaces Using a High-Dimensional Neural Network Potential: CHD 3 + Cu(111). J Phys Chem Lett 2019; 10:1763-1768. [PMID: 30922058 PMCID: PMC6477808 DOI: 10.1021/acs.jpclett.9b00560] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
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 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.
Collapse
Affiliation(s)
- N. Gerrits
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- E-mail:
| | - Khosrow Shakouri
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Behler
- Institut
für Physikalische Chemie, Theoretische Chemie, Universität Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Geert-Jan Kroes
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- E-mail:
| |
Collapse
|
22
|
Ghosh S, Ray D, Tiwari AK. Effects of alloying on mode-selectivity in H2O dissociation on Cu/Ni bimetallic surfaces. J Chem Phys 2019; 150:114702. [DOI: 10.1063/1.5085696] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Smita Ghosh
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
| | - Dhiman Ray
- Department of Chemistry, University of California Irvine, Irvine, California 92617, USA
| | - Ashwani K. Tiwari
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
| |
Collapse
|
23
|
Gerrits N, Kroes GJ. An AIMD study of dissociative chemisorption of methanol on Cu(111) with implications for formaldehyde formation. J Chem Phys 2019; 150:024706. [DOI: 10.1063/1.5070129] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nick Gerrits
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
24
|
Spiering P, Wijzenbroek M, Somers MF. An improved static corrugation model. J Chem Phys 2018; 149:234702. [DOI: 10.1063/1.5058271] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Affiliation(s)
- P. Spiering
- Leiden University, Leiden,
Zuid-holland 2300 RA, The Netherlands
| | | | - M. F. Somers
- Leiden University, Leiden,
Zuid-holland 2300 RA, The Netherlands
| |
Collapse
|
25
|
Gerrits N, Migliorini D, Kroes GJ. Dissociation of CHD3 on Cu(111), Cu(211), and single atom alloys of Cu(111). J Chem Phys 2018; 149:224701. [DOI: 10.1063/1.5053990] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nick Gerrits
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Davide Migliorini
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
26
|
Mandal S, Ghosh S, Sardar S, Adhikari S. The TDDVR approach for molecular photoexcitation, molecule–surface and triatomic reactive scattering processes. INT REV PHYS CHEM 2018. [DOI: 10.1080/0144235x.2018.1548103] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Souvik Mandal
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, India
| | - Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, India
| | | | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, India
| |
Collapse
|
27
|
Spiering P, Meyer J. Testing Electronic Friction Models: Vibrational De-excitation in Scattering of H 2 and D 2 from Cu(111). J Phys Chem Lett 2018; 9:1803-1808. [PMID: 29528648 PMCID: PMC5890313 DOI: 10.1021/acs.jpclett.7b03182] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 03/12/2018] [Indexed: 05/29/2023]
Abstract
At present, molecular dynamics with electronic friction (MDEF) is the workhorse model to go beyond the Born-Oppenheimer approximation in modeling dynamics of molecules at metal surfaces. Concomitant friction coefficients can be calculated with either the local density friction approximation (LDFA) or orbital-dependent friction (ODF), which, unlike LDFA, accounts for anisotropy while relying on other approximations. Due to the computational cost of ODF, extensive high-dimensional MDEF trajectory calculations of experimentally measurable observables have hitherto only been performed based on LDFA. We overcome this limitation with a continuous neural-network-based representation. In our first application to the scattering of vibrationally excited H2 and D2 from Cu(111), we predict up to three times higher vibrational de-excitation probabilities with ODF than with LDFA. These results indicate that anisotropic electronic friction can be important for specific molecular observables. Future experiments can test for this "fingerprint" of different approximations underlying state-of-the-art MDEF.
Collapse
Affiliation(s)
- Paul Spiering
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
28
|
Doblhoff-Dier K, Meyer J, Hoggan PE, Kroes GJ. Quantum Monte Carlo Calculations on a Benchmark Molecule-Metal Surface Reaction: H 2 + Cu(111). J Chem Theory Comput 2017; 13:3208-3219. [PMID: 28514594 PMCID: PMC5508338 DOI: 10.1021/acs.jctc.7b00344] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
Accurate
modeling of heterogeneous catalysis requires the availability
of highly accurate potential energy surfaces. Within density functional
theory, these can—unfortunately—depend heavily on the
exchange-correlation functional. High-level ab initio calculations, on the other hand, are challenging due to the system
size and the metallic character of the metal slab. Here, we present
a quantum Monte Carlo (QMC) study for the benchmark system H2 + Cu(111), focusing on the dissociative chemisorption barrier height.
These computationally extremely challenging ab initio calculations agree to within 1.6 ± 1.0 kcal/mol with a chemically
accurate semiempirical value. Remaining errors, such as time-step
errors and locality errors, are analyzed in detail in order to assess
the reliability of the results. The benchmark studies presented here
are at the cutting edge of what is computationally feasible at the
present time. Illustrating not only the achievable accuracy but also
the challenges arising within QMC in such a calculation, our study
presents a clear picture of where we stand at the moment and which
approaches might allow for even more accurate results in the future.
Collapse
Affiliation(s)
- Katharina Doblhoff-Dier
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , Post Office Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , Post Office Box 9502, 2300 RA Leiden, The Netherlands
| | - Philip E Hoggan
- Institute Pascal, UMR 6602 CNRS, University Blaise Pascal , 4 avenue Blaise Pascal, TSA 60026, CS 60026, 63178 Aubiere Cedex, France
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , Post Office Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
29
|
Shuai Q, Kaufmann S, Auerbach DJ, Schwarzer D, Wodtke AM. Evidence for Electron-Hole Pair Excitation in the Associative Desorption of H 2 and D 2 from Au(111). J Phys Chem Lett 2017; 8:1657-1663. [PMID: 28338329 DOI: 10.1021/acs.jpclett.7b00265] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The dissociative adsorption reaction of hydrogen on noble metals is believed to be well-described within the Born-Oppenheimer approximation. In this work, we have experimentally derived translational energy distributions for selected quantum states of H2 and D2 formed in associative desorption reactions at a Au(111) surface. Using the principle of detailed balance, we compare our results to theory carried out at the same level of sophistication as was done for the reaction on copper. The theory predicts translational excitation that is much higher than is seen in experiment and fails to reproduce the experimentally observed isotope effect. The large deviations between experiment and theory are surprising because, for the same reactions occurring on Cu(111), a similar theoretical strategy agreed with experiment, yielding "chemical accuracy". We argue that electron-hole pair excitation is more important for the reaction on gold, an effect that may be related to the reaction's later transition state.
Collapse
Affiliation(s)
- Quan Shuai
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry , Am Faßberg 11, 37077 Göttingen, Germany
| | - Sven Kaufmann
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry , Am Faßberg 11, 37077 Göttingen, Germany
- Institute for Physical Chemistry, Georg-August University of Goettingen , Tammannstraße 6, 37077 Göttingen, Germany
| | - Daniel J Auerbach
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry , Am Faßberg 11, 37077 Göttingen, Germany
- Institute for Physical Chemistry, Georg-August University of Goettingen , Tammannstraße 6, 37077 Göttingen, Germany
| | - Dirk Schwarzer
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry , Am Faßberg 11, 37077 Göttingen, Germany
- Institute for Physical Chemistry, Georg-August University of Goettingen , Tammannstraße 6, 37077 Göttingen, Germany
| | - Alec M Wodtke
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry , Am Faßberg 11, 37077 Göttingen, Germany
- Institute for Physical Chemistry, Georg-August University of Goettingen , Tammannstraße 6, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Georg-August University of Goettingen , Tammannstraße 6, 37077 Göttingen, Germany
| |
Collapse
|
30
|
Seenivasan H, Jackson B, Tiwari AK. Water dissociation on Ni(100), Ni(110), and Ni(111) surfaces: Reaction path approach to mode selectivity. J Chem Phys 2017; 146:074705. [DOI: 10.1063/1.4976133] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- H. Seenivasan
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
| | - Bret Jackson
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Ashwani K. Tiwari
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
| |
Collapse
|
31
|
Wijzenbroek M, Helstone D, Meyer J, Kroes GJ. Dynamics of H2 dissociation on the close-packed (111) surface of the noblest metal: H2 + Au(111). J Chem Phys 2016; 145:144701. [DOI: 10.1063/1.4964486] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mark Wijzenbroek
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Darcey Helstone
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
32
|
Kroes GJ, Díaz C. Quantum and classical dynamics of reactive scattering of H2 from metal surfaces. Chem Soc Rev 2016; 45:3658-700. [DOI: 10.1039/c5cs00336a] [Citation(s) in RCA: 120] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
State-of-the-art theoretical models allow nowadays an accurate description of H2/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry
- Gorlaeus Laboratories
- Leiden University
- 2300 RA Leiden
- The Netherlands
| | - Cristina Díaz
- Departamento de Química
- Módulo 13
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
| |
Collapse
|
33
|
Kroes GJ. Toward a Database of Chemically Accurate Barrier Heights for Reactions of Molecules with Metal Surfaces. J Phys Chem Lett 2015; 6:4106-14. [PMID: 26722785 DOI: 10.1021/acs.jpclett.5b01344] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
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.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
34
|
Wijzenbroek M, Klein DM, Smits B, Somers MF, Kroes GJ. Performance of a Non-Local van der Waals Density Functional on the Dissociation of H2 on Metal Surfaces. J Phys Chem A 2015; 119:12146-58. [DOI: 10.1021/acs.jpca.5b06008] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mark Wijzenbroek
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA, Leiden, The Netherlands
| | - David M. Klein
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA, Leiden, The Netherlands
| | - Bauke Smits
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA, Leiden, The Netherlands
| | - Mark F. Somers
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA, Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry,
Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA, Leiden, The Netherlands
| |
Collapse
|
35
|
Mandal S, Sahoo T, Ghosh S, Adhikari S. The effect of phonon modes and electron–hole pair couplings on molecule–surface scattering processes. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2015. [DOI: 10.1142/s0219633615500285] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The effect of phonon modes and electron–hole pair (elhp) couplings at different surface temperature on D 2(v = 0, 1; j = 0)– Cu (111) collision has been explored by assuming weakly correlated interactions between molecular Degrees of Freedoms (DOFs) with surface modes and elhp excitations through a Hartree product type wavefunction, where the initial state distributions for the phonon modes and the elhp couplings are incorporated by using Bose–Einstein and Fermi–Dirac probability factors, respectively. We carry out four (4D⊗2D)- and six (6D)- dimensional quantum dynamics on such an effective Hamiltonian, and depict the calculated sticking/transition probabilities and energy transfer from molecule to the surface. The phonon modes slightly affect the sticking probability by broadening the profile, but the transition probability are substantially changed with respect to the rigid surface. On the contrary, the inclusion of elhp coupling along with phonon modes does not change the results much compared to the only phonon case.
Collapse
Affiliation(s)
- Souvik Mandal
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Tapas Sahoo
- Department of Chemical Physics, Wiseman Institute of Science, Rehovot 76100, Israel
| | - Sandip Ghosh
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Satrajit Adhikari
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| |
Collapse
|
36
|
Wang W, Zhao Y. The direct and precursor mediated dissociation rates of H2 on a Ni(111) surface. Phys Chem Chem Phys 2015; 17:5901-12. [PMID: 25630487 DOI: 10.1039/c4cp05624h] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dissociation and recombination rates of physisorbed H2, and the direct and steady state dissociation (i.e., the precursor mediated dissociation) rates of gas phase H2 on Ni(111), as well as the corresponding kinetic isotope effects, are calculated using the quantum instanton method, together with path integral Monte Carlo and adaptive umbrella sampling techniques. All these rates except the recombination one first decrease and then increase with the increasing temperature, and their minimum values appear at about 250, 300 and 250 K, respectively. These non-monotonic behaviors reveal that the quantum effect of H2 should be very remarkable at low temperatures. The steady state rates are smaller than the direct rates at low temperatures, however, they become larger than the direct ones at high temperatures, these two kinds of rates become equal at about 400 and 300 K on the rigid and quantum lattices, respectively. The quantum motion of the lattice can enhance the direct and steady state rates, and it increases the steady state rate much more than the direct one, for instance, the direct and steady state rates on the quantum lattice are 1.30 and 2.08 times larger than that on the rigid one at 300 K. The calculated kinetic isotope effects are much larger than 1, which reveals that H2 always has a larger rate than that of D2, and the direct process predicts much larger kinetic isotope effects than the steady state process at low temperatures. In addition, the kinetic isotope effects are not affected by the lattice motion.
Collapse
Affiliation(s)
- Wenji Wang
- College of Science, Northwest A&F University, Yangling, 712100, Shaanxi Province, P. R. China.
| | | |
Collapse
|
37
|
Nattino F, Genova A, Guijt M, Muzas AS, Díaz C, Auerbach DJ, Kroes GJ. Dissociation and recombination of D2 on Cu(111): Ab initio molecular dynamics calculations and improved analysis of desorption experiments. J Chem Phys 2014; 141:124705. [DOI: 10.1063/1.4896058] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Francesco Nattino
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Alessandro Genova
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Marieke Guijt
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Alberto S. Muzas
- Departamento de Química Módulo 13, Universitad Autónoma de Madrid, 28049 Madrid, Spain
| | - Cristina Díaz
- Departamento de Química Módulo 13, Universitad Autónoma de Madrid, 28049 Madrid, Spain
| | - Daniel J. Auerbach
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
38
|
Nave S, Tiwari AK, Jackson B. Dissociative Chemisorption of Methane on Ni and Pt Surfaces: Mode-Specific Chemistry and the Effects of Lattice Motion. J Phys Chem A 2014; 118:9615-31. [DOI: 10.1021/jp5063644] [Citation(s) in RCA: 98] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sven Nave
- Institut des Sciences
Moléculaires d’Orsay, Université Paris-Sud 11/CNRS-UMR
8214, Université Paris-Sud, Bât. 351, 91405 Orsay Cedex, France
| | - Ashwani K. Tiwari
- Indian Institute of Science Education and Research Kolkata, Mohanpur 741252, India
| | - Bret Jackson
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, United States
| |
Collapse
|
39
|
Seenivasan H, Tiwari AK. Water adsorption and dissociation on Ni(110): How is it different from its close packed counterparts? J Chem Phys 2014; 140:174704. [DOI: 10.1063/1.4873898] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Wijzenbroek M, Kroes GJ. The effect of the exchange-correlation functional on H2 dissociation on Ru(0001). J Chem Phys 2014; 140:084702. [DOI: 10.1063/1.4865946] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Wang W, Zhao Y. Dissociation rates of H2on a Ni(100) surface: the role of the physisorbed state. Phys Chem Chem Phys 2014; 16:13318-28. [DOI: 10.1039/c4cp01705f] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
42
|
Boereboom JM, Wijzenbroek M, Somers MF, Kroes GJ. Towards a specific reaction parameter density functional for reactive scattering of H2 from Pd(111). J Chem Phys 2013; 139:244707. [DOI: 10.1063/1.4851355] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Seenivasan H, Tiwari AK. Water dissociation on Ni(100) and Ni(111): Effect of surface temperature on reactivity. J Chem Phys 2013; 139:174707. [DOI: 10.1063/1.4827641] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|