1
|
Saito S. Unraveling the dynamic slowdown in supercooled water: The role of dynamic disorder in jump motions. J Chem Phys 2024; 160:194506. [PMID: 38767263 DOI: 10.1063/5.0209713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Accepted: 05/01/2024] [Indexed: 05/22/2024] Open
Abstract
When a liquid is rapidly cooled below its melting point without inducing crystallization, its dynamics slow down significantly without noticeable structural changes. Elucidating the origin of this slowdown has been a long-standing challenge. Here, we report a theoretical investigation into the mechanism of the dynamic slowdown in supercooled water, a ubiquitous yet extraordinary substance characterized by various anomalous properties arising from local density fluctuations. Using molecular dynamics simulations, we found that the jump dynamics, which are elementary structural change processes, deviate from Poisson statistics with decreasing temperature. This deviation is attributed to slow variables competing with the jump motions, i.e., dynamic disorder. The present analysis of the dynamic disorder showed that the primary slow variable is the displacement of the fourth nearest oxygen atom of a jumping molecule, which occurs in an environment created by the fluctuations of molecules outside the first hydration shell. As the temperature decreases, the jump dynamics become slow and intermittent. These intermittent dynamics are attributed to the prolonged trapping of jumping molecules within extended and stable low-density domains. As the temperature continues to decrease, the number of slow variables increases due to the increased cooperative motions. Consequently, the jump dynamics proceed in a higher-dimensional space consisting of multiple slow variables, becoming slower and more intermittent. It is then conceivable that with further decreasing temperature, the slowing and intermittency of the jump dynamics intensify, eventually culminating in a glass transition.
Collapse
Affiliation(s)
- Shinji Saito
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan and The Graduate University for Advanced Studies (SOKENDAI), Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
2
|
Niinomi H, Kouch A, Hama T, Nada H, Yamazaki T, Kimura Y. Low- and High-Density Unknown Waters at Ice-Water Interfaces. J Phys Chem Lett 2022; 13:4251-4256. [PMID: 35543729 DOI: 10.1021/acs.jpclett.2c00660] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Experimental confirmation of liquid polymorphs of water, high-density liquid (HDL) and low-density liquid (LDL), is desired for understanding not only the liquid state of matter but also the origin of the mysterious properties of water. However, this remains challenging because the liquid-liquid critical point of water lies in experimentally inaccessible supercooling conditions known as "no man's land". Here, we show by in situ optical microscopy that droplets and layers of low- and high-density unknown waters (LDUW and HDUW) appear macroscopically depending upon ice polymorphs at non-equilibrium interfaces between water and ices under experimentally accessible (de)pressurization conditions. These unknown waters were found to have characteristic velocities (about 20 and 100 m/s for LDUW and HDUW, respectively) different from water (about 40 m/s) and quasi-liquid layers (QLLs) (about 2 and 0.2 m/s for droplet and layer forms of QLLs, respectively). Our discoveries provide insight on liquid polymorphism of water.
Collapse
Affiliation(s)
- Hiromasa Niinomi
- Institute of Low Temperature Science, Hokkaido University, Kita-19, Nishi-8, Kita-ku, Sapporo 060-0819, Japan
- Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Akira Kouch
- Institute of Low Temperature Science, Hokkaido University, Kita-19, Nishi-8, Kita-ku, Sapporo 060-0819, Japan
| | - Tetsuya Hama
- Komaba Institute for Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
| | - Hiroki Nada
- National Institute of Advanced Industrial Science and Technology (AIST), 16-1 Onogawa, Tsukuba 305-8569, Japan
| | - Tomoya Yamazaki
- Institute of Low Temperature Science, Hokkaido University, Kita-19, Nishi-8, Kita-ku, Sapporo 060-0819, Japan
| | - Yuki Kimura
- Institute of Low Temperature Science, Hokkaido University, Kita-19, Nishi-8, Kita-ku, Sapporo 060-0819, Japan
| |
Collapse
|
3
|
Li P, Lu H, Fu YQ. Phase transition of supercooled water confined in cooperative two-state domain. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:165403. [PMID: 35114647 DOI: 10.1088/1361-648x/ac519b] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Accepted: 02/03/2022] [Indexed: 06/14/2023]
Abstract
The question of 'what is the structure of water?' has been regarded as one of the major scientific conundrums in condensed-matter physics due to the complex phase behavior and condensed structure of supercooled water. Great effort has been made so far using both theoretical analysis based on various mathematical models and computer simulations such as molecular dynamics and first-principle. However, these theoretical and simulation studies often do not have strong evidences of condensed-matter physics to support. In this study, a cooperative domain model is formulated to describe the dynamic phase transition of supercooled water between supercooled water and amorphous ice, both of which are composed of low- and high-density liquid water. Free volume theory is initially employed to identify the working principle of dynamic phase transition and its connection to glass transition in the supercooled water. Then a cooperative two-state model is developed to characterize the dynamic anomalies of supercooled water, including density, viscosity and self-diffusion coefficient. Finally, the proposed model is verified using the experimental results reported in literature.
Collapse
Affiliation(s)
- Peizhao Li
- Science and Technology on Advanced Composites in Special Environments Laboratory, Harbin Institute of Technology, Harbin 150080, People's Republic of China
| | - Haibao Lu
- Science and Technology on Advanced Composites in Special Environments Laboratory, Harbin Institute of Technology, Harbin 150080, People's Republic of China
| | - Yong-Qing Fu
- Faculty of Engineering and Environment, Northumbria University, Newcastle upon Tyne, NE1 8ST, United Kingdom
| |
Collapse
|
4
|
Oliveira ASF, Ciccotti G, Haider S, Mulholland AJ. Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:144. [PMID: 34720710 PMCID: PMC8549953 DOI: 10.1140/epjb/s10051-021-00157-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 07/05/2021] [Indexed: 05/03/2023]
Abstract
ABSTRACT A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in β -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol, BS8 1TQ UK
| | - Giovanni Ciccotti
- Institute for Applied Computing “Mauro Picone” (IAC), CNR, Via dei Taurini 19, 00185 Rome, Italy
- School of Physics, University College of Dublin, UCD-Belfield, Dublin 4, Ireland
- Università di Roma La Sapienza, Ple. A. Moro 5, 00185 Rome, Italy
| | - Shozeb Haider
- School of Pharmacy, University College London, London, WC1N 1AX UK
| | - Adrian J. Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
| |
Collapse
|
5
|
Montero de Hijes P, Sanz E, Joly L, Valeriani C, Caupin F. Viscosity and self-diffusion of supercooled and stretched water from molecular dynamics simulations. J Chem Phys 2018; 149:094503. [DOI: 10.1063/1.5042209] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Pablo Montero de Hijes
- Departamento de Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, Madrid 28040, Spain
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Laurent Joly
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Lyon, France
| | - Chantal Valeriani
- Departamento de Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, Madrid 28040, Spain
| | - Frédéric Caupin
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Lyon, France
| |
Collapse
|
6
|
Galamba N. On the hydrogen-bond network and the non-Arrhenius transport properties of water. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:015101. [PMID: 27831934 DOI: 10.1088/0953-8984/29/1/015101] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We study the structural and dynamic transformations of SPC/E water with temperature, through molecular dynamics (MD), and discuss the non-Arrhenius behavior of the transport properties and orientational dynamics, and the magnitude of the breakdown of the Stokes-Einstein (SE) and the Stokes-Einstein-Debye (SED) relations, in the light of these transformations. Our results show that deviations from Arrhenius behavior of the self-diffusion at low temperatures cannot be exclusively explained by the reduction of water defects (interstitial waters) and the increase of the local tetrahedrality, thus, suggesting the importance of the slowdown of collective rearrangements. Interestingly we find that at high temperatures (T ⩾ 340 K) water defects lead to a slight increase of the tetrahedrality and a decrease of the self-diffusion, opposite to water at low temperatures. The relative magnitude of the breakdown of the SE and the SED relations is found to be in accord with recent experiments (Dehaoui et al 2015 Proc. Natl Acad. Sci. USA 112 12020) resolving the discrepancy with previous MD results. Further, we show that SPC/E hydrogen-bond (HB) lifetimes deviate from Arrhenious behaviour at low temperatures in contrast with some previous MD studies. This deviation is nevertheless much smaller than that observed for the orientational dynamics and the transport properties of water, consistent with the relaxation times measured by several experimental methods. The HB acceptor exchange dynamics defined here by the acceptor switch and reform (librational dynamics) frequencies exhibit similar Arrhenius deviations, thus explaining to some extent the non-Arrhenius behavior of the transport properties and of the orientational dynamics of water. Our results also show that the fraction of HB switches through a bifurcated pathway follow a power law with the temperature decrease. Thus, at low temperatures HB acceptor switches are less frequent but occur on a faster time scale consistent with the temperature dependence of the ratio of the rotational relaxation times for the different Legendre polynomial ranks.
Collapse
Affiliation(s)
- N Galamba
- Centro de Química e Bioquímica da Faculdade de Ciências da Universidade de Lisboa, Edifício C8, Campo Grande, 1749-016 Lisboa, Portugal
| |
Collapse
|
7
|
Rogers TR, Leong KY, Wang F. Possible Evidence for a New Form of Liquid Buried in the Surface Tension of Supercooled Water. Sci Rep 2016; 6:33284. [PMID: 27615518 PMCID: PMC5018856 DOI: 10.1038/srep33284] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2016] [Accepted: 08/24/2016] [Indexed: 01/29/2023] Open
Abstract
Contrary to the historical data, several recent experiments indicate that the surface tension of supercooled water follows a smooth extrapolation of the IAPWS equation in the supercooled regime. It can be seen, however, that a small deviation from the IAPWS equation is present in the recent experimental measurements. It is shown with simulations using the WAIL water potential that the small deviation in the experimental data is consistent with the tail of an exponential growth in surface tension as temperature decreases. The emergence temperature, Te, of a substantial deviation from the IAPWS equation is shown to be 227 K for the WAIL water and 235 K for real water. Since the 227 K Te is close to the Widom line in WAIL water, we argue that real water at 235 K approaches a similar crossover line at one atmospheric pressure.
Collapse
Affiliation(s)
- T. Ryan Rogers
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, USA
| | - Kai-Yang Leong
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, USA
| | - Feng Wang
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, USA
| |
Collapse
|