1
|
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
|
2
|
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
|
3
|
Hou Z, Lin X, Wu K, Chi H, Zhang W, Ma L, Xi Y. A density functional theory study on the adsorption of different organic sulfides on boron nitride nanosheet. RSC Adv 2023; 13:31622-31631. [PMID: 37908653 PMCID: PMC10614038 DOI: 10.1039/d3ra05718f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 10/23/2023] [Indexed: 11/02/2023] Open
Abstract
The adsorption of methanethiol (MT), thiophene (T), benzothiophene (BT), dibenzothiophene (DBT) on hexagonal boron nitride (h-BN) has been investigated by the framework of the density functional theory (DFT) calculations in this work. The prefer adsorption sites and interfacial angles of different sulfur compounds on the surface of the h-BN are investigated and analyzed. The adsorption energy results indicated that the adsorption of MT (Ead ≈ -6 kcal mol-1), T (Ead ≈ -10 kcal mol-1), BT (Ead ≈ -15 kcal mol-1), and DBT (Ead ≈ -21 kcal mol-1) on monolayer h-BN is physical interaction, and the value of Ead on bilayer h-BN is more than that on monolayer h-BN 0.05%. Adsorptive conformations show that sulfides prefer to be adsorbed on center B atoms rather than N atoms. Meanwhile, thiophene and its analogues tend to be adsorbed parallel on h-BN plane. Energy decomposition, natural population analysis (NPA), and electrostatic potential (ESP) analysis used to better understand the nature of adsorption on h-BN. van der Waals force plays a dominant role in adsorption process. Due to the π-π interactions, T, BT, and DBT have larger van der Waals forces than MT and the value of adsorption energy is negative correlated to the number of benzene rings. These findings are helpful for deeper understanding the adsorptive desulfurization mechanism and help develop better adsorbents for desulfurization in the future.
Collapse
Affiliation(s)
- Zhengjian Hou
- College of Chemistry and Chemical Engineering, China University of Petroleum (East China) Qingdao 266580 China +86-532-86981975
| | - Xufeng Lin
- College of Chemistry and Chemical Engineering, China University of Petroleum (East China) Qingdao 266580 China +86-532-86981975
- State Key Laboratory of Heavy Oil Processing, China University of Petroleum (East China) Qingdao 266580 China
| | - Ke Wu
- Changqing Engineering Design Co. Ltd., PetroChina Changqing Oilfield Company Xi'an 710000 Shanxi China
| | - Hua Chi
- College of Chemistry and Chemical Engineering, China University of Petroleum (East China) Qingdao 266580 China +86-532-86981975
| | - Wumin Zhang
- College of Chemistry and Chemical Engineering, China University of Petroleum (East China) Qingdao 266580 China +86-532-86981975
| | - Lishuang Ma
- College of Chemistry and Chemical Engineering, China University of Petroleum (East China) Qingdao 266580 China +86-532-86981975
| | - Yanyan Xi
- Advanced Chemical Engineering and Energy Materials Research Center, China University of Petroleum (East China) Qingdao 266580 China
| |
Collapse
|
4
|
Iyer GR, Rubenstein BM. Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry. J Phys Chem A 2022; 126:4636-4646. [PMID: 35820033 DOI: 10.1021/acs.jpca.2c01957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H2O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.
Collapse
Affiliation(s)
- Gopal R Iyer
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
5
|
Puzzarini C, Spada L, Alessandrini S, Barone V. The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:343002. [PMID: 32203942 DOI: 10.1088/1361-648x/ab8253] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2019] [Accepted: 03/23/2020] [Indexed: 06/10/2023]
Abstract
In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
| | - Lorenzo Spada
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Silvia Alessandrini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
6
|
Patkowski K. Recent developments in symmetry‐adapted perturbation theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1452] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Konrad Patkowski
- Department of Chemistry and Biochemistry Auburn University Auburn Alabama
| |
Collapse
|
7
|
Alessio M, Usvyat D, Sauer J. Chemically Accurate Adsorption Energies: CO and H2O on the MgO(001) Surface. J Chem Theory Comput 2018; 15:1329-1344. [DOI: 10.1021/acs.jctc.8b01122] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maristella Alessio
- Institut für Chemie, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
| | - Denis Usvyat
- 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
|
8
|
Dubecký M. Noncovalent Interactions by Fixed-Node Diffusion Monte Carlo: Convergence of Nodes and Energy Differences vs Gaussian Basis-Set Size. J Chem Theory Comput 2017; 13:3626-3635. [DOI: 10.1021/acs.jctc.7b00537] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701
03 Ostrava, Czech Republic
- ATRI, Faculty of Materials
Science and Technology, Slovak University of Technology, Paulínska
16, 917 24 Trnava, Slovakia
| |
Collapse
|
9
|
Tsatsoulis T, Hummel F, Usvyat D, Schütz M, Booth GH, Binnie SS, Gillan MJ, Alfè D, Michaelides A, Grüneis A. A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface. J Chem Phys 2017; 146:204108. [PMID: 28571392 PMCID: PMC5446292 DOI: 10.1063/1.4984048] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2017] [Accepted: 05/11/2017] [Indexed: 12/03/2022] Open
Abstract
We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.
Collapse
Affiliation(s)
- Theodoros Tsatsoulis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Felix Hummel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Denis Usvyat
- Institut für Chemie, Humboldt-Universität zu Berlin, Brook-Taylor-St. 2, D-12489 Berlin, Germany
| | - Martin Schütz
- Institut für Chemie, Humboldt-Universität zu Berlin, Brook-Taylor-St. 2, D-12489 Berlin, Germany
| | - George H Booth
- Department of Physics, King's College London, Strand, London WC2R 2LS, United Kingdom
| | - Simon S Binnie
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Michael J Gillan
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Dario Alfè
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Andreas Grüneis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| |
Collapse
|
10
|
Kebede GG, Spångberg D, Mitev PD, Broqvist P, Hermansson K. Comparing van der Waals DFT methods for water on NaCl(001) and MgO(001). J Chem Phys 2017; 146:064703. [PMID: 28201901 DOI: 10.1063/1.4971790] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
In this work, a range of van der Waals type density functionals are applied to the H2O/NaCl(001) and H2O/MgO(001) interface systems to explore the effect of an explicit dispersion treatment. The functionals we use are the self-consistent vdW functionals vdW-DF, vdW-DF2, optPBE-vdW, optB88-vdW, optB86b-vdW, and vdW-DF-cx, as well as the dispersion-corrected PBE-TS and PBE-D2 methods; they are all compared with the standard PBE functional. For both NaCl(001) and MgO(001), we find that the dispersion-flavoured functionals stabilize the water-surface interface by approximately 20%-40% compared to the PBE results. For NaCl(001), where the water molecules remain intact for all overlayers, the dominant contribution to the adsorption energy from "density functional theory dispersion" stems from the water-surface interactions rather than the water-water interactions. The optPBE-vdW and vdW-DF-cx functionals yield adsorption energies in good agreement with available experimental values for both NaCl and MgO. To probe the strengths of the perturbations of the adsorbed water molecules, we also calculated water dipole moments and found an increase up to 85% for water at the MgO(001) surface and 70% at the NaCl(001) surface, compared to the gas-phase dipole moment.
Collapse
Affiliation(s)
- Getachew G Kebede
- Department of Chemistry-Ångström, Uppsala University, Box 538, SE-751 21Uppsala, Sweden
| | - Daniel Spångberg
- Department of Chemistry-Ångström, Uppsala University, Box 538, SE-751 21Uppsala, Sweden
| | - Pavlin D Mitev
- Department of Chemistry-Ångström, Uppsala University, Box 538, SE-751 21Uppsala, Sweden
| | - Peter Broqvist
- Department of Chemistry-Ångström, Uppsala University, Box 538, SE-751 21Uppsala, Sweden
| | - Kersti Hermansson
- Department of Chemistry-Ångström, Uppsala University, Box 538, SE-751 21Uppsala, Sweden
| |
Collapse
|
11
|
|
12
|
Moc J. Adsorption, Dissociation, and Dehydrogenation of Water Monomer and Water Dimer on the Smallest 3D Aluminum Particle. The O-H Dissociation Barrier Disappears for the Dimer. J Phys Chem A 2016; 120:8725-8737. [PMID: 27715043 DOI: 10.1021/acs.jpca.6b08278] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a detailed mechanistic study on the interaction and reaction of water monomer and water dimer with the smallest 3D aluminum particle (Al6) by employing density functional and explicitly correlated coupled cluster CCSD(T)-F12 theories. Water adsorption, dissociation, and dehydrogenation are considered. For the monomer reaction, where core-valence correlation and an extrapolation to the complete basis set limit is allowed for, our coupled cluster calculations predict the O-H dissociation barrier of about 2 kcal/mol. For the dimer reaction, two distinct reaction paths are identified, initiated by forming separate dimer complexes wherein (H2O)2 adsorbs mainly via the oxygen atom of the donor H2O molecule. The key O-H dissociation transition states of the dimer reaction involve a concerted migration of two H atoms resulting in the dissociation of the donor molecule and formation of the OH-water complex adsorbed on the metal cluster's surface. The most remarkable feature of both dimer reaction energy profiles is the lack of the overall energy barrier for the (rate-determining) O-H dissociation. The hydrogen bond acceptor molecule is suggested to have an extra catalytic effect on the O-H dissociation barrier of the hydrogen bond donor molecule by removing this barrier. A similar effect on the dehydrogenation step is indicated.
Collapse
Affiliation(s)
- Jerzy Moc
- Faculty of Chemistry, Wroclaw University , F. Joliot-Curie 14, 50-383 Wroclaw, Poland
| |
Collapse
|
13
|
Affiliation(s)
- Matúš Dubecký
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
14
|
Abstract
Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
15
|
Hofer TS, Tirler AO. Combining 2d-Periodic Quantum Chemistry with Molecular Force Fields: A Novel QM/MM Procedure for the Treatment of Solid-State Surfaces and Interfaces. J Chem Theory Comput 2015; 11:5873-87. [DOI: 10.1021/acs.jctc.5b00548] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division,
Institute for General Inorganic and Theoretical Chemistry, Center
for Chemistry and Biomedicine, University of Innsbruck, Innrain
80-82, A-6020 Innsbruck, Austria
| | - Andreas O. Tirler
- Theoretical Chemistry Division,
Institute for General Inorganic and Theoretical Chemistry, Center
for Chemistry and Biomedicine, University of Innsbruck, Innrain
80-82, A-6020 Innsbruck, Austria
| |
Collapse
|
16
|
de Lara-Castells MP, Mitrushchenkov AO, Stoll H. Combining density functional and incremental post-Hartree-Fock approaches for van der Waals dominated adsorbate-surface interactions: Ag2/graphene. J Chem Phys 2015; 143:102804. [DOI: 10.1063/1.4919397] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Affiliation(s)
| | - Alexander O. Mitrushchenkov
- Université Paris-Est, Laboratoire Modélisation et Simulation Multi Echelle, MSME UMR 8208 CNRS, 5 bd Descartes, 77454 Marne-la-Vallée, France
| | - Hermann Stoll
- Institut für Theoretische Chemie, Universität Stuttgart, D-70550 Stuttgart, Germany
| |
Collapse
|
17
|
Huang Y, Goldey M, Head-Gordon M, Beran GJO. Achieving High-Accuracy Intermolecular Interactions by Combining Coulomb-Attenuated Second-Order Møller–Plesset Perturbation Theory with Coupled Kohn–Sham Dispersion. J Chem Theory Comput 2014; 10:2054-63. [DOI: 10.1021/ct5002329] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Yuanhang Huang
- Department
of Chemistry, University of California, Riverside, California 92521, United States
| | - Matthew Goldey
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Gregory J. O. Beran
- Department
of Chemistry, University of California, Riverside, California 92521, United States
| |
Collapse
|
18
|
de Lara-Castells MP, Stoll H, Mitrushchenkov AO. Assessing the Performance of Dispersionless and Dispersion-Accounting Methods: Helium Interaction with Cluster Models of the TiO2(110) Surface. J Phys Chem A 2014; 118:6367-84. [DOI: 10.1021/jp412765t] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
| | - Hermann Stoll
- Institut für Theoretische Chemie, Universität Stuttgart, D-70550 Stuttgart, Germany
| | - Alexander O. Mitrushchenkov
- Laboratoire Modélisation et Simulation Multi Echelle, MSME
UMR 8208 CNRS, Université Paris-Est, 5 bd Descartes, 77454 Marne-la-Vallée, France
| |
Collapse
|
19
|
Huang Y, Shao Y, Beran GJO. Accelerating MP2C dispersion corrections for dimers and molecular crystals. J Chem Phys 2014; 138:224112. [PMID: 23781788 DOI: 10.1063/1.4809981] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
The MP2C dispersion correction of Pitonak and Hesselmann [J. Chem. Theory Comput. 6, 168 (2010)] substantially improves the performance of second-order Møller-Plesset perturbation theory for non-covalent interactions, albeit with non-trivial computational cost. Here, the MP2C correction is computed in a monomer-centered basis instead of a dimer-centered one. When applied to a single dimer MP2 calculation, this change accelerates the MP2C dispersion correction several-fold while introducing only trivial new errors. More significantly, in the context of fragment-based molecular crystal studies, combination of the new monomer basis algorithm and the periodic symmetry of the crystal reduces the cost of computing the dispersion correction by two orders of magnitude. This speed-up reduces the MP2C dispersion correction calculation from a significant computational expense to a negligible one in crystals like aspirin or oxalyl dihydrazide, without compromising accuracy.
Collapse
Affiliation(s)
- Yuanhang Huang
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | | | | |
Collapse
|
20
|
Bučko T, Lebègue S, Hafner J, Ángyán JG. Improved Density Dependent Correction for the Description of London Dispersion Forces. J Chem Theory Comput 2013; 9:4293-9. [DOI: 10.1021/ct400694h] [Citation(s) in RCA: 149] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Tomáš Bučko
- Department
of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University,
Mlynská Dolina, SK-84215 Bratislava, Slovakia
- Institute of Inorganic Chemistry,
Slovak Academy of Sciences, Dubravska cesta 9, SK-84236 Bratislava,
Slovakia
- Fakultät
für Physik and Center for Computational
Materials Science, Universität Wien, Sensengasse, Wien 1090, Austria
| | - Sébastien Lebègue
- Université de Lorraine, CRM2, UMR 7036, Vandoeuvre-lès-Nancy, F-54506, France
- CNRS, CRM2, UMR 7036, Vandoeuvre-lès-Nancy, F-54506, France
| | - Jürgen Hafner
- Fakultät
für Physik and Center for Computational
Materials Science, Universität Wien, Sensengasse, Wien 1090, Austria
| | - János G. Ángyán
- Université de Lorraine, CRM2, UMR 7036, Vandoeuvre-lès-Nancy, F-54506, France
- CNRS, CRM2, UMR 7036, Vandoeuvre-lès-Nancy, F-54506, France
| |
Collapse
|
21
|
Beran GJO, Wen S, Nanda K, Huang Y, Heit Y. Accurate and Robust Molecular Crystal Modeling Using Fragment-Based Electronic Structure Methods. Top Curr Chem (Cham) 2013; 345:59-93. [DOI: 10.1007/128_2013_502] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|