1
|
Förster M, Ukoji N, Sahle CJ, Niskanen J, Sakrowski R, Surmeier G, Weis C, Irifune T, Imoto S, Yavas H, Huotari S, Marx D, Sternemann C, Tse JS. Generating interstitial water within the persisting tetrahedral H-bond network explains density increase upon compressing liquid water. Proc Natl Acad Sci U S A 2024; 121:e2403662121. [PMID: 39284048 PMCID: PMC11441526 DOI: 10.1073/pnas.2403662121] [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: 02/23/2024] [Accepted: 07/16/2024] [Indexed: 10/02/2024] Open
Abstract
Despite its ubiquitous nature, the atomic structure of water in its liquid state is still controversially debated. We use a combination of X-ray Raman scattering spectroscopy in conjunction with ab initio and path integral molecular dynamics simulations to study the local atomic and electronic structure of water under high pressure conditions. Systematically increasing fingerprints of non-hydrogen-bonded H[Formula: see text]O molecules in the first hydration shell are identified in the experimental and computational oxygen K-edge excitation spectra. This provides evidence for a compaction mechanism in terms of a continuous collapse of the second hydration shell with increasing pressure via generation of interstitial water within locally tetrahedral hydrogen-bonding environments.
Collapse
Affiliation(s)
- Mirko Förster
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Nnanna Ukoji
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SKS7N 5E2, Canada
| | | | - Johannes Niskanen
- Department of Physics and Astronomy, University of Turku, Turun yliopistoFI-20014, Finland
| | - Robin Sakrowski
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Göran Surmeier
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Christopher Weis
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Tetsuo Irifune
- Geodynamics Research Center, Ehime University, Matsuyama790, Japan
| | - Sho Imoto
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum44780, Germany
| | - Hasan Yavas
- Deutsches Elektronen-Synchrotron, Hamburg22607, Germany
- Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, CA94025
| | - Simo Huotari
- Department of Physics, University of Helsinki, HelsinkiFI-00014, Finland
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum44780, Germany
| | | | - John S. Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SKS7N 5E2, Canada
| |
Collapse
|
2
|
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
|
3
|
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
|
4
|
Eltareb A, Lopez GE, Giovambattista N. A continuum of amorphous ices between low-density and high-density amorphous ice. Commun Chem 2024; 7:36. [PMID: 38378859 PMCID: PMC10879119 DOI: 10.1038/s42004-024-01117-2] [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: 10/11/2023] [Accepted: 02/01/2024] [Indexed: 02/22/2024] Open
Abstract
Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY, 11210, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| | - Gustavo E Lopez
- Department of Chemistry, Lehman College of the City University of New York, Bronx, NY, 10468, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| | - Nicolas Giovambattista
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY, 11210, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| |
Collapse
|
5
|
Daidone I, Foffi R, Amadei A, Zanetti-Polzi L. A statistical mechanical model of supercooled water based on minimal clusters of correlated molecules. J Chem Phys 2023; 159:094502. [PMID: 37655770 DOI: 10.1063/5.0157505] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 08/14/2023] [Indexed: 09/02/2023] Open
Abstract
In this paper, we apply a theoretical model for fluid state thermodynamics to investigate simulated water in supercooled conditions. This model, which we recently proposed and applied to sub- and super-critical fluid water [Zanetti-Polzi et al., J. Chem. Phys. 156(4), 44506 (2022)], is based on a combination of the moment-generating functions of the enthalpy and volume fluctuations as provided by two gamma distributions and provides the free energy of the system as well as other relevant thermodynamic quantities. The application we make here provides a thermodynamic description of supercooled water fully consistent with that expected by crossing the liquid-liquid Widom line, indicating the presence of two distinct liquid states. In particular, the present model accurately reproduces the Widom line temperatures estimated with other two-state models and well describes the heat capacity anomalies. Differently from previous models, according to our description, a cluster of molecules that extends beyond the first hydration shell is necessary to discriminate between the statistical fluctuation regimes typical of the two liquid states.
Collapse
Affiliation(s)
- Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| | - Riccardo Foffi
- Department of Civil, Environmental and Geomatic Engineering, Institute of Environmental Engineering, ETH Zürich, Laura-Hezner-Weg 7, 8093 Zürich, Switzerland
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome "Tor Vergata," Via della Ricerca Scientifica, I-00185 Rome, Italy
| | - Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| |
Collapse
|
6
|
Dhabal D, Molinero V. Kinetics and Mechanisms of Pressure-Induced Ice Amorphization and Polyamorphic Transitions in a Machine-Learned Coarse-Grained Water Model. J Phys Chem B 2023; 127:2847-2862. [PMID: 36920450 DOI: 10.1021/acs.jpcb.3c00434] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Water glasses have attracted considerable attention due to their potential connection to a liquid-liquid transition in supercooled water. Here we use molecular simulations to investigate the formation and phase behavior of water glasses using the machine-learned bond-order parameter (ML-BOP) water model. We produce glasses through hyperquenching of water, pressure-induced amorphization (PIA) of ice, and pressure-induced polyamorphic transformations. We find that PIA of polycrystalline ice occurs at a lower pressure than that of monocrystalline ice and through a different mechanism. The temperature dependence of the amorphization pressure of polycrystalline ice for ML-BOP agrees with that in experiments. We also find that ML-BOP accurately reproduces the density, coordination number, and structural features of low-density (LDA), high-density (HDA), and very high-density (VHDA) amorphous water glasses. ML-BOP accurately reproduces the experimental radial distribution function of LDA but overpredicts the minimum between the first two shells in high-density glasses. We examine the kinetics and mechanism of the transformation between low-density and high-density glasses and find that the sharp nature of these transitions in ML-BOP is similar to that in experiments and all-atom water models with a liquid-liquid transition. Transitions between ML-BOP glasses occur through a spinodal-like mechanism, similar to ice crystallization from LDA. Both glass-to-glass and glass-to-ice transformations have Avrami-Kolmogorov kinetics with exponent n = 1.5 ± 0.2 in experiments and simulations. Importantly, ML-BOP reproduces the competition between crystallization and HDA→LDA transition above the glass transition temperature Tg, and separation of their time scales below Tg, observed also in experiments. These findings demonstrate the ability of ML-BOP to accurately reproduce water properties across various regimes, making it a promising model for addressing the competition between polyamorphic transitions and crystallization in water and solutions.
Collapse
Affiliation(s)
- Debdas Dhabal
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
Martelli F. Steady-like topology of the dynamical hydrogen bond network in supercooled water. PNAS NEXUS 2022; 1:pgac090. [PMID: 36741425 PMCID: PMC9896910 DOI: 10.1093/pnasnexus/pgac090] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 06/08/2022] [Indexed: 02/07/2023]
Abstract
We investigate the link between topology of the hydrogen bond network (HBN) and large-scale density fluctuations in water from ambient conditions to the glassy state. We observe a transition from a temperature-dependent topology at high temperatures, to a steady-like topology below the Widom temperature TW ∼ 220 K signaling the fragile-to-strong crossover and the maximum in structural fluctuations. As a consequence of the steady topology, the network suppresses large-scale density fluctuations much more efficiently than at higher temperatures. Below TW , the contribution of coordination defects of the kind A 2 D 1 (two acceptors and one donor) to the kinetics of the HBN becomes progressively more pronounced, suggesting that A 2 D 1 configurations may represent the main source of dynamical heterogeneities. Below the vitrification temperature, the freezing of rotational and translational degrees of freedom allow for an enhanced suppression of large-scale density fluctuations and the sample reaches the edges of nearly hyperuniformity. The formed network still hosts coordination defects, hence implying that nearly hyperuniformity goes beyond the classical continuous random network paradigm of tetrahedral networks and can emerge in scenarios much more complex than previously assumed. Our results unveil a hitherto undisclosed link between network topology and properties of water essential for better understanding water's rich and complex nature. Beyond implications for water, our findings pave the way to a better understanding of the physics of supercooled liquids and disordered hyperuniform networks at large.
Collapse
|
10
|
Faccio C, Benzi M, Zanetti-Polzi L, Daidone I. Low- and high-density forms of liquid water revealed by a new medium-range order descriptor. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118922] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|