1
|
Korobov A. A possibility to infer frustrations of supported catalytic clusters from macro-scale observations. Sci Rep 2024; 14:3801. [PMID: 38361133 PMCID: PMC10869823 DOI: 10.1038/s41598-024-54485-z] [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: 11/27/2023] [Accepted: 02/12/2024] [Indexed: 02/17/2024] Open
Abstract
Recent experimental and theoretical studies suggest that dynamic active centres of supported heterogeneous catalysts may, under certain conditions, be frustrated. Such out-of-equilibrium materials are expected to possess unique catalytic properties and also higher level of functionality. The latter is associated with the navigation through the free energy landscapes with energetically close local minima. The lack of common approaches to the study of out-of-equilibrium materials motivates the search for specific ones. This paper suggests a way to infer some valuable information from the interplay between the intensity of reagent supply and regularities of product formation.
Collapse
Affiliation(s)
- Alexander Korobov
- Materials Chemistry Department, V. N. Karazin Kharkiv National University, Kharkiv, 61022, Ukraine.
| |
Collapse
|
2
|
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
|
3
|
Van Speybroeck V, Bocus M, Cnudde P, Vanduyfhuys L. Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations. ACS Catal 2023; 13:11455-11493. [PMID: 37671178 PMCID: PMC10476167 DOI: 10.1021/acscatal.3c01945] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2023] [Revised: 07/27/2023] [Indexed: 09/07/2023]
Abstract
Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.
Collapse
Affiliation(s)
| | - Massimo Bocus
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Pieter Cnudde
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Louis Vanduyfhuys
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| |
Collapse
|
4
|
Rajan A, Pushkar AP, Dharmalingam BC, Varghese JJ. Iterative multiscale and multi-physics computations for operando catalyst nanostructure elucidation and kinetic modeling. iScience 2023; 26:107029. [PMID: 37360694 PMCID: PMC10285649 DOI: 10.1016/j.isci.2023.107029] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023] Open
Abstract
Modern heterogeneous catalysis has benefitted immensely from computational predictions of catalyst structure and its evolution under reaction conditions, first-principles mechanistic investigations, and detailed kinetic modeling, which are rungs on a multiscale workflow. Establishing connections across these rungs and integration with experiments have been challenging. Here, operando catalyst structure prediction techniques using density functional theory simulations and ab initio thermodynamics calculations, molecular dynamics, and machine learning techniques are presented. Surface structure characterization by computational spectroscopic and machine learning techniques is then discussed. Hierarchical approaches in kinetic parameter estimation involving semi-empirical, data-driven, and first-principles calculations and detailed kinetic modeling via mean-field microkinetic modeling and kinetic Monte Carlo simulations are discussed along with methods and the need for uncertainty quantification. With these as the background, this article proposes a bottom-up hierarchical and closed loop modeling framework incorporating consistency checks and iterative refinements at each level and across levels.
Collapse
Affiliation(s)
- Ajin Rajan
- Department of Chemical Engineering, Indian Institute of Technology Madras, Chennai, Tamil Nadu 600036, India
| | - Anoop P. Pushkar
- Department of Chemical Engineering, Indian Institute of Technology Madras, Chennai, Tamil Nadu 600036, India
| | - Balaji C. Dharmalingam
- Department of Chemical Engineering, Indian Institute of Technology Madras, Chennai, Tamil Nadu 600036, India
| | - Jithin John Varghese
- Department of Chemical Engineering, Indian Institute of Technology Madras, Chennai, Tamil Nadu 600036, India
| |
Collapse
|
5
|
Zhang Y, Ren D. Mechanisms for Catalytic CO Oxidation on SiAu n ( n = 1-5) Cluster. Molecules 2023; 28:1917. [PMID: 36838905 PMCID: PMC9962203 DOI: 10.3390/molecules28041917] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Revised: 02/14/2023] [Accepted: 02/14/2023] [Indexed: 02/19/2023] Open
Abstract
Significant progress has been made in understanding the reactivity and catalytic activity of gas-phase and loaded gold clusters for CO oxidation. However, little research has focused on mixed silicon/gold clusters (SiAun) for CO oxidation. In the present work, we performed density function theory (DFT) calculations for a SiAun (n = 1-5) cluster at the CAM-B3LYP/aug-cc-pVDZ-PP level and investigated the effects on the reactivity and catalytic activity of the SiAun cluster for CO oxidation. The calculated results show that the effect is very low for the activation barriers for the formation of OOCO intermediates on SiAu clusters, SiAu3 clusters, and SiAu5 clusters in the catalytic oxidation of CO and the activation energy barriers for the formation of OCO intermediates on OSiAu3, OSiAu4, and OSiAu5. Our calculations show that, compared with the conventional small Au cluster, the incorporation of Si enhances the catalytic performance towards CO oxidation.
Collapse
Affiliation(s)
| | - Dasen Ren
- College of Chemical Engineering, Guizhou Minzu University, Guiyang 550025, China
| |
Collapse
|
6
|
Samadi P, Binczarski MJ, Maniukiewicz W, Pawlaczyk A, Rogowski J, Szubiakiewicz E, Szynkowska-Jozwik MI, Witonska IA. Zn Modification of Pd/TiO 2/Ti Catalyst for CO Oxidation. MATERIALS (BASEL, SWITZERLAND) 2023; 16:1216. [PMID: 36770219 PMCID: PMC9921276 DOI: 10.3390/ma16031216] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Revised: 01/22/2023] [Accepted: 01/28/2023] [Indexed: 06/18/2023]
Abstract
The main goal of this study was to modify the activity of Pd/TiO2/Ti catalyst in the reaction of CO oxidation by the addition of Zn. Plasma electrolytic oxidation (PEO) of Ti wire was conducted to produce a uniform porous layer of TiO2. A mixture of Pd and Zn was then introduced by means of adsorption. After reduction treatment, the activity of the samples was examined by oxidation of 5% CO in a temperature range from 80-350 °C. Model catalysts with sufficient amounts of the metals for physico-chemical investigation were prepared to further investigate the reaction between Pd and Zn during CO oxidation. The structures and compositions of the samples were investigated using scanning electron microscopy with energy dispersive spectroscopy (SEM-EDS), time-of-flight secondary ion mass spectrometry (TOF-SIMS), inductively coupled plasma mass spectrometry (ICP-MS), Fourier transform infrared (FTIR), and X-ray diffraction (XRD). Modification of Pd/TiO2/Ti catalyst by Zn with a Pd:Zn atomic ratio of 2:1 decreased the temperature of complete CO oxidation from 220 °C for Pd/TiO2/Ti to 180 °C for Pd-Zn/TiO2/Ti. The temperature of 50% CO conversion on Pd-Zn(2:1)/TiO2/Ti was around 55 °C lower than in the reaction on monometallic Pd catalyst. The addition of Zn to the Pd catalyst lowered the binding energy of CO on the surface and improved the dissociative adsorption of oxygen, facilitating the oxidation of CO. FTIR showed that the bridging form of adsorbed CO is preferred on bimetallic systems. Analysis of the surface compositions of the samples (SEM-EDS, TOF-SIMS) showed higher amounts of oxygen on the bimetallic systems.
Collapse
|
7
|
|
8
|
Kennes K, Kubarev A, Demaret C, Treps L, Delpoux O, Rivallan M, Guillon E, Méthivier A, de Bruin T, Gomez A, Harbuzaru B, Roeffaers MB, Chizallet C. Multiscale Visualization and Quantification of the Effect of Binders on the Acidity of Shaped Zeolites. ACS Catal 2022. [DOI: 10.1021/acscatal.2c02152] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Koen Kennes
- Centre for Membrane Separations, Adsorption, Catalysis and Spectroscopy for Sustainable Solutions, Faculty of Bioscience Engineering, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Alexey Kubarev
- Centre for Membrane Separations, Adsorption, Catalysis and Spectroscopy for Sustainable Solutions, Faculty of Bioscience Engineering, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Coralie Demaret
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Laureline Treps
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Olivier Delpoux
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Mickael Rivallan
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Emmanuelle Guillon
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Alain Méthivier
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Theodorus de Bruin
- IFP Energies nouvelles, 1 et 4 Avenue de Bois-Préau, BP3, 92852 Rueil-Malmaison, France
| | - Axel Gomez
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
- Département de Chimie, École Normale Supérieure, PSL University, 75005 Paris, France
| | - Bogdan Harbuzaru
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| | - Maarten B.J. Roeffaers
- Centre for Membrane Separations, Adsorption, Catalysis and Spectroscopy for Sustainable Solutions, Faculty of Bioscience Engineering, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Céline Chizallet
- IFP Energies nouvelles, Rond-Point de L’Echangeur de Solaize, BP 3 69360 Solaize, France
| |
Collapse
|
9
|
Kellouai W, Judeinstein P, Plazanet M, Baudoin S, Drobek M, Julbe A, Coasne B. Gas Adsorption in Zeolite and Thin Zeolite Layers: Molecular Simulation, Experiment, and Adsorption Potential Theory. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:5428-5438. [PMID: 35486814 DOI: 10.1021/acs.langmuir.1c03420] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Molecular simulations and experiments are used to investigate methane adsorption in bulk and thin layers of MFI zeolite (silicalite-1). After comparing the theoretical adsorption data obtained using Grand Canonical Monte Carlo simulations for bulk MFI at various temperatures against experiments, zeolite layers with different crystalline orientations and levels of surface flexibility are considered. The data obtained for such prototypical systems allow us to rationalize both the qualitative and quantitative impact of external surface in nanoporous solids. In particular, due to strong confinement in zeolite pores, methane is found to adsorb at low pressures in the core of the zeolite while external surface adsorption occurs at pressures where the internal porosity of zeolite is saturated. Using Polanyi's adsorption potential theory, which is derived here from Hill's general scheme for adsorption, we provide a simple thermodynamic formalism to predict consistently adsorption both in the internal porosity and at the external surface of nanoporous solids. While this seminal theory has been already applied for gases in nanoporous solids, its extension to describe both surface and volume adsorption is important to provide a general rational framework for fluid adsorption in finely divided materials. We also discuss the applicability of this formalism for gas adsorption data under supercritical conditions.
Collapse
Affiliation(s)
- Wanda Kellouai
- Université Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| | | | - Marie Plazanet
- Université Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| | - Simon Baudoin
- Institut Laue Langevin, 71 Av. des Martyrs, 38000 Grenoble, France
| | - Martin Drobek
- Institut Européen des Membranes, CNRS, ENSCM, Univ Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Anne Julbe
- Institut Européen des Membranes, CNRS, ENSCM, Univ Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Benoit Coasne
- Université Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| |
Collapse
|
10
|
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
|
11
|
Trunschke A. Prospects and challenges for autonomous catalyst discovery viewed from an experimental perspective. Catal Sci Technol 2022. [DOI: 10.1039/d2cy00275b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Autonomous catalysis research requires elaborate integration of operando experiments into automated workflows. Suitable experimental data for analysis by artificial intelligence can be measured more readily according to standard operating procedures.
Collapse
Affiliation(s)
- Annette Trunschke
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Department of Inorganic Chemistry, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
12
|
Pigeon T, Chizallet C, Raybaud P. Revisiting γ-alumina surface models through the topotactic transformation of boehmite surfaces. J Catal 2022. [DOI: 10.1016/j.jcat.2021.11.011] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|