1
|
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
|
2
|
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]
|
3
|
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: 42] [Impact Index Per Article: 14.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
|
4
|
Rivero Santamaría A, Ramos M, Alducin M, Busnengo HF, Díez Muiño R, Juaristi JI. High-Dimensional Atomistic Neural Network Potential to Study the Alignment-Resolved O 2 Scattering from Highly Oriented Pyrolytic Graphite. J Phys Chem A 2021; 125:2588-2600. [PMID: 33734696 DOI: 10.1021/acs.jpca.1c00835] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A high dimensional and accurate atomistic neural network potential energy surface (ANN-PES) that describes the interaction between one O2 molecule and a highly oriented pyrolytic graphite (HOPG) surface has been constructed using the open-source package (aenet). The validation of the PES is performed by paying attention to static characteristics as well as by testing its performance in reproducing previous ab initio molecular dynamics simulation results. Subsequently, the ANN-PES is used to perform quasi-classical molecular dynamics calculations of the alignment-dependent scattering of O2 from HOPG. The results are obtained for 200 meV O2 molecules with different initial alignments impinging with a polar incidence angle with respect to the surface normal of 22.5° on a thermalized (110 and 300 K) graphite surface. The choice of these initial conditions in our simulations is made to perform comparisons to recent experimental results on this system. Our results show that the scattering of O2 from the HOPG surface is a rather direct process, that the angular distributions are alignment dependent, and that the final translational energy of end-on molecules is around 20% lower than that of side-on molecules. Upon collision with the surface, the molecules that are initially aligned perpendicular to the surface become highly rotationally excited, whereas a very small change in the rotational state of the scattered molecules is observed for the initial parallel alignments. The latter confirms the energy transfer dependence on the stereodynamics for the present system. The results of our simulations are in overall agreement with the experimental observations regarding the shape of the angular distributions and the alignment dependence of the in-plane reflected molecules.
Collapse
Affiliation(s)
- Alejandro Rivero Santamaría
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain.,Donostia International Physics Center DIPC, Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Maximiliano Ramos
- Instituto de Física Rosario, CONICET and Universidad Nacional de Rosario, Bv. 27 de Febrero 210 bis, 2000 Rosario, Argentina.,Facultad de Ciencias Exactas, Ingeniera y Agrimensura, Universidad Nacional de Rosario, Av. Pellegrini 250, S2000BTP Rosario, Argentina
| | - Maite Alducin
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain.,Donostia International Physics Center DIPC, Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Heriberto Fabio Busnengo
- Instituto de Física Rosario, CONICET and Universidad Nacional de Rosario, Bv. 27 de Febrero 210 bis, 2000 Rosario, Argentina
| | - Ricardo Díez Muiño
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain.,Donostia International Physics Center DIPC, Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - J Iñaki Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain.,Donostia International Physics Center DIPC, Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain.,Departamento de Polímeros y Materiales Avanzados: Física, Química y Tecnología, Facultad de Químicas, Universidad del País Vasco (UPV/EHU), Apartado 1072, 20080 Donostia-San Sebastián, Spain
| |
Collapse
|
5
|
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
|
6
|
Wijzenbroek M, Kroes GJ. An ab initio molecular dynamics study of D2 dissociation on CO-precovered Ru(0001). Phys Chem Chem Phys 2016; 18:21190-201. [DOI: 10.1039/c6cp00291a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In dynamics studies of hydrogen dissociation on CO pre-covered Ru(0001) the simulation cell size is important for describing energy exchange.
Collapse
Affiliation(s)
- M. Wijzenbroek
- Leiden Institute of Chemistry
- Gorlaeus Laboratories
- Leiden University
- Leiden
- The Netherlands
| | - G. J. Kroes
- Leiden Institute of Chemistry
- Gorlaeus Laboratories
- Leiden University
- Leiden
- The Netherlands
| |
Collapse
|
7
|
Kroes GJ, Pavanello M, Blanco-Rey M, Alducin M, Auerbach DJ. Ab initio molecular dynamics calculations on scattering of hyperthermal H atoms from Cu(111) and Au(111). J Chem Phys 2014; 141:054705. [DOI: 10.1063/1.4891483] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Michele Pavanello
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - María Blanco-Rey
- Departamento de Física de Materiales, Facultad de Químicas UPV/EHU, Apartado 1072, 20080 Donostia-San Sebastián, Spain
- Donostia International Physics Center, Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Maite Alducin
- Donostia International Physics Center, Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Centro de Física de Materiales, Centro Mixto CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Daniel J. Auerbach
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
- Institute for Physical Chemistry, Georg-August University of Göttingen, Göttingen, Germany
| |
Collapse
|
8
|
Pétuya R, Larrégaray P, Crespos C, Busnengo HF, Martínez AE. Dynamics of H2 Eley-Rideal abstraction from W(110): Sensitivity to the representation of the molecule-surface potential. J Chem Phys 2014; 141:024701. [DOI: 10.1063/1.4885139] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- R. Pétuya
- Université de Bordeaux, ISM, CNRS UMR 5255, 33405 Talence Cedex, France
- CNRS, ISM, UMR5255, F-33400 Talence, France
| | - P. Larrégaray
- Université de Bordeaux, ISM, CNRS UMR 5255, 33405 Talence Cedex, France
- CNRS, ISM, UMR5255, F-33400 Talence, France
| | - C. Crespos
- Université de Bordeaux, ISM, CNRS UMR 5255, 33405 Talence Cedex, France
- CNRS, ISM, UMR5255, F-33400 Talence, France
| | - H. F. Busnengo
- Instituto de Física Rosario (IFIR) CONICET-UNR. Ocampo y Esmeralda (2000) Rosario, Argentina
| | - A. E. Martínez
- Instituto de Física Rosario (IFIR) CONICET-UNR. Ocampo y Esmeralda (2000) Rosario, Argentina
| |
Collapse
|
9
|
Meng Q, Chen J, Kilin D. Proton reduction at surface of transition metal nanocatalysts. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.911871] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
10
|
Blanco-Rey M, Juaristi JI, Díez Muiño R, Busnengo HF, Kroes GJ, Alducin M. Electronic friction dominates hydrogen hot-atom relaxation on Pd(100). PHYSICAL REVIEW LETTERS 2014; 112:103203. [PMID: 24679290 DOI: 10.1103/physrevlett.112.103203] [Citation(s) in RCA: 64] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2013] [Indexed: 05/25/2023]
Abstract
We study the dynamics of transient hot H atoms on Pd(100) that originated from dissociative adsorption of H2. The methodology developed here, denoted AIMDEF, consists of ab initio molecular dynamics simulations that include a friction force to account for the energy transfer to the electronic system. We find that the excitation of electron-hole pairs is the main channel for energy dissipation, which happens at a rate that is five times faster than energy transfer into Pd lattice motion. Our results show that electronic excitations may constitute the dominant dissipation channel in the relaxation of hot atoms on surfaces.
Collapse
Affiliation(s)
- M Blanco-Rey
- Departamento de Física de Materiales, Facultad de Químicas UPV/EHU, Apartado 1072, 20018 Donostia-San Sebastián, Spain and Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - J I Juaristi
- Departamento de Física de Materiales, Facultad de Químicas UPV/EHU, Apartado 1072, 20018 Donostia-San Sebastián, Spain and Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain and Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - R Díez Muiño
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain and Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - H F Busnengo
- Instituto de Física Rosario and Universidad Nacional de Rosario, 2000 Rosario, Argentina
| | - G J Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - M Alducin
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain and Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| |
Collapse
|
11
|
|
12
|
Groß A. Coverage effects in the adsorption of H2 on Pd(100) studied by ab initio molecular dynamics simulations. J Chem Phys 2011; 135:174707. [DOI: 10.1063/1.3656765] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Groß A. Ab Initio Molecular Dynamics Simulations of the Adsorption of H2 on Palladium Surfaces. Chemphyschem 2010; 11:1374-81. [DOI: 10.1002/cphc.200900818] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
14
|
Xiao Y, Dong W, Busnengo HF. Reactive force fields for surface chemical reactions: A case study with hydrogen dissociation on Pd surfaces. J Chem Phys 2010; 132:014704. [DOI: 10.1063/1.3265854] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|