1
|
Creazzo F, Luber S. Water-air interface revisited by means of path-integral ab initio molecular dynamics. Phys Chem Chem Phys 2024; 26:21290-21302. [PMID: 39078670 PMCID: PMC11305098 DOI: 10.1039/d4cp02500h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2024] [Accepted: 07/21/2024] [Indexed: 07/31/2024]
Abstract
Although nuclear quantum effects (NQEs) have been considered on bulk liquid water, the impact of these latter on the air-water interface has not yet been reported. Herein, by performing and comparing ab initio molecular dynamics (AIMD) and path integral AIMD (PI-AIMD) simulations, we reveal the impact of NQEs on structural, dynamical and electronic properties as well as IR spectra of the air-water interface at room temperature. NQEs, being able to describe a more accurate proton delocalization in H-bonded system than AIMD, reveal a different structural arrangement and dynamical behaviour of both bulk and interfacial water molecules in comparison to AIMD results. A more de-structured and de-bound water arrangement and coordination are identified when the quantum nature of nuclei are considered for both bulk and interfacial water molecules. Structural properties, such as inter-/intra-molecular bond lengths, coordination numbers and H-bonding angles of bulk and interfacial water molecules here calculated, are affected by NQEs mitigating the overstructured description given by AIMD. Further evidences of an AIMD overstructured description of bulk water are in the computed IR spectra, where an increased absorption peak intensity and an increased strength of the hydrogen-bond network are alleviated by NQEs. In addition, NQEs show a valuable impact on the electronic structure of the air-water interface, reducing the total valence bandwidth and the electronic energy band-gap when passing from bulk to interfacial water. This work proves how NQEs significantly affect properties and features of the air-water interface, that are essential to accurately describe H-bonded interfacial systems.
Collapse
Affiliation(s)
- Fabrizio Creazzo
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| |
Collapse
|
2
|
Hartke B. On the brink of self-hydration: the water heptadecamer. Phys Chem Chem Phys 2024; 26:15445-15451. [PMID: 38747364 DOI: 10.1039/d4cp00816b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
For pure, neutral, isolated molecular clusters, (H2O)17 marks the transition from structures with all water molecules on the cluster surface to water self-hydration, i.e., cluster structures around one central water molecule. Getting this right with water model potentials turns out to be challenging. Even the best water potentials currently available, which reproduce collective properties very well, still deliver contradicting results for (H2O)17, when different low-energy isomers from global structure optimizations are examined. Interestingly, ab initio quantum chemistry also struggles with the only seemingly simple question if (H2O)17 is all-surface or water-centered. Hence, although the long history of water potential development may be entering its final phase, it is not quite finished yet.
Collapse
Affiliation(s)
- Bernd Hartke
- Institute for Physical Chemistry, Kiel University, 24118 Kiel, Germany.
| |
Collapse
|
3
|
Ma R, Baradwaj N, Nomura KI, Krishnamoorthy A, Kalia RK, Nakano A, Vashishta P. Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results. J Chem Phys 2024; 160:134309. [PMID: 38568947 DOI: 10.1063/5.0186058] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2023] [Accepted: 03/18/2024] [Indexed: 04/05/2024] Open
Abstract
Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+
Collapse
Affiliation(s)
- Ruru Ma
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
| | - Nitish Baradwaj
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
| | - Ken-Ichi Nomura
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
| | - Aravind Krishnamoorthy
- Department of Mechanical Engineering, Texas A&M University, College Station, Texas 77843-3123, USA
| | - Rajiv K Kalia
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
| | - Aiichiro Nakano
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
| | - Priya Vashishta
- Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
| |
Collapse
|
4
|
Herman KM, Stone AJ, Xantheas SS. Accurate Calculation of Many-Body Energies in Water Clusters Using a Classical Geometry-Dependent Induction Model. J Chem Theory Comput 2023; 19:6805-6815. [PMID: 37703063 DOI: 10.1021/acs.jctc.3c00575] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/14/2023]
Abstract
We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The moment expansion is carried out up to the hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA classical interaction potentials for water. The classical induction model with no adjustable parameters reproduces the ab initio 3-/4-body terms contained in these two databases with a root-mean-square error (RMSE) of 0.104/0.058 and a mean-absolute error (MAE) of 0.054/0.026 kcal/mol, respectively. These results are on par with the ones obtained by fitting the same data using over 14,000 (for the 3-body) and 200 (for the 4-body) parameters via Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy, included in the calculation of the 3-body energy, was found to be small but not quite negligible. The model represents a practical, efficient, and transferable approach for obtaining accurate nonadditive interactions for multicomponent systems without the need to perform tens of thousands of high level electronic structure calculations and fitting them with PIPs.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98185, United States
| | - Anthony J Stone
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, U.K
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98185, United States
- Advanced Computing Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MSIN J7-10, Richland, Washington 99352, United States
| |
Collapse
|
5
|
Hasegawa T. Nuclear Quantum Dynamics of Three-Dimensional Condensed-Phase Systems by Constant Uncertainty Molecular Dynamics. J Phys Chem Lett 2023; 14:8043-8049. [PMID: 37651720 DOI: 10.1021/acs.jpclett.3c01254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Constant uncertainty molecular dynamics (CUMD), which was developed to include nuclear quantum effects (NQEs) in molecular dynamics (MD) simulations, is extended to three-dimensional condensed-phase systems. Its applicability was verified via CUMD simulations of bulk water and ice Ih with the q-TIP4P/F potential model. The simulated radial distribution functions, proton momentum distributions, and infrared spectra were compared with those of previous studies in which the simulations were carried out on the basis of the imaginary-time path-integral formalism. We found that CUMD can function as an alternative trajectory-based simulation method for including static and dynamic NQEs in MD simulations of three-dimensional condensed-phase systems.
Collapse
Affiliation(s)
- Taisuke Hasegawa
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| |
Collapse
|
6
|
Yu Q, Qu C, Houston PL, Nandi A, Pandey P, Conte R, Bowman JM. A Status Report on "Gold Standard" Machine-Learned Potentials for Water. J Phys Chem Lett 2023; 14:8077-8087. [PMID: 37656898 PMCID: PMC10510435 DOI: 10.1021/acs.jpclett.3c01791] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 08/28/2023] [Indexed: 09/03/2023]
Abstract
Owing to the central importance of water to life as well as its unusual properties, potentials for water have been the subject of extensive research over the past 50 years. Recently, five potentials based on different machine learning approaches have been reported that are at or near the "gold standard" CCSD(T) level of theory. The development of such high-level potentials enables efficient and accurate simulations of water systems using classical and quantum dynamical approaches. This Perspective serves as a status report of these potentials, focusing on their methodology and applications to water systems across different phases. Their performances on the energies of gas phase water clusters, as well as condensed phase structural and dynamical properties, are discussed.
Collapse
Affiliation(s)
- Qi Yu
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chen Qu
- Independent
Researcher, Toronto, Ontario M9B 0E3, Canada
| | - Paul L. Houston
- Department
of Chemistry and Chemical Biology, Cornell
University, Ithaca, New York 14853, United States
- Department of Chemistry
and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Apurba Nandi
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511, Luxembourg City, Luxembourg
| | - Priyanka Pandey
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Riccardo Conte
- Dipartimento
di Chimica, Università degli Studi
di Milano, via Golgi 19, 20133 Milano, Italy
| | - Joel M. Bowman
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
7
|
Riera M, Knight C, Bull-Vulpe EF, Zhu X, Agnew H, Smith DGA, Simmonett AC, Paesani F. MBX: A many-body energy and force calculator for data-driven many-body simulations. J Chem Phys 2023; 159:054802. [PMID: 37526156 PMCID: PMC10550339 DOI: 10.1063/5.0156036] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 07/11/2023] [Indexed: 08/02/2023] Open
Abstract
Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.
Collapse
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Argonne National Laboratory, Computational Science Division, Lemont, Illinois 60439, USA
| | - Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Xuanyu Zhu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | | | - Andrew C. Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | | |
Collapse
|
8
|
Chen MS, Lee J, Ye HZ, Berkelbach TC, Reichman DR, Markland TE. Data-Efficient Machine Learning Potentials from Transfer Learning of Periodic Correlated Electronic Structure Methods: Liquid Water at AFQMC, CCSD, and CCSD(T) Accuracy. J Chem Theory Comput 2023; 19:4510-4519. [PMID: 36730728 DOI: 10.1021/acs.jctc.2c01203] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Obtaining the atomistic structure and dynamics of disordered condensed-phase systems from first-principles remains one of the forefront challenges of chemical theory. Here we exploit recent advances in periodic electronic structure and provide a data-efficient approach to obtain machine-learned condensed-phase potential energy surfaces using AFQMC, CCSD, and CCSD(T) from a very small number (≤200) of energies by leveraging a transfer learning scheme starting from lower-tier electronic structure methods. We demonstrate the effectiveness of this approach for liquid water by performing both classical and path integral molecular dynamics simulations on these machine-learned potential energy surfaces. By doing this, we uncover the interplay of dynamical electron correlation and nuclear quantum effects across the entire liquid range of water while providing a general strategy for efficiently utilizing periodic correlated electronic structure methods to explore disordered condensed-phase systems.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California94305, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York10027, United States
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York10027, United States
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York10010, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York10027, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California94305, United States
| |
Collapse
|
9
|
Heindel JP, Herman KM, Xantheas SS. Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development. Annu Rev Phys Chem 2023; 74:337-360. [PMID: 37093659 DOI: 10.1146/annurev-physchem-062422-023532] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born-Oppenheimer potential energy surface. In particular, we focus on ab initio-based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through (a) parameterization of distributed multipoles, (b) explicit fitting of the MBE, (c) inclusion of many-atom features in a neural network, and (d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington, USA; ,
| |
Collapse
|
10
|
Herman KM, Xantheas SS. An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2-25. Phys Chem Chem Phys 2023; 25:7120-7143. [PMID: 36853239 DOI: 10.1039/d2cp03241d] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2-11, 16-17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (De's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2-25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83-119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (De's) throughout the entire cluster range (n = 2-25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ∼4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6-3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA.
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA. .,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, WA, 99352, USA.
| |
Collapse
|
11
|
Bowman JM, Qu C, Conte R, Nandi A, Houston PL, Yu Q. Δ-Machine Learned Potential Energy Surfaces and Force Fields. J Chem Theory Comput 2023; 19:1-17. [PMID: 36527383 DOI: 10.1021/acs.jctc.2c01034] [Citation(s) in RCA: 22] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
There has been great progress in developing machine-learned potential energy surfaces (PESs) for molecules and clusters with more than 10 atoms. Unfortunately, this number of atoms generally limits the level of electronic structure theory to less than the "gold standard" CCSD(T) level. Indeed, for the well-known MD17 dataset for molecules with 9-20 atoms, all of the energies and forces were obtained with DFT calculations (PBE). This Perspective is focused on a Δ-machine learning method that we recently proposed and applied to bring DFT-based PESs to close to CCSD(T) accuracy. This is demonstrated for hydronium, N-methylacetamide, acetyl acetone, and ethanol. For 15-atom tropolone, it appears that special approaches (e.g., molecular tailoring, local CCSD(T)) are needed to obtain the CCSD(T) energies. A new aspect of this approach is the extension of Δ-machine learning to force fields. The approach is based on many-body corrections to polarizable force field potentials. This is examined in detail using the TTM2.1 water potential. The corrections make use of our recent CCSD(T) datasets for 2-b, 3-b, and 4-b interactions for water. These datasets were used to develop a new fully ab initio potential for water, termed q-AQUA.
Collapse
Affiliation(s)
- Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chen Qu
- Independent Researcher, Toronto, Canada 66777
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
12
|
Jónsson EÖ, Rasti S, Galynska M, Meyer J, Jónsson H. Transferable Potential Function for Flexible H 2O Molecules Based on the Single-Center Multipole Expansion. J Chem Theory Comput 2022; 18:7528-7543. [PMID: 36395502 DOI: 10.1021/acs.jctc.2c00598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
A potential function is presented for describing a system of flexible H2O molecules based on the single-center multipole expansion (SCME) of the electrostatic interaction. The model, referred to as SCME/f, includes the variation of the molecular quadrupole moment as well as the dipole moment with changes in bond length and angle so as to reproduce results of high-level electronic structure calculations. The multipole expansion also includes fixed octupole and hexadecapole moments, as well as anisotropic dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole polarizability tensors. The model contains five adjustable parameters related to the repulsive interaction and damping functions in the electrostatic and dispersion interactions. Their values are adjusted to reproduce the lowest energy isomers of small clusters, (H2O)n with n = 2-6, as well as measured properties of the ice Ih crystal. Subsequent calculations of the energy difference between the various isomer configurations of the clusters show that SCME/f gives good agreement with results of electronic structure calculations and represents a significant improvement over the previously presented rigid SCME potential function. Analysis of the vibrational frequencies of the clusters and structural properties of ice Ih crystal show the importance of accurately describing the variation of the quadrupole moment with molecular structures.
Collapse
Affiliation(s)
- Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| | - Soroush Rasti
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RALeiden, The Netherlands
| | - Marta Galynska
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RALeiden, The Netherlands
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| |
Collapse
|
13
|
Rodríguez-Segundo R, Gijón A, Prosmiti R. Quantum molecular simulations of micro-hydrated halogen anions. Phys Chem Chem Phys 2022; 24:14964-14974. [PMID: 35686995 DOI: 10.1039/d2cp01396g] [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 report the results of a detailed and accurate investigation focused on structures and energetics of poly-hydrated halides employing first-principles polarizable halide-water potentials to describe the underlying forces. Following a bottom-up data-driven potential approach, we initially looked into the classical behavior of higher-order X-(H2O)N clusters. We have located several low-lying energies, such as global and local minima, structures for each cluster, with various water molecules (up to N = 8) surrounding the halide anion (X- = F-, Cl-, Br-, I-), employing an evolutionary programming method. It is found that the F--water clusters exhibit different structural configurations than the heavier halides, however independently of the halide anion, all clusters show in general a selective growth with the anion preferring to be connected to the outer shell of the water molecule arrangements. In turn, path-integral molecular dynamics simulations are performed to incorporate explicitly nuclear quantum and thermal effects in describing the nature of halide ion microsolvation in such prototypical model systems. Our data reveal that at low finite temperatures, nuclear quantum effects affect certain structural properties, such as weakening hydrogen bonding between the halide anion and water molecules, with minor distortions in the water network beyond the first hydration shell, indicating local structure rearrangements. Such structural characteristics and the promising cluster size trends observed in the single-ion solvation energies motivated us to draw connections of small size cluster data to the limits of continuum bulk values, toward the investigation of the challenging computational modeling of bulk single ion hydration.
Collapse
Affiliation(s)
- Raúl Rodríguez-Segundo
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain. .,Atelgraphics S.L., Mota de Cuervo 42, 28043, Madrid, Spain
| | - Alfonso Gijón
- Materials Science Institute of Madrid (ICMM-CSIC), CSIC, Sor Juana Inés de la Cruz 3, 28049 Madrid, Spain
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
14
|
Yu Q, Qu C, Houston PL, Conte R, Nandi A, Bowman JM. q-AQUA: A Many-Body CCSD(T) Water Potential, Including Four-Body Interactions, Demonstrates the Quantum Nature of Water from Clusters to the Liquid Phase. J Phys Chem Lett 2022; 13:5068-5074. [PMID: 35652912 DOI: 10.1021/acs.jpclett.2c00966] [Citation(s) in RCA: 37] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Many model potential energy surfaces (PESs) have been reported for water; however, none are strictly from "first-principles". Here we report such a potential, based on a many-body representation at the CCSD(T) level of theory up to the four-body interaction. The new PES is benchmarked for the isomers of the water hexamer for dissociation energies, harmonic frequencies, and unrestricted diffusion Monte Carlo (DMC) calculations of the zero-point energies of the Prism, Book, and Cage isomers. Dissociation energies of several isomers of the 20-mer agree well with recent benchmark energies. Exploratory DMC calculations on this cluster verify the robustness of the new PES for quantum simulations. The accuracy and speed of the new PES are demonstrated for standard condensed phase properties, i.e., the radial distribution function and the self-diffusion constant. Quantum effects are shown to be substantial for these observables and also needed to bring theory into excellent agreement with experiment.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Chen Qu
- Independent Researcher, Toronto, Ontario, Canada
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
15
|
A Benchmark Protocol for DFT Approaches and Data-Driven Models for Halide-Water Clusters. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27051654. [PMID: 35268757 PMCID: PMC8924895 DOI: 10.3390/molecules27051654] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Revised: 02/18/2022] [Accepted: 02/26/2022] [Indexed: 11/17/2022]
Abstract
Dissolved ions in aqueous media are ubiquitous in many physicochemical processes, with a direct impact on research fields, such as chemistry, climate, biology, and industry. Ions play a crucial role in the structure of the surrounding network of water molecules as they can either weaken or strengthen it. Gaining a thorough understanding of the underlying forces from small clusters to bulk solutions is still challenging, which motivates further investigations. Through a systematic analysis of the interaction energies obtained from high-level electronic structure methodologies, we assessed various dispersion-corrected density functional approaches, as well as ab initio-based data-driven potential models for halide ion-water clusters. We introduced an active learning scheme to automate the generation of optimally weighted datasets, required for the development of efficient bottom-up anion-water models. Using an evolutionary programming procedure, we determined optimized and reference configurations for such polarizable and first-principles-based representation of the potentials, and we analyzed their structural characteristics and energetics in comparison with estimates from DF-MP2 and DFT+D quantum chemistry computations. Moreover, we presented new benchmark datasets, considering both equilibrium and non-equilibrium configurations of higher-order species with an increasing number of water molecules up to 54 for each F, Cl, Br, and I anions, and we proposed a validation protocol to cross-check methods and approaches. In this way, we aim to improve the predictive ability of future molecular computer simulations for determining the ongoing conflicting distribution of different ions in aqueous environments, as well as the transition from nanoscale clusters to macroscopic condensed phases.
Collapse
|
16
|
Yue S, Riera M, Ghosh R, Panagiotopoulos AZ, Paesani F. Transferability of data-driven, many-body models for CO2 simulations in the vapor and liquid phases. J Chem Phys 2022; 156:104503. [DOI: 10.1063/5.0080061] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shuwen Yue
- Princeton University, United States of America
| | - Marc Riera
- Chemistry and Biochemistry, University of California San Diego Department of Chemistry and Biochemistry, United States of America
| | - Raja Ghosh
- University of California San Diego, United States of America
| | | | | |
Collapse
|
17
|
Omodemi O, Sprouse S, Herbert D, Kaledin M, Kaledin AL. On the Cartesian Representation of the Molecular Polarizability Tensor Surface by Polynomial Fitting to Ab Initio Data. J Chem Theory Comput 2021; 18:37-45. [PMID: 34958587 DOI: 10.1021/acs.jctc.1c01015] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We describe an approach to constructing an analytic Cartesian representation of the molecular dipole polarizability tensor surface in terms of polynomials in interatomic distances with a training set of ab initio data points obtained from a molecular dynamics (MD) simulation or by any other available means. The proposed formulation is based on a perturbation treatment of the unmodified point dipole polarizability model of Applequist [ J. Am. Chem. Soc. 1972, 94, 2952] and is shown here to be, by construction (i) free of short-range or other singularities or discontinuities, (ii) symmetric and translationally invariant, and (iii) nonreliant on a body-fixed coordinate system. Permutational invariance of like nuclei is demonstrated to be readily applicable, making this approach useful for highly fluxional and reactive systems. Derivation of the method is described in detail, adding brief didactic numerical examples of H2 and H2O and concluding with an MD simulation of the Raman spectrum of H5O2+ at 300 K with the polarizability tensor fitted to CCSD(T)/aug-cc-pVTZ data obtained using the HBB-4B potential [ J. Chem. Phys. 2005, 122, 044308].
Collapse
Affiliation(s)
- Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Sarah Sprouse
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Destyni Herbert
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Alexey L Kaledin
- Cherry L. Emerson Center for Scientific Computation, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| |
Collapse
|
18
|
Heindel JP, Xantheas SS. Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD). J Chem Theory Comput 2021; 17:7341-7352. [PMID: 34723531 DOI: 10.1021/acs.jctc.1c00780] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-molecular dynamics (MD) classical and nuclear quantum dynamical simulations, in which the energies and forces of the full system are approximated by the two-, three-, and four-body terms of the MBE, and compare the average potential and the vibrational density of states with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the four-body term converges with near-identical behavior to the one from the full simulation. The three-body makes a substantial contribution (∼20%) to the energy, whereas the four-body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (∼2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm-1) at the four-body level and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a four-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N4, N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that the MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States.,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| |
Collapse
|
19
|
Nandi A, Qu C, Houston PL, Conte R, Yu Q, Bowman JM. A CCSD(T)-Based 4-Body Potential for Water. J Phys Chem Lett 2021; 12:10318-10324. [PMID: 34662138 DOI: 10.1021/acs.jpclett.1c03152] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
High-level, ab initio calculations find that the 4-body (4-b) interaction is needed to account for near-100% of the total interaction energy for water clusters as large as the 21-mer. Motivated by this, we report a permutationally invariant polynomial potential energy surface (PES) for the 4-body interaction. This machine-learned PES is a fit to 2119 symmetry-unique, CCSD(T)-F12a/haTZ 4-b interaction energies. Configurations for these come from tetramer direct-dynamics calculations, fragments from an MD water simulation at 300 K, and tetramer fragments in a variety of water clusters. The PIP basis is purified to ensure that the PES goes rigorously to zero in monomer+trimer and dimer+dimer dissociations. The 4-b energies of isomers of the hexamer calculated with the new PES are shown to be in better agreement with benchmark CCSD(T) results than those from the MB-pol potential. Tests on larger clusters further validate the high-fidelity of the PES. The PES is shown to be fast to evaluate, taking 2.4 s for 105 evaluations on a single core of 2.4 GHz Intel Xeon processor, and significantly faster using a parallel version of the PES.
Collapse
Affiliation(s)
- Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chen Qu
- Department of Chemistry & Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
20
|
Bull-Vulpe EF, Riera M, Götz AW, Paesani F. MB-Fit: Software infrastructure for data-driven many-body potential energy functions. J Chem Phys 2021; 155:124801. [PMID: 34598567 DOI: 10.1063/5.0063198] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.
Collapse
Affiliation(s)
- Ethan F Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - 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
| |
Collapse
|
21
|
Lambros E, Dasgupta S, Palos E, Swee S, Hu J, Paesani F. General Many-Body Framework for Data-Driven Potentials with Arbitrary Quantum Mechanical Accuracy: Water as a Case Study. J Chem Theory Comput 2021; 17:5635-5650. [PMID: 34370954 DOI: 10.1021/acs.jctc.1c00541] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a general framework for the development of data-driven many-body (MB) potential energy functions (MB-QM PEFs) that represent the interactions between small molecules at an arbitrary quantum-mechanical (QM) level of theory. As a demonstration, a family of MB-QM PEFs for water is rigorously derived from density functionals belonging to different rungs across Jacob's ladder of approximations within density functional theory (MB-DFT) and from Møller-Plesset perturbation theory (MB-MP2). Through a systematic analysis of individual MB contributions to the interaction energies of water clusters, we demonstrate that all MB-QM PEFs preserve the same accuracy as the corresponding ab initio calculations, with the exception of those derived from density functionals within the generalized gradient approximation (GGA). The differences between the DFT and MB-DFT results are traced back to density-driven errors that prevent GGA functionals from accurately representing the underlying molecular interactions for different cluster sizes and hydrogen-bonding arrangements. We show that this shortcoming may be overcome, within the MB formalism, by using density-corrected functionals (DC-DFT) that provide a more consistent representation of each individual MB contribution. This is demonstrated through the development of a MB-DFT PEF derived from DC-PBE-D3 data, which more accurately reproduce the corresponding ab initio results.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Steven Swee
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Jie Hu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
22
|
Muniz MC, Gartner TE, Riera M, Knight C, Yue S, Paesani F, Panagiotopoulos AZ. Vapor-liquid equilibrium of water with the MB-pol many-body potential. J Chem Phys 2021; 154:211103. [PMID: 34240989 DOI: 10.1063/5.0050068] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. These results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.
Collapse
Affiliation(s)
- Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Thomas E Gartner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Shuwen Yue
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | | |
Collapse
|
23
|
Cruzeiro VWD, Lambros E, Riera M, Roy R, Paesani F, Götz AW. Highly Accurate Many-Body Potentials for Simulations of N 2O 5 in Water: Benchmarks, Development, and Validation. J Chem Theory Comput 2021; 17:3931-3945. [PMID: 34029079 DOI: 10.1021/acs.jctc.1c00069] [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/14/2022]
Abstract
Dinitrogen pentoxide (N2O5) is an important intermediate in the atmospheric chemistry of nitrogen oxides. Although there has been much research, the processes that govern the physical interactions between N2O5 and water are still not fully understood at a molecular level. Gaining a quantitative insight from computer simulations requires going beyond the accuracy of classical force fields while accessing length scales and time scales that are out of reach for high-level quantum-chemical approaches. To this end, we present the development of MB-nrg many-body potential energy functions for nonreactive simulations of N2O5 in water. This MB-nrg model is based on electronic structure calculations at the coupled cluster level of theory and is compatible with the successful MB-pol model for water. It provides a physically correct description of long-range many-body interactions in combination with an explicit representation of up to three-body short-range interactions in terms of multidimensional permutationally invariant polynomials. In order to further investigate the importance of the underlying interactions in the model, a TTM-nrg model was also devised. TTM-nrg is a more simplistic representation that contains only two-body short-range interactions represented through Born-Mayer functions. In this work, an active learning approach was employed to efficiently build representative training sets of monomer, dimer, and trimer structures, and benchmarks are presented to determine the accuracy of our new models in comparison to a range of density functional theory methods. By assessing the binding curves, distortion energies of N2O5, and interaction energies in clusters of N2O5 and water, we evaluate the importance of two-body and three-body short-range potentials. The results demonstrate that our MB-nrg model has high accuracy with respect to the coupled cluster reference, outperforms current density functional theory models, and thus enables highly accurate simulations of N2O5 in aqueous environments.
Collapse
Affiliation(s)
- Vinícius Wilian D Cruzeiro
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States.,Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Ronak Roy
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States.,Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
24
|
Riera M, Hirales A, Ghosh R, Paesani F. Data-Driven Many-Body Models with Chemical Accuracy for CH4/H2O Mixtures. J Phys Chem B 2020; 124:11207-11221. [DOI: 10.1021/acs.jpcb.0c08728] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Alan Hirales
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Raja Ghosh
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
25
|
Lynch C, Rao S, Sansom MSP. Water in Nanopores and Biological Channels: A Molecular Simulation Perspective. Chem Rev 2020; 120:10298-10335. [PMID: 32841020 PMCID: PMC7517714 DOI: 10.1021/acs.chemrev.9b00830] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Indexed: 12/18/2022]
Abstract
This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.
Collapse
Affiliation(s)
- Charlotte
I. Lynch
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Shanlin Rao
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| |
Collapse
|
26
|
Lambros E, Paesani F. How good are polarizable and flexible models for water: Insights from a many-body perspective. J Chem Phys 2020; 153:060901. [DOI: 10.1063/5.0017590] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, 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
- Materials Science and Engineering, 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
|
27
|
Pan C, Liu C, Peng J, Ren P, Huang X. Three-site and five-site fixed-charge water models compatible with AMOEBA force field. J Comput Chem 2020; 41:1034-1044. [PMID: 31976572 DOI: 10.1002/jcc.26151] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2019] [Accepted: 01/01/2020] [Indexed: 11/06/2022]
Abstract
In a typical biomolecular simulation using Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field, the vast majority molecules in the simulation box consist of water, and these water molecules consume the most CPU power due to the explicit mutual induction effect. To improve the computational efficiency, we here develop two new nonpolarizable water models (with flexible bonds and fixed charges) that are compatible with AMOEBA solute: the 3-site AW3C and 5-site AW5C. To derive the force-field parameters for AW3C and AW5C, we fit to six experimental liquid thermodynamic properties: liquid density, enthalpy of vaporization, dielectric constant, isobaric heat capacity, isothermal compressibility and thermal expansion coefficient, at a broad range of temperatures from 261.15 to 353.15 K under 1.0 atm pressure. We further validate our AW3C and AW5C water models by showing that they can well reproduce the radial distribution function g(r), self-diffusion constant D, and hydration free energy from the AMOEBA03 water model and the experimental observations. Furthermore, we show that our AW3C and AW5C water models can greatly accelerate (>5 times) the bulk water as well as biomolecular simulations when compared to AMOEBA water. Specifically, we demonstrate that the applications of AW3C and AW5C water models to simulate a DNA duplex lead to a threefold acceleration, and in the meanwhile well maintain the structural properties as the fully polarizable AMOEBA water. We expect that our AW3C and AW5C water models hold great promise to be widely applied to simulate complex bio-molecules using the AMOEBA force field.
Collapse
Affiliation(s)
- Cong Pan
- Department of Chemistry, The Hong Kong University of Science and Technology, Hong Kong
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas
| | - Junhui Peng
- Department of Chemistry, The Hong Kong University of Science and Technology, Hong Kong
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Hong Kong.,Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| |
Collapse
|
28
|
Zhai Y, Caruso A, Gao S, Paesani F. Active learning of many-body configuration space: Application to the Cs+–water MB-nrg potential energy function as a case study. J Chem Phys 2020; 152:144103. [DOI: 10.1063/5.0002162] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Yaoguang Zhai
- Department of Computer Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
| | - Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Sicun Gao
- Department of Computer Science and Engineering, 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
- Materials Science and Engineering, 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
|
29
|
Riera M, Yeh EP, Paesani F. Data-Driven Many-Body Models for Molecular Fluids: CO2/H2O Mixtures as a Case Study. J Chem Theory Comput 2020; 16:2246-2257. [DOI: 10.1021/acs.jctc.9b01175] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Eric P. Yeh
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
30
|
Yu Q, Bowman JM. Tracking Hydronium/Water Stretches in Magic H3O+(H2O)20 Clusters through High-level Quantum VSCF/VCI Calculations. J Phys Chem A 2020; 124:1167-1175. [DOI: 10.1021/acs.jpca.9b11983] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Qi Yu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
31
|
Liu C, Piquemal JP, Ren P. Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential. J Phys Chem Lett 2020; 11:419-426. [PMID: 31865706 PMCID: PMC7384396 DOI: 10.1021/acs.jpclett.9b03489] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Molecular dynamics (MD) simulations employing classical force fields (FFs) have been widely used to model molecular systems. The important ingredient of the current FFs, atomic charge, remains fixed during MD simulations despite the atomic environment or local geometry changes. This approximation hinders the transferability of the potential being used in multiple phases. Here we implement a geometry-dependent charge flux (GDCF) model into the multipole-based AMOEBA+ polarizable potential. The CF in the current work explicitly depends on the local geometry (bond and angle) of the molecule. To our knowledge, this is the first study that derives energy and force expressions due to GDCF in a multipole-based polarizable FF framework. Due to the inclusion of GDCF, the AMOEBA+ water model is noticeably improved in terms of describing the monomer properties, cluster binding/interaction energy, and a variety of liquid properties, including the infrared spectra that previous flexible water models were not able to capture.
Collapse
Affiliation(s)
- Chengwen Liu
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
- Laboratoire de Chimie Théorique , Sorbonne Université, UMR7616 CNRS , 75252 Paris , France
- Institut Universitaire de France , 75005 , Paris , France
| | - Pengyu Ren
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
32
|
Wu Y, Car R. Quantum momentum distribution and quantum entanglement in the deep tunneling regime. J Chem Phys 2020; 152:024106. [PMID: 31941303 DOI: 10.1063/1.5133053] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we consider the momentum operator of a quantum particle directed along the displacement of two of its neighbors. A modified open-path path integral molecular dynamics is presented to sample the distribution of this directional momentum distribution, where we derive and use a new estimator for this distribution. Variationally enhanced sampling is used to obtain this distribution for an example molecule, malonaldehyde, in the very low temperature regime where deep tunneling happens. We find no secondary feature in the directional momentum distribution and that its absence is due to quantum entanglement through a further study of the reduced density matrix.
Collapse
Affiliation(s)
- Yantao Wu
- The Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
| | - Roberto Car
- The Department of Chemistry and the Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
33
|
Yong X, Burnham CJ, English NJ, Tse JS. Classical and path-integral molecular-dynamics study on liquid water and ice melting using non-empirical TTM2.1-F model. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1652774] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Xue Yong
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Canada
| | - Christian J. Burnham
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | - John S. Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Canada
| |
Collapse
|
34
|
Bizzarro BB, Egan CK, Paesani F. Nature of Halide–Water Interactions: Insights from Many-Body Representations and Density Functional Theory. J Chem Theory Comput 2019; 15:2983-2995. [DOI: 10.1021/acs.jctc.9b00064] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
35
|
Lao KU, Jia J, Maitra R, DiStasio RA. On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions. J Chem Phys 2018; 149:204303. [PMID: 30501247 DOI: 10.1063/1.5051458] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigate how geometric changes influence the static dipole polarizability (α) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-ζ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this ab initio surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order C 6 dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H2O. Using first-order perturbation theory, we compute the zero-point vibrational correction to α at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic (α iso) and nearly 50% to the anisotropic (α aniso) polarizability values. In doing so, we find that α iso = 9.8307 bohr3, which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.
Collapse
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Rahul Maitra
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
36
|
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
|
37
|
Heßelmann A. Correlation effects and many-body interactions in water clusters. Beilstein J Org Chem 2018; 14:979-991. [PMID: 29977369 PMCID: PMC6009095 DOI: 10.3762/bjoc.14.83] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2018] [Accepted: 04/11/2018] [Indexed: 11/23/2022] Open
Abstract
Background: The quantum-chemical description of the interactions in water clusters is an essential basis for deriving accurate and physically sound models of the interaction potential for water to be used in molecular simulations. In particular, the role of many-body interactions beyond the two-body interactions, which are often not explicitly taken into account by empirical force fields, can be accurately described by quantum chemistry methods on an adequate level, e.g., random-phase approximation electron correlation methods. The relative magnitudes of the different interaction energy contributions obtained by accurate ab initio calculations can therefore provide useful insights that can be exploited to develop enhanced force field methods. Results: In line with earlier theoretical studies of the interactions in water clusters, it has been found that the main contribution to the many-body interactions in clusters with a size of up to N = 13 molecules are higher-order polarisation interaction terms. Compared to this, many-body dispersion interactions are practically negligible for all studied sytems. The two-body dispersion interaction, however, plays a significant role in the formation of the structures of the water clusters and their stability, since it leads to a distinct compression of the cluster sizes compared to the structures optimized on an uncorrelated level. Overall, the many-body interactions amount to about 13% of the total interaction energy, irrespective of the cluster size. The electron correlation contribution to these, however, amounts to only about 30% to the total many-body interactions for the largest clusters studied and is repulsive for all structures considered in this work. Conclusion: While this shows that three- and higher-body interactions can not be neglected in the description of water complexes, the electron correlation contributions to these are much smaller in comparison to the two-body electron correlation effects. Efficient quantum chemistry approaches for describing intermolecular interactions between water molecules may therefore describe higher-body interactions on an uncorrelated Hartree-Fock level without a serious loss in accuracy.
Collapse
Affiliation(s)
- Andreas Heßelmann
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, 91058 Erlangen, Germany
| |
Collapse
|
38
|
Affiliation(s)
- Xinzijian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| | - Jian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| |
Collapse
|
39
|
León-Merino I, Rodríguez-Segundo R, Arismendi-Arrieta DJ, Prosmiti R. Assessing Intermolecular Interactions in Guest-Free Clathrate Hydrate Systems. J Phys Chem A 2018; 122:1479-1487. [PMID: 29328645 DOI: 10.1021/acs.jpca.7b12107] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Recently, empty hydrate structures sI, sII, sH, and others have been proposed as low-density ice structures by both experimental observations and computer simulations. Some of them have been synthesized in the laboratory, which motivates further investigations on the stability of such guest-free clathrate structures. Using semiempirical and ab initio-based water models, as well as dispersion-corrected density functional theory approaches, we predict their stability, including cooperative many-body effects, in comparison with reference data from converged wave function-based DF-MP2 electronic structure calculations. We show that large basis sets and counterpoise corrections are required to improve convergence in the interaction/binding energies for such systems. Therefore, extrapolation schemes based on triple/quadruple and quadruple/quintuple ζ quality basis sets are used to reach high accuracy. Eleven different water structures corresponding to dodecahedron, edge sharing, face sharing, and fused cubes, as a part of the WATER27 database, as well as cavities from the sI, sII, and sH clathrate hydrates formed by 20, 24, 28, and 36 water molecules, are employed, and new benchmark energies are reported. Using these benchmark sets of interaction energies, we assess the performance of both analytical models and direct DFT calculations for such clathrate-like systems. In particular, seven popular water models (TIP4P/ice, TIP4P/2005, q-TIP4P/F, TTM2-F, TTM3-F, TTM4-F, and MB-pol) available in the literature, and nine density functional approximations (3 meta-GGAs, 3 hybrids, and 3 range separated functionals) are used to investigate their accuracy. By including dispersion corrections, our results show that errors in the interaction energies are reduced for most of the DFT functionals. Despite the difficulties faced by current water models and DFT functionals to accurately describe the interactions in such water systems, we found some general trends that could serve to extend their applicability to larger systems.
Collapse
Affiliation(s)
- Iván León-Merino
- Institute of Fundamental Physics (IFF-CSIC), CSIC , Serrano 123, 28006 Madrid, Spain
| | | | | | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC , Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
40
|
Demerdash O, Wang L, Head‐Gordon T. Advanced models for water simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1355] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Omar Demerdash
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California, Davis Davis CA USA
| | - Teresa Head‐Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
- Department of Bioengineering University of California Berkeley CA USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley CA USA
| |
Collapse
|
41
|
Riera M, Mardirossian N, Bajaj P, Götz AW, Paesani F. Toward chemical accuracy in the description of ion–water interactions through many-body representations. Alkali-water dimer potential energy surfaces. J Chem Phys 2017; 147:161715. [DOI: 10.1063/1.4993213] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Narbe Mardirossian
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - 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
| |
Collapse
|
42
|
Futera Z, English NJ. Communication: Influence of external static and alternating electric fields on water from long-time non-equilibrium ab initio molecular dynamics. J Chem Phys 2017; 147:031102. [DOI: 10.1063/1.4994694] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Zdenek Futera
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield Dublin 4, Ireland
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield Dublin 4, Ireland
| |
Collapse
|
43
|
Galib M, Duignan TT, Misteli Y, Baer MD, Schenter GK, Hutter J, Mundy CJ. Mass density fluctuations in quantum and classical descriptions of liquid water. J Chem Phys 2017; 146:244501. [DOI: 10.1063/1.4986284] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Mirza Galib
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Timothy T. Duignan
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Yannick Misteli
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Marcel D. Baer
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Gregory K. Schenter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Jürg Hutter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Christopher J. Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| |
Collapse
|
44
|
Affiliation(s)
- Emilie B. Guidez
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Mark S. Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| |
Collapse
|
45
|
Reddy SK, Straight SC, Bajaj P, Huy Pham C, Riera M, Moberg DR, Morales MA, Knight C, Götz AW, Paesani F. On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice. J Chem Phys 2016; 145:194504. [DOI: 10.1063/1.4967719] [Citation(s) in RCA: 161] [Impact Index Per Article: 20.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Sandeep K. Reddy
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Shelby C. Straight
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - C. Huy Pham
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Daniel R. Moberg
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, USA
| | - Chris Knight
- Leadership Computing Facility, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, USA
| | - 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
| |
Collapse
|
46
|
Paesani F. Getting the Right Answers for the Right Reasons: Toward Predictive Molecular Simulations of Water with Many-Body Potential Energy Functions. Acc Chem Res 2016; 49:1844-51. [PMID: 27548325 DOI: 10.1021/acs.accounts.6b00285] [Citation(s) in RCA: 101] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The central role played by water in fundamental processes relevant to different disciplines, including chemistry, physics, biology, materials science, geology, and climate research, cannot be overemphasized. It is thus not surprising that, since the pioneering work by Stillinger and Rahman, many theoretical and computational studies have attempted to develop a microscopic description of the unique properties of water under different thermodynamic conditions. Consequently, numerous molecular models based on either molecular mechanics or ab initio approaches have been proposed over the years. However, despite continued progress, the correct prediction of the properties of water from small gas-phase clusters to the liquid phase and ice through a single molecular model remains challenging. To large extent, this is due to the difficulties encountered in the accurate modeling of the underlying hydrogen-bond network in which both number and strength of the hydrogen bonds vary continuously as a result of a subtle interplay between energetic, entropic, and nuclear quantum effects. In the past decade, the development of efficient algorithms for correlated electronic structure calculations of small molecular complexes, accompanied by tremendous progress in the analytical representation of multidimensional potential energy surfaces, opened the doors to the design of highly accurate potential energy functions built upon rigorous representations of the many-body expansion (MBE) of the interaction energies. This Account provides a critical overview of the performance of the MB-pol many-body potential energy function through a systematic analysis of energetic, structural, thermodynamic, and dynamical properties as well as of vibrational spectra of water from the gas to the condensed phase. It is shown that MB-pol achieves unprecedented accuracy across all phases of water through a quantitative description of each individual term of the MBE, with a physically correct representation of both short- and long-range many-body contributions. Comparisons with experimental data probing different regions of the water potential energy surface from clusters to bulk demonstrate that MB-pol represents a major step toward the long-sought-after "universal model" capable of accurately describing the molecular properties of water under different conditions and in different environments. Along this path, future challenges include the extension of the many-body scheme adopted by MB-pol to the description of generic solutes as well as the integration of MB-pol in an efficient theoretical and computational framework to model acid-base reactions in aqueous environments. In this context, given the nontraditional form of the MB-pol energy and force expressions, synergistic efforts by theoretical/computational chemists/physicists and computer scientists will be critical for the development of high-performance software for many-body molecular dynamics simulations.
Collapse
Affiliation(s)
- Francesco Paesani
- Department of Chemistry and
Biochemistry, University of California—San Diego, La Jolla, California 92093, United States
| |
Collapse
|
47
|
|
48
|
Verstraelen T, Vandenbrande S, Heidar-Zadeh F, Vanduyfhuys L, Van Speybroeck V, Waroquier M, Ayers PW. Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development. J Chem Theory Comput 2016; 12:3894-912. [DOI: 10.1021/acs.jctc.6b00456] [Citation(s) in RCA: 83] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Toon Verstraelen
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Steven Vandenbrande
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Farnaz Heidar-Zadeh
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 West
Main Street, Hamilton, Ontario L8S 4M1, Canada
| | - Louis Vanduyfhuys
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Veronique Van Speybroeck
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Michel Waroquier
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Paul W. Ayers
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 West
Main Street, Hamilton, Ontario L8S 4M1, Canada
| |
Collapse
|
49
|
Cisneros G, Wikfeldt KT, Ojamäe L, Lu J, Xu Y, Torabifard H, Bartók AP, Csányi G, Molinero V, Paesani F. Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions. Chem Rev 2016; 116:7501-28. [PMID: 27186804 PMCID: PMC5450669 DOI: 10.1021/acs.chemrev.5b00644] [Citation(s) in RCA: 278] [Impact Index Per Article: 34.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2015] [Indexed: 12/17/2022]
Abstract
Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.
Collapse
Affiliation(s)
| | - Kjartan Thor Wikfeldt
- Science
Institute, University of Iceland, VR-III, 107, Reykjavik, Iceland
- Department
of Physics, Albanova, Stockholm University, S-106 91 Stockholm, Sweden
| | - Lars Ojamäe
- Department
of Chemistry, Linköping University, SE-581 83 Linköping, Sweden
| | - Jibao Lu
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Yao Xu
- Lehrstuhl
Physikalische Chemie II, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Hedieh Torabifard
- Department
of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Albert P. Bartók
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Gábor Csányi
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
50
|
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
|