Water Potential from Adaptive Force Matching for Ice and Liquid with Revised Dispersion Predicts Supercooled Liquid Anomalies in Good Agreement with Two Independent Experimental Fits

Overview

Journal J Phys Chem B

Publisher American Chemical Society

Specialty Chemistry

Date 2024 Mar 27

PMID 38536126

Authors

Raymond Weldon

Feng Wang

Affiliations

Soon will be listed here.

Abstract

A revised version of the Water potential from Adaptive force matching for Ice and Liquid (WAIL) was developed by using the previous data set for fitting the WAIL model but with a dispersion term calculated using symmetry adapted perturbation theory (SAPT). The model has no adjustable parameters and relies solely on fitting first-principles information. The new model, named revised WAIL (rWAIL), shows improved predictions of most properties of water when compared to the previously published WAIL model. The rWAIL model also compares favorably to other first-principles-derived water models, such as MB-Pol, at only a fraction of the computational cost. The rWAIL model is used to study the properties of supercooled water. The model shows evidence of a liquid-liquid phase transition (LLPT) in the supercooled regimes with the liquid-liquid critical point (LLCP) at 203 K and 90 MPa. This estimate is in good agreement with a recent polynomial fit to the experimental density of water. Also, the fit to the surface tension of supercooled water based on the rWAIL model shows excellent agreement with the corresponding fit to the experimental data. Consistent with previously published molecular dynamics and experimental data, the surface tension of water exhibits exponential growth in the supercooled regime, which is likely a result of the emergence of a low-density liquid form of water. The simulation thus unites two separate experimental fits with one first-principles-based model, lending strong evidence of an LLPT in real water.

References

Akin-Ojo O, Song Y, Wang F . Developing ab initio quality force fields from condensed phase quantum-mechanics/molecular-mechanics calculations through the adaptive force matching method. J Chem Phys. 2008; 129(6):064108. DOI: 10.1063/1.2965882. View

Gallo P, Amann-Winkel K, Angell C, Anisimov M, Caupin F, Chakravarty C . Water: A Tale of Two Liquids. Chem Rev. 2016; 116(13):7463-500. PMC: 5424717. DOI: 10.1021/acs.chemrev.5b00750. View

Habershon S, Markland T, Manolopoulos D . Competing quantum effects in the dynamics of a flexible water model. J Chem Phys. 2009; 131(2):024501. DOI: 10.1063/1.3167790. View

Muniz M, Gartner 3rd T, Riera M, Knight C, Yue S, Paesani F . Vapor-liquid equilibrium of water with the MB-pol many-body potential. J Chem Phys. 2021; 154(21):211103. DOI: 10.1063/5.0050068. View

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. 2015; 9(12):5395-403. DOI: 10.1021/ct400863t. View

Hoover . Canonical dynamics: Equilibrium phase-space distributions. Phys Rev A Gen Phys. 1985; 31(3):1695-1697. DOI: 10.1103/physreva.31.1695. View

Weldon R, Wang F . Simulating a flexible water model as rigid: Best practices and lessons learned. J Chem Phys. 2023; 158(13):134506. PMC: 10076064. DOI: 10.1063/5.0143836. View

Vega C, Sanz E, Abascal J . The melting temperature of the most common models of water. J Chem Phys. 2005; 122(11):114507. DOI: 10.1063/1.1862245. View

Grimme S, Antony J, Ehrlich S, Krieg H . A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J Chem Phys. 2010; 132(15):154104. DOI: 10.1063/1.3382344. View

10.

Mishima O, Sumita T . Equation of State of Liquid Water Written by Simple Experimental Polynomials and the Liquid-Liquid Critical Point. J Phys Chem B. 2023; 127(6):1414-1421. DOI: 10.1021/acs.jpcb.2c08342. View

11.

Li J, Wang F . Accurate Prediction of the Hydration Free Energies of 20 Salts through Adaptive Force Matching and the Proper Comparison with Experimental References. J Phys Chem B. 2017; 121(27):6637-6645. PMC: 5665695. DOI: 10.1021/acs.jpcb.7b04618. View

12.

Bussi G, Donadio D, Parrinello M . Canonical sampling through velocity rescaling. J Chem Phys. 2007; 126(1):014101. DOI: 10.1063/1.2408420. View

13.

Foffi R, Russo J, Sciortino F . Structural and topological changes across the liquid-liquid transition in water. J Chem Phys. 2021; 154(18):184506. DOI: 10.1063/5.0049299. View

14.

Rogers T, Leong K, Wang F . Possible Evidence for a New Form of Liquid Buried in the Surface Tension of Supercooled Water. Sci Rep. 2016; 6:33284. PMC: 5018856. DOI: 10.1038/srep33284. View

15.

Song Y, Akin-Ojo O, Wang F . Correcting for dispersion interaction and beyond in density functional theory through force matching. J Chem Phys. 2010; 133(17):174115. DOI: 10.1063/1.3503656. View

16.

Akin-Ojo O, Wang F . The quest for the best nonpolarizable water model from the adaptive force matching method. J Comput Chem. 2010; 32(3):453-62. DOI: 10.1002/jcc.21634. View

17.

Abascal J, Sanz E, Garcia Fernandez R, Vega C . A potential model for the study of ices and amorphous water: TIP4P/Ice. J Chem Phys. 2005; 122(23):234511. DOI: 10.1063/1.1931662. View

18.

Bore S, Paesani F . Realistic phase diagram of water from "first principles" data-driven quantum simulations. Nat Commun. 2023; 14(1):3349. PMC: 10250386. DOI: 10.1038/s41467-023-38855-1. View

19.

Gartner 3rd T, Zhang L, Piaggi P, Car R, Panagiotopoulos A, Debenedetti P . Signatures of a liquid-liquid transition in an ab initio deep neural network model for water. Proc Natl Acad Sci U S A. 2020; 117(42):26040-26046. PMC: 7584908. DOI: 10.1073/pnas.2015440117. View

20.

Hu H, Wang F . The liquid-vapor equilibria of TIP4P/2005 and BLYPSP-4F water models determined through direct simulations of the liquid-vapor interface. J Chem Phys. 2015; 142(21):214507. PMC: 4457660. DOI: 10.1063/1.4922166. View