1
|
Romanelli G, Andreani C, Bocedi A, Senesi R. Quantum motion of oxygen and hydrogen in water: Atomic and total kinetic energy across melting from neutron scattering measurements. J Chem Phys 2024; 160:234503. [PMID: 38884402 DOI: 10.1063/5.0211165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2024] [Accepted: 05/22/2024] [Indexed: 06/18/2024] Open
Abstract
We provide a concurrent measurement of the hydrogen and oxygen nuclear kinetic energies in the water molecule across melting at 270 K in the solid phase and 276 K in the liquid phase. Experimental values are obtained by analyzing the neutron Compton profiles of each atomic species in a deep inelastic neutron scattering experiment. The concurrent measurement of the atom kinetic energy of both hydrogen and oxygen allows the estimate of the total kinetic energy per molecule due to the motion of nuclei, specifically 35.3 ± 0.8 and 34.8 ± 0.8 kJ/mol for the solid and liquid phases, respectively. Such a small difference supports results from ab initio simulations and phenomenological models from the literature on the mechanism of competing quantum effects across the phase change. Despite the experimental uncertainties, the results are consistent with the trend from state-of-the-art computer simulations, whereby the atom and molecule kinetic energies in the liquid phase would be slightly lower than in the solid phase. Moreover, the small change of nuclear kinetic energy across melting can be used to simplify the calculation of neutron-related environmental dose in complex locations, such as high altitude or polar neutron radiation research stations where liquid water and ice are both present: for neutron energies between hundreds of meV and tens of keV, the total scattering cross section per molecule in the two phases can be considered the same, with the macroscopic cross section only depending upon the density changes of water near the melting point.
Collapse
Affiliation(s)
- Giovanni Romanelli
- Dipartimento di Fisica and NAST Centre, Università degli Studi di Roma "Tor Vergata," via della Ricerca Scientifica 1, 00133 Rome, Italy
| | - Carla Andreani
- Dipartimento di Fisica and NAST Centre, Università degli Studi di Roma "Tor Vergata," via della Ricerca Scientifica 1, 00133 Rome, Italy
| | - Alessio Bocedi
- Dipartimento di Scienze e Tecnologie Chimiche and NAST Centre, Università degli Studi di Roma "Tor Vergata," via della Ricerca Scientifica 1, 00133 Rome, Italy
| | - Roberto Senesi
- Dipartimento di Fisica and NAST Centre, Università degli Studi di Roma "Tor Vergata," via della Ricerca Scientifica 1, 00133 Rome, Italy
| |
Collapse
|
2
|
Lambros E, Fetherolf JH, Hammes-Schiffer S, Li X. A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters. J Phys Chem Lett 2024; 15:4070-4075. [PMID: 38587257 DOI: 10.1021/acs.jpclett.4c00439] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. By performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jonathan H Fetherolf
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
3
|
Montero de Hijes P, Dellago C, Jinnouchi R, Schmiedmayer B, Kresse G. Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks. J Chem Phys 2024; 160:114107. [PMID: 38506284 DOI: 10.1063/5.0197105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 03/03/2024] [Indexed: 03/21/2024] Open
Abstract
In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests.
Collapse
Affiliation(s)
- Pablo Montero de Hijes
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- University of Vienna, Faculty of Earth Sciences, Geography and Astronomy, Josef-Holaubuek-Platz 2, 1090 Vienna, Austria
| | - Christoph Dellago
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
| | - Ryosuke Jinnouchi
- Toyota Central R&D Labs., Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | | | - Georg Kresse
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- VASP Software GmbH, Berggasse 21, A-1090 Vienna, Austria
| |
Collapse
|
4
|
Reinhardt A, Chew PY, Cheng B. A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals. J Chem Phys 2023; 159:184110. [PMID: 37962445 DOI: 10.1063/5.0173341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 10/20/2023] [Indexed: 11/15/2023] Open
Abstract
Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances.
Collapse
Affiliation(s)
- Aleks Reinhardt
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Pin Yu Chew
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Bingqing Cheng
- Institute of Science and Technology Austria, Am Campus 1, 3400 Klosterneuburg, Austria
| |
Collapse
|
5
|
Feldman YMY, Hirshberg B. Quadratic scaling bosonic path integral molecular dynamics. J Chem Phys 2023; 159:154107. [PMID: 37855315 DOI: 10.1063/5.0173749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Accepted: 09/15/2023] [Indexed: 10/20/2023] Open
Abstract
Bosonic exchange symmetry leads to fascinating quantum phenomena, from exciton condensation in quantum materials to the superfluidity of liquid 4He. Unfortunately, path integral molecular dynamics (PIMD) simulations of bosons are computationally prohibitive beyond ∼100 particles, due to a cubic scaling with the system size. We present an algorithm that reduces the complexity from cubic to quadratic, allowing the first simulations of thousands of bosons using PIMD. Our method is orders of magnitude faster, with a speedup that scales linearly with the number of particles and the number of imaginary time slices (beads). Simulations that would have otherwise taken decades can now be done in days. In practice, the new algorithm eliminates most of the added computational cost of including bosonic exchange effects, making them almost as accessible as PIMD simulations of distinguishable particles.
Collapse
Affiliation(s)
- Yotam M Y Feldman
- School of Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
- The Ratner Center for Single Molecule Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Barak Hirshberg
- School of Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
- The Ratner Center for Single Molecule Science, Tel Aviv University, Tel Aviv 6997801, Israel
- The Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
6
|
Daru J, Forbert H, Behler J, Marx D. Coupled Cluster Molecular Dynamics of Condensed Phase Systems Enabled by Machine Learning Potentials: Liquid Water Benchmark. PHYSICAL REVIEW LETTERS 2022; 129:226001. [PMID: 36493459 DOI: 10.1103/physrevlett.129.226001] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/17/2022] [Revised: 09/05/2022] [Accepted: 10/05/2022] [Indexed: 06/17/2023]
Abstract
Coupled cluster theory is a general and systematic electronic structure method, but in particular the highly accurate "gold standard" coupled cluster singles, doubles and perturbative triples, CCSD(T), can only be applied to small systems. To overcome this limitation, we introduce a framework to transfer CCSD(T) accuracy of finite molecular clusters to extended condensed phase systems using a high-dimensional neural network potential. This approach, which is automated, allows one to perform high-quality coupled cluster molecular dynamics, CCMD, as we demonstrate for liquid water including nuclear quantum effects. The machine learning strategy is very efficient, generic, can be systematically improved, and is applicable to a variety of complex systems.
Collapse
Affiliation(s)
- János Daru
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
7
|
Reinhardt A, Bethkenhagen M, Coppari F, Millot M, Hamel S, Cheng B. Thermodynamics of high-pressure ice phases explored with atomistic simulations. Nat Commun 2022; 13:4707. [PMID: 35948550 PMCID: PMC9365810 DOI: 10.1038/s41467-022-32374-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 07/27/2022] [Indexed: 11/25/2022] Open
Abstract
Most experimentally known high-pressure ice phases have a body-centred cubic (bcc) oxygen lattice. Our large-scale molecular-dynamics simulations with a machine-learning potential indicate that, amongst these bcc ice phases, ices VII, VII′ and X are the same thermodynamic phase under different conditions, whereas superionic ice VII″ has a first-order phase boundary with ice VII′. Moreover, at about 300 GPa, the transformation between ice X and the Pbcm phase has a sharp structural change but no apparent activation barrier, whilst at higher pressures the barrier gradually increases. Our study thus clarifies the phase behaviour of the high-pressure ices and reveals peculiar solid–solid transition mechanisms not known in other systems. Many experimentally known high-pressure ice phase are structurally very similar. Here authors elucidate the phase behaviour of the high-pressure insulating ices and reveal solid-solid transition mechanisms not known in other systems.
Collapse
Affiliation(s)
- Aleks Reinhardt
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK
| | - Mandy Bethkenhagen
- École Normale Supérieure de Lyon, Université Lyon 1, Laboratoire de Géologie de Lyon, CNRS UMR 5276, 69364, Lyon Cedex 07, France
| | - Federica Coppari
- Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
| | - Marius Millot
- Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
| | - Sebastien Hamel
- Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
| | - Bingqing Cheng
- Institute of Science and Technology Austria, Am Campus 1, 3400, Klosterneuburg, Austria.
| |
Collapse
|
8
|
Eckhoff M, Behler J. Insights into lithium manganese oxide-water interfaces using machine learning potentials. J Chem Phys 2021; 155:244703. [PMID: 34972388 DOI: 10.1063/5.0073449] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Unraveling the atomistic and the electronic structure of solid-liquid interfaces is the key to the design of new materials for many important applications, from heterogeneous catalysis to battery technology. Density functional theory (DFT) calculations can, in principle, provide a reliable description of such interfaces, but the high computational costs severely restrict the accessible time and length scales. Here, we report machine learning-driven simulations of various interfaces between water and lithium manganese oxide (LixMn2O4), an important electrode material in lithium ion batteries and a catalyst for the oxygen evolution reaction. We employ a high-dimensional neural network potential to compute the energies and forces several orders of magnitude faster than DFT without loss in accuracy. In addition, a high-dimensional neural network for spin prediction is utilized to analyze the electronic structure of the manganese ions. Combining these methods, a series of interfaces is investigated by large-scale molecular dynamics. The simulations allow us to gain insights into a variety of properties, such as the dissociation of water molecules, proton transfer processes, and hydrogen bonds, as well as the geometric and electronic structure of the solid surfaces, including the manganese oxidation state distribution, Jahn-Teller distortions, and electron hopping.
Collapse
Affiliation(s)
- Marco Eckhoff
- Institut für Physikalische Chemie, Theoretische Chemie, Universität Göttingen, Tammannstraße 6, 37077 Göttingen, Germany
| | - Jörg Behler
- Institut für Physikalische Chemie, Theoretische Chemie, Universität Göttingen, Tammannstraße 6, 37077 Göttingen, Germany
| |
Collapse
|
9
|
Campos-Villalobos G, Boattini E, Filion L, Dijkstra M. Machine learning many-body potentials for colloidal systems. J Chem Phys 2021; 155:174902. [PMID: 34742191 DOI: 10.1063/5.0063377] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Simulations of colloidal suspensions consisting of mesoscopic particles and smaller species such as ions or depletants are computationally challenging as different length and time scales are involved. Here, we introduce a machine learning (ML) approach in which the degrees of freedom of the microscopic species are integrated out and the mesoscopic particles interact with effective many-body potentials, which we fit as a function of all colloid coordinates with a set of symmetry functions. We apply this approach to a colloid-polymer mixture. Remarkably, the ML potentials can be assumed to be effectively state-independent and can be used in direct-coexistence simulations. We show that our ML method reduces the computational cost by several orders of magnitude compared to a numerical evaluation and accurately describes the phase behavior and structure, even for state points where the effective potential is largely determined by many-body contributions.
Collapse
Affiliation(s)
- Gerardo Campos-Villalobos
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 1, 3584 CC Utrecht, The Netherlands
| | - Emanuele Boattini
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 1, 3584 CC Utrecht, The Netherlands
| | - Laura Filion
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 1, 3584 CC Utrecht, The Netherlands
| | - Marjolein Dijkstra
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 1, 3584 CC Utrecht, The Netherlands
| |
Collapse
|
10
|
Keith JA, Vassilev-Galindo V, Cheng B, Chmiela S, Gastegger M, Müller KR, Tkatchenko A. Combining Machine Learning and Computational Chemistry for Predictive Insights Into Chemical Systems. Chem Rev 2021; 121:9816-9872. [PMID: 34232033 PMCID: PMC8391798 DOI: 10.1021/acs.chemrev.1c00107] [Citation(s) in RCA: 190] [Impact Index Per Article: 63.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Indexed: 12/23/2022]
Abstract
Machine learning models are poised to make a transformative impact on chemical sciences by dramatically accelerating computational algorithms and amplifying insights available from computational chemistry methods. However, achieving this requires a confluence and coaction of expertise in computer science and physical sciences. This Review is written for new and experienced researchers working at the intersection of both fields. We first provide concise tutorials of computational chemistry and machine learning methods, showing how insights involving both can be achieved. We follow with a critical review of noteworthy applications that demonstrate how computational chemistry and machine learning can be used together to provide insightful (and useful) predictions in molecular and materials modeling, retrosyntheses, catalysis, and drug design.
Collapse
Affiliation(s)
- John A. Keith
- Department
of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Valentin Vassilev-Galindo
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Bingqing Cheng
- Accelerate
Programme for Scientific Discovery, Department
of Computer Science and Technology, 15 J. J. Thomson Avenue, Cambridge CB3 0FD, United Kingdom
| | - Stefan Chmiela
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Michael Gastegger
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587, Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea
- Max-Planck-Institut für Informatik, 66123 Saarbrücken, Germany
- Google Research, Brain Team, 10117 Berlin, Germany
| | - Alexandre Tkatchenko
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
11
|
Zhang L, Wang H, Car R, E W. Phase Diagram of a Deep Potential Water Model. PHYSICAL REVIEW LETTERS 2021; 126:236001. [PMID: 34170175 DOI: 10.1103/physrevlett.126.236001] [Citation(s) in RCA: 85] [Impact Index Per Article: 28.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Accepted: 04/28/2021] [Indexed: 06/13/2023]
Abstract
Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. The fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII^{''} and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.
Collapse
Affiliation(s)
- Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Fenghao East Road 2, Beijing 100094, People's Republic of China
| | - Roberto Car
- Department of Chemistry, Department of Physics, Program in Applied and Computational Mathematics, Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, USA
| | - Weinan E
- Department of Mathematics and Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA and Beijing Institute of Big Data Research, Beijing 100871, People's Republic of China
| |
Collapse
|
12
|
Rossi M. Progress and challenges in ab initio simulations of quantum nuclei in weakly bonded systems. J Chem Phys 2021; 154:170902. [PMID: 34241065 DOI: 10.1063/5.0042572] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Atomistic simulations based on the first-principles of quantum mechanics are reaching unprecedented length scales. This progress is due to the growth in computational power allied with the development of new methodologies that allow the treatment of electrons and nuclei as quantum particles. In the realm of materials science, where the quest for desirable emergent properties relies increasingly on soft weakly bonded materials, such methods have become indispensable. In this Perspective, an overview of simulation methods that are applicable for large system sizes and that can capture the quantum nature of electrons and nuclei in the adiabatic approximation is given. In addition, the remaining challenges are discussed, especially regarding the inclusion of nuclear quantum effects (NQEs) beyond a harmonic or perturbative treatment, the impact of NQEs on electronic properties of weakly bonded systems, and how different first-principles potential energy surfaces can change the impact of NQEs on the atomic structure and dynamics of weakly bonded systems.
Collapse
Affiliation(s)
- Mariana Rossi
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
13
|
Drużbicki K, Gaboardi M, Fernandez-Alonso F. Dynamics & Spectroscopy with Neutrons-Recent Developments & Emerging Opportunities. Polymers (Basel) 2021; 13:1440. [PMID: 33947108 PMCID: PMC8125526 DOI: 10.3390/polym13091440] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Accepted: 04/27/2021] [Indexed: 12/19/2022] Open
Abstract
This work provides an up-to-date overview of recent developments in neutron spectroscopic techniques and associated computational tools to interrogate the structural properties and dynamical behavior of complex and disordered materials, with a focus on those of a soft and polymeric nature. These have and continue to pave the way for new scientific opportunities simply thought unthinkable not so long ago, and have particularly benefited from advances in high-resolution, broadband techniques spanning energy transfers from the meV to the eV. Topical areas include the identification and robust assignment of low-energy modes underpinning functionality in soft solids and supramolecular frameworks, or the quantification in the laboratory of hitherto unexplored nuclear quantum effects dictating thermodynamic properties. In addition to novel classes of materials, we also discuss recent discoveries around water and its phase diagram, which continue to surprise us. All throughout, emphasis is placed on linking these ongoing and exciting experimental and computational developments to specific scientific questions in the context of the discovery of new materials for sustainable technologies.
Collapse
Affiliation(s)
- Kacper Drużbicki
- Materials Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastian, Spain;
- Polish Academy of Sciences, Center of Molecular and Macromolecular Studies, Sienkiewicza 112, 90-363 Lodz, Poland
| | - Mattia Gaboardi
- Elettra—Sincrotrone Trieste S.C.p.A., S.S. 14 km 163.5 in Area Science Park, 34149 Trieste, Italy;
| | - Felix Fernandez-Alonso
- Materials Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastian, Spain;
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastian, Spain
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, 48009 Bilbao, Spain
| |
Collapse
|
14
|
Affiliation(s)
- Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstraße 6, 37077 Göttingen, Germany
| |
Collapse
|
15
|
Reinhardt A, Cheng B. Quantum-mechanical exploration of the phase diagram of water. Nat Commun 2021; 12:588. [PMID: 33500405 PMCID: PMC7838264 DOI: 10.1038/s41467-020-20821-w] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Accepted: 12/21/2020] [Indexed: 11/10/2022] Open
Abstract
The set of known stable phases of water may not be complete, and some of the phase boundaries between them are fuzzy. Starting from liquid water and a comprehensive set of 50 ice structures, we compute the phase diagram at three hybrid density-functional-theory levels of approximation, accounting for thermal and nuclear fluctuations as well as proton disorder. Such calculations are only made tractable because we combine machine-learning methods and advanced free-energy techniques. The computed phase diagram is in qualitative agreement with experiment, particularly at pressures ≲ 8000 bar, and the discrepancy in chemical potential is comparable with the subtle uncertainties introduced by proton disorder and the spread between the three hybrid functionals. None of the hypothetical ice phases considered is thermodynamically stable in our calculations, suggesting the completeness of the experimental water phase diagram in the region considered. Our work demonstrates the feasibility of predicting the phase diagram of a polymorphic system from first principles and provides a thermodynamic way of testing the limits of quantum-mechanical calculations.
Collapse
Affiliation(s)
- Aleks Reinhardt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.
| | - Bingqing Cheng
- Accelerate Programme for Scientific Discovery, Department of Computer Science and Technology, 15 J.J. Thomson Avenue, Cambridge, CB3 0FD, UK. .,Cavendish Laboratory, University of Cambridge, J.J. Thomson Avenue, Cambridge, CB3 0HE, UK.
| |
Collapse
|
16
|
Yue S, Muniz MC, Calegari Andrade MF, Zhang L, Car R, Panagiotopoulos AZ. When do short-range atomistic machine-learning models fall short? J Chem Phys 2021; 154:034111. [DOI: 10.1063/5.0031215] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Shuwen Yue
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | | |
Collapse
|
17
|
Engel EA. Identification of synthesisable crystalline phases of water – a prototype for the challenges of computational materials design. CrystEngComm 2021. [DOI: 10.1039/d0ce01260b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We discuss the identification of experimentally realisable crystalline phases of water to outline and contextualise some of the diverse building blocks of a computational materials design process.
Collapse
Affiliation(s)
- Edgar A. Engel
- TCM Group
- Cavendish Laboratory
- University of Cambridge
- Cambridge CB3 0HE
- UK
| |
Collapse
|
18
|
Ulpiani P, Romanelli G, Onorati D, Krzystyniak M, Andreani C, Senesi R. The effective isotropy of the hydrogen local potential in biphenyl and other hydrocarbons. J Chem Phys 2020; 153:234306. [PMID: 33353342 DOI: 10.1063/5.0029578] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an experimental investigation of the hydrogen nuclear momentum distribution in biphenyl using deep inelastic neutron scattering. Our experimental results suggest that the local potential affecting hydrogen is both harmonic and isotropic within experimental uncertainties. This feature is interpreted as a consequence of the central limit theorem, whereby the three-dimensional momentum distribution is expected to become a purely Gaussian function as the number of independent vibrational modes in a system increases. We also performed ab initio phonon calculations on biphenyl and other saturated hydrocarbons, from methane to decane. From the results of the simulations, one can observe that the nuclear momentum distribution becomes more isotropic as the number of atoms and normal modes in the molecule increases. Moreover, the predicted theoretical anisotropy in biphenyl is clearly larger than in the experiment. The reason is that the total number of normal modes necessary to reproduce the experimental results is much larger than the number of normal modes encompassed by a single unit cell due to the presence of structural disorder and intermolecular interactions in the real crystal, as well as coupling of different normal modes. Finally, experimental data were collected, over a subset of detectors on the VESUVIO spectrometer at ISIS, with a novel setup to increase the count rate and signal-to-background ratio. We envision that such an optimized experimental setup can provide faster measurements and more stringent constraints for phonon calculations.
Collapse
Affiliation(s)
- Pierfrancesco Ulpiani
- Università degli Studi di Roma "Tor Vergata," Dipartimento di Scienze e Tecnologie Chimiche, Via della Ricerca Scientifica 1, Rome 00133, Italy
| | - Giovanni Romanelli
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11OQX, United Kingdom
| | - Dalila Onorati
- Università degli Studi di Roma "Tor Vergata," Dipartimento di Fisica and NAST Center, Via della Ricerca Scientifica 1, Rome 00133, Italy
| | - Matthew Krzystyniak
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11OQX, United Kingdom
| | - Carla Andreani
- Università degli Studi di Roma "Tor Vergata," Dipartimento di Fisica and NAST Center, Via della Ricerca Scientifica 1, Rome 00133, Italy
| | - Roberto Senesi
- Università degli Studi di Roma "Tor Vergata," Dipartimento di Fisica and NAST Center, Via della Ricerca Scientifica 1, Rome 00133, Italy
| |
Collapse
|
19
|
Boattini E, Bezem N, Punnathanam SN, Smallenburg F, Filion L. Modeling of many-body interactions between elastic spheres through symmetry functions. J Chem Phys 2020; 153:064902. [DOI: 10.1063/5.0015606] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Emanuele Boattini
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Utrecht, The Netherlands
| | - Nina Bezem
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Utrecht, The Netherlands
| | - Sudeep N. Punnathanam
- Department of Chemical Engineering, Indian Institute of Science, Bangalore 560012, Karnataka, India
| | - Frank Smallenburg
- Université Paris-Saclay, CNRS, Laboratoire de Physique des Solides, 91405 Orsay, France
| | - Laura Filion
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Utrecht, The Netherlands
| |
Collapse
|
20
|
Pattnaik P, Raghunathan S, Kalluri T, Bhimalapuram P, Jawahar CV, Priyakumar UD. Machine Learning for Accurate Force Calculations in Molecular Dynamics Simulations. J Phys Chem A 2020; 124:6954-6967. [DOI: 10.1021/acs.jpca.0c03926] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Punyaslok Pattnaik
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Shampa Raghunathan
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Tarun Kalluri
- Center for Visual Information Technology, KCIS, International Institute of Information Technology, Hyderabad 500 032, India
| | - Prabhakar Bhimalapuram
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - C. V. Jawahar
- Center for Visual Information Technology, KCIS, International Institute of Information Technology, Hyderabad 500 032, India
| | - U. Deva Priyakumar
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| |
Collapse
|
21
|
Cheng Z, Zhao D, Ma J, Li W, Li S. An On-the-Fly Approach to Construct Generalized Energy-Based Fragmentation Machine Learning Force Fields of Complex Systems. J Phys Chem A 2020; 124:5007-5014. [DOI: 10.1021/acs.jpca.0c04526] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Zheng Cheng
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Dongbo Zhao
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Jing Ma
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Wei Li
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Shuhua Li
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| |
Collapse
|
22
|
Wang L, Ceriotti M, Markland TE. Quantum kinetic energy and isotope fractionation in aqueous ionic solutions. Phys Chem Chem Phys 2020; 22:10490-10499. [PMID: 31942581 DOI: 10.1039/c9cp06483d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
At room temperature, the quantum contribution to the kinetic energy of a water molecule exceeds the classical contribution by an order of magnitude. The quantum kinetic energy (QKE) of a water molecule is modulated by its local chemical environment and leads to uneven partitioning of isotopes between different phases in thermal equilibrium, which would not occur if the nuclei behaved classically. In this work, we use ab initio path integral simulations to show that QKEs of the water molecules and the equilibrium isotope fractionation ratios of the oxygen and hydrogen isotopes are sensitive probes of the hydrogen bonding structures in aqueous ionic solutions. In particular, we demonstrate how the QKE of water molecules in path integral simulations can be decomposed into translational, rotational and vibrational degrees of freedom, and use them to determine the impact of solvation on different molecular motions. By analyzing the QKEs and isotope fractionation ratios, we show how the addition of the Na+, Cl- and HPO42- ions perturbs the competition between quantum effects in liquid water and impacts their local solvation structures.
Collapse
Affiliation(s)
- Lu Wang
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA.
| | | | | |
Collapse
|
23
|
Karandashev K, Vaníček J. Accelerating equilibrium isotope effect calculations. II. Stochastic implementation of direct estimators. J Chem Phys 2019; 151:134116. [PMID: 31594323 DOI: 10.1063/1.5124995] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Path integral calculations of equilibrium isotope effects and isotopic fractionation are expensive due to the presence of path integral discretization errors, statistical errors, and thermodynamic integration errors. Whereas the discretization errors can be reduced by high-order factorization of the path integral and statistical errors by using centroid virial estimators, two recent papers proposed alternative ways to completely remove the thermodynamic integration errors: Cheng and Ceriotti [J. Chem. Phys. 141, 244112 (2015)] employed a variant of free-energy perturbation called "direct estimators," while Karandashev and Vaníček [J. Chem. Phys. 143, 194104 (2017)] combined the thermodynamic integration with a stochastic change of mass and piecewise-linear umbrella biasing potential. Here, we combine the former approach with the stochastic change in mass in order to decrease its statistical errors when applied to larger isotope effects and perform a thorough comparison of different methods by computing isotope effects first on a harmonic model and then on methane and methanium, where we evaluate all isotope effects of the form CH4-xDx/CH4 and CH5-xDx +/CH5 +, respectively. We discuss the reasons for a surprising behavior of the original method of direct estimators, which performed well for a much larger range of isotope effects than what had been expected previously, as well as some implications of our work for the more general problem of free energy difference calculations.
Collapse
Affiliation(s)
- Konstantin Karandashev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
24
|
Ulpiani P, Romanelli G, Onorati D, Parmentier A, Festa G, Schooneveld E, Cazzaniga C, Arcidiacono L, Andreani C, Senesi R. Optimization of detection strategies for epithermal neutron spectroscopy using photon-sensitive detectors. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2019; 90:073901. [PMID: 31370488 DOI: 10.1063/1.5091084] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2019] [Accepted: 06/15/2019] [Indexed: 06/10/2023]
Abstract
In this work, we discuss an improved detection procedure for the photon-sensitive yttrium-aluminum-perovskite detectors installed on the VESUVIO spectrometer at the ISIS pulsed neutron and muon source. By decreasing the low-level energy threshold of detected photons, we observe an increased count rate up to a factor ∼3, and a decrease of relative error bars and noise of ∼40% and 35%, respectively, for deep inelastic neutron scattering measurements. In addition, we demonstrate how the reported optimization may increase the accuracy in the line shape analysis of neutron Compton profiles, as well as in the application of the mean-force approach to detect the anisotropy and anharmonicity in the single-particle local potential. We envisage that such an upgrade of the detection procedure would have a substantial impact on the VESUVIO scientific programme based on deep inelastic neutron scattering investigations.
Collapse
Affiliation(s)
- Pierfrancesco Ulpiani
- Università degli studi di Roma "Tor Vergata", Dipartimento di Scienze e Tecnologie Chimiche, Via della Ricerca Scientifica 1, Rome, 00133, Italy
| | - Giovanni Romanelli
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 OQX, United Kingdom
| | - Dalila Onorati
- Università degli studi di Roma "Tor Vergata", Centro NAST, Via della Ricerca Scientifica 1, Rome, 00133, Italy
| | - Alexandra Parmentier
- INFN Sezione di Roma Tor Vergata, Via della Ricerca Scientifica 1, I-00133 Rome, Italy
| | - Giulia Festa
- Centro Fermi - Museo Storico della Fisica e Centro Studi e Ricerche "Enrico Fermi", Piazza del Viminale 1, Rome, 00184, Italy
| | - Erik Schooneveld
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 OQX, United Kingdom
| | - Carlo Cazzaniga
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 OQX, United Kingdom
| | - Laura Arcidiacono
- Università degli studi di Roma "Tor Vergata", Centro NAST, Via della Ricerca Scientifica 1, Rome, 00133, Italy
| | - Carla Andreani
- Università degli studi di Roma "Tor Vergata", Centro NAST, Via della Ricerca Scientifica 1, Rome, 00133, Italy
| | - Roberto Senesi
- Università degli studi di Roma "Tor Vergata", Centro NAST, Via della Ricerca Scientifica 1, Rome, 00133, Italy
| |
Collapse
|
25
|
Eckhoff M, Behler J. From Molecular Fragments to the Bulk: Development of a Neural Network Potential for MOF-5. J Chem Theory Comput 2019; 15:3793-3809. [PMID: 31091097 DOI: 10.1021/acs.jctc.8b01288] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The development of first-principles-quality reactive atomistic potentials for organic-inorganic hybrid materials is still a substantial challenge because of the very different physics of the atomic interactions-from covalent via ionic bonding to dispersion-that have to be described in an accurate and balanced way. In this work we used a prototypical metal-organic framework, MOF-5, as a benchmark case to investigate the applicability of high-dimensional neural network potentials (HDNNPs) to this class of materials. In HDNNPs, which belong to the class of machine learning potentials, the energy is constructed as a sum of environment-dependent atomic energy contributions. We demonstrate that by the use of this approach it is possible to obtain a high-quality potential for the periodic MOF-5 crystal using density functional theory (DFT) reference calculations of small molecular fragments only. The resulting HDNNP, which has a root-mean-square error (RMSE) of 1.6 meV/atom for the energies of molecular fragments not included in the training set, is able to provide the equilibrium lattice constant of the bulk MOF-5 structure with an error of about 0.1% relative to DFT, and also, the negative thermal expansion behavior is accurately predicted. The total energy RMSE of periodic structures that are completely absent in the training set is about 6.5 meV/atom, with errors on the order of 2 meV/atom for energy differences. We show that in contrast to energy differences, achieving a high accuracy for total energies requires careful variation of the stoichiometries of the training structures to avoid energy offsets, as atomic energies are not physical observables. The forces, which have RMSEs of about 94 meV/ a0 for the molecular fragments and 130 meV/ a0 for bulk structures not included in the training set, are insensitive to such offsets. Therefore, forces, which are the relevant properties for molecular dynamics simulations, provide a realistic estimate of the accuracy of atomistic potentials.
Collapse
Affiliation(s)
- Marco Eckhoff
- Universität Göttingen , Institut für Physikalische Chemie, Theoretische Chemie , Tammannstraße 6 , D-37077 Göttingen , Germany
| | - Jörg Behler
- Universität Göttingen , Institut für Physikalische Chemie, Theoretische Chemie , Tammannstraße 6 , D-37077 Göttingen , Germany
| |
Collapse
|
26
|
Affiliation(s)
- Wei Fang
- School of Physics and Collaborative Innovation Centre of Quantum Matter, Peking University, Beijing, People's Republic of China
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Physics and Astronomy, University College London, London, UK
- Laboratory of Physical Chemistry, ETH Zurich, Zurich, Switzerland
| | - Ji Chen
- Department of Electronic Structure Theory, Max Plank Institute for Solid State Research, Stuttgart, Germany
| | - Yexin Feng
- School of Physics and Electronics, Hunan University, Changsha, People's Republic of China
| | - Xin-Zheng Li
- School of Physics and Collaborative Innovation Centre of Quantum Matter, Peking University, Beijing, People's Republic of China
- State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Peking University, Beijing, People's Republic of China
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Physics and Astronomy, University College London, London, UK
| |
Collapse
|
27
|
Singraber A, Behler J, Dellago C. Library-Based LAMMPS Implementation of High-Dimensional Neural Network Potentials. J Chem Theory Comput 2019; 15:1827-1840. [DOI: 10.1021/acs.jctc.8b00770] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Andreas Singraber
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstraße 6, 37077 Göttingen, Germany
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| |
Collapse
|
28
|
Cheng B, Engel EA, Behler J, Dellago C, Ceriotti M. Ab initio thermodynamics of liquid and solid water. Proc Natl Acad Sci U S A 2019; 116:1110-1115. [PMID: 30610171 PMCID: PMC6347673 DOI: 10.1073/pnas.1815117116] [Citation(s) in RCA: 152] [Impact Index Per Article: 30.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Thermodynamic properties of liquid water as well as hexagonal (Ih) and cubic (Ic) ice are predicted based on density functional theory at the hybrid-functional level, rigorously taking into account quantum nuclear motion, anharmonic fluctuations, and proton disorder. This is made possible by combining advanced free-energy methods and state-of-the-art machine-learning techniques. The ab initio description leads to structural properties in excellent agreement with experiments and reliable estimates of the melting points of light and heavy water. We observe that nuclear-quantum effects contribute a crucial [Formula: see text] to the stability of ice Ih, making it more stable than ice Ic. Our computational approach is general and transferable, providing a comprehensive framework for quantitative predictions of ab initio thermodynamic properties using machine-learning potentials as an intermediate step.
Collapse
Affiliation(s)
- Bingqing Cheng
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland;
| | - Edgar A Engel
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Universität Göttingen, 37073 Göttingen, Germany
| | | | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
29
|
Hellström M, Ceriotti M, Behler J. Nuclear Quantum Effects in Sodium Hydroxide Solutions from Neural Network Molecular Dynamics Simulations. J Phys Chem B 2018; 122:10158-10171. [DOI: 10.1021/acs.jpcb.8b06433] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Matti Hellström
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstr. 6, 37077 Göttingen, Germany
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstr. 6, 37077 Göttingen, Germany
| |
Collapse
|
30
|
Nguyen TT, Székely E, Imbalzano G, Behler J, Csányi G, Ceriotti M, Götz AW, Paesani F. Comparison of permutationally invariant polynomials, neural networks, and Gaussian approximation potentials in representing water interactions through many-body expansions. J Chem Phys 2018; 148:241725. [DOI: 10.1063/1.5024577] [Citation(s) in RCA: 118] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Thuong T. Nguyen
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - Eszter Székely
- Engineering Department, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| | - Giulio Imbalzano
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstr. 6, 37077 Göttingen, Germany
| | - Gábor Csányi
- Engineering Department, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
31
|
Imbalzano G, Anelli A, Giofré D, Klees S, Behler J, Ceriotti M. Automatic selection of atomic fingerprints and reference configurations for machine-learning potentials. J Chem Phys 2018; 148:241730. [DOI: 10.1063/1.5024611] [Citation(s) in RCA: 163] [Impact Index Per Article: 27.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Affiliation(s)
- Giulio Imbalzano
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Andrea Anelli
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Daniele Giofré
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Sinja Klees
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstr. 6, 37077 Göttingen, Germany
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
32
|
Kapil V, Cuzzocrea A, Ceriotti M. Anisotropy of the Proton Momentum Distribution in Water. J Phys Chem B 2018; 122:6048-6054. [DOI: 10.1021/acs.jpcb.8b03896] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Alice Cuzzocrea
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| |
Collapse
|
33
|
Rossi M, Kapil V, Ceriotti M. Fine tuning classical and quantum molecular dynamics using a generalized Langevin equation. J Chem Phys 2018; 148:102301. [DOI: 10.1063/1.4990536] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mariana Rossi
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
34
|
Schran C, Uhl F, Behler J, Marx D. High-dimensional neural network potentials for solvation: The case of protonated water clusters in helium. J Chem Phys 2018; 148:102310. [DOI: 10.1063/1.4996819] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Felix Uhl
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
- Theoretische Chemie, Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
35
|
|
36
|
Behler J. First Principles Neural Network Potentials for Reactive Simulations of Large Molecular and Condensed Systems. Angew Chem Int Ed Engl 2017; 56:12828-12840. [PMID: 28520235 DOI: 10.1002/anie.201703114] [Citation(s) in RCA: 317] [Impact Index Per Article: 45.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2017] [Indexed: 11/06/2022]
Abstract
Modern simulation techniques have reached a level of maturity which allows a wide range of problems in chemistry and materials science to be addressed. Unfortunately, the application of first principles methods with predictive power is still limited to rather small systems, and despite the rapid evolution of computer hardware no fundamental change in this situation can be expected. Consequently, the development of more efficient but equally reliable atomistic potentials to reach an atomic level understanding of complex systems has received considerable attention in recent years. A promising new development has been the introduction of machine learning (ML) methods to describe the atomic interactions. Once trained with electronic structure data, ML potentials can accelerate computer simulations by several orders of magnitude, while preserving quantum mechanical accuracy. This Review considers the methodology of an important class of ML potentials that employs artificial neural networks.
Collapse
Affiliation(s)
- Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, Tammannstrasse 6, 37077, Göttingen, Germany
| |
Collapse
|
37
|
Behler J. Hochdimensionale neuronale Netze für Potentialhyperflächen großer molekularer und kondensierter Systeme. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201703114] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Jörg Behler
- Universität Göttingen; Institut für Physikalische Chemie, Theoretische Chemie; Tammannstraße 6 37077 Göttingen Deutschland
| |
Collapse
|
38
|
Affiliation(s)
- Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum, Germany
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
39
|
Andreani C, Romanelli G, Senesi R. Direct Measurements of Quantum Kinetic Energy Tensor in Stable and Metastable Water near the Triple Point: An Experimental Benchmark. J Phys Chem Lett 2016; 7:2216-2220. [PMID: 27214268 DOI: 10.1021/acs.jpclett.6b00926] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
This study presents the first direct and quantitative measurement of the nuclear momentum distribution anisotropy and the quantum kinetic energy tensor in stable and metastable (supercooled) water near its triple point, using deep inelastic neutron scattering (DINS). From the experimental spectra, accurate line shapes of the hydrogen momentum distributions are derived using an anisotropic Gaussian and a model-independent framework. The experimental results, benchmarked with those obtained for the solid phase, provide the state of the art directional values of the hydrogen mean kinetic energy in metastable water. The determinations of the direction kinetic energies in the supercooled phase, provide accurate and quantitative measurements of these dynamical observables in metastable and stable phases, that is, key insight in the physical mechanisms of the hydrogen quantum state in both disordered and polycrystalline systems. The remarkable findings of this study establish novel insight into further expand the capacity and accuracy of DINS investigations of the nuclear quantum effects in water and represent reference experimental values for theoretical investigations.
Collapse
Affiliation(s)
- Carla Andreani
- Università degli Studi di Roma "Tor Vergata" , Dipartimento di Fisica e Centro NAST, Via della Ricerca Scientifica 1, 00133 Roma, Italy
- Consiglio Nazionale delle Ricerche, CNR-IPCF, Sezione di Messina 98122, Italy
| | - Giovanni Romanelli
- ISIS Neutron Source, Science Technology Facility Council, Chilton, Oxfordshire OX11 0QX, United Kingdom
| | - Roberto Senesi
- Università degli Studi di Roma "Tor Vergata" , Dipartimento di Fisica e Centro NAST, Via della Ricerca Scientifica 1, 00133 Roma, Italy
- Consiglio Nazionale delle Ricerche, CNR-IPCF, Sezione di Messina 98122, Italy
| |
Collapse
|
40
|
Natarajan SK, Behler J. Neural network molecular dynamics simulations of solid–liquid interfaces: water at low-index copper surfaces. Phys Chem Chem Phys 2016; 18:28704-28725. [DOI: 10.1039/c6cp05711j] [Citation(s) in RCA: 114] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Molecular dynamics simulation of the water–copper interface have been carried out using high-dimensional neural network potential based on density functional theory.
Collapse
Affiliation(s)
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie
- Ruhr-Universität Bochum
- D-44780 Bochum
- Germany
| |
Collapse
|