1
|
Yang C, Ladd-Parada M, Nam K, Jeong S, You S, Eklund T, Späh A, Pathak H, Lee JH, Eom I, Kim M, Perakis F, Nilsson A, Kim KH, Amann-Winkel K. Unveiling a common phase transition pathway of high-density amorphous ices through time-resolved x-ray scattering. J Chem Phys 2024; 160:244503. [PMID: 38916268 DOI: 10.1063/5.0216904] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2024] [Accepted: 06/07/2024] [Indexed: 06/26/2024] Open
Abstract
Here, we investigate the hypothesis that despite the existence of at least two high-density amorphous ices, only one high-density liquid state exists in water. We prepared a very-high-density amorphous ice (VHDA) sample and rapidly increased its temperature to around 205 ± 10 K using laser-induced isochoric heating. This temperature falls within the so-called "no-man's land" well above the glass-liquid transition, wherein the IR laser pulse creates a metastable liquid state. Subsequently, this high-density liquid (HDL) state of water decompresses over time, and we examined the time-dependent structural changes using short x-ray pulses from a free electron laser. We observed a liquid-liquid transition to low-density liquid water (LDL) over time scales ranging from 20 ns to 3 μs, consistent with previous experimental results using expanded high-density amorphous ice (eHDA) as the initial state. In addition, the resulting LDL derived both from VHDA and eHDA displays similar density and degree of inhomogeneity. Our observation supports the idea that regardless of the initial annealing states of the high-density amorphous ices, the same HDL and final LDL states are reached at temperatures around 205 K.
Collapse
Affiliation(s)
- Cheolhee Yang
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk 37673, Republic of Korea
| | - Marjorie Ladd-Parada
- Chemistry Department, Glycoscience Division, Kungliga Tekniska Högskola, Roslagstullsbacken 21, 11421 Stockholm, Sweden
| | - Kyeongmin Nam
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk 37673, Republic of Korea
| | - Sangmin Jeong
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk 37673, Republic of Korea
| | - Seonju You
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk 37673, Republic of Korea
| | - Tobias Eklund
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
- Max-Planck-Institute for Polymer Research, 55128 Mainz, Germany
- Institute for Physics, Johannes Gutenberg University Mainz, 55128 Mainz, Germany
| | - Alexander Späh
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Harshad Pathak
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Jae Hyuk Lee
- Pohang Accelerator Laboratory, POSTECH, Pohang, Gyeongbuk 37673, Republic of Korea
| | - Intae Eom
- Pohang Accelerator Laboratory, POSTECH, Pohang, Gyeongbuk 37673, Republic of Korea
| | - Minseok Kim
- Pohang Accelerator Laboratory, POSTECH, Pohang, Gyeongbuk 37673, Republic of Korea
| | - Fivos Perakis
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Anders Nilsson
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Kyung Hwan Kim
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk 37673, Republic of Korea
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
- Max-Planck-Institute for Polymer Research, 55128 Mainz, Germany
- Institute for Physics, Johannes Gutenberg University Mainz, 55128 Mainz, Germany
| |
Collapse
|
2
|
Dhabal D, Kumar R, Molinero V. Liquid-liquid transition and ice crystallization in a machine-learned coarse-grained water model. Proc Natl Acad Sci U S A 2024; 121:e2322853121. [PMID: 38709921 PMCID: PMC11098087 DOI: 10.1073/pnas.2322853121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2023] [Accepted: 03/27/2024] [Indexed: 05/08/2024] Open
Abstract
Mounting experimental evidence supports the existence of a liquid-liquid transition (LLT) in high-pressure supercooled water. However, fast crystallization of supercooled water has impeded identification of the LLT line TLL(p) in experiments. While the most accurate all-atom (AA) water models display a LLT, their computational cost limits investigations of its interplay with ice formation. Coarse-grained (CG) models provide over 100-fold computational efficiency gain over AA models, enabling the study of water crystallization, but have not yet shown to have a LLT. Here, we demonstrate that the CG machine-learned water model Machine-Learned Bond-Order Potential (ML-BOP) has a LLT that ends in a critical point at pc = 170 ± 10 MPa and Tc = 181 ± 3 K. The TLL(p) of ML-BOP is almost identical to the one of TIP4P/2005, adding to the similarity in the equation of state of liquid water in both models. Cooling simulations reveal that ice crystallization is fastest at the LLT and its supercritical continuation of maximum heat capacity, supporting a mechanistic relationship between the structural transformation of water to a low-density liquid (LDL) and ice formation. We find no signature of liquid-liquid criticality in the ice crystallization temperatures. ML-BOP replicates the competition between formation of LDL and ice observed in ultrafast experiments of decompression of the high-density liquid (HDL) into the region of stability of LDL. The simulations reveal that crystallization occurs prior to the coarsening of the HDL and LDL domains, obscuring the distinction between the highly metastable first-order LLT and pronounced structural fluctuations along its supercritical continuation.
Collapse
Affiliation(s)
- Debdas Dhabal
- Department of Chemistry, The University of Utah, Salt Lake City, UT84112-0850
| | - Rajat Kumar
- Department of Chemistry, The University of Utah, Salt Lake City, UT84112-0850
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, UT84112-0850
| |
Collapse
|
3
|
Eltareb A, Lopez GE, Giovambattista N. Potential energy landscape of a flexible water model: Equation of state, configurational entropy, and Adam-Gibbs relationship. J Chem Phys 2024; 160:154510. [PMID: 38639318 PMCID: PMC11184974 DOI: 10.1063/5.0200306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Accepted: 03/28/2024] [Indexed: 04/20/2024] Open
Abstract
The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid-liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam-Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully).
Collapse
Affiliation(s)
- Ali Eltareb
- Authors to whom correspondence should be addressed: ; ; and
| | | | | |
Collapse
|
4
|
Weldon R, Wang F. 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. J Phys Chem B 2024; 128:3398-3407. [PMID: 38536126 PMCID: PMC11017247 DOI: 10.1021/acs.jpcb.3c06495] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2023] [Revised: 02/15/2024] [Accepted: 03/11/2024] [Indexed: 04/12/2024]
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.
Collapse
Affiliation(s)
- Raymond Weldon
- Department of Chemistry and
Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Feng Wang
- Department of Chemistry and
Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| |
Collapse
|
5
|
Sciortino F, Gartner TE, Debenedetti PG. Free-energy landscape and spinodals for the liquid-liquid transition of the TIP4P/2005 and TIP4P/Ice models of water. J Chem Phys 2024; 160:104501. [PMID: 38456528 DOI: 10.1063/5.0196964] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 02/20/2024] [Indexed: 03/09/2024] Open
Abstract
Continued increases in computational power now make it possible to evaluate the free-energy landscape associated with the first-order liquid-liquid transition in realistic models of water for which an accurate estimate of the liquid-liquid critical point exists, and to explore its change with pressure near the coexistence line. We report the results of 50 μs-long NPT umbrella sampling simulations for two realistic models for water, TIP4P/2005 and TIP4P/ice, 3-9 K below their critical temperatures. The free energy profile at different pressures clearly shows the presence of two well-defined free energy basins and makes it possible to identify the liquid-liquid spinodal points, the limits of stability that define the (temperature dependent) pressure range within which two distinct free energy basins exist. The results show that for temperatures less than 10 K below the critical temperature, metastable states are possible across a very limited pressure interval, information that is relevant to the interpretation of experiments probing the metastable phase behavior of deeply supercooled water in the so-called no-man's land.
Collapse
Affiliation(s)
- Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Rome, Italy
| | - Thomas E Gartner
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, USA
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
6
|
Maturi F, Raposo Filho RS, Brites CDS, Fan J, He R, Zhuang B, Liu X, Carlos LD. Deciphering Density Fluctuations in the Hydration Water of Brownian Nanoparticles via Upconversion Thermometry. J Phys Chem Lett 2024; 15:2606-2615. [PMID: 38420927 PMCID: PMC10926164 DOI: 10.1021/acs.jpclett.4c00044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 02/14/2024] [Accepted: 02/16/2024] [Indexed: 03/02/2024]
Abstract
We investigate the intricate relationship among temperature, pH, and Brownian velocity in a range of differently sized upconversion nanoparticles (UCNPs) dispersed in water. These UCNPs, acting as nanorulers, offer insights into assessing the relative proportion of high-density and low-density liquid in the surrounding hydration water. The study reveals a size-dependent reduction in the onset temperature of liquid-water fluctuations, indicating an augmented presence of high-density liquid domains at the nanoparticle surfaces. The observed upper-temperature threshold is consistent with a hypothetical phase diagram of water, validating the two-state model. Moreover, an increase in pH disrupts the organization of water molecules, similar to external pressure effects, allowing simulation of the effects of temperature and pressure on hydrogen bonding networks. The findings underscore the significance of the surface of suspended nanoparticles for understanding high- to low-density liquid fluctuations and water behavior at charged interfaces.
Collapse
Affiliation(s)
- Fernando
E. Maturi
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
- Institute
of Chemistry, São Paulo State University
(UNESP), 14800-060 Araraquara, SP, Brazil
| | - Ramon S. Raposo Filho
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Carlos D. S. Brites
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Jingyue Fan
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Ruihua He
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Bilin Zhuang
- Harvey
Mudd College, 301 Platt
Boulevard, Claremont, California 91711, United States
| | - Xiaogang Liu
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Luís D. Carlos
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| |
Collapse
|
7
|
Foffi R, Sciortino F. Identification of local structures in water from supercooled to ambient conditions. J Chem Phys 2024; 160:094504. [PMID: 38436442 DOI: 10.1063/5.0188764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 02/09/2024] [Indexed: 03/05/2024] Open
Abstract
Studies of water thermodynamics have long been tied to the identification of two distinct families of local structures, whose competition could explain the origin of the many thermodynamic anomalies and the hypothesized liquid-liquid critical point in water. Despite the many successes and insights gained, the structural indicators proposed throughout the years were not able to unequivocally identify these two families over a wide range of conditions. We show that a recently introduced indicator, Ψ, which exploits information on the hydrogen bond network connectivity, can reliably identify these two distinct local environments over a wide range of thermodynamic conditions (188-300 K and 0-13 kbar) and that close to the liquid-liquid critical point, the spatial correlations of density fluctuations are identical to those of the Ψ indicator. Our results strongly support the idea that water thermodynamic properties arise from the competition between two distinct and identifiable local environments.
Collapse
Affiliation(s)
- Riccardo Foffi
- Department of Civil, Environmental and Geomatic Engineering, Institute for Environmental Engineering, ETH Zürich, Laura-Hezner-Weg 7, 8093 Zürich, Switzerland
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Rome, Italy
| |
Collapse
|
8
|
Palmer JC, Sarupria S, Truskett TM. Tribute to Pablo G. Debenedetti. J Phys Chem B 2023; 127:8075-8078. [PMID: 37766640 DOI: 10.1021/acs.jpcb.3c06020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/29/2023]
Affiliation(s)
- Jeremy C Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, United States
| | - Sapna Sarupria
- Department of Chemistry, Chemical Theory Center, University of Minnesota Twin Cities, Minneapolis, Minnesota 55455, United States
| | - Thomas M Truskett
- McKetta Department of Chemical Engineering and Department of Physics, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
9
|
Piaggi PM, Gartner TE, Car R, Debenedetti PG. Melting curves of ice polymorphs in the vicinity of the liquid-liquid critical point. J Chem Phys 2023; 159:054502. [PMID: 37531247 DOI: 10.1063/5.0159288] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Accepted: 07/14/2023] [Indexed: 08/04/2023] Open
Abstract
The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164-168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. We explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.
Collapse
Affiliation(s)
- Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Thomas E Gartner
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
- Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
10
|
Guidarelli Mattioli F, Sciortino F, Russo J. Are Neural Network Potentials Trained on Liquid States Transferable to Crystal Nucleation? A Test on Ice Nucleation in the mW Water Model. J Phys Chem B 2023; 127:3894-3901. [PMID: 37075256 PMCID: PMC10165654 DOI: 10.1021/acs.jpcb.3c00693] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 04/06/2023] [Indexed: 04/21/2023]
Abstract
Neural network potentials (NNPs) are increasingly being used to study processes that happen on long time scales. A typical example is crystal nucleation, which rate is controlled by the occurrence of a rare fluctuation, i.e., the appearance of the critical nucleus. Because the properties of this nucleus are far from those of the bulk crystal, it is yet unclear whether NN potentials trained on equilibrium liquid states can accurately describe nucleation processes. So far, nucleation studies on NNPs have been limited to ab initio models whose nucleation properties are unknown, preventing an accurate comparison. Here we train a NN potential on the mW model of water─a classical three-body potential whose nucleation time scale is accessible in standard simulations. We show that a NNP trained only on a small number of liquid state points can reproduce with great accuracy the nucleation rates and free energy barriers of the original model, computed from both spontaneous and biased trajectories, strongly supporting the use of NNPs to study nucleation events.
Collapse
Affiliation(s)
| | | | - John Russo
- Sapienza University of Rome, Piazzale Aldo Moro 2, 00185 Rome, Italy
| |
Collapse
|
11
|
Cerdeiriña CA, González-Salgado D, Troncoso J. Liquid-Liquid Criticality in TIP4P/2005 and Three-State Models of Water. J Phys Chem B 2023; 127:3902-3910. [PMID: 37097210 PMCID: PMC10165646 DOI: 10.1021/acs.jpcb.3c00696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Molecular dynamics simulations leading to the isothermal compressibility, the isobaric thermal expansivity, and the isobaric heat capacity of TIP4P/2005 water are found to be consistent with the coordinates of its second, liquid-liquid critical point reported recently by Debenedetti et al. [ Science 2020, 369, 289-292]. In accord with the theory of critical phenomena, we encounter that the rise in the magnitude of these response functions as temperature is lowered is especially marked along the critical isochore. Furthermore, response-function ratios provide a test for thermodynamic consistency at the critical point and manifest nonuniversal features sharply distinguishing liquid-liquid from standard gas-liquid criticality. The whole pattern of behavior revealed by simulations is qualitatively the same as the one of a three-state Ising model of water exhibiting a low-temperature liquid-liquid critical point. Exact solutions for the two-state components of such a three-state model are also provided.
Collapse
Affiliation(s)
- Claudio A Cerdeiriña
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQFR, Ourense 32004, Spain
| | - Diego González-Salgado
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQFR, Ourense 32004, Spain
| | - Jacobo Troncoso
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQFR, Ourense 32004, Spain
| |
Collapse
|
12
|
Guidarelli Mattioli F, Sciortino F, Russo J. A neural network potential with self-trained atomic fingerprints: A test with the mW water potential. J Chem Phys 2023; 158:104501. [PMID: 36922151 DOI: 10.1063/5.0139245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023] Open
Abstract
We present a neural network (NN) potential based on a new set of atomic fingerprints built upon two- and three-body contributions that probe distances and local orientational order, respectively. Compared with the existing NN potentials, the atomic fingerprints depend on a small set of tunable parameters that are trained together with the NN weights. In addition to simplifying the selection of the atomic fingerprints, this strategy can also considerably increase the overall accuracy of the network representation. To tackle the simultaneous training of the atomic fingerprint parameters and NN weights, we adopt an annealing protocol that progressively cycles the learning rate, significantly improving the accuracy of the NN potential. We test the performance of the network potential against the mW model of water, which is a classical three-body potential that well captures the anomalies of the liquid phase. Trained on just three state points, the NN potential is able to reproduce the mW model in a very wide range of densities and temperatures, from negative pressures to several GPa, capturing the transition from an open random tetrahedral network to a dense interpenetrated network. The NN potential also reproduces very well properties for which it was not explicitly trained, such as dynamical properties and the structure of the stable crystalline phases of mW.
Collapse
Affiliation(s)
| | | | - John Russo
- Sapienza University of Rome, Piazzale Aldo Moro 2, 00185 Rome, Italy
| |
Collapse
|
13
|
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:1414-1421. [PMID: 36744620 DOI: 10.1021/acs.jpcb.2c08342] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Experimental and theoretical evidence has been accumulating to support the liquid-liquid critical point (LLCP) hypothesis for water. However, no agreement has yet been reached on the pressure and temperature of LLCP. Here we made simple experimental equations that reproduced the volume of liquid water measured over a wide pressure-temperature range. They were polynomials that calculate pressure using volume and temperature as variables, and coefficients were determined by the method of least-squares. We analyzed the polynomials by changing the volume data and the number of terms in the polynomials and extrapolated them slightly to low temperatures. Consequently, the available experimental volume of liquid water indicated (but did not prove) the existence of LLCP at low temperature. Representative polynomials suggested that LLCP locates around 105 ± ∼ 9 MPa, 207 ± ∼ 5 K, and 0.993 ± ∼ 0.009 cm3/g.
Collapse
Affiliation(s)
- Osamu Mishima
- Research Center for Advanced Measurement and Characterization, National Institute for Materials Science (NIMS), Tsukuba, Ibaraki305-0044, Japan
| | - Tatsuya Sumita
- Geological Survey of Japan, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki305-8567, Japan
| |
Collapse
|
14
|
Amann-Winkel K, Kim KH, Giovambattista N, Ladd-Parada M, Späh A, Perakis F, Pathak H, Yang C, Eklund T, Lane TJ, You S, Jeong S, Lee JH, Eom I, Kim M, Park J, Chun SH, Poole PH, Nilsson A. Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice. Nat Commun 2023; 14:442. [PMID: 36707522 PMCID: PMC9883474 DOI: 10.1038/s41467-023-36091-1] [Citation(s) in RCA: 16] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Accepted: 01/13/2023] [Indexed: 01/28/2023] Open
Abstract
Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale < 100 ns and widely separated from ice formation, which begins at times >1 μs.
Collapse
Affiliation(s)
- Katrin Amann-Winkel
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden ,grid.419547.a0000 0001 1010 1663Max Planck Institute for Polymer Research and Johannes Gutenberg University, 55128 Mainz, Germany
| | - Kyung Hwan Kim
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Nicolas Giovambattista
- grid.183006.c0000 0001 0671 7844Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY 11210 USA ,grid.253482.a0000 0001 0170 7903The Graduate Center of the City University of New York, New York, NY 10016 USA
| | - Marjorie Ladd-Parada
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden ,grid.411313.50000 0004 0512 3288Present Address: Division of Glycoscience, School of Biotechnology, KTH Royal Institute of Technology, AlbaNova University Center, SE-10691 Stockholm, Sweden
| | - Alexander Späh
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Fivos Perakis
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Harshad Pathak
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Cheolhee Yang
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Tobias Eklund
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Thomas J. Lane
- grid.445003.60000 0001 0725 7771SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025 USA
| | - Seonju You
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Sangmin Jeong
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Jae Hyuk Lee
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Intae Eom
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Minseok Kim
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Jaeku Park
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Sae Hwan Chun
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Peter H. Poole
- grid.264060.60000 0004 1936 7363Department of Physics, St. Francis Xavier University, Antigonish, NS B2G2W5 Canada
| | - Anders Nilsson
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
15
|
Foffi R, Sciortino F. Correlated Fluctuations of Structural Indicators Close to the Liquid-Liquid Transition in Supercooled Water. J Phys Chem B 2022; 127:378-386. [PMID: 36538764 PMCID: PMC9841516 DOI: 10.1021/acs.jpcb.2c07169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Multiple numerical studies have unambiguously shown the existence of a liquid-liquid critical point in supercooled states for different numerical models of water, and various structural indicators have been put forward to describe the transformation associated with this phase transition. Here we analyze numerical simulations of near-critical supercooled water to compare the behavior of several of such indicators with critical density fluctuations. We show that close to the critical point most indicators are strongly correlated to density, and some of them even display identical distributions of fluctuations. These indicators probe the exact same free energy landscape, therefore providing a thermodynamic description of critical supercooled water which is identical to that provided by the density order parameter. This implies that close to the critical point, there is a tight coupling between many, only apparently distinct, structural degrees of freedom.
Collapse
Affiliation(s)
- Riccardo Foffi
- Institute
for Environmental Engineering, Department of Civil, Environmental
and Geomatic Engineering, ETH Zürich, 8093Zürich, Switzerland
| | - Francesco Sciortino
- Dipartimento
di Fisica, Sapienza Università di
Roma, Piazzale Aldo Moro
5, I-00185Rome, Italy,E-mail:
| |
Collapse
|
16
|
Gartner TE, Piaggi PM, Car R, Panagiotopoulos AZ, Debenedetti PG. Liquid-Liquid Transition in Water from First Principles. PHYSICAL REVIEW LETTERS 2022; 129:255702. [PMID: 36608224 DOI: 10.1103/physrevlett.129.255702] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Accepted: 11/14/2022] [Indexed: 06/17/2023]
Abstract
A long-standing question in water research is the possibility that supercooled liquid water can undergo a liquid-liquid phase transition (LLT) into high- and low-density liquids. We used several complementary molecular simulation techniques to evaluate the possibility of an LLT in an ab initio neural network model of water trained on density functional theory calculations with the SCAN exchange correlation functional. We conclusively show the existence of a first-order LLT and an associated critical point in the SCAN description of water, representing the first definitive computational evidence for an LLT in water from first principles.
Collapse
Affiliation(s)
- Thomas E Gartner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
- Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
- Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
17
|
Mondal A, Ramesh G, Singh RS. Manifestations of the structural origin of supercooled water’s anomalies in the heterogeneous relaxation on the potential energy landscape. J Chem Phys 2022; 157:184503. [DOI: 10.1063/5.0124041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
Liquid water is well-known for its intriguing thermodynamic anomalies in the supercooled state. The phenomenological two-state models—based on the assumption of the existence of two types of competing local states (or, structures) in liquid water—have been extremely successful in describing water’s thermodynamic anomalies. However, the precise structural features of these competing local states in liquid water still remain elusive. Here, we have employed a predefined structural order parameter-free approach to unambiguously identify two types of competing local states—entropically and energetically favored—with significantly different structural and energetic features in the TIP4P/2005 liquid water. This identification is based on the heterogeneous structural relaxation of the system in the potential energy landscape (PEL) during the steepest-descent energy minimization. This heterogeneous relaxation is characterized using order parameters inspired by the spin-glass transition in frustrated magnetic systems. We have further established a direct relationship between the population fluctuation of the two states and the anomalous behavior of the heat capacity in supercooled water. The composition-dependent spatial distribution of the entropically favored local states shows an interesting crossover from a spanning network-like single cluster to the spatially delocalized clusters in the close vicinity of the Widom line. Additionally, this study establishes a direct relationship between the topographic features of the PEL and the water’s thermodynamic anomalies in the supercooled state and provides alternate markers (in addition to the locus of maxima of thermodynamic response functions) for the Widom line in the phase plane.
Collapse
Affiliation(s)
- Arijit Mondal
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| | - Gadha Ramesh
- Department of Physics, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| | - Rakesh S. Singh
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| |
Collapse
|
18
|
Cerdeiriña CA. Water's Unusual Thermodynamics in the Realm of Physical Chemistry. J Phys Chem B 2022; 126:6608-6613. [PMID: 36001372 PMCID: PMC9797112 DOI: 10.1021/acs.jpcb.2c05274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
While it is known since the early work by Edsall, Frank and Evans, Kauzmann, and others that the thermodynamics of solvation of nonpolar solutes in water is unusual and has implications for the thermodynamics of protein folding, only recently have its connections with the unusual temperature dependence of the density of solvent water been illuminated. Such density behavior is, in turn, one of the manifestations of a nonstandard thermodynamic pattern contemplating a second, liquid-liquid critical point at conditions of temperature and pressure at which water exists as a deeply supercooled liquid. Recent experimental and computational work unambiguously points toward the existence of such a critical point, thereby providing concrete answers to the questions posed by the 1976 pioneering experiments by Speedy and Angell and the associated "liquid-liquid transition hypothesis" posited in 1992 by Stanley and co-workers. Challenges of this phenomenology to the branch of Statistical Mechanics remain.
Collapse
|