1
|
Aasen A, Hammer M, Reguera D, Wilhelmsen Ø. Estimating metastable thermodynamic properties by isochoric extrapolation from stable states. J Chem Phys 2024; 161:044113. [PMID: 39051829 DOI: 10.1063/5.0220207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2024] [Accepted: 07/01/2024] [Indexed: 07/27/2024] Open
Abstract
The description of metastable fluids, those in local but not global equilibrium, remains an important problem of thermodynamics, and it is crucial for many industrial applications and all first order phase transitions. One way to estimate their properties is by extrapolation from nearby stable states. This is often done isothermally, in terms of a virial expansion for gases or a Taylor expansion in density for liquids. This work presents evidence that an isochoric expansion of pressure at a given temperature is superior to an isothermal density expansion. Two different isochoric extrapolation strategies are evaluated, one best suited for vapors and one for liquids. Both are exact for important model systems, including the van der Waals equation of state. Moreover, we present a simple method to evaluate all the coefficients of the isochoric expansion directly from a simulation in the canonical ensemble. Using only the properties of stable states, the isochoric extrapolation methods reproduce simulation results with Lennard-Jones potentials, mostly within their uncertainties. The isochoric extrapolation methods are able to predict deeply metastable pressures accurately even from temperatures well above the critical. Isochoric extrapolation also predicts a mechanical stability limit, i.e., the thermodynamic spinodal. For water, the liquid spinodal pressure is predicted to be monotonically decreasing with decreasing temperature, in contrast to the re-entrant behavior predicted by the direct extension of the reference equation of state.
Collapse
Affiliation(s)
- Ailo Aasen
- Department of Gas Technology, SINTEF Energy Research, NO-7465 Trondheim, Norway
| | - Morten Hammer
- Department of Gas Technology, SINTEF Energy Research, NO-7465 Trondheim, Norway
- Porelab, Department of Chemistry, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| | - David Reguera
- Departament de Física de la Matèria Condensada, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
- Universitat de Barcelona Institute of Complex Systems (UBICS), Martí i Franquès 1, 08028 Barcelona, Spain
| | - Øivind Wilhelmsen
- Department of Gas Technology, SINTEF Energy Research, NO-7465 Trondheim, Norway
- Porelab, Department of Chemistry, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| |
Collapse
|
2
|
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
|
3
|
Guo C, Yang K, Qin H, Zhu Y, Chen M, Lü Y. Abnormal condensation of water vapour at ambient temperature. Phys Chem Chem Phys 2024; 26:8784-8793. [PMID: 38420852 DOI: 10.1039/d3cp05628g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2024]
Abstract
The homogeneous condensation of water vapor at ambient temperature is studied using molecular dynamics simulation. We reveal that there is a droplet size at the nanoscale where water droplets can be stabilized in the condensation process. Our simulations show that the growth of water droplets is dominated by collision and coagulation between small water droplets after nucleation. This process is found to be accompanied by exceptionally fast evaporation such that droplet growth is balanced by evaporation when water droplets grow to a critical size, approximately 12.5 Å in radius, reaching a stable size distribution. The extremely high evaporation rate is attributed to the curvature dependence of surface tension. Surface tension shows a significant decrease with decreasing droplet size below 20 Å, which causes the total free energy of nanoscaled water droplets to rise after collision and coagulation. Consequently, water droplets have to shrink via fast evaporation. The curvature dependence of surface tension is related to the dielectric ordering of water molecules near the surface of water droplets. Owing to fast evaporation, secondary condensation occurs, and many small water clusters form, ultimately exhibiting a bimodal distribution of water-droplet size.
Collapse
Affiliation(s)
- Chenchen Guo
- School of Physics, Beijing Institute of Technology, Beijing 100081, China.
- Changping School Attached to Beijing Normal University, Beijing 102206, China
| | - Kun Yang
- School of Physics, Beijing Institute of Technology, Beijing 100081, China.
| | - Hairong Qin
- School of Physics, Beijing Institute of Technology, Beijing 100081, China.
| | - Yong Zhu
- Science and Technology on Electromagnetic Scattering Laboratory, Beijing 100854, China
| | - Min Chen
- Department of Engineering Mechanics, Center for Nano and Micro Mechanics, Tsinghua University, Beijing 100084, China
| | - Yongjun Lü
- School of Physics, Beijing Institute of Technology, Beijing 100081, China.
| |
Collapse
|
4
|
Fijan D, Wilson M. Thermodynamic anomalies, polyamorphism and all that. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220336. [PMID: 37634531 PMCID: PMC10460645 DOI: 10.1098/rsta.2022.0336] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2023] [Accepted: 05/18/2023] [Indexed: 08/29/2023]
Abstract
The appearance and evolution of thermodynamics anomalies, and related properties, are studied for two classes of system, modelling those dominated by covalent and ionic interactions, respectively. Such anomalies are most familiar in the density but are also present in other thermodynamics variables such as the compressibility and heat capacity. By systematically varying key model parameters the emergence and evolution of these anomalies can be tracked across the phase space. The interaction of the anomalies can often be rationalized by thermodynamics 'rules'. The emergence of these anomalies may also be correlated with the appearance of polyamorphism, the existence of multiple amorphous states which differ in density and entropy. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.
Collapse
Affiliation(s)
- Domagoj Fijan
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford, UK
| | - Mark Wilson
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford, UK
| |
Collapse
|
5
|
Perin L, Gallo P. Phase Diagram of Aqueous Solutions of LiCl: a Study of Concentration Effects on the Anomalies of Water. J Phys Chem B 2023; 127:4613-4622. [PMID: 37167579 DOI: 10.1021/acs.jpcb.3c00703] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
We perform molecular dynamics simulations in order to study thermodynamics and the structure of supercooled aqueous solutions of lithium chloride (LiCl) at concentrations c = 0.678 and 2.034 mol/kg. We model the solvent using the TIP4P/2005 potential and the ions using the Madrid-2019 force field, a force field particularly suited for studying this solution. We find that, for c = 0.678 mol/kg, the behavior of the equation of state, studied in the P-T plane, indicates the presence of a liquid-liquid phase transition, similar to what was previously found for bulk water. We estimate the position of the liquid-liquid critical point to be at Tc ≈ 174 K, Pc ≈ 1775 bar, and ρc ≈ 1.065 g/cm3. When the concentration is tripled to c = 2.034 mol/kg, no critical point is observed, indicating its possible disappearance at this concentration. We also study the water-water and water-ions structure in the two solutions, and we find that at the concentrations examined the effect of ions on the water-water structure is not strong, and all the features found in bulk water are preserved. We also calculate the hydration number of the Li and Cl ions, and in line with experiments, we find the value of 4 for Li+ and between 5.5 and 6 for Cl-, confirming the good performances of the Madrid-2019 force field.
Collapse
Affiliation(s)
- Leonardo Perin
- Dipartimento di Fisica, Università Roma Tre, Via della Vasca Navale 84, I-00146 Roma, Italy
| | - Paola Gallo
- Dipartimento di Fisica, Università Roma Tre, Via della Vasca Navale 84, I-00146 Roma, Italy
| |
Collapse
|
6
|
Škvára J, Nezbeda I. Thermodynamics and structure of supercooled water. II. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120508] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
7
|
Agrawal A, Gravelle S, Kamal C, Botto L. Viscous peeling of a nanosheet. SOFT MATTER 2022; 18:3967-3980. [PMID: 35551304 PMCID: PMC9131316 DOI: 10.1039/d1sm01743h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 05/02/2022] [Indexed: 06/15/2023]
Abstract
Combining molecular dynamics (MD) and continuum simulations, we study the dynamics of propagation of a peeling front in a system composed of multilayered graphene nanosheets completely immersed in water. Peeling is induced by lifting one of the nanosheet edges with an assigned pulling velocity normal to the flat substrate. Using MD, we compute the pulling force as a function of the pulling velocity, and quantify the viscous resistance to the advancement of the peeling front. We compare the MD results to a 1D continuum model of a sheet loaded with modelled hydrodynamic loads. Our results show that the viscous dependence of the force on the velocity is negligible below a threshold velocity. Above this threshold, the hydrodynamics is mainly controlled by the viscous resistance associated to the flow near the crack opening, while lubrication forces are negligible owing to the large hydrodynamic slip at the liquid-solid boundary. Two dissipative mechanisms are identified: a drag resistance to the upward motion of the edge, and a resistance to the gap opening associated to the curvature of the flow streamlines near the entrance. Surprisingly, the shape of the sheet was found to be approximately independent of the pulling velocity even for the largest velocities considered.
Collapse
Affiliation(s)
- Adyant Agrawal
- School of Engineering and Material Science, Queen Mary University of London, London, UK
| | - Simon Gravelle
- School of Engineering and Material Science, Queen Mary University of London, London, UK
| | - Catherine Kamal
- School of Engineering and Material Science, Queen Mary University of London, London, UK
| | - Lorenzo Botto
- Process and Energy Department, 3ME Faculty of Mechanical, Maritime and Materials Engineering, TU Delft, Delft, The Netherlands.
| |
Collapse
|
8
|
Eltareb A, Lopez GE, Giovambattista N. Evidence of a liquid-liquid phase transition in H[Formula: see text]O and D[Formula: see text]O from path-integral molecular dynamics simulations. Sci Rep 2022; 12:6004. [PMID: 35397618 PMCID: PMC8994788 DOI: 10.1038/s41598-022-09525-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2021] [Accepted: 02/23/2022] [Indexed: 01/22/2023] Open
Abstract
We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H[Formula: see text]O and D[Formula: see text]O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density [Formula: see text], isothermal compressibility [Formula: see text], and self-diffusion coefficients D(T) of H[Formula: see text]O and D[Formula: see text]O are in excellent agreement with available experimental data; the isobaric heat capacity [Formula: see text] obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H[Formula: see text]O and D[Formula: see text]O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H[Formula: see text]O and D[Formula: see text]O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H[Formula: see text]O, from PIMD simulations, is located at [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D[Formula: see text]O is estimated to be [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of [Formula: see text] for D[Formula: see text]O and, particularly, H[Formula: see text]O suggest that improved water models are needed for the study of supercooled water.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, New York 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, New York 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
|
9
|
Longo TJ, Anisimov MA. Phase transitions affected by natural and forceful molecular interconversion. J Chem Phys 2022; 156:084502. [DOI: 10.1063/5.0081180] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
If a binary liquid mixture, composed of two alternative species with equal amounts, is quenched from a high temperature to a low temperature, below the critical point of demixing, then the mixture will phase separate through a process known as spinodal decomposition. However, if the two alternative species are allowed to interconvert, either naturally (e.g., the equilibrium interconversion of enantiomers) or forcefully (e.g., via an external source of energy or matter), then the process of phase separation may drastically change. In this case, depending on the nature of interconversion, two phenomena could be observed: either phase amplification, the growth of one phase at the expense of another stable phase, or microphase separation, the formation of nongrowing (steady-state) microphase domains. In this work, we phenomenologically generalize the Cahn–Hilliard theory of spinodal decomposition to include the molecular interconversion of species and describe the physical properties of systems undergoing either phase amplification or microphase separation. We apply the developed phenomenology to accurately describe the simulation results of three atomistic models that demonstrate phase amplification and/or microphase separation. We also discuss the application of our approach to phase transitions in polyamorphic liquids. Finally, we describe the effects of fluctuations of the order parameter in the critical region on phase amplification and microphase separation.
Collapse
Affiliation(s)
- Thomas J. Longo
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - Mikhail A. Anisimov
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
10
|
Li P, Lu H, Fu YQ. Phase transition of supercooled water confined in cooperative two-state domain. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:165403. [PMID: 35114647 DOI: 10.1088/1361-648x/ac519b] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Accepted: 02/03/2022] [Indexed: 06/14/2023]
Abstract
The question of 'what is the structure of water?' has been regarded as one of the major scientific conundrums in condensed-matter physics due to the complex phase behavior and condensed structure of supercooled water. Great effort has been made so far using both theoretical analysis based on various mathematical models and computer simulations such as molecular dynamics and first-principle. However, these theoretical and simulation studies often do not have strong evidences of condensed-matter physics to support. In this study, a cooperative domain model is formulated to describe the dynamic phase transition of supercooled water between supercooled water and amorphous ice, both of which are composed of low- and high-density liquid water. Free volume theory is initially employed to identify the working principle of dynamic phase transition and its connection to glass transition in the supercooled water. Then a cooperative two-state model is developed to characterize the dynamic anomalies of supercooled water, including density, viscosity and self-diffusion coefficient. Finally, the proposed model is verified using the experimental results reported in literature.
Collapse
Affiliation(s)
- Peizhao Li
- Science and Technology on Advanced Composites in Special Environments Laboratory, Harbin Institute of Technology, Harbin 150080, People's Republic of China
| | - Haibao Lu
- Science and Technology on Advanced Composites in Special Environments Laboratory, Harbin Institute of Technology, Harbin 150080, People's Republic of China
| | - Yong-Qing Fu
- Faculty of Engineering and Environment, Northumbria University, Newcastle upon Tyne, NE1 8ST, United Kingdom
| |
Collapse
|
11
|
de Almeida Ribeiro I, de Koning M, Molinero V. Is It Possible to Follow the Structural Evolution of Water in "No-Man's Land" Using a Pulsed-Heating Procedure? J Phys Chem Lett 2022; 13:1085-1089. [PMID: 35080178 DOI: 10.1021/acs.jpclett.1c04106] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The anomalous increase in compressibility and heat capacity of supercooled water has been attributed to its structural transformation of into a four-coordinated liquid. Experiments revealed that κT and Cp peak at TWthermo ≈ 229 K [Kim et al. Science 2017, 358, 1589; Pathak et al. Proc. Natl. Acad. Sci. 2021, 118, e2018379118]. Recently, a pulsed heating procedure (PHP) was employed to interrogate the structure of water, reporting a steep increase in tetrahedrality around TWPHP = 210 ± 3 K [Kringle et al. Science 2020, 369, 1490]. This discrepancy questions whether water structure and thermodynamics are decoupled, or if the shift in TW is an artifact of PHP. Here we implement PHP in molecular simulations. We find that the stationary states captured at the bottom of the pulse are not representative of the thermalized liquid or its inherent structure. Our analysis reveals a temperature-dependent distortion that shifts TWPHP to ∼20 K below TWthermo. We conclude that 2 orders of magnitude faster rates are required to sample water's inherent structure with PHP.
Collapse
Affiliation(s)
- Ingrid de Almeida Ribeiro
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Maurice de Koning
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861, Campinas, São Paulo, Brazil
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
12
|
Russo J, Leoni F, Martelli F, Sciortino F. The physics of empty liquids: from patchy particles to water. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2022; 85:016601. [PMID: 34905739 DOI: 10.1088/1361-6633/ac42d9] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/16/2021] [Accepted: 12/14/2021] [Indexed: 06/14/2023]
Abstract
Empty liquids represent a wide class of materials whose constituents arrange in a random network through reversible bonds. Many key insights on the physical properties of empty liquids have originated almost independently from the study of colloidal patchy particles on one side, and a large body of theoretical and experimental research on water on the other side. Patchy particles represent a family of coarse-grained potentials that allows for a precise control of both the geometric and the energetic aspects of bonding, while water has arguably the most complex phase diagram of any pure substance, and a puzzling amorphous phase behavior. It was only recently that the exchange of ideas from both fields has made it possible to solve long-standing problems and shed new light on the behavior of empty liquids. Here we highlight the connections between patchy particles and water, focusing on the modelling principles that make an empty liquid behave like water, including the factors that control the appearance of thermodynamic and dynamic anomalies, the possibility of liquid-liquid phase transitions, and the crystallization of open crystalline structures.
Collapse
Affiliation(s)
- John Russo
- Department of Physics, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Fabio Leoni
- Department of Physics, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Fausto Martelli
- IBM Research Europe, Hartree Centre, Daresbury WA4 4AD, United Kingdom
| | - Francesco Sciortino
- Department of Physics, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185 Rome, Italy
| |
Collapse
|
13
|
Malek SMA, Kwan V, Saika-Voivod I, Consta S. Low Density Interior in Supercooled Aqueous Nanodroplets Expels Ions to the Subsurface. J Am Chem Soc 2021; 143:13113-13123. [PMID: 34375522 DOI: 10.1021/jacs.1c04142] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The interaction between water and ions within droplets plays a key role in the chemical reactivity of atmospheric and man-made aerosols. Here we report direct computational evidence that in supercooled aqueous nanodroplets a lower density core of tetrahedrally coordinated water expels the cosmotropic ions to the denser and more disordered subsurface. In contrast, at room temperature, depending on the nature of the ion, the radial distribution in the droplet core is nearly uniform or elevated toward the center. We analyze the spatial distribution of a single ion in terms of a reference electrostatic model. The energy of the system in the analytical model is expressed as the sum of the electrostatic and surface energy of a deformable droplet. The model predicts that the ion is subject to a harmonic potential centered at the droplet's center of mass. We name this effect "electrostatic confinement". The model's predictions are consistent with the simulation findings for a single ion at room temperature but not at supercooling. We anticipate this study to be the starting point for investigating the structure of supercooled (electro)sprayed droplets that are used to preserve the conformations of macromolecules originating from the bulk solution.
Collapse
Affiliation(s)
- Shahrazad M A Malek
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's A1B 3X7, Canada
| | - Victor Kwan
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Ivan Saika-Voivod
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's A1B 3X7, Canada.,Department of Applied Mathematics, Western University, London, Ontario N6A 3K7, Canada
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
14
|
Eltareb A, Lopez GE, Giovambattista N. Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water. Phys Chem Chem Phys 2021; 23:6914-6928. [PMID: 33729222 PMCID: PMC9311234 DOI: 10.1039/d0cp04325g] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We perform path-integral molecular dynamics (PIMD) simulations of H2O and D2O using the q-TIP4P/F model. Simulations are performed at P = 1 bar and over a wide range of temperatures that include the equilibrium (T≥ 273 K) and supercooled (210 ≤T < 273 K) liquid states of water. The densities of both H2O and D2O calculated from PIMD simulations are in excellent agreement with experiments in the equilibrium and supercooled regimes. We also evaluate important thermodynamic response functions, specifically, the thermal expansion coefficient αP(T), isothermal compressibility κT(T), isobaric heat capacity CP(T), and static dielectric constant ε(T). While these properties are in excellent [αP(T) and κT(T)] or semi-quantitative agreement [CP(T) and ε(T)] with experiments in the equilibrium regime, they are increasingly underestimated upon further cooling. It follows that the inclusion of nuclear quantum effects in PIMD simulations of (q-TIP4P/F) water is not sufficient to reproduce the anomalous large fluctuations in density, entropy, and electric dipole moment characteristic of supercooled water. It has been hypothesized that water may exhibit a liquid-liquid critical point (LLCP) in the supercooled regime at P > 1 bar and that such a LLCP generates a maximum in CP(T) and κT(T) at 1 bar. Consistent with this hypothesis and in particular, with experiments, we find a maximum in the κT(T) of q-TIP4P/F light and heavy water at T≈ 230-235 K. No maximum in CP(T) could be detected down to T≥ 210 K. We also calculate the diffusion coefficient D(T) of H2O and D2O using the ring-polymer molecular dynamics (RPMD) technique and find that computer simulations are in remarkable good agreement with experiments at all temperatures studied. The results from RPMD/PIMD simulations are also compared with the corresponding results obtained from classical MD simulations of q-TIP4P/F water where atoms are represented by single interacting sites. Surprisingly, we find minor differences in most of the properties studied, with CP(T), D(T), and structural properties being the only (expected) exceptions.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY 11210, USA.
| | | | | |
Collapse
|
15
|
Abstract
Many important processes affecting the earth's climate are determined by the physical properties of seawater. In addition, desalination of seawater is a significant source of drinking water for the human population living in coastal areas. Since the physical properties of seawater governing these processes depend on the molecular interactions among its components, a deeper knowledge of seawater at the molecular level would contribute to a better understanding of these phenomena. However, in strong contrast with the situation in other areas such as biomolecules or materials science, molecular simulation studies reporting the physical properties of seawater are currently lacking. This is probably due to the usual perception of the seawater composition being too complex to approach. This point of view ignores the fact that physical properties of seawater are dependent on a single parameter representing the composition, namely the salinity. This is because the relative proportions of any two major constituents of seasalt are always the same. Another obstacle to performing molecular simulations of seawater could have been the unavailability of a satisfactory force field representing the interactions between water molecules and dissolved substances. However, this drawback has recently been overcome with the proposal of the Madrid-2019 force field. In this work we show for the first time that molecular simulation of seawater is feasible. We have performed molecular dynamics simulations of a system, the composition of which is close to the average composition of standard seawater and with the molecular interactions given by the Madrid-2019 force field. In this way we are able to provide quantitative or semiquantitative predictions for a number of relevant physical properties of seawater for temperatures and salinities from the oceanographic range to those relevant to desalination processes. The computed magnitudes include static (density), dynamical (viscosity and diffusion coefficients), structural (ionic hydration, ion-ion distribution functions), and interfacial (surface tension) properties.
Collapse
Affiliation(s)
- Ivan M Zeron
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Miguel A Gonzalez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Edoardo Errani
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L F Abascal
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
16
|
Puibasset J, Judeinstein P, Zanotti JM. Bulk supercooled water versus adsorbed films on silica surfaces: specific heat by Monte Carlo simulation. Phys Chem Chem Phys 2021; 23:2275-2285. [PMID: 33443254 DOI: 10.1039/d0cp05387b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Between 150 and 230.6 K, bulk supercooled water freezes upon cooling, and amorphous ice crystallizes upon heating: bulk water thus exists only in its stable ice form. To circumvent this problem, experiments are generally performed on water adsorbed in SiO2 based porous systems. In this work, we take advantage of Monte Carlo simulations to explore this metastable supercooled region inaccessible to experiments. Using three rigid, non-polarizable water models, namely SPC, TIP4P and TIP4P/2005, we investigate the isobaric specific heat capacity (Cp), between 100 and 300 K, of bulk water and water films of few monolayers adsorbed on different SiO2 surfaces: a smooth surface, a non-hydroxylated (0001) surface of quartz, and a fully hydroxylated (001) surface of cristobalite. As Cp is directly related to the entropy fluctuations and we focus on low temperatures, the convergence of the Monte Carlo simulations is a critical point of this work. Also, due to the small mass of the hydrogen atoms, quantum corrections are taken into account, and lead to an excellent agreement of the simulated and experimental Cp values at low temperature (100 K region). Altogether, we conclude that, in bulk, Cp is shown to exhibit a broad peak around 225 K for the SPC and TIP4P models, and around 250 K for the TIP4P/2005 model, in qualitative agreement with the experimentally observed features in Cp measurements. For interfacial water, in all cases, the broad Cp peak disappears. This result, at odds with experimental observations, suggests that disorder and hydrogen bonding at the interface (not yet taken into account) have a fundamental role in confined water transitions.
Collapse
Affiliation(s)
- J Puibasset
- ICMN, UMR7374, CNRS, Université d'Orléans, 1b, Rue de la Férollerie, 45071, Orléans Cedex 2, France.
| | | | | |
Collapse
|
17
|
Cerdeiriña CA, González-Salgado D. Temperature, Pressure, and Length-Scale Dependence of Solvation in Water-like Solvents. I. Small Solvophobic Solutes. J Phys Chem B 2020; 125:297-306. [DOI: 10.1021/acs.jpcb.0c09569] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Claudio A. Cerdeiriña
- Departamento de Física Aplicada, Universidad de Vigo, Campus del Agua, Ourense 32004, Spain
| | - Diego González-Salgado
- Departamento de Física Aplicada, Universidad de Vigo, Campus del Agua, Ourense 32004, Spain
| |
Collapse
|
18
|
Dubey V, Daschakraborty S. Breakdown of the Stokes-Einstein Relation in Supercooled Water/Methanol Binary Mixtures: Explanation Using the Translational Jump-Diffusion Approach. J Phys Chem B 2020; 124:10398-10408. [PMID: 33153260 DOI: 10.1021/acs.jpcb.0c07318] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A recent experiment has directly checked the validity of the Stokes-Einstein (SE) relation for pure water, pure methanol, and their binary mixtures of three different compositions at different temperatures. The effect of composition on the nature of breakdown of the SE relation is interesting. While in the majority of the systems, an increasing SE breakdown is observed with decreasing temperature, the breakdown is already significant at higher temperatures for the equimolar mixture. Violations of the SE relation in pure supercooled water at different temperatures and pressures have been previously explained using the translational jump-diffusion (TJD) approach, which provides a fundamental molecular basis, by directly connecting the SE breakdown with jump-diffusion of the molecules. We have used the same TJD approach for explaining the SE breakdown for the methanol/water binary mixtures of compositions studied in the experiment over a wide range of temperatures between 220 K and 300 K. We have understood that the jump-diffusion is the key responsible factor for the SE breakdown. The maximum jump-diffusion contribution gives rise to the early SE breakdown observed for the equimolar mixture observed in the experiment. This study, therefore, provides molecular insight into the SE breakdown for the supercooled water/methanol binary mixture, as found in the experiment.
Collapse
Affiliation(s)
- Vikas Dubey
- Department of Chemistry, Indian Institute of Technology Patna, Patna, Bihar 801106, India
| | | |
Collapse
|
19
|
Gartner TE, Zhang L, Piaggi PM, Car R, Panagiotopoulos AZ, Debenedetti PG. 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:26040-26046. [PMID: 33008883 PMCID: PMC7584908 DOI: 10.1073/pnas.2015440117] [Citation(s) in RCA: 55] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The possible existence of a metastable liquid-liquid transition (LLT) and a corresponding liquid-liquid critical point (LLCP) in supercooled liquid water remains a topic of much debate. An LLT has been rigorously proved in three empirically parametrized molecular models of water, and evidence consistent with an LLT has been reported for several other such models. In contrast, experimental proof of this phenomenon has been elusive due to rapid ice nucleation under deeply supercooled conditions. In this work, we combined density functional theory (DFT), machine learning, and molecular simulations to shed additional light on the possible existence of an LLT in water. We trained a deep neural network (DNN) model to represent the ab initio potential energy surface of water from DFT calculations using the Strongly Constrained and Appropriately Normed (SCAN) functional. We then used advanced sampling simulations in the multithermal-multibaric ensemble to efficiently explore the thermophysical properties of the DNN model. The simulation results are consistent with the existence of an LLCP, although they do not constitute a rigorous proof thereof. We fit the simulation data to a two-state equation of state to provide an estimate of the LLCP's location. These combined results-obtained from a purely first-principles approach with no empirical parameters-are strongly suggestive of the existence of an LLT, bolstering the hypothesis that water can separate into two distinct liquid forms.
Collapse
Affiliation(s)
- Thomas E Gartner
- Department of Chemistry, Princeton University, Princeton, NJ 08544
| | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ 08544
| | - Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, NJ 08544
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ 08544
- Department of Physics, Princeton University, Princeton, NJ 08544
- Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, NJ 08544
| | - Athanassios Z Panagiotopoulos
- Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, NJ 08544;
- Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ 08544
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ 08544
| |
Collapse
|
20
|
Zaragoza A, Tripathi CSP, Gonzalez MA, Abascal JLF, Caupin F, Valeriani C. Effect of dissolved salt on the anomalies of water at negative pressure. J Chem Phys 2020; 152:194501. [PMID: 33687252 DOI: 10.1063/5.0002745] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Adding salt to water at ambient pressure affects its thermodynamic properties. At low salt concentration, anomalies such as the density maximum are shifted to lower temperature, while at large enough salt concentration, they cannot be observed any more. Here, we investigate the effect of salt on an anomaly recently observed in pure water at negative pressure: the existence of a sound velocity minimum along isochores. We compare experiments and simulations for an aqueous solution of sodium chloride with molality around 1.2 mol kg-1, reaching pressures beyond -100 MPa. We also discuss the origin of the minima in the sound velocity and emphasize the importance of the relative position of the temperatures of sound velocity and density anomalies.
Collapse
Affiliation(s)
- Alberto Zaragoza
- Departamento Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | | | - Miguel A Gonzalez
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - José Luis F Abascal
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Frédéric Caupin
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Chantal Valeriani
- Departamento Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
21
|
Montazeri K, Hao S, Abdolhosseini Qomi MJ, Won Y. Molecular Dynamics Investigation of Liquid and Vapor Interactions Near an Evaporating Interface: A Theoretical Genetics Perspective. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000017] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Kimia Montazeri
- Department of Mechanical and Aerospace EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| | - Shuai Hao
- Department of Mechanical and Aerospace EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| | - Mohammad Javad Abdolhosseini Qomi
- Department of Civil and Environmental EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| | - Yoonjin Won
- Department of Mechanical and Aerospace EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| |
Collapse
|
22
|
Taschin A, Bartolini P, Fanetti S, Lapini A, Citroni M, Righini R, Bini R, Torre R. Pressure Effects on Water Dynamics by Time-Resolved Optical Kerr Effect. J Phys Chem Lett 2020; 11:3063-3068. [PMID: 32212735 DOI: 10.1021/acs.jpclett.0c00363] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Despite water being the most common and most widely studied substance in the world, it still presents unknown aspects. In particular, water shows several thermodynamic and dynamical anomalies in the liquid and supercooled metastable phases, and the natures of these phases are still hotly debated. Here, we report measurements of water using the optical Kerr effect as a function of pressure along two isotherms, at 273 K from 0.1 to 750 MPa and at 297 K from 0.1 to 1350 MPa, reaching the supercooled metastable phase. The structural relaxation and the low frequency vibrational dynamics of water show a peculiar pressure dependence similar to that of other dynamical properties. The data analysis suggests the presence in the water phase diagram of a crossover area that divides two regions characterized by different dynamic regimes, which appear to be related to two liquid forms, one dominated by the high density water and the other by the low density water.
Collapse
Affiliation(s)
- Andrea Taschin
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
| | - Paolo Bartolini
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
| | - Samuele Fanetti
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
- ICCOM-CNR, Istituto di Chimica dei Composti OrganoMetallici, Consiglio Nazionale delle Ricerche, Via Madonna del Piano 10, I-50019 Sesto Fiorentino, Firenze, Italy
| | - Andrea Lapini
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
- INRIM, Istituto Nazionale di Ricerca Metrologica, Strada delle Cacce 91, I-10135 Torino, Italy
| | - Margherita Citroni
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
| | - Roberto Righini
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
- Dipartimento di Chimica "Ugo Schiff", Università di Firenze, via Della Lastruccia 13, I-50019 Sesto Fiorentino, Firenze, Italy
| | - Roberto Bini
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
- ICCOM-CNR, Istituto di Chimica dei Composti OrganoMetallici, Consiglio Nazionale delle Ricerche, Via Madonna del Piano 10, I-50019 Sesto Fiorentino, Firenze, Italy
- Dipartimento di Chimica "Ugo Schiff", Università di Firenze, via Della Lastruccia 13, I-50019 Sesto Fiorentino, Firenze, Italy
| | - Renato Torre
- European Laboratory for Non-Linear Spectroscopy (LENS), Università di Firenze, via N. Carrara 1, I-50019 Sesto Fiorentino, Firenze, Italy
- Dipartimento di Fisica e Astronomia, Università di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Firenze, Italy
| |
Collapse
|
23
|
Wang Q, Tang Q, Tian S. Molecular dynamics simulation of sI methane hydrate under compression and tension. OPEN CHEM 2020. [DOI: 10.1515/chem-2020-0008] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
AbstractMolecular dynamics (MD) analysis of methane hydrate is important for the application of methane hydrate technology. This study investigated the microstructure changes of sI methane hydrate and the laws of stress–strain evolution under the condition of compression and tension by using MD simulation. This study further explored the mechanical property and stability of sI methane hydrate under different stress states. Results showed that tensile and compressive failures produced an obvious size effect under a certain condition. At low temperature and high pressure, most of the clathrate hydrate maintained a stable structure in the tensile fracture process, during which only a small amount of unstable methane broke the structure, thereby, presenting a free-motion state. The methane hydrate cracked when the system reached the maximum stress in the loading process, in which the maximum compressive stress is larger than the tensile stress under the same experimental condition. This study provides a basis for understanding the microscopic stress characteristics of methane hydrate.
Collapse
Affiliation(s)
- Qiang Wang
- State Key Laboratory of Coal Mine Disaster Dynamics and Control, School of Resources and Safety Engineering, Chongqing University, Chongqing400044, China
- Key Laboratory of Low-grade Energy Utilization Technology and System, Ministry of Education, College of Energy and Power Engineering, Chongqing University, Chongqing400044, China
| | - Qizhong Tang
- Chongqing Key Laboratory of Heterogeneous Material Mechanics, College of Aerospace Engineering, Chongqing University, Chongqing400044, China
| | - Sen Tian
- State Key Laboratory of Coal Mine Disaster Dynamics and Control, School of Resources and Safety Engineering, Chongqing University, Chongqing400044, China
- School of Civil Engineering, The University of Queensland, Brisbane, QLD 4072, Australia
| |
Collapse
|
24
|
Pathak H, Späh A, Amann-Winkel K, Perakis F, Kim KKH, Nilsson A. Temperature dependent anomalous fluctuations in water: shift of ≈1 kbar between experiment and classical force field simulations. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1649486] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Harshad Pathak
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Alexander Späh
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Fivos Perakis
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Kyung Kyung Hwan Kim
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Anders Nilsson
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| |
Collapse
|
25
|
Dueby S, Dubey V, Daschakraborty S. Decoupling of Translational Diffusion from the Viscosity of Supercooled Water: Role of Translational Jump Diffusion. J Phys Chem B 2019; 123:7178-7189. [DOI: 10.1021/acs.jpcb.9b01719] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shivam Dueby
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India
| | - Vikas Dubey
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India
| | | |
Collapse
|
26
|
Caupin F, Anisimov MA. Thermodynamics of supercooled and stretched water: Unifying two-structure description and liquid-vapor spinodal. J Chem Phys 2019; 151:034503. [PMID: 31325919 DOI: 10.1063/1.5100228] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
We have applied a two-structure approach to the description of the thermodynamic properties of supercooled and stretched water, metastable toward vapor, ice, or both, by incorporating the stability limit of liquid with respect to vapor at negative pressures. In addition to the properties of water considered in previous studies, we include new data recently obtained in deeply supercooled and stretched regions. Our model reproduces the experimentally observed anomalies in metastable water up to 400 MPa and down to -140 MPa, and can provide a physically based extrapolation in regions where no measurements are available yet. Moreover, we are able to elucidate the thermodynamic nature of the alternative "states" of liquid water, namely, high-temperature denser water (state A) and "mother-of-ice" lighter water (state B). Based on the internal consistency of the described anomalies and new data on the isothermal compressibility, we exclude the critical-point-free scenario in which the first-order liquid-liquid transition line would continue into the stretched liquid state (doubly metastable) crossing the vapor-liquid spinodal. A "singularity-free" scenario remains an option for explaining supercooled water's anomalies within the framework of two-state thermodynamics; however, the extreme case of the singularity-free scenario, ideal mixing of A and B, seems improbable. We have also clarified the concept of fast interconversion of alternative states in supercooled water as a phenomenological representation of distribution of short-ranged local structures.
Collapse
Affiliation(s)
- Frédéric Caupin
- Institut Lumière Matière, Université Claude Bernard Lyon 1, CNRS, Université de Lyon, F-69622, Villeurbanne, France
| | - Mikhail A Anisimov
- Department of Chemical and Biomolecular Engineering and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
27
|
Spurious violation of the Stokes-Einstein-Debye relation in supercooled water. Sci Rep 2019; 9:8118. [PMID: 31148561 PMCID: PMC6544661 DOI: 10.1038/s41598-019-44517-4] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2019] [Accepted: 05/17/2019] [Indexed: 11/12/2022] Open
Abstract
The theories of Brownian motion, the Debye rotational diffusion model, and hydrodynamics together provide us with the Stokes–Einstein–Debye (SED) relation between the rotational relaxation time of the \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\boldsymbol{\ell }}$$\end{document}ℓ-th degree Legendre polynomials \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ, and viscosity divided by temperature, η/T. Experiments on supercooled liquids are frequently performed to measure the SED relations, \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓkBT/η and Dt\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ, where Dt is the translational diffusion constant. However, the SED relations break down, and its molecular origin remains elusive. Here, we assess the validity of the SED relations in TIP4P/2005 supercooled water using molecular dynamics simulations. Specifically, we demonstrate that the higher-order \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ values exhibit a temperature dependence similar to that of η/T, whereas the lowest-order \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ values are decoupled with η/T, but are coupled with the translational diffusion constant Dt. We reveal that the SED relations are so spurious that they significantly depend on the degree of Legendre polynomials.
Collapse
|
28
|
Uralcan B, Latinwo F, Debenedetti PG, Anisimov MA. Pattern of property extrema in supercooled and stretched water models and a new correlation for predicting the stability limit of the liquid state. J Chem Phys 2019; 150:064503. [DOI: 10.1063/1.5078446] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Betul Uralcan
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Folarin Latinwo
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Mikhail A. Anisimov
- Department of Chemical and Biomolecular Engineering and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
29
|
Abstract
Liquid water is considered poorly understood. How are water's physical properties encoded in its molecular structure? We introduce a statistical mechanical model (CageWater) of water's hydrogen-bonding (HB) and Lennard-Jones (LJ) interactions. It predicts the energetic and volumetric and anomalous properties accurately. Yet, because the model is analytical, it is essentially instantaneous to compute. This model advances our understanding beyond current molecular simulations and experiments. Water has long been regarded as a "2-density liquid": a dense LJ liquid and a looser HB one. Instead, we find here a different antagonism underlying water structure-property relations: HBs in water-water pairs drive density, while HBs in cooperative cages drive openness. The balance shifts strongly with temperature and pressure. This model interprets the molecular structures underlying the liquid-liquid phase transition in supercooled water. It may have value in geophysics, biomolecular modeling, and engineering of materials for water purification and green chemistry.
Collapse
Affiliation(s)
- Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, Ljubljana SI-1000, Slovenia
| | - Ken A. Dill
- Laufer Center for Physical and Quantitative Biology and Departments of Chemistry and of Physics & Astronomy, Stony Brook University, Stony Brook, New York 11794-5252, United States
| |
Collapse
|
30
|
Montero de Hijes P, Sanz E, Joly L, Valeriani C, Caupin F. Viscosity and self-diffusion of supercooled and stretched water from molecular dynamics simulations. J Chem Phys 2018; 149:094503. [DOI: 10.1063/1.5042209] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Pablo Montero de Hijes
- Departamento de Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, Madrid 28040, Spain
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Laurent Joly
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Lyon, France
| | - Chantal Valeriani
- Departamento de Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, Madrid 28040, Spain
| | - Frédéric Caupin
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Lyon, France
| |
Collapse
|
31
|
Silletta EV, Tuckerman ME, Jerschow A. Unusual Proton Transfer Kinetics in Water at the Temperature of Maximum Density. PHYSICAL REVIEW LETTERS 2018; 121:076001. [PMID: 30169046 DOI: 10.1103/physrevlett.121.076001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2018] [Revised: 06/26/2018] [Indexed: 06/08/2023]
Abstract
Water exhibits numerous anomalous properties, many of which remain poorly understood. One of its intriguing behaviors is that it exhibits a temperature of maximum density (TMD) at 4 °C. We provide here new experimental evidence for hitherto unknown abrupt changes in proton transfer kinetics at the TMD. In particular, we show that the lifetime of OH^{-} ions has a maximum at this temperature, in contrast to hydronium ions. Furthermore, base-catalyzed proton transfer shows a sharp local minimum at this temperature, and activation energies change abruptly as well. The measured lifetimes agree with earlier theoretical predictions as the temperature approaches the TMD. Similar results are also found for heavy water at its own TMD. These findings point to a high propensity of forming fourfold coordinated OH^{-} solvation complexes at the TMD, underlining the asymmetry between hydroxide and hydronium transport. These results could help to further elucidate the unusual properties of water and related liquids.
Collapse
Affiliation(s)
- Emilia V Silletta
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Mark E Tuckerman
- Department of Chemistry, New York University, New York, New York 10003, USA
- Courant Institute of Mathematical Science, New York University, New York, New York 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| | - Alexej Jerschow
- Department of Chemistry, New York University, New York, New York 10003, USA
| |
Collapse
|
32
|
|
33
|
Singraber A, Morawietz T, Behler J, Dellago C. Density anomaly of water at negative pressures from first principles. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:254005. [PMID: 29762140 DOI: 10.1088/1361-648x/aac4f4] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Using molecular dynamics simulations based on ab initio trained high-dimensional neural network potentials, we study the equation of state of liquid water at negative pressures. From density isobars computed for various pressures down to p = -230 MPa we determine the line of density maxima for two potentials based on the BLYP and the RPBE functionals, respectively. In both cases, dispersion corrections are included to account for non-local long-range correlations that give rise to van der Waals forces. We have followed the density maximum down to negative pressures close to the spinodal instability. For both functionals, the temperature of maximum density increases with decreasing pressure under moderate stretching, but changes slope at [Formula: see text] MPa and [Formula: see text] MPa for BLYP and RPBE, respectively. Our calculations confirm qualitatively the retracing shape of the line of density maxima found for empirical water models, indicating that the spinodal line maintains a positive slope even at strongly negative pressures.
Collapse
Affiliation(s)
- Andreas Singraber
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | | | | | | |
Collapse
|
34
|
Abstract
We report a numerical study of the statistical properties of the potential energy landscape of TIP4P/2005, one of the most accurate rigid water models. We show that, in the region where equilibrated configurations can be generated, a Gaussian landscape description is able to properly describe the model properties. We also find that the volume dependence of the landscape properties is consistent with the existence of a locus of density maxima in the phase diagram. The landscape-based equation of state accurately reproduces the TIP4P/2005 pressure-vs-volume curves, providing a sound extrapolation of the free-energy at low T. A positive-pressure liquid-liquid critical point is predicted by the resulting free-energy.
Collapse
Affiliation(s)
- Philip H Handle
- Department of Physics, Sapienza-University of Rome, Piazzale Aldo Moro 5, I-00185 Roma, Italy
| | - Francesco Sciortino
- Department of Physics, Sapienza-University of Rome, Piazzale Aldo Moro 5, I-00185 Roma, Italy
| |
Collapse
|
35
|
Abstract
Water is the most common and yet least understood material on Earth. Despite its simplicity, water tends to form tetrahedral order locally by directional hydrogen bonding. This structuring is known to be responsible for a vast array of unusual properties, e.g., the density maximum at 4 °C, which play a fundamental role in countless natural and technological processes, with the Earth’s climate being one of the most important examples. By systematically tuning the degree of tetrahedrality, we succeed in continuously interpolating between water-like behavior and simple liquid-like behavior. Our approach reveals what physical factors make water so anomalous and special even compared with other tetrahedral liquids. Tetrahedral interactions describe the behavior of the most abundant and technologically important materials on Earth, such as water, silicon, carbon, germanium, and countless others. Despite their differences, these materials share unique common physical behaviors, such as liquid anomalies, open crystalline structures, and extremely poor glass-forming ability at ambient pressure. To reveal the physical origin of these anomalies and their link to the shape of the phase diagram, we systematically study the properties of the Stillinger–Weber potential as a function of the strength of the tetrahedral interaction λ. We uncover a unique transition to a reentrant spinodal line at low values of λ, accompanied with a change in the dynamical behavior, from non-Arrhenius to Arrhenius. We then show that a two-state model can provide a comprehensive understanding on how the thermodynamic and dynamic anomalies of this important class of materials depend on the strength of the tetrahedral interaction. Our work establishes a deep link between the shape of the phase diagram and the thermodynamic and dynamic properties through local structural ordering in liquids and hints at why water is so special among all substances.
Collapse
|
36
|
Tonauer CM, Seidl-Nigsch M, Loerting T. High-density amorphous ice: nucleation of nanosized low-density amorphous ice. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:034002. [PMID: 29189205 DOI: 10.1088/1361-648x/aa9e76] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The pressure dependence of the crystallization temperature of different forms of expanded high-density amorphous ice (eHDA) was scrutinized. Crystallization at pressures 0.05-0.30 GPa was followed using volumetry and powder x-ray diffraction. eHDA samples were prepared via isothermal decompression of very high-density amorphous ice at 140 K to different end pressures between 0.07-0.30 GPa (eHDA0.07-0.3). At 0.05-0.17 GPa the crystallization line T x (p) of all eHDA variants is the same. At pressures >0.17 GPa, all eHDA samples decompressed to pressures <0.20 GPa exhibit significantly lower T x values than eHDA0.2 and eHDA0.3. We rationalize our findings with the presence of nanoscaled low-density amorphous ice (LDA) seeds that nucleate in eHDA when it is decompressed to pressures <0.20 GPa at 140 K. Below ~0.17 GPa, these nanosized LDA domains are latent within the HDA matrix, exhibiting no effect on T x of eHDA<0.2. Upon heating at pressures ⩾0.17 GPa, these nanosized LDA nuclei transform to ice IX nuclei. They are favored sites for crystallization and, hence, lower T x . By comparing crystallization experiments of bulk LDA with the ones involving nanosized LDA we are able to estimate the Laplace pressure and radius of ~0.3-0.8 nm for the nanodomains of LDA. The nucleation of LDA in eHDA revealed here is evidence for the first-order-like nature of the HDA → LDA transition, supporting water's liquid-liquid transition scenarios.
Collapse
Affiliation(s)
- Christina M Tonauer
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria
| | | | | |
Collapse
|
37
|
Holten V, Qiu C, Guillerm E, Wilke M, Rička J, Frenz M, Caupin F. Compressibility Anomalies in Stretched Water and Their Interplay with Density Anomalies. J Phys Chem Lett 2017; 8:5519-5522. [PMID: 29043801 DOI: 10.1021/acs.jpclett.7b02563] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Water keeps puzzling scientists because of its numerous properties which behave oppositely to those of usual liquids: for instance, water expands upon cooling, and liquid water is denser than ice. To explain this anomalous behavior, several theories have been proposed, with different predictions for the properties of supercooled water (liquid at conditions where ice is stable). However, discriminating between those theories with experiments has remained elusive because of spontaneous ice nucleation. Here we measure the sound velocity in liquid water stretched to negative pressure and derive an experimental equation of state, which reveals compressibility anomalies. We show by rigorous thermodynamic relations how these anomalies are intricately linked with the density anomaly. Some features we observe are necessary conditions for the validity of two theories of water.
Collapse
Affiliation(s)
- Vincent Holten
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Chen Qiu
- Institute of Applied Physics, University of Bern , Sidlerstrasse 5, 3012 Bern, Switzerland
| | - Emmanuel Guillerm
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Max Wilke
- Deutsches GeoForschungsZentrum and Universität Potsdam, Erd- und Umweltwissenschaften, Karl-Liebknecht-Str. 24-25, 14473 Potsdam, Germany
| | - Jaroslav Rička
- Institute of Applied Physics, University of Bern , Sidlerstrasse 5, 3012 Bern, Switzerland
| | - Martin Frenz
- Institute of Applied Physics, University of Bern , Sidlerstrasse 5, 3012 Bern, Switzerland
| | - Frédéric Caupin
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| |
Collapse
|
38
|
Altabet YE, Singh RS, Stillinger FH, Debenedetti PG. Thermodynamic Anomalies in Stretched Water. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11771-11778. [PMID: 28885852 DOI: 10.1021/acs.langmuir.7b02339] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Via molecular dynamics simulations of the TIP4P/2005 water model, we study liquid water's anomalous behavior at large negative pressure produced through isochoric cooling. We find that isochores without a pressure minimum can display "reentrant" behavior whereby a system that cavitates upon cooling can then rehomogenize upon further cooling. This behavior is a consequence of the underlying density maximum along the spinodal, but its actual manifestation in simulations is strongly influenced by finite size effects. These observations suggest that water under strong hydrophilic confinement may display richer phase behavior than hitherto assumed. This also suggests that propensity toward cavitation does not always correlate with greater tension, contrary to the prevailing assumption for interpreting water stretching experiments. We also show that a maximum spinodal density in water results in a locus of maximum compressibility and a minimum speed of sound that are independent from any influence of a liquid-liquid critical point (LLCP). However, we demonstrate that structural signatures of a Widom line, which likely emanates from an LLCP at elevated pressure, extend to large negative pressure, but such signatures are only observed upon sampling water's underlying potential energy landscape, rather than the thermalized metastable liquid.
Collapse
Affiliation(s)
- Y Elia Altabet
- Department of Chemical and Biological Engineering, Princeton University , Princeton, New Jersey 08544, United States
| | - Rakesh S Singh
- Department of Chemical and Biological Engineering, Princeton University , Princeton, New Jersey 08544, United States
| | - Frank H Stillinger
- Department of Chemistry, Princeton University , Princeton, New Jersey 08544, United States
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University , Princeton, New Jersey 08544, United States
| |
Collapse
|
39
|
Horstmann R, Vogel M. Common behaviors associated with the glass transitions of water-like models. J Chem Phys 2017; 147:034505. [DOI: 10.1063/1.4993445] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- R. Horstmann
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - M. Vogel
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| |
Collapse
|
40
|
Singh LP, Issenmann B, Caupin F. Pressure dependence of viscosity in supercooled water and a unified approach for thermodynamic and dynamic anomalies of water. Proc Natl Acad Sci U S A 2017; 114:4312-4317. [PMID: 28404733 PMCID: PMC5410815 DOI: 10.1073/pnas.1619501114] [Citation(s) in RCA: 53] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The anomalous decrease of the viscosity of water with applied pressure has been known for over a century. It occurs concurrently with major structural changes: The second coordination shell around a molecule collapses onto the first shell. Viscosity is thus a macroscopic witness of the progressive breaking of the tetrahedral hydrogen bond network that makes water so peculiar. At low temperature, water at ambient pressure becomes more tetrahedral and the effect of pressure becomes stronger. However, surprisingly, no data are available for the viscosity of supercooled water under pressure, in which dramatic anomalies are expected based on interpolation between ambient pressure data for supercooled water and high pressure data for stable water. Here we report measurements with a time-of-flight viscometer down to [Formula: see text] and up to [Formula: see text], revealing a reduction of viscosity by pressure by as much as 42%. Inspired by a previous attempt [Tanaka H (2000) J Chem Phys 112:799-809], we show that a remarkably simple extension of a two-state model [Holten V, Sengers JV, Anisimov MA (2014) J Phys Chem Ref Data 43:043101], initially developed to reproduce thermodynamic properties, is able to accurately describe dynamic properties (viscosity, self-diffusion coefficient, and rotational correlation time) as well. Our results support the idea that water is a mixture of a high density, "fragile" liquid, and a low density, "strong" liquid, the varying proportion of which explains the anomalies and fragile-to-strong crossover in water.
Collapse
Affiliation(s)
- Lokendra P Singh
- Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumiére Matiére and Institut Universitaire de France, F-69622 Villeurbanne, France
| | - Bruno Issenmann
- Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumiére Matiére and Institut Universitaire de France, F-69622 Villeurbanne, France
| | - Frédéric Caupin
- Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumiére Matiére and Institut Universitaire de France, F-69622 Villeurbanne, France
| |
Collapse
|
41
|
Rovigatti L, Bianco V, Tavares JM, Sciortino F. Communication: Re-entrant limits of stability of the liquid phase and the Speedy scenario in colloidal model systems. J Chem Phys 2017; 146:041103. [DOI: 10.1063/1.4974830] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lorenzo Rovigatti
- Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, Oxford OX1 3NP, United Kingdom
| | - Valentino Bianco
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090 Vienna, Austria
| | - José Maria Tavares
- Instituto Superior de Engenharia de Lisboa-ISEL, Rua Conselheiro Emídio Navarro 1, P-1950-062 Lisbon, Portugal and Centro de Física Teórica e Computacional, Universidade de Lisboa, Campo Grande, P-1749-016 Lisbon, Portugal
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza-Universitá di Roma, Piazzale A. Moro 5, 00185 Roma, Italy and Istituto Sistemi Complessi (CNR-ISC), Via dei Taurini 19, 00185 Roma, Italy
| |
Collapse
|
42
|
Biddle JW, Singh RS, Sparano EM, Ricci F, González MA, Valeriani C, Abascal JLF, Debenedetti PG, Anisimov MA, Caupin F. Two-structure thermodynamics for the TIP4P/2005 model of water covering supercooled and deeply stretched regions. J Chem Phys 2017; 146:034502. [DOI: 10.1063/1.4973546] [Citation(s) in RCA: 86] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- John W. Biddle
- Institute of Physical Science and Technology and Department of
Chemical and Biomolecular Engineering, University of Maryland,
College Park, Maryland 20742, USA
| | - Rakesh S. Singh
- Department of Chemical and Biological Engineering,
Princeton University, Princeton, New Jersey 08544,
USA
| | - Evan M. Sparano
- Department of Chemical and Biological Engineering,
Princeton University, Princeton, New Jersey 08544,
USA
| | - Francesco Ricci
- Department of Chemical and Biological Engineering,
Princeton University, Princeton, New Jersey 08544,
USA
| | - Miguel A. González
- Departamento Química Física I, Facultad Ciencias Químicas,
Universidad Complutense de Madrid, 28040 Madrid,
Spain
| | - Chantal Valeriani
- Departamento Química Física I, Facultad Ciencias Químicas,
Universidad Complutense de Madrid, 28040 Madrid,
Spain
| | - José L. F. Abascal
- Departamento Química Física I, Facultad Ciencias Químicas,
Universidad Complutense de Madrid, 28040 Madrid,
Spain
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering,
Princeton University, Princeton, New Jersey 08544,
USA
| | - Mikhail A. Anisimov
- Institute of Physical Science and Technology and Department of
Chemical and Biomolecular Engineering, University of Maryland,
College Park, Maryland 20742, USA
- Oil and Gas Research Institute of the Russian Academy of Sciences, Moscow 119333, Russia
| | - Frédéric Caupin
- Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne,
France
| |
Collapse
|
43
|
Dhabal D, Chakravarty C, Molinero V, Kashyap HK. Comparison of liquid-state anomalies in Stillinger-Weber models of water, silicon, and germanium. J Chem Phys 2016; 145:214502. [DOI: 10.1063/1.4967939] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Debdas Dhabal
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Charusita Chakravarty
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Hemant K. Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|