1
|
Cerdeiriña CA, Troncoso J. Ising Paradigm in Isobaric Ensembles. ENTROPY (BASEL, SWITZERLAND) 2024; 26:438. [PMID: 38920447 PMCID: PMC11203012 DOI: 10.3390/e26060438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Revised: 05/17/2024] [Accepted: 05/20/2024] [Indexed: 06/27/2024]
Abstract
We review recent work on Ising-like models with "compressible cells" of fluctuating volume that, as such, are naturally treated in NpT and μpT ensembles. Besides volumetric phenomena, local entropic effects crucially underlie the models. We focus on "compressible cell gases" (CCG), namely, lattice gases with fluctuating cell volumes, and "compressible cell liquids" (CCL) with singly occupied cells and fluctuating cell volumes. CCGs contemplate singular diameters and "Yang-Yang features" predicted by the "complete scaling" formulation of asymmetric fluid criticality, with a specific version incorporating "ice-like" hydrogen bonding further describing the "singularity-free scenario" for the low-temperature unusual thermodynamics of supercooled water. In turn, suitable CCL variants constitute adequate prototypes of water-like liquid-liquid criticality and the freezing transition of a system of hard spheres. On incorporating vacant cells to such two-state CCL variants, one obtains three-state, BEG-like models providing a satisfactory description of water's "second-critical-point scenario" and the whole phase behavior of a simple substance like argon. Future challenges comprise water's crystal-fluid phase behavior and metastable states.
Collapse
|
2
|
Kimmel GA. Isotope effects in supercooled H2O and D2O and a corresponding-states-like rescaling of the temperature and pressure. J Chem Phys 2024; 160:194508. [PMID: 38767262 DOI: 10.1063/5.0207719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2024] [Accepted: 04/30/2024] [Indexed: 05/22/2024] Open
Abstract
Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H2O and D2O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations-with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a "corresponding-states-like" rescaling for water properties in three classical water models that all exhibit a liquid-liquid transition and critical point [Uralcan et al., J. Chem. Phys. 150, 064503 (2019)], the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H2O and D2O is investigated here. Utilizing previously published data and equations-of-state for H2O and D2O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid-liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water's anomalies.
Collapse
Affiliation(s)
- Greg A Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| |
Collapse
|
3
|
Saito S. Unraveling the dynamic slowdown in supercooled water: The role of dynamic disorder in jump motions. J Chem Phys 2024; 160:194506. [PMID: 38767263 DOI: 10.1063/5.0209713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Accepted: 05/01/2024] [Indexed: 05/22/2024] Open
Abstract
When a liquid is rapidly cooled below its melting point without inducing crystallization, its dynamics slow down significantly without noticeable structural changes. Elucidating the origin of this slowdown has been a long-standing challenge. Here, we report a theoretical investigation into the mechanism of the dynamic slowdown in supercooled water, a ubiquitous yet extraordinary substance characterized by various anomalous properties arising from local density fluctuations. Using molecular dynamics simulations, we found that the jump dynamics, which are elementary structural change processes, deviate from Poisson statistics with decreasing temperature. This deviation is attributed to slow variables competing with the jump motions, i.e., dynamic disorder. The present analysis of the dynamic disorder showed that the primary slow variable is the displacement of the fourth nearest oxygen atom of a jumping molecule, which occurs in an environment created by the fluctuations of molecules outside the first hydration shell. As the temperature decreases, the jump dynamics become slow and intermittent. These intermittent dynamics are attributed to the prolonged trapping of jumping molecules within extended and stable low-density domains. As the temperature continues to decrease, the number of slow variables increases due to the increased cooperative motions. Consequently, the jump dynamics proceed in a higher-dimensional space consisting of multiple slow variables, becoming slower and more intermittent. It is then conceivable that with further decreasing temperature, the slowing and intermittency of the jump dynamics intensify, eventually culminating in a glass transition.
Collapse
Affiliation(s)
- Shinji Saito
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan and The Graduate University for Advanced Studies (SOKENDAI), Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
4
|
Abraham GR, Chaderjian AS, N Nguyen AB, Wilken S, Saleh OA. Nucleic acid liquids. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2024; 87:066601. [PMID: 38697088 DOI: 10.1088/1361-6633/ad4662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Accepted: 05/02/2024] [Indexed: 05/04/2024]
Abstract
The confluence of recent discoveries of the roles of biomolecular liquids in living systems and modern abilities to precisely synthesize and modify nucleic acids (NAs) has led to a surge of interest in liquid phases of NAs. These phases can be formed primarily from NAs, as driven by base-pairing interactions, or from the electrostatic combination (coacervation) of negatively charged NAs and positively charged molecules. Generally, the use of sequence-engineered NAs provides the means to tune microsopic particle properties, and thus imbue specific, customizable behaviors into the resulting liquids. In this way, researchers have used NA liquids to tackle fundamental problems in the physics of finite valence soft materials, and to create liquids with novel structured and/or multi-functional properties. Here, we review this growing field, discussing the theoretical background of NA liquid phase separation, quantitative understanding of liquid material properties, and the broad and growing array of functional demonstrations in these materials. We close with a few comments discussing remaining open questions and challenges in the field.
Collapse
Affiliation(s)
- Gabrielle R Abraham
- Physics Department,University of California, Santa Barbara, CA 93106, United States of America
| | - Aria S Chaderjian
- Physics Department,University of California, Santa Barbara, CA 93106, United States of America
| | - Anna B N Nguyen
- Biomolecular Science and Engineering Program, University of California, Santa Barbara, CA 93106, United States of America
| | - Sam Wilken
- Physics Department,University of California, Santa Barbara, CA 93106, United States of America
- Materials Department, University of California, Santa Barbara, CA 93106, United States of America
| | - Omar A Saleh
- Physics Department,University of California, Santa Barbara, CA 93106, United States of America
- Biomolecular Science and Engineering Program, University of California, Santa Barbara, CA 93106, United States of America
- Materials Department, University of California, Santa Barbara, CA 93106, United States of America
| |
Collapse
|
5
|
Dhabal D, Kumar R, Molinero V. Liquid-liquid transition and ice crystallization in a machine-learned coarse-grained water model. Proc Natl Acad Sci U S A 2024; 121:e2322853121. [PMID: 38709921 PMCID: PMC11098087 DOI: 10.1073/pnas.2322853121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2023] [Accepted: 03/27/2024] [Indexed: 05/08/2024] Open
Abstract
Mounting experimental evidence supports the existence of a liquid-liquid transition (LLT) in high-pressure supercooled water. However, fast crystallization of supercooled water has impeded identification of the LLT line TLL(p) in experiments. While the most accurate all-atom (AA) water models display a LLT, their computational cost limits investigations of its interplay with ice formation. Coarse-grained (CG) models provide over 100-fold computational efficiency gain over AA models, enabling the study of water crystallization, but have not yet shown to have a LLT. Here, we demonstrate that the CG machine-learned water model Machine-Learned Bond-Order Potential (ML-BOP) has a LLT that ends in a critical point at pc = 170 ± 10 MPa and Tc = 181 ± 3 K. The TLL(p) of ML-BOP is almost identical to the one of TIP4P/2005, adding to the similarity in the equation of state of liquid water in both models. Cooling simulations reveal that ice crystallization is fastest at the LLT and its supercritical continuation of maximum heat capacity, supporting a mechanistic relationship between the structural transformation of water to a low-density liquid (LDL) and ice formation. We find no signature of liquid-liquid criticality in the ice crystallization temperatures. ML-BOP replicates the competition between formation of LDL and ice observed in ultrafast experiments of decompression of the high-density liquid (HDL) into the region of stability of LDL. The simulations reveal that crystallization occurs prior to the coarsening of the HDL and LDL domains, obscuring the distinction between the highly metastable first-order LLT and pronounced structural fluctuations along its supercritical continuation.
Collapse
Affiliation(s)
- Debdas Dhabal
- Department of Chemistry, The University of Utah, Salt Lake City, UT84112-0850
| | - Rajat Kumar
- Department of Chemistry, The University of Utah, Salt Lake City, UT84112-0850
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, UT84112-0850
| |
Collapse
|
6
|
Loubet NA, Verde AR, Appignanesi GA. A water structure indicator suitable for generic contexts: Two-liquid behavior at hydration and nanoconfinement conditions and a molecular approach to hydrophobicity and wetting. J Chem Phys 2024; 160:144502. [PMID: 38587223 DOI: 10.1063/5.0203989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2024] [Accepted: 03/25/2024] [Indexed: 04/09/2024] Open
Abstract
In a recent work, we have briefly introduced a new structural index for water that, unlike previous indicators, was devised specifically for generic contexts beyond bulk conditions, making it suitable for hydration and nanoconfinement settings. In this work, we shall study this metric in detail, demonstrating its ability to reveal the existence of a fine-tuned interplay between the local structure and energetics in liquid water. This molecular principle enables the establishment of an extended hydrogen bond network, while simultaneously allowing for the existence of network defects by compensating for uncoordinated sites. By studying different water models and different temperatures encompassing both the normal liquid and the supercooled regime, this molecular mechanism will be shown to underlie the two-state behavior of bulk water. In addition, by studying functionalized self-assembled monolayers and diverse graphene-like surfaces, we shall show that this principle is also operative at hydration and nanoconfinement conditions, thus generalizing the validity of the two-liquid scenario of water to these contexts. This approach will allow us to define conditions for wettability, providing an accurate measure of hydrophobicity and a reliable predictor of filling and drying transitions. Hence, it might open the possibility of elucidating the active role of water in the broad fields of biophysics and materials science. As a preliminary step, we shall study the hydration structure and hydrophilicity of graphene-like systems (parallel graphene sheets and carbon nanotubes) as a function of the confinement dimensionality.
Collapse
Affiliation(s)
- Nicolás A Loubet
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Alejandro R Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
7
|
Bachler J, Daidone I, Zanetti-Polzi L, Loerting T. Tuning the low-temperature phase behavior of aqueous ionic liquids. Phys Chem Chem Phys 2024; 26:9741-9753. [PMID: 38470827 DOI: 10.1039/d3cp06101a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2024]
Abstract
Water's anomalous behavior is often explained using a two-liquid model, where two types of water, high-density liquid (HDL) and low-density liquid (LDL), can be separated via a liquid-liquid phase transition (LLPT) at low temperature. Mixtures of water and the ionic liquid hydrazinium trifluoroacetate were suggested to also show an LLPT but with the advantage that there is no rapid ice crystallization hampering its observation. It remains controversial whether these solutions exhibit an LLPT or are instead associated with complex phase separation phenomena. We here show detailed low-temperature calorimetry and diffraction experiments on aqueous solutions containing hydrazinium trifluoroacetate and other similar ionic liquids, all at a solute mole fraction of x = 0.175. Hydrazinium trifluoroacetate, ammonium trifluoroacetate, ethylammonium trifluoroacetate and hydrazinium pentafluoropropionate all boast exothermic transitions unrelated to crystallization as well as remarkable structural changes upon cooling into the glassy state. We propose a model inspired by micelle formation and decomposition in surfactant solutions, which is complemented by MD simulations and allows rationalizing the rich phase behavior of our mixtures during cooling. The fundamental aspect of the model is the hydrophobic nature of fluorinated anions that enables aggregation, which is reversed upon cooling and culminates in the remarkable exothermic first-order transition observed at low temperature. That is, we assign the first-order transition not to an LLPT but to phase-separations similar to the ones when falling below the Krafft temperature. All other solutions merely show simple vitrification behavior. Still, they exhibit distinct differences in liquid fragility, which is decreased continuously with decreasing hydrophobicity of the anions. This might enable the systematic tuning of ionic liquids with the goal of designing aqueous solutions of specific fragility.
Collapse
Affiliation(s)
- Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, Innsbruck A-6020, Austria.
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, L'Aquila 67010, Italy
| | | | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, Innsbruck A-6020, Austria.
| |
Collapse
|
8
|
Sciortino F, Gartner TE, Debenedetti PG. Free-energy landscape and spinodals for the liquid-liquid transition of the TIP4P/2005 and TIP4P/Ice models of water. J Chem Phys 2024; 160:104501. [PMID: 38456528 DOI: 10.1063/5.0196964] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 02/20/2024] [Indexed: 03/09/2024] Open
Abstract
Continued increases in computational power now make it possible to evaluate the free-energy landscape associated with the first-order liquid-liquid transition in realistic models of water for which an accurate estimate of the liquid-liquid critical point exists, and to explore its change with pressure near the coexistence line. We report the results of 50 μs-long NPT umbrella sampling simulations for two realistic models for water, TIP4P/2005 and TIP4P/ice, 3-9 K below their critical temperatures. The free energy profile at different pressures clearly shows the presence of two well-defined free energy basins and makes it possible to identify the liquid-liquid spinodal points, the limits of stability that define the (temperature dependent) pressure range within which two distinct free energy basins exist. The results show that for temperatures less than 10 K below the critical temperature, metastable states are possible across a very limited pressure interval, information that is relevant to the interpretation of experiments probing the metastable phase behavior of deeply supercooled water in the so-called no-man's land.
Collapse
Affiliation(s)
- Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Rome, Italy
| | - Thomas E Gartner
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, USA
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
9
|
Maturi F, Raposo Filho RS, Brites CDS, Fan J, He R, Zhuang B, Liu X, Carlos LD. Deciphering Density Fluctuations in the Hydration Water of Brownian Nanoparticles via Upconversion Thermometry. J Phys Chem Lett 2024; 15:2606-2615. [PMID: 38420927 PMCID: PMC10926164 DOI: 10.1021/acs.jpclett.4c00044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 02/14/2024] [Accepted: 02/16/2024] [Indexed: 03/02/2024]
Abstract
We investigate the intricate relationship among temperature, pH, and Brownian velocity in a range of differently sized upconversion nanoparticles (UCNPs) dispersed in water. These UCNPs, acting as nanorulers, offer insights into assessing the relative proportion of high-density and low-density liquid in the surrounding hydration water. The study reveals a size-dependent reduction in the onset temperature of liquid-water fluctuations, indicating an augmented presence of high-density liquid domains at the nanoparticle surfaces. The observed upper-temperature threshold is consistent with a hypothetical phase diagram of water, validating the two-state model. Moreover, an increase in pH disrupts the organization of water molecules, similar to external pressure effects, allowing simulation of the effects of temperature and pressure on hydrogen bonding networks. The findings underscore the significance of the surface of suspended nanoparticles for understanding high- to low-density liquid fluctuations and water behavior at charged interfaces.
Collapse
Affiliation(s)
- Fernando
E. Maturi
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
- Institute
of Chemistry, São Paulo State University
(UNESP), 14800-060 Araraquara, SP, Brazil
| | - Ramon S. Raposo Filho
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Carlos D. S. Brites
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Jingyue Fan
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Ruihua He
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Bilin Zhuang
- Harvey
Mudd College, 301 Platt
Boulevard, Claremont, California 91711, United States
| | - Xiaogang Liu
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Luís D. Carlos
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| |
Collapse
|
10
|
Eltareb A, Lopez GE, Giovambattista N. A continuum of amorphous ices between low-density and high-density amorphous ice. Commun Chem 2024; 7:36. [PMID: 38378859 PMCID: PMC10879119 DOI: 10.1038/s42004-024-01117-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 02/01/2024] [Indexed: 02/22/2024] Open
Abstract
Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY, 11210, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| | - Gustavo E Lopez
- Department of Chemistry, Lehman College of the City University of New York, Bronx, NY, 10468, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| | - Nicolas Giovambattista
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY, 11210, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| |
Collapse
|
11
|
Finney JL. The structure of water: A historical perspective. J Chem Phys 2024; 160:060901. [PMID: 38341786 DOI: 10.1063/5.0182665] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/11/2024] [Indexed: 02/13/2024] Open
Abstract
Attempts to understand the molecular structure of water were first made well over a century ago. Looking back at the various attempts, it is illuminating to see how these were conditioned by the state of knowledge of chemistry and physics at the time and the experimental and theoretical tools then available. Progress in the intervening years has been facilitated by not only conceptual and theoretical advances in physics and chemistry but also the development of experimental techniques and instrumentation. Exploitation of powerful computational methods in interpreting what at first sight may seem impenetrable experimental data has led us to the consistent and detailed picture we have today of not only the structure of liquid water itself and how it changes with temperature and pressure but also its interactions with other molecules, in particular those relevant to water's role in important chemical and biological processes. Much remains to be done in the latter areas, but the experimental and computational techniques that now enable us to do what might reasonably be termed "liquid state crystallography" have opened the door to make possible further advances. Consequently, we now have the tools to explore further the role of water in those processes that underpin life itself-the very prospect that inspired Bernal to develop his ideas on the structure of liquids in general and of water in particular.
Collapse
Affiliation(s)
- John L Finney
- Department of Physics and Astronomy and London Centre for Nanotechnology, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
12
|
Fan Z, Tanaka H. Microscopic mechanisms of pressure-induced amorphous-amorphous transitions and crystallisation in silicon. Nat Commun 2024; 15:368. [PMID: 38228606 DOI: 10.1038/s41467-023-44332-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Accepted: 12/08/2023] [Indexed: 01/18/2024] Open
Abstract
Some low-coordination materials, including water, silica, and silicon, exhibit polyamorphism, having multiple amorphous forms. However, the microscopic mechanism and kinetic pathway of amorphous-amorphous transition (AAT) remain largely unknown. Here, we use a state-of-the-art machine-learning potential and local structural analysis to investigate the microscopic kinetics of AAT in silicon after a rapid pressure change. We find that the transition from low-density-amorphous (LDA) to high-density-amorphous (HDA) occurs through nucleation and growth, resulting in non-spherical interfaces that underscore the mechanical nature of AAT. In contrast, the reverse transition occurs through spinodal decomposition. Further pressurisation transforms LDA into very-high-density amorphous (VHDA), with HDA serving as an intermediate state. Notably, the final amorphous states are inherently unstable, transitioning into crystals. Our findings demonstrate that AAT and crystallisation are driven by joint thermodynamic and mechanical instabilities, assisted by preordering, occurring without diffusion. This unique mechanical and diffusion-less nature distinguishes AAT from liquid-liquid transitions.
Collapse
Affiliation(s)
- Zhao Fan
- Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8904, Japan
| | - Hajime Tanaka
- Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8904, Japan.
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan.
| |
Collapse
|
13
|
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
|
14
|
Palmer JC, Sarupria S, Truskett TM. Tribute to Pablo G. Debenedetti. J Phys Chem B 2023; 127:8075-8078. [PMID: 37766640 DOI: 10.1021/acs.jpcb.3c06020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/29/2023]
Affiliation(s)
- Jeremy C Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, United States
| | - Sapna Sarupria
- Department of Chemistry, Chemical Theory Center, University of Minnesota Twin Cities, Minneapolis, Minnesota 55455, United States
| | - Thomas M Truskett
- McKetta Department of Chemical Engineering and Department of Physics, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
15
|
Daidone I, Foffi R, Amadei A, Zanetti-Polzi L. A statistical mechanical model of supercooled water based on minimal clusters of correlated molecules. J Chem Phys 2023; 159:094502. [PMID: 37655770 DOI: 10.1063/5.0157505] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 08/14/2023] [Indexed: 09/02/2023] Open
Abstract
In this paper, we apply a theoretical model for fluid state thermodynamics to investigate simulated water in supercooled conditions. This model, which we recently proposed and applied to sub- and super-critical fluid water [Zanetti-Polzi et al., J. Chem. Phys. 156(4), 44506 (2022)], is based on a combination of the moment-generating functions of the enthalpy and volume fluctuations as provided by two gamma distributions and provides the free energy of the system as well as other relevant thermodynamic quantities. The application we make here provides a thermodynamic description of supercooled water fully consistent with that expected by crossing the liquid-liquid Widom line, indicating the presence of two distinct liquid states. In particular, the present model accurately reproduces the Widom line temperatures estimated with other two-state models and well describes the heat capacity anomalies. Differently from previous models, according to our description, a cluster of molecules that extends beyond the first hydration shell is necessary to discriminate between the statistical fluctuation regimes typical of the two liquid states.
Collapse
Affiliation(s)
- Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| | - Riccardo Foffi
- Department of Civil, Environmental and Geomatic Engineering, Institute of Environmental Engineering, ETH Zürich, Laura-Hezner-Weg 7, 8093 Zürich, Switzerland
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome "Tor Vergata," Via della Ricerca Scientifica, I-00185 Rome, Italy
| | - Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| |
Collapse
|
16
|
Loubet NA, Verde AR, Lockhart JA, Appignanesi GA. Turning an energy-based defect detector into a multi-molecule structural indicator for water. J Chem Phys 2023; 159:064512. [PMID: 37578063 DOI: 10.1063/5.0159060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2023] [Accepted: 07/27/2023] [Indexed: 08/15/2023] Open
Abstract
Recent studies have provided conclusive evidence for the existence of a liquid-liquid critical point in numerical models of water. Such a scenario implies the competition between two local molecular arrangements of different densities: a high-density liquid (HDL) and a low-density liquid (LDL). Within this context, the development of accurate structural indicators to properly characterize the two interconverting local structures is demanded. In a previous study, we introduced a reliable energy-based structural descriptor that properly discriminates water molecules into tetrahedrally arranged molecules (T molecules) and distorted molecules (D molecules). The latter constitute defects in terms of hydrogen bond (HB) coordination and have been shown to represent a minority component, even at high temperatures above the melting point. In addition, the D molecules tend to form high-quality HBs with three T molecules and to be surrounded by T and D molecules at further distances. Thus, it became evident that, while the LDL state might consist of a virtually pure T state, the HDL state would comprise mixed molecular arrangements including the D molecules. Such a need to abandon the single-molecule description requires the investigation of the degree of structural information to be incorporated in order to build an appropriate multi-molecule indicator. Hence, in this work, we shall study the effect of the local structural constraints on the water molecules in order to discriminate the different molecular arrangements into two disjoint classes. This will enable us to build a multi-molecule structural indicator for water whose performance will then be investigated within the water's supercooled regime.
Collapse
Affiliation(s)
- Nicolás A Loubet
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Alejandro R Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Jano A Lockhart
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
17
|
Piaggi PM, Gartner TE, Car R, Debenedetti PG. Melting curves of ice polymorphs in the vicinity of the liquid-liquid critical point. J Chem Phys 2023; 159:054502. [PMID: 37531247 DOI: 10.1063/5.0159288] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Accepted: 07/14/2023] [Indexed: 08/04/2023] Open
Abstract
The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164-168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. We explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.
Collapse
Affiliation(s)
- Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Thomas E Gartner
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
- Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
18
|
Li Z, Turner J, Snurr RQ. Computational investigation of hysteresis and phase equilibria of n-alkanes in a metal-organic framework with both micropores and mesopores. Commun Chem 2023; 6:90. [PMID: 37156883 PMCID: PMC10167368 DOI: 10.1038/s42004-023-00889-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Accepted: 04/19/2023] [Indexed: 05/10/2023] Open
Abstract
Adsorption hysteresis is a phenomenon related to phase transitions that can impact applications such as gas storage and separations in porous materials. Computational approaches can greatly facilitate the understanding of phase transitions and phase equilibria in porous materials. In this work, adsorption isotherms for methane, ethane, propane, and n-hexane were calculated from atomistic grand canonical Monte Carlo (GCMC) simulations in a metal-organic framework having both micropores and mesopores to better understand hysteresis and phase equilibria between connected pores of different size and the external bulk fluid. At low temperatures, the calculated isotherms exhibit sharp steps accompanied by hysteresis. As a complementary simulation method, canonical (NVT) ensemble simulations with Widom test particle insertions are demonstrated to provide additional information about these systems. The NVT+Widom simulations provide the full van der Waals loop associated with the sharp steps and hysteresis, including the locations of the spinodal points and points within the metastable and unstable regions that are inaccessible to GCMC simulations. The simulations provide molecular-level insight into pore filling and equilibria between high- and low-density states within individual pores. The effect of framework flexibility on adsorption hysteresis is also investigated for methane in IRMOF-1.
Collapse
Affiliation(s)
- Zhao Li
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL, 60208, USA
| | - Jake Turner
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL, 60208, USA
| | - Randall Q Snurr
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL, 60208, USA.
| |
Collapse
|
19
|
Guidarelli Mattioli F, Sciortino F, Russo J. Are Neural Network Potentials Trained on Liquid States Transferable to Crystal Nucleation? A Test on Ice Nucleation in the mW Water Model. J Phys Chem B 2023; 127:3894-3901. [PMID: 37075256 PMCID: PMC10165654 DOI: 10.1021/acs.jpcb.3c00693] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 04/06/2023] [Indexed: 04/21/2023]
Abstract
Neural network potentials (NNPs) are increasingly being used to study processes that happen on long time scales. A typical example is crystal nucleation, which rate is controlled by the occurrence of a rare fluctuation, i.e., the appearance of the critical nucleus. Because the properties of this nucleus are far from those of the bulk crystal, it is yet unclear whether NN potentials trained on equilibrium liquid states can accurately describe nucleation processes. So far, nucleation studies on NNPs have been limited to ab initio models whose nucleation properties are unknown, preventing an accurate comparison. Here we train a NN potential on the mW model of water─a classical three-body potential whose nucleation time scale is accessible in standard simulations. We show that a NNP trained only on a small number of liquid state points can reproduce with great accuracy the nucleation rates and free energy barriers of the original model, computed from both spontaneous and biased trajectories, strongly supporting the use of NNPs to study nucleation events.
Collapse
Affiliation(s)
| | | | - John Russo
- Sapienza University of Rome, Piazzale Aldo Moro 2, 00185 Rome, Italy
| |
Collapse
|
20
|
Piskulich ZA, Borkowski AK, Thompson WH. A Maxwell relation for dynamical timescales with application to the pressure and temperature dependence of water self-diffusion and shear viscosity. Phys Chem Chem Phys 2023; 25:12820-12832. [PMID: 37129891 DOI: 10.1039/d3cp01386c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
A Maxwell relation for a reaction rate constant (or other dynamical timescale) obtained under constant pressure, p, and temperature, T, is introduced and discussed. Examination of this relationship in the context of fluctuation theory provides insight into the p and T dependence of the timescale and the underlying molecular origins. This Maxwell relation motivates a suggestion for the general form of the timescale as a function of pressure and temperature. This is illustrated by accurately fitting simulation results and existing experimental data on the self-diffusion coefficient and shear viscosity of liquid water. A key advantage of this approach is that each fitting parameter is physically meaningful.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, KS 66045, USA.
| | | | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, KS 66045, USA.
| |
Collapse
|
21
|
Kringle L, Kay BD, Kimmel GA. Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water. J Phys Chem B 2023; 127:3919-3930. [PMID: 37097190 DOI: 10.1021/acs.jpcb.3c01465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.
Collapse
Affiliation(s)
- Loni Kringle
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Bruce D Kay
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Greg A Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|
22
|
Cerdeiriña CA, González-Salgado D, Troncoso J. Liquid-Liquid Criticality in TIP4P/2005 and Three-State Models of Water. J Phys Chem B 2023; 127:3902-3910. [PMID: 37097210 PMCID: PMC10165646 DOI: 10.1021/acs.jpcb.3c00696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Molecular dynamics simulations leading to the isothermal compressibility, the isobaric thermal expansivity, and the isobaric heat capacity of TIP4P/2005 water are found to be consistent with the coordinates of its second, liquid-liquid critical point reported recently by Debenedetti et al. [ Science 2020, 369, 289-292]. In accord with the theory of critical phenomena, we encounter that the rise in the magnitude of these response functions as temperature is lowered is especially marked along the critical isochore. Furthermore, response-function ratios provide a test for thermodynamic consistency at the critical point and manifest nonuniversal features sharply distinguishing liquid-liquid from standard gas-liquid criticality. The whole pattern of behavior revealed by simulations is qualitatively the same as the one of a three-state Ising model of water exhibiting a low-temperature liquid-liquid critical point. Exact solutions for the two-state components of such a three-state model are also provided.
Collapse
Affiliation(s)
- Claudio A Cerdeiriña
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQFR, Ourense 32004, Spain
| | - Diego González-Salgado
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQFR, Ourense 32004, Spain
| | - Jacobo Troncoso
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQFR, Ourense 32004, Spain
| |
Collapse
|
23
|
Singh Y, Santra M, Singh RS. Anomalous Vapor and Ice Nucleation in Water at Negative Pressures: A Classical Density Functional Theory Study. J Phys Chem B 2023; 127:3312-3324. [PMID: 36989467 DOI: 10.1021/acs.jpcb.2c09136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/31/2023]
Abstract
In contrast to the abundance of work on the anomalous behavior of water, the relationship between the water's thermodynamic anomalies and kinetics of phase transition from metastable water is relatively unexplored. In this work, we have employed classical density functional theory to provide a unified and coherent picture of nucleation (both vapor and ice) from metastable water at negative pressure conditions. Our results suggest a peculiar nonmonotonic temperature dependence of vapor-liquid surface tension at temperatures where vapor-liquid coexistence is metastable with respect to the ice phase. The vapor nucleation barrier on isochoric cooling also shows a nonmonotonic temperature dependence. We further report that, for low density isochores, the temperature of the minimum vapor nucleation barrier (TΔΩv/min*) does not coincide with the temperature of maximum density (TMD) where metastability is maximum. The difference between the TΔΩv/min* and the TMD, however, decreases with increasing the density of the isochore. The vapor nucleation barrier along isobars shows an interesting crossover behavior in the vicinity of the Widom line on lowering the temperature. Our results on the ice nucleation suggest an anomalous retracing behavior of the nucleation barrier along isotherms at negative pressures and theoretically validate the recent findings that the reentrant ice(Ih)-liquid coexistence line can induce a drastic change in the kinetics of ice nucleation. Thus, this study establishes a direct connection between the metastable water's thermodynamic anomalies and the (vapor and ice) nucleation kinetics. In addition, this study provides deeper insights into the origin of the isothermal compressibility maximum on isochoric cooling.
Collapse
Affiliation(s)
- Yuvraj Singh
- Department of Physics, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| | - Mantu Santra
- School of Chemical and Materials Sciences, Indian Institute of Technology Goa, Ponda, Goa 403401, India
| | - Rakesh S Singh
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| |
Collapse
|
24
|
Gorfer A, Dellago C, Sega M. High-density liquid (HDL) adsorption at the supercooled water/vapor interface and its possible relation to the second surface tension inflection point. J Chem Phys 2023; 158:054503. [PMID: 36754827 DOI: 10.1063/5.0132985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
We investigate the properties of water along the liquid/vapor coexistence line in the supercooled regime down to the no-man's land. Extensive molecular dynamics simulations of the TIP4P/2005 liquid/vapor interface in the range 198-348 K allow us to locate the second surface tension inflection point with a high accuracy at 283 ± 5 K, close to the temperature of maximum density. This temperature also coincides with the appearance of a density anomaly at the interface known as the apophysis. We relate the emergence of the apophysis to the observation of high-density liquid (HDL) water adsorption in the proximity of the liquid/vapor interface.
Collapse
Affiliation(s)
- Alexander Gorfer
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, Wien A-1090, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, Wien A-1090, Austria
| | - Marcello Sega
- Department of Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
| |
Collapse
|
25
|
Takahashi O, Pettersson LGM. Dynamical and interference effects in X-ray emission spectroscopy of H-bonded water – origin of the split lone-pair peaks. Mol Phys 2023. [DOI: 10.1080/00268976.2023.2170686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Osamu Takahashi
- Basic Chemistry Program, Graduate School of Advanced Science and Engineering, Hiroshima University, Higashi-Hiroshima, Japan
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| |
Collapse
|
26
|
Amann-Winkel K, Kim KH, Giovambattista N, Ladd-Parada M, Späh A, Perakis F, Pathak H, Yang C, Eklund T, Lane TJ, You S, Jeong S, Lee JH, Eom I, Kim M, Park J, Chun SH, Poole PH, Nilsson A. Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice. Nat Commun 2023; 14:442. [PMID: 36707522 PMCID: PMC9883474 DOI: 10.1038/s41467-023-36091-1] [Citation(s) in RCA: 16] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Accepted: 01/13/2023] [Indexed: 01/28/2023] Open
Abstract
Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale < 100 ns and widely separated from ice formation, which begins at times >1 μs.
Collapse
Affiliation(s)
- Katrin Amann-Winkel
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden ,grid.419547.a0000 0001 1010 1663Max Planck Institute for Polymer Research and Johannes Gutenberg University, 55128 Mainz, Germany
| | - Kyung Hwan Kim
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Nicolas Giovambattista
- grid.183006.c0000 0001 0671 7844Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY 11210 USA ,grid.253482.a0000 0001 0170 7903The Graduate Center of the City University of New York, New York, NY 10016 USA
| | - Marjorie Ladd-Parada
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden ,grid.411313.50000 0004 0512 3288Present Address: Division of Glycoscience, School of Biotechnology, KTH Royal Institute of Technology, AlbaNova University Center, SE-10691 Stockholm, Sweden
| | - Alexander Späh
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Fivos Perakis
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Harshad Pathak
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Cheolhee Yang
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Tobias Eklund
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Thomas J. Lane
- grid.445003.60000 0001 0725 7771SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025 USA
| | - Seonju You
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Sangmin Jeong
- grid.49100.3c0000 0001 0742 4007Department of Chemistry, POSTECH, Pohang, 37673 Republic of Korea
| | - Jae Hyuk Lee
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Intae Eom
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Minseok Kim
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Jaeku Park
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Sae Hwan Chun
- grid.49100.3c0000 0001 0742 4007Pohang Accelerator Laboratory, Pohang, Gyeongbuk 37673 Republic of Korea
| | - Peter H. Poole
- grid.264060.60000 0004 1936 7363Department of Physics, St. Francis Xavier University, Antigonish, NS B2G2W5 Canada
| | - Anders Nilsson
- grid.10548.380000 0004 1936 9377Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
27
|
Jedrecy A, Saitta AM, Pietrucci F. Free energy calculations and unbiased molecular dynamics targeting the liquid-liquid transition in water no man's land. J Chem Phys 2023; 158:014502. [PMID: 36610960 DOI: 10.1063/5.0120789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The existence of a first-order phase transition between a low-density liquid (LDL) and a high-density liquid (HDL) form of supercooled water has been a central and highly debated issue of physics and chemistry for the last three decades. We present a computational study that allows us to determine the free-energy landscapes of supercooled water over a wide range of pressure and temperature conditions using the TIP4P/2005 force field. Our approach combines topology-based structural transformation coordinates, state-of-the-art free-energy calculation methods, and extensive unbiased molecular dynamics. All our diverse simulations cannot detect any barrier within the investigated timescales and system size, for a discontinuous transition between the LDL and HDL forms throughout the so-called "no man's land," until the onset of the solid, non-diffusive amorphous forms.
Collapse
Affiliation(s)
- Alexandre Jedrecy
- Insitut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, CNRS, MNHN, UMR 7590, Paris, France
| | - A Marco Saitta
- Insitut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, CNRS, MNHN, UMR 7590, Paris, France
| | - Fabio Pietrucci
- Insitut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, CNRS, MNHN, UMR 7590, Paris, France
| |
Collapse
|
28
|
Foffi R, Sciortino F. Correlated Fluctuations of Structural Indicators Close to the Liquid-Liquid Transition in Supercooled Water. J Phys Chem B 2022; 127:378-386. [PMID: 36538764 PMCID: PMC9841516 DOI: 10.1021/acs.jpcb.2c07169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Multiple numerical studies have unambiguously shown the existence of a liquid-liquid critical point in supercooled states for different numerical models of water, and various structural indicators have been put forward to describe the transformation associated with this phase transition. Here we analyze numerical simulations of near-critical supercooled water to compare the behavior of several of such indicators with critical density fluctuations. We show that close to the critical point most indicators are strongly correlated to density, and some of them even display identical distributions of fluctuations. These indicators probe the exact same free energy landscape, therefore providing a thermodynamic description of critical supercooled water which is identical to that provided by the density order parameter. This implies that close to the critical point, there is a tight coupling between many, only apparently distinct, structural degrees of freedom.
Collapse
Affiliation(s)
- Riccardo Foffi
- Institute
for Environmental Engineering, Department of Civil, Environmental
and Geomatic Engineering, ETH Zürich, 8093Zürich, Switzerland
| | - Francesco Sciortino
- Dipartimento
di Fisica, Sapienza Università di
Roma, Piazzale Aldo Moro
5, I-00185Rome, Italy,E-mail:
| |
Collapse
|
29
|
Gartner TE, Piaggi PM, Car R, Panagiotopoulos AZ, Debenedetti PG. Liquid-Liquid Transition in Water from First Principles. PHYSICAL REVIEW LETTERS 2022; 129:255702. [PMID: 36608224 DOI: 10.1103/physrevlett.129.255702] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Accepted: 11/14/2022] [Indexed: 06/17/2023]
Abstract
A long-standing question in water research is the possibility that supercooled liquid water can undergo a liquid-liquid phase transition (LLT) into high- and low-density liquids. We used several complementary molecular simulation techniques to evaluate the possibility of an LLT in an ab initio neural network model of water trained on density functional theory calculations with the SCAN exchange correlation functional. We conclusively show the existence of a first-order LLT and an associated critical point in the SCAN description of water, representing the first definitive computational evidence for an LLT in water from first principles.
Collapse
Affiliation(s)
- Thomas E Gartner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
- Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
- Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
30
|
Camisasca G, Tenuzzo L, Gallo P. Protein hydration water: focus on low density and high density local structures upon cooling. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
31
|
Mondal A, Ramesh G, Singh RS. Manifestations of the structural origin of supercooled water’s anomalies in the heterogeneous relaxation on the potential energy landscape. J Chem Phys 2022; 157:184503. [DOI: 10.1063/5.0124041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
Liquid water is well-known for its intriguing thermodynamic anomalies in the supercooled state. The phenomenological two-state models—based on the assumption of the existence of two types of competing local states (or, structures) in liquid water—have been extremely successful in describing water’s thermodynamic anomalies. However, the precise structural features of these competing local states in liquid water still remain elusive. Here, we have employed a predefined structural order parameter-free approach to unambiguously identify two types of competing local states—entropically and energetically favored—with significantly different structural and energetic features in the TIP4P/2005 liquid water. This identification is based on the heterogeneous structural relaxation of the system in the potential energy landscape (PEL) during the steepest-descent energy minimization. This heterogeneous relaxation is characterized using order parameters inspired by the spin-glass transition in frustrated magnetic systems. We have further established a direct relationship between the population fluctuation of the two states and the anomalous behavior of the heat capacity in supercooled water. The composition-dependent spatial distribution of the entropically favored local states shows an interesting crossover from a spanning network-like single cluster to the spatially delocalized clusters in the close vicinity of the Widom line. Additionally, this study establishes a direct relationship between the topographic features of the PEL and the water’s thermodynamic anomalies in the supercooled state and provides alternate markers (in addition to the locus of maxima of thermodynamic response functions) for the Widom line in the phase plane.
Collapse
Affiliation(s)
- Arijit Mondal
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| | - Gadha Ramesh
- Department of Physics, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| | - Rakesh S. Singh
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| |
Collapse
|
32
|
Cerdeiriña CA. Water's Unusual Thermodynamics in the Realm of Physical Chemistry. J Phys Chem B 2022; 126:6608-6613. [PMID: 36001372 PMCID: PMC9797112 DOI: 10.1021/acs.jpcb.2c05274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
While it is known since the early work by Edsall, Frank and Evans, Kauzmann, and others that the thermodynamics of solvation of nonpolar solutes in water is unusual and has implications for the thermodynamics of protein folding, only recently have its connections with the unusual temperature dependence of the density of solvent water been illuminated. Such density behavior is, in turn, one of the manifestations of a nonstandard thermodynamic pattern contemplating a second, liquid-liquid critical point at conditions of temperature and pressure at which water exists as a deeply supercooled liquid. Recent experimental and computational work unambiguously points toward the existence of such a critical point, thereby providing concrete answers to the questions posed by the 1976 pioneering experiments by Speedy and Angell and the associated "liquid-liquid transition hypothesis" posited in 1992 by Stanley and co-workers. Challenges of this phenomenology to the branch of Statistical Mechanics remain.
Collapse
|
33
|
Bachler J, Giebelmann J, Amann-Winkel K, Loerting T. Pressure-annealed high-density amorphous ice made from vitrified water droplets: A systematic calorimetry study on water's second glass transition. J Chem Phys 2022; 157:064502. [DOI: 10.1063/5.0100571] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In previous work, water's second glass transition was investigated based on an amorphous sample made from crystalline ice (Amann-Winkel et al., Proc. Natl. Acad. Sci. U.S.A. 110 (44) 17720-17725). In the present work, we investigate water's second glass transition based on the genuine glassy state of high-density water as prepared from micron-sized liquid water droplets, avoiding crystallinity at all stages. All the calorimetric features of water's second glass transition observed in the previous work are also observed here on the genuine glassy samples. This suggests that the glass transition indeed thermodynamically links amorphous ices continuously with deeply supercooled water. We proceed to extend the earlier study by investigating the effect of preparation history on the calorimetric glass transition temperature. The best samples prepared here feature both a lower glass transition temperature Tg,2 and a higher polyamorphic transition temperature Tons, thereby extending the range of thermal stability in which the deeply supercooled liquid can be observed by about 4 K. Just before the polyamorphic transition, we observe a spike-like increase of heat capacity that we interpret in terms of nucleation of low-density water. Without this spike, the width of water's second glass transition is 15 K, and the Δcp amounts to 3{plus minus}1 J K-1 mol-1, making the case for HDL being a strong liquid. We suggest that samples annealed at 1.9 GPa to 175 K and decompressed at 140 K to {greater than or equal to}0.10 GPa are free from such nuclei and represent the most ideal HDA glasses.
Collapse
Affiliation(s)
- Johannes Bachler
- University of Innsbruck Institute of Physical Chemistry, Austria
| | | | | | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck, Austria
| |
Collapse
|
34
|
Yun Y, Khaliullin RZ, Jung Y. Correlated Local Fluctuations in the Hydrogen Bond Network of Liquid Water. J Am Chem Soc 2022; 144:13127-13136. [PMID: 35820142 DOI: 10.1021/jacs.2c02362] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The hypothesis that liquid water can separate into two phases in the supercooled state has been supported by recent experimental and theoretical studies. However, whether such structural inhomogeneity extends to ambient conditions is under intense debate. Due to the dynamic nature of the hydrogen bond network of liquid water, exploring its structure requires detailed insight into the collective motion of neighboring water molecules, a missing link that has not been examined so far. Here, highly sensitive quantum mechanical calculations detect that the time evolution of nearby hydrogen bonds is strongly correlated, revealing a direct mechanism for the appearance of short-range structural fluctuations in the hydrogen bond network of liquid water for the first time. This correlated dynamics is found to be closely connected to the static structural picture. The distortions from the tetrahedral structure do not occur independently but are correlated due to the preference of nearby donors and acceptors to be in similar environments. The existence of such cooperative fluctuations is further supported by the temperature dependence of the local structural evolution and explained by conventional analysis of localized orbitals. It was found that such correlated structural fluctuations are only observed on a short length scale in simulations at ambient conditions. The correlations of the nearby hydrogen bond pairs of liquid water unveiled here are expected to offer a new insight into connecting the dynamics of individual water molecules and the local structure of the hydrogen bond network.
Collapse
Affiliation(s)
- Yonghwan Yun
- Department of Chemical and Biomolecular Engineering (BK21 four), Korea Advanced Institute of Science and Technology (KAIST), 291 Daehak-ro, Daejeon 34141, South Korea
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, 801 Sherbrooke Street West, Montreal, Quebec H3A 0B8, Canada
| | - Yousung Jung
- Department of Chemical and Biomolecular Engineering (BK21 four), Korea Advanced Institute of Science and Technology (KAIST), 291 Daehak-ro, Daejeon 34141, South Korea
| |
Collapse
|
35
|
Structural changes across thermodynamic maxima in supercooled liquid tellurium: A water-like scenario. Proc Natl Acad Sci U S A 2022; 119:e2202044119. [PMID: 35867742 DOI: 10.1073/pnas.2202044119] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Liquid polymorphism is an intriguing phenomenon that has been found in a few single-component systems, the most famous being water. By supercooling liquid Te to more than 130 K below its melting point and performing simultaneous small-angle and wide-angle X-ray scattering measurements, we observe clear maxima in its thermodynamic response functions around 615 K, suggesting the possible existence of liquid polymorphism. A close look at the underlying structural evolution shows the development of intermediate-range order upon cooling, most strongly around the thermodynamic maxima, which we attribute to bond-orientational ordering. The striking similarities between our results and those of water, despite the lack of hydrogen-bonding and tetrahedrality in Te, indicate that water-like anomalies may be a general phenomenon among liquid systems with competing bond- and density-ordering.
Collapse
|
36
|
Coexistence of vitreous and crystalline phases of H 2O at ambient temperature. Proc Natl Acad Sci U S A 2022; 119:e2117281119. [PMID: 35763575 PMCID: PMC9271169 DOI: 10.1073/pnas.2117281119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Understanding the phase behavior of H2O is essential in geoscience, extreme biology, biological imaging, chemistry, and physics. Vitreous phases of H2O are of particular importance since these phases avoid the typical expansion of H2O during ice-crystal formation. Here, we confirm the existence of vitreous ice and ice VI in mixed-phase samples of H2O at room temperature and high pressure. We show how Raman scattering and X-ray diffraction alone lead to misleading characterization and understanding of the mixed-phase material, a conclusion supported by molecular dynamics simulations. The coexistence of vitreous and crystalline components of H2O under these conditions is crucial for experimental studies of biological systems. The results have implications for related metastable transitions in other materials under pressure. Formation of vitreous ice during rapid compression of water at room temperature is important for biology and the study of biological systems. Here, we show that Raman spectra of rapidly compressed water at greater than 1 GPa at room temperature exhibits the signature of high-density amorphous ice, whereas the X-ray diffraction (XRD) pattern is dominated by crystalline ice VI. To resolve this apparent contradiction, we used molecular dynamics simulations to calculate full vibrational spectra and diffraction patterns of mixtures of vitreous ice and ice VI, including embedded interfaces between the two phases. We show quantitatively that Raman spectra, which probe the local polarizability with respect to atomic displacements, are dominated by the vitreous phase, whereas a small amount of the crystalline component is readily apparent by XRD. The results of our combined experimental and theoretical studies have implications for detecting vitreous phases of water, survival of biological systems under extreme conditions, and biological imaging. The results provide additional insight into the stable and metastable phases of H2O as a function of pressure and temperature, as well as of other materials undergoing pressure-induced amorphization and other metastable transitions.
Collapse
|
37
|
Liu J, Lan J, He X. Toward High-level Machine Learning Potential for Water Based on Quantum Fragmentation and Neural Networks. J Phys Chem A 2022; 126:3926-3936. [PMID: 35679610 DOI: 10.1021/acs.jpca.2c00601] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller-Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China.,Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Jinggang Lan
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,New York University-East China Normal University Center for Computational Chemistry, NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
38
|
Faccio C, Benzi M, Zanetti-Polzi L, Daidone I. Low- and high-density forms of liquid water revealed by a new medium-range order descriptor. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118922] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
39
|
Yoshida K, Nishimoto S, Yamaguchi T. Structural analysis of hydrazinium trifluoroacetate aqueous solution by X-ray diffraction and empirical potential structure refinement modeling in the temperature range of 25 to −125 °C. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118802] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
40
|
Gartner TE, Hunter KM, Lambros E, Caruso A, Riera M, Medders GR, Panagiotopoulos AZ, Debenedetti PG, Paesani F. Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model. J Phys Chem Lett 2022; 13:3652-3658. [PMID: 35436129 DOI: 10.1021/acs.jpclett.2c00567] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
For the past 50 years, researchers have sought molecular models that can accurately reproduce water's microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ∼208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. This work suggests that MB-pol may provide predictive capability for studies of liquid water's physical properties across broad ranges of thermodynamic states, including the so-called water's "no man's land" which is difficult to probe experimentally.
Collapse
Affiliation(s)
- Thomas E Gartner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Kelly M Hunter
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Gregory R Medders
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | | | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California, San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
41
|
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: 1] [Impact Index Per Article: 0.5] [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
|
42
|
Martelli F, Palmer JC. Signatures of sluggish dynamics and local structural ordering during ice nucleation. J Chem Phys 2022; 156:114502. [DOI: 10.1063/5.0083638] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems.
Collapse
Affiliation(s)
- Fausto Martelli
- IBM Research Europe, Hartree Centre, Daresbury WA4 4AD, United Kingdom
| | - Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, USA
| |
Collapse
|
43
|
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
|
44
|
Kringle L, Thornley WA, Kay BD, Kimmel GA. Isotope effects on the structural transformation and relaxation of deeply supercooled water. J Chem Phys 2022; 156:084501. [DOI: 10.1063/5.0078796] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
We have examined the structure of supercooled liquid D2O as a function of temperature between 185 and 255 K using pulsed laser heating to rapidly heat and cool the sample on a nanosecond timescale. The liquid structure can be represented as a linear combination of two structural motifs, with a transition between them described by a logistic function centered at 218 K with a width of 10 K. The relaxation to a metastable state, which occurred prior to crystallization, exhibited nonexponential kinetics with a rate that was dependent on the initial structural configuration. When the temperature is scaled by the temperature of maximum density, which is an isostructural point of the isotopologues, the structural transition and the non-equilibrium relaxation kinetics of D2O agree remarkably well with those for H2O.
Collapse
Affiliation(s)
- Loni Kringle
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Wyatt A. Thornley
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Bruce D. Kay
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Greg A. Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| |
Collapse
|
45
|
Multiple Melting Temperatures in Glass-Forming Melts. SUSTAINABILITY 2022. [DOI: 10.3390/su14042351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
All materials are vitrified by fast quenching even monoatomic substances. Second melting temperatures accompanied by weak exothermic or endothermic heat are often observed at Tn+ after remelting them above the equilibrium thermodynamic melting transition at Tm. These temperatures, Tn+, are due to the breaking of bonds (configurons formation) or antibonds depending on the thermal history, which is explained by using a nonclassical nucleation equation. Their multiple existence in monoatomic elements is now demonstrated by molecular dynamics simulations and still predicted. Proposed equations show that crystallization enthalpy is reduced at the temperature Tx due to new vitrification of noncrystallized parts and their melting at Tn+. These glassy parts, being equal above Tx to singular values or to their sum, are melted at various temperatures Tn+ and attain 100% in Cu46Zr46Al8 and 86.7% in bismuth. These first order transitions at Tn+ are either reversible or irreversible, depending on the formation of super atoms, either solid or liquid.
Collapse
|
46
|
Zanetti-Polzi L, Daidone I, Amadei A. A general statistical mechanical model for fluid system thermodynamics: Application to sub- and super-critical water. J Chem Phys 2022; 156:044506. [DOI: 10.1063/5.0079206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L’Aquila, via Vetoio (Coppito 1), 67010 L’Aquila, Italy
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome “Tor Vergata”, Via della Ricerca Scientifica, I-00185 Rome, Italy
| |
Collapse
|
47
|
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
|
48
|
Ibagon I, Furlan AP, Oliveira TJ, Dickman R. Phase diagram and critical properties of a two-dimensional associating lattice gas. Phys Rev E 2021; 104:064120. [PMID: 35030842 DOI: 10.1103/physreve.104.064120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Accepted: 12/01/2021] [Indexed: 06/14/2023]
Abstract
We revisit the associating lattice gas (ALG) introduced by Henriques et al. [Phys. Rev. E 71, 031504 (2005)PLEEE81539-375510.1103/PhysRevE.71.031504] in its symmetric version. In this model, defined on the triangular lattice, interaction between molecules occupying nearest-neighbor sites depends on their relative orientation, mimicking the formation of hydrogen bonds in network-forming fluids. Although all previous studies of this model agree that it has a disordered fluid (DF), a low-density liquid (LDL), and a high-density liquid (HDL) phase, quite different forms have been reported for its phase diagram. Here, we present a thorough investigation of its phase behavior using both transfer matrix calculations and Monte Carlo (MC) simulations, along with finite-size scaling extrapolations. Results in striking agreement are found using these methods. The critical point associated with the DF-HDL transition at full occupancy, identified by Furlan et al. [Phys. Rev. E 100, 022109 (2019)2470-004510.1103/PhysRevE.100.022109] is shown to be one terminus of a critical line separating these phases. In opposition to previous simulation studies, we find that the transition between the DF and LDL phases is always discontinuous, similar to the LDL-HDL transition. The associated coexistence lines meet at the point where the DF-HDL critical line ends, making it a critical-end-point. Overall, the form of the phase diagram observed in our simulations is very similar to that found in the exact solution of the model on a Husimi lattice. Our results confirm that, despite the existence of some waterlike anomalies in this model, it is unable to reproduce key features of the phase behavior of liquid water.
Collapse
Affiliation(s)
- I Ibagon
- Departamento de Física, ICEx, Universidade Federal de Minas Gerais, C.P. 702, 30123-970, Belo Horizonte, Minas Gerais-Brazil
| | - A P Furlan
- Departamento de Física, ICEx, Universidade Federal de Minas Gerais, C.P. 702, 30123-970, Belo Horizonte, Minas Gerais-Brazil
| | - T J Oliveira
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, Minas Gerais, Brazil
| | - R Dickman
- Departamento de Física and National Institute of Science and Technology for Complex Systems, ICEx, Universidade Federal de Minas Gerais, C.P. 702, 30123-970 Belo Horizonte, Minas Gerais-Brazil
| |
Collapse
|
49
|
Liu J, He X. Ab initio molecular dynamics simulation of liquid water with fragment-based quantum mechanical approach under periodic boundary conditions. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
50
|
Gallo P, Bachler J, Bove LE, Böhmer R, Camisasca G, Coronas LE, Corti HR, de Almeida Ribeiro I, de Koning M, Franzese G, Fuentes-Landete V, Gainaru C, Loerting T, de Oca JMM, Poole PH, Rovere M, Sciortino F, Tonauer CM, Appignanesi GA. Advances in the study of supercooled water. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:143. [PMID: 34825973 DOI: 10.1140/epje/s10189-021-00139-1] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 10/17/2021] [Indexed: 06/13/2023]
Abstract
In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.
Collapse
Affiliation(s)
- Paola Gallo
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy.
| | - Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Livia E Bove
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
- Sorbonne Université, CNRS UMR 7590, IMPMC, 75005, Paris, France
| | - Roland Böhmer
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Gaia Camisasca
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Luis E Coronas
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Horacio R Corti
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina
| | - Ingrid de Almeida Ribeiro
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859, Campinas, São Paulo, Brazil
| | - 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
| | - Giancarlo Franzese
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Violeta Fuentes-Landete
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Catalin Gainaru
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | | | - Peter H Poole
- Department of Physics, St. Francis Xavier University, Antigonish, NS, B2G 2W5, Canada
| | - Mauro Rovere
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
| | - Christina M Tonauer
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
| |
Collapse
|