1
|
Ye HZ, Berkelbach TC. Adsorption and vibrational spectroscopy of CO on the surface of MgO from periodic local coupled-cluster theory. Faraday Discuss 2024. [PMID: 39049598 DOI: 10.1039/d4fd00041b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of +14.7 cm-1, in excellent agreement with the experimental shift of +14.3 cm-1. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by 5-15 kJ mol-1 and predict a negative frequency shift of about -20 cm-1, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, NY 10027, USA.
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, NY 10027, USA.
- Initiative for Computational Catalysis, Flatiron Institute, New York, NY 10010, USA
| |
Collapse
|
2
|
Kanhaiya K, Nathanson M, In 't Veld PJ, Zhu C, Nikiforov I, Tadmor EB, Choi YK, Im W, Mishra RK, Heinz H. Accurate Force Fields for Atomistic Simulations of Oxides, Hydroxides, and Organic Hybrid Materials up to the Micrometer Scale. J Chem Theory Comput 2023; 19:8293-8322. [PMID: 37962992 DOI: 10.1021/acs.jctc.3c00750] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
The simulation of metals, oxides, and hydroxides can accelerate the design of therapeutics, alloys, catalysts, cement-based materials, ceramics, bioinspired composites, and glasses. Here we introduce the INTERFACE force field (IFF) and surface models for α-Al2O3, α-Cr2O3, α-Fe2O3, NiO, CaO, MgO, β-Ca(OH)2, β-Mg(OH)2, and β-Ni(OH)2. The force field parameters are nonbonded, including atomic charges for Coulomb interactions, Lennard-Jones (LJ) potentials for van der Waals interactions with 12-6 and 9-6 options, and harmonic bond stretching for hydroxide ions. The models outperform DFT calculations and earlier atomistic models (Pedone, ReaxFF, UFF, CLAYFF) up to 2 orders of magnitude in reliability, compatibility, and interpretability due to a quantitative representation of chemical bonding consistent with other compounds across the periodic table and curated experimental data for validation. The IFF models exhibit average deviations of 0.2% in lattice parameters, <10% in surface energies (to the extent known), and 6% in bulk moduli relative to experiments. The parameters and models can be used with existing parameters for solvents, inorganic compounds, organic compounds, biomolecules, and polymers in IFF, CHARMM, CVFF, AMBER, OPLS-AA, PCFF, and COMPASS, to simulate bulk oxides, hydroxides, electrolyte interfaces, and multiphase, biological, and organic hybrid materials at length scales from atoms to micrometers. The nonbonded character of the models also enables the analysis of mixed oxides, glasses, and certain chemical reactions, and well-performing nonbonded models for silica phases, SiO2, are introduced. Automated model building is available in the CHARMM-GUI Nanomaterial Modeler. We illustrate applications of the models to predict the structure of mixed oxides, and energy barriers of ion migration, as well as binding energies of water and organic molecules in outstanding agreement with experimental data and calculations at the CCSD(T) level. Examples of model building for hydrated, pH-sensitive oxide surfaces to simulate solid-electrolyte interfaces are discussed.
Collapse
Affiliation(s)
- Krishan Kanhaiya
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| | - Michael Nathanson
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| | - Pieter J In 't Veld
- BASF SE, Molecular Modeling & Drug Discovery, Carl Bosch Str. 38, 67056 Ludwigshafen, Germany
| | - Cheng Zhu
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| | - Ilia Nikiforov
- Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Ellad B Tadmor
- Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Yeol Kyo Choi
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Ratan K Mishra
- BASF SE, Molecular Modeling & Drug Discovery, Carl Bosch Str. 38, 67056 Ludwigshafen, Germany
| | - Hendrik Heinz
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
3
|
Shi B, Zen A, Kapil V, Nagy PR, Grüneis A, Michaelides A. Many-Body Methods for Surface Chemistry Come of Age: Achieving Consensus with Experiments. J Am Chem Soc 2023; 145:25372-25381. [PMID: 37948071 PMCID: PMC10683001 DOI: 10.1021/jacs.3c09616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Revised: 10/15/2023] [Accepted: 10/17/2023] [Indexed: 11/12/2023]
Abstract
The adsorption energy of a molecule onto the surface of a material underpins a wide array of applications, spanning heterogeneous catalysis, gas storage, and many more. It is the key quantity where experimental measurements and theoretical calculations meet, with agreement being necessary for reliable predictions of chemical reaction rates and mechanisms. The prototypical molecule-surface system is CO adsorbed on MgO, but despite intense scrutiny from theory and experiment, there is still no consensus on its adsorption energy. In particular, the large cost of accurate many-body methods makes reaching converged theoretical estimates difficult, generating a wide range of values. In this work, we address this challenge, leveraging the latest advances in diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] to obtain accurate predictions for CO on MgO. These reliable theoretical estimates allow us to evaluate the inconsistencies in published temperature-programed desorption experiments, revealing that they arise from variations in employed pre-exponential factors. Utilizing this insight, we derive new experimental estimates of the (electronic) adsorption energy with a (more) precise pre-exponential factor. As a culmination of all of this effort, we are able to reach a consensus between multiple theoretical calculations and multiple experiments for the first time. In addition, we show that our recently developed cluster-based CCSD(T) approach provides a low-cost route toward achieving accurate adsorption energies. This sets the stage for affordable and reliable theoretical predictions of chemical reactions on surfaces to guide the realization of new catalysts and gas storage materials.
Collapse
Affiliation(s)
- Benjamin
X. Shi
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Andrea Zen
- Dipartimento
di Fisica Ettore Pancini, Università
di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
- Department
of Earth Sciences, University College London, Gower Street, WC1E 6BT London, U.K.
| | - Venkat Kapil
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Andreas Grüneis
- Institute
for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| |
Collapse
|
4
|
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
|
5
|
Röckert A, Kullgren J, Hermansson K. Predicting Frequency from the External Chemical Environment: OH Vibrations on Hydrated and Hydroxylated Surfaces. J Chem Theory Comput 2022; 18:7683-7694. [PMID: 36458913 DOI: 10.1021/acs.jctc.2c00135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
Robust correlation curves are essential to decipher structural information from IR-vibrational spectra. However, for surface-adsorbed water and hydroxides, few such correlations have been presented in the literature. In this paper, OH vibrational frequencies are correlated against 12 structural descriptors representing the quantum mechanical or geometrical environment, focusing on those external to the vibrating molecule. A nonbiased fitting procedure based on Gaussian process regression (GPR) was used alongside simple analytical functional forms. The training data consist of 217 structurally unique OH groups from 38 water/metal oxide interface systems for MgO, CaO and CeO2, all optimized at the DFT level, and the fully anharmonic and uncoupled OH vibrational signatures were calculated. Among our results, we find the following: (i) The intermolecular R(H···O) hydrogen bond distance is particularly strong, indicating the primary cause of the frequency shift. (ii) Similarly, the electric field along the H-bond vector is also a good descriptor. (iii) Highly detailed machine learning descriptors (ACSF, SOAP) are less intuitive but were found to be more capable descriptors. (iv) Combinations of geometric and QM descriptors give the best predictions, supplying complementary information.
Collapse
Affiliation(s)
- Andreas Röckert
- Department of Chemistry, Ångström Laboratory, Uppsala University, Uppsala751 21, Sweden
| | - Jolla Kullgren
- Department of Chemistry, Ångström Laboratory, Uppsala University, Uppsala751 21, Sweden
| | - Kersti Hermansson
- Department of Chemistry, Ångström Laboratory, Uppsala University, Uppsala751 21, Sweden
| |
Collapse
|
6
|
Adsorption energies on transition metal surfaces: towards an accurate and balanced description. Nat Commun 2022; 13:6853. [PMID: 36369277 PMCID: PMC9652424 DOI: 10.1038/s41467-022-34507-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Accepted: 10/27/2022] [Indexed: 11/13/2022] Open
Abstract
Density functional theory predictions of binding energies and reaction barriers provide invaluable data for analyzing chemical transformations in heterogeneous catalysis. For high accuracy, effects of band structure and coverage, as well as the local bond strength in both covalent and non-covalent interactions, must be reliably described and much focus has been put on improving functionals to this end. Here, we show that a correction from higher-level calculations on small metal clusters can be applied to improve periodic band structure adsorption energies and barriers. We benchmark against 38 reliable experimental covalent and non-covalent adsorption energies and five activation barriers with mean absolute errors of 2.2 kcal mol-1, 2.7 kcal mol-1, and 1.1 kcal mol-1, respectively, which are lower than for functionals widely used and tested for surface science evaluations, such as BEEF-vdW and RPBE.
Collapse
|
7
|
Lu Y, Zheng F, Lan Q, Schnedler M, Ebert P, Dunin-Borkowski RE. Counting Point Defects at Nanoparticle Surfaces by Electron Holography. NANO LETTERS 2022; 22:6936-6941. [PMID: 36041122 DOI: 10.1021/acs.nanolett.2c01510] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Metal oxide nanoparticles exhibit outstanding catalytic properties, believed to be related to the presence of oxygen vacancies at the particle's surface. However, little quantitative information is known about concentrations of point defects inside and at surfaces of these nanoparticles, due to the challenges in achieving an atomically resolved experimental access. By employing off-axis electron holography, we demonstrate, using MgO nanoparticles as an example, a methodology that discriminates between mobile charge induced by electron beam irradiation and immobile charge associated with deep traps induced by point defects as well as distinguishes between bulk and surface point defects. Counting the immobile charge provides a quantification of the concentration of F2+ centers induced by oxygen vacancies at the MgO nanocube surfaces.
Collapse
Affiliation(s)
- Yan Lu
- Ernst Ruska-Centre for Microscopy and Spectroscopy with Electrons (ER-C 1) and Peter Grünberg Institute (PGI 5), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
- Beijing Key Lab and Institute of Microstructure and Properties of Advanced Materials, Beijing University of Technology, Beijing 100124, China
| | - Fengshan Zheng
- Ernst Ruska-Centre for Microscopy and Spectroscopy with Electrons (ER-C 1) and Peter Grünberg Institute (PGI 5), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Qianqian Lan
- Ernst Ruska-Centre for Microscopy and Spectroscopy with Electrons (ER-C 1) and Peter Grünberg Institute (PGI 5), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Michael Schnedler
- Ernst Ruska-Centre for Microscopy and Spectroscopy with Electrons (ER-C 1) and Peter Grünberg Institute (PGI 5), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Philipp Ebert
- Ernst Ruska-Centre for Microscopy and Spectroscopy with Electrons (ER-C 1) and Peter Grünberg Institute (PGI 5), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Rafal E Dunin-Borkowski
- Ernst Ruska-Centre for Microscopy and Spectroscopy with Electrons (ER-C 1) and Peter Grünberg Institute (PGI 5), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| |
Collapse
|
8
|
Mitra A, Hermes MR, Cho M, Agarawal V, Gagliardi L. Periodic Density Matrix Embedding for CO Adsorption on the MgO(001) Surface. J Phys Chem Lett 2022; 13:7483-7489. [PMID: 35939641 PMCID: PMC9393885 DOI: 10.1021/acs.jpclett.2c01915] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Accepted: 08/01/2022] [Indexed: 05/19/2023]
Abstract
The adsorption of simple gas molecules to metal oxide surfaces is a primary step in many heterogeneous catalysis applications. Quantum chemical modeling of these reactions is a challenge in terms of both cost and accuracy, and quantum-embedding methods are promising, especially for localized chemical phenomena. In this work, we employ density matrix embedding theory (DMET) for periodic systems to calculate the adsorption energy of CO to the MgO(001) surface. Using coupled-cluster theory with single and double excitations and second-order Møller-Plesset perturbation theory as quantum chemical solvers, we perform calculations with embedding clusters up to 266 electrons in 306 orbitals, with the largest embedding models agreeing to within 1.2 kcal/mol of the non-embedding references. Moreover, we present a memory-efficient procedure of storing and manipulating electron repulsion integrals in the embedding space within the framework of periodic DMET.
Collapse
Affiliation(s)
- Abhishek Mitra
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Minsik Cho
- Department
of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Valay Agarawal
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Laura Gagliardi
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Argonne
National Laboratory 9700
South Cass Avenue Lemont, Illinois 60439, United
States
- E-mail:
| |
Collapse
|
9
|
Gao Y, Dong L, Huang L, Huang Z, Li F, Zhang H, Zhang S. A First-Principles Study on the Hydration Behavior of (MgO)n Clusters and the Effect Mechanism of Anti-Hydration Agents. MATERIALS 2022; 15:ma15103521. [PMID: 35629547 PMCID: PMC9146619 DOI: 10.3390/ma15103521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Revised: 05/03/2022] [Accepted: 05/09/2022] [Indexed: 02/04/2023]
Abstract
Magnesia-based refractory is widely used in high-temperature industries; its easy hydration is, however, a key concern in refractory processing. Understanding the hydration mechanism of MgO will help in solving its hydration problem. Herein, the hydration behavior of (MgO)n (n = 1–6) at the molecular level and the effect mechanisms of several anti-hydration agents on the hydration of (MgO)4 were investigated with first-principles calculations. The results indicated that the following: (1) The smaller the (MgO)n cluster size, the more favorable the hydration of MgO and the tendency to convert into Mg(OH)2 crystal; (2) Anti-hydration agents can coordinate with the unsaturated Mg atom of (MgO)4 to form a bond, increasing the coordination number of Mg, thus reducing its activity when reacting with H2O; (3) The greater the number of −COOH groups and the longer the chain length in the anti-hydration agents, the better its effect of inhibiting the hydration of MgO. These findings could enhance the understanding of the mechanism of hydration of MgO and provide theoretical guidance for the design of novel anti-hydration agents.
Collapse
Affiliation(s)
- Yu Gao
- The State Key Laboratory of Refractories and Metallurgy, College of Materials and Metallurgy, Wuhan University of Science and Technology, Wuhan 430081, China; (Y.G.); (L.D.); (Z.H.); (F.L.)
| | - Long Dong
- The State Key Laboratory of Refractories and Metallurgy, College of Materials and Metallurgy, Wuhan University of Science and Technology, Wuhan 430081, China; (Y.G.); (L.D.); (Z.H.); (F.L.)
| | - Liang Huang
- The State Key Laboratory of Refractories and Metallurgy, College of Materials and Metallurgy, Wuhan University of Science and Technology, Wuhan 430081, China; (Y.G.); (L.D.); (Z.H.); (F.L.)
- Correspondence: (L.H.); (H.Z.)
| | - Zhong Huang
- The State Key Laboratory of Refractories and Metallurgy, College of Materials and Metallurgy, Wuhan University of Science and Technology, Wuhan 430081, China; (Y.G.); (L.D.); (Z.H.); (F.L.)
| | - Faliang Li
- The State Key Laboratory of Refractories and Metallurgy, College of Materials and Metallurgy, Wuhan University of Science and Technology, Wuhan 430081, China; (Y.G.); (L.D.); (Z.H.); (F.L.)
| | - Haijun Zhang
- The State Key Laboratory of Refractories and Metallurgy, College of Materials and Metallurgy, Wuhan University of Science and Technology, Wuhan 430081, China; (Y.G.); (L.D.); (Z.H.); (F.L.)
- Correspondence: (L.H.); (H.Z.)
| | - Shaowei Zhang
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter EX4 4QF, UK;
| |
Collapse
|
10
|
Gardeh MG, Kistanov AA, Nguyen H, Manzano H, Cao W, Kinnunen P. Exploring Mechanisms of Hydration and Carbonation of MgO and Mg(OH) 2 in Reactive Magnesium Oxide-Based Cements. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:6196-6206. [PMID: 35449521 PMCID: PMC9014411 DOI: 10.1021/acs.jpcc.1c10590] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/15/2021] [Revised: 03/16/2022] [Indexed: 06/14/2023]
Abstract
Reactive magnesium oxide (MgO)-based cement (RMC) can play a key role in carbon capture processes. However, knowledge on the driving forces that control the degree of carbonation and hydration and rate of reactions in this system remains limited. In this work, density functional theory-based simulations are used to investigate the physical nature of the reactions taking place during the fabrication of RMCs under ambient conditions. Parametric indicators such as adsorption energies, charge transfer, electron localization function, adsorption/dissociation energy barriers, and the mechanisms of interaction of H2O and CO2 molecules with MgO and brucite (Mg(OH)2) clusters are considered. The following hydration and carbonation interactions relevant to RMCs are evaluated: (i) carbonation of MgO, (ii) hydration of MgO, carbonation of hydrated MgO, (iii) carbonation of Mg(OH)2, (iv) hydration of Mg(OH)2, and (v) hydration of carbonated Mg(OH)2. A comparison of the energy barriers and reaction pathways of these mechanisms shows that the carbonation of MgO is hindered by the presence of H2O molecules, while the carbonation of Mg(OH)2 is hindered by the formation of initial carbonate and hydrate layers as well as presence of excessed H2O molecules. To compare these finding to bulk mineral surfaces, the interactions of the CO2 and H2O molecules with the MgO(001) and Mg(OH)2 (001) surfaces are studied. Therefore, this work presents deep insights into the physical nature of the reactions and the mechanisms involved in hydrated magnesium carbonates production that can be beneficial for its development.
Collapse
Affiliation(s)
- Mina Ghane Gardeh
- Fibre
and Particle Engineering Research Unit, University of Oulu, Pentti Kaiteran katu 1, 90014 Oulu, Finland
| | - Andrey A. Kistanov
- Nano
and Molecular Systems Research Unit, University
of Oulu, Pentti Kaiteran katu 1, 90014 Oulu, Finland
| | - Hoang Nguyen
- Fibre
and Particle Engineering Research Unit, University of Oulu, Pentti Kaiteran katu 1, 90014 Oulu, Finland
| | - Hegoi Manzano
- Departament
of Condensed Matter Physics, University
of the Basque Country (UPV/EHU), Barrio Sarriena, s/n, 48940 Leioa, Spain
| | - Wei Cao
- Nano
and Molecular Systems Research Unit, University
of Oulu, Pentti Kaiteran katu 1, 90014 Oulu, Finland
| | - Paivo Kinnunen
- Fibre
and Particle Engineering Research Unit, University of Oulu, Pentti Kaiteran katu 1, 90014 Oulu, Finland
| |
Collapse
|
11
|
Christlmaier EM, Kats D, Alavi A, Usvyat D. Full Configuration Interaction Quantum Monte Carlo treatment of fragments embedded in a periodic mean field. J Chem Phys 2022; 156:154107. [DOI: 10.1063/5.0084040] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an embedded fragment approach for high-level quantum chemical calculations on local features in periodic systems. The fragment is defined as a set of localized orbitals (occupied and virtual) corresponding to a converged periodic Hartree-Fock solution. These orbitals serve as the basis for the in-fragment post-Hartree Fock treatment. The embedding field for the fragment, consisting of the Coulomb and exchange potential from the rest of the crystal, is included in the fragment's one-electron Hamiltonian. As an application of the embedded fragment approach we investigate the performanceof full configuration interaction quantum Monte Carlo (FCIQMC) with the adaptive shift. As the orbital choice we use the natural orbitals from the distinguishable cluster method with singles and doubles. FCIQMC is a stochastic approximation to the full CI method and can be routinely applied to much larger active spaces than the latter. This makes this method especially attractive in the context of open shell defects in crystals, where fragments of adequate size can be ratherlarge. As a test case we consider dissociation of a fluorine atom from a fluorographane surface. This process poses a challenge for high-level electronic structure models as both the static and dynamic correlations are essential here. Furthermore the active space for an adequate fragment (32 electrons in 173 orbitals) is already quite large even for FCIQMC. Despite this, FCIQMC delivers accurate dissociation and total energies.
Collapse
Affiliation(s)
| | - Daniel Kats
- Max-Planck-Institute for Solid State Research, Germany
| | - Ali Alavi
- Max-Planck-Institute for Solid State Research, Germany
| | - Denis Usvyat
- Institute of Chemistry, Humboldt University of Berlin, Germany
| |
Collapse
|
12
|
Mullan T, Maschio L, Saalfrank P, Usvyat D. Reaction barriers on non-conducting surfaces beyond periodic local MP2: Diffusion of hydrogen on \ce{\alpha-Al2O3}(0001) as a test case. J Chem Phys 2022; 156:074109. [DOI: 10.1063/5.0082805] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Lorenzo Maschio
- Dipartimento di Chimica, Università degli Studi di Torino, Italy
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam Institut für Chemie, Germany
| | - Denis Usvyat
- Institute of Chemistry, Humboldt University of Berlin, Germany
| |
Collapse
|
13
|
Hansen AS, Aurbakken E, Pedersen TB. Smooth potential-energy surfaces in fragmentation-based local correlation methods for periodic systems. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1896046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- A. S. Hansen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - E. Aurbakken
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - T. B. Pedersen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| |
Collapse
|
14
|
Amsler J, Sarma BB, Agostini G, Prieto G, Plessow PN, Studt F. Prospects of Heterogeneous Hydroformylation with Supported Single Atom Catalysts. J Am Chem Soc 2020; 142:5087-5096. [DOI: 10.1021/jacs.9b12171] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Jonas Amsler
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Bidyut B. Sarma
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Giovanni Agostini
- ALBA Synchrotron Light Source, Carrer de la Llum 2-26, 08290 Cerdanyola del Vallès, Barcelona, Spain
| | - Gonzalo Prieto
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
- ITQ Instituto de Tecnologı́a Quı́mica, Universitat Politècnica de València-Consejo Superior de Investigaciones Científicas (UPV-CSIC), Av. Los Naranjos s/n, 46022 Valencia, Spain
| | - 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
| |
Collapse
|
15
|
Rehak FR, Piccini G, Alessio M, Sauer J. Including dispersion in density functional theory for adsorption on flat oxide surfaces, in metal-organic frameworks and in acidic zeolites. Phys Chem Chem Phys 2020; 22:7577-7585. [PMID: 32227013 DOI: 10.1039/d0cp00394h] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We examine the performance of nine commonly used methods for including dispersion interactions in density functional theory (DFT): three different parametrizations of damped 1/Rn terms (n = 6, 8, …) added to the DFT energy (Grimme's D2 and D3 parameterizations as well as that of Tkatchenko and Scheffler), three different implementations of the many-body dispersion approach (MBD, MBD/HI and MBD/FI), the density-dependent energy correction, called dDsC, and two "first generation" van der Waals density functionals, revPBE-vdW and optB86b-vdW. As test set we use eight molecule-surface systems for which agreement has been reached between experiment and hybrid QM:QM calculations within chemical accuracy limits (±4.2 kJ mol-1). It includes adsorption of carbon monoxide and dioxide in the Mg2(2,5-dioxido-1,4-benzenedicarboxylate) metal-organic framework (Mg-MOF-74, CPO-27-Mg), adsorption of carbon monoxide as well as of monolayers of methane and ethane on the MgO(001) surface, as well as adsorption of methane, ethane and propane in H-chabazite (H-CHA). D2 with Ne parameters for Mg2+, D2(Ne), MBD/HI and MBD/FI perform best. With the PBE functional, the mean unsigned errors are 6.1, 5.6 and 5.4 kJ mol-1, respectively.
Collapse
Affiliation(s)
- Florian R Rehak
- Institut für Chemie, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany.
| | - GiovanniMaria Piccini
- Institut für Chemie, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany.
| | - Maristella Alessio
- Institut für Chemie, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany.
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany.
| |
Collapse
|
16
|
Sauer J. Ab Initio Calculations for Molecule-Surface Interactions with Chemical Accuracy. Acc Chem Res 2019; 52:3502-3510. [PMID: 31765121 DOI: 10.1021/acs.accounts.9b00506] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Atomistic understanding of complex surface phenomena such as heterogeneous catalysis or storage and separation of energy-relevant gases in nanoporous materials (zeolites; metal-organic frameworks, MOFs) requires knowledge about reaction energies and energy barriers for elementary steps. This is difficult to obtain from experiment since the number of possible chemical, adsorption/desorption, and diffusion steps coupled to complex reaction networks is large, and so is the number of possible surface sites. Here is an important role of quantum chemistry which can provide rate and equilibrium constants for elementary steps "ab initio." To be useful, the predictions have to reach chemical accuracy (4 kJ/mol) which is difficult to achieve because realistic models of the surface systems may comprise of the order of a thousand atoms. While density functional theory (DFT) as a rule cannot be trusted to yield results within chemical accuracy limits, methods that are accurate enough (Coupled Cluster with Single, Double, and perturbative Triple Substitution, CCSD(T)) cannot be applied because of their exponential scaling with system size. This Account presents a hybrid high-level-low-level quantum method that combines DFT with dispersion for the full periodic system with second order Møller-Plesset perturbation theory (MP2) for the reaction site within a mechanical embedding scheme. In addition, to check if MP2 is accurate enough, we calculate Coupled Cluster (CC) corrections with Single, Double, and perturbatively treated Triple substitutions (CCSD(T)) for sufficiently small models of the reaction site. This multilevel hybrid MP2:DFT-D+ΔCC method is shown to yield chemical accuracy for a set of 12 molecule-surface interaction systems for which reliable experimental data are available. For CO/MgO(001), the history of the experiment-theory comparison illustrates two problems: (i) Do experiment and theory look at the same surface site? (ii) Does theory calculate the same quantity as derived from experiment? The hybrid MP2:DFT-D+ΔCC data set generated includes the MgO(001) surface, the Mg2(dobdc) metal-organic framework, and the proton forms of the CHA and MFI zeolites interacting with the H2, N2, CO, CO2, CH4, and C2H6 molecules. It serves two purposes. First, it will be useful for testing density functionals with respect to their performance for molecule-surface interactions. Second, it establishes the hybrid MP2:DFT-D+ΔCC method as a reliable and powerful tool for ab initio predictions of adsorption and reaction energies as well as energy barriers when testing reaction mechanisms. For adsorption of small molecules in MOFs, isotherm predictions have reached a level of accuracy that deviations between theoretical predictions and experiments indicate sample imperfections. For elementary steps of the industrially important methanol-to-olefin process, our hybrid MP2:PBE+D+ΔCC calculations yield rate constants in agreement with experiment within chemical accuracy limits, finally achieving for molecule-surface reactions which was possible so hitherto only for gas phase reactions involving not more than 10 atoms.
Collapse
Affiliation(s)
- Joachim Sauer
- Institute of Chemistry, Humboldt University, 10117 Berlin, Germany
| |
Collapse
|
17
|
Maul J, Spoto G, Mino L, Erba A. Elucidating the structure and dynamics of CO ad-layers on MgO surfaces. Phys Chem Chem Phys 2019; 21:26279-26283. [PMID: 31687690 DOI: 10.1039/c9cp05418a] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The combination of quantum-mechanical simulations and infrared absorption spectroscopy measurements provides a clear picture for a long standing puzzle in surface science: the actual structure and vibrational dynamics of the low-temperature ordered CO monolayer adsorbed on (001) MgO surfaces. The equilibrium structure of the commensurate (4 × 2) adsorbed phase consists of three CO molecules per primitive cell (surface coverage of 75%) located at two inequivalent sites: one molecule seats upright on top of a Mg site while two molecules, tilted off the normal to the surface, are symmetrically positioned relative to the upright one with anti-parallel projections on the surface. This configuration, long believed to be incompatible with measured polarization infrared spectra, is shown to reproduce all observed spectral features, including a new, unexpected one: the vanishing anharmonicity of CO stretching modes in the monolayer.
Collapse
Affiliation(s)
- Jefferson Maul
- Dipartimento di Chimica, Università di Torino, via Giuria 7, 10125, Torino, Italy.
| | | | | | | |
Collapse
|
18
|
Kundu A, Sillar K, Sauer J. Predicting adsorption selectivities from pure gas isotherms for gas mixtures in metal-organic frameworks. Chem Sci 2019; 11:643-655. [PMID: 34123036 PMCID: PMC8146500 DOI: 10.1039/c9sc03008e] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
We perform Grand Canonical Monte Carlo simulations on a lattice of Mg2+ sites (GCMC) for adsorption of four binary A/B mixtures, CH4/N2, CO/N2, CO2/N2, and CO2/CH4, in the metal–organic framework Mg2(2,5-dioxidobenzedicarboxylate), also known as CPO-27–Mg or Mg–MOF-74. We present a mean field co-adsorption isotherm model and show that its predictions agree with the GCMC results if the same quantum chemical ab initio data are used for Gibbs free energies of adsorption at the individual sites and for lateral interaction energies between the same, A⋯A and B⋯B, and unlike, A⋯B, adsorbed molecules. We use both approaches to test the assumption underlying Ideal Adsorbed Solution Theory (IAST), namely approximating A⋯B interaction energies as the arithmetic mean of A⋯A and B⋯B interaction energies. While IAST works well for mixtures with weak lateral interactions, CH4/N2 and CO/N2, the deviations are large for mixtures with stronger lateral interactions, CO2/N2 and CO2/CH4. Motivated by the theory of London dispersion forces, we propose use of the geometric mean instead of the arithmetic mean and achieve substantial improvements. For CO2/CH4, the lateral interactions become anisotropic. To include this in the geometric mean co-adsorption model, we introduce an anisotropy factor. We propose a protocol, named co-adsorption mean field theory (CAMT), for co-adsorption selectivity prediction from known (experiment or simulation) pure component isotherms which is similar to the IAST protocol but uses the geometric mean to approximate mixed pair interaction energies and yields improved results for non-ideal mixtures. A new mixing rule (geometric mean) is proposed with substantial improvements compared to the widely used ideal adsorbed solution theory for adsorbates with strong lateral interactions.![]()
Collapse
Affiliation(s)
- Arpan Kundu
- Humboldt Universität zu Berlin, Institut für Chemie Unter den Linden 6 10099 Berlin Germany
| | - Kaido Sillar
- Humboldt Universität zu Berlin, Institut für Chemie Unter den Linden 6 10099 Berlin Germany .,University of Tartu, Institute of Chemistry Ravila 14a 50411 Tartu Estonia
| | - Joachim Sauer
- Humboldt Universität zu Berlin, Institut für Chemie Unter den Linden 6 10099 Berlin Germany
| |
Collapse
|
19
|
Kropp T, Mavrikakis M. Brønsted–Evans–Polanyi relation for CO oxidation on metal oxides following the Mars–van Krevelen mechanism. J Catal 2019. [DOI: 10.1016/j.jcat.2019.08.002] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
20
|
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
|