1
|
Blöndal K, Badger K, Sargsyan K, Bross DH, Ruscic B, Goldsmith CF. Importance sampling within configuration space integration for adsorbate thermophysical properties: a case study for CH 3/Ni(111). Phys Chem Chem Phys 2024; 26:17265-17273. [PMID: 38856369 DOI: 10.1039/d4cp01197j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2024]
Abstract
A new strategy is presented for computing anharmonic partition functions for the motion of adsorbates relative to a catalytic surface. Importance sampling is compared with conventional Monte Carlo. The importance sampling is significantly more efficient. This new approach is applied to CH3* on Ni(111) as a test case. The motion of methyl relative to the nickel surface is found to be anharmonic, with significantly higher entropy compared to the standard harmonic oscillator model. The new method is freely available as part of the Minima-Preserving Neural Network within the ADTHERM package.
Collapse
Affiliation(s)
- Katrín Blöndal
- Chemical Engineering Group, School of Engineering, Brown University, Providence, RI 02912, USA.
| | - Kirk Badger
- Chemical Engineering Group, School of Engineering, Brown University, Providence, RI 02912, USA.
| | | | - David H Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - C Franklin Goldsmith
- Chemical Engineering Group, School of Engineering, Brown University, Providence, RI 02912, USA.
| |
Collapse
|
2
|
Jaykhedkar N, Bystrický R, Sýkora M, Bučko T. Investigating the role of dispersion corrections and anharmonic effects on the phase transition in SrZrS3: A systematic analysis from AIMD free energy calculations. J Chem Phys 2024; 160:014710. [PMID: 38180257 DOI: 10.1063/5.0185319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Accepted: 12/13/2023] [Indexed: 01/06/2024] Open
Abstract
A thermally driven needle-like (NL) to distorted perovskite (DP) phase transition in SrZrS3 was investigated by means of ab initio free energy calculations accelerated by machine learning. As a first step, a systematic screening of the methods to include long-range interactions in semilocal density functional theory Perdew-Burke-Ernzerhof calculations was performed. Out of the ten correction schemes tested, the Tkatchenko-Scheffler method with iterative Hirshfeld partitioning method was found to yield the best match between calculated and experimental lattice geometries, while predicting the correct order of stability of NL and DP phases at zero temperature. This method was then used in free energy calculations, performed using several approaches, so as to determine the effect of various anharmonicity contributions, such as the anisotropic thermal lattice expansion or the thermally induced internal structure changes, on the phase transition temperature (TNP→DP). Accounting for the full anharmonicity by combining the NPT molecular dynamics data with thermodynamic integration with harmonic reference provided our best estimate of TNL→DP = 867 K. Although this result is ∼150 K lower than the experimental value, it still provides an improvement by nearly 300 K compared to the previous theoretical report by Koocher et al. [Inorg. Chem. 62, 11134-11141 (2023)].
Collapse
Affiliation(s)
- Namrata Jaykhedkar
- Laboratory of Advanced Materials, Comenius University, Ilkovičova 6, 84104 Bratislava, Slovakia
| | - Roman Bystrický
- Laboratory of Advanced Materials, Comenius University, Ilkovičova 6, 84104 Bratislava, Slovakia
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská Cesta 9, 84236 Bratislava, Slovakia
| | - Milan Sýkora
- Laboratory of Advanced Materials, Comenius University, Ilkovičova 6, 84104 Bratislava, Slovakia
| | - Tomáš Bučko
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská Cesta 9, 84236 Bratislava, Slovakia
- Department of Physical and Theoretical Chemistry, Comenius University, Ilkovičova 6, 84104 Bratislava, Slovakia
| |
Collapse
|
3
|
Yook H, Hwang J, Yeo W, Bang J, Kim J, Kim TY, Choi JS, Han JW. Design Strategies for Hydroxyapatite-Based Materials to Enhance Their Catalytic Performance and Applicability. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2204938. [PMID: 35917488 DOI: 10.1002/adma.202204938] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Revised: 07/16/2022] [Indexed: 06/15/2023]
Abstract
Hydroxyapatite (HAP) is a green catalyst that has a wide range of applications in catalysis due to its high flexibility and multifunctionality. These properties allow HAP to accommodate a large number of catalyst modifications that can selectively improve the catalytic performance in target reactions. To date, many studies have been conducted to elucidate the effect of HAP modification on the catalytic activities for various reactions. However, systematic design strategies for HAP catalysts are not established yet due to an incomplete understanding of underlying structure-activity relationships. In this review, tuning methods of HAP for improving the catalytic performance are discussed: 1) ionic composition change, 2) morphology control, 3) incorporation of other metal species, and 4) catalytic support engineering. Detailed mechanisms and effects of structural modulations on the catalytic performances for attaining the design insights of HAP catalysts are investigated. In addition, computational studies to understand catalytic reactions on HAP materials are also introduced. Finally, important areas for future research are highlighted.
Collapse
Affiliation(s)
- Hyunwoo Yook
- School of Interdisciplinary Bioscience and Bioengineering, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk, 37673, Republic of Korea
| | - Jinwoo Hwang
- Department of Chemical Engineering, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk, 37673, Republic of Korea
| | - Woonsuk Yeo
- Department of Chemical Engineering, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk, 37673, Republic of Korea
| | - Jungup Bang
- Catalyst R&D Division, LG Chem Ltd, 188, Munji-ro, Yuseong-gu, Daejeon, 34122, Republic of Korea
| | - Jaeyoung Kim
- Catalyst R&D Division, LG Chem Ltd, 188, Munji-ro, Yuseong-gu, Daejeon, 34122, Republic of Korea
| | - Tae Yong Kim
- Department of Chemical Engineering, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk, 37673, Republic of Korea
| | - Jae-Soon Choi
- Catalyst R&D Division, LG Chem Ltd, 188, Munji-ro, Yuseong-gu, Daejeon, 34122, Republic of Korea
| | - Jeong Woo Han
- School of Interdisciplinary Bioscience and Bioengineering, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk, 37673, Republic of Korea
- Department of Chemical Engineering, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk, 37673, Republic of Korea
| |
Collapse
|
4
|
Amsler J, Plessow PN, Studt F, Bučko T. Anharmonic Correction to Free Energy Barriers from DFT-Based Molecular Dynamics Using Constrained Thermodynamic Integration. J Chem Theory Comput 2023; 19:2455-2468. [PMID: 37043693 DOI: 10.1021/acs.jctc.3c00169] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/14/2023]
Abstract
For the calculation of anharmonic contributions to free energy barriers, constrained thermodynamic λ-path integration (λ-TI) from a harmonic reference force field to density functional theory is presented as an alternative to the established Blue Moon ensemble method (ξ-TI), in which free energy gradients along the reaction coordinate ξ are integrated. With good agreement in all cases, the λ-TI method is benchmarked against the ξ-TI method for several reactions, including the internal CH3 group rotation in ethane, a nucleophilic substitution of CH3Cl, a retro-Diels-Alder reaction, and a proton transfer in zeolite H-SSZ-13. An advantage of λ-TI is that one can use virtually any reference state to compute anharmonic contributions to reaction free energies or free energy barriers. This is particularly relevant for catalysis, where it is now possible to compute anharmonic corrections to the free energy of a transition state relative to any reference, for example, the most stable state of the active site and the reactants in the gas phase. This is in contrast to ξ-TI, where free energy barriers can only be computed relative to an initial state with all reactants coadsorbed. Finally, the Bennett acceptance ratio method combined with λ-TI is demonstrated to reduce the number of required integration grid points with tolerable accuracy, favoring thus λ-TI over ξ-TI in terms of computational efficiency.
Collapse
Affiliation(s)
- Jonas Amsler
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Philipp N Plessow
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Felix Studt
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institute for Chemical Technology and Polymer Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - 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
| |
Collapse
|
5
|
Berger F, Rybicki M, Sauer J. Molecular Dynamics with Chemical Accuracy─Alkane Adsorption in Acidic Zeolites. ACS Catal 2023. [DOI: 10.1021/acscatal.2c05493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Fabian Berger
- Institut für Chemie, Humboldt-Universität zu Berlin, D-10099Berlin, Germany
| | - Marcin Rybicki
- Institut für Chemie, Humboldt-Universität zu Berlin, D-10099Berlin, Germany
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität zu Berlin, D-10099Berlin, Germany
- Department of Physical and Macromolecular Chemistry & Charles University Center of Advanced Materials, Charles University, Hlavova 8, 128 43Prague 2, Czech Republic
| |
Collapse
|
6
|
Yang Y, Meng X, Zhu L, Yang J, Zhang G, Shen H, Cao X. Rapid Synthesis of Si-Rich SSZ-13 Zeolite under Fluoride-Free Conditions. Inorg Chem 2022; 61:21115-21122. [PMID: 36521022 DOI: 10.1021/acs.inorgchem.2c03749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Rapid synthesis of Si-rich (SiO2/Al2O3 > 100) SSZ-13 zeolite under fluoride-free conditions is highly desirable but still challenging. Herein, we for the first time report a rapid synthesis of all silica and aluminosilicate (SiO2/Al2O3 > 100) SSZ-13 zeolite without the addition of fluoride species. The crystallization could be fully completed at 160 °C for 4 h when the aging of the starting gel is 3 h at room temperature after the addition of a zeolite seed. The key to success is the formation of more basic building units (4- and 6-membered rings) in the initial gel with the aging time of 3 h after the addition of a zeolite seed, leading to the successful rapid synthesis of Si-rich SSZ-13 zeolite. The obtained Si-rich SSZ-13 zeolite displays high crystallinity, uniform cubic morphology with a nanoparticle feature, and a large surface area. More importantly, the obtained Si-rich SSZ-13 zeolite displays excellent performance in the adsorption of ethanol and methanol-to-olefin reaction.
Collapse
Affiliation(s)
- Yichang Yang
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou 310014, China.,College of Biological, Chemical Science and Engineering, Jiaxing University, Jiaxing 314001, China
| | - Xiangyu Meng
- College of Biological, Chemical Science and Engineering, Jiaxing University, Jiaxing 314001, China
| | - Longfeng Zhu
- College of Biological, Chemical Science and Engineering, Jiaxing University, Jiaxing 314001, China
| | - Jinghuai Yang
- College of Biological, Chemical Science and Engineering, Jiaxing University, Jiaxing 314001, China
| | - Guoliang Zhang
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou 310014, China
| | - Hongxia Shen
- College of Biological, Chemical Science and Engineering, Jiaxing University, Jiaxing 314001, China
| | - Xuebo Cao
- College of Biological, Chemical Science and Engineering, Jiaxing University, Jiaxing 314001, China
| |
Collapse
|
7
|
Blöndal K, Sargsyan K, Bross DH, Ruscic B, Goldsmith CF. Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3OH–Cu(111). ACS Catal 2022. [DOI: 10.1021/acscatal.2c04246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Katrín Blöndal
- Chemical Engineering Group, School of Engineering, Brown University, Providence, Rhode Island02912, United States
| | - Khachik Sargsyan
- Sandia National Laboratories, Livermore, California94550, United States
| | - David H. Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - C. Franklin Goldsmith
- Chemical Engineering Group, School of Engineering, Brown University, Providence, Rhode Island02912, United States
| |
Collapse
|
8
|
Bocus M, Van Speybroeck V. Insights into the Mechanism and Reactivity of Zeolite-Catalyzed Alkylphenol Dealkylation. ACS Catal 2022. [DOI: 10.1021/acscatal.2c03844] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Massimo Bocus
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052Zwijnaarde, Belgium
| | | |
Collapse
|
9
|
Li SC, Lin YC, Li YP. Comparative Analysis of Uncoupled Mode Approximations for Molecular Thermochemistry and Kinetics. J Chem Theory Comput 2022; 18:6866-6877. [PMID: 36269729 DOI: 10.1021/acs.jctc.2c00664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate prediction of thermochemistry and kinetic parameters is an important task for reaction modeling. Unfortunately, the commonly used harmonic oscillator model is often not accurate enough due to the absence of anharmonic effects. In this work, we improve the representation of an anharmonic potential energy surface (PES) using uncoupled mode (UM) approximations, which model the full-dimensional PES as a sum of one-dimensional potentials of each mode. We systematically analyze different PES sampling schemes and coordinate systems for constructing the one-dimensional potentials, and benchmark the performance of UM methods on data sets of molecular thermochemistry and kinetic properties. The results show that the accuracy of the UM approach strongly depends on how the one-dimensional potentials are defined. If one-dimensional potentials are constructed by sampling along normal mode directions (UM-N) or along the directions that minimize intermode coupling (E- and E'-optimized), the accuracies of the predicted properties are not significantly improved compared to the harmonic oscillator model. However, significant improvements can be achieved by sampling the torsional modes separately from the vibrational modes (UM-T and UM-VT). In this work, three types of coordinate systems are examined, including redundant internal coordinates (RIC), hybrid internal coordinates (HIC), and translation-rotation-internal coordinates (TRIC). The HIC and TRIC coordinate systems can outperform RIC since transition state species may contain large-amplitude interfragmentary motions that regular internal coordinates can not describe adequately. Among all the methods we examined, the activation energies and pre-exponential factors calculated using UM-VT with either TRIC or HIC best agree with the reference values. Since UM-VT requires only a number of additional single point energy calculations for each independent mode, the scaling of computational costs of UM-VT is the same as that of the standard harmonic oscillator model, making UM-VT an appealing way of calculating the thermochemistry and kinetic properties for large-size systems.
Collapse
Affiliation(s)
- Shih-Cheng Li
- Department of Chemical Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan
| | - Yen-Chun Lin
- Department of Chemical Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan
| | - Yi-Pei Li
- Department of Chemical Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan.,Taiwan International Graduate Program on Sustainable Chemical Science and Technology (TIGP-SCST), Academia Sinica, No. 128, Sec. 2, Academia Road, Taipei 11529, Taiwan
| |
Collapse
|
10
|
De Wispelaere K, Plessow PN, Studt F. Toward Computing Accurate Free Energies in Heterogeneous Catalysis: a Case Study for Adsorbed Isobutene in H-ZSM-5. ACS PHYSICAL CHEMISTRY AU 2022; 2:399-406. [PMID: 36855690 PMCID: PMC9955322 DOI: 10.1021/acsphyschemau.2c00020] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Herein, we propose a novel computational protocol that enables calculating free energies with improved accuracy by combining the best available techniques for enthalpy and entropy calculation. While the entropy is described by enhanced sampling molecular dynamics techniques, the energy is calculated using ab initio methods. We apply the method to assess the stability of isobutene adsorption intermediates in the zeolite H-SSZ-13, a prototypical problem that is computationally extremely challenging in terms of calculating enthalpy and entropy. We find that at typical operating conditions for zeolite catalysis (400 °C), the physisorbed π-complex, and not the tertiary carbenium ion as often reported, is the most stable intermediate. This method paves the way for sampling-based techniques to calculate the accurate free energies in a broad range of chemistry-related disciplines, thus presenting a big step forward toward predictive modeling.
Collapse
Affiliation(s)
- Kristof De Wispelaere
- Center
for Molecular Modeling, Ghent University, Technologiepark 46, B-9052 Ghent, Belgium,
| | - Philipp N. Plessow
- Institute
of Catalysis Research and Technology, Karlsruhe
Institute of Technology, Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany,
| | - Felix Studt
- Institute
of Catalysis Research and Technology, Karlsruhe
Institute of Technology, Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany,Institute
for Chemical Technology and Polymer Chemistry, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany,
| |
Collapse
|
11
|
Rybicki M, Sauer J. Rigid Body Approximation for the Anharmonic Description of Molecule-Surface Vibrations. J Chem Theory Comput 2022; 18:5618-5635. [PMID: 35913469 DOI: 10.1021/acs.jctc.2c00597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present an anharmonic approach for molecule-surface vibrations that employs rigid body coordinates, based on the Rodrigues rotation formula, to describe curvilinear displacements (rotations and translations) of the molecule along normal modes. These displacements are used to calculate energy data points from which one-dimensional polynomial potentials are fitted using cubic splines. In these potentials, for each of the six rigid body modes separately, one-dimensional Schrödinger equations are solved with harmonic oscillator or Fourier functions (for pure rotations) as basis sets. The anharmonic vibrational energies obtained are used to calculate partition functions and from them enthalpies, entropies, and Gibbs free energies of adsorption. Our numerical implementation has been successfully tested for Morse and cosine potentials with known analytical solutions. The methods have been applied to adsorption of CH4 on the hydroxyl group of the proton form of the chabazite zeolite (H-CHA), as well as to adsorption of CH4 and CO on the Mg2+ ions of the metal-organic framework (MOF) Mg2(dobdc). To obtain the best estimates for thermodynamic functions, we include the coupling between molecule-surface vibrations and intrasystem vibrations at the harmonic level. The calculated Gibbs free energies show that the coupling is small for CH4/H-CHA and CO/MOF (between -0.7 and +0.1 kJ/mol) but substantial for CH4/MOF (-3.4 kJ/mol). The predicted anharmonic effect on the Gibbs free energy of adsorption for CH4/H-CHA, CH4/MOF, and CO/MOF is -4.7, 0.3 ± 0.7, and -2.4 ± 0.6 kJ/mol, respectively, which results in +4.2, +0.9 ± 0.7, and -0.4 ± 0.6 kJ/mol, respectively, for the deviation from experiment. This is well within chemical accuracy limits (±4.2 kJ/mol) for the adsorption of CH4 and CO in the MOF. The larger deviation for CH4/H-CHA, at the edge of the chemical accuracy range, is most likely due to contributions from soft zeolite modes which are neglected in our approach.
Collapse
Affiliation(s)
- Marcin Rybicki
- Institut für Chemie, Humboldt-Universität zu Berlin, 10117 Berlin, Germany
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität zu Berlin, 10117 Berlin, Germany.,Department of Physical and Macromolecular Chemistry/Charles University Center of Advanced Materials, Charles University, Hlavova 8, 128 43 Prague 2, Czech Republic
| |
Collapse
|
12
|
|
13
|
Mortén M, Cordero-Lanzac T, Cnudde P, Redekop EA, Svelle S, van Speybroeck V, Olsbye U. Acidity effect on benzene methylation kinetics over substituted H-MeAlPO-5 catalysts. J Catal 2021. [DOI: 10.1016/j.jcat.2021.11.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
14
|
Borodin D, Rahinov I, Galparsoro O, Fingerhut J, Schwarzer M, Golibrzuch K, Skoulatakis G, Auerbach DJ, Kandratsenka A, Schwarzer D, Kitsopoulos TN, Wodtke AM. Kinetics of NH 3 Desorption and Diffusion on Pt: Implications for the Ostwald Process. J Am Chem Soc 2021; 143:18305-18316. [PMID: 34672570 PMCID: PMC8569812 DOI: 10.1021/jacs.1c09269] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
We report accurate time-resolved measurements of NH3 desorption from Pt(111) and Pt(332) and use these results to determine elementary rate constants for desorption from steps, from (111) terrace sites and for diffusion on (111) terraces. Modeling the extracted rate constants with transition state theory, we find that conventional models for partition functions, which rely on uncoupled degrees of freedom (DOFs), are not able to reproduce the experimental observations. The results can be reproduced using a more sophisticated partition function, which couples DOFs that are most sensitive to NH3 translation parallel to the surface; this approach yields accurate values for the NH3 binding energy to Pt(111) (1.13 ± 0.02 eV) and the diffusion barrier (0.71 ± 0.04 eV). In addition, we determine NH3's binding energy preference for steps over terraces on Pt (0.23 ± 0.03 eV). The ratio of the diffusion barrier to desorption energy is ∼0.65, in violation of the so-called 12% rule. Using our derived diffusion/desorption rates, we explain why established rate models of the Ostwald process incorrectly predict low selectivity and yields of NO under typical reactor operating conditions. Our results suggest that mean-field kinetics models have limited applicability for modeling the Ostwald process.
Collapse
Affiliation(s)
- Dmitriy Borodin
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, 37077 Goettingen, Germany.,Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany
| | - Igor Rahinov
- Department of Natural Sciences, The Open University of Israel, 4353701 Raanana, Israel
| | - Oihana Galparsoro
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain.,Kimika Fakultatea, Euskal Herriko Unibertsitatea UPV/EHU, P.K. 1072 Donostia-San Sebastián, Spain
| | - Jan Fingerhut
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, 37077 Goettingen, Germany
| | - Michael Schwarzer
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, 37077 Goettingen, Germany
| | - Kai Golibrzuch
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany
| | - Georgios Skoulatakis
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany
| | - Daniel J Auerbach
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany
| | - Alexander Kandratsenka
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany
| | - Dirk Schwarzer
- Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany
| | - Theofanis N Kitsopoulos
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, 37077 Goettingen, Germany.,Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany.,Department of Chemistry, University of Crete, 71003 Heraklion, Greece.,Institute of Electronic Structure and Laser - FORTH, 71110 Heraklion, Greece
| | - Alec M Wodtke
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, 37077 Goettingen, Germany.,Department of Dynamics at Surfaces, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Goettingen, Germany.,International Center for Advanced Studies of Energy Conversion, Georg-August University of Goettingen, Tammannstraße 6, 37077 Goettingen, Germany
| |
Collapse
|
15
|
Galimberti DR, Sauer J. Chemically Accurate Vibrational Free Energies of Adsorption from Density Functional Theory Molecular Dynamics: Alkanes in Zeolites. J Chem Theory Comput 2021; 17:5849-5862. [PMID: 34459582 PMCID: PMC8444336 DOI: 10.1021/acs.jctc.1c00519] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present a methodology
to compute, at reduced computational cost,
Gibbs free energies, enthalpies, and entropies of adsorption from
molecular dynamics. We calculate vibrational partition functions from
vibrational energies, which we obtain from the vibrational density
of states by projection on the normal modes. The use of a set of well-chosen
reference structures along the trajectories accounts for the anharmonicities
of the modes. For the adsorption of methane, ethane, and propane in
the H-CHA zeolite, we limit our treatment to a set of vibrational
modes localized at the adsorption site (zeolitic OH group) and the
alkane molecule interacting with it. Only two short trajectories (1–20
ps) are required to reach convergence (<1 kJ/mol) for the thermodynamic
functions. The mean absolute deviations from the experimentally measured
values are 2.6, 2.8, and 4.7 kJ/mol for the Gibbs free energy, the
enthalpy, and the entropy term (−TΔS),
respectively. In particular, the entropy terms show a major improvement
compared to the harmonic approximation and almost reach the accuracy
of the previous use of anharmonic frequencies obtained with curvilinear
distortions of individual modes. The thermodynamic functions so obtained
follow the trend of the experimental values for methane, ethane, and
propane, and the Gibbs free energy of adsorption at experimental conditions
is correctly predicted to change from positive for methane (5.9 kJ/mol)
to negative for ethane (−4.8 kJ/mol) and propane (−7.1
kJ/mol).
Collapse
Affiliation(s)
- Daria Ruth Galimberti
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany.,Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany
| |
Collapse
|
16
|
Lee EMY, Ludwig T, Yu B, Singh AR, Gygi F, Nørskov JK, de Pablo JJ. Neural Network Sampling of the Free Energy Landscape for Nitrogen Dissociation on Ruthenium. J Phys Chem Lett 2021; 12:2954-2962. [PMID: 33729797 DOI: 10.1021/acs.jpclett.1c00195] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction-the dissociation of molecular nitrogen on ruthenium-at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an "on-top" adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. These findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.
Collapse
Affiliation(s)
- Elizabeth M Y Lee
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Thomas Ludwig
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Boyuan Yu
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Aayush R Singh
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
| | - François Gygi
- Department of Computer Science, University of California, Davis, California 95616, United States
| | - Jens K Nørskov
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
- Department of Physics, Technical University of Denmark, Lyngby 2800, Denmark
| | - Juan J de Pablo
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
- Argonne National Laboratory, 9700 Cass Avenue, Lemont, Illinois 60439, United States
| |
Collapse
|
17
|
Dahri M, Abolmaali SS, Abedanzadeh M, Salmanpour M, Maleki R. Composition and surface chemistry engineering of graphene grafting chitosan for stimuli-responsive cancer therapy: An in-silico study. INFORMATICS IN MEDICINE UNLOCKED 2021. [DOI: 10.1016/j.imu.2021.100627] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
|