1
|
Wagner J, Sibener SJ. Coverage-Dependent Site-Specific Placement and Correlated Diffusion of Atomic Oxygen on Moiré-Patterned Graphene on Ru(0001). J Phys Chem Lett 2024:2936-2943. [PMID: 38451507 DOI: 10.1021/acs.jpclett.4c00274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
Nano-periodic arrays of atomic oxygen are visualized on epitaxial graphene on Ru(0001) via STM following supersonic beam exposure to non-equilibrium fluxes of atomic oxygen. Self-organization of atomic oxygen on graphene is directed by the intrinsic moiré pattern of the ruthenium-graphene interface. Atom-resolved STM imaging reveals the richness of multiparticle interactions, leading to correlated atomic diffusion and placement. Pair-distribution functions demonstrate that repulsive oxygen-oxygen interactions play an increasingly important role in the site specificity and diffusivity of atomic oxygen on the moiré lattice with increasing coverage. Atomic visualization shows the number of oxygen atoms in a local region changes overall diffusion rates and promotes the correlated motion of oxygen atoms. Understanding the site specificity of oxygen adsorption and diffusive behavior of atomic oxygen on epitaxial graphene on Ru(0001) provides insight for both the synthesis and stability of moiré-templated two-dimensional materials which show promise as platforms for next-generation quantum materials and catalysts.
Collapse
Affiliation(s)
- Joshua Wagner
- The James Franck Institute and Department of Chemistry, The University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, United States
| | - Steven J Sibener
- The James Franck Institute and Department of Chemistry, The University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, United States
| |
Collapse
|
2
|
Gerrits N, Jackson B, Bogaerts A. Accurate Reaction Probabilities for Translational Energies on Both Sides of the Barrier of Dissociative Chemisorption on Metal Surfaces. J Phys Chem Lett 2024; 15:2566-2572. [PMID: 38416779 PMCID: PMC10926167 DOI: 10.1021/acs.jpclett.3c03408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Revised: 02/16/2024] [Accepted: 02/26/2024] [Indexed: 03/01/2024]
Abstract
Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule-metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.
Collapse
Affiliation(s)
- Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, Post Office
Box 9502, 2300 RA Leiden, Netherlands
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, BE-2610, Wilrijk, Antwerp, Belgium
| | - Bret Jackson
- Department
of Chemistry, University of Massachusetts
Amherst, Amherst, Massachusetts 01003, United States
| | - Annemie Bogaerts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, BE-2610, Wilrijk, Antwerp, Belgium
| |
Collapse
|
3
|
Wang J, Wang GC. Mechanisms of CH 4 activation over oxygen-preadsorbed transition metals by ReaxFF and AIMD simulations. J Comput Chem 2024; 45:238-246. [PMID: 37746925 DOI: 10.1002/jcc.27233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Revised: 08/26/2023] [Accepted: 09/09/2023] [Indexed: 09/26/2023]
Abstract
The chemisorbed oxygen usually promotes the CH bond activation over less active metals like IB group metals but has no effect or even an inhibition effect over more active metals like Pd based on the static electronic structure study. However, the understanding in terms of dynamics knowledge is far from complete. In the present work, methane dissociation on the oxygen-preadsorbed transition metals including Au, Cu, Ni, Pt, and Pd is systemically studied by reactive force field (ReaxFF). The ReaxFF simulation results indicate that CH4 molecules mainly undergo the direct dissociation on Ni, Pt, and Pd surfaces, while undergo the oxygen-assisted dissociation on Au and Cu surfaces. Additionally, the ab initio molecular dynamics (AIMD) simulations with the umbrella sampling are employed to study the free-energy changes of CH4 dissociation, and the results further support the CH4 dissociation pathway during the ReaxFF simulations. The present results based on ReaxFF and AIMD will provide a deeper dynamic understanding of the effects of pre-adsorbed oxygen species on the CH bond activation compared to that of static DFT.
Collapse
Affiliation(s)
- Jie Wang
- Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education) and the Tianjin key Lab and Molecule-based Material Chemistry, College of Chemistry, Nankai University, Tianjin, China
| | - Gui-Chang Wang
- Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education) and the Tianjin key Lab and Molecule-based Material Chemistry, College of Chemistry, Nankai University, Tianjin, China
| |
Collapse
|
4
|
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
|
5
|
Tchakoua T, Jansen T, van Nies Y, van den Elshout RFA, van Boxmeer BAB, Poort SP, Ackermans MG, Beltrão GS, Hildebrand SA, Beekman SEJ, van der Drift T, Kaart S, Šantić A, Spuijbroek EE, Gerrits N, Somers MF, Kroes GJ. Constructing Mixed Density Functionals for Describing Dissociative Chemisorption on Metal Surfaces: Basic Principles. J Phys Chem A 2023; 127:10481-10498. [PMID: 38051300 PMCID: PMC10726370 DOI: 10.1021/acs.jpca.3c01932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
The production of a majority of chemicals involves heterogeneous catalysis at some stage, and the rates of many heterogeneously catalyzed processes are governed by transition states for dissociative chemisorption on metals. Accurate values of barrier heights for dissociative chemisorption on metals are therefore important to benchmarking electronic structure theory in general and density functionals in particular. Such accurate barriers can be obtained using the semiempirical specific reaction parameter (SRP) approach to density functional theory. However, this approach has thus far been rather ad hoc in its choice of the generic expression of the SRP functional to be used, and there is a need for better heuristic approaches to determining the mixing parameters contained in such expressions. Here we address these two issues. We investigate the ability of several mixed, parametrized density functional expressions combining exchange at the generalized gradient approximation (GGA) level with either GGA or nonlocal correlation to reproduce barrier heights for dissociative chemisorption on metal surfaces. For this, seven expressions of such mixed density functionals are tested on a database consisting of results for 16 systems taken from a recently published slightly larger database called SBH17. Three expressions are derived that exhibit high tunability and use correlation functionals that are either of the PBE GGA form or of one of two limiting nonlocal forms also describing the attractive van der Waals interaction in an approximate way. We also find that, for mixed density functionals incorporating GGA correlation, the optimum fraction of repulsive RPBE GGA exchange obtained with a specific GGA density functional is correlated with the charge-transfer parameter, which is equal to the difference in the work function of the metal surface and the electron affinity of the molecule. However, the correlation is generally not large and not large enough to obtain accurate guesses of the mixing parameter for the systems considered, suggesting that it does not give rise to a very effective search strategy.
Collapse
Affiliation(s)
- Théophile Tchakoua
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Tim Jansen
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Youri van Nies
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | | | - Bart A B van Boxmeer
- Faculty of Applied Engineering, University of Antwerp, Groenenborgerlaan 171, 2020 Antwerpen, Belgium
| | - Saskia P Poort
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Michelle G Ackermans
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Gabriel Spiller Beltrão
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Stefan A Hildebrand
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Steijn E J Beekman
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Thijs van der Drift
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Sam Kaart
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Anthonie Šantić
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Esmee E Spuijbroek
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, 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
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
6
|
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
|
7
|
Jackson B. An examination of phonon-inelastic molecule-metal scattering using reduced density matrix and stochastic wave packet methods. J Chem Phys 2023; 158:024701. [PMID: 36641393 DOI: 10.1063/5.0133638] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
We explore the application of reduced density matrix-based approaches to molecules interacting with the lattice vibrations of metals, an interaction responsible for the temperature dependence of many of the fundamental steps of catalysis. We avoid the use of simple models for the bath and instead use density functional theory to compute all molecule-phonon interactions and the properties of the lattice phonons, for methane scattering from Ir(111). We find that while the large metal mass leads to long bath correlation times, these are not significantly longer than the time over which the reduced density matrix changes due to interactions with the bath. We show that the neglect of memory is reasonable and the use of the Redfield equation is justified. We also show how the commonly used rotating wave approximation is far too severe for this scattering problem. A less restrictive approximation that is nearly exact for our system gives an equation of motion in the Lindblad form. As a result, the Monte Carlo wave packet methods can be used to describe gas-phonon scattering, guaranteeing positivity, and with all couplings derived from first-principles.
Collapse
Affiliation(s)
- Bret Jackson
- Department of Chemistry, University of Massachusetts Amherst, Amherst, Massachusetts 01003, USA
| |
Collapse
|
8
|
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
|
9
|
Hertl N, Krüger K, Bünermann O. Electronically Nonadiabatic H Atom Scattering from Low Miller Index Surfaces of Silver. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:14162-14171. [PMID: 36350763 DOI: 10.1021/acs.langmuir.2c02140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The reactivity of a surface depends strongly on the surface structure. To study the influence of surface structure on H atom adsorption, we performed inelastic scattering experiments and complementary electronically nonadiabatic molecular dynamics (MD) simulations for H atoms colliding with the three low Miller index surface facets of silver. Experiment reveals very similar energy loss distributions for all three investigated facets. However, for the (100) facet a dependence on the surface orientation is observed that is absent for the other two facets. The nonadiabatic MD simulations manage to describe the experiments well. Despite the observed insignificant influence of the surface geometry on the energy loss distributions, our simulations predict that the capability of the H atoms to penetrate the surface critically depends on the surface structure. The observed crystal orientation dependence of the energy loss distributions in the experiment for Ag(100) cannot be explained with our simulations, and we provide a discussion for a better theoretical description of this system to stimulate future computational investigations.
Collapse
Affiliation(s)
- Nils Hertl
- Max-Planck-Institut für Multidisziplinäre Naturwissenschaften, Am Faßberg 11, 37077Göttingen, Germany
| | - Kerstin Krüger
- Institut für Physikalische Chemie, Georg-August-Universität, Tammannstraße 6, 37073Göttingen, Germany
| | - Oliver Bünermann
- Institut für Physikalische Chemie, Georg-August-Universität, Tammannstraße 6, 37073Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, 37077Göttingen, Germany
| |
Collapse
|
10
|
Meng Q, Chen J, Ma J, Zhang X, Chen J. Adiabatic models for the quantum dynamics of surface scattering with lattice effects. Phys Chem Chem Phys 2022; 24:16415-16436. [PMID: 35766107 DOI: 10.1039/d2cp01560a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Junbo Chen
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China. .,Xi'an Modern Chemistry Research Institute, China North Industries Group Corp., Ltd., East Zhangba Road 168, 710065 Xi'an, China
| | - Jianxing Ma
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Jun Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Yangqiao Road West 155, 350002 Fuzhou, China.,Fujian Science & Technology Innovation Laboratory for Optoelectronic Information of China, Optoelectronic Industry Base at High-tech Zone, 350108 Fuzhou, China
| |
Collapse
|
11
|
Jackson B. Quantum studies of methane-metal inelastic diffraction and trapping: the variation with molecular orientation and phonon coupling. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111516] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
12
|
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
|
13
|
Guo H, Poths P, Sautet P, Alexandrova AN. Oxidation Dynamics of Supported Catalytic Cu Clusters: Coupling to Fluxionality. ACS Catal 2021. [DOI: 10.1021/acscatal.1c04100] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Han Guo
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Patricia Poths
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Philippe Sautet
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
- California NanoSystems Institute, Los Angeles, California 90095, United States
| | - Anastassia N. Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
- California NanoSystems Institute, Los Angeles, California 90095, United States
| |
Collapse
|
14
|
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
|
15
|
Putra SEM, Morikawa Y, Hamada I. Isotope effect of methane adsorbed on fcc metal (1 1 1) surfaces. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138943] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
16
|
Gerrits N, Geweke J, Auerbach DJ, Beck RD, Kroes GJ. Highly Efficient Activation of HCl Dissociation on Au(111) via Rotational Preexcitation. J Phys Chem Lett 2021; 12:7252-7260. [PMID: 34313445 PMCID: PMC8350909 DOI: 10.1021/acs.jpclett.1c02093] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 07/23/2021] [Indexed: 06/13/2023]
Abstract
The probability for dissociation of molecules on metal surfaces, which often controls the rate of industrially important catalytic processes, can depend strongly on how energy is partitioned in the incident molecule. There are many example systems where the addition of vibrational energy promotes reaction more effectively than the addition of translational energy, but for rotational pre-excitation similar examples have not yet been discovered. Here, we make an experimentally testable theoretical prediction that adding energy to the rotation of HCl can promote its dissociation on Au(111) 20 times more effectively than increasing its translational energy. In the underlying mechanism, the molecule's initial rotational motion allows it to pass through a critical region of the reaction path, where this path shows a strong and nonmonotonic dependence on the molecular orientation.
Collapse
Affiliation(s)
- Nick Gerrits
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jan Geweke
- Department
of Dynamics at Surfaces, Max Planck Institute
for Biophysical Chemistry, Göttingen, Am Fassberg 11, 37077 Göttingen, Germany
- Institute
for Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Daniel J. Auerbach
- Department
of Dynamics at Surfaces, Max Planck Institute
for Biophysical Chemistry, Göttingen, Am Fassberg 11, 37077 Göttingen, Germany
| | - Rainer D. Beck
- Laboratoire
de Chimie Physique Moléculaire, École
Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Geert-Jan Kroes
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
17
|
Jackson B. The trapping of methane on Ir(111): A first-principles quantum study. J Chem Phys 2021; 155:044705. [PMID: 34340380 DOI: 10.1063/5.0058672] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We implement a fully quantum mechanical study of methane trapping on Ir(111), where the phonons, the molecule-surface interaction, and the molecule-phonon coupling are all computed from first-principles. We find that both the surface corrugation and the phonon coupling vary strongly with molecular orientation and that there is a "chemical" aspect to this due to the catalytic nature of the metal. For example, molecules with reactive orientations can approach close to surface sites with low barriers to dissociation. Moreover, lattice motion can lower the barrier to dissociation, leading to unusual behavior for the phonon coupling. We find good agreement with experiment and two recent classical studies if we average our potential energy surface over several orientations of the molecule. We also find reasonable agreement with a recent study of methane diffraction, although we show that diffraction does not play a major role in trapping on the smooth Ir(111) surface and that trapping obeys normal energy scaling, consistent with experiment. We show that the trapping probability can be sensitive to the temperature at both high and low incidence energies. Relaxation and desorption of trapped particles are examined.
Collapse
Affiliation(s)
- Bret Jackson
- Department of Chemistry, University of Massachusetts Amherst, Amherst, Massachusetts 01003, USA
| |
Collapse
|
18
|
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
|
19
|
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
|
20
|
Zhou X, Zhang Y, Guo H, Jiang B. Towards bridging the structure gap in heterogeneous catalysis: the impact of defects in dissociative chemisorption of methane on Ir surfaces. Phys Chem Chem Phys 2021; 23:4376-4385. [DOI: 10.1039/d0cp06535h] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The negatively activated region in CH4 dissociation is attributed to a precursor-mediated mechanism involving surface defects.
Collapse
Affiliation(s)
- Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Hua Guo
- Department of Chemistry and Chemical Biology
- University of New Mexico
- Albuquerque
- USA
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| |
Collapse
|
21
|
Gerrits N, Smeets EWF, Vuckovic S, Powell AD, Doblhoff-Dier K, Kroes GJ. Density Functional Theory for Molecule-Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not. J Phys Chem Lett 2020; 11:10552-10560. [PMID: 33295770 PMCID: PMC7751010 DOI: 10.1021/acs.jpclett.0c02452] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
While density functional theory (DFT) is perhaps the most used electronic structure theory in chemistry, many of its practical aspects remain poorly understood. For instance, DFT at the generalized gradient approximation (GGA) tends to fail miserably at describing gas-phase reaction barriers, while it performs surprisingly well for many molecule-metal surface reactions. GGA-DFT also fails for many systems in the latter category, and up to now it has not been clear when one may expect it to work. We show that GGA-DFT tends to work if the difference between the work function of the metal and the molecule's electron affinity is greater than ∼7 eV and to fail if this difference is smaller, with sticking of O2 on Al(111) being a spectacular example. Using dynamics calculations we show that, for this system, the DFT problem may be solved as done for gas-phase reactions, i.e., by resorting to hybrid functionals, but using screening at long-range to obtain a correct description of the metal. Our results suggest the GGA error in the O2 + Al(111) barrier height to be functional driven. Our results also suggest the possibility to compute potential energy surfaces for the difficult-to-treat systems with computationally cheap nonself-consistent calculations in which a hybrid functional is applied to a GGA density.
Collapse
Affiliation(s)
- Nick Gerrits
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Egidius W. F. Smeets
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Stefan Vuckovic
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Andrew D. Powell
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
22
|
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
|
23
|
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
|
24
|
Zhang L, Jiang B. State-to-state quantum dynamics of H2O/HOD scattering from Cu(111): Mode- and bond-selective vibrational energy transfer. J Chem Phys 2020; 153:214702. [DOI: 10.1063/5.0030490] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Liang Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
25
|
Ding X, Zhu H, Ren H, Liu D, Yu Z, Shi N, Guo W. Adsorption and dehydrogenation of C 2-C 6n-alkanes over a Pt catalyst: a theoretical study on the size effects of alkane molecules and Pt substrates. Phys Chem Chem Phys 2020; 22:21835-21843. [PMID: 32966439 DOI: 10.1039/d0cp03194a] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Adsorption and dehydrogenation of C2-C6n-alkanes are investigated on a Pt substrate using density functional theory (DFT) calculations, and the size effects of alkane molecules and Pt substrates are discussed in detail. The Pt(111) surface and Pt55 cluster are chosen to represent large and small Pt nanoparticles, respectively. The C2-C6 straight-chain alkanes show no site preference on Pt(111) drifting over the surface, but prefer to locate along the edge sites of Pt55. Our results suggest that a linear relationship holds for the adsorption energies of n-alkanes against the chain length on Pt(111), in accordance with the experimental observations. Pt55 also exhibits a similar linear relationship for n-alkanes but with larger adsorption energies due to the low-coordinated Pt atoms at the edge site. For the two-step dehydrogenation from alkanes to alkenes, the first dehydrogenation reaction is the rate-determining step (RDS) on Pt(111), and a larger size of alkane molecule will lead to a lower dehydrogenation activity. While on Pt55, no RDS is present and the dehydrogenation activity oscillates slightly as the chain length of n-alkane increases. Generally, Pt55 involves lower energy barriers for most dehydrogenation steps compared to Pt(111), indicating that small Pt particles with more low-coordinated Pt atoms are more active towards alkane dehydrogenation. In addition, a clear BEP relationship is identified for all the dehydrogenation reactions of C2-C6n-alkanes on Pt substrates, and this linear relationship is independent of the particle size of the Pt substrate and the chain length of alkanes.
Collapse
Affiliation(s)
- Xuefei Ding
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, P. R. China.
| | - Houyu Zhu
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, P. R. China.
| | - Hao Ren
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, P. R. China.
| | - Dongyuan Liu
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, P. R. China.
| | - Zehua Yu
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, P. R. China.
| | - Naiyou Shi
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, P. R. China.
| | - Wenyue Guo
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, P. R. China.
| |
Collapse
|
26
|
Jackson B. Direct and trapping-mediated pathways to dissociative chemisorption: CH4 dissociation on Ir(111) with step defects. J Chem Phys 2020; 153:034704. [DOI: 10.1063/5.0012252] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Bret Jackson
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, USA
| |
Collapse
|
27
|
Roy S, K. J. N, Tiwari N, Tiwari AK. Energetics and dynamics of CH4 and H2O dissociation on metal surfaces. INT REV PHYS CHEM 2020. [DOI: 10.1080/0144235x.2020.1765598] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Sudipta Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| | - Nayanthara K. J.
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| | - Nidhi Tiwari
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| | - Ashwani K. Tiwari
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| |
Collapse
|
28
|
Guo H, Sautet P, Alexandrova AN. Reagent-Triggered Isomerization of Fluxional Cluster Catalyst via Dynamic Coupling. J Phys Chem Lett 2020; 11:3089-3094. [PMID: 32227852 DOI: 10.1021/acs.jpclett.0c00548] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Metallic cluster catalysts have many thermodynamically accessible isomers with diverse active sites and low reaction barriers, and lately a strong hypothesis emerged that the many catalyst states collectively drive the catalysis. However, it remained a hypothesis that catalyst isomerization is actually kinetically feasible under the current reaction conditions. Using high-temperature dynamics simulations and sampling, a range of orientations, and vibrational energy distributions, we probe how thermal effects and molecular events affect cluster catalyst dynamics. We show that even such a delicate affair as the dissociation or scattering of a methane molecule on the heavy and thus slow Pt13 cluster triggers substantial isomerization of the catalyst, far beyond thermal at 700 K. A kinetic coupling between the methane activity and cluster catalyst dynamics is observed. In return, the thermal dynamics of the cluster affects the methane reaction and scattering probabilities. Hence, molecular events at the surfaces of fluxional cluster catalysts should facilitate the population of an ensemble of catalyst states under the current reaction conditions, with implications for available active sites, reaction mechanisms, and apparent rates.
Collapse
Affiliation(s)
| | - Philippe Sautet
- California NanoSystems Institute, Los Angeles, California 90095, United States
| | | |
Collapse
|
29
|
Gutiérrez-González A, Beck RD. Quantum state and surface-site-resolved studies of methane chemisorption by vibrational spectroscopies. Phys Chem Chem Phys 2020; 22:17448-17459. [DOI: 10.1039/d0cp03134h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Infrared spectroscopic methods enable quantum-state-specific and surface-site-selective studies of methane chemisorption on stepped platinum surfaces.
Collapse
Affiliation(s)
- Ana Gutiérrez-González
- Laboratoire de Chimie Physique Moléculaire (LCPM)
- École Polytechnique Fédérale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| | - Rainer D. Beck
- Laboratoire de Chimie Physique Moléculaire (LCPM)
- École Polytechnique Fédérale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| |
Collapse
|
30
|
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
|
31
|
Liu T, Fu B, Zhang DH. Six-dimensional potential energy surfaces for the dissociative chemisorption of HCl on rigid Ag(100) and Ag(110) surfaces. J Chem Phys 2019; 151:144707. [DOI: 10.1063/1.5122218] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tianhui Liu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| |
Collapse
|
32
|
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
|
33
|
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
|
34
|
Guo H, Jackson B. Methane dissociation on stepped Ni surfaces resolved by impact site, collision energy, vibrational state, and lattice distortion. J Chem Phys 2019; 150:204703. [DOI: 10.1063/1.5095145] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Han Guo
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Bret Jackson
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, USA
| |
Collapse
|
35
|
Jiang B, Guo H. Dynamics in reactions on metal surfaces: A theoretical perspective. J Chem Phys 2019; 150:180901. [PMID: 31091904 DOI: 10.1063/1.5096869] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.
Collapse
Affiliation(s)
- Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
36
|
Ghassemi E, Somers MF, Kroes GJ. Assessment of Two Problems of Specific Reaction Parameter Density Functional Theory: Sticking and Diffraction of H 2 on Pt(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:10406-10418. [PMID: 31049122 PMCID: PMC6488140 DOI: 10.1021/acs.jpcc.9b00981] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/31/2019] [Revised: 03/26/2019] [Indexed: 06/09/2023]
Abstract
It is important that theory is able to accurately describe dissociative chemisorption reactions on metal surfaces, as such reactions are often rate-controlling in heterogeneously catalyzed processes. Chemically accurate theoretical descriptions have recently been obtained on the basis of the specific reaction parameter (SRP) approach to density functional (DF) theory (DFT), allowing reaction barriers to be obtained with chemical accuracy. However, being semiempirical, this approach suffers from two basic problems. The first is that sticking probabilities (to which SRP density functionals (DFs) are usually fitted) might show differences across experiments, of which the origins are not always clear. The second is that it has proven hard to use experiments on diffractive scattering of H2 from metals for validation purposes, as dynamics calculations using a SRP-DF may yield a rather poor description of the measured data, especially if the potential used contains a van der Waals well. We address the first problem by performing dynamics calculations on three sets of molecular beam experiments on D2 + Pt(111), using four sets of molecular beam parameters to obtain sticking probabilities, and the SRP-DF recently fitted to one set of experiments on D2 + Pt(111). It is possible to reproduce all three sets of experiments with chemical accuracy with the aid of two sets of molecular beam parameters. The theoretical simulations with the four different sets of beam parameters allow one to determine for which range of incidence conditions the experiments should agree well and for which conditions they should show specific differences. This allows one to arrive at conclusions about the quality of the experiments and about problems that might affect the experiments. Our calculations on diffraction of H2 scattering from Pt(111) show both quantitative and qualitative differences with previously measured diffraction probabilities, which were Debye-Waller (DW)-extrapolated to 0 K. We suggest that DW extrapolation, which is appropriate for direct scattering, might fail if the scattering is affected by the presence of a van der Waals well and that theory should attempt to model surface atom motion for reproducing diffraction experiments performed for surface temperatures of 500 K and higher.
Collapse
|
37
|
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
|
38
|
|
39
|
Al Taleb A, Anemone G, Zhou L, Guo H, Farías D. Diffraction of CH 4 from a Metal Surface. J Phys Chem Lett 2019; 10:1574-1580. [PMID: 30855971 DOI: 10.1021/acs.jpclett.9b00158] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Diffraction with matter waves has been reported since the beginning of quantum mechanics. In free space, diffraction effects have been observed even with objects as large as C60 molecules. However, in scattering from a solid surface, pure elastic diffraction features have never been observed with molecules larger than D2. Here we report the observation of pure molecular diffraction for CH4 scattered off of an Ir(111) surface. These results prove that quantum coherence is preserved, despite the small separation between rotational levels and the interaction with surface phonons. Density functional theory calculations of the potential energy surface provide some clues to understand the larger corrugation sampled by CH4 molecules in comparison to Ne atoms. Our results show that isotope separation of polyatomic molecules may be possible using gas-surface diffraction.
Collapse
Affiliation(s)
- Amjad Al Taleb
- Departamento de Física de la Materia Condensada , Universidad Autónoma de Madrid , 28049 Madrid , Spain
- Instituto Madrileño de Estudios Avanzados en Nanociencia (IMDEA-Nanociencia) , Cantoblanco, 28049 Madrid , Spain
| | - Gloria Anemone
- Departamento de Física de la Materia Condensada , Universidad Autónoma de Madrid , 28049 Madrid , Spain
| | - Linsen Zhou
- Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , United States
| | - Daniel Farías
- Departamento de Física de la Materia Condensada , Universidad Autónoma de Madrid , 28049 Madrid , Spain
- Instituto "Nicolás Cabrera" , Universidad Autónoma de Madrid , 28049 Madrid , Spain
- Condensed Matter Physics Center (IFIMAC) , Universidad Autónoma de Madrid , 28049 Madrid , Spain
| |
Collapse
|
40
|
Chadwick H, Gutiérrez-González A, Beck RD, Kroes GJ. Transferability of the SRP32-vdW specific reaction parameter functional to CHD3 dissociation on Pt(110)-(2 × 1). J Chem Phys 2019; 150:124702. [DOI: 10.1063/1.5081005] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Helen Chadwick
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Ana Gutiérrez-González
- Laboratoire de Chimie Physique Moléculaire, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Rainer D. Beck
- Laboratoire de Chimie Physique Moléculaire, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
41
|
Zhang Y, Zhou X, Jiang B. Bridging the Gap between Direct Dynamics and Globally Accurate Reactive Potential Energy Surfaces Using Neural Networks. J Phys Chem Lett 2019; 10:1185-1191. [PMID: 30802067 DOI: 10.1021/acs.jpclett.9b00085] [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/18/2023]
Abstract
Direct dynamics simulations become increasingly popular in studying reaction dynamics for complex systems where analytical potential energy surfaces (PESs) are unavailable. Yet, the number and/or the propagation time of trajectories are often limited by high computational costs, and numerous energies and forces generated on-the-fly become wasted after simulations. We demonstrate here an example of reusing only a very small portion of existing direct dynamics data to reconstruct a 90-dimensional globally accurate reactive PES describing the interaction of CO2 with a movable Ni(100) surface based on a machine learning approach. In addition to reproducing previous results with much better statistics, we predict scattering probabilities of CO2 at the state-to-state level, which is extremely demanding for direct dynamics. We propose this unified way to investigate gaseous and gas-surface reactions of medium size, initiating with hundreds of preliminary direct dynamics trajectories, followed by low-cost and high-quality simulations on full-dimensional analytical PESs.
Collapse
Affiliation(s)
- Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes , University of Science and Technology of China , Hefei , Anhui 230026 , China
| | - Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes , University of Science and Technology of China , Hefei , Anhui 230026 , China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes , University of Science and Technology of China , Hefei , Anhui 230026 , China
| |
Collapse
|
42
|
Ghassemi EN, Smeets EWF, Somers MF, Kroes GJ, Groot IMN, Juurlink LBF, Füchsel G. Transferability of the Specific Reaction Parameter Density Functional for H 2 + Pt(111) to H 2 + Pt(211). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:2973-2986. [PMID: 30792827 PMCID: PMC6376921 DOI: 10.1021/acs.jpcc.8b11018] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2018] [Revised: 12/20/2018] [Indexed: 05/25/2023]
Abstract
The accurate description of heterogeneously catalyzed reactions may require chemically accurate evaluation of barriers for reactions of molecules at the edges of metal nanoparticles. It was recently shown that a semiempirical density functional describing the interaction of a molecule dissociating on a flat metal surface (CHD3 + Pt(111)) is transferable to the same molecule reacting on a stepped surface of the same metal (Pt(211)). However, validation of the method for additional systems is desirable. To address the question whether the specific reaction parameter (SRP) functional that describes H2 + Pt(111) with chemical accuracy is also capable of accurately describing H2 + Pt(211), we have performed molecular beam simulations with the quasi-classical trajectory (QCT) method, using the SRP functional developed for H2 + Pt(111). Our calculations used the Born-Oppenheimer static surface model. The accuracy of the QCT method was assessed by comparison with quantum dynamics results for reaction of the ro-vibrational ground state of H2. The theoretical results for sticking of H2 and D2 on Pt(211) are in quite good agreement with the experiment, but uncertainties remain because of a lack of accuracy of the QCT simulations at low incidence energies and possible inaccuracies in the reported experimental incidence energies at high energies. We also investigated the nonadiabatic effect of electron-hole pair excitation on the reactivity using the molecular dynamics with the electron friction (MDEF) method, employing the local density friction approximation (LDFA). Only small effects of electron-hole pair excitation on sticking are found.
Collapse
Affiliation(s)
- Elham Nour Ghassemi
- Leiden Institute of Chemistry, Gorlaeus Laboratories,
Leiden University, P.O. Box 9502, 2300 RA Leiden,
The Netherlands
| | - Egidius W. F. Smeets
- 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
| | - Irene M. N. Groot
- Leiden Institute of Chemistry, Gorlaeus Laboratories,
Leiden University, P.O. Box 9502, 2300 RA Leiden,
The Netherlands
| | - Ludo B. F. Juurlink
- Leiden Institute of Chemistry, Gorlaeus Laboratories,
Leiden University, P.O. Box 9502, 2300 RA Leiden,
The Netherlands
| | - Gernot Füchsel
- Leiden Institute of Chemistry, Gorlaeus Laboratories,
Leiden University, P.O. Box 9502, 2300 RA Leiden,
The Netherlands
- Institut für Chemie und
Biochemie—Physikalische und Theoretische Chemie, Freie
Universität Berlin, Takustraße 3, 14195 Berlin,
Germany
| |
Collapse
|
43
|
Füchsel G, Zhou X, Jiang B, Juaristi JI, Alducin M, Guo H, Kroes GJ. Reactive and Nonreactive Scattering of HCl from Au(111): An Ab Initio Molecular Dynamics Study. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:2287-2299. [PMID: 30740194 PMCID: PMC6366682 DOI: 10.1021/acs.jpcc.8b10686] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 12/19/2018] [Indexed: 05/20/2023]
Abstract
The HCl + Au(111) system has recently become a benchmark for highly activated dissociative chemisorption, which presumably is strongly affected by electron-hole pair excitation. Previous dynamics calculations, which were based on density functional theory at the generalized gradient approximation level (GGA-DFT) for the molecule-surface interaction, have all overestimated measured reaction probabilities by at least an order of magnitude. Here, we perform ab initio molecular dynamics (AIMD) and AIMD with electronic friction (AIMDEF) calculations employing a density functional that includes the attractive van der Waals interaction. Our calculations model the simultaneous and possibly synergistic effects of surface temperature, surface atom motion, electron-hole pair excitation, the molecular beam conditions of the experiments, and the van der Waals interaction on the reactivity. We find that reaction probabilities computed with AIMDEF and the SRP32-vdW functional still overestimate the measured reaction probabilities, by a factor 18 for the highest incidence energy at which measurements were performed (≈2.5 eV). Even granting that the experiment could have underestimated the sticking probability by about a factor three, this still translates into a considerable overestimation of the reactivity by the current theory. Likewise, scaled transition probabilities for vibrational excitation from ν = 1, j = 1 to ν = 2 are overestimated by the AIMDEF theory, by factors 3-8 depending on the initial conditions modeled. Energy losses to the surface and translational energy losses are, however, in good agreement with experimental values.
Collapse
Affiliation(s)
- Gernot Füchsel
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- Institut
für Chemie und Biochemie—Physikalische und Theoretische
Chemie, Freie Universität Berlin, Takustraße3, 14195 Berlin, Germany
- E-mail: (G.F.)
| | - Xueyao Zhou
- Hefei
National Laboratory for Physical Science at the Microscale, Department
of Chemical Physics, School of Chemistry and Materials, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Bin Jiang
- Hefei
National Laboratory for Physical Science at the Microscale, Department
of Chemical Physics, School of Chemistry and Materials, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - J. Iñaki Juaristi
- Departamento
de Física de Materiales, Facultad
de Químicas (UPV/EHU), Apartado 1072, 20080 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
- Donostia
International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - 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
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University
of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- E-mail: . Phone: +31 (0)71 527
4396 (G.-J.K.)
| |
Collapse
|
44
|
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
|
45
|
Guo H, Menzel JP, Jackson B. Quantum dynamics studies of the dissociative chemisorption of CH4 on the steps and terraces of Ni(211). J Chem Phys 2018; 149:244704. [DOI: 10.1063/1.5066553] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Affiliation(s)
- Han Guo
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Jan Paul Menzel
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Bret Jackson
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, USA
| |
Collapse
|
46
|
Migliorini D, Nattino F, Tiwari AK, Kroes GJ. HOD on Ni(111): Ab Initio molecular dynamics prediction of molecular beam experiments. J Chem Phys 2018; 149:244706. [DOI: 10.1063/1.5059357] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Davide Migliorini
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Francesco Nattino
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Ashwani K. Tiwari
- Indian Institute of Science Education and Research Kolkata, Mohanpur 741246 West Bengal, India
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
47
|
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
|
48
|
Liu T, Fu B, Zhang DH. Six-dimensional quantum dynamics for the dissociative chemisorption of HCl on rigid Ag(111) on three potential energy surfaces with different density functionals. J Chem Phys 2018; 149:174702. [DOI: 10.1063/1.5053827] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tianhui Liu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| |
Collapse
|
49
|
Nour Ghassemi E, Somers M, Kroes GJ. Test of the Transferability of the Specific Reaction Parameter Functional for H 2 + Cu(111) to D 2 + Ag(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2018; 122:22939-22952. [PMID: 30344838 PMCID: PMC6189907 DOI: 10.1021/acs.jpcc.8b05658] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Revised: 09/04/2018] [Indexed: 06/08/2023]
Abstract
The accurate description of the dissociative chemisorption of a molecule on a metal surface requires a chemically accurate description of the molecule-surface interaction. Previously, it was shown that the specific reaction parameter approach to density functional theory (SRP-DFT) enables accurate descriptions of the reaction of dihydrogen with metal surfaces in, for instance, H2 + Pt(111), H2 + Cu(111), and H2 + Cu(100). SRP-DFT likewise allowed a chemically accurate description of dissociation of methane on Ni(111) and Pt(111), and the SRP functional for CH4 + Ni(111) was transferable to CH4 + Pt(111), where Ni and Pt belong to the same group. Here, we investigate whether the SRP density functional derived for H2 + Cu(111) also gives chemically accurate results for H2 + Ag(111), where Ag belongs to the same group as Cu. To do this, we have performed quasi-classical trajectory calculations using the six-dimensional potential energy surface of H2 + Ag(111) within the Born-Oppenheimer static surface approximation. The computed reaction probabilities are compared with both state-resolved associative desorption and molecular beam sticking experiments. Our results do not yet show transferability, as the computed sticking probabilities and initial-state selected reaction probabilities are shifted relative to experiment to higher energies by about 2-3 kcal/mol. The lack of transferability may be due to the different character of the SRP functionals for H2 + Cu and CH4 + group 10 metals, the latter containing a van der Waals correlation functional and the former not.
Collapse
|
50
|
Chen J, Zhou X, Zhang Y, Jiang B. Vibrational control of selective bond cleavage in dissociative chemisorption of methanol on Cu(111). Nat Commun 2018; 9:4039. [PMID: 30279479 PMCID: PMC6168487 DOI: 10.1038/s41467-018-06478-6] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2018] [Accepted: 08/24/2018] [Indexed: 11/25/2022] Open
Abstract
Controlling product branching ratios in a chemical reaction represents a desired but difficult achievement in chemistry. In this work, we demonstrate the first example of altering the branching ratios in a multichannel reaction, i.e., methanol dissociative chemisorption on Cu(111), via selectively exciting specific vibrational modes. To this end, we develop a globally accurate full-dimensional potential energy surface for the CH3OH/Cu(111) system and perform extensive vibrational state-selected molecular dynamics simulations. Our results show that O-H/C-H/C-O stretching vibrational excitations substantially enhance the respective bond scission processes, representing extraordinary bond selectivity. At a given total energy, the branching ratio of C-O/C-H dissociation can increase by as large as 100 times by exciting the C-O stretching mode which possesses an unprecedentedly strong vibrational efficacy on reactivity. This vibrational control can be realized by the well-designed experiment using a linearly polarized laser.
Collapse
Affiliation(s)
- Jialu Chen
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China.
| |
Collapse
|