1
|
Malosso C, Manko N, Izzo MG, Baroni S, Hassanali A. Evidence of ferroelectric features in low-density supercooled water from ab initio deep neural-network simulations. Proc Natl Acad Sci U S A 2024; 121:e2407295121. [PMID: 39083416 PMCID: PMC11317578 DOI: 10.1073/pnas.2407295121] [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: 04/12/2024] [Accepted: 06/30/2024] [Indexed: 08/02/2024] Open
Abstract
Over the last decade, an increasing body of evidence has emerged, supporting the existence of a metastable liquid-liquid critical point in supercooled water whereby two distinct liquid phases of different densities coexist. Analyzing long molecular dynamics simulations performed using deep neural-network force fields trained to accurate quantum mechanical data, we demonstrate that the low-density liquid phase displays a strong propensity toward spontaneous polarization, as witnessed by large and long-lived collective dipole fluctuations. Our findings suggest that the dynamical stability of the low-density phase, and hence the transition from high-density to low-density liquid, is triggered by a collective process involving an accumulation of rotational angular jumps, which could ignite large dipole fluctuations. This dynamical transition involves subtle changes in the electronic polarizability of water molecules which affects their rotational mobility within the two phases. These findings hold the potential for catalyzing activity in the search for dielectric-based probes of the putative second critical point.
Collapse
Affiliation(s)
- Cesare Malosso
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
| | - Natalia Manko
- Condensed Matter and Statistical Physics (CMSP), The Abdus Salam Centre for Theoretical Physics, Trieste34151, Italy
| | - Maria Grazia Izzo
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
| | - Stefano Baroni
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
- Consiglio Nazionale delle Ricerche-Istituto Officina dei Materiali, Scuola Internazionale Superiore di Studi Avanzati Unit, Trieste34136, Italy
| | - Ali Hassanali
- Condensed Matter and Statistical Physics (CMSP), The Abdus Salam Centre for Theoretical Physics, Trieste34151, Italy
| |
Collapse
|
2
|
Di Fonte N, Faccio C, Zanetti-Polzi L, Daidone I. Early prediction of spinodal-like relaxation events in supercooled liquid water. J Chem Phys 2024; 161:034505. [PMID: 39012812 DOI: 10.1063/5.0211031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Accepted: 07/01/2024] [Indexed: 07/18/2024] Open
Abstract
Several computational studies on different water models reported evidence of a phase transition in supercooled conditions between two liquid states of water differing in density: the high-density liquid (HDL) and the low-density liquid (LDL). Yet, conclusive experimental evidence of the existence of a phase transition between the two liquid water phases could not be obtained due to fast crystallization in the region where the phase transition should occur. For the same reason, the investigation of possible transition mechanisms between the two phases is committed to computational investigations. In this work, we simulate an out-of-equilibrium temperature-induced transition from the LDL to the HDL-like state in the TIP4P/2005 water model. To structurally characterize the system relaxation, we use the node total communicability (NTC) we recently proposed as an effective order parameter to discriminate the two liquid phases differing in density. We find that the relaxation process is compatible with a spinodal-like scenario. We observe the formation of HDL-like domains in the LDL phase and we characterize their fluctuating behavior and subsequent coarsening and stabilization. Furthermore, we find that the formation of stable HDL-like domains is favored in the regions where the early formation of small patches of highly connected HDL-like molecules (i.e., with very high NTC values) is observed. Besides characterizing the LDL- to HDL-like relaxation from a structural point of view, these results also show that the NTC order parameter can serve as an early-time predictor of the regions from which the transition process initiates.
Collapse
Affiliation(s)
- Nico Di Fonte
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| | - Chiara Faccio
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
- Department of Mathematics "Tullio Levi-Civita," University of Padova, Via Trieste 63, Padova 35121, Italy
| | - Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, Modena 41125, Italy
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| |
Collapse
|
3
|
Ong CY, Lee CS, Gao XY, Zhai Q, Yu Z, Shi R, Deng HY, Lam CH. Relating fragile-to-strong transition to fragile glass via lattice model simulations. Phys Rev E 2024; 109:054124. [PMID: 38907442 DOI: 10.1103/physreve.109.054124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Accepted: 04/11/2024] [Indexed: 06/24/2024]
Abstract
Glass formers are, in general, classified as strong or fragile depending on whether their relaxation rates follow Arrhenius or super-Arrhenius temperature dependence. There are, however, notable exceptions, such as water, which exhibit a fragile-to-strong (FTS) transition and behave as fragile and strong, respectively, at high and low temperatures. In this work, the FTS transition is studied using a distinguishable-particle lattice model previously demonstrated to be capable of simulating both strong and fragile glasses [C.-S. Lee, M. Lulli, L.-H. Zhang, H.-Y. Deng, and C.-H. Lam, Phys. Rev. Lett. 125, 265703 (2020)0031-900710.1103/PhysRevLett.125.265703]. Starting with a bimodal pair-interaction distribution appropriate for fragile glasses, we show that by narrowing down the energy dispersion in the low-energy component of the distribution, a FTS transition is observed. The transition occurs at a temperature at which the stretching exponent of the relaxation is minimized, in agreement with previous molecular dynamics simulations.
Collapse
Affiliation(s)
- Chin-Yuan Ong
- Department of Applied Physics, Hong Kong Polytechnic University, Hong Kong, China
| | - Chun-Shing Lee
- Department of Applied Physics, Hong Kong Polytechnic University, Hong Kong, China
- School of Science, Harbin Institute of Technology (Shenzhen), Shenzhen 518055, China
| | - Xin-Yuan Gao
- Department of Applied Physics, Hong Kong Polytechnic University, Hong Kong, China
| | - Qiang Zhai
- School of Physics, MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Xi'an Jiaotong University, Xi'an 710049, China
| | - Zhenhao Yu
- Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, School of Physics, Zhejiang University, Hangzhou 310027, China
| | - Rui Shi
- Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, School of Physics, Zhejiang University, Hangzhou 310027, China
| | - Hai-Yao Deng
- School of Physics and Astronomy, Cardiff University, 5 The Parade, Cardiff CF24 3AA, Wales, United Kingdom
| | - Chi-Hang Lam
- Department of Applied Physics, Hong Kong Polytechnic University, Hong Kong, China
| |
Collapse
|
4
|
Muthachikavil AV, Sun G, Peng B, Tanaka H, Kontogeorgis GM, Liang X. Unraveling thermodynamic anomalies of water: A molecular simulation approach to probe the two-state theory with atomistic and coarse-grained water models. J Chem Phys 2024; 160:154505. [PMID: 38624123 DOI: 10.1063/5.0194036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2023] [Accepted: 03/19/2024] [Indexed: 04/17/2024] Open
Abstract
Thermodynamic and dynamic anomalies of water play a crucial role in supporting life on our planet. The two-state theory attributes these anomalies to a dynamic equilibrium between locally favored tetrahedral structures (LFTSs) and disordered normal liquid structures. This theory provides a straightforward, phenomenological explanation for water's unique thermodynamic and dynamic characteristics. To validate this two-state feature, it is critical to unequivocally identify these structural motifs in a dynamically fluctuating disordered liquid. In this study, we employ a recently introduced structural parameter (θavg) that characterizes the local angular order within the first coordination shell to identify these LFTSs through molecular dynamics simulations. We employ both realistic water models with a liquid-liquid critical point (LLCP) and a coarse-grained water model without an LLCP to study water's anomalies in low-pressure regions below 2 kbar. The two-state theory consistently describes water's thermodynamic anomalies in these models, both with and without an LLCP. This suggests that the anomalies predominantly result from the two-state features rather than criticality, particularly within experimentally accessible temperature-pressure regions.
Collapse
Affiliation(s)
- Aswin V Muthachikavil
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Building 229, Lyngby DK-2800, Denmark
| | - Gang Sun
- Department of Physics, Center for Advanced Quantum Studies, Beijing Normal University, Beijing 100875, China
| | - Baoliang Peng
- Research Institute of Petroleum Exploration & Development (RIPED), PetroChina, Beijing 100083, China
| | - Hajime Tanaka
- Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Georgios M Kontogeorgis
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Building 229, Lyngby DK-2800, Denmark
| | - Xiaodong Liang
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Building 229, Lyngby DK-2800, Denmark
| |
Collapse
|
5
|
Donkor ED, Offei-Danso A, Rodriguez A, Sciortino F, Hassanali A. Beyond Local Structures in Critical Supercooled Water through Unsupervised Learning. J Phys Chem Lett 2024; 15:3996-4005. [PMID: 38574274 DOI: 10.1021/acs.jpclett.4c00383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/06/2024]
Abstract
The presence of a second critical point in water has been a topic of intense investigation for the last few decades. The molecular origins underlying this phenomenon are typically rationalized in terms of the competition between local high-density (HD) and low-density (LD) structures. Their identification often requires designing parameters that are subject to human intervention. Herein, we use unsupervised learning to discover structures in atomistic simulations of water close to the liquid-liquid critical point (LLCP). Encoding the information on the environment using local descriptors, we do not find evidence for two distinct thermodynamic structures. In contrast, when we deploy nonlocal descriptors that probe instead heterogeneities on the nanometer length scale, this leads to the emergence of LD and HD domains rationalizing the microscopic origins of the density fluctuations close to criticality.
Collapse
Affiliation(s)
- Edward Danquah Donkor
- The Abdus Salam International Center for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Adu Offei-Danso
- The Abdus Salam International Center for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Alex Rodriguez
- The Abdus Salam International Center for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
- Dipartimento di Matematica, Informatica e Geoscienze, Università degli studi di Trieste, via Valerio 12/1, 34127 Trieste, Italy
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, P. le Aldo Moro 5, 00185 Rome, Italy
| | - Ali Hassanali
- The Abdus Salam International Center for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
6
|
Maturi F, Raposo Filho RS, Brites CDS, Fan J, He R, Zhuang B, Liu X, Carlos LD. Deciphering Density Fluctuations in the Hydration Water of Brownian Nanoparticles via Upconversion Thermometry. J Phys Chem Lett 2024; 15:2606-2615. [PMID: 38420927 PMCID: PMC10926164 DOI: 10.1021/acs.jpclett.4c00044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 02/14/2024] [Accepted: 02/16/2024] [Indexed: 03/02/2024]
Abstract
We investigate the intricate relationship among temperature, pH, and Brownian velocity in a range of differently sized upconversion nanoparticles (UCNPs) dispersed in water. These UCNPs, acting as nanorulers, offer insights into assessing the relative proportion of high-density and low-density liquid in the surrounding hydration water. The study reveals a size-dependent reduction in the onset temperature of liquid-water fluctuations, indicating an augmented presence of high-density liquid domains at the nanoparticle surfaces. The observed upper-temperature threshold is consistent with a hypothetical phase diagram of water, validating the two-state model. Moreover, an increase in pH disrupts the organization of water molecules, similar to external pressure effects, allowing simulation of the effects of temperature and pressure on hydrogen bonding networks. The findings underscore the significance of the surface of suspended nanoparticles for understanding high- to low-density liquid fluctuations and water behavior at charged interfaces.
Collapse
Affiliation(s)
- Fernando
E. Maturi
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
- Institute
of Chemistry, São Paulo State University
(UNESP), 14800-060 Araraquara, SP, Brazil
| | - Ramon S. Raposo Filho
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Carlos D. S. Brites
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Jingyue Fan
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Ruihua He
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Bilin Zhuang
- Harvey
Mudd College, 301 Platt
Boulevard, Claremont, California 91711, United States
| | - Xiaogang Liu
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Luís D. Carlos
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| |
Collapse
|
7
|
Zielkiewicz J. Mechanism of antifreeze protein functioning and the "anchored clathrate water" concept. J Chem Phys 2023; 159:085101. [PMID: 37622597 DOI: 10.1063/5.0158590] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Accepted: 08/04/2023] [Indexed: 08/26/2023] Open
Abstract
In liquid water, there is a natural tendency to form aggregates that consist of water molecules linked by hydrogen bonds. Such spontaneously formed aggregates are surrounded by a "sea" of disordered water molecules, with both forms remaining in equilibrium. The process of creating water aggregates also takes place in the solvation water of proteins, but in this case, the interactions of water molecules with the protein surface shift the equilibrium of the process. In this paper, we analyze the structural properties of the solvation water in antifreeze proteins (AFPs). The results of molecular dynamics analysis with the use of various parameters related to the structure of solvation water on the protein surface are presented. We found that in the vicinity of the active region responsible for the binding of AFPs to ice, the equilibrium is clearly shifted toward the formation of "ice-like aggregates," and the solvation water has a more ordered ice-like structure. We have demonstrated that a reduction in the tendency to create "ice-like aggregates" results in a significant reduction in the antifreeze activity of the protein. We conclude that shifting the equilibrium in favor of the formation of "ice-like aggregates" in the solvation water in the active region is a prerequisite for the biological functionality of AFPs, at least for AFPs having a well-defined ice binding area. In addition, our results fully confirm the validity of the "anchored clathrate water" concept, formulated by Garnham et al. [Proc. Natl. Acad. Sci. U. S. A. 108, 7363 (2011)].
Collapse
Affiliation(s)
- Jan Zielkiewicz
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| |
Collapse
|
8
|
Zhang H, Wang X, Zhang J, Yu HB, Douglas JF. Approach to hyperuniformity in a metallic glass-forming material exhibiting a fragile to strong glass transition. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2023; 46:50. [PMID: 37380868 DOI: 10.1140/epje/s10189-023-00308-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2023] [Accepted: 06/05/2023] [Indexed: 06/30/2023]
Abstract
We investigate a metallic glass-forming (GF) material (Al90Sm10) exhibiting a fragile-strong (FS) glass-formation by molecular dynamics simulation to better understand this highly distinctive pattern of glass-formation in which many of the usual phenomenological relations describing relaxation times and diffusion of ordinary GF liquids no longer apply, and where instead genuine thermodynamic features are observed in response functions and little thermodynamic signature is exhibited at the glass transition temperature, Tg. Given the many unexpected similarities between the thermodynamics and dynamics of this metallic GF material with water, we first focus on the anomalous static scattering in this liquid, following recent studies on water, silicon and other FS GF liquids. We quantify the "hyperuniformity index" H of our liquid, which provides a quantitative measure of molecular "jamming". To gain insight into the T-dependence and magnitude of H, we also estimate another more familiar measure of particle localization, the Debye-Waller parameter 〈u2〉 describing the mean-square particle displacement on a timescale on the order of the fast relaxation time, and we also calculate H and 〈u2〉 for heated crystalline Cu. This comparative analysis between H and 〈u2〉 for crystalline and metallic glass materials allows us to understand the critical value of H on the order of 10-3 as being analogous to the Lindemann criterion for both the melting of crystals and the "softening" of glasses. We further interpret the emergence of FS GF and liquid-liquid phase separation in this class of liquids to arise from a cooperative self-assembly process in the GF liquid.
Collapse
Affiliation(s)
- Hao Zhang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, AB, T6G 1H9, Canada.
| | - Xinyi Wang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, AB, T6G 1H9, Canada
| | - Jiarui Zhang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, AB, T6G 1H9, Canada
| | - Hai-Bin Yu
- Wuhan National High Magnetic Field Center, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Jack F Douglas
- Material Measurement Laboratory, Material Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, MD, 20899, USA.
| |
Collapse
|
9
|
Malik S, Karmakar S, Debnath A. Quantifying dynamical heterogeneity length scales of interface water across model membrane phase transitions. J Chem Phys 2023; 158:091103. [PMID: 36889951 DOI: 10.1063/5.0137727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023] Open
Abstract
All-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine lipid membranes reveal a drastic growth in the heterogeneity length scales of interface water (IW) across fluid to ripple to gel phase transitions. It acts as an alternate probe to capture the ripple size of the membrane and follows an activated dynamical scaling with the relaxation time scale solely within the gel phase. The results quantify the mostly unknown correlations between the spatiotemporal scales of the IW and membranes at various phases under physiological and supercooled conditions.
Collapse
Affiliation(s)
- Sheeba Malik
- Department of Chemistry, IIT Jodhpur, Jodhpur, Rajasthan, India
| | - Smarajit Karmakar
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad, India
| | - Ananya Debnath
- Department of Chemistry, IIT Jodhpur, Jodhpur, Rajasthan, India
| |
Collapse
|
10
|
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
|
11
|
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
|
12
|
Jin Z, Zhao J, Chen G, Chen G, Luo Z, Xu L. Revealing the three-component structure of water with principal component analysis (PCA) of X-ray spectra. SOFT MATTER 2022; 18:7486-7496. [PMID: 36000526 DOI: 10.1039/d2sm00576j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Combining principal component analysis (PCA) of X-ray spectra with MD simulations, we experimentally reveal the existence of three basic components in water. These components exhibit distinct structures, densities, and temperature dependencies. Among the three, the two major components correspond to the low-density liquid (LDL) and the high-density liquid (HDL) predicted by the two-component model, and the third component exhibits a unique 5-hydrogen-bond configuration with ultra-high local density. As the temperature increases, the LDL component decreases and the HDL component increases, while the third component varies non-monotonically with a peak around 20 °C to 30 °C. The 3D structure of the third component is further illustrated as the uniform distribution of five hydrogen-bonded neighbors on a spherical surface. Our study reveals experimental evidence for water's possible three-component structure, which provides a fundamental basis for understanding water's special properties and anomalies.
Collapse
Affiliation(s)
- Zhipeng Jin
- Department of Physics, the Chinese University of Hong Kong, Hong Kong, China.
| | - Jiangtao Zhao
- National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230026, China.
| | - Gang Chen
- School of Physical Science and Technology, ShanghaiTech University, Shanghai 201210, China
| | - Guo Chen
- Chongqing Key Laboratory of Soft Condensed Matter Physics and Smart Materials, College of Physics, Chongqing University, Chongqing 400044, China.
| | - Zhenlin Luo
- National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230026, China.
| | - Lei Xu
- Department of Physics, the Chinese University of Hong Kong, Hong Kong, China.
| |
Collapse
|
13
|
Toffano A, Russo J, Rescigno M, Ranieri U, Bove LE, Martelli F. Temperature- and pressure-dependence of the hydrogen bond network in plastic ice VII. J Chem Phys 2022; 157:094502. [DOI: 10.1063/5.0111189] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We model, via classical molecular dynamics simulations, the plastic phase of ice VII across a wide range of the phase diagram of interest for planetary investigations. Although structural and dynamical properties of plastic ice VII are mostly independent on the thermodynamic conditions, the hydrogen bond network (HBN) acquires a diverse spectrum of topologies distinctly different from that of liquid water and of ice VII simulated at the same pressure. We observe that the HBN topology of plastic ice carries some degree of similarity with the crystal phase, stronger at thermodynamic conditions proximal to ice VII, and gradually lessening upon approaching the liquid state. Our results enrich our understanding of the properties of water at high pressure and high temperature, and may help in rationalizing the geology of
Collapse
Affiliation(s)
| | | | - Maria Rescigno
- Physics, Università degli Studi di Roma La Sapienza, Italy
| | | | | | | |
Collapse
|
14
|
Martelli F. Steady-like topology of the dynamical hydrogen bond network in supercooled water. PNAS NEXUS 2022; 1:pgac090. [PMID: 36741425 PMCID: PMC9896910 DOI: 10.1093/pnasnexus/pgac090] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 06/08/2022] [Indexed: 02/07/2023]
Abstract
We investigate the link between topology of the hydrogen bond network (HBN) and large-scale density fluctuations in water from ambient conditions to the glassy state. We observe a transition from a temperature-dependent topology at high temperatures, to a steady-like topology below the Widom temperature TW ∼ 220 K signaling the fragile-to-strong crossover and the maximum in structural fluctuations. As a consequence of the steady topology, the network suppresses large-scale density fluctuations much more efficiently than at higher temperatures. Below TW , the contribution of coordination defects of the kind A 2 D 1 (two acceptors and one donor) to the kinetics of the HBN becomes progressively more pronounced, suggesting that A 2 D 1 configurations may represent the main source of dynamical heterogeneities. Below the vitrification temperature, the freezing of rotational and translational degrees of freedom allow for an enhanced suppression of large-scale density fluctuations and the sample reaches the edges of nearly hyperuniformity. The formed network still hosts coordination defects, hence implying that nearly hyperuniformity goes beyond the classical continuous random network paradigm of tetrahedral networks and can emerge in scenarios much more complex than previously assumed. Our results unveil a hitherto undisclosed link between network topology and properties of water essential for better understanding water's rich and complex nature. Beyond implications for water, our findings pave the way to a better understanding of the physics of supercooled liquids and disordered hyperuniform networks at large.
Collapse
|
15
|
Herrero C, Pauletti M, Tocci G, Iannuzzi M, Joly L. Connection between water's dynamical and structural properties: Insights from ab initio simulations. Proc Natl Acad Sci U S A 2022; 119:e2121641119. [PMID: 35588447 PMCID: PMC9173753 DOI: 10.1073/pnas.2121641119] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 04/12/2022] [Indexed: 01/25/2023] Open
Abstract
SignificanceFirst-principles calculations, which explicitly account for the electronic structure of matter, can shed light on the molecular structure and dynamics of water in its supercooled state. In this work, we use density functional theory, which relies on a functional to describe electronic exchange and correlations, to evaluate which functional best describes the temperature evolution of bulk water transport coefficients. We also assess the validity of the Stokes-Einstein relation for all the functionals in the temperature range studied, and explore the link between structure and dynamics. Based on these results, we show how transport coefficients can be computed from structural descriptors, which require shorter simulation times to converge, and we point toward strategies to develop better functionals.
Collapse
Affiliation(s)
- Cecilia Herrero
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Michela Pauletti
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
- Institut Universitaire de France (IUF), 75005 Paris, France
| |
Collapse
|
16
|
Gao Y, Fang H, Ni K, Feng Y. Water clusters and density fluctuations in liquid water based on extended hierarchical clustering methods. Sci Rep 2022; 12:8036. [PMID: 35577839 PMCID: PMC9110331 DOI: 10.1038/s41598-022-11947-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2022] [Accepted: 04/25/2022] [Indexed: 11/29/2022] Open
Abstract
The microscopic structures of liquid water at ambient temperatures remain a hot debate, which relates with structural and density fluctuations in the hydrogen bond network. Here, we use molecular dynamics simulations of liquid water to study the properties of three-dimensional cage-like water clusters, which we investigate using extended graph-based hierarchical clustering methods. The water clusters can cover over 95% of hydrogen bond network, among which some clusters maximally encompass thousands of molecules extending beyond 3.0 nm. The clusters imply fractal behaviors forming percolating networks and the morphologies of small and large clusters show different scaling rules. The local favored clusters and the preferred connections between adjacent clusters correspond to lower energy and conformational entropy depending on cluster topologies. Temperature can destroy large clusters into small ones. We show further that the interior of clusters favors high-density patches. The water molecules in the small clusters, inside which are the void regarded as hydrophobic objects, have a preference for being more tetrahedral. Our results highlight the properties and changes of water clusters as the fundamental building blocks of hydrogen bond networks. In addition, the water clusters can elucidate structural and density fluctuations on different length scales in liquid water.
Collapse
|
17
|
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
|
18
|
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
|
19
|
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
|
20
|
Berkowicz S, Perakis F. Exploring the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes. Phys Chem Chem Phys 2021; 23:25490-25499. [PMID: 34494639 DOI: 10.1039/d1cp02866a] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The breakdown of Stokes-Einstein relation in liquid water is one of the many anomalies that take place upon cooling and indicates the decoupling of diffusion and viscosity. It is hypothesized that these anomalies manifest due to the appearance of nanometer-scale spatial fluctuations, which become increasingly pronounced in the supercooled regime. Here, we explore the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes. We capture the diffusive dynamics of the probes using dynamic light scattering and target dynamics at different length scales by varying the probe size, from ≈100 nm silica spheres to molecular-sized polyhydroxylated fullerenes (≈1 nm). We find that all the studied probes, independent of size, display similar diffusive dynamics with an Arrhenius activation energy of ≈23 kJ mol-1. Analysis of the diffusion coefficient further indicates that the probes, independent of their size, experience similar dynamic environment, which coincides with the macroscopic viscosity, while single water molecules effectively experience a comparatively lower viscosity. Finally, we conclude that our results indicate that the Stokes-Einstein relation is preserved for diffusion of probes in supercooled water T ≥ 260 K with size as small as ≈1 nm.
Collapse
Affiliation(s)
- Sharon Berkowicz
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden.
| | - Fivos Perakis
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden.
| |
Collapse
|
21
|
Corsaro C, Fazio E. From Critical Point to Critical Point: The Two-States Model Describes Liquid Water Self-Diffusion from 623 to 126 K. Molecules 2021; 26:molecules26195899. [PMID: 34641442 PMCID: PMC8512083 DOI: 10.3390/molecules26195899] [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] [Received: 09/03/2021] [Revised: 09/23/2021] [Accepted: 09/26/2021] [Indexed: 12/02/2022] Open
Abstract
Liquid’s behaviour, when close to critical points, is of extreme importance both for fundamental research and industrial applications. A detailed knowledge of the structural–dynamical correlations in their proximity is still today a target to reach. Liquid water anomalies are ascribed to the presence of a second liquid–liquid critical point, which seems to be located in the very deep supercooled regime, even below 200 K and at pressure around 2 kbar. In this work, the thermal behaviour of the self-diffusion coefficient for liquid water is analyzed, in terms of a two-states model, for the first time in a very wide thermal region (126 K < T < 623 K), including those of the two critical points. Further, the corresponding configurational entropy and isobaric-specific heat have been evaluated within the same interval. The two liquid states correspond to high and low-density water local structures that play a primary role on water dynamical behavior over 500 K.
Collapse
|
22
|
Dubey V, Dueby S, Daschakraborty S. Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective. Phys Chem Chem Phys 2021; 23:19964-19986. [PMID: 34515269 DOI: 10.1039/d1cp02202d] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.
Collapse
Affiliation(s)
- Vikas Dubey
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | - Shivam Dueby
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | | |
Collapse
|
23
|
Kikutsuji T, Kim K, Matubayasi N. Transition pathway of hydrogen bond switching in supercooled water analyzed by the Markov state model. J Chem Phys 2021; 154:234501. [PMID: 34241244 DOI: 10.1063/5.0055531] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
In this work, we examine hydrogen-bond (H-bond) switching by employing the Markov State Model (MSM). During the H-bond switching, a water hydrogen initially H-bonded with water oxygen becomes H-bonded to a different water oxygen. MSM analysis was applied to trajectories generated from molecular dynamics simulations of the TIP4P/2005 model from a room-temperature state to a supercooled state. We defined four basis states to characterize the configuration between two water molecules: H-bonded ("H"), unbound ("U"), weakly H-bonded ("w"), and alternative H-bonded ("a") states. A 16 × 16 MSM matrix was constructed, describing the transition probability between states composed of three water molecules. The mean first-passage time of the H-bond switching was estimated by calculating the total flux from the HU to UH states. It is demonstrated that the temperature dependence of the mean first-passage time is in accordance with that of the H-bond lifetime determined from the H-bond correlation function. Furthermore, the flux for the H-bond switching is decomposed into individual pathways that are characterized by different forms of H-bond configurations of trimers. The dominant pathway of the H-bond switching is found to be a direct one without passing through such intermediate states as "w" and "a," the existence of which becomes evident in supercooled water. The pathway through "w" indicates a large reorientation of the donor molecule. In contrast, the pathway through "a" utilizes the tetrahedral H-bond network, which is revealed by the further decomposition based on the H-bond number of the acceptor molecule.
Collapse
Affiliation(s)
- Takuma Kikutsuji
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kang Kim
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
24
|
Foffi R, Russo J, Sciortino F. Structural and topological changes across the liquid-liquid transition in water. J Chem Phys 2021; 154:184506. [PMID: 34241034 DOI: 10.1063/5.0049299] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
It has recently been shown that the TIP4P/Ice model of water can be studied numerically in metastable equilibrium at and below its liquid-liquid critical temperature. We report here simulations along a subcritical isotherm, for which two liquid states with the same pressure and temperature but different density can be equilibrated. This allows for a clear visualization of the structural changes taking place across the transition. We specifically focus on how the topological properties of the H-bond network change across the liquid-liquid transition. Our results demonstrate that the structure of the high-density liquid, characterized by the existence of interstitial molecules and commonly explained in terms of the collapse of the second neighbor shell, actually originates from the folding back of long rings, bringing pairs of molecules separated by several hydrogen-bonds close by in space.
Collapse
Affiliation(s)
- Riccardo Foffi
- Department of Physics, Sapienza Università di Roma, Piazzale Aldo Moro, 2, 00185 Rome, Italy
| | - John Russo
- Department of Physics, Sapienza Università di Roma, Piazzale Aldo Moro, 2, 00185 Rome, Italy
| | - Francesco Sciortino
- Department of Physics, Sapienza Università di Roma, Piazzale Aldo Moro, 2, 00185 Rome, Italy
| |
Collapse
|
25
|
|
26
|
Zhang H, Wang X, Yu HB, Douglas JF. Dynamic heterogeneity, cooperative motion, and Johari-Goldstein [Formula: see text]-relaxation in a metallic glass-forming material exhibiting a fragile-to-strong transition. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:56. [PMID: 33871722 DOI: 10.1140/epje/s10189-021-00060-7] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Accepted: 03/21/2021] [Indexed: 06/12/2023]
Abstract
We investigate the Johari-Goldstein (JG) [Formula: see text]-relaxation process in a model metallic glass-forming (GF) material ([Formula: see text]), previously studied extensively by both frequency-dependent mechanical measurements and simulation studies devoted to equilibrium properties, by molecular dynamics simulations based on validated and optimized interatomic potentials with the primary aim of better understanding the nature of this universal relaxation process from a dynamic heterogeneity (DH) perspective. The present relatively low temperature and long-time simulations reveal a direct correspondence between the JG [Formula: see text]-relaxation time [Formula: see text] and the lifetime of the mobile particle clusters [Formula: see text], defined as in previous DH studies, a relationship dual to the corresponding previously observed relationship between the [Formula: see text]-relaxation time [Formula: see text] and the lifetime of immobile particle clusters [Formula: see text]. Moreover, we find that the average diffusion coefficient D nearly coincides with [Formula: see text] of the smaller atomic species (Al) and that the 'hopping time' associated with D coincides with [Formula: see text] to within numerical uncertainty, both trends being in accord with experimental studies. This indicates that the JG [Formula: see text]-relaxation is dominated by the smaller atomic species and the observation of a direct relation between this relaxation process and rate of molecular diffusion in GF materials at low temperatures where the JG [Formula: see text]-relaxation becomes the prevalent mode of structural relaxation. As an unanticipated aspect of our study, we find that [Formula: see text] exhibits fragile-to-strong (FS) glass formation, as found in many other metallic GF liquids, but this fact does not greatly alter the geometrical nature of DH in this material and the relation of DH to dynamical properties. On the other hand, the temperature dependence of the DH and dynamical properties, such as the structural relaxation time, can be significantly altered from 'ordinary' GF liquids.
Collapse
Affiliation(s)
- Hao Zhang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta, T6G 1H9, Canada.
| | - Xinyi Wang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta, T6G 1H9, Canada
| | - Hai-Bin Yu
- Wuhan National High Magnetic Field Center, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Jack F Douglas
- Material Measurement Laboratory, Materials Science and Engineering Division, National Institute of Standards and Technology(NIST), Gaithersburg, MD, 20899, USA.
| |
Collapse
|
27
|
Zakhvataev VE, Kompaniets LA. On the existence of soliton-like collective modes in liquid water at the viscoelastic crossover. Sci Rep 2021; 11:5417. [PMID: 33686146 PMCID: PMC7940660 DOI: 10.1038/s41598-021-84277-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Accepted: 02/04/2021] [Indexed: 11/11/2022] Open
Abstract
The problem of large-density variations in supercooled and ambient water has been widely discussed in the past years. Recent studies have indicated the possibility of nanometer-sized density variations on the subpicosecond and picosecond time scales. The nature of fluctuating density heterogeneities remains a highly debated issue. In the present work, we address the problem of possible association of such density variations with the dynamics of terahertz longitudinal acoustic-like modes in liquid water. Our study is based on the fact that the subpicosecond dynamics of liquid water are essentially governed by the structural relaxation. Using a mode coupling theory approach, we found that for typical values of parameters of liquid water, the dynamic mechanism coming from the combination of the structural relaxation process and the finiteness of the amplitude of terahertz longitudinal acoustic-like mode gives rise to a soliton-like collective mode on a temperature-dependent nanometer length scale. The characteristics of this mode are consistent with the estimates of the amplitudes and temperature-dependent correlation lengths of density fluctuations in liquid water obtained in experiments and simulations. Thus, the fully dynamic mechanism could contribute to the formation and dynamics of fluctuating density heterogeneities. The soliton-like collective excitations suggested by our analysis may be relevant to different phenomena connected with supercooled water and can be expected to be associated with some ultrafast biological processes.
Collapse
Affiliation(s)
- V E Zakhvataev
- Federal Research Center "Krasnoyarsk Scientific Center of the Siberian Branch of the Russian Academy of Sciences", 660036, Krasnoyarsk, Russia.
- Siberian Federal University, 660041, Krasnoyarsk, Russia.
| | - L A Kompaniets
- Institute of Computational Modelling of the Siberian Branch of the Russian Academy of Sciences, 660036, Krasnoyarsk, Russia
| |
Collapse
|
28
|
Abstract
Water determines the properties of biological systems. Therefore, understanding the nature of the mutual interaction between water and biosystems is of primary importance for a proper assessment of any biological activity, e.g., the efficacy of new drugs or vaccines. A convenient way to characterize the interactions between biosystems and water is to analyze their impact on water density and dynamics in the proximity of the interfaces. It is commonly accepted that water bulk density and dynamical properties are recovered at distances of the order of 1 nm away from the surface of biological systems. This notion leads to the definition of hydration or biological water as the nanoscopic layer of water covering the surface of biosystems and to the expectation that all the effects of the water-interface interaction are limited to this thin region. Here, we review some of our latest contributions, showing that phospholipid membranes affect the water dynamics, structural properties, and hydrogen bond network at a distance that is more than twice as large as the commonly evoked ∼1nm thick layer and of the order of 2.4 nm. Furthermore, we unveil that at a shorter distance ∼0.5nm from the membrane, instead, there is an additional interface between lipid-bound and unbound water. Bound water has a structural role in the stability of the membrane. Our results imply that the concept of hydration water should be revised or extended and pave the way to a deeper understanding of the mutual interactions between water and biological systems.
Collapse
|
29
|
Zhang H, Wang X, Yu HB, Douglas JF. Fast dynamics in a model metallic glass-forming material. J Chem Phys 2021; 154:084505. [DOI: 10.1063/5.0039162] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Hao Zhang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 1H9, Canada
| | - Xinyi Wang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 1H9, Canada
| | - Hai-Bin Yu
- Wuhan National High Magnetic Field Center, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
| | - Jack F. Douglas
- Material Measurement Laboratory, Material Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
30
|
Mendonça BHS, Ternes P, Salcedo E, de Oliveira AB, Barbosa MC. Water diffusion in carbon nanotubes: Interplay between confinement, surface deformation, and temperature. J Chem Phys 2020; 153:244504. [PMID: 33380079 DOI: 10.1063/5.0031084] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
In this article, we investigate, through molecular dynamics simulations, the diffusion behavior of the TIP4P/2005 water confined in pristine and deformed carbon nanotubes (armchair and zigzag). To analyze different diffusive mechanisms, the water temperature was varied as 210 ≤ T ≤ 380 K. The results of our simulations reveal that water presents a non-Arrhenius to Arrhenius diffusion crossover. The confinement shifts the diffusion transition to higher temperatures when compared with the bulk system. In addition, for narrower nanotubes, water diffuses in a single line, which leads to its mobility independent of the activation energy.
Collapse
Affiliation(s)
- Bruno H S Mendonça
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Porto Alegre, RS 91501-970, Brazil
| | - Patricia Ternes
- School of Geography, University of Leeds, Leeds LS2 9NL, United Kingdom
| | - Evy Salcedo
- Coordenadoria Especial de Física, Química e Matemática, Universidade Federal de Santa Catarina, Araranguá, SC 88905-120, Brazil
| | - Alan B de Oliveira
- Departamento de Física, Universidade Federal de Ouro Preto, Ouro Preto, MG 35400-000, Brazil
| | - Marcia C Barbosa
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Porto Alegre, RS 91501-970, Brazil
| |
Collapse
|
31
|
Cheng YH, Yang HC, Chou PT. Could Chemical Reaction at the Molecular Level Show Distinction between Two Liquid-Water States? Study of the Excited-State Water-Catalyzed Proton Transfer Reaction Provides a Clue. J Phys Chem Lett 2020; 11:9468-9475. [PMID: 33108192 DOI: 10.1021/acs.jpclett.0c02896] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The two liquid-water states, which lead to some anomalies when temperature crosses over 50 ± 10 °C at the atmospheric pressure, have been continuously catching popular attention. In this study, using the excited-state proton transfer (ESPT) catalyzed by water molecules as a prototypical reaction, we demonstrate that the kinetics of ESPT indeed is influenced by the two liquid-water states. In the water-catalyzed ESPT of 3-cyano-7-azaindole (3CAI), a repetitive and comprehensive temperature-dependent study of ESPT in H2O from 0 to 90 °C shows anomalous behavior. The plot of the logarithm of ESPT rate constant as a function of inverse of absolute temperature deviates from a straight line. The convex-Arrhenius behavior manifests the activation free energy for water-assisted ESPT being dependent on temperature and hence the liquid water structure. To simplify the discussion, the plot is well fitted by using two straight lines that are crossed over in the vicinity of 40 °C. The free energy difference between water-solvated 3CAI and the 1:1 H2O:3CAI complex is deduced to be 2.29 ± 0.04 and 1.96 ± 0.04 kcal·mol-1 in the regions of 0-40 and 40-90 °C water, respectively, which also results in different frequency factors, i.e., the proton transfer/tunneling rates of (5.83 ± 0.36) × 1010 and (3.48 ± 0.27) × 1010 s-1, respectively. In a qualitative manner, the results are then rationalized by the different types of H-bonding configuration as proposed for two liquid-water phases, rendering experimental evidence to support the different water phases in ambient temperatures at 1 bar.
Collapse
Affiliation(s)
- Yu-Hsuan Cheng
- Department of Chemistry, National Taiwan University, Taipei 10617, Taiwan
| | - Hsiao-Ching Yang
- Department of Chemistry, Fu Jen Catholic University, New Taipei City 24205, Taiwan
| | - Pi-Tai Chou
- Department of Chemistry, National Taiwan University, Taipei 10617, Taiwan
| |
Collapse
|
32
|
Abstract
The origin of water's anomalies has been a matter of long-standing debate. A two-state model, dating back to Röntgen, relies on the dynamical coexistence of two types of local structures-locally favored tetrahedral structure (LFTS) and disordered normal-liquid structure (DNLS)-in liquid water. Phenomenologically, this model not only explains water's thermodynamic anomalies but also can rationalize the existence of a liquid-liquid critical point (LLCP) if there is a cooperative formation of LFTS. We recently found direct evidence for the coexistence of LFTS and DNLS in the experimental structure factor of liquid water. However, the existence of the LLCP and its impact on water's properties has remained elusive, leaving the origin of water's anomalies unclear. Here we propose a unique strategy to locate the LLCP of liquid water. First, we make a comprehensive analysis of a large set of experimental structural, thermodynamic, and dynamic data based on our hierarchical two-state model. This model predicts that the two thermodynamic and dynamical fluctuation maxima lines should cross at the LLCP if it exists, which we confirm by hundred-microsecond simulations for model waters. Based on recent experimental results of the compressibility and diffusivity measurements in the no man's land, we reveal that the two lines cross around 184 K and 173 MPa for real water, suggesting the presence of the LLCP around there. Nevertheless, we find that the criticality is almost negligible in the experimentally accessible region of liquid water because it is too far from the LLCP. Our findings would provide a clue to settle the long-standing debate.
Collapse
|
33
|
Affiliation(s)
- Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
34
|
Martelli F, Leoni F, Sciortino F, Russo J. Connection between liquid and non-crystalline solid phases in water. J Chem Phys 2020; 153:104503. [DOI: 10.1063/5.0018923] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Fausto Martelli
- IBM Research Europe, Hartree Centre, Daresbury WA4 4AD, United Kingdom
- School of Mathematics, University of Bristol, Bristol BS8 1UG, United Kingdom
| | - Fabio Leoni
- School of Mathematics, University of Bristol, Bristol BS8 1UG, United Kingdom
- Department of Physics, Sapienza University of Rome, P.le Aldo Moro 5, 00185 Rome, Italy
| | - Francesco Sciortino
- Department of Physics, Sapienza University of Rome, P.le Aldo Moro 5, 00185 Rome, Italy
| | - John Russo
- School of Mathematics, University of Bristol, Bristol BS8 1UG, United Kingdom
- Department of Physics, Sapienza University of Rome, P.le Aldo Moro 5, 00185 Rome, Italy
| |
Collapse
|
35
|
Brites CDS, Zhuang B, Debasu ML, Ding D, Qin X, Maturi FE, Lim WWY, Soh DW, Rocha J, Yi Z, Liu X, Carlos LD. Decoding a Percolation Phase Transition of Water at ∼330 K with a Nanoparticle Ruler. J Phys Chem Lett 2020; 11:6704-6711. [PMID: 32672973 DOI: 10.1021/acs.jpclett.0c02147] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Liquid water, despite its simple molecular structure, remains one of the most fascinating and complex substances. Most notably, many questions continue to exist regarding the phase transitions and anomalous properties of water, which are subtle to observe experimentally. Here, we report a sharp transition in water at 330 K unveiled through experimental measurements of the instantaneous Brownian velocity of NaYF4:Yb/Er upconversion nanoparticles in water. Our experimental investigations, corroborated by molecular dynamics simulations, elucidate a geometrical phase transition where a low-density liquid (LDL) clusters become percolated below 330 K. Around this critical temperature, we find the sizes of the LDL clusters to be similar to those of the nanoparticles, confirming the role of the upconversion nanoparticle as a powerful ruler for measuring the extensiveness of the LDL hydrogen-bond network and nanometer-scale spatial changes (20-100 nm) in liquids. Additionally, a new order parameter that unequivocally classifies water molecules into two local geometric states is introduced, providing a new tool for understanding and modeling water's many anomalous properties and phase transitions.
Collapse
Affiliation(s)
- Carlos D S Brites
- Phantom-g, CICECO - Aveiro Institute of Materials, Department of Physics, Universidade de Aveiro, 3810-193 Aveiro, Portugal
| | - Bilin Zhuang
- Institute of High Performance Computing, Singapore 138632, Singapore
- Yale-NUS College, Singapore 138527, Singapore
| | - Mengistie L Debasu
- Phantom-g, CICECO - Aveiro Institute of Materials, Department of Physics, Universidade de Aveiro, 3810-193 Aveiro, Portugal
- Department of Chemistry and CICECO - Aveiro Institute of Materials, Universidade de Aveiro, 3810-193 Aveiro, Portugal
| | - Ding Ding
- Institute of Materials Research and Engineering, Singapore 138634, Singapore
| | - Xian Qin
- Department of Chemistry, National University of Singapore, Singapore 117543, Singapore
| | - Fernando E Maturi
- Phantom-g, CICECO - Aveiro Institute of Materials, Department of Physics, Universidade de Aveiro, 3810-193 Aveiro, Portugal
| | - Winnie W Y Lim
- Institute of High Performance Computing, Singapore 138632, Singapore
| | - De Wen Soh
- Institute of High Performance Computing, Singapore 138632, Singapore
| | - J Rocha
- Department of Chemistry and CICECO - Aveiro Institute of Materials, Universidade de Aveiro, 3810-193 Aveiro, Portugal
| | - Zhigao Yi
- Department of Chemistry, National University of Singapore, Singapore 117543, Singapore
| | - Xiaogang Liu
- Department of Chemistry, National University of Singapore, Singapore 117543, Singapore
- Joint School of National University of Singapore and Tianjin University, International Campus of Tianjin University, Fuzhou 350207, P.R. China
| | - Luís D Carlos
- Phantom-g, CICECO - Aveiro Institute of Materials, Department of Physics, Universidade de Aveiro, 3810-193 Aveiro, Portugal
| |
Collapse
|
36
|
Martelli F, Crain J, Franzese G. Network Topology in Water Nanoconfined between Phospholipid Membranes. ACS NANO 2020; 14:8616-8623. [PMID: 32578978 DOI: 10.1021/acsnano.0c02984] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Water provides the driving force for the assembly and stability of many cellular components. Despite its impact on biological functions, a nanoscale understanding of the relationship between its structure and dynamics under soft confinement has remained elusive. As expected, water in contact with biological membranes recovers its bulk density and dynamics at ∼1 nm from phospholipid headgroups but surprisingly enhances its intermediate range order (IRO) over a distance, at least, twice as large. Here, we explore how the IRO is related to the water's hydrogen-bond network (HBN) and its coordination defects. We characterize the increased IRO by an alteration of the HBN up to more than eight coordination shells of hydration water. The HBN analysis emphasizes the existence of a bound-unbound water interface at ∼0.8 nm from the membrane. The unbound water has a distribution of defects intermediate between bound and bulk water, but with density and dynamics similar to bulk, while bound water has reduced thermal energy and many more HBN defects than low-temperature water. This observation could be fundamental for developing nanoscale models of biological interactions and for understanding how alteration of the water structure and topology, for example, due to changes in extracellular ions concentration, could affect diseases and signaling. More generally, it gives us a different perspective to study nanoconfined water.
Collapse
Affiliation(s)
- Fausto Martelli
- Hartree Centre, IBM Research Europe, Daresbury WA4 4AD, United Kingdom
| | - Jason Crain
- Hartree Centre, IBM Research Europe, Daresbury WA4 4AD, United Kingdom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Giancarlo Franzese
- Secció de Física Estadística i Interdisciplinària, Departament de Física de la Matèria Condensada, and Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028 Barcelona, Spain
| |
Collapse
|
37
|
Rosa ACP, Cruz C, Santana WS, Brito E, Moret MA. Non-Arrhenius behavior and fragile-to-strong transition of glass-forming liquids. Phys Rev E 2020; 101:042131. [PMID: 32422727 DOI: 10.1103/physreve.101.042131] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2020] [Accepted: 04/03/2020] [Indexed: 11/07/2022]
Abstract
Characterization of the non-Arrhenius behavior of glass-forming liquids is a broad avenue for research toward the understanding of the formation mechanisms of noncrystalline materials. In this context, this paper explores the main properties of the viscosity of glass-forming systems, considering super-Arrhenius diffusive processes. We establish the viscous activation energy as a function of the temperature, measure the degree of fragility of the system, and characterize the fragile-to-strong transition through the standard Angell's plot. Our results show that the non-Arrhenius behavior observed in fragile liquids can be understood through the non-Markovian dynamics that characterize the diffusive processes of these systems. Moreover, the fragile-to-strong transition corresponds to a change in the spatiotemporal range of correlations during the glass transition process.
Collapse
Affiliation(s)
- A C P Rosa
- Grupo de Informação Quântica e Física Estatística, Centro de Ciências Exatas e das Tecnologias, Universidade Federal do Oeste da Bahia, Rua Bertioga, 892, Morada Nobre I, 47810-059 Barreiras, Bahia, Brazil
| | - C Cruz
- Grupo de Informação Quântica e Física Estatística, Centro de Ciências Exatas e das Tecnologias, Universidade Federal do Oeste da Bahia, Rua Bertioga, 892, Morada Nobre I, 47810-059 Barreiras, Bahia, Brazil
| | - W S Santana
- Grupo de Informação Quântica e Física Estatística, Centro de Ciências Exatas e das Tecnologias, Universidade Federal do Oeste da Bahia, Rua Bertioga, 892, Morada Nobre I, 47810-059 Barreiras, Bahia, Brazil
| | - E Brito
- Grupo de Informação Quântica e Física Estatística, Centro de Ciências Exatas e das Tecnologias, Universidade Federal do Oeste da Bahia, Rua Bertioga, 892, Morada Nobre I, 47810-059 Barreiras, Bahia, Brazil
| | - M A Moret
- Programa de Modelagem Computacional-SENAI-CIMATEC, 41650-010 Salvador, Bahia, Brazil.,Universidade do Estado da Bahia (UNEB), 41150-000 Salvador, Bahia, Brazil
| |
Collapse
|
38
|
Role of hydrodynamics in liquid-liquid transition of a single-component substance. Proc Natl Acad Sci U S A 2020; 117:4471-4479. [PMID: 32051252 DOI: 10.1073/pnas.1911544117] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Liquid-liquid transition (LLT) is an unconventional transition between two liquid states in a single-component system. This phenomenon has recently attracted considerable attention not only because of its counterintuitive nature but also since it is crucial for our fundamental understanding of the liquid state. However, its physical understanding has remained elusive, particularly of the critical dynamics and phase-ordering kinetics. So far, the hydrodynamic degree of freedom, which is the most intrinsic kinetic feature of liquids, has been neglected in its theoretical description. Here we develop a Ginzburg-Landau-type kinetic theory of LLT taking it into account, based on a two-order parameter model. We examine slow critical fluctuations of the nonconserved order parameter coupled to the hydrodynamic degree of freedom in equilibrium. We also study the nonequilibrium process of LLT. We show both analytically and numerically that domain growth becomes faster (slower), depending upon the density decrease (increase) upon the transition, as a consequence of hydrodynamic flow induced by the density change. The coupling between nonconserved order parameter and hydrodynamic interaction results in anomalous domain growth in both nucleation-growth-type and spinodal-decomposition-type LLT. Our study highlights the characteristic features of hydrodynamic fluctuations and phase ordering during LLT under complex interplay among conserved and nonconserved order parameters and the hydrodynamic transport intrinsic to the liquid state.
Collapse
|
39
|
Shi R, Tanaka H. Direct Evidence in the Scattering Function for the Coexistence of Two Types of Local Structures in Liquid Water. J Am Chem Soc 2020; 142:2868-2875. [DOI: 10.1021/jacs.9b11211] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Rui Shi
- Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
40
|
Dubey V, Erimban S, Indra S, Daschakraborty S. Understanding the Origin of the Breakdown of the Stokes-Einstein Relation in Supercooled Water at Different Temperature-Pressure Conditions. J Phys Chem B 2019; 123:10089-10099. [PMID: 31702917 DOI: 10.1021/acs.jpcb.9b08309] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A recent experiment has measured the viscosity of water down to approximately 244 K and up to 300 MPa. The correct viscosity and translational diffusivity data at various temperature-pressure (T-P) state points allowed for checking the validity of the Stokes-Einstein (SE) relation, which accounts for the coupling between translational self-diffusion and medium viscosity. The diffusion-viscosity decoupling increases with decreasing temperature, but the increasing pressure reduces the extent of the decoupling. Earlier simulation studies explained the breakdown of the SE relation in terms of the location of the Widom line, emanating from the liquid-liquid critical point (LLCP). Although these studies made a significant contribution to the current understanding of the above phenomena, a detailed molecular picture is still lacking. Recently, our group has explained the diffusion-viscosity decoupling from a jump-diffusion perspective. The jump-diffusion coefficient, emanating from the jump translation of water molecules, is calculated using a quantitative approach for different temperatures at ambient pressure. It has been observed that jump-diffusion is the key factor for diffusion-viscosity decoupling in supercooled water. The same method is adopted in the present work to estimate the jump-diffusion coefficient for different T-P state points and, thereby, explains the role of jump-diffusion for the different extents of the SE relation breakdown at different pressures. The residual diffusion coefficient, the other component of the total diffusion that originates from small step displacement and that is calculated by subtracting the jump-diffusion coefficient from the total diffusion, is seen to be fairly coupled to the viscosity at the entire range of temperature and pressure. Furthermore, we have calculated the average number of H-bonds per water molecule and the tetrahedral order for different T-P state points and investigated an approximate correlation between the average local structure and the contribution of the jump-diffusion to the total diffusion of water. This study, therefore, puts forward a new perspective for explaining the SE relation breakdown in supercooled water under different pressure conditions.
Collapse
Affiliation(s)
- Vikas Dubey
- Department of Chemistry , Indian Institute of Technology Patna , Patna , Bihar 801106 , India
| | - Shakkira Erimban
- Department of Chemistry , Indian Institute of Technology Patna , Patna , Bihar 801106 , India
| | - Sandipa Indra
- Department of Chemistry , Indian Institute of Technology Patna , Patna , Bihar 801106 , India
| | - Snehasis Daschakraborty
- Department of Chemistry , Indian Institute of Technology Patna , Patna , Bihar 801106 , India
| |
Collapse
|
41
|
Leoni F, Shi R, Tanaka H, Russo J. Crystalline clusters in mW water: Stability, growth, and grain boundaries. J Chem Phys 2019; 151:044505. [DOI: 10.1063/1.5100812] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Fabio Leoni
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | - Rui Shi
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - John Russo
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
42
|
Caupin F, Anisimov MA. Thermodynamics of supercooled and stretched water: Unifying two-structure description and liquid-vapor spinodal. J Chem Phys 2019; 151:034503. [PMID: 31325919 DOI: 10.1063/1.5100228] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
We have applied a two-structure approach to the description of the thermodynamic properties of supercooled and stretched water, metastable toward vapor, ice, or both, by incorporating the stability limit of liquid with respect to vapor at negative pressures. In addition to the properties of water considered in previous studies, we include new data recently obtained in deeply supercooled and stretched regions. Our model reproduces the experimentally observed anomalies in metastable water up to 400 MPa and down to -140 MPa, and can provide a physically based extrapolation in regions where no measurements are available yet. Moreover, we are able to elucidate the thermodynamic nature of the alternative "states" of liquid water, namely, high-temperature denser water (state A) and "mother-of-ice" lighter water (state B). Based on the internal consistency of the described anomalies and new data on the isothermal compressibility, we exclude the critical-point-free scenario in which the first-order liquid-liquid transition line would continue into the stretched liquid state (doubly metastable) crossing the vapor-liquid spinodal. A "singularity-free" scenario remains an option for explaining supercooled water's anomalies within the framework of two-state thermodynamics; however, the extreme case of the singularity-free scenario, ideal mixing of A and B, seems improbable. We have also clarified the concept of fast interconversion of alternative states in supercooled water as a phenomenological representation of distribution of short-ranged local structures.
Collapse
Affiliation(s)
- Frédéric Caupin
- Institut Lumière Matière, Université Claude Bernard Lyon 1, CNRS, Université de Lyon, F-69622, Villeurbanne, France
| | - Mikhail A Anisimov
- Department of Chemical and Biomolecular Engineering and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
43
|
Camisasca G, Schlesinger D, Zhovtobriukh I, Pitsevich G, Pettersson LGM. A proposal for the structure of high- and low-density fluctuations in liquid water. J Chem Phys 2019; 151:034508. [PMID: 31325915 DOI: 10.1063/1.5100875] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Based on recent experimental data that can be interpreted as indicating the presence of specific structures in liquid water, we build and optimize two structural models which we compare with the available experimental data. To represent the proposed high-density liquid structures, we use a model consisting of chains of water molecules, and for low-density liquid, we investigate fused dodecahedra as templates for tetrahedral fluctuations. The computed infrared spectra of the models are in very good agreement with the extracted experimental spectra for the two components, while the extracted structures from molecular dynamics (MD) simulations give spectra that are intermediate between the experimentally derived spectra. Computed x-ray absorption and emission spectra as well as the O-O radial distribution functions of the proposed structures are not contradicted by experiment. The stability of the proposed dodecahedral template structures is investigated in MD simulations by seeding the starting structure, and remnants found to persist on an ∼30 ps time scale. We discuss the possible significance of such seeds in simulations and whether they can be viable candidates as templates for structural fluctuations below the compressibility minimum of liquid water.
Collapse
Affiliation(s)
- Gaia Camisasca
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Daniel Schlesinger
- Department of Environmental Science and Analytical Chemistry & Bolin Centre for Climate Research, Stockholm University, 114 18 Stockholm, Sweden
| | - Iurii Zhovtobriukh
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - George Pitsevich
- Belarusian State University, Nezavisimosti Ave., 4, 220030 Minsk, Belarus
| | - Lars G M Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| |
Collapse
|
44
|
Shi R, Tanaka H. Homogeneous nucleation of ferroelectric ice crystal driven by spontaneous dipolar ordering in supercooled TIP5P water. J Chem Phys 2019; 151:024501. [DOI: 10.1063/1.5100634] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Rui Shi
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
45
|
Abstract
The idea that water is a mixture of two distinct states is analyzed in some detail. It is shown that the known compressibility of water is in fact sufficiently small that for a volume of water of size 1 nm3, the density fluctuations are of order 4% of the average density. This is much smaller than the ≈25% density fluctuations that would be required for significant regions of high and low density water to occur on this volume scale. It is also pointed out that the density fluctuations in water are, if anything, smaller than those that occur in other common liquids which do not have the anomalous properties of water. It is shown that if the distribution of density fluctuations is unimodal, the system is in the one-phase region, and if bimodal, it is in the two-phase region. None of the liquid or amorphous phases of water explored in this work give any sign of being in the two-phase region. Existing neutron and X-ray scattering data on water in the amorphous phases, and in the stable liquid phases as a function pressure and temperature, are subject to a new set of empirical potential structure refinement simulations. These simulations are interrogated for their configurational entropy, using a spherical harmonic reconstruction of the full orientational pair correlation function. It is shown that the excess pair entropy derived from this function, plus the known perfect gas contributions, give a reasonable account of the total entropy of water, within the likely errors. This estimated entropy follows the expected declining trend with decreasing temperature. Evidence that higher density water will have higher entropy than lower density water emerges, in accordance with what is expected from the negative thermal expansion coefficient of water at low temperatures. However, this entropy increase is not large and goes through a maximum before declining at yet higher densities and pressures, in a manner reminiscent of what has been previously observed in the diffusion coefficient as a function of pressure. There is no evidence that ambient water can be regarded as patches of high density, high entropy and low density, low entropy liquid, as some have claimed, since high density water has a similar entropy to low density water. There is some evidence that the distinction between these two states will become more pronounced as the temperature is lowered. Extensive discussion of the use of order parameters to describe water structure is given, and it is pointed out that these indices generally cannot be used to infer two-state behavior.
Collapse
Affiliation(s)
- A K Soper
- ISIS Facility, UKRI-STFC Rutherford Appleton Laboratory, Harwell Campus, Didcot, Oxfordshire OX11 0QX, United Kingdom
| |
Collapse
|
46
|
Kikutsuji T, Kim K, Matubayasi N. Diffusion dynamics of supercooled water modeled with the cage-jump motion and hydrogen-bond rearrangement. J Chem Phys 2019; 150:204502. [DOI: 10.1063/1.5095978] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Takuma Kikutsuji
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kang Kim
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
47
|
Martelli F. Unravelling the contribution of local structures to the anomalies of water: The synergistic action of several factors. J Chem Phys 2019; 150:094506. [PMID: 30849899 DOI: 10.1063/1.5087471] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the microscopic origin of water's anomalies by inspecting the hydrogen bond network (HBN) and the spatial organization of low-density-liquid (LDL) like and high-density-liquid (HDL) like environments. Specifically, we simulate-via classical molecular dynamics simulations-the isobaric cooling of a sample composed of 512 water molecules from ambient to deeply undercooled conditions at three pressures, namely, 1 bar, 400 bars, and 1000 bars. In correspondence with the Widom line (WL), (i) the HDL-like dominating cluster undergoes fragmentation caused by the percolation of LDL-like aggregates following a spinodal-like kinetics; (ii) such fragmentation always occurs at a "critical" concentration of ∼20%-30% in LDL; (iii) the HBN within LDL-like environments is characterized by an equal number of pentagonal and hexagonal rings that create a state of maximal frustration between a configuration that promotes crystallization (hexagonal ring) and a configuration that hinders it (pentagonal ring); (iv) the spatial organization of HDL-like environments shows a marked variation. Moreover, the inspection of the global symmetry shows that the intermediate-range order decreases in correspondence with the WL and such a decrease becomes more pronounced upon increasing the pressure, hence supporting the hypothesis of a liquid-liquid critical point. Our results reveal and rationalize the complex microscopic origin of water's anomalies as the cooperative effect of several factors acting synergistically. Beyond implications for water, our findings may be extended to other materials displaying anomalous behaviours.
Collapse
Affiliation(s)
- Fausto Martelli
- IBM Research, Hartree Centre, Daresbury WA4 4AD, United Kingdom
| |
Collapse
|
48
|
Pettersson LGM. A Two-State Picture of Water and the Funnel of Life. SPRINGER PROCEEDINGS IN PHYSICS 2019. [DOI: 10.1007/978-3-030-21755-6_1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|