1
|
Novak N, Liang X, Kontogeorgis GM. Prediction of water anomalous properties by introducing the two-state theory in SAFT. J Chem Phys 2024; 160:104505. [PMID: 38465683 DOI: 10.1063/5.0186752] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Accepted: 02/14/2024] [Indexed: 03/12/2024] Open
Abstract
Water is one of the most abundant substances on earth, but it is still not entirely understood. It shows unusual behavior, and its properties present characteristic extrema unlike any other fluid. This unusual behavior has been linked to the two-state theory of water, which proposes that water forms different clusters, one with a high density and one with a low density, which may even form two distinct phases at low temperatures. Models incorporating the two-state theory manage to capture the unusual extrema of water, unlike traditional equations of state, which fail. In this work, we have derived the framework to incorporate the two-state theory of water into the Statistical-Associating-Fluid-Theory (SAFT). More specifically, we have assumed that water is an ideal solution of high density water molecules and low density water molecules that are in chemical equilibrium. Using this assumption, we have generalized the association term SAFT to allow for the simultaneous existence of the two water types, which have the same physical parameters but different association properties. We have incorporated the newly derived association term in the context of the Perturbed Chain-SAFT (PC-SAFT). The new model is referred to as PC-SAFT-Two-State (PC-SAFT-TS). Using PC-SAFT-TS, we have succeeded in predicting the characteristic extrema of water, such as its density and speed of sound maximum, etc., without loss of accuracy compared to the original PC-SAFT. This new framework is readily extended to mixtures, and PC-SAFT-TS manages to capture the solubility minimum of hydrocarbons in water in a straightforward manner.
Collapse
Affiliation(s)
- Nefeli Novak
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, 2800 Kgs Lyngby, Denmark
| | - Xiaodong Liang
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, 2800 Kgs Lyngby, Denmark
| | - Georgios M Kontogeorgis
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, 2800 Kgs Lyngby, Denmark
| |
Collapse
|
2
|
Mokshin AV, Vlasov RV. Liquid-Liquid Crossover in Water Model: Local Structure vs Kinetics of Hydrogen Bonds. J Phys Chem B 2024. [PMID: 38411102 DOI: 10.1021/acs.jpcb.3c07650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/28/2024]
Abstract
In equilibrium and supercooled liquids, polymorphism is manifested by thermodynamic regions defined in the phase diagram, which are predominantly of different short- and medium-range order (local structure). It is found that on the phase diagram of the water model, the thermodynamic region corresponding to the equilibrium liquid phase is divided by a line of the smooth liquid-liquid crossover. In the case of the water model TIP4P/2005, this crossover is revealed by various local order parameters and corresponds to pressures on the order of 3150 ± 350 atm at ambient temperature. In the vicinity of the crossover, the dynamics of water molecules change significantly, which is reflected, in particular, in the fact that the self-diffusion coefficient reaches its maximum values. In addition, changes in the structure also manifest themselves in changes in the kinetics of hydrogen bonding, which are captured by values of such quantities as the average lifetime of hydrogen bonding, the average lifetimes of different local coordination numbers, and the frequencies of changes in different local coordination numbers. An interpretation of the hydrogen bond kinetics in terms of the free energy landscape concept in the space of possible coordination numbers is proposed.
Collapse
Affiliation(s)
- Anatolii V Mokshin
- Department of Computational Physics, Kazan (Volga Region) Federal University, Kazan 420008, Russia
| | - Roman V Vlasov
- Department of Computational Physics, Kazan (Volga Region) Federal University, Kazan 420008, Russia
| |
Collapse
|
3
|
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
|
4
|
Muthachikavil AV, Kontogeorgis GM, Liang X, Lei Q, Peng B. Structural characteristics of low-density environments in liquid water. Phys Rev E 2022; 105:034604. [PMID: 35428046 DOI: 10.1103/physreve.105.034604] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 02/24/2022] [Indexed: 06/14/2023]
Abstract
The existence of two structural forms in liquid water has been a point of discussion for a long time. A phase transition between these two forms of liquid water has been proposed based on evidence from molecular simulations, and experiments have also been very recently able to track the proposed transition of the low-density liquid form to the high-density liquid form. We propose to use the average angle an oxygen atom makes with its neighbors to describe the structural environment of a water molecule. The distribution of this order parameter is observed to have two peaks with one peak at ∼109.5^{∘}, corresponding to the internal angle of a regular tetrahedron, indicating tetrahedral arrangement. The other peak corresponds to an environment with a tighter arrangement of neighboring molecules. The distribution of O-O-O angles is decomposed into two skewed distributions to estimate the fractions of the two liquid forms in water. A good similarity is observed between the temperature and pressure trends of fractions of locally favored tetrahedral structure (LFTS) form estimated using the new order parameter and the reports in the literature, over a range of temperatures and pressures. We also compare the structural environments indicated by different order parameters and find that the order parameter proposed in this paper captures the structure of first solvation shell of the LFTS accurately.
Collapse
Affiliation(s)
- Aswin V Muthachikavil
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| | - Georgios M Kontogeorgis
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| | - Xiaodong Liang
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Kongens Lyngby 2800, Denmark
| | - Qun Lei
- Research Institute of Petroleum Exploration and Development (RIPED), PetroChina, Beijing 100083, China
| | - Baoliang Peng
- Research Institute of Petroleum Exploration and Development (RIPED), PetroChina, Beijing 100083, China
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Dubey V, Dueby S, Daschakraborty S. Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective. Phys Chem Chem Phys 2021; 23:19964-19986. [PMID: 34515269 DOI: 10.1039/d1cp02202d] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.
Collapse
Affiliation(s)
- Vikas Dubey
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | - Shivam Dueby
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | | |
Collapse
|
8
|
Parui S, Jana B. Cold denaturation induced helix-to-helix transition and its implication to activity of helical antifreeze protein. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116627] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
9
|
Chalyi AV. Synergetic dialogue “physics – medicine”: Hexagons in living and inanimate nature. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
10
|
Horstmann R, Vogel M. Relations between thermodynamics, structures, and dynamics for modified water models in their supercooled regimes. J Chem Phys 2021; 154:054502. [DOI: 10.1063/5.0037080] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- R. Horstmann
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - M. Vogel
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| |
Collapse
|
11
|
Liquid-liquid transition as a perturbation on the van der Waals' equation. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
12
|
Chalyi AV. Synergetic dialogue “physics – medicine”: Hexagons in living and inanimate nature. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
13
|
Cheng YH, Yang HC, Chou PT. Could Chemical Reaction at the Molecular Level Show Distinction between Two Liquid-Water States? Study of the Excited-State Water-Catalyzed Proton Transfer Reaction Provides a Clue. J Phys Chem Lett 2020; 11:9468-9475. [PMID: 33108192 DOI: 10.1021/acs.jpclett.0c02896] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The two liquid-water states, which lead to some anomalies when temperature crosses over 50 ± 10 °C at the atmospheric pressure, have been continuously catching popular attention. In this study, using the excited-state proton transfer (ESPT) catalyzed by water molecules as a prototypical reaction, we demonstrate that the kinetics of ESPT indeed is influenced by the two liquid-water states. In the water-catalyzed ESPT of 3-cyano-7-azaindole (3CAI), a repetitive and comprehensive temperature-dependent study of ESPT in H2O from 0 to 90 °C shows anomalous behavior. The plot of the logarithm of ESPT rate constant as a function of inverse of absolute temperature deviates from a straight line. The convex-Arrhenius behavior manifests the activation free energy for water-assisted ESPT being dependent on temperature and hence the liquid water structure. To simplify the discussion, the plot is well fitted by using two straight lines that are crossed over in the vicinity of 40 °C. The free energy difference between water-solvated 3CAI and the 1:1 H2O:3CAI complex is deduced to be 2.29 ± 0.04 and 1.96 ± 0.04 kcal·mol-1 in the regions of 0-40 and 40-90 °C water, respectively, which also results in different frequency factors, i.e., the proton transfer/tunneling rates of (5.83 ± 0.36) × 1010 and (3.48 ± 0.27) × 1010 s-1, respectively. In a qualitative manner, the results are then rationalized by the different types of H-bonding configuration as proposed for two liquid-water phases, rendering experimental evidence to support the different water phases in ambient temperatures at 1 bar.
Collapse
Affiliation(s)
- Yu-Hsuan Cheng
- Department of Chemistry, National Taiwan University, Taipei 10617, Taiwan
| | - Hsiao-Ching Yang
- Department of Chemistry, Fu Jen Catholic University, New Taipei City 24205, Taiwan
| | - Pi-Tai Chou
- Department of Chemistry, National Taiwan University, Taipei 10617, Taiwan
| |
Collapse
|
14
|
Affiliation(s)
- Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
15
|
|
16
|
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
|
17
|
Abstract
The idea that water is a mixture of two distinct states is analyzed in some detail. It is shown that the known compressibility of water is in fact sufficiently small that for a volume of water of size 1 nm3, the density fluctuations are of order 4% of the average density. This is much smaller than the ≈25% density fluctuations that would be required for significant regions of high and low density water to occur on this volume scale. It is also pointed out that the density fluctuations in water are, if anything, smaller than those that occur in other common liquids which do not have the anomalous properties of water. It is shown that if the distribution of density fluctuations is unimodal, the system is in the one-phase region, and if bimodal, it is in the two-phase region. None of the liquid or amorphous phases of water explored in this work give any sign of being in the two-phase region. Existing neutron and X-ray scattering data on water in the amorphous phases, and in the stable liquid phases as a function pressure and temperature, are subject to a new set of empirical potential structure refinement simulations. These simulations are interrogated for their configurational entropy, using a spherical harmonic reconstruction of the full orientational pair correlation function. It is shown that the excess pair entropy derived from this function, plus the known perfect gas contributions, give a reasonable account of the total entropy of water, within the likely errors. This estimated entropy follows the expected declining trend with decreasing temperature. Evidence that higher density water will have higher entropy than lower density water emerges, in accordance with what is expected from the negative thermal expansion coefficient of water at low temperatures. However, this entropy increase is not large and goes through a maximum before declining at yet higher densities and pressures, in a manner reminiscent of what has been previously observed in the diffusion coefficient as a function of pressure. There is no evidence that ambient water can be regarded as patches of high density, high entropy and low density, low entropy liquid, as some have claimed, since high density water has a similar entropy to low density water. There is some evidence that the distinction between these two states will become more pronounced as the temperature is lowered. Extensive discussion of the use of order parameters to describe water structure is given, and it is pointed out that these indices generally cannot be used to infer two-state behavior.
Collapse
Affiliation(s)
- A K Soper
- ISIS Facility, UKRI-STFC Rutherford Appleton Laboratory, Harwell Campus, Didcot, Oxfordshire OX11 0QX, United Kingdom
| |
Collapse
|
18
|
Martelli F. Unravelling the contribution of local structures to the anomalies of water: The synergistic action of several factors. J Chem Phys 2019; 150:094506. [PMID: 30849899 DOI: 10.1063/1.5087471] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the microscopic origin of water's anomalies by inspecting the hydrogen bond network (HBN) and the spatial organization of low-density-liquid (LDL) like and high-density-liquid (HDL) like environments. Specifically, we simulate-via classical molecular dynamics simulations-the isobaric cooling of a sample composed of 512 water molecules from ambient to deeply undercooled conditions at three pressures, namely, 1 bar, 400 bars, and 1000 bars. In correspondence with the Widom line (WL), (i) the HDL-like dominating cluster undergoes fragmentation caused by the percolation of LDL-like aggregates following a spinodal-like kinetics; (ii) such fragmentation always occurs at a "critical" concentration of ∼20%-30% in LDL; (iii) the HBN within LDL-like environments is characterized by an equal number of pentagonal and hexagonal rings that create a state of maximal frustration between a configuration that promotes crystallization (hexagonal ring) and a configuration that hinders it (pentagonal ring); (iv) the spatial organization of HDL-like environments shows a marked variation. Moreover, the inspection of the global symmetry shows that the intermediate-range order decreases in correspondence with the WL and such a decrease becomes more pronounced upon increasing the pressure, hence supporting the hypothesis of a liquid-liquid critical point. Our results reveal and rationalize the complex microscopic origin of water's anomalies as the cooperative effect of several factors acting synergistically. Beyond implications for water, our findings may be extended to other materials displaying anomalous behaviours.
Collapse
Affiliation(s)
- Fausto Martelli
- IBM Research, Hartree Centre, Daresbury WA4 4AD, United Kingdom
| |
Collapse
|
19
|
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
|
20
|
Hestand NJ, Skinner JL. Perspective: Crossing the Widom line in no man’s land: Experiments, simulations, and the location of the liquid-liquid critical point in supercooled water. J Chem Phys 2018; 149:140901. [DOI: 10.1063/1.5046687] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nicholas J. Hestand
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - J. L. Skinner
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
21
|
Zhao G, Wang H, Hu DM, Ding MC, Zhao XG, Yan JL. Anomalous phase behavior of first-order fluid-liquid phase transition in phosphorus. J Chem Phys 2018; 147:204501. [PMID: 29195280 DOI: 10.1063/1.4999009] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Although the existence of liquid-liquid phase transition has become more and more convincing, whether it will terminate at a critical point and what is the order parameter are still open. To explore these questions, we revisit the fluid-liquid phase transition (FLPT) in phosphorus (P) and study its phase behavior by performing extensive first-principles molecular dynamics simulations. The FLPT observed in experiments is well reproduced, and a fluid-liquid critical point (FLCP) at T = 3000 ∼ 3500 K, P = 1.5-2.0 Kbar is found. With decreasing temperature from the FLCP along the transition line, the density difference (Δρ) between two coexisting phases first increases from zero and then anomalously decreases; however, the entropy difference (ΔS) continuously increases from zero. These features suggest that an order parameter containing contributions from both the density and the entropy is needed to describe the FLPT in P, and at least at low temperatures, the entropy, instead of the density, governs the FLPT.
Collapse
Affiliation(s)
- G Zhao
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - H Wang
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - D M Hu
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - M C Ding
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - X G Zhao
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - J L Yan
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| |
Collapse
|
22
|
Ni Y, Hestand NJ, Skinner JL. Communication: Diffusion constant in supercooled water as the Widom line is crossed in no man’s land. J Chem Phys 2018; 148:191102. [DOI: 10.1063/1.5029822] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yicun Ni
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Nicholas J. Hestand
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - J. L. Skinner
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
23
|
Anomalous propagation and scattering of sound in 2-propanol water solution near its singular point. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.01.040] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
24
|
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
|
25
|
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
|
26
|
Galamba N. On the hydrogen-bond network and the non-Arrhenius transport properties of water. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:015101. [PMID: 27831934 DOI: 10.1088/0953-8984/29/1/015101] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We study the structural and dynamic transformations of SPC/E water with temperature, through molecular dynamics (MD), and discuss the non-Arrhenius behavior of the transport properties and orientational dynamics, and the magnitude of the breakdown of the Stokes-Einstein (SE) and the Stokes-Einstein-Debye (SED) relations, in the light of these transformations. Our results show that deviations from Arrhenius behavior of the self-diffusion at low temperatures cannot be exclusively explained by the reduction of water defects (interstitial waters) and the increase of the local tetrahedrality, thus, suggesting the importance of the slowdown of collective rearrangements. Interestingly we find that at high temperatures (T ⩾ 340 K) water defects lead to a slight increase of the tetrahedrality and a decrease of the self-diffusion, opposite to water at low temperatures. The relative magnitude of the breakdown of the SE and the SED relations is found to be in accord with recent experiments (Dehaoui et al 2015 Proc. Natl Acad. Sci. USA 112 12020) resolving the discrepancy with previous MD results. Further, we show that SPC/E hydrogen-bond (HB) lifetimes deviate from Arrhenious behaviour at low temperatures in contrast with some previous MD studies. This deviation is nevertheless much smaller than that observed for the orientational dynamics and the transport properties of water, consistent with the relaxation times measured by several experimental methods. The HB acceptor exchange dynamics defined here by the acceptor switch and reform (librational dynamics) frequencies exhibit similar Arrhenius deviations, thus explaining to some extent the non-Arrhenius behavior of the transport properties and of the orientational dynamics of water. Our results also show that the fraction of HB switches through a bifurcated pathway follow a power law with the temperature decrease. Thus, at low temperatures HB acceptor switches are less frequent but occur on a faster time scale consistent with the temperature dependence of the ratio of the rotational relaxation times for the different Legendre polynomial ranks.
Collapse
Affiliation(s)
- N Galamba
- Centro de Química e Bioquímica da Faculdade de Ciências da Universidade de Lisboa, Edifício C8, Campo Grande, 1749-016 Lisboa, Portugal
| |
Collapse
|
27
|
Urbic T. Liquid-liquid critical point in a simple analytical model of water. Phys Rev E 2016; 94:042126. [PMID: 27841542 DOI: 10.1103/physreve.94.042126] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2016] [Indexed: 11/07/2022]
Abstract
A statistical model for a simple three-dimensional Mercedes-Benz model of water was used to study phase diagrams. This model on a simple level describes the thermal and volumetric properties of waterlike molecules. A molecule is presented as a soft sphere with four directions in which hydrogen bonds can be formed. Two neighboring waters can interact through a van der Waals interaction or an orientation-dependent hydrogen-bonding interaction. For pure water, we explored properties such as molar volume, density, heat capacity, thermal expansion coefficient, and isothermal compressibility and found that the volumetric and thermal properties follow the same trends with temperature as in real water and are in good general agreement with Monte Carlo simulations. The model exhibits also two critical points for liquid-gas transition and transition between low-density and high-density fluid. Coexistence curves and a Widom line for the maximum and minimum in thermal expansion coefficient divides the phase space of the model into three parts: in one part we have gas region, in the second a high-density liquid, and the third region contains low-density liquid.
Collapse
Affiliation(s)
- Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, 1000 Lubljana, Slovenia
| |
Collapse
|
28
|
Gallo P, Amann-Winkel K, Angell CA, Anisimov MA, Caupin F, Chakravarty C, Lascaris E, Loerting T, Panagiotopoulos AZ, Russo J, Sellberg JA, Stanley HE, Tanaka H, Vega C, Xu L, Pettersson LGM. Water: A Tale of Two Liquids. Chem Rev 2016; 116:7463-500. [PMID: 27380438 PMCID: PMC5424717 DOI: 10.1021/acs.chemrev.5b00750] [Citation(s) in RCA: 421] [Impact Index Per Article: 52.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Water is the most abundant liquid
on earth and also the substance
with the largest number of anomalies in its properties. It is a prerequisite
for life and as such a most important subject of current research
in chemical physics and physical chemistry. In spite of its simplicity
as a liquid, it has an enormously rich phase diagram where different
types of ices, amorphous phases, and anomalies disclose a path that
points to unique thermodynamics of its supercooled liquid state that
still hides many unraveled secrets. In this review we describe the
behavior of water in the regime from ambient conditions to the deeply
supercooled region. The review describes simulations and experiments
on this anomalous liquid. Several scenarios have been proposed to
explain the anomalous properties that become strongly enhanced in
the supercooled region. Among those, the second critical-point scenario
has been investigated extensively, and at present most experimental
evidence point to this scenario. Starting from very low temperatures,
a coexistence line between a high-density amorphous phase and a low-density
amorphous phase would continue in a coexistence line between a high-density
and a low-density liquid phase terminating in a liquid–liquid
critical point, LLCP. On approaching this LLCP from the one-phase
region, a crossover in thermodynamics and dynamics can be found. This
is discussed based on a picture of a temperature-dependent balance
between a high-density liquid and a low-density liquid favored by,
respectively, entropy and enthalpy, leading to a consistent picture
of the thermodynamics of bulk water. Ice nucleation is also discussed,
since this is what severely impedes experimental investigation of
the vicinity of the proposed LLCP. Experimental investigation of stretched
water, i.e., water at negative pressure, gives access to a different
regime of the complex water diagram. Different ways to inhibit crystallization
through confinement and aqueous solutions are discussed through results
from experiments and simulations using the most sophisticated and
advanced techniques. These findings represent tiles of a global picture
that still needs to be completed. Some of the possible experimental
lines of research that are essential to complete this picture are
explored.
Collapse
Affiliation(s)
- Paola Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre , Via della Vasca Navale 84, 00146 Rome, Italy
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Charles Austen Angell
- Department of Chemistry and Biochemistry, Arizona State University , Tempe, Arizona 85287, United States
| | - Mikhail Alexeevich Anisimov
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland , College Park, Maryland 20742, United States
| | - Frédéric Caupin
- Institut Lumière Matière, UMR5306 Université Claude Bernard Lyon 1-CNRS, Université de Lyon, Institut Universitaire de France , 69622 Villeurbanne, France
| | - Charusita Chakravarty
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas , New Delhi 110016, India
| | - Erik Lascaris
- Center for Polymer Studies and Department of Physics, Boston University , Boston, Massachusetts 02215, United States
| | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck , 6020 Innsbruck, Austria
| | | | - John Russo
- Institute of Industrial Science, University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan.,School of Mathematics, University of Bristol , Bristol BS8 1TW, United Kingdom
| | - Jonas Alexander Sellberg
- Biomedical and X-ray Physics, Department of Applied Physics, AlbaNova University Center, KTH Royal Institute of Technology , SE-106 91 Stockholm, Sweden
| | - Harry Eugene Stanley
- Center for Polymer Studies and Department of Physics, Boston University , Boston, Massachusetts 02215, United States
| | - Hajime Tanaka
- Institute of Industrial Science, University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Carlos Vega
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid , 28040 Madrid, Spain
| | - Limei Xu
- International Centre for Quantum Materials and School of Physics, Peking University , Beijing 100871, China.,Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
| | | |
Collapse
|
29
|
Singh RS, Biddle JW, Debenedetti PG, Anisimov MA. Two-state thermodynamics and the possibility of a liquid-liquid phase transition in supercooled TIP4P/2005 water. J Chem Phys 2016; 144:144504. [DOI: 10.1063/1.4944986] [Citation(s) in RCA: 121] [Impact Index Per Article: 15.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Rakesh S. Singh
- Department of Chemical & Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - John W. Biddle
- Institute of Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Pablo G. Debenedetti
- Department of Chemical & 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
| |
Collapse
|
30
|
Lascaris E. Tunable Liquid-Liquid Critical Point in an Ionic Model of Silica. PHYSICAL REVIEW LETTERS 2016; 116:125701. [PMID: 27058086 DOI: 10.1103/physrevlett.116.125701] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/08/2015] [Indexed: 06/05/2023]
Abstract
Recently, it was shown that the Woodcock-Angell-Cheeseman model for liquid silica [L. V. Woodcock, C. A. Angell, and P. Cheeseman, J. Chem. Phys. 65, 1565 (1976)] is remarkably close to having a liquid-liquid critical point (LLCP). We demonstrate that increasing the ion charge separates the global maxima of the response functions, while reducing the charge smoothly merges them into a LLCP, a phenomenon that might be experimentally observable with charged colloids. An analysis of the Si and O coordination numbers suggests that a sufficiently low Si/O coordination number ratio is needed to attain a LLCP.
Collapse
Affiliation(s)
- Erik Lascaris
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
| |
Collapse
|
31
|
Luo J, Xu L, Angell CA, Stanley HE, Buldyrev SV. Physics of the Jagla model as the liquid-liquid coexistence line slope varies. J Chem Phys 2015; 142:224501. [DOI: 10.1063/1.4921559] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jiayuan Luo
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
| | - Limei Xu
- International Center for Quantum Materials, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Quantum Matter, Beijing, China
| | - C. Austen Angell
- Department of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287, USA
| | - H. Eugene Stanley
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
| | - Sergey V. Buldyrev
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
- Department of Physics, Yeshiva University, 500 West 185th Street, New York, New York 10033, USA
| |
Collapse
|
32
|
Huš M, Urbic T. Existence of a liquid-liquid phase transition in methanol. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:062306. [PMID: 25615092 DOI: 10.1103/physreve.90.062306] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/26/2014] [Indexed: 06/04/2023]
Abstract
A simple model is constructed to study the phase diagram and thermodynamic properties of methanol, which is described as a dimer of an apolar sphere mimicking the methyl group and a sphere with core-softened potential as the hydroxyl group. Performing classical Monte Carlo simulations, we obtained the phase diagram, showing a second critical point between two different liquid phases. Evaluating systems with a different number of particles, we extrapolate to infinite size in accordance with Ising universality class to obtain bulk values for critical temperature, pressure, and density. Strong evidence that the structure of the liquid changes upon transition from high- to low-density phase was provided. From the experimentally determined hydrogen bond strength and length in methanol and water, we propose where the second critical point of methanol should be.
Collapse
Affiliation(s)
- Matej Huš
- Chair of Physical Chemistry, Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Chair of Physical Chemistry, Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
33
|
Holten V, Palmer JC, Poole PH, Debenedetti PG, Anisimov MA. Two-state thermodynamics of the ST2 model for supercooled water. J Chem Phys 2014; 140:104502. [PMID: 24628177 DOI: 10.1063/1.4867287] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Thermodynamic properties of the ST2 model for supercooled liquid water exhibit anomalies similar to those observed in real water. A possible explanation of these anomalies is the existence of a metastable, liquid-liquid transition terminated by a critical point. This phenomenon, whose possible existence in real water is the subject of much current experimental work, has been unambiguously demonstrated for this particular model by most recent simulations. In this work, we reproduce the anomalies of two versions of the ST2 model with an equation of state describing water as a non-ideal "mixture" of two different types of local molecular order. We show that the liquid-liquid transition in the ST2 water is energy-driven. This is in contrast to another popular model, mW, in which non-ideality in mixing of two alternative local molecular orders is entropy-driven, and is not sufficiently strong to induce a liquid-liquid transition.
Collapse
Affiliation(s)
- Vincent Holten
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Jeremy C Palmer
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Peter H Poole
- Department of Physics, St. Francis Xavier University, Antigonish, Nova Scotia B2G 2W5, Canada
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Mikhail A Anisimov
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
34
|
Biddle JW, Holten V, Anisimov MA. Behavior of supercooled aqueous solutions stemming from hidden liquid–liquid transition in water. J Chem Phys 2014; 141:074504. [DOI: 10.1063/1.4892972] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- John W. Biddle
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Vincent Holten
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Mikhail A. Anisimov
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
35
|
Luo J, Xu L, Lascaris E, Stanley HE, Buldyrev SV. Behavior of the Widom line in critical phenomena. PHYSICAL REVIEW LETTERS 2014; 112:135701. [PMID: 24745439 DOI: 10.1103/physrevlett.112.135701] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/06/2013] [Indexed: 05/15/2023]
Abstract
Using linear scaling theory, we study the behavior of response functions extrema in the vicinity of the critical point. We investigate how the speed of convergence of the loci of response function extrema to the Widom line depends on the parameters of the linear scaling theory. We find that when the slope of the coexistence line is near zero, the line of specific heat maxima does not follow the Widom line but instead follows the coexistence line. This has relevance for the detection of liquid-liquid critical points, which can exhibit a near-horizontal coexistence line. Our theoretical predictions are confirmed by computer simulations of a family of spherically symmetric potentials.
Collapse
Affiliation(s)
- Jiayuan Luo
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
| | - Limei Xu
- International Center for Quantum Materials, Peking University, Beijing 100871, China and Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
| | - Erik Lascaris
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
| | - H Eugene Stanley
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
| | - Sergey V Buldyrev
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA and Department of Physics, Yeshiva University, 500 West 185th Street, New York, New York 10033, USA
| |
Collapse
|
36
|
Bianco V, Franzese G. Critical behavior of a water monolayer under hydrophobic confinement. Sci Rep 2014; 4:4440. [PMID: 24699181 PMCID: PMC3975237 DOI: 10.1038/srep04440] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2013] [Accepted: 03/06/2014] [Indexed: 12/16/2022] Open
Abstract
The properties of water can have a strong dependence on the confinement. Here, we consider a water monolayer nanoconfined between hydrophobic parallel walls under conditions that prevent its crystallization. We investigate, by simulations of a many-body coarse-grained water model, how the properties of the liquid are affected by the confinement. We show, by studying the response functions and the correlation length and by performing finite-size scaling of the appropriate order parameter, that at low temperature the monolayer undergoes a liquid-liquid phase transition ending in a critical point in the universality class of the two-dimensional (2D) Ising model. Surprisingly, by reducing the linear size L of the walls, keeping the walls separation h constant, we find a 2D-3D crossover for the universality class of the liquid-liquid critical point for L=h ≃ 50, i.e. for a monolayer thickness that is small compared to its extension. This result is drastically different from what is reported for simple liquids, where the crossover occurs for L=h ≃ 5, and is consistent with experimental results and atomistic simulations. We shed light on these findings showing that they are a consequence of the strong cooperativity and the low coordination number of the hydrogen bond network that characterizes water.
Collapse
Affiliation(s)
- Valentino Bianco
- Departament de Física Fonamental, Universitat de Barcelona, Martí i Franquès 1, ES-08028 Barcelona, Spain
| | - Giancarlo Franzese
- Departament de Física Fonamental, Universitat de Barcelona, Martí i Franquès 1, ES-08028 Barcelona, Spain
| |
Collapse
|
37
|
Overduin SD, Patey GN. An analysis of fluctuations in supercooled TIP4P/2005 water. J Chem Phys 2013; 138:184502. [PMID: 23676051 DOI: 10.1063/1.4803868] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Large-scale, long-time molecular dynamics simulations are used to investigate fluctuations in the TIP4P/2005 water model in the supercooled region (240-190 K). Particular attention is focused in the vicinity of a previously reported liquid-liquid critical point [J. L. F. Abascal and C. Vega, J. Chem. Phys. 133, 234502 (2010)]. Water is viewed as an equimolar binary mixture with "species" defined based on a local tetrahedral order parameter. A Bhatia-Thornton fluctuation analysis is used to show that species concentration fluctuations couple to density fluctuations and completely account for the anomalous increase in the structure factor at small wave number observed under supercooled conditions. Although we find that both concentration and density fluctuations increase with decreasing temperature along the proposed critical isochore, we cannot confirm the existence of a liquid-liquid critical point. Our simulations suggest that the parameters previously reported are not a true liquid-liquid critical point and we find no evidence of two-phase coexistence in its vicinity. It is shown that very long simulations (on the order of 8 μs for 500 molecules) are necessary to obtain well converged density distributions for deeply supercooled water and this is especially important if one is seeking direct evidence of a two-phase region.
Collapse
Affiliation(s)
- S D Overduin
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | | |
Collapse
|
38
|
Soper AK. Radical re-appraisal of water structure in hydrophilic confinement. Chem Phys Lett 2013; 590:1-15. [PMID: 25843963 PMCID: PMC4376068 DOI: 10.1016/j.cplett.2013.10.075] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2013] [Accepted: 10/22/2013] [Indexed: 11/29/2022]
Abstract
The structure of water confined in MCM41 silica cylindrical pores is studied to determine whether confined water is simply a version of the bulk liquid which can be substantially supercooled without crystallisation. A combination of total neutron scattering from the porous silica, both wet and dry, and computer simulation using a realistic model of the scattering substrate is used. The water in the pore is divided into three regions: core, interfacial and overlap. The average local densities of water in these simulations are found to be about 20% lower than bulk water density, while the density in the core region is below, but closer to, the bulk density. There is a decrease in both local and core densities when the temperature is lowered from 298 K to 210 K. The radical proposal is made here that water in hydrophilic confinement is under significant tension, around -100 MPa, inside the pore.
Collapse
|
39
|
Nagoe A, Oguni M. Hydroxylamine-doping effect on the Tg of 160 K for water confined in silica-gel nanopores. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2013; 25:465110. [PMID: 24162332 DOI: 10.1088/0953-8984/25/46/465110] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
The glass transition behavior of hydroxylamine (HA) aqueous solutions in bulk and confined in silica-gel nanopores with average width of 1.1 nm was studied by means of differential scanning calorimetry measurements and adiabatic calorimetry. The glass transition temperature (Tg) of the confined solution with high HA mole-fraction (xHA) was essentially the same as the value of the bulk. This suggests that the nano-size confinement affects the Tg of HA aqueous solution little. Meanwhile, the bulk solution with xHA < 0.3 revealed partial crystallization on cooling and, on the other hand, the confined solution with the same xHA did not crystallize. The Tg of the xHA = 0.076 confined solution was 174 K which is higher than the value of 160 K for pure water confined in the same silica-gel pores. This demonstrates that HA doping leads to no abrupt Tg-decrease, unlike doping of all the other second components reported so far, suggesting that HA is set neatly in a hydrogen-bond network formed by water molecules. We discuss the xHA dependence of Tg for the HA aqueous solutions from a viewpoint related to peculiar phase-behavior of pure water. Considering that the xHA = 0.076 aqueous solution revealed no anomaly compared with pure water, it was recognized as corresponding to the high-temperature phase of pure water.
Collapse
|
40
|
Holten V, Limmer DT, Molinero V, Anisimov MA. Nature of the anomalies in the supercooled liquid state of the mW model of water. J Chem Phys 2013; 138:174501. [PMID: 23656138 DOI: 10.1063/1.4802992] [Citation(s) in RCA: 90] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The thermodynamic properties of the supercooled liquid state of the mW model of water show anomalous behavior. Like in real water, the heat capacity and compressibility sharply increase upon supercooling. One of the possible explanations of these anomalies, the existence of a second (liquid-liquid) critical point, is not supported by simulations for this model. In this work, we reproduce the anomalies of the mW model with two thermodynamic scenarios: one based on a non-ideal "mixture" with two different types of local order of the water molecules, and one based on weak crystallization theory. We show that both descriptions accurately reproduce the model's basic thermodynamic properties. However, the coupling constant required for the power laws implied by weak crystallization theory is too large relative to the regular backgrounds, contradicting assumptions of weak crystallization theory. Fluctuation corrections outside the scope of this work would be necessary to fit the forms predicted by weak crystallization theory. For the two-state approach, the direct computation of the low-density fraction of molecules in the mW model is in agreement with the prediction of the phenomenological equation of state. The non-ideality of the "mixture" of the two states never becomes strong enough to cause liquid-liquid phase separation, also in agreement with simulation results.
Collapse
Affiliation(s)
- Vincent Holten
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | | | | | | |
Collapse
|
41
|
|
42
|
Biddle JW, Holten V, Sengers JV, Anisimov MA. Thermal conductivity of supercooled water. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:042302. [PMID: 23679409 DOI: 10.1103/physreve.87.042302] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/25/2013] [Indexed: 06/02/2023]
Abstract
The heat capacity of supercooled water, measured down to -37°C, shows an anomalous increase as temperature decreases. The thermal diffusivity, i.e., the ratio of the thermal conductivity and the heat capacity per unit volume, shows a decrease. These anomalies may be associated with a hypothesized liquid-liquid critical point in supercooled water below the line of homogeneous nucleation. However, while the thermal conductivity is known to diverge at the vapor-liquid critical point due to critical density fluctuations, the thermal conductivity of supercooled water, calculated as the product of thermal diffusivity and heat capacity, does not show any sign of such an anomaly. We have used mode-coupling theory to investigate the possible effect of critical fluctuations on the thermal conductivity of supercooled water and found that indeed any critical thermal-conductivity enhancement would be too small to be measurable at experimentally accessible temperatures. Moreover, the behavior of thermal conductivity can be explained by the observed anomalies of the thermodynamic properties. In particular, we show that thermal conductivity should go through a minimum when temperature is decreased, as Kumar and Stanley observed in the TIP5P model of water. We discuss physical reasons for the striking difference between the behavior of thermal conductivity in water near the vapor-liquid and liquid-liquid critical points.
Collapse
Affiliation(s)
- John W Biddle
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | | | | | | |
Collapse
|
43
|
Anisimov MA. Cold and supercooled water: A novel supercritical-fluid solvent. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY B 2013. [DOI: 10.1134/s199079311208009x] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
44
|
Bertrand CE, Zhang Y, Chen SH. Deeply-cooled water under strong confinement: neutron scattering investigations and the liquid–liquid critical point hypothesis. Phys Chem Chem Phys 2013. [DOI: 10.1039/c2cp43235h] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
45
|
Liu Y, Palmer JC, Panagiotopoulos AZ, Debenedetti PG. Liquid-liquid transition in ST2 water. J Chem Phys 2012; 137:214505. [DOI: 10.1063/1.4769126] [Citation(s) in RCA: 128] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
46
|
Mazza MG, Stokely K, Stanley HE, Franzese G. Effect of pressure on the anomalous response functions of a confined water monolayer at low temperature. J Chem Phys 2012. [DOI: 10.1063/1.4767355] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
47
|
Holten V, Anisimov MA. Entropy-driven liquid-liquid separation in supercooled water. Sci Rep 2012; 2:713. [PMID: 23056905 PMCID: PMC3465811 DOI: 10.1038/srep00713] [Citation(s) in RCA: 176] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2012] [Accepted: 09/19/2012] [Indexed: 11/09/2022] Open
Abstract
Twenty years ago Poole et al. suggested that the anomalous properties of supercooled water may be caused by a critical point that terminates a line of liquid–liquid separation of lower-density and higher-density water. Here we present a thermodynamic model based on this hypothesis, which describes all available experimental data for supercooled water with better quality and fewer adjustable parameters than any other model. Liquid water at low temperatures is viewed as an ‘athermal solution' of two molecular structures with different entropies and densities. Alternatively to popular models for water, in which liquid–liquid separation is driven by energy, the phase separation in the athermal two-state water is driven by entropy upon increasing the pressure, while the critical temperature is defined by the ‘reaction' equilibrium constant. The model predicts the location of density maxima at the locus of a near-constant fraction of the lower-density structure.
Collapse
Affiliation(s)
- V Holten
- Institute for Physical Science & Technology, University of Maryland, College Park, Maryland 20742, USA
| | | |
Collapse
|
48
|
Holten V, Bertrand CE, Anisimov MA, Sengers JV. Thermodynamics of supercooled water. J Chem Phys 2012; 136:094507. [PMID: 22401452 DOI: 10.1063/1.3690497] [Citation(s) in RCA: 174] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We review the available experimental information on the thermodynamic properties of supercooled water and demonstrate the possibility of modeling these thermodynamic properties on a theoretical basis. We show that by assuming the existence of a liquid-liquid critical point in supercooled water, the theory of critical phenomena can give an accurate account of the experimental thermodynamic-property data up to a pressure of 150 MPa. In addition, we show that a phenomenological extension of the theoretical model can account for all currently available experimental data in the supercooled region, up to 400 MPa. The stability limit of the liquid state and possible coupling between crystallization and liquid-liquid separation are also discussed. It is concluded that critical-point thermodynamics describes the available thermodynamic data for supercooled water within experimental accuracy, thus establishing a benchmark for further developments in this area.
Collapse
Affiliation(s)
- V Holten
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | | | | | | |
Collapse
|
49
|
Nanoscale dynamics of phase flipping in water near its hypothesized liquid-liquid critical point. Sci Rep 2012; 2:474. [PMID: 22761987 PMCID: PMC3386518 DOI: 10.1038/srep00474] [Citation(s) in RCA: 96] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2012] [Accepted: 05/30/2012] [Indexed: 11/25/2022] Open
Abstract
One hypothesized explanation for water's anomalies imagines the existence of a liquid-liquid (LL) phase transition line separating two liquid phases and terminating at a LL critical point. We simulate the classic ST2 model of water for times up to 1000 ns and system size up to N = 729. We find that for state points near the LL transition line, the entire system flips rapidly between liquid states of high and low density. Our finite-size scaling analysis accurately locates both the LL transition line and its associated LL critical point. We test the stability of the two liquids with respect to the crystal and find that of the 350 systems simulated, only 3 of them crystallize and these 3 for the relatively small system size N = 343 while for all other simulations the incipient crystallites vanish on a time scales smaller than ≈ 100 ns.
Collapse
|
50
|
Suzuki H, Matsuzaki Y, Muraoka A, Tachikawa M. Raman spectroscopy of optically levitated supercooled water droplet. J Chem Phys 2012; 136:234508. [DOI: 10.1063/1.4729476] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|