1
|
Shi H, Liu T, Fu Y, Wu H, Fu B, Zhang DH. Fundamental invariant-neural network potential energy surface and dissociative chemisorption dynamics of N2 on rigid Ni(111). COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
2
|
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
|
3
|
Rodríguez-Fernández A, Bonnet L, Larrégaray P, Muiño RD. Ab initio molecular dynamics of hydrogen on tungsten surfaces. Phys Chem Chem Phys 2021; 23:7919-7925. [PMID: 33347525 DOI: 10.1039/d0cp05423b] [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 dissociation process of hydrogen molecules on W(110) was studied using density functional theory and classical molecular dynamics. We have calculated the dissociation probability for molecules with energies below 300 meV and analyzed the dynamics of the adsorption process. Our results show that the fate of each trajectory is determined at distances relatively far from the surface, at roughly 2-2.5 Å. This distance varies slightly with the initial kinetic energy of the molecule. Part of our simulations include van der Waals dispersion effects in the interaction between molecule and surface. We present a comparison between these results and other theoretical and experimental results previously published. The inclusion of the van der Waals term provokes an increase in the far-distance attraction that is compensated by a stronger repulsion at short distances. The combination of both effects appreciably decreases the value of the dissociation probability. The successful comparison of our results with experimental information confirms that the methodology employed can be considered as a rich and accurate instrument to study the dissociation of hydrogen on surfaces.
Collapse
|
4
|
Pétuya C, Martin-Gondre L, Aurel P, Damay F, Desmedt A. Unraveling the metastability of the SI and SII carbon monoxide hydrate with a combined DFT-neutron diffraction investigation. J Chem Phys 2019; 150:184705. [PMID: 31091912 DOI: 10.1063/1.5093202] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Clathrate hydrates are crystalline compounds consisting of water molecules forming cages (so-called "host") inside of which "guest" molecules are encapsulated depending on the thermodynamic conditions of formation (systems stable at low temperature and high pressure). These icelike systems are naturally abundant on Earth and are generally expected to exist on icy celestial bodies. Carbon monoxide hydrate might be considered an important component of the carbon cycle in the solar system since CO gas is one of the predominant forms of carbon. Intriguing fundamental properties have also been reported: the CO hydrate initially forms in the sI structure (kinetically favored) and transforms into the sII structure (thermodynamically stable). Understanding and predicting the gas hydrate structural stability then become essential. The aim of this work is, thereby, to study the structural and energetic properties of the CO hydrate using density functional theory (DFT) calculations together with neutron diffraction measurements. In addition to the comparison of DFT-derived structural properties with those from experimental neutron diffraction, the originality of this work lies in the DFT-derived energy calculations performed on a complete unit cell (sI and sII) and not only by considering guest molecules confined in an isolated water cage (as usually performed for extracting the binding energies). Interestingly, an excellent agreement (within less than 1% error) is found between the measured and DFT-derived unit cell parameters by considering the Perdew-Burke-Ernzerhof (denoted PBE) functional. Moreover, a strategy is proposed for evaluating the hydrate structural stability on the basis of potential energy analysis of the total nonbonding energies (i.e., binding energy and water substructure nonbonding energy). It is found that the sII structure is the thermodynamically stable hydrate phase. In addition, increasing the CO content in the large cages has a stabilizing effect on the sII structure, while it destabilizes the sI structure. Such findings are in agreement with the recent experimental results evidencing the structural metastability of the CO hydrate.
Collapse
Affiliation(s)
- Claire Pétuya
- Institut des Sciences Moléculaires, UMR5255 CNRS-Université de Bordeaux, Talence, France
| | | | - Philippe Aurel
- Institut des Sciences Moléculaires, UMR5255 CNRS-Université de Bordeaux, Talence, France
| | - Françoise Damay
- Laboratoire Léon Brillouin UMR 12 CEA-CNRS, Gif-sur-Yvette, France
| | - Arnaud Desmedt
- Institut des Sciences Moléculaires, UMR5255 CNRS-Université de Bordeaux, Talence, France
| |
Collapse
|
5
|
Peña-Torres A, Busnengo HF, Juaristi JI, Larregaray P, Crespos C. Dynamics of N2 sticking on W(100): the decisive role of van der Waals interactions. Phys Chem Chem Phys 2018; 20:19326-19331. [DOI: 10.1039/c8cp03515f] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The reactive dynamics of N2 on W(100) has been investigated by means of quasi-classical trajectory calculations using an interpolated six-dimensional potential energy surface (PES) based on density functional theory energies obtained employing the vdW-DF2 functional.
Collapse
Affiliation(s)
| | - H. Fabio Busnengo
- Instituto de Física de Rosario (CONICET-UNR) and Facultad de Ciencias Exactas
- Ingeniería y Agrimensura
- Universidad Nacional de Rosario
- 2000 Rosario
- Argentina
| | - J. Iñaki Juaristi
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-San Sebastián
- Spain
- Departamento de Física de Materiales
- Facultad de Químicas (UPV/EHU)
| | | | | |
Collapse
|
6
|
Novko D, Tremblay JC, Blanco-Rey M. On the tautomerisation of porphycene on copper (111): Finding the subtle balance between van der Waals interactions and hybridisation. J Chem Phys 2016; 145:244701. [DOI: 10.1063/1.4972213] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Dino Novko
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Jean Christophe Tremblay
- Institute for Chemistry and Biochemistry, Freie Universität Berlin, Takustraße 3, D-14195 Berlin, Germany
| | - María Blanco-Rey
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Departamento de Física de Materiales, Facultad de Químicas UPV/EHU, Apartado 1072, 20018 Donostia-San Sebastián, Spain
| |
Collapse
|
7
|
Nattino F, Galparsoro O, Costanzo F, Díez Muiño R, Alducin M, Kroes GJ. Modeling surface motion effects in N2 dissociation on W(110): Ab initio molecular dynamics calculations and generalized Langevin oscillator model. J Chem Phys 2016; 144:244708. [DOI: 10.1063/1.4954773] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [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
| | - Oihana Galparsoro
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Université de Bordeaux, F-33400 Talence, France
- CNRS, ISM, UMR 5255, F-33400 Talence, France
| | - Francesca Costanzo
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Ricardo Díez Muiño
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Maite Alducin
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
8
|
Migliorini D, Nattino F, Kroes GJ. Application of van der Waals functionals to the calculation of dissociative adsorption of N2 on W(110) for static and dynamic systems. J Chem Phys 2016; 144:084702. [DOI: 10.1063/1.4942198] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Davide Migliorini
- Gorlaeus Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden,The Netherlands
| | - Francesco Nattino
- 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
|
9
|
Nattino F, Migliorini D, Bonfanti M, Kroes GJ. Methane dissociation on Pt(111): Searching for a specific reaction parameter density functional. J Chem Phys 2016; 144:044702. [DOI: 10.1063/1.4939520] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [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
| | - Davide Migliorini
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Matteo Bonfanti
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
10
|
Allouche A. First principles calculations on nitrogen reactivity on tungsten surfaces. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:015001. [PMID: 26569510 DOI: 10.1088/0953-8984/28/1/015001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Using the spin-polarized gradient-corrected density functional theory we investigate the adsorption and dissociation of nitrogen molecule and hydrogen nitride radicals on two model surfaces of tungsten, W{1 0 0} and W{1 1 0}. The goal of the investigations is to predict the number and nature of the nitrogenized moieties that could be found on a tungsten surface after nitrogen-seeded hydrogen (deuterium) plasma discharge. The surfaces are considered clean or saturated in hydrogen or nitrogen atoms as they are expected to be after nitrogen seeded hydrogen plasma irradiation. We find that the radicals NH, NH2 and NH3 are dissociated by the catalytic action of the metallic surfaces if they are animated with an initial kinetic energy of less than 2.5 eV. The products of the reaction are a nitrogen atom incorporated into the metal surface and an adsorbed proton. The general conclusion is that nitrogen reduces hydrogen retention on the surface and more generally that the nitrided surfaces should be less reactive than the clean tungsten even if the two surfaces we consider behave quite differently.
Collapse
Affiliation(s)
- A Allouche
- Physique des Interactions Ioniques et Moléculaires, CNRS and Aix-Marseille Université (UMR7345), Campus Scientifique de Saint Jérôme, service 242, 13397 Marseille Cedex 20-, France
| |
Collapse
|
11
|
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
|
12
|
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
|
13
|
Nattino F, Costanzo F, Kroes GJ. N2 dissociation on W(110): An ab initio molecular dynamics study on the effect of phonons. J Chem Phys 2015; 142:104702. [DOI: 10.1063/1.4913979] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Francesco Nattino
- Gorlaeus Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Francesca Costanzo
- 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
|