1
|
Shen H, Shen X, Wu Z. Simulating the isotropic Raman spectra of O-H stretching mode in liquid H 2O based on a machine learning potential: the influence of vibrational couplings. Phys Chem Chem Phys 2023; 25:28180-28188. [PMID: 37819214 DOI: 10.1039/d3cp03035k] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/13/2023]
Abstract
In this study, we trained a deep potential (DP) for H2O, an accurate machine learning (ML) potential. We performed molecular dynamics (MD) simulations of liquid water using the DP model (or DeePMD simulations). Our results showed that the DP model exhibits DFT-level accuracy, and the DeePMD simulation is a promising approach for modeling the structural properties of liquid water. Based on the DeePMD simulation trajectories, we calculated the isotropic Raman spectra of the O-H stretching mode using the surface-specific velocity-velocity correlation function (ssVVCF), showing that the DeePMD/ssVVCF approach can correctly capture the bimodal characteristics of the experimental Raman spectra, with one peak located near 3400 cm-1 and the other near 3250 cm-1. The success of the DeePMD/ssVVCF approach should be credited to (1) the DFT-level accuracy of the DP model for H2O, (2) the ssVVCF formulation considering the coupling between vibrational modes, and (3) non-Condon effects. Furthermore, the DeePMD simulations revealed that the anharmonic interactions between the coupled water molecules in the first and second hydration shells should play an essential role in the strong mixing of the H-O-H bending mode and the O-H stretching mode, leading to the delocalization of the O-H stretching band. In particular, increasing the strength of hydrogen bonds would enhance the bend-stretch coupling, leading to the red-shifting of the O-H vibrational spectra and the increase in the intensity of the shoulder around 3250 cm-1.
Collapse
Affiliation(s)
- Hujun Shen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Education University, Guiyang 550018, China.
| | - Xu Shen
- National Center of Technology Innovation for Intelligent Design and Numerical Control, Huazhong University of Science and Technology, Wuhan, 430074, China
| | - Zhenhua Wu
- Department of Big Data and Artificial Intelligence, Guizhou Vocational Technology College of Electronics & Information, Kaili, 556000, China
| |
Collapse
|
2
|
Zhai Y, Caruso A, Bore SL, Luo Z, Paesani F. A "short blanket" dilemma for a state-of-the-art neural network potential for water: Reproducing experimental properties or the physics of the underlying many-body interactions? J Chem Phys 2023; 158:084111. [PMID: 36859071 DOI: 10.1063/5.0142843] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Deep neural network (DNN) potentials have recently gained popularity in computer simulations of a wide range of molecular systems, from liquids to materials. In this study, we explore the possibility of combining the computational efficiency of the DeePMD framework and the demonstrated accuracy of the MB-pol data-driven, many-body potential to train a DNN potential for large-scale simulations of water across its phase diagram. We find that the DNN potential is able to reliably reproduce the MB-pol results for liquid water, but provides a less accurate description of the vapor-liquid equilibrium properties. This shortcoming is traced back to the inability of the DNN potential to correctly represent many-body interactions. An attempt to explicitly include information about many-body effects results in a new DNN potential that exhibits the opposite performance, being able to correctly reproduce the MB-pol vapor-liquid equilibrium properties, but losing accuracy in the description of the liquid properties. These results suggest that DeePMD-based DNN potentials are not able to correctly "learn" and, consequently, represent many-body interactions, which implies that DNN potentials may have limited ability to predict the properties for state points that are not explicitly included in the training process. The computational efficiency of the DeePMD framework can still be exploited to train DNN potentials on data-driven many-body potentials, which can thus enable large-scale, "chemically accurate" simulations of various molecular systems, with the caveat that the target state points must have been adequately sampled by the reference data-driven many-body potential in order to guarantee a faithful representation of the associated properties.
Collapse
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
| | - Sigbjørn Løland Bore
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Zhishang Luo
- 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
| |
Collapse
|
3
|
Jankowski P, Grabowska E, Szalewicz K. On the role of coupled-clusters' full triple and perturbative quadruple excitations on rovibrational spectra of van der Waals complexes. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1955989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | - Ewelina Grabowska
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE, USA
| |
Collapse
|
4
|
Konovalov A, Symons BCB, Popelier PLA. On the many-body nature of intramolecular forces in FFLUX and its implications. J Comput Chem 2020; 42:107-116. [PMID: 33107993 DOI: 10.1002/jcc.26438] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Revised: 09/28/2020] [Accepted: 09/28/2020] [Indexed: 12/24/2022]
Abstract
FFLUX is a biomolecular force field under construction, based on Quantum Chemical Topology (QCT) and machine learning (kriging), with a minimalistic and physically motivated design. A detailed analysis of the forces within the kriging models as treated in FFLUX is presented, taking as a test example a liquid water model. The energies of topological atoms are modeled as 3Natoms -6 dimensional potential energy surfaces, using atomic local frames to represent the internal degrees of freedom. As a result, the forces within the kriging models in FFLUX are inherently N-body in nature where N refers to Natoms . This provides a fuller picture that is closer to a true quantum mechanical representation of interactions between atoms. The presented computational example quantitatively showcases the non-negligible (as much as 9%) three-body nature of bonded forces and angular forces in a water molecule. We discuss the practical impact on the pressure calculation with N-body forces and periodic boundary conditions (PBC) in molecular dynamics, as opposed to classical force fields with two-body forces. The equivalence between the PBC-related correction terms in the general virial equation is shown mathematically.
Collapse
Affiliation(s)
- Anton Konovalov
- Manchester Institute of Biotechnology (MIB), Manchester, United Kingdom.,Department of Chemistry, University of Manchester, Manchester, United Kingdom
| | - Benjamin C B Symons
- Manchester Institute of Biotechnology (MIB), Manchester, United Kingdom.,Department of Chemistry, University of Manchester, Manchester, United Kingdom
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB), Manchester, United Kingdom.,Department of Chemistry, University of Manchester, Manchester, United Kingdom
| |
Collapse
|
5
|
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
|
6
|
Garcia J, Podeszwa R, Szalewicz K. SAPT codes for calculations of intermolecular interaction energies. J Chem Phys 2020; 152:184109. [PMID: 32414261 DOI: 10.1063/5.0005093] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.
Collapse
Affiliation(s)
- Javier Garcia
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia at Katowice, Szkolna 9, Katowice, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
7
|
Metz MP, Szalewicz K. Automatic Generation of Flexible-Monomer Intermolecular Potential Energy Surfaces. J Chem Theory Comput 2020; 16:2317-2339. [PMID: 32240593 DOI: 10.1021/acs.jctc.9b01241] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
A method is developed for automatic generation of nonreactive intermolecular two-body potential energy surfaces (PESs) including intramonomer degrees of freedom. This method, called flex-autoPES, is an extension of the autoPES method developed earlier, which assumes rigid monomers. In both cases, the whole PES development proceeds without any human intervention. The functional form used is a sum of products of site-site functions (both atomic and off-atomic sites can be used). The leading terms with sites involving different monomers are of physically motivated form. The long-range part of a PES is computed from monomer properties without using any dimer information. The close-range part is fitted to dimer interaction energies computed using electronic structure methods. Virtually any method can be used in such calculations, but the use of symmetry-adapted perturbation theory provides a seamless connection to the long-range part of the PES. The performance of the flex-autoPES code was tested by developing a full-dimensional PES for the water dimer and PESs including only some soft intramonomer degrees of freedom for the ethylene glycol dimer and for the ethylene glycol-water dimer. In the case of the water dimer, the root-mean-square error (RMSE) of the PES from the data points with negative total energies is 0.03 kcal/mol, and we expect this PES to be more accurate than any previously published PES of this type. For the ethylene glycol dimer and the ethylene glycol-water dimers, the analogous RMSEs are 0.25 and 0.1 kcal/mol, respectively.
Collapse
Affiliation(s)
- Michael P Metz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
8
|
Zhang XL, Ma YT, Zhai Y, Li H. Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements. J Chem Phys 2019; 151:074301. [PMID: 31438702 DOI: 10.1063/1.5115496] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The rovibrational energy levels and intensities of the CH3F-H2 dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H2O-H2 dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H2 are more strongly bound than dimers containing para-H2. Energies and wavefunctions of the discrete rovibrational levels of CH3F-paraH2 complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH3F-paraH2 complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH3F-orthoH2 dimer has been predicted for the first time.
Collapse
Affiliation(s)
- Xiao-Long Zhang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, China
| | - Yong-Tao Ma
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, China
| | - Yu Zhai
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, China
| | - Hui Li
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, China
| |
Collapse
|
9
|
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]
|
10
|
Zhuang D, Riera M, Schenter GK, Fulton JL, Paesani F. Many-Body Effects Determine the Local Hydration Structure of Cs + in Solution. J Phys Chem Lett 2019; 10:406-412. [PMID: 30629438 DOI: 10.1021/acs.jpclett.8b03829] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
A systematic analysis of the hydration structure of Cs+ ions in solution is derived from simulations carried out using a series of molecular models built upon a hierarchy of approximate representations of many-body effects in ion-water interactions. It is found that a pairwise-additive model, commonly used in biomolecular simulations, provides poor agreement with experimental X-ray spectra, indicating an incorrect description of the underlying hydration structure. Although the agreement with experiment improves in simulations with a polarizable model, the predicted hydration structure is found to lack the correct sequence of water shells. Progressive inclusion of explicit many-body effects in the representation of Cs+-water interactions as well as accounting for nuclear quantum effects is shown to be necessary for quantitatively reproducing the experimental X-ray spectra. Besides emphasizing the importance of many-body effects, these results suggest that molecular models rigorously derived from many-body expansions hold promise for realistic simulations of aqueous solutions.
Collapse
Affiliation(s)
- Debbie Zhuang
- 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
| | - Gregory K Schenter
- Physical and Computational Sciences Directorate , Pacific Northwest National Laboratory , Richland , Washington 99352 , United States
| | - John L Fulton
- Physical and Computational Sciences Directorate , Pacific Northwest National Laboratory , Richland , Washington 99352 , 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
|
11
|
Metz MP, Szalewicz K, Sarka J, Tóbiás R, Császár AG, Mátyus E. Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states. Phys Chem Chem Phys 2019; 21:13504-13525. [DOI: 10.1039/c9cp00993k] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H2O)2, CH4·H2O, and (CH4)2.
Collapse
Affiliation(s)
- Michael P. Metz
- Department of Physics and Astronomy
- University of Delaware
- Newark
- USA
| | | | - János Sarka
- Department of Chemistry and Biochemistry
- Texas Tech University
- Lubbock
- USA
| | - Roland Tóbiás
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
- MTA-ELTE Complex Chemical Systems Research Group
| | - Attila G. Császár
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
- MTA-ELTE Complex Chemical Systems Research Group
| | - Edit Mátyus
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
| |
Collapse
|
12
|
Garberoglio G, Jankowski P, Szalewicz K, Harvey AH. Fully quantum calculation of the second and third virial coefficients of water and its isotopologues from ab initio potentials. Faraday Discuss 2018; 212:467-497. [PMID: 30302450 PMCID: PMC6561489 DOI: 10.1039/c8fd00092a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Path-Integral Monte Carlo methods were applied to calculate the second, B(T), and the third, C(T), virial coefficients for water. A fully quantum approach and state-of-the-art flexible-monomer pair and three-body potentials were used. Flexible-monomer potentials allow calculations for any isotopologue; we performed calculations for both H2O and D2O. For B(T) of H2O, the quantum effect contributes 25% of the value at 300 K and is not entirely negligible even at 1000 K, in accordance with recent literature findings. The effect of monomer flexibility, while not as large as some claims in the literature, is significant compared to the experimental uncertainty. It is of opposite sign to the quantum effect, smaller in magnitude than the latter below 500 K, and varies from 2% at 300 K to 10% at 700 K. When monomer flexibility is accounted for, results from the CCpol-8sf pair potential are in excellent agreement with the available experimental data and provide reliable B(T) values at temperatures outside the range of experimental data. The flexible-monomer MB-pol pair potential yields B(T) values that are slightly too high compared to experiment. For C(T), our calculations confirm earlier findings that the use of three-body potential is necessary for meaningful predictions. However, due to various uncertainties of the potentials used, especially the three-body ones, we were not able to establish benchmark values of C(T), although our results are in qualitative agreement with available experimental data. The quantum effect, never before included for water, reduces the magnitude of the classical value for H2O by a factor of 2.5 at 300 K and is not entirely negligible even at 1000 K.
Collapse
Affiliation(s)
- Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear Physics and Related Areas (FBK-ECT*), strada delle Tabarelle 286, I-38123 Trento, Italy. and Trento Institute for Fundamental Physics and Applications (TIFPA-INFN), via Sommarive 18, I-38213 Trento, Italy
| | - Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University in Torun, Gagarina 7, PL-87-100 Torun, Poland.
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA.
| | - Allan H Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305-3337, USA.
| |
Collapse
|
13
|
Hunter KM, Shakib FA, Paesani F. Disentangling Coupling Effects in the Infrared Spectra of Liquid Water. J Phys Chem B 2018; 122:10754-10761. [DOI: 10.1021/acs.jpcb.8b09910] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Kelly M. Hunter
- Department of Chemistry and Biochemistry, University of California, San Diego La Jolla, California 92093, United States
| | - Farnaz A. Shakib
- 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
| |
Collapse
|
14
|
Riera M, Brown SE, Paesani F. Isomeric Equilibria, Nuclear Quantum Effects, and Vibrational Spectra of M+(H2O)n=1–3 Clusters, with M = Li, Na, K, Rb, and Cs, through Many-Body Representations. J Phys Chem A 2018; 122:5811-5821. [DOI: 10.1021/acs.jpca.8b04106] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
15
|
|
16
|
Analysis of polarization in hydrogen bonded complexes: An asymptotic projection approach. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.02.034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
17
|
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
|
18
|
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
|
19
|
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
|
20
|
Brown SE, Götz AW, Cheng X, Steele RP, Mandelshtam VA, Paesani F. Monitoring Water Clusters “Melt” Through Vibrational Spectroscopy. J Am Chem Soc 2017; 139:7082-7088. [DOI: 10.1021/jacs.7b03143] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
| | | | - Xiaolu Cheng
- Department
of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Ryan P. Steele
- Department
of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | | | | |
Collapse
|
21
|
Hellmann R. Intermolecular potential energy surface and thermophysical properties of propane. J Chem Phys 2017; 146:114304. [PMID: 28330364 DOI: 10.1063/1.4978412] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A six-dimensional potential energy surface (PES) for the interaction of two rigid propane molecules was determined from supermolecular ab initio calculations up to the coupled cluster with single, double, and perturbative triple excitations level of theory for 9452 configurations. An analytical site-site potential function with 14 sites per molecule was fitted to the calculated interaction energies. To validate the analytical PES, the second virial coefficient and the dilute gas shear viscosity and thermal conductivity of propane were computed. The dispersion part of the potential function was slightly adjusted such that quantitative agreement with the most accurate experimental data for the second virial coefficient at room temperature was achieved. The adjusted PES yields values for the three properties that are in very good agreement with the best experimental data at all temperatures.
Collapse
Affiliation(s)
- Robert Hellmann
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| |
Collapse
|
22
|
Pham CH, Reddy SK, Chen K, Knight C, Paesani F. Many-Body Interactions in Ice. J Chem Theory Comput 2017; 13:1778-1784. [DOI: 10.1021/acs.jctc.6b01248] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- C. Huy Pham
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Sandeep K. Reddy
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Karl Chen
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Chris Knight
- Leadership
Computing Facility, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, United States
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| |
Collapse
|
23
|
Affiliation(s)
- Ramachandran Subramanian
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York, USA
| | - Andrew J. Schultz
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York, USA
| | - David A. Kofke
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York, USA
| |
Collapse
|
24
|
Hellmann R. Nonadditive three-body potential and third to eighth virial coefficients of carbon dioxide. J Chem Phys 2017; 146:054302. [DOI: 10.1063/1.4974995] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Robert Hellmann
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
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
|
28
|
Salim MA, Willow SY, Hirata S. Ice Ih anomalies: Thermal contraction, anomalous volume isotope effect, and pressure-induced amorphization. J Chem Phys 2016; 144:204503. [DOI: 10.1063/1.4951687] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael A. Salim
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - Soohaeng Yoo Willow
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| |
Collapse
|
29
|
Zhu W, Zhao WH, Wang L, Yin D, Jia M, Yang J, Zeng XC, Yuan LF. Two-dimensional interlocked pentagonal bilayer ice: how do water molecules form a hydrogen bonding network? Phys Chem Chem Phys 2016; 18:14216-21. [DOI: 10.1039/c5cp07524f] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The tradeoff between the conditions of an ideal hydrogen bonding network can serve as a generic guidance to understand the rich phase behaviors of nanoconfined water.
Collapse
Affiliation(s)
- Weiduo Zhu
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Wen-Hui Zhao
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Lu Wang
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Di Yin
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Min Jia
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Xiao Cheng Zeng
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Lan-Feng Yuan
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| |
Collapse
|
30
|
Crusius JP, Hellmann R, Hassel E, Bich E. Ab initio intermolecular potential energy surface and thermophysical properties of nitrous oxide. J Chem Phys 2015; 142:244307. [PMID: 26133428 DOI: 10.1063/1.4922830] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an analytical intermolecular potential energy surface (PES) for two rigid nitrous oxide (N2O) molecules derived from high-level quantum-chemical ab initio calculations. Interaction energies for 2018 N2O-N2O configurations were computed utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory using basis sets up to aug-cc-pVQZ supplemented with bond functions. A site-site potential function with seven sites per N2O molecule was fitted to the pair interaction energies. We validated our PES by computing the second virial coefficient as well as shear viscosity and thermal conductivity in the dilute-gas limit. The values of these properties are substantiated by the best experimental data.
Collapse
Affiliation(s)
| | - Robert Hellmann
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| | - Egon Hassel
- Lehrstuhl für Technische Thermodynamik, Universität Rostock, 18059 Rostock, Germany
| | - Eckard Bich
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| |
Collapse
|
31
|
Miliordos E, Xantheas SS. An accurate and efficient computational protocol for obtaining the complete basis set limits of the binding energies of water clusters at the MP2 and CCSD(T) levels of theory: Application to (H2O)m, m = 2-6, 8, 11, 16, and 17. J Chem Phys 2015; 142:234303. [DOI: 10.1063/1.4922262] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Evangelos Miliordos
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| | - Sotiris S. Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| |
Collapse
|
32
|
Jankowski P, Murdachaew G, Bukowski R, Akin-Ojo O, Leforestier C, Szalewicz K. Ab initio water pair potential with flexible monomers. J Phys Chem A 2015; 119:2940-64. [PMID: 25687650 DOI: 10.1021/jp512847z] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects.
Collapse
Affiliation(s)
- Piotr Jankowski
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States.,§Faculty of Chemistry, Nicolaus Copernicus University in Torun, Gagarina 7, 87-100 Torun, Poland
| | - Garold Murdachaew
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Robert Bukowski
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Omololu Akin-Ojo
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Claude Leforestier
- ‡Institut Charles Gerhardt (CTMM)-UMR 5253, CC 15.01, Université Montpellier II-CNRS, Montpellier, Cedex 05, France 34095
| | - Krzysztof Szalewicz
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
33
|
Medders GR, Paesani F. Infrared and Raman Spectroscopy of Liquid Water through "First-Principles" Many-Body Molecular Dynamics. J Chem Theory Comput 2015; 11:1145-54. [PMID: 26579763 DOI: 10.1021/ct501131j] [Citation(s) in RCA: 161] [Impact Index Per Article: 17.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Vibrational spectroscopy is a powerful technique to probe the structure and dynamics of water. However, deriving an unambiguous molecular-level interpretation of the experimental spectral features remains a challenge due to the complexity of the underlying hydrogen-bonding network. In this contribution, we present an integrated theoretical and computational framework (named many-body molecular dynamics or MB-MD) that, by systematically removing uncertainties associated with existing approaches, enables a rigorous modeling of vibrational spectra of water from quantum dynamical simulations. Specifically, we extend approaches used to model the many-body expansion of interaction energies to develop many-body representations of the dipole moment and polarizability of water. The combination of these "first-principles" representations with centroid molecular dynamics simulations enables the simulation of infrared and Raman spectra of liquid water under ambient conditions that, without relying on any ad hoc parameters, are in good agreement with the corresponding experimental results. Importantly, since the many-body energy, dipole, and polarizability surfaces employed in the simulations are derived independently from accurate fits to correlated electronic structure data, MB-MD allows for a systematic analysis of the calculated spectra in terms of both electronic and dynamical contributions. The present analysis suggests that, while MB-MD correctly reproduces both the shifts and the shapes of the main spectroscopic features, an improved description of quantum dynamical effects possibly combined with a dissociable water potential may be necessary for a quantitative representation of the OH stretch band.
Collapse
Affiliation(s)
- Gregory R Medders
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92037, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92037, United States
| |
Collapse
|
34
|
Góra U, Cencek W, Podeszwa R, van der Avoird A, Szalewicz K. Predictions for water clusters from a first-principles two- and three-body force field. J Chem Phys 2014; 140:194101. [PMID: 24852524 DOI: 10.1063/1.4875097] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new rigid-monomer three-body potential has been developed for water by fitting it to more than 70 thousand trimer interaction energies computed ab initio using coupled-cluster methods and augmented triple-zeta-quality basis sets. This potential was used together with a modified form of a previously developed two-body potential and with a polarization model of four- and higher-body interactions to predict the energetics of the water trimer, hexamer, and 24-mer. Despite using the rigid-monomer approximation, these predictions agree better with flexible-monomer benchmarks than published results obtained with flexible-monomer force fields. An unexpected finding of our work is that simple polarization models predict four-body interactions to within a few percent, whereas for three-body interactions these models are known to have errors on the order of 50%.
Collapse
Affiliation(s)
- Urszula Góra
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Wojciech Cencek
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Ad van der Avoird
- Theoretical Chemistry, Institute for Molecules and Materials, Radboud University Nijmegen, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
35
|
Arismendi-Arrieta D, Medina JS, Fanourgakis GS, Prosmiti R, Delgado-Barrio G. Simulating liquid water for determining its structural and transport properties. Appl Radiat Isot 2014; 83 Pt B:115-21. [DOI: 10.1016/j.apradiso.2013.01.020] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2012] [Revised: 01/10/2013] [Accepted: 01/10/2013] [Indexed: 11/29/2022]
|
36
|
Babin V, Leforestier C, Paesani F. Development of a “First Principles” Water Potential with Flexible Monomers: Dimer Potential Energy Surface, VRT Spectrum, and Second Virial Coefficient. J Chem Theory Comput 2013; 9:5395-403. [DOI: 10.1021/ct400863t] [Citation(s) in RCA: 318] [Impact Index Per Article: 28.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Volodymyr Babin
- Department
of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Claude Leforestier
- ICG-CTMM
CC 15.01, Universite Montpellier
II , 34095 Montpellier, Cedex 05, France
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
37
|
Medders GR, Paesani F. Many-Body Convergence of the Electrostatic Properties of Water. J Chem Theory Comput 2013; 9:4844-52. [DOI: 10.1021/ct400696d] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gregory R. Medders
- 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
| |
Collapse
|
38
|
|
39
|
Affiliation(s)
- Branko Ruscic
- Chemical
Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States, and Computation Institute, University of Chicago, Chicago, Illinois 60637, United
States
| |
Collapse
|
40
|
Hellmann R. Ab initiopotential energy surface for the nitrogen molecule pair and thermophysical properties of nitrogen gas. Mol Phys 2013. [DOI: 10.1080/00268976.2012.726379] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
41
|
Medders GR, Babin V, Paesani F. A Critical Assessment of Two-Body and Three-Body Interactions in Water. J Chem Theory Comput 2013; 9:1103-14. [PMID: 26588754 DOI: 10.1021/ct300913g] [Citation(s) in RCA: 114] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The microscopic behavior of water under different conditions and in different environments remains the subject of intense debate. A great number of the controversies arise due to the contradictory predictions obtained within different theoretical models. Relative to conclusions derived from force fields or density functional theory, there is comparably less room to dispute highly correlated electronic structure calculations. Unfortunately, such ab initio calculations are severely limited by system size. In this study, a detailed analysis of the two- and three-body water interactions evaluated at the CCSD(T) level is carried out to quantitatively assess the accuracy of several force fields, DFT models, and ab initio based interaction potentials that are commonly used in molecular simulations. On the basis of this analysis, a new model, HBB2-pol, is introduced which is capable of accurately mapping CCSD(T) results for water dimers and trimers into an efficient analytical function. The accuracy of HBB2-pol is further established through comparison with the experimentally determined second and third virial coefficients.
Collapse
Affiliation(s)
- Gregory R Medders
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92103, United States
| | - Volodymyr Babin
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92103, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92103, United States
| |
Collapse
|
42
|
Akin-Ojo O, Szalewicz K. How well can polarization models of pairwise nonadditive forces describe liquid water? J Chem Phys 2013; 138:024316. [DOI: 10.1063/1.4773821] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Omololu Akin-Ojo
- Theoretical Physics Department, African University of Science and Technology, Km 10 Airport Road, Galadimawa, Abuja, Nigeria
| | | |
Collapse
|
43
|
Schnell M, Erlekam U, Bunker PR, von Helden G, Grabow JU, Meijer G, van der Avoird A. Unraveling the internal dynamics of the benzene dimer: a combined theoretical and microwave spectroscopy study. Phys Chem Chem Phys 2013; 15:10207-23. [DOI: 10.1039/c3cp51181b] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
44
|
Leforestier C, Szalewicz K, van der Avoird A. Spectra of water dimer from a new ab initio potential with flexible monomers. J Chem Phys 2012; 137:014305. [PMID: 22779646 DOI: 10.1063/1.4722338] [Citation(s) in RCA: 92] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
We report the definition and testing of a new ab initio 12-dimensional potential for the water dimer with flexible monomers. Using our recent accurate CCpol-8s rigid water pair potential [W. Cencek, K. Szalewicz, C. Leforestier, R. van Harrevelt, and A. van der Avoird, Phys. Chem. Chem. Phys. 10, 4716 (2008)] as a reference for the undistorted monomers' geometries, a distortion correction has been added, which was taken from a former flexible-monomer ab initio potential. This correction allows us to retrieve the correct binding energy D(e)=21.0 kJ mol(-1), and leads to an equilibrium geometry in close agreement with the one obtained from benchmark calculations. The kinetic energy operator describing the flexible-monomer water dimer has been expressed in terms of Radau coordinates for each monomer and a recent general cluster polyspherical formulation describing their relative motions. Within this formulation, an adiabatic scheme has been invoked in order to decouple fast (intramolecular) modes and slow (intermolecular) ones. Different levels of approximation were tested, which differ in the way in which the residual potential coupling between the intramolecular modes located on different monomers and the dependence of the monomer rotational constants on the dimer geometry are handled. Accurate calculations of the vibration-rotation-tunneling levels of (H(2)O)(2) and (D(2)O)(2) were performed, which show the best agreement with experiments achieved so far for any water potential. Intramolecular excitations of the two monomers were calculated within two limiting cases, to account for the lack of non-adiabatic coupling between intramolecular modes due to the intermolecular motion. In the first model, the excitation was assumed to stay either on the donor or the acceptor molecule, and to hop between the two moieties upon donor-acceptor interchange. In the second model, the excitation remains on the same molecule whatever is the dimer geometry. Marginal frequency differences, less than 2 cm(-1), were obtained for all modes, and the resulting infrared shifts are in good agreement with experiments.
Collapse
Affiliation(s)
- Claude Leforestier
- Institut Charles Gerhardt (CTMM)-UMR 5253, CC 1501, Université Montpellier II-CNRS, 34095 Montpellier Cedex 05, France.
| | | | | |
Collapse
|
45
|
Ziemkiewicz MP, Pluetzer C, Nesbitt DJ, Scribano Y, Faure A, van der Avoird A. Overtone vibrational spectroscopy in H2-H2O complexes: A combined high level theoretical ab initio, dynamical and experimental study. J Chem Phys 2012; 137:084301. [DOI: 10.1063/1.4732581] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Wang J, Cieplak P, Cai Q, Hsieh MJ, Wang J, Duan Y, Luo R. Development of polarizable models for molecular mechanical calculations. 3. Polarizable water models conforming to Thole polarization screening schemes. J Phys Chem B 2012; 116:7999-8008. [PMID: 22712654 PMCID: PMC3493849 DOI: 10.1021/jp212117d] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
As an integrated step toward a coherent polarizable force field for biomolecular modeling, we analyzed four polarizable water models to evaluate their consistencies with the Thole polarization screening schemes utilized in our latest Amber polarizable force field. Specifically, we studied the performance of both the Thole linear and exponential schemes in these water models to assess their abilities to reproduce experimental water properties. The analysis shows that the tested water models reproduce most of the room-temperature properties of liquid water reasonably well but fall short of reproducing the dynamic properties and temperature-dependent properties. This study demonstrates the necessity to further fine-tune water polarizable potentials for more robust polarizable force fields for biomolecular simulations.
Collapse
Affiliation(s)
- Jun Wang
- Molecular Biology and Biochemistry, University of California-Irvine, Irvine, California 92697, United States
| | | | | | | | | | | | | |
Collapse
|
47
|
Tennyson J, Barber MJ, Kelly REA. An adiabatic model for calculating overtone spectra of dimers such as (H(2)O)(2). PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2012; 370:2656-2674. [PMID: 22547237 DOI: 10.1098/rsta.2011.0243] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
The near-infrared and visible wavelength spectrum of the water dimer is considered to be the major contributor to the so-called water continuum at these wavelengths. However, theoretical models of this spectrum require the simultaneous treatment of both monomer and dimer excitations. A model for treating this problem is proposed which is based upon a Franck-Condon-like separation between the monomer and dimer vibrational motions. In this model, one of the monomers is treated as the chromophore and its absorption is assumed to be given by its, possibly perturbed, vibrational band intensity. The main computational issue is the treatment of separate monomer and dimer motions. Various approaches for obtaining dimer vibration-rotation tunnelling spectra that allow for monomer motion are explored. These approaches include ways of treating the adiabatic separation of dimer vibrational modes from monomer vibrational modes. We classify the adiabatic separation methods under four main approaches: namely fixed-geometry, free-monomer, perturbed-monomer and coupled-monomer methods. The latter being the most computationally expensive as the monomer wave functions are dependent on the dimer coordinates. For each of these approaches, expectation values over the full potential are calculated for the given monomer vibrational wave functions. Various full (named VAP 2pD in the text) and partial (VAP (+p)D) averaging techniques are outlined to calculate the vibrationally averaged, monomer state-dependent, dimer interaction potentials. The computational costs associated with application of these techniques to the water dimer are estimated and the prospects for full calculations based on this approach are assessed.
Collapse
Affiliation(s)
- J Tennyson
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.
| | | | | |
Collapse
|
48
|
|
49
|
Musso F, Mignon P, Ugliengo P, Sodupe M. Cooperative effects at water–crystalline silica interfaces strengthen surface silanol hydrogen bonding. An ab initio molecular dynamics study. Phys Chem Chem Phys 2012; 14:10507-14. [DOI: 10.1039/c2cp40756f] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
50
|
Park K, Lin W, Paesani F. A Refined MS-EVB Model for Proton Transport in Aqueous Environments. J Phys Chem B 2011; 116:343-52. [DOI: 10.1021/jp208946p] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Kyoyeon Park
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| | - Wei Lin
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| |
Collapse
|