1
|
van Bree RAB, Gerrits N, Kroes GJ. Dissociative chemisorption of O 2 on Al(111): dynamics on a potential energy surface computed with a non-self-consistent screened hybrid density functional approach. Faraday Discuss 2024. [PMID: 38787655 DOI: 10.1039/d3fd00165b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2024]
Abstract
Density functional theory (DFT) at the generalized gradient approximation (GGA) level is often considered the best compromise between feasibility and accuracy for reactions of molecules on metal surfaces. Recent work, however, strongly suggests that density functionals (DFs) based on GGA exchange are not able to describe molecule-metal surface reactions for which the work function of the metal surface minus the electron affinity of the molecule is less than 7 eV. Systems for which this is true exhibit an increased charge transfer from the metal to the molecule at the transition state, increasing the delocalisation of the electron density. This enlarged delocalisation can cause GGA-DFT to underestimate energy values relative to the gas-phase and thus underestimate the barrier height, similar to what has been observed for several gas-phase reactions. An example of such a molecule-metal surface system is O2 + Al(111). Following a similar strategy as for gas-phase reactions, previous work showed results of increased accuracy when using a screened hybrid DF for O2 + Al(111). However, even screened hybrid DFs are computationally expensive to use for metal surfaces. To resolve this, we test a non-self-consistent field (NSCF) screened hybrid DF approach. This approach computes screened hybrid DFT energies based on self-consistent-field (SCF) GGA electronic densities. Here, we explore the accuracy of the NSCF screened hybrid DF approach by implementing the NSCF HSE03-1/3x@RPBE DF for O2 + Al(111). We compute and analyse molecular beam sticking probabilities as well as a set of sticking probabilities for rotationally aligned O2. Our results show that the NSCF approach results in reaction probability curves that reproduce SCF results with near-chemical accuracy, suggesting that the NSCF approach can be used advantageously for exploratory purposes. An analysis of the potential energy surface and the barriers gives insight into the cause of the disagreement between the SCF and NSCF reaction probabilities and into the changes needed in theoretical modelling to further improve the description of the O2 + Al(111) system. Finally, the hole model yields fair agreement with dynamics results for the reaction probability curve, but results in an increased slope of the reaction probability curve compared to the molecular dynamics, with a shift to lower or higher energies depending on whether the vibrational energy of the molecule is included in the initial energy of the molecule or not.
Collapse
Affiliation(s)
- Robert A B van Bree
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| | - Nick Gerrits
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| | - Geert-Jan Kroes
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| |
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
|
Luo W, Yan X, Pan X, Jiao J, Mai L. What Makes On-Chip Microdevices Stand Out in Electrocatalysis? SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2024; 20:e2305020. [PMID: 37875658 DOI: 10.1002/smll.202305020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Revised: 09/03/2023] [Indexed: 10/26/2023]
Abstract
Clean and sustainable energy conversion and storage through electrochemistry shows great promise as an alternative to traditional fuel or fossil-consumption energy systems. With regards to practical and high-efficient electrochemistry application, the rational design of active sites and the accurate description of mechanism remain a challenge. Toward this end, in this Perspective, a unique on-chip micro/nano device coupling nanofabrication and low-dimensional electrochemical materials is presented, in which material structure analysis, field-effect regulation, in situ monitoring, and simulation modeling are highlighted. The critical mechanisms that influence electrochemical response are discussed, and how on-chip micro/nano device distinguishes itself is emphasized. The key challenges and opportunities of on-chip electrochemical platforms are also provided through the Perspective.
Collapse
Affiliation(s)
- Wen Luo
- Department of Physics, School of Science, Wuhan University of Technology, Wuhan, 430070, China
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan, 430070, China
| | - Xin Yan
- Department of Physics, School of Science, Wuhan University of Technology, Wuhan, 430070, China
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan, 430070, China
| | - Xuelei Pan
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan, 430070, China
- Wolfson Catalysis Centre, Department of Chemistry, University of Oxford, Oxford, OX1 3QR, UK
| | - Jinying Jiao
- Department of Physics, School of Science, Wuhan University of Technology, Wuhan, 430070, China
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan, 430070, China
| | - Liqiang Mai
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan, 430070, China
| |
Collapse
|
4
|
Rey J, Chizallet C, Rocca D, Bučko T, Badawi M. Reference-Quality Free Energy Barriers in Catalysis from Machine Learning Thermodynamic Perturbation Theory. Angew Chem Int Ed Engl 2024; 63:e202312392. [PMID: 38055209 DOI: 10.1002/anie.202312392] [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: 08/23/2023] [Revised: 11/11/2023] [Accepted: 12/06/2023] [Indexed: 12/07/2023]
Abstract
For the first time, we report calculations of the free energies of activation of cracking and isomerization reactions of alkenes that combine several different electronic structure methods with molecular dynamics simulations. We demonstrate that the use of a high level of theory (here Random Phase Approximation-RPA) is necessary to bridge the gap between experimental and computed values. These transformations, catalyzed by zeolites and proceeding via cationic intermediates and transition states, are building blocks of many chemical transformations for valorization of long chain paraffins originating, e.g., from plastic waste, vegetable oils, Fischer-Tropsch waxes or crude oils. Compared with the free energy barriers computed at the PBE+D2 production level of theory via constrained ab initio molecular dynamics, the barriers computed at the RPA level by the application of Machine Learning thermodynamic Perturbation Theory (MLPT) show a significant decrease for isomerization reaction and an increase of a similar magnitude for cracking, yielding an unprecedented agreement with the results obtained by experiments and kinetic modeling.
Collapse
Affiliation(s)
- Jérôme Rey
- Laboratoire de Physique et Chimie Théoriques LPCT UMR 7019-CNRS, Université de Lorraine, Vandœuvre-lés-Nancy, France
| | - Céline Chizallet
- IFP Energies nouvelles, Rond-Point de l'Ēchangeur de Solaize, BP3, 69360, Solaize, France
| | - Dario Rocca
- Laboratoire de Physique et Chimie Théoriques LPCT UMR 7019-CNRS, Université de Lorraine, Vandœuvre-lés-Nancy, France
| | - Tomáš Bučko
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, SK-84215, Bratislava, Slovakia
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská cesta 9, SK-84236, Bratislava, Slovakia
| | - Michael Badawi
- Laboratoire de Physique et Chimie Théoriques LPCT UMR 7019-CNRS, Université de Lorraine, Vandœuvre-lés-Nancy, France
- Laboratoire Lorrain de Chimie Moléculaire L2CM UMR 7053-CNRS, Université de Lorraine, Metz, France
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Hasegawa T, Hagiwara S, Otani M, Maeda S. A Combined Reaction Path Search and Hybrid Solvation Method for the Systematic Exploration of Elementary Reactions at the Solid-Liquid Interface. J Phys Chem Lett 2023; 14:8796-8804. [PMID: 37747821 DOI: 10.1021/acs.jpclett.3c02233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
We present a combined simulation method of single-component artificial force induced reaction (SC-AFIR) and effective screening medium combined with the reference interaction site model (ESM-RISM), termed SC-AFIR+ESM-RISM. SC-AFIR automatically and systematically explores the chemical reaction pathway, and ESM-RISM directly simulates the precise electronic structure at the solid-liquid interface. Hence, SC-AFIR+ESM-RISM enables us to explore reliable reaction pathways at the solid-liquid interface. We applied it to explore the dissociation pathway of an H2O molecule at the Cu(111)/water interface. The reaction path networks of the whole reaction and the minimum energy paths from H2O to H2 + O depend on the interfacial environment. The qualitative difference in the energy diagrams and the resulting change in the kinematically favored dissociation pathway upon changing the solvation environments are discussed. We believe that SC-AFIR+ESM-RISM will be a powerful tool to reveal the details of chemical reactions in surface catalysis and electrochemistry.
Collapse
Affiliation(s)
- Taisuke Hasegawa
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10 Nishi 8, Kita-ku, Sapporo 060-0810, Japan
| | - Satoshi Hagiwara
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tenno-dai, Tsukuba 305-8577, Japan
| | - Minoru Otani
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tenno-dai, Tsukuba 305-8577, Japan
| | - Satoshi Maeda
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10 Nishi 8, Kita-ku, Sapporo 060-0810, Japan
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Kita 13, Nishi 8, Sapporo 060-8628 Japan
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21 Nishi 10, Kita-ku, Sapporo 001-0021, Japan
| |
Collapse
|
8
|
Lin Q, Jiang B. Modeling Equilibration Dynamics of Hyperthermal Products of Surface Reactions Using Scalable Neural Network Potential with First-Principles Accuracy. J Phys Chem Lett 2023; 14:7513-7518. [PMID: 37582162 DOI: 10.1021/acs.jpclett.3c01708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/17/2023]
Abstract
Equilibration dynamics of hot oxygen atoms following the dissociation of O2 on Pd(100) and Pd(111) surfaces are investigated by molecular dynamics simulations based on a scalable neural network potential enabling first-principles description of the interaction of O2 and O interacting with variable Pd supercells. By analyzing hundreds of trajectories with appropriate initial sampling, the measured distance distribution of equilibrated atom pairs on Pd(111) is well reproduced. However, our results on Pd(100) suggest that the ballistic motion of hot atoms predicted previously is a rare event under ideal conditions, while initial molecular orientation and surface thermal fluctuation could significantly affect the overall postdissociation dynamics. On both surfaces, dissociated hyperthermal oxygen atoms primarily locate their nascent positions and experience a similar random walk motion nearby.
Collapse
Affiliation(s)
- Qidong Lin
- Key Laboratory of Precision and Intelligent Chemistry, 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, Anhui 230026, China
| | - Bin Jiang
- Key Laboratory of Precision and Intelligent Chemistry, 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, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, 230088, China
| |
Collapse
|
9
|
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
|
10
|
Foppa L, Rüther F, Geske M, Koch G, Girgsdies F, Kube P, Carey SJ, Hävecker M, Timpe O, Tarasov AV, Scheffler M, Rosowski F, Schlögl R, Trunschke A. Data-Centric Heterogeneous Catalysis: Identifying Rules and Materials Genes of Alkane Selective Oxidation. J Am Chem Soc 2023; 145:3427-3442. [PMID: 36745555 PMCID: PMC9936587 DOI: 10.1021/jacs.2c11117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Artificial intelligence (AI) can accelerate catalyst design by identifying key physicochemical descriptive parameters correlated with the underlying processes triggering, favoring, or hindering the performance. In analogy to genes in biology, these parameters might be called "materials genes" of heterogeneous catalysis. However, widely used AI methods require big data, and only the smallest part of the available data meets the quality requirement for data-efficient AI. Here, we use rigorous experimental procedures, designed to consistently take into account the kinetics of the catalyst active states formation, to measure 55 physicochemical parameters as well as the reactivity of 12 catalysts toward ethane, propane, and n-butane oxidation reactions. These materials are based on vanadium or manganese redox-active elements and present diverse phase compositions, crystallinities, and catalytic behaviors. By applying the sure-independence-screening-and-sparsifying-operator symbolic-regression approach to the consistent data set, we identify nonlinear property-function relationships depending on several key parameters and reflecting the intricate interplay of processes that govern the formation of olefins and oxygenates: local transport, site isolation, surface redox activity, adsorption, and the material dynamical restructuring under reaction conditions. These processes are captured by parameters derived from N2 adsorption, X-ray photoelectron spectroscopy (XPS), and near-ambient-pressure in situ XPS. The data-centric approach indicates the most relevant characterization techniques to be used for catalyst design and provides "rules" on how the catalyst properties may be tuned in order to achieve the desired performance.
Collapse
Affiliation(s)
- Lucas Foppa
- The
NOMAD Laboratory at the Fritz-Haber-Institut of the Max-Planck-Gesellschaft
and IRIS-Adlershof of the Humboldt-Universität zu Berlin, Faradayweg 4-6, D-14195 Berlin, Germany,
| | - Frederik Rüther
- BasCat
- UniCat BASF JointLab, Hardenbergstraße 36, D-10623 Berlin, Germany
| | - Michael Geske
- BasCat
- UniCat BASF JointLab, Hardenbergstraße 36, D-10623 Berlin, Germany
| | - Gregor Koch
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Frank Girgsdies
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Pierre Kube
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Spencer J. Carey
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Michael Hävecker
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany,Max
Planck Institute for Chemical Energy Conversion, 45470 Mülheim, Germany
| | - Olaf Timpe
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Andrey V. Tarasov
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Matthias Scheffler
- The
NOMAD Laboratory at the Fritz-Haber-Institut of the Max-Planck-Gesellschaft
and IRIS-Adlershof of the Humboldt-Universität zu Berlin, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Frank Rosowski
- BasCat
- UniCat BASF JointLab, Hardenbergstraße 36, D-10623 Berlin, Germany,BASF
SE, Catalysis Research, Carl-Bosch-Straße 38, D-67065 Ludwigshafen, Germany
| | - Robert Schlögl
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Annette Trunschke
- Department
of Inorganic Chemistry, Fritz-Haber-Institut
of the Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany,
| |
Collapse
|
11
|
Xu W, Yang B. Microkinetic modeling with machine learning predicted binding energies of reaction intermediates of ethanol steam reforming: The limitations. MOLECULAR CATALYSIS 2023. [DOI: 10.1016/j.mcat.2023.112940] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
|
12
|
Zha S, Sharapa DI, Liu S, Zhao ZJ, Studt F. Modeling CoCu Nanoparticles Using Neural Network-Accelerated Monte Carlo Simulations. J Phys Chem A 2022; 126:9440-9446. [PMID: 36512375 DOI: 10.1021/acs.jpca.2c07888] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The correct description of catalytic reactions happening on bimetallic particles is not feasible without proper accounting of the segregation process. In this study, we tried to shed light on the structure of large CoCu particles, for which quite controversial results were published before. However, density functional theory (DFT) is challenging to be directly used for the systematic study of nanometer-sized particles. Therefore, we constructed a neural network-based potential and further applied it to the Monte Carlo simulations for the description of the segregation phenomenon. The resulting approach shows high efficiency and can be used in systems with thousands of atoms. The accuracy and transferability of the model to other sizes and compositions make this methodology useful for solving segregation problems.
Collapse
Affiliation(s)
- Shenjun Zha
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz Platz 1, 76344Eggenstein-Leopoldshafen, Germany
| | - Dmitry I Sharapa
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz Platz 1, 76344Eggenstein-Leopoldshafen, Germany
| | - Sihang Liu
- Key Laboratory for Green Chemical Technology of Ministry of Education, School of Chemical Engineering and Technology, Collaborative Innovation Center of Chemical Science and Engineering, Tianjin University, Tianjin300072, China
| | - Zhi-Jian Zhao
- Key Laboratory for Green Chemical Technology of Ministry of Education, School of Chemical Engineering and Technology, Collaborative Innovation Center of Chemical Science and Engineering, Tianjin University, Tianjin300072, China
| | - Felix Studt
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz Platz 1, 76344Eggenstein-Leopoldshafen, Germany.,Institute for Chemical Technology and Polymer Chemistry, Karlsruhe Institute of Technology, Engesserstrasse 18, 76131Karlsruhe, Germany
| |
Collapse
|
13
|
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
|
14
|
Fuel Cell Reactors for the Clean Cogeneration of Electrical Energy and Value-Added Chemicals. ELECTROCHEM ENERGY R 2022. [DOI: 10.1007/s41918-022-00168-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
AbstractFuel cell reactors can be tailored to simultaneously cogenerate value-added chemicals and electrical energy while releasing negligible CO2 emissions or other pollution; moreover, some of these reactors can even “breathe in” poisonous gas as feedstock. Such clean cogeneration favorably offsets the fast depletion of fossil fuel resources and eases growing environmental concerns. These unique reactors inherit advantages from fuel cells: a high energy conversion efficiency and high selectivity. Compared with similar energy conversion devices with sandwich structures, fuel cell reactors have successfully “hit three birds with one stone” by generating power, producing chemicals, and maintaining eco-friendliness. In this review, we provide a systematic summary on the state of the art regarding fuel cell reactors and key components, as well as the typical cogeneration reactions accomplished in these reactors. Most strategies fall short in reaching a win–win situation that meets production demand while concurrently addressing environmental issues. The use of fuel cells (FCs) as reactors to simultaneously produce value-added chemicals and electrical power without environmental pollution has emerged as a promising direction. The FC reactor has been well recognized due to its “one stone hitting three birds” merit, namely, efficient chemical production, electrical power generation, and environmental friendliness. Fuel cell reactors for cogeneration provide multidisciplinary perspectives on clean chemical production, effective energy utilization, and even pollutant treatment, with far-reaching implications for the wider scientific community and society. The scope of this review focuses on unique reactors that can convert low-value reactants and/or industrial wastes to value-added chemicals while simultaneously cogenerating electrical power in an environmentally friendly manner.
Graphical Abstract
A schematic diagram for the concept of fuel cell reactors for cogeneration of electrical energy and value-added chemicals
Collapse
|
15
|
|
16
|
Pineda M, Stamatakis M. Kinetic Monte Carlo simulations for heterogeneous catalysis: Fundamentals, current status, and challenges. J Chem Phys 2022; 156:120902. [DOI: 10.1063/5.0083251] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Kinetic Monte Carlo (KMC) simulations in combination with first-principles (1p)-based calculations are rapidly becoming the gold-standard computational framework for bridging the gap between the wide range of length scales and time scales over which heterogeneous catalysis unfolds. 1p-KMC simulations provide accurate insights into reactions over surfaces, a vital step toward the rational design of novel catalysts. In this Perspective, we briefly outline basic principles, computational challenges, successful applications, as well as future directions and opportunities of this promising and ever more popular kinetic modeling approach.
Collapse
Affiliation(s)
- M. Pineda
- Thomas Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London WC1E 7JE, United Kingdom
| | - M. Stamatakis
- Thomas Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London WC1E 7JE, United Kingdom
| |
Collapse
|
17
|
Steiner M, Reiher M. Autonomous Reaction Network Exploration in Homogeneous and Heterogeneous Catalysis. Top Catal 2022; 65:6-39. [PMID: 35185305 PMCID: PMC8816766 DOI: 10.1007/s11244-021-01543-9] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/17/2021] [Indexed: 12/11/2022]
Abstract
Autonomous computations that rely on automated reaction network elucidation algorithms may pave the way to make computational catalysis on a par with experimental research in the field. Several advantages of this approach are key to catalysis: (i) automation allows one to consider orders of magnitude more structures in a systematic and open-ended fashion than what would be accessible by manual inspection. Eventually, full resolution in terms of structural varieties and conformations as well as with respect to the type and number of potentially important elementary reaction steps (including decomposition reactions that determine turnover numbers) may be achieved. (ii) Fast electronic structure methods with uncertainty quantification warrant high efficiency and reliability in order to not only deliver results quickly, but also to allow for predictive work. (iii) A high degree of autonomy reduces the amount of manual human work, processing errors, and human bias. Although being inherently unbiased, it is still steerable with respect to specific regions of an emerging network and with respect to the addition of new reactant species. This allows for a high fidelity of the formalization of some catalytic process and for surprising in silico discoveries. In this work, we first review the state of the art in computational catalysis to embed autonomous explorations into the general field from which it draws its ingredients. We then elaborate on the specific conceptual issues that arise in the context of autonomous computational procedures, some of which we discuss at an example catalytic system.
Collapse
Affiliation(s)
- Miguel Steiner
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
18
|
Li Z. First-principles-based microkinetic rate equation theory for oxygen carrier reduction in chemical looping. Chem Eng Sci 2022. [DOI: 10.1016/j.ces.2021.117042] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
19
|
Gerrits N. Accurate Simulations of the Reaction of H 2 on a Curved Pt Crystal through Machine Learning. J Phys Chem Lett 2021; 12:12157-12164. [PMID: 34918518 PMCID: PMC8724818 DOI: 10.1021/acs.jpclett.1c03395] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Accepted: 12/14/2021] [Indexed: 06/14/2023]
Abstract
Theoretical studies on molecule-metal surface reactions have so far been limited to small surface unit cells due to computational costs. Here, for the first time molecular dynamics simulations on very large surface unit cells at the level of density functional theory are performed, allowing a direct comparison to experiments performed on a curved crystal. Specifically, the reaction of D2 on a curved Pt crystal is investigated with a neural network potential (NNP). The developed NNP is also accurate for surface unit cells considerably larger than those that have been included in the training data, allowing dynamical simulations on very large surface unit cells that otherwise would have been intractable. Important and complex aspects of the reaction mechanism are discovered such as diffusion and a shadow effect of the step. Furthermore, conclusions from simulations on smaller surface unit cells cannot always be transfered to larger surface unit cells, limiting the applicability of theoretical studies of smaller surface unit cells to heterogeneous catalysts with small defect densities.
Collapse
Affiliation(s)
- Nick Gerrits
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, BE-2610 Wilrijk, Antwerp, Belgium
| |
Collapse
|
20
|
Zhang L, Cai J, Chen Y, Huang J. Modelling electrocatalytic reactions with a concerted treatment of multistep electron transfer kinetics and local reaction conditions. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33. [PMID: 34525456 DOI: 10.1088/1361-648x/ac26fb] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Accepted: 09/15/2021] [Indexed: 05/17/2023]
Abstract
A physicochemical model is developed for electrocatalytic reactions involving multiple electron transfer steps occurring in the electric double layer (EDL). The local reaction conditions are calculated using a mean-field EDL model, which is derived from a comprehensive grand potential that considers the steric effects, solvent polarization, and chemisorption-induced surface dipoles. Macroscopic mass transport in the so-called diffusion layer is controlled by the same set of controlling equations of the EDL model, without imposing the electroneutrality assumption as usual. The Gerischer's formulation of electron transfer theory, corrected with local reaction conditions, is used to describe the kinetics of elementary steps. Multistep kinetics of the electrocatalytic reaction is treated using microkinetics modelling, without resorting to the usual rate-determining step approximation. In formal analysis of the model, we retrieve canonical models with additional assumptions. Self-consistent numerical implementation of the model is demonstrated for oxygen reduction reaction (ORR) at Pt(111) in acidic solution, and the aptness of the model is verified by comparison with experimental data. A comparative study of the full model and its simplified versions allows us to examine how the ORR is influenced by asymmetric steric effects, finite concentration of ions, solvent polarization, surface charge effects, and metal electronic structure effects. We find that the difference in terms of the overpotential between the full model and the simplest model can be up to ∼0.1 V at a current density of -6 mAcm-2.
Collapse
Affiliation(s)
- Lulu Zhang
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei 230026, People's Republic of China
| | - Jun Cai
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei 230026, People's Republic of China
| | - Yanxia Chen
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei 230026, People's Republic of China
| | - Jun Huang
- College of Chemistry and Chemical Engineering, Central South University, Changsha 410083, People's Republic of China
| |
Collapse
|
21
|
Huang J, Zhu X, Eikerling M. The rate-determining term of electrocatalytic reactions with first-order kinetics. Electrochim Acta 2021. [DOI: 10.1016/j.electacta.2021.139019] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
22
|
Achievements and Expectations in the Field of Computational Heterogeneous Catalysis in an Innovation Context. Top Catal 2021. [DOI: 10.1007/s11244-021-01489-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
23
|
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: 5] [Impact Index Per Article: 1.7] [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
|
24
|
Bhattacharjee H, Anesiadis N, Vlachos DG. Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies. Sci Rep 2021; 11:14372. [PMID: 34257362 PMCID: PMC8277863 DOI: 10.1038/s41598-021-93854-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2021] [Accepted: 06/22/2021] [Indexed: 11/25/2022] Open
Abstract
A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.
Collapse
Affiliation(s)
- Himaghna Bhattacharjee
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, DE, 19716, USA.,RAPID Manufacturing Institute and Delaware Energy Institute (DEI), 221 Academy Street, Newark, DE, 19716, USA
| | - Nikolaos Anesiadis
- Department of Chemical Engineering and Applied Chemistry, University of Toronto, 200 College St., Toronto, ON, M5S 3E5, Canada
| | - Dionisios G Vlachos
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, DE, 19716, USA. .,RAPID Manufacturing Institute and Delaware Energy Institute (DEI), 221 Academy Street, Newark, DE, 19716, USA.
| |
Collapse
|
25
|
Fajín JLC, Moura AS, Cordeiro MNDS. First-principles-based kinetic Monte Carlo simulations of CO oxidation on catalytic Au(110) and Ag(110) surfaces. Phys Chem Chem Phys 2021; 23:14037-14050. [PMID: 34151916 DOI: 10.1039/d1cp00729g] [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/21/2022]
Abstract
At the core of the development of more efficient and reliable fuel cells (FCs), there are several essential chemical reactions, namely carbon monoxide (CO) oxidation. This reaction is a keystone in the cleaning of hydrogen fuel used in fuel cells due to strong poisoning by this species of the platinum catalyst used in these devices. The present work aims to provide insight regarding the activation of CO oxidation by gold or silver microfacets possessing low coordinated atoms. To achieve this, density functional theory (DFT) quantum calculations, which determined two competing reaction pathways for CO oxidation, i.e., by molecularly adsorbed oxygen, and by dissociated oxygen, are combined with first-principles kinetic Monte Carlo (1p-kMC) simulations, which employed the resulting DFT parameters in order to address the effect of temperature and partial pressures and the interplay of the elementary reaction events. The use of 1p-kMC is a step further from available works regarding the CO oxidation on gold- and silver-based catalysts for cleansing of hydrogen that is used as a fuel in FCs. Indeed, this research contributes to the conclusion that CO oxidation should preferentially occur on silver microfacets, while the obtained turnover frequencies (TOFs) reinforced such a conclusion.
Collapse
Affiliation(s)
- Jose L C Fajín
- LAQV@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, 4169-007 Porto, Portugal.
| | | | | |
Collapse
|
26
|
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
|
27
|
Réocreux R, Fampiou I, Stamatakis M. The role of oxygenated species in the catalytic self-coupling of MeOH on O pre-covered Au(111). Faraday Discuss 2021; 229:251-266. [PMID: 33646205 DOI: 10.1039/c9fd00134d] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The oxidation of alcohols plays a central role in the valorisation of biomass, in particular when performed with a non-toxic oxidant such as O2. Aerobic oxidation of methanol on gold has attracted attention lately and the main steps of its mechanism have been described experimentally. However, the exact role of O and OH on each elementary step and the effect of the interactions between adsorbates are still not completely understood. Here we investigate the mechanism of methanol oxidation to HCOOCH3 and CO2. We use Density Functional Theory (DFT) to assess the energetics of the underlying pathways, and subsequently build lattice kinetic Monte Carlo (kMC) models of increasing complexity, to elucidate the role of different oxygenates. Detailed comparisons of our simulation results with experimental temperature programmed desorption (TPD) spectra enable us to validate the mechanism and identify rate determining steps. Crucially, taking into account dispersion (van der Waals forces) and adsorbate-adsorbate lateral interactions are both important for reproducing the experimental data.
Collapse
Affiliation(s)
- R Réocreux
- Thomas Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London, WC1E 7JE, UK.
| | - I Fampiou
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - M Stamatakis
- Thomas Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London, WC1E 7JE, UK.
| |
Collapse
|
28
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (Eb, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate Eb can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for Eb for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH4. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of Eb for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate Eb for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O2, H2O, NH3, CO2, and CH3OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
29
|
Nezam I, Zhou W, Gusmão GS, Realff MJ, Wang Y, Medford AJ, Jones CW. Direct aromatization of CO2 via combined CO2 hydrogenation and zeolite-based acid catalysis. J CO2 UTIL 2021. [DOI: 10.1016/j.jcou.2020.101405] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
30
|
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
|
31
|
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
|
32
|
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
|
33
|
Auras SV, van Lent R, Bashlakov D, Piñeiros Bastidas JM, Roorda T, Spierenburg R, Juurlink LBF. Scaling Platinum-Catalyzed Hydrogen Dissociation on Corrugated Surfaces. Angew Chem Int Ed Engl 2020; 59:20973-20979. [PMID: 32749736 PMCID: PMC7692953 DOI: 10.1002/anie.202005616] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 07/17/2020] [Indexed: 11/07/2022]
Abstract
We determine absolute reactivities for dissociation at low coordinated Pt sites. Two curved Pt(111) single-crystal surfaces allow us to probe either straight or highly kinked step edges with molecules impinging at a low impact energy. A model extracts the average reactivity of inner and outer kink atoms, which is compared to the reactivity of straight A- and B-type steps. Local surface coordination numbers do not adequately capture reactivity trends for H2 dissociation. We utilize the increase of reactivity with step density to determine the area over which a step causes increased dissociation. This step-type specific reactive area extends beyond the step edge onto the (111) terrace. It defines the reaction cross-section for H2 dissociation at the step, bypassing assumptions about contributions of individual types of surface atoms. Our results stress the non-local nature of H2 interaction with a surface and provide insight into reactivity differences for nearly identical step sites.
Collapse
Affiliation(s)
- Sabine V Auras
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Richard van Lent
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Dima Bashlakov
- ILTPE, National Academy of Sciences of Ukraine, 47 Nauky Ave., Kharkiv, 61103, Ukraine
| | | | - Tycho Roorda
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Rick Spierenburg
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Ludo B F Juurlink
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| |
Collapse
|
34
|
Auras SV, Lent R, Bashlakov D, Piñeiros Bastidas JM, Roorda T, Spierenburg R, Juurlink LBF. Scaling Platinum‐Catalyzed Hydrogen Dissociation on Corrugated Surfaces. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202005616] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Sabine V. Auras
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Richard Lent
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Dima Bashlakov
- ILTPE, National Academy of Sciences of Ukraine 47 Nauky Ave. Kharkiv 61103 Ukraine
| | | | - Tycho Roorda
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Rick Spierenburg
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Ludo B. F. Juurlink
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| |
Collapse
|
35
|
Ravipati S, d'Avezac M, Nielsen J, Hetherington J, Stamatakis M. A Caching Scheme To Accelerate Kinetic Monte Carlo Simulations of Catalytic Reactions. J Phys Chem A 2020; 124:7140-7154. [PMID: 32786994 DOI: 10.1021/acs.jpca.0c03571] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Kinetic Monte Carlo (KMC) simulations have been instrumental in advancing our fundamental understanding of heterogeneously catalyzed reactions, with particular emphasis on structure sensitivity, ensemble effects, and the interplay between adlayer structure and adsorbate-adsorbate lateral interactions in shaping the observed kinetics. Yet, the computational cost of KMC remains high, thereby motivating the development of acceleration schemes that would improve the simulation efficiency. We present an exact such scheme, which implements a caching algorithm along with shared-memory parallelization to improve the computational performance of simulations incorporating long-range adsorbate-adsorbate lateral interactions. This scheme is based on caching information about the energetic interaction patterns associated with the products of each possible lattice process (adsorption, desorption, reaction etc.). Thus, every time a reaction occurs ("ongoing reaction"), it enables fast updates of the rate constants of "affected reactions", i.e., possible reactions in the region of influence of the "ongoing reaction". Benchmarks on KMC simulations of NOx oxidation/reduction, yielded acceleration factors of up to 20, when comparing single-thread runs without caching to runs on 16 threads with caching, for simulations with a cluster expansion Hamiltonian that incorporates up to 8th-nearest-neighbor interactions.
Collapse
Affiliation(s)
- Srikanth Ravipati
- Thomas Young Centre and Department of Chemical Engineering, University College London, Torrington Place, London WC1E 7JE, United Kingdom
| | - Mayeul d'Avezac
- Research Software Development Group, Research IT Services, University College London, Torrington Place, London WC1E 6BT, United Kingdom
| | - Jens Nielsen
- Research Software Development Group, Research IT Services, University College London, Torrington Place, London WC1E 6BT, United Kingdom
| | - James Hetherington
- Research Software Development Group, Research IT Services, University College London, Torrington Place, London WC1E 6BT, United Kingdom
| | - Michail Stamatakis
- Thomas Young Centre and Department of Chemical Engineering, University College London, Torrington Place, London WC1E 7JE, United Kingdom
| |
Collapse
|
36
|
Trends and Outlook of Computational Chemistry and Microkinetic Modeling for Catalytic Synthesis of Methanol and DME. Catalysts 2020. [DOI: 10.3390/catal10060655] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
The first-principle modeling of heterogeneous catalysts is a revolutionarily approach, as the electronic structure of a catalyst is closely related to its reactivity on the surface with reactant molecules. In the past, detailed reaction mechanisms could not be understood, however, computational chemistry has made it possible to analyze a specific elementary reaction of a reaction system. Microkinetic modeling is a powerful tool for investigating elementary reactions and reaction mechanisms for kinetics. Using a microkinetic model, the dominant pathways and rate-determining steps can be elucidated among the competitive reactions, and the effects of operating conditions on the reaction mechanisms can be determined. Therefore, the combination of computational chemistry and microkinetic modeling can significantly improve computational catalysis research. In this study, we reviewed the trends and outlook of this combination technique as applied to the catalytic synthesis of methanol (MeOH) and dimethyl ether (DME), whose detailed mechanisms are still controversial. Although the scope is limited to the catalytic synthesis of limited species, this study is expected to provide a foundation for future works in the field of catalysis research based on computational catalysis.
Collapse
|
37
|
Sengar A, van Santen RA, Kuipers JA. Deactivation Kinetics of the Catalytic Alkylation Reaction. ACS Catal 2020. [DOI: 10.1021/acscatal.0c00932] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Aditya Sengar
- Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
| | - Rutger A. van Santen
- Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
- Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
| | - Johannes A.M. Kuipers
- Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
| |
Collapse
|
38
|
Cheula R, Maestri M, Mpourmpakis G. Modeling Morphology and Catalytic Activity of Nanoparticle Ensembles Under Reaction Conditions. ACS Catal 2020. [DOI: 10.1021/acscatal.0c01005] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Raffaele Cheula
- Laboratory of Catalysis and Catalytic Processes, Dipartimento di Energia, Politecnico di Milano, Via La Masa, 34, 20156 Milano, Italy
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, 4200 Fifth Avenue, Pittsburgh, Pennsylvania 15260, United States
| | - Matteo Maestri
- Laboratory of Catalysis and Catalytic Processes, Dipartimento di Energia, Politecnico di Milano, Via La Masa, 34, 20156 Milano, Italy
| | - Giannis Mpourmpakis
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, 4200 Fifth Avenue, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
39
|
Abstract
Modern computational chemistry has reached a stage at which massive exploration into chemical reaction space with unprecedented resolution with respect to the number of potentially relevant molecular structures has become possible. Various algorithmic advances have shown that such structural screenings must and can be automated and routinely carried out. This will replace the standard approach of manually studying a selected and restricted number of molecular structures for a chemical mechanism. The complexity of the task has led to many different approaches. However, all of them address the same general target, namely to produce a complete atomistic picture of the kinetics of a chemical process. It is the purpose of this overview to categorize the problems that should be targeted and to identify the principal components and challenges of automated exploration machines so that the various existing approaches and future developments can be compared based on well-defined conceptual principles.
Collapse
Affiliation(s)
- Jan P. Unsleber
- Laboratory for Physical Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory for Physical Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| |
Collapse
|
40
|
Chizallet C. Toward the Atomic Scale Simulation of Intricate Acidic Aluminosilicate Catalysts. ACS Catal 2020. [DOI: 10.1021/acscatal.0c01136] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Céline Chizallet
- IFP Energies nouvelles Solaize, Rond-Point de l’Echangeur de Solaize, BP 3, 69360 Solaize, France
| |
Collapse
|
41
|
Organic Carbonate Production Utilizing Crude Glycerol Derived as By-Product of Biodiesel Production: A Review. ENERGIES 2020. [DOI: 10.3390/en13061483] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
As a promising alternative renewable liquid fuel, biodiesel production has increased and eventually led to an increase in the production of its by-product, crude glycerol. The vast generation of glycerol has surpassed the market demand. Hence, the crude glycerol produced should be utilized effectively to increase the viability of biodiesel production. One of them is through crude glycerol upgrading, which is not economical. A good deal of attention has been dedicated to research for alternative material and chemicals derived from sustainable biomass resources. It will be more valuable if the crude glycerol is converted into glycerol derivatives, and so, increase the economic possibility of the biodiesel production. Studies showed that glycerol carbonate plays an important role, as a building block, in synthesizing the glycerol oligomers at milder conditions under microwave irradiation. This review presents a brief outline of the physio-chemical, thermodynamic, toxicological, production methods, reactivity, and application of organic carbonates derived from glycerol with a major focus on glycerol carbonate and dimethyl carbonate (DMC), as a green chemical, for application in the chemical and biotechnical field. Research gaps and further improvements have also been discussed.
Collapse
|
42
|
Andersson S, Radl S, Svenum IH, Shevlin SA, Guo ZX, Amini S. Towards rigorous multiscale flow models of nanoparticle reactivity in chemical looping applications. Catal Today 2019. [DOI: 10.1016/j.cattod.2019.06.024] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
43
|
Smeets EF, Füchsel G, Kroes GJ. Quantum Dynamics of Dissociative Chemisorption of H 2 on the Stepped Cu(211) Surface. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:23049-23063. [PMID: 31565113 PMCID: PMC6757508 DOI: 10.1021/acs.jpcc.9b06539] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/10/2019] [Revised: 08/22/2019] [Indexed: 06/10/2023]
Abstract
Reactions on stepped surfaces are relevant to heterogeneous catalysis, in which a reaction often takes place at the edges of nanoparticles where the edges resemble steps on single-crystal stepped surfaces. Previous results on H2 + Cu(211) showed that, in this system, steps do not enhance the reactivity and raised the question of whether this effect could be, in any way, related to the neglect of quantum dynamical effects in the theory. To investigate this, we present full quantum dynamical molecular beam simulations of sticking of H2 on Cu(211), in which all important rovibrational states populated in a molecular beam experiment are taken into account. We find that the reaction of H2 with Cu(211) is very well described with quasi-classical dynamics when simulating molecular beam sticking experiments, in which averaging takes place over a large number of rovibrational states and over translational energy distributions. Our results show that the stepped Cu(211) surface is distinct from its component Cu(111) terraces and Cu(100) steps and cannot be described as a combination of its component parts with respect to the reaction dynamics when considering the orientational dependence. Specifically, we present evidence that, at translational energies close to the reaction threshold, vibrationally excited molecules show a negative rotational quadrupole alignment parameter on Cu(211), which is not found on Cu(111) and Cu(100). The effect arises because these molecules react with a site-specific reaction mechanism at the step, that is, inelastic rotational enhancement, which is only effective for molecules with a small absolute value of the magnetic rotation quantum number. From a comparison to recent associative desorption experiments as well as Born-Oppenheimer molecular dynamics calculations, it follows that the effects of surface atom motion and electron-hole pair excitation on the reactivity fall within chemical accuracy, that is, modeling these effect shifts extracted reaction probability curves by less than 1 kcal/mol translational energy. We found no evidence in our fully state-resolved calculations for the "slow" reaction channel that was recently reported for associative desorption of H2 from Cu(111) and Cu(211), but our results for the fast channel are in good agreement with the experiments on H2 + Cu(211).
Collapse
Affiliation(s)
- Egidius
W. F. Smeets
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Gernot Füchsel
- Institut
für Chemie und Biochemie - Physikalische und Theoretische Chemie, Freie Universität Berlin, Takustraße 3, 14195 Berlin, Germany
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
44
|
Calderón LA, Montoya A, Soon A, Stampfl C. Non-dissociative adsorption of glycerol on the (111) surface of Ni and Pt-based metallic systems: Hints on reforming activity from d-band center. MOLECULAR CATALYSIS 2019. [DOI: 10.1016/j.mcat.2019.110412] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
45
|
Tian H, Rangarajan S. Predicting Adsorption Energies Using Multifidelity Data. J Chem Theory Comput 2019; 15:5588-5600. [DOI: 10.1021/acs.jctc.9b00336] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Huijie Tian
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem 18015, United States
| | - Srinivas Rangarajan
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem 18015, United States
| |
Collapse
|
46
|
Bruix A, Margraf JT, Andersen M, Reuter K. First-principles-based multiscale modelling of heterogeneous catalysis. Nat Catal 2019. [DOI: 10.1038/s41929-019-0298-3] [Citation(s) in RCA: 131] [Impact Index Per Article: 26.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
47
|
Smeets EWF, 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] [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
|
48
|
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
|
49
|
Andersen M, Panosetti C, Reuter K. A Practical Guide to Surface Kinetic Monte Carlo Simulations. Front Chem 2019; 7:202. [PMID: 31024891 PMCID: PMC6465329 DOI: 10.3389/fchem.2019.00202] [Citation(s) in RCA: 83] [Impact Index Per Article: 16.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2018] [Accepted: 03/15/2019] [Indexed: 11/26/2022] Open
Abstract
This review article is intended as a practical guide for newcomers to the field of kinetic Monte Carlo (KMC) simulations, and specifically to lattice KMC simulations as prevalently used for surface and interface applications. We will provide worked out examples using the kmos code, where we highlight the central approximations made in implementing a KMC model as well as possible pitfalls. This includes the mapping of the problem onto a lattice and the derivation of rate constant expressions for various elementary processes. Example KMC models will be presented within the application areas surface diffusion, crystal growth and heterogeneous catalysis, covering both transient and steady-state kinetics as well as the preparation of various initial states of the system. We highlight the sensitivity of KMC models to the elementary processes included, as well as to possible errors in the rate constants. For catalysis models in particular, a recurrent challenge is the occurrence of processes at very different timescales, e.g., fast diffusion processes and slow chemical reactions. We demonstrate how to overcome this timescale disparity problem using recently developed acceleration algorithms. Finally, we will discuss how to account for lateral interactions between the species adsorbed to the lattice, which can play an important role in all application areas covered here.
Collapse
Affiliation(s)
- Mie Andersen
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Garching, Germany
| | | | | |
Collapse
|
50
|
Konwar M, Chetia M, Sarma D. A Low-Cost, Well-Designed Catalytic System Derived from Household Waste "Egg Shell": Applications in Organic Transformations. Top Curr Chem (Cham) 2019; 377:6. [PMID: 30675643 DOI: 10.1007/s41061-018-0230-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2018] [Accepted: 12/29/2018] [Indexed: 01/22/2023]
Abstract
A waste feedstock-derived economical basic alternative catalyst is described in this review. Eggshell is one of the household wastes created in tons of weight daily. Therefore, in order to reduce the environmental pollution-related problems, its use in heterogeneous catalysis can be attributed as a great contribution for the chemical and material science society to carry out several known reactions and for the much-needed energy alternative biodiesel production as low-cost catalytic system. Keeping green chemistry in mind, industrial use of these catalysts may also reduce the use of other traditionally used high-cost chemical catalytic systems.
Collapse
Affiliation(s)
- Manashjyoti Konwar
- Department of Chemistry, Dibrugarh University, Dibrugarh, Assam, 786004, India
| | - Mitali Chetia
- Department of Chemistry, Dibrugarh University, Dibrugarh, Assam, 786004, India
| | - Diganta Sarma
- Department of Chemistry, Dibrugarh University, Dibrugarh, Assam, 786004, India.
| |
Collapse
|