51
|
Zen A, Roch LM, Cox SJ, Hu XL, Sorella S, Alfè D, Michaelides A. Toward Accurate Adsorption Energetics on Clay Surfaces. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2016; 120:26402-26413. [PMID: 27917256 PMCID: PMC5126707 DOI: 10.1021/acs.jpcc.6b09559] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Revised: 10/31/2016] [Indexed: 05/31/2023]
Abstract
Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.
Collapse
Affiliation(s)
- Andrea Zen
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Loïc M. Roch
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Stephen J. Cox
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Xiao Liang Hu
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Sandro Sorella
- SISSA−International
School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- INFM
Democritos National Simulation Center, 34151 Trieste, Italy
| | - Dario Alfè
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United
Kingdom
| | - Angelos Michaelides
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
52
|
Arabi AA. Evaluating dispersion forces for optimization of van der Waals complexes using a non-empirical functional. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2016; 374:rsta.2016.0145. [PMID: 27698041 PMCID: PMC5052729 DOI: 10.1098/rsta.2016.0145] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 07/05/2016] [Indexed: 05/27/2023]
Abstract
Modelling dispersion interactions with traditional density functional theory (DFT) is a challenge that has been extensively addressed in the past decade. The exchange-dipole moment (XDM), among others, is a non-empirical add-on dispersion correction model in DFT. The functional PW86+PBE+XDM for exchange, correlation and dispersion, respectively, compromises an accurate functional for thermochemistry and for van der Waals (vdW) complexes at equilibrium and non-equilibrium geometries. To use this functional in optimizing vdW complexes, rather than computing single point energies, it is necessary to evaluate accurate forces. The purpose of this paper is to validate that, along the potential energy surface, the distance at which the energy is minimum is commensurate with the distance at which the forces vanish to zero. This test was validated for 10 rare gas diatomic molecules using various integration grids and different convergence criteria. It was found that the use of either convergence criterion, 10-6 or 10-8, in Gaussian09, does not affect the accuracy of computed optimal distances and binding energies. An ultra-fine grid needs to be used when computing accurate energies using generalized gradient approximation functionals.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.
Collapse
Affiliation(s)
- Alya A Arabi
- College of Natural and Health Sciences, Zayed University, Abu Dhabi, PO Box 144534, United Arab Emirates
| |
Collapse
|
53
|
Liu Y, Ojamäe L. Raman and IR Spectra of Ice Ih and Ice XI with an Assessment of DFT Methods. J Phys Chem B 2016; 120:11043-11051. [PMID: 27690444 DOI: 10.1021/acs.jpcb.6b07001] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
IR and Raman spectroscopic technology can be directly used to identify the occurrence of ferroelectric ice XI in laboratory or extraterrestrial settings. The performance of 16 different DFT methods applied on the ice Ih, VIII, IX, and XI crystal phases is evaluated. Based on a selected DFT computational scheme, the IR and Raman spectra of ice Ih and XI are derived and compared. When the spectra, both IR and Raman, of ice Ih and ice XI are compared, the librational vibrations are found to be the most affected by the proton ordering. The spectroscopic fingerprint of ice XI can be used to distinguish ferroelectric ice XI from ice Ih in the universe. Furthermore, the existence of only one kind of H-bond in ice Ih is demonstrated from the overlapping subspectra for different types of H-bonded pair configurations in 16 isomers of ice Ih, which provides an illustration to the historic debate on whether one or two kinds of H-bonds existed in ice.
Collapse
Affiliation(s)
- Yuan Liu
- Department of Chemistry, IFM, Linköping University , SE-581 83 Linköping, Sweden
| | - Lars Ojamäe
- Department of Chemistry, IFM, Linköping University , SE-581 83 Linköping, Sweden
| |
Collapse
|
54
|
Klimeš J. Lattice energies of molecular solids from the random phase approximation with singles corrections. J Chem Phys 2016; 145:094506. [DOI: 10.1063/1.4962188] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jiří Klimeš
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-18223 Prague 8, Czech Republic and Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
55
|
Curtis F, Wang X, Marom N. Effect of packing motifs on the energy ranking and electronic properties of putative crystal structures of tricyano-1,4-dithiino[c]-isothiazole. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2016; 72:562-570. [PMID: 27484377 DOI: 10.1107/s2052520616009227] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2016] [Accepted: 06/07/2016] [Indexed: 06/06/2023]
Abstract
We present an analysis of putative structures of tricyano-1,4-dithiino[c]-isothiazole (TCS3), generated within the sixth crystal structure prediction blind test. Typical packing motifs are identified and characterized in terms of distinct patterns of close contacts and regions of electrostatic and dispersion interactions. We find that different dispersion-inclusive density functional theory (DFT) methods systematically favor specific packing motifs, which may affect the outcome of crystal structure prediction efforts. The effect of crystal packing on the electronic and optical properties of TCS3 is investigated using many-body perturbation theory within the GW approximation and the Bethe-Salpeter equation (BSE). We find that a structure with Pna21 symmetry and a bilayer packing motif exhibits intermolecular bonding patterns reminiscent of π-π stacking and has markedly different electronic and optical properties than the experimentally observed P21/n structure with a cyclic dimer motif, including a narrower band gap, enhanced band dispersion and broader optical absorption. The Pna21 bilayer structure is close in energy to the observed structure and may be feasible to grow.
Collapse
Affiliation(s)
- Farren Curtis
- Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Xiaopeng Wang
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Noa Marom
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| |
Collapse
|
56
|
Engel EA, Monserrat B, Needs RJ. Vibrational effects on surface energies and band gaps in hexagonal and cubic ice. J Chem Phys 2016; 145:044703. [DOI: 10.1063/1.4959283] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edgar A. Engel
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Bartomeu Monserrat
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
- Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019, USA
| | - Richard J. Needs
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
57
|
Ambrosio F, Miceli G, Pasquarello A. Structural, Dynamical, and Electronic Properties of Liquid Water: A Hybrid Functional Study. J Phys Chem B 2016; 120:7456-70. [DOI: 10.1021/acs.jpcb.6b03876] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Francesco Ambrosio
- Chaire de Simulation à
l’Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Giacomo Miceli
- Chaire de Simulation à
l’Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Alfredo Pasquarello
- Chaire de Simulation à
l’Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
58
|
Otero-de-la-Roza A, DiLabio GA, Johnson ER. Exchange–Correlation Effects for Noncovalent Interactions in Density Functional Theory. J Chem Theory Comput 2016; 12:3160-75. [DOI: 10.1021/acs.jctc.6b00298] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- A. Otero-de-la-Roza
- National Institute for Nanotechnology, National Research Council of Canada, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Gino A. DiLabio
- National Institute for Nanotechnology, National Research Council of Canada, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Erin R. Johnson
- Department
of Chemistry, Dalhousie University, 6274 Coburg Road, Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
59
|
Fritz M, Fernández-Serra M, Soler JM. Optimization of an exchange-correlation density functional for water. J Chem Phys 2016; 144:224101. [DOI: 10.1063/1.4953081] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Michelle Fritz
- Departamento de Física de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
| | - Marivi Fernández-Serra
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794-3800, USA
| | - José M. Soler
- Departamento e Instituto de Física de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
| |
Collapse
|
60
|
Sun J, Remsing RC, Zhang Y, Sun Z, Ruzsinszky A, Peng H, Yang Z, Paul A, Waghmare U, Wu X, Klein ML, Perdew JP. Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional. Nat Chem 2016; 8:831-6. [DOI: 10.1038/nchem.2535] [Citation(s) in RCA: 538] [Impact Index Per Article: 67.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2015] [Accepted: 04/26/2016] [Indexed: 11/09/2022]
|
61
|
Sakong S, Forster-Tonigold K, Groß A. The structure of water at a Pt(111) electrode and the potential of zero charge studied from first principles. J Chem Phys 2016; 144:194701. [DOI: 10.1063/1.4948638] [Citation(s) in RCA: 104] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
| | | | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
- Helmholtz Institute Ulm (HIU), Electrochemical Energy Storage, 89069 Ulm, Germany
| |
Collapse
|
62
|
Gillan MJ, Alfè D, Michaelides A. Perspective: How good is DFT for water? J Chem Phys 2016; 144:130901. [DOI: 10.1063/1.4944633] [Citation(s) in RCA: 478] [Impact Index Per Article: 59.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Michael J. Gillan
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Dario Alfè
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
63
|
Sun H, Kang D, Dai J, Ma W, Zhou L, Zeng J. First-principles study on equation of states and electronic structures of shock compressed Ar up to warm dense regime. J Chem Phys 2016; 144:124503. [DOI: 10.1063/1.4943767] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Huayang Sun
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Dongdong Kang
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Jiayu Dai
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Wen Ma
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Liangyuan Zhou
- Academy of Ocean Science and Engineering, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Jiaolong Zeng
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
- IFSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China
| |
Collapse
|
64
|
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
|
65
|
Abstract
A combination of periodic, self-consistent density functional theory (DFT-GGA-PW91) calculations, reaction kinetics experiments on a SiO2-supported Pd catalyst, and mean-field microkinetic modeling are used to probe key aspects of H2O2 decomposition on Pd in the absence of cofeeding H2 We conclude that both Pd(111) and OH-partially covered Pd(100) surfaces represent the nature of the active site for H2O2 decomposition on the supported Pd catalyst reasonably well. Furthermore, all reaction flux in the closed catalytic cycle is predicted to flow through an O-O bond scission step in either H2O2 or OOH, followed by rapid H-transfer steps to produce the H2O and O2 products. The barrier for O-O bond scission is sensitive to Pd surface structure and is concluded to be the central parameter governing H2O2 decomposition activity.
Collapse
|
66
|
Brandenburg JG, Maas T, Grimme S. Benchmarking DFT and semiempirical methods on structures and lattice energies for ten ice polymorphs. J Chem Phys 2016; 142:124104. [PMID: 25833562 DOI: 10.1063/1.4916070] [Citation(s) in RCA: 71] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water in different phases under various external conditions is very important in bio-chemical systems and for material science at surfaces. Density functional theory methods and approximations thereof have to be tested system specifically to benchmark their accuracy regarding computed structures and interaction energies. In this study, we present and test a set of ten ice polymorphs in comparison to experimental data with mass densities ranging from 0.9 to 1.5 g/cm(3) and including explicit corrections for zero-point vibrational and thermal effects. London dispersion inclusive density functionals at the generalized gradient approximation (GGA), meta-GGA, and hybrid level as well as alternative low-cost molecular orbital methods are considered. The widely used functional of Perdew, Burke and Ernzerhof (PBE) systematically overbinds and overall provides inconsistent results. All other tested methods yield reasonable to very good accuracy. BLYP-D3(atm) gives excellent results with mean absolute errors for the lattice energy below 1 kcal/mol (7% relative deviation). The corresponding optimized structures are very accurate with mean absolute relative deviations (MARDs) from the reference unit cell volume below 1%. The impact of Axilrod-Teller-Muto (atm) type three-body dispersion and of non-local Fock exchange is small but on average their inclusion improves the results. While the density functional tight-binding model DFTB3-D3 performs well for low density phases, it does not yield good high density structures. As low-cost alternative for structure related problems, we recommend the recently introduced minimal basis Hartree-Fock method HF-3c with a MARD of about 3%.
Collapse
Affiliation(s)
- Jan Gerit Brandenburg
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Tilo Maas
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
67
|
Lin TJ, Hsing CR, Wei CM, Kuo JL. Structure prediction of the solid forms of methanol: an ab initio random structure searching approach. Phys Chem Chem Phys 2016; 18:2736-46. [DOI: 10.1039/c5cp06583f] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Liquid methanol and methanol clusters have been comprehensively studied to reveal their local structure and hydrogen bond networks.
Collapse
Affiliation(s)
- Tzu-Jen Lin
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| | - Cheng-Rong Hsing
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| | - Ching-Ming Wei
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| |
Collapse
|
68
|
Engel EA, Monserrat B, Needs RJ. Vibrational renormalisation of the electronic band gap in hexagonal and cubic ice. J Chem Phys 2015; 143:244708. [DOI: 10.1063/1.4938029] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edgar A. Engel
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Bartomeu Monserrat
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
- Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019, USA
| | - Richard J. Needs
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
69
|
Fischer M. Water adsorption in SAPO-34: elucidating the role of local heterogeneities and defects using dispersion-corrected DFT calculations. Phys Chem Chem Phys 2015; 17:25260-71. [PMID: 26352329 DOI: 10.1039/c5cp04189a] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The chabazite-type silicoaluminophosphate SAPO-34 is a promising adsorbent for applications in thermal energy storage using water adsorption-desorption cycles. In order to develop a microscopic understanding of the impact of local heterogeneities and defects on the water adsorption properties, the interaction of different models of SAPO-34 with water was studied using dispersion-corrected density-functional theory (DFT-D) calculations. In addition to SAPO-34 with isolated silicon atoms, the calculations considered models incorporating two types of heterogeneities (silicon islands, aluminosilicate domains), and two defect-containing (partially and fully desilicated) systems. DFT-D optimisations were performed for systems with small amounts of adsorbed water, in which all H2O molecules can interact with framework protons, and systems with large amounts of adsorbed water (30 H2O molecules per unit cell). At low loadings, the host-guest interaction energy calculated for SAPO-34 with isolated Si atoms amounts to approximately -90 kJ mol(-1). While the presence of local heterogeneities leads to the creation of some adsorption sites that are energetically slightly more favourable, the interaction strength is drastically reduced in systems with defects. At high water loadings, energies in the range of -70 kJ mol(-1) are obtained for all models. The DFT-D interaction energies are in good agreement with experimentally measured heats of water adsorption. A detailed analysis of the equilibrium structures was used to gain insights into the binding modes at low coverages, and to assess the extent of framework deprotonation and changes in the coordination environment of aluminium atoms at high water loadings.
Collapse
Affiliation(s)
- Michael Fischer
- Fachgebiet Kristallographie, Fachbereich Geowissenschaften, Universität Bremen, Klagenfurter Straße 2, 28359 Bremen, Germany.
| |
Collapse
|
70
|
Del Ben M, Hutter J, VandeVondele J. Probing the structural and dynamical properties of liquid water with models including non-local electron correlation. J Chem Phys 2015; 143:054506. [PMID: 26254660 DOI: 10.1063/1.4927325] [Citation(s) in RCA: 75] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaric-isothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance.
Collapse
Affiliation(s)
- Mauro Del Ben
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Joost VandeVondele
- Department of Materials, ETH Zurich, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland
| |
Collapse
|
71
|
Grimme S, Brandenburg JG, Bannwarth C, Hansen A. Consistent structures and interactions by density functional theory with small atomic orbital basis sets. J Chem Phys 2015; 143:054107. [PMID: 26254642 DOI: 10.1063/1.4927476] [Citation(s) in RCA: 523] [Impact Index Per Article: 58.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
A density functional theory (DFT) based composite electronic structure approach is proposed to efficiently compute structures and interaction energies in large chemical systems. It is based on the well-known and numerically robust Perdew-Burke-Ernzerhoff (PBE) generalized-gradient-approximation in a modified global hybrid functional with a relatively large amount of non-local Fock-exchange. The orbitals are expanded in Ahlrichs-type valence-double zeta atomic orbital (AO) Gaussian basis sets, which are available for many elements. In order to correct for the basis set superposition error (BSSE) and to account for the important long-range London dispersion effects, our well-established atom-pairwise potentials are used. In the design of the new method, particular attention has been paid to an accurate description of structural parameters in various covalent and non-covalent bonding situations as well as in periodic systems. Together with the recently proposed three-fold corrected (3c) Hartree-Fock method, the new composite scheme (termed PBEh-3c) represents the next member in a hierarchy of "low-cost" electronic structure approaches. They are mainly free of BSSE and account for most interactions in a physically sound and asymptotically correct manner. PBEh-3c yields good results for thermochemical properties in the huge GMTKN30 energy database. Furthermore, the method shows excellent performance for non-covalent interaction energies in small and large complexes. For evaluating its performance on equilibrium structures, a new compilation of standard test sets is suggested. These consist of small (light) molecules, partially flexible, medium-sized organic molecules, molecules comprising heavy main group elements, larger systems with long bonds, 3d-transition metal systems, non-covalently bound complexes (S22 and S66×8 sets), and peptide conformations. For these sets, overall deviations from accurate reference data are smaller than for various other tested DFT methods and reach that of triple-zeta AO basis set second-order perturbation theory (MP2/TZ) level at a tiny fraction of computational effort. Periodic calculations conducted for molecular crystals to test structures (including cell volumes) and sublimation enthalpies indicate very good accuracy competitive to computationally more involved plane-wave based calculations. PBEh-3c can be applied routinely to several hundreds of atoms on a single processor and it is suggested as a robust "high-speed" computational tool in theoretical chemistry and physics.
Collapse
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Jan Gerit Brandenburg
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
72
|
Delle Piane M, Corno M, Orlando R, Dovesi R, Ugliengo P. Elucidating the fundamental forces in protein crystal formation: the case of crambin. Chem Sci 2015; 7:1496-1507. [PMID: 29899894 PMCID: PMC5963673 DOI: 10.1039/c5sc03447g] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2015] [Accepted: 11/23/2015] [Indexed: 12/05/2022] Open
Abstract
This study demonstrates the feasibility of periodic all-electron hybrid density functional theory calculations in the description of protein crystals, using crambin as a test case.
Molecular simulations of proteins have been usually accomplished through empirical or semi-empirical potentials, due to the large size and inherent complexity of these biological systems. On the other hand, a theoretical description of proteins based on quantum-mechanical methods would however provide an unbiased characterization of their electronic properties, possibly offering a link between these and the ultimate biological activity. Yet, such approaches have been historically hindered by the large amount of requested computational power. Here we demonstrate the feasibility of periodic all-electron density functional theory calculations in the description of the crystal of the protein crambin (46 aminoacids), which is determined with exceptional structural accuracy. We have employed the hybrid B3LYP functional, coupled to an empirical description of London interactions (D*) to simulate the crambin crystal with an increasing amount of lattice water molecules in the cell (up to 172H2O per cell). The agreement with the experiment is good for both protein geometry and protein–water interactions. The energetics was computed to predict crystal formation energies, protein–water and protein–protein interaction energies. We studied the role of dispersion interactions which are crucial for holding the crambin crystal in place. B3LYP-D* electrostatic potential and dipole moment of crambin as well as the electronic charge flow from crambin to the solvating water molecules (0.0015e per H2O) have also been predicted. These results proved that quantum-mechanical simulations of small proteins, both free and in their crystalline state, are now feasible in a reasonable amount of time, by programs capable of exploiting high performance computing architectures, allowing the study of protein properties not easily amenable through classical force fields.
Collapse
Affiliation(s)
- Massimo Delle Piane
- Department of Chemistry and NIS Centre , University of Torino , via Pietro Giuria 7 , 10125 , Torino , Italy .
| | - Marta Corno
- Department of Chemistry and NIS Centre , University of Torino , via Pietro Giuria 7 , 10125 , Torino , Italy .
| | - Roberto Orlando
- Department of Chemistry and NIS Centre , University of Torino , via Pietro Giuria 7 , 10125 , Torino , Italy .
| | - Roberto Dovesi
- Department of Chemistry and NIS Centre , University of Torino , via Pietro Giuria 7 , 10125 , Torino , Italy .
| | - Piero Ugliengo
- Department of Chemistry and NIS Centre , University of Torino , via Pietro Giuria 7 , 10125 , Torino , Italy .
| |
Collapse
|
73
|
Al-Hamdani YS, Alfè D, von Lilienfeld OA, Michaelides A. Water on BN doped benzene: a hard test for exchange-correlation functionals and the impact of exact exchange on weak binding. J Chem Phys 2015; 141:18C530. [PMID: 25399195 DOI: 10.1063/1.4898356] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces, whereas the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B3N3H6) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry, Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| |
Collapse
|
74
|
Teeratchanan P, Hermann A. Computational phase diagrams of noble gas hydrates under pressure. J Chem Phys 2015; 143:154507. [DOI: 10.1063/1.4933371] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Pattanasak Teeratchanan
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| | - Andreas Hermann
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| |
Collapse
|
75
|
Azadi S, Foulkes WMC. Systematic study of finite-size effects in quantum Monte Carlo calculations of real metallic systems. J Chem Phys 2015; 143:102807. [DOI: 10.1063/1.4922619] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Sam Azadi
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - W. M. C. Foulkes
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| |
Collapse
|
76
|
Gaiduk AP, Gygi F, Galli G. Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals. J Phys Chem Lett 2015; 6:2902-2908. [PMID: 26267178 DOI: 10.1021/acs.jpclett.5b00901] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We determined the equilibrium density and compressibility of water and ice from first-principles molecular dynamics simulations using gradient-corrected (PBE) and hybrid (PBE0) functionals. Both functionals predicted the density of ice to be larger than that of water, by 15 (PBE) and 35% (PBE0). The PBE0 functional yielded a lower density of both ice and water with respect to PBE, leading to better agreement with experiment for ice but not for liquid water. Approximate inclusion of dispersion interactions on computed molecular-dynamics trajectories led to a substantial improvement of the PBE0 results for the density of liquid water, which, however, resulted to be slightly lower than that of ice.
Collapse
Affiliation(s)
- Alex P Gaiduk
- †Institute for Molecular Engineering, The University of Chicago, 5747 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - François Gygi
- ‡Department of Computer Science, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Giulia Galli
- †Institute for Molecular Engineering, The University of Chicago, 5747 South Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
77
|
Bankura A, Santra B, DiStasio RA, Swartz CW, Klein ML, Wu X. A systematic study of chloride ion solvation in water using van der Waals inclusive hybrid density functional theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1059959] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Arindam Bankura
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | | | - Charles W. Swartz
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Michael L. Klein
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Xifan Wu
- Department of Physics, Temple University , Philadelphia, PA, USA
| |
Collapse
|
78
|
Affiliation(s)
- Biswajit Santra
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | | | - Fausto Martelli
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | - Roberto Car
- Department of Chemistry, Princeton University , Princeton, NJ, USA
- Department of Physics, Princeton University , Princeton, NJ, USA
- Princeton Institute for the Science and Technology of Materials, Princeton University , Princeton, NJ, USA
- Program in Applied and Computational Mathematics, Princeton University , Princeton, NJ, USA
| |
Collapse
|
79
|
Soubiran F, Militzer B. MISCIBILITY CALCULATIONS FOR WATER AND HYDROGEN IN GIANT PLANETS. ACTA ACUST UNITED AC 2015. [DOI: 10.1088/0004-637x/806/2/228] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
80
|
Murphy CJ, Carrasco J, Lawton TJ, Liriano ML, Baber AE, Lewis EA, Michaelides A, Sykes ECH. Structure and energetics of hydrogen-bonded networks of methanol on close packed transition metal surfaces. J Chem Phys 2015; 141:014701. [PMID: 25005297 DOI: 10.1063/1.4882863] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Methanol is a versatile chemical feedstock, fuel source, and energy storage material. Many reactions involving methanol are catalyzed by transition metal surfaces, on which hydrogen-bonded methanol overlayers form. As with water, the structure of these overlayers is expected to depend on a delicate balance of hydrogen bonding and adsorbate-substrate bonding. In contrast to water, however, relatively little is known about the structures methanol overlayers form and how these vary from one substrate to another. To address this issue, herein we analyze the hydrogen bonded networks that methanol forms as a function of coverage on three catalytically important surfaces, Au(111), Cu(111), and Pt(111), using a combination of scanning tunneling microscopy and density functional theory. We investigate the effect of intermolecular interactions, surface coverage, and adsorption energies on molecular assembly and compare the results to more widely studied water networks on the same surfaces. Two main factors are shown to direct the structure of methanol on the surfaces studied: the surface coverage and the competition between the methanol-methanol and methanol-surface interactions. Additionally, we report a new chiral form of buckled hexamer formed by surface bound methanol that maximizes the interactions between methanol monomers by sacrificing interactions with the surface. These results serve as a direct comparison of interaction strength, assembly, and chirality of methanol networks on Au(111), Cu(111), and Pt(111) which are catalytically relevant for methanol oxidation, steam reforming, and direct methanol fuel cells.
Collapse
Affiliation(s)
- Colin J Murphy
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, USA
| | - Javier Carrasco
- CIC Energigune, Albert Einstein 48, 01510 Miñano (Álava), Spain
| | - Timothy J Lawton
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, USA
| | - Melissa L Liriano
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, USA
| | - Ashleigh E Baber
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, USA
| | - Emily A Lewis
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, USA
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Chemistry, University College London, London WC1E 6BT, United Kingdom
| | - E Charles H Sykes
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, USA
| |
Collapse
|
81
|
DiStasio RA, Santra B, Li Z, Wu X, Car R. The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water. J Chem Phys 2015; 141:084502. [PMID: 25173016 DOI: 10.1063/1.4893377] [Citation(s) in RCA: 183] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, POOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.
Collapse
Affiliation(s)
- Robert A DiStasio
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Zhaofeng Li
- Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
82
|
Cox SJ, Towler MD, Alfè D, Michaelides A. Benchmarking the performance of density functional theory and point charge force fields in their description of sI methane hydrate against diffusion Monte Carlo. J Chem Phys 2015; 140:174703. [PMID: 24811651 DOI: 10.1063/1.4871873] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
High quality reference data from diffusion Monte Carlo calculations are presented for bulk sI methane hydrate, a complex crystal exhibiting both hydrogen-bond and dispersion dominated interactions. The performance of some commonly used exchange-correlation functionals and all-atom point charge force fields is evaluated. Our results show that none of the exchange-correlation functionals tested are sufficient to describe both the energetics and the structure of methane hydrate accurately, while the point charge force fields perform badly in their description of the cohesive energy but fair well for the dissociation energetics. By comparing to ice Ih, we show that a good prediction of the volume and cohesive energies for the hydrate relies primarily on an accurate description of the hydrogen bonded water framework, but that to correctly predict stability of the hydrate with respect to dissociation to ice Ih and methane gas, accuracy in the water-methane interaction is also required. Our results highlight the difficulty that density functional theory faces in describing both the hydrogen bonded water framework and the dispersion bound methane.
Collapse
Affiliation(s)
- Stephen J Cox
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Michael D Towler
- Department of Earth Sciences, University College London Gower Street, London WC1E 6BT, United Kingdom
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| |
Collapse
|
83
|
Lobanova O, Avendaño C, Lafitte T, Müller EA, Jackson G. SAFT-γ force field for the simulation of molecular fluids: 4. A single-site coarse-grained model of water applicable over a wide temperature range. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1004804] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
84
|
Miceli G, de Gironcoli S, Pasquarello A. Isobaric first-principles molecular dynamics of liquid water with nonlocal van der Waals interactions. J Chem Phys 2015; 142:034501. [DOI: 10.1063/1.4905333] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
85
|
Fang C, Li WF, Koster RS, Klimeš J, van Blaaderen A, van Huis MA. The accurate calculation of the band gap of liquid water by means of GW corrections applied to plane-wave density functional theory molecular dynamics simulations. Phys Chem Chem Phys 2015; 17:365-75. [DOI: 10.1039/c4cp04202f] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Knowledge about the intrinsic electronic properties of water is imperative for understanding the behaviour of aqueous solutions that are used throughout biology, chemistry, physics, and industry.
Collapse
Affiliation(s)
- Changming Fang
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Wun-Fan Li
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Rik S. Koster
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Jiří Klimeš
- University of Vienna
- Faculty of Physics and Center for Computational Materials Science
- A-1090 Vienna
- Austria
| | - Alfons van Blaaderen
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| | - Marijn A. van Huis
- Debye Institute for Nanomaterials Science and Center for Extreme Matter and Emergent Phenomena
- Utrecht University
- 3584 CC Utrecht
- The Netherlands
| |
Collapse
|
86
|
Pedersen A, Wikfeldt KT, Karssemeijer L, Cuppen H, Jónsson H. Molecular reordering processes on ice (0001) surfaces from long timescale simulations. J Chem Phys 2014; 141:234706. [DOI: 10.1063/1.4903812] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
|
87
|
Gillan MJ. Many-body exchange-overlap interactions in rare gases and water. J Chem Phys 2014; 141:224106. [DOI: 10.1063/1.4903240] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
88
|
Del Ben M, VandeVondele J, Slater B. Periodic MP2, RPA, and Boundary Condition Assessment of Hydrogen Ordering in Ice XV. J Phys Chem Lett 2014; 5:4122-8. [PMID: 26278943 DOI: 10.1021/jz501985w] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Ice XV is the hydrogen-ordered form of the ice VI phase whose structure was predicted to be Cc and ferroelectric using periodic DFT approaches. However, neutron diffraction and Raman spectroscopy data show the structure to have P1̅ symmetry and to be antiferroelectric. Recent work1 using fragment-based MP2 and CCSD(T) approaches predicts the experimental structure as the ground state. We have reconsidered this problem using fully periodic MP2 and RPA approaches and find that the ferroelectric Cc structure is the lowest energy configuration. However, ubiquitously employed tinfoil boundary conditions stabilize polar structures. We suggest that ferroelectric Cc crystals can grow within a paraelectric ice VI matrix but may become unstable once a fraction of the matrix has become hydrogen-ordered. The reduction in dielectric constant causes P1̅ and other structures with small polarization to become favored, providing a possible resolution between observation and theoretical predictions.
Collapse
Affiliation(s)
- Mauro Del Ben
- †Department of Chemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Joost VandeVondele
- ‡Department of Materials, ETH Zürich, Wolfgang-Pauli-Strasse 27, CH-8093 Zürich, Switzerland
| | - Ben Slater
- §Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
89
|
Alfè D, Bartók AP, Csányi G, Gillan MJ. Analyzing the errors of DFT approximations for compressed water systems. J Chem Phys 2014; 141:014104. [DOI: 10.1063/1.4885440] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D. Alfè
- Department of Earth Sciences, UCL, London WC1E 6BT, United Kingdom
- London Centre for Nanotechnology, UCL, London WC1H 0AH, United Kingdom
- Thomas Young Centre, UCL, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, UCL, London WC1E 6BT, United Kingdom
| | - A. P. Bartók
- Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - G. Csányi
- Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - M. J. Gillan
- London Centre for Nanotechnology, UCL, London WC1H 0AH, United Kingdom
- Thomas Young Centre, UCL, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, UCL, London WC1E 6BT, United Kingdom
| |
Collapse
|
90
|
Slater B, Quigley D. Crystal nucleation: Zeroing in on ice. NATURE MATERIALS 2014; 13:670-671. [PMID: 24947779 DOI: 10.1038/nmat4017] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Affiliation(s)
- Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London W C1H 0AJ, UK
| | - David Quigley
- Department of Physics and Centre for Scientific Computing, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK
| |
Collapse
|
91
|
Macher M, Klimeš J, Franchini C, Kresse G. The random phase approximation applied to ice. J Chem Phys 2014; 140:084502. [DOI: 10.1063/1.4865748] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
92
|
Karssemeijer LJ, de Wijs GA, Cuppen HM. Interactions of adsorbed CO2 on water ice at low temperatures. Phys Chem Chem Phys 2014; 16:15630-9. [DOI: 10.1039/c4cp01622j] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
93
|
Gillan MJ, Alfè D, Bartók AP, Csányi G. First-principles energetics of water clusters and ice: A many-body analysis. J Chem Phys 2013; 139:244504. [DOI: 10.1063/1.4852182] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|