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
|
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
|
3
|
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
|
4
|
Madl P, Renati P. Quantum Electrodynamics Coherence and Hormesis: Foundations of Quantum Biology. Int J Mol Sci 2023; 24:14003. [PMID: 37762305 PMCID: PMC10530466 DOI: 10.3390/ijms241814003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2023] [Revised: 07/12/2023] [Accepted: 07/21/2023] [Indexed: 09/29/2023] Open
Abstract
BACKGROUND "Quantum biology" (QB) is a promising theoretical approach addressing questions about how living systems are able to unfold dynamics that cannot be solved on a chemical basis or seem to violate some fundamental laws (e.g., thermodynamic yield, morphogenesis, adaptation, autopoiesis, memory, teleology, biosemiotics). Current "quantum" approaches in biology are still very basic and "corpuscular", as these rely on a semi-classical and approximated view. We review important considerations of theory and experiments of the recent past in the field of condensed matter, water, physics of living systems, and biochemistry to join them by creating a consistent picture applicable for life sciences. Within quantum field theory (QFT), the field (also in the matter field) has the primacy whereby the particle, or "quantum", is a derivative of it. The phase of the oscillation and not the number of quanta is the most important observable of the system. Thermodynamics of open systems, symmetry breaking, fractals, and quantum electrodynamics (QED) provide a consistent picture of condensed matter, liquid water, and living matter. Coherence, resonance-driven biochemistry, and ion cyclotron resonance (Liboff-Zhadin effect) emerge as crucial hormetic phenomena. We offer a paradigmatic approach when dealing with living systems in order to enrich and ultimately better understand the implications of current research activities in the field of life sciences.
Collapse
Affiliation(s)
- Pierre Madl
- Department of Biosciences & Medical Biology, University of Salzburg, A-5020 Salzburg, Austria
- Prototyping Unit, Edge-Institute, ER-System Mechatronics, A-5440 Golling, Austria
| | - Paolo Renati
- Prototyping Unit, Edge-Institute, ER-System Mechatronics, A-5440 Golling, Austria
- The World Water Community, Marconistraat 16, 3029 AK Rotterdam, The Netherlands;
| |
Collapse
|
5
|
Daidone I, Foffi R, Amadei A, Zanetti-Polzi L. A statistical mechanical model of supercooled water based on minimal clusters of correlated molecules. J Chem Phys 2023; 159:094502. [PMID: 37655770 DOI: 10.1063/5.0157505] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 08/14/2023] [Indexed: 09/02/2023] Open
Abstract
In this paper, we apply a theoretical model for fluid state thermodynamics to investigate simulated water in supercooled conditions. This model, which we recently proposed and applied to sub- and super-critical fluid water [Zanetti-Polzi et al., J. Chem. Phys. 156(4), 44506 (2022)], is based on a combination of the moment-generating functions of the enthalpy and volume fluctuations as provided by two gamma distributions and provides the free energy of the system as well as other relevant thermodynamic quantities. The application we make here provides a thermodynamic description of supercooled water fully consistent with that expected by crossing the liquid-liquid Widom line, indicating the presence of two distinct liquid states. In particular, the present model accurately reproduces the Widom line temperatures estimated with other two-state models and well describes the heat capacity anomalies. Differently from previous models, according to our description, a cluster of molecules that extends beyond the first hydration shell is necessary to discriminate between the statistical fluctuation regimes typical of the two liquid states.
Collapse
Affiliation(s)
- Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| | - Riccardo Foffi
- Department of Civil, Environmental and Geomatic Engineering, Institute of Environmental Engineering, ETH Zürich, Laura-Hezner-Weg 7, 8093 Zürich, Switzerland
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome "Tor Vergata," Via della Ricerca Scientifica, I-00185 Rome, Italy
| | - Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| |
Collapse
|
6
|
Kringle L, Kay BD, Kimmel GA. Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water. J Phys Chem B 2023; 127:3919-3930. [PMID: 37097190 DOI: 10.1021/acs.jpcb.3c01465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.
Collapse
Affiliation(s)
- Loni Kringle
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Bruce D Kay
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Greg A Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|
7
|
de Oliveira PMC, de Souza JIR, da Silva JAB, Longo RL. Temperature Dependence of Hydrogen Bond Networks of Liquid Water: Thermodynamic Properties and Structural Heterogeneity from Topological Descriptors. J Phys Chem B 2023; 127:2250-2257. [PMID: 36877152 DOI: 10.1021/acs.jpcb.2c08873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Abstract
Topological analyses of hydrogen bond networks were performed based on the complex network and island statistics of liquid water at different temperatures. The influence of temperature on the liquid water structures and the topological properties of the hydrogen bond networks was investigated by Metropolis Monte Carlo simulations with the TIP4P/2005 potential model. The bilinear behavior of the second peak in the radial distribution function with the temperature was properly reproduced by these simulations. The average connectivity also displayed a bilinear behavior consistent with being a local descriptor. The semiglobal average path length (or geodesic distance) descriptor showed an unprecedented trimodal distribution, whose areas were dependent on the temperature. Considering equilibrium between these three sets of networks, standard enthalpy and entropy of equilibrium were determined for the first time, providing new insights into the structural heterogeneities of liquid water with interesting perspectives for modeling these quantitative properties of hydrogen bond networks.
Collapse
Affiliation(s)
- Paulo M C de Oliveira
- Núcleo Interdisciplinar de Ciências Exatas e da Natureza, Campus do Agreste. Universidade Federal de Pernambuco, 55.014-900 Caruaru, Pernambuco, Brazil
| | - Jéssica I R de Souza
- Programa de Pós-Graduação em Ciência de Materiais, Universidade Federal de Pernambuco, 50740-560 Recife, Pernambuco, Brazil
| | - Juliana A B da Silva
- Núcleo Interdisciplinar de Ciências Exatas e da Natureza, Campus do Agreste. Universidade Federal de Pernambuco, 55.014-900 Caruaru, Pernambuco, Brazil.,Programa de Pós-Graduação em Química, Universidade Federal Rural de Pernambuco, 52.171-900 Recife, Pernambuco, Brazil
| | - Ricardo L Longo
- Programa de Pós-Graduação em Ciência de Materiais, Universidade Federal de Pernambuco, 50740-560 Recife, Pernambuco, Brazil.,Departamento de Química Fundamental, Universidade Federal de Pernambuco, 50740-540 Recife, Pernambuco, Brazil
| |
Collapse
|
8
|
Lanza G. Water model for hydrophobic cavities: structure and energy from quantum-chemical calculations. Phys Chem Chem Phys 2023; 25:6902-6913. [PMID: 36799662 DOI: 10.1039/d2cp05195h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
This ab initio study aims to design a series of large water clusters having a hollow clathrate-like cage able to host hydrophobic solutes of various sizes. Starting from the (H2O)n (n = 18, 20, 24 and 28) hollow cages, water layers have been added in a stepwise manner in order to model the configuration of water molecules beyond the primary shell. The large (H2O)100, (H2O)120 and (H2O)140 clusters complete the hydrogen bonding network of the cage with optimal and regular tiling of the do-, tetra-decahedron and hexa-decahedron, respectively. This study is corroborated by an investigation of dense water clusters up to the (H2O)123 one, being highly consistent with experimental data on ice concerning the electronic and zero-point energies for aggregate formation at 0 K and enthalpy and entropy at 273 K. The cavity creation profoundly alters the orientation of water molecules compared with those found in dense clusters. Nevertheless, such a large reorganization is necessary to maximize the water-water attraction by making it similar to the one found in dense clusters. The cage formation is an endothermic process; however, the computed values are large compared with previous reports for hydrocarbon aqueous solutions. Larger clusters are required for a more fruitful comparison.
Collapse
Affiliation(s)
- Giuseppe Lanza
- Dipartimento di Scienze del Farmaco e della Salute, Università di Catania, Viale A. Doria 6, Catania, 95125, Italy.
| |
Collapse
|
9
|
Foffi R, Sciortino F. Correlated Fluctuations of Structural Indicators Close to the Liquid-Liquid Transition in Supercooled Water. J Phys Chem B 2022; 127:378-386. [PMID: 36538764 PMCID: PMC9841516 DOI: 10.1021/acs.jpcb.2c07169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Multiple numerical studies have unambiguously shown the existence of a liquid-liquid critical point in supercooled states for different numerical models of water, and various structural indicators have been put forward to describe the transformation associated with this phase transition. Here we analyze numerical simulations of near-critical supercooled water to compare the behavior of several of such indicators with critical density fluctuations. We show that close to the critical point most indicators are strongly correlated to density, and some of them even display identical distributions of fluctuations. These indicators probe the exact same free energy landscape, therefore providing a thermodynamic description of critical supercooled water which is identical to that provided by the density order parameter. This implies that close to the critical point, there is a tight coupling between many, only apparently distinct, structural degrees of freedom.
Collapse
Affiliation(s)
- Riccardo Foffi
- Institute
for Environmental Engineering, Department of Civil, Environmental
and Geomatic Engineering, ETH Zürich, 8093Zürich, Switzerland
| | - Francesco Sciortino
- Dipartimento
di Fisica, Sapienza Università di
Roma, Piazzale Aldo Moro
5, I-00185Rome, Italy,E-mail:
| |
Collapse
|
10
|
Camisasca G, Tenuzzo L, Gallo P. Protein hydration water: focus on low density and high density local structures upon cooling. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
11
|
Škvára J, Nezbeda I. Thermodynamics and structure of supercooled water. II. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120508] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
12
|
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
|
13
|
Yun Y, Khaliullin RZ, Jung Y. Correlated Local Fluctuations in the Hydrogen Bond Network of Liquid Water. J Am Chem Soc 2022; 144:13127-13136. [PMID: 35820142 DOI: 10.1021/jacs.2c02362] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The hypothesis that liquid water can separate into two phases in the supercooled state has been supported by recent experimental and theoretical studies. However, whether such structural inhomogeneity extends to ambient conditions is under intense debate. Due to the dynamic nature of the hydrogen bond network of liquid water, exploring its structure requires detailed insight into the collective motion of neighboring water molecules, a missing link that has not been examined so far. Here, highly sensitive quantum mechanical calculations detect that the time evolution of nearby hydrogen bonds is strongly correlated, revealing a direct mechanism for the appearance of short-range structural fluctuations in the hydrogen bond network of liquid water for the first time. This correlated dynamics is found to be closely connected to the static structural picture. The distortions from the tetrahedral structure do not occur independently but are correlated due to the preference of nearby donors and acceptors to be in similar environments. The existence of such cooperative fluctuations is further supported by the temperature dependence of the local structural evolution and explained by conventional analysis of localized orbitals. It was found that such correlated structural fluctuations are only observed on a short length scale in simulations at ambient conditions. The correlations of the nearby hydrogen bond pairs of liquid water unveiled here are expected to offer a new insight into connecting the dynamics of individual water molecules and the local structure of the hydrogen bond network.
Collapse
Affiliation(s)
- Yonghwan Yun
- Department of Chemical and Biomolecular Engineering (BK21 four), Korea Advanced Institute of Science and Technology (KAIST), 291 Daehak-ro, Daejeon 34141, South Korea
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, 801 Sherbrooke Street West, Montreal, Quebec H3A 0B8, Canada
| | - Yousung Jung
- Department of Chemical and Biomolecular Engineering (BK21 four), Korea Advanced Institute of Science and Technology (KAIST), 291 Daehak-ro, Daejeon 34141, South Korea
| |
Collapse
|
14
|
Offei-Danso A, Hassanali A, Rodriguez A. High-Dimensional Fluctuations in Liquid Water: Combining Chemical Intuition with Unsupervised Learning. J Chem Theory Comput 2022; 18:3136-3150. [PMID: 35472272 DOI: 10.1021/acs.jctc.1c01292] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The microscopic description of the local structure of water remains an open challenge. Here, we adopt an agnostic approach to understanding water's hydrogen bond network using data harvested from molecular dynamics simulations of an empirical water model. A battery of state-of-the-art unsupervised data-science techniques are used to characterize the free-energy landscape of water starting from encoding the water environment using local atomic descriptors, through dimensionality reduction and finally the use of advanced clustering techniques. Analysis of the free energy under ambient conditions was found to be consistent with a rough single basin and independent of the choice of the water model. We find that the fluctuations of the water network occur in a high-dimensional space, which we characterize using a combination of both atomic descriptors and chemical-intuition-based coordinates. We demonstrate that a combination of both types of variables is needed in order to adequately capture the complexity of the fluctuations in the hydrogen bond network at different length scales both at room temperature and also close to the critical point of water. Our results provide a general framework for examining fluctuations in water under different conditions.
Collapse
Affiliation(s)
- Adu Offei-Danso
- The Abdus Salam International Center for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy.,SISSA─International School for Advanced Studies, Via Bonomea 265, 34136 Trieste, Italy
| | - Ali Hassanali
- The Abdus Salam International Center for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| | - Alex Rodriguez
- The Abdus Salam International Center for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
15
|
Malik S, Debnath A. Structural Changes of Interfacial Water upon Fluid-Ripple-Gel Phase Transitions of Bilayers. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139613] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
16
|
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
|
17
|
Kringle L, Thornley WA, Kay BD, Kimmel GA. Isotope effects on the structural transformation and relaxation of deeply supercooled water. J Chem Phys 2022; 156:084501. [DOI: 10.1063/5.0078796] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
We have examined the structure of supercooled liquid D2O as a function of temperature between 185 and 255 K using pulsed laser heating to rapidly heat and cool the sample on a nanosecond timescale. The liquid structure can be represented as a linear combination of two structural motifs, with a transition between them described by a logistic function centered at 218 K with a width of 10 K. The relaxation to a metastable state, which occurred prior to crystallization, exhibited nonexponential kinetics with a rate that was dependent on the initial structural configuration. When the temperature is scaled by the temperature of maximum density, which is an isostructural point of the isotopologues, the structural transition and the non-equilibrium relaxation kinetics of D2O agree remarkably well with those for H2O.
Collapse
Affiliation(s)
- Loni Kringle
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Wyatt A. Thornley
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Bruce D. Kay
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Greg A. Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| |
Collapse
|
18
|
Takahashi O, Yamamura R, Tokushima T, Harada Y. Interpretation of the X-Ray Emission Spectra of Liquid Water through Temperature and Isotope Dependence. PHYSICAL REVIEW LETTERS 2022; 128:086002. [PMID: 35275678 DOI: 10.1103/physrevlett.128.086002] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Accepted: 01/31/2022] [Indexed: 06/14/2023]
Abstract
The interpretation of x-ray emission spectroscopy (XES) spectra in terms of their sensitivity to the hydrogen bonding and the consequent microheterogeneity in liquid water has been debated over a decade. To shed a light on this problem, we report the theoretical reproduction of the debated 1b_{1} peaks observed in the XES spectra of liquid water using semiclassical Kramers-Heisenberg formula. The essence of the temperature and isotope dependence of the 1b_{1} double peaks is explained by molecular dynamics simulations including full vibrational (O─H stretching, bending, and) modes, rotational combined with the density functional theory and core-hole induced dynamics. Some inconsistencies exist with the experimental XES profile, which illustrates the need to employ a more precise theoretical calculations for both geometry sampling and electronic structure using a more sophisticated procedure.
Collapse
Affiliation(s)
- Osamu Takahashi
- Basic Chemistry Program, Graduate School of Advanced Science and Engineering, Hiroshima University, Higashi-Hiroshima 739-8526, Japan
| | - Ryosuke Yamamura
- Department of Chemistry, Hiroshima University, Higashi-Hiroshima 739-8526, Japan
| | - Takashi Tokushima
- MAX IV Laboratory, Lund University, Fotongatan 2, 224 84 Lund, Sweden
| | - Yoshihisa Harada
- Institute for Solid State Physics (ISSP), University of Tokyo, Kashiwanoha, Kashiwa, Chiba 277-8581, Japan
- Synchrotron Radiation Research Organization, University of Tokyo, Sayo-cho, Sayo, Hyogo 679-5198, Japan
| |
Collapse
|
19
|
de Almeida Ribeiro I, de Koning M, Molinero V. Is It Possible to Follow the Structural Evolution of Water in "No-Man's Land" Using a Pulsed-Heating Procedure? J Phys Chem Lett 2022; 13:1085-1089. [PMID: 35080178 DOI: 10.1021/acs.jpclett.1c04106] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The anomalous increase in compressibility and heat capacity of supercooled water has been attributed to its structural transformation of into a four-coordinated liquid. Experiments revealed that κT and Cp peak at TWthermo ≈ 229 K [Kim et al. Science 2017, 358, 1589; Pathak et al. Proc. Natl. Acad. Sci. 2021, 118, e2018379118]. Recently, a pulsed heating procedure (PHP) was employed to interrogate the structure of water, reporting a steep increase in tetrahedrality around TWPHP = 210 ± 3 K [Kringle et al. Science 2020, 369, 1490]. This discrepancy questions whether water structure and thermodynamics are decoupled, or if the shift in TW is an artifact of PHP. Here we implement PHP in molecular simulations. We find that the stationary states captured at the bottom of the pulse are not representative of the thermalized liquid or its inherent structure. Our analysis reveals a temperature-dependent distortion that shifts TWPHP to ∼20 K below TWthermo. We conclude that 2 orders of magnitude faster rates are required to sample water's inherent structure with PHP.
Collapse
Affiliation(s)
- Ingrid de Almeida Ribeiro
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Maurice de Koning
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861, Campinas, São Paulo, Brazil
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
20
|
Esmaeildoost N, Pathak H, Späh A, Lane TJ, Kim KH, Yang C, Amann-Winkel K, Ladd-Parada M, Perakis F, Koliyadu J, Oggenfuss AR, Johnson PJM, Deng Y, Zerdane S, Mankowsky R, Beaud P, Lemke HT, Nilsson A, Sellberg JA. Anomalous temperature dependence of the experimental x-ray structure factor of supercooled water. J Chem Phys 2021; 155:214501. [PMID: 34879659 DOI: 10.1063/5.0075499] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The structural changes of water upon deep supercooling were studied through wide-angle x-ray scattering at SwissFEL. The experimental setup had a momentum transfer range of 4.5 Å-1, which covered the principal doublet of the x-ray structure factor of water. The oxygen-oxygen structure factor was obtained for temperatures down to 228.5 ± 0.6 K. Similar to previous studies, the second diffraction peak increased strongly in amplitude as the structural change accelerated toward a local tetrahedral structure upon deep supercooling. We also observed an anomalous trend for the second peak position of the oxygen-oxygen structure factor (q2). We found that q2 exhibits an unprecedented positive partial derivative with respect to temperature for temperatures below 236 K. Based on Fourier inversion of our experimental data combined with reference data, we propose that the anomalous q2 shift originates from that a repeat spacing in the tetrahedral network, associated with all peaks in the oxygen-oxygen pair-correlation function, gives rise to a less dense local ordering that resembles that of low-density amorphous ice. The findings are consistent with that liquid water consists of a pentamer-based hydrogen-bonded network with low density upon deep supercooling.
Collapse
Affiliation(s)
- Niloofar Esmaeildoost
- Biomedical and X-ray Physics, Department of Applied Physics, AlbaNova University Center, KTH Royal Institute of Technology, S-106 91 Stockholm, Sweden
| | - Harshad Pathak
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Alexander Späh
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Thomas J Lane
- Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Kyung Hwan Kim
- Department of Chemistry, Pohang University of Science and Technology, Pohang 37673, Republic of Korea
| | - Cheolhee Yang
- Department of Chemistry, Pohang University of Science and Technology, Pohang 37673, Republic of Korea
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Marjorie Ladd-Parada
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Fivos Perakis
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | | | | | | | - Yunpei Deng
- SwissFEL, Paul Scherrer Institute, CH-5232 Villigen, Switzerland
| | - Serhane Zerdane
- SwissFEL, Paul Scherrer Institute, CH-5232 Villigen, Switzerland
| | - Roman Mankowsky
- SwissFEL, Paul Scherrer Institute, CH-5232 Villigen, Switzerland
| | - Paul Beaud
- SwissFEL, Paul Scherrer Institute, CH-5232 Villigen, Switzerland
| | - Henrik T Lemke
- SwissFEL, Paul Scherrer Institute, CH-5232 Villigen, Switzerland
| | - Anders Nilsson
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Jonas A Sellberg
- Biomedical and X-ray Physics, Department of Applied Physics, AlbaNova University Center, KTH Royal Institute of Technology, S-106 91 Stockholm, Sweden
| |
Collapse
|
21
|
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
|
22
|
Tetrahedral structure of supercooled water at ambient pressure and its influence on dynamic relaxation: Comparative study of water models. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117269] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|
23
|
Gao Y, Fang H, Ni K. A hierarchical clustering method of hydrogen bond networks in liquid water undergoing shear flow. Sci Rep 2021; 11:9542. [PMID: 33953246 PMCID: PMC8100111 DOI: 10.1038/s41598-021-88810-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2020] [Accepted: 04/14/2021] [Indexed: 02/02/2023] Open
Abstract
Many properties of water, such as turbulent flow, are closely related to water clusters, whereas how water clusters form and transform in bulk water remains unclear. A hierarchical clustering method is introduced to search out water clusters in hydrogen bonded network based on modified Louvain algorithm of graph community. Hydrogen bonds, rings and fragments are considered as 1st-, 2nd-, and 3rd-level structures, respectively. The distribution, dynamics and structural characteristics of 4th- and 5th-level clusters undergoing non-shear- and shear-driven flow are also analyzed at various temperatures. At low temperatures, nearly 50% of water molecules are included in clusters. Over 60% of clusters remain unchanged between neighboring configurations. Obvious collective translational motion of clusters is observed. The topological difference for clusters is elucidated between the inner layer, which favors 6-membered rings, and the external surface layer, which contains more 5-membered rings. Temperature and shearing can not only accelerate the transformation or destruction of clusters at all levels but also change cluster structures. The assembly of large clusters can be used to discretize continuous liquid water to elucidate the properties of liquid water.
Collapse
Affiliation(s)
- Yitian Gao
- grid.12527.330000 0001 0662 3178State Key Laboratory of Hydro-Science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing, 100084 China
| | - Hongwei Fang
- grid.12527.330000 0001 0662 3178State Key Laboratory of Hydro-Science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing, 100084 China
| | - Ke Ni
- grid.12527.330000 0001 0662 3178State Key Laboratory of Hydro-Science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing, 100084 China
| |
Collapse
|
24
|
|
25
|
Moid M, Sastry S, Dasgupta C, Pascal TA, Maiti PK. Dimensionality dependence of the Kauzmann temperature: A case study using bulk and confined water. J Chem Phys 2021; 154:164510. [PMID: 33940812 DOI: 10.1063/5.0047656] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Kauzmann temperature (TK) of a supercooled liquid is defined as the temperature at which the liquid entropy becomes equal to that of the crystal. The excess entropy, the difference between liquid and crystal entropies, is routinely used as a measure of the configurational entropy, whose vanishing signals the thermodynamic glass transition. The existence of the thermodynamic glass transition is a widely studied subject, and of particular recent interest is the role of dimensionality in determining the presence of a glass transition at a finite temperature. The glass transition in water has been investigated intensely and is challenging as the experimental glass transition appears to occur at a temperature where the metastable liquid is strongly prone to crystallization and is not stable. To understand the dimensionality dependence of the Kauzmann temperature in water, we study computationally bulk water (three-dimensions), water confined in the slit pore of the graphene sheet (two-dimensions), and water confined in the pore of the carbon nanotube of chirality (11,11) having a diameter of 14.9 Å (one-dimension), which is the lowest diameter where amorphous water does not always crystallize into nanotube ice in the supercooled region. Using molecular dynamics simulations, we compute the entropy of water in bulk and under reduced dimensional nanoscale confinement to investigate the variation of the Kauzmann temperature with dimension. We obtain a value of TK (133 K) for bulk water in good agreement with experiments [136 K (C. A. Angell, Science 319, 582-587 (2008) and K. Amann-Winkel et al., Proc. Natl. Acad. Sci. U. S. A. 110, 17720-17725 (2013)]. However, for confined water, in two-dimensions and one-dimension, we find that there is no finite temperature Kauzmann point (in other words, the Kauzmann temperature is 0 K). Analysis of the fluidicity factor, a measure of anharmonicity in the oscillation of normal modes, reveals that the Kauzmann temperature can also be computed from the difference in the fluidicity factor between amorphous and ice phases.
Collapse
Affiliation(s)
- Mohd Moid
- Department of Physics, Centre for Condensed Matter Theory, Indian Institute of Science, Bangalore 560012, India
| | - Srikanth Sastry
- Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore 560064, India
| | - Chandan Dasgupta
- Department of Physics, Centre for Condensed Matter Theory, Indian Institute of Science, Bangalore 560012, India
| | - Tod A Pascal
- Department of Nanoengineering and Chemical Engineering, University of California San Diego, La Jolla, California 92023, USA
| | - Prabal K Maiti
- Department of Physics, Centre for Condensed Matter Theory, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
26
|
Structural relaxation and crystallization in supercooled water from 170 to 260 K. Proc Natl Acad Sci U S A 2021; 118:2022884118. [PMID: 33790015 DOI: 10.1073/pnas.2022884118] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The origin of water's anomalous properties has been debated for decades. Resolution of the problem is hindered by a lack of experimental data in a crucial region of temperatures, T, and pressures where supercooled water rapidly crystallizes-a region often referred to as "no man's land." A recently developed technique where water is heated and cooled at rates greater than 109 K/s now enables experiments in this region. Here, it is used to investigate the structural relaxation and crystallization of deeply supercooled water for 170 K < T < 260 K. Water's relaxation toward a new equilibrium structure depends on its initial structure with hyperquenched glassy water (HQW) typically relaxing more quickly than low-density amorphous solid water (LDA). For HQW and T > 230 K, simple exponential relaxation kinetics is observed. For HQW at lower temperatures, increasingly nonexponential relaxation is observed, which is consistent with the dynamics expected on a rough potential energy landscape. For LDA, approximately exponential relaxation is observed for T > 230 K and T < 200 K, with nonexponential relaxation only at intermediate temperatures. At all temperatures, water's structure can be reproduced by a linear combination of two, local structural motifs, and we show that a simple model accounts for the complex kinetics within this context. The relaxation time, τ rel , is always shorter than the crystallization time, τ xtal For HQW, the ratio, τ xtal /τ rel , goes through a minimum at ∼198 K where the ratio is about 60.
Collapse
|
27
|
Verde AR, de Oca JMM, Accordino SR, Alarcón LM, Appignanesi GA. Structural aspects of an energy-based water classification index and the structure-dynamics link in glassy relaxation. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:47. [PMID: 33783648 DOI: 10.1140/epje/s10189-021-00057-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Accepted: 03/15/2021] [Indexed: 06/12/2023]
Abstract
An energy-based structural indicator for water, [Formula: see text], has been recently introduced by our group. In turn, in this work we aim at: (1) demonstrating that [Formula: see text] is indeed able to correctly classify water molecules between locally structured tetrahedral (T) and locally distorted (D) ones, circumventing the usual problem of certain previous indicators of overestimating the distorted state; (2) correlating [Formula: see text] with dynamic propensity, a measure of the molecular mobility tendency, in order to seek for the existence of a connection between structure and dynamics within the supercooled regime. More specifically, in the first part of this work we will show that [Formula: see text] accurately discriminates between merely thermally deformed local molecular arrangements and truly distorted molecules (defects). This fact will be made evident not only from radial distribution function results but also from the dynamic propensity distributions of the different kinds of molecules. In turn, we shall devote the second part of this work to finding correlations between T and D molecules with low- and high-dynamic-propensity molecules, respectively, thus revealing the existence of a link between local structure and dynamics, while also making evident the dominant role of the D molecules (defects) in the structural relaxation. Moreover, the availability of a proper molecular classification technique will enable us to study the timescale of such influence of structure on dynamics by defining a modified dynamic propensity measure and by applying it to the structured and unstructured water molecular states.
Collapse
Affiliation(s)
- Alejandro R Verde
- Departamento de Química, INQUISUR, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
| | - Joan Manuel Montes de Oca
- Departamento de Química, INQUISUR, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, 60637, USA
| | - Sebastián R Accordino
- Departamento de Química, INQUISUR, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
| | - Laureano M Alarcón
- Departamento de Química, INQUISUR, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- Departamento de Química, INQUISUR, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina.
| |
Collapse
|
28
|
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
|
29
|
Horstmann R, Vogel M. Relations between thermodynamics, structures, and dynamics for modified water models in their supercooled regimes. J Chem Phys 2021; 154:054502. [DOI: 10.1063/5.0037080] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- R. Horstmann
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - M. Vogel
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| |
Collapse
|
30
|
Jin J, Pak AJ, Han Y, Voth GA. A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). II. Temperature transferability and structural properties at low temperature. J Chem Phys 2021; 154:044105. [PMID: 33514078 PMCID: PMC7826166 DOI: 10.1063/5.0026652] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Accepted: 12/14/2020] [Indexed: 11/14/2022] Open
Abstract
A number of studies have constructed coarse-grained (CG) models of water to understand its anomalous properties. Most of these properties emerge at low temperatures, and an accurate CG model needs to be applicable to these low-temperature ranges. However, direct use of CG models parameterized from other temperatures, e.g., room temperature, encounters a problem known as transferability, as the CG potential essentially follows the form of the many-body CG free energy function. Therefore, temperature-dependent changes to CG interactions must be accounted for. The collective behavior of water at low temperature is generally a many-body process, which often motivates the use of expensive many-body terms in the CG interactions. To surmount the aforementioned problems, we apply the Bottom-Up Many-Body Projected Water (BUMPer) CG model constructed from Paper I to study the low-temperature behavior of water. We report for the first time that the embedded three-body interaction enables BUMPer, despite its pairwise form, to capture the growth of ice at the ice/water interface with corroborating many-body correlations during the crystal growth. Furthermore, we propose temperature transferable BUMPer models that are indirectly constructed from the free energy decomposition scheme. Changes in CG interactions and corresponding structures are faithfully recapitulated by this framework. We further extend BUMPer to examine its ability to predict the structure, density, and diffusion anomalies by employing an alternative analysis based on structural correlations and pairwise potential forms to predict such anomalies. The presented analysis highlights the existence of these anomalies in the low-temperature regime and overcomes potential transferability problems.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Alexander J. Pak
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Yining Han
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
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
|
Kringle L, Thornley WA, Kay BD, Kimmel GA. Reversible structural transformations in supercooled liquid water from 135 to 245 K. Science 2020; 369:1490-1492. [DOI: 10.1126/science.abb7542] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Accepted: 07/28/2020] [Indexed: 11/02/2022]
Abstract
A fundamental understanding of the unusual properties of water remains elusive because of the limited data at the temperatures and pressures needed to decide among competing theories. We investigated the structural transformations of transiently heated supercooled water films, which evolved for several nanoseconds per pulse during fast laser heating before quenching to 70 kelvin (K). Water’s structure relaxed from its initial configuration to a steady-state configuration before appreciable crystallization. Over the full temperature range investigated, all structural changes were reversible and reproducible by a linear combination of high- and low-temperature structural motifs. The fraction of the liquid with the high-temperature motif decreased rapidly as the temperature decreased from 245 to 190 K, consistent with the predictions of two-state “mixture” models for supercooled water in the supercritical regime.
Collapse
Affiliation(s)
- Loni Kringle
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, WA 99352, USA
| | - Wyatt A. Thornley
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, WA 99352, USA
| | - Bruce D. Kay
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, WA 99352, USA
| | - Greg A. Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, WA 99352, USA
| |
Collapse
|
34
|
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
|
35
|
Kim KH, Späh A, Pathak H, Yang C, Bonetti S, Amann-Winkel K, Mariedahl D, Schlesinger D, Sellberg JA, Mendez D, van der Schot G, Hwang HY, Clark J, Shigeki O, Tadashi T, Harada Y, Ogasawara H, Katayama T, Nilsson A, Perakis F. Anisotropic X-Ray Scattering of Transiently Oriented Water. PHYSICAL REVIEW LETTERS 2020; 125:076002. [PMID: 32857536 DOI: 10.1103/physrevlett.125.076002] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Accepted: 07/08/2020] [Indexed: 06/11/2023]
Abstract
We study the structural dynamics of liquid water by time-resolved anisotropic x-ray scattering under the optical Kerr effect condition. In this way, we can separate the anisotropic scattering decay of 160 fs from the delayed temperature increase of ∼0.1 K occurring at 1 ps and quantify transient changes in the O-O pair distribution function. Polarizable molecular dynamics simulations reproduce well the experiment, indicating transient alignment of molecules along the electric field, which shortens the nearest-neighbor distances. In addition, analysis of the simulated water local structure provides evidence that two hypothesized fluctuating water configurations exhibit different polarizability.
Collapse
Affiliation(s)
- Kyung Hwan Kim
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
- Department of Chemistry, POSTECH, Pohang 37673, Republic of Korea
| | - Alexander Späh
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Harshad Pathak
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Cheolhee Yang
- Department of Chemistry, POSTECH, Pohang 37673, Republic of Korea
| | - Stefano Bonetti
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
- Department of Molecular Sciences and Nanosystems, Ca' Foscari University of Venice, 30172 Venice-Mestre, Italy
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Daniel Mariedahl
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Daniel Schlesinger
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
- Department of Environmental Science and Bolin Centre for Climate Research, Stockholm University, 114 18 Stockholm, Sweden
| | - Jonas A Sellberg
- Biomedical and X-Ray Physics, Department of Applied Physics, AlbaNova University Center, KTH Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Derek Mendez
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Gijs van der Schot
- Department of Cell and Molecular Biology, Laboratory of Molecular Biophysics, Uppsala University, SE-75124 Uppsala, Sweden
| | - Harold Y Hwang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Jesse Clark
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Owada Shigeki
- Japan Synchrotron Radiation Research Institute, Kouto 1-1-1, Sayo, Hyogo 679-5198, Japan
| | - Togashi Tadashi
- Japan Synchrotron Radiation Research Institute, Kouto 1-1-1, Sayo, Hyogo 679-5198, Japan
| | - Yoshihisa Harada
- Institute for Solid State Physics, The University of Tokyo, Kouto 1-1-1, Sayo, Hyogo 679-5148, Japan
| | | | - Tetsuo Katayama
- Japan Synchrotron Radiation Research Institute, Kouto 1-1-1, Sayo, Hyogo 679-5198, Japan
| | - Anders Nilsson
- 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
|
36
|
Tang PH, Wu TM. Instantaneous normal mode analysis for OKE reduced spectra of liquid and supercooled water: Contributions of low-density and high-density liquids. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2019.112363] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
37
|
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
|
38
|
Bachler J, Handle PH, Giovambattista N, Loerting T. Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations. Phys Chem Chem Phys 2019; 21:23238-23268. [PMID: 31556899 DOI: 10.1039/c9cp02953b] [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/21/2022]
Abstract
One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.
Collapse
Affiliation(s)
- Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria.
| | | | | | | |
Collapse
|
39
|
Pathak H, Späh A, Amann-Winkel K, Perakis F, Kim KKH, Nilsson A. Temperature dependent anomalous fluctuations in water: shift of ≈1 kbar between experiment and classical force field simulations. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1649486] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Harshad Pathak
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Alexander Späh
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Fivos Perakis
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Kyung Kyung Hwan Kim
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| | - Anders Nilsson
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| |
Collapse
|
40
|
Camisasca G, Pathak H, Wikfeldt KT, Pettersson LGM. Radial distribution functions of water: Models vs experiments. J Chem Phys 2019; 151:044502. [DOI: 10.1063/1.5100871] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Affiliation(s)
- Gaia Camisasca
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10609 Stockholm, Sweden
| | - Harshad Pathak
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10609 Stockholm, Sweden
| | - Kjartan Thor Wikfeldt
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10609 Stockholm, Sweden
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10609 Stockholm, Sweden
| |
Collapse
|
41
|
Montes de Oca JM, Accordino SR, Verde AR, Alarcón LM, Appignanesi GA. Structural features of high-local-density water molecules: Insights from structure indicators based on the translational order between the first two molecular shells. Phys Rev E 2019; 99:062601. [PMID: 31330696 DOI: 10.1103/physreve.99.062601] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2019] [Indexed: 11/07/2022]
Abstract
The two-liquids scenario for liquid water assumes the existence of two competing preferential local molecular structural states characterized by either low or high local density. While the former is expected to present good local order thus involving privileged structures, the latter is usually regarded as conforming a high-entropy unstructured state. A main difference in the local arrangement of such "classes" of water molecules can be inferred from the degree of translational order between the first and second molecular shells. This is so, since the low-local-density molecules present a clear gap between the first two shells while in the case of the high-local-density ones, one or more molecules from the second shell have collapsed toward the first one, thus populating the intershell region. Some structural indicators, like the widely employed local structure index and the recently introduced ζ index, have been devised precisely on the basis of this observation, being successful in detecting well-structured low-local-density molecules. However, the nature of the high-local-density state has been mainly disregarded over the years. In this work we employ molecular dynamics simulations for two water models (the extended simple point charge model and the five-site model) at the liquid and supercooled regimes combined with the inherent dynamics approach (energy minimizations of the instantaneous configurations) in order to both rationalize the detailed structural and topological information that these indicators provide and to advance in our understanding of the high-density state.
Collapse
Affiliation(s)
- Joan Manuel Montes de Oca
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Sebastián R Accordino
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Alejandro R Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Laureano M Alarcón
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
42
|
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
|
43
|
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
|
44
|
Handle PH, Sciortino F, Giovambattista N. Glass polymorphism in TIP4P/2005 water: A description based on the potential energy landscape formalism. J Chem Phys 2019; 150:244506. [PMID: 31255050 DOI: 10.1063/1.5100346] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The potential energy landscape (PEL) formalism is a statistical mechanical approach to describe supercooled liquids and glasses. Here, we use the PEL formalism to study the pressure-induced transformations between low-density amorphous ice (LDA) and high-density amorphous ice (HDA) using computer simulations of the TIP4P/2005 molecular model of water. We find that the properties of the PEL sampled by the system during the LDA-HDA transformation exhibit anomalous behavior. In particular, at conditions where the change in density during the LDA-HDA transformation is approximately discontinuous, reminiscent of a first-order phase transition, we find that (i) the inherent structure (IS) energy, eIS(V), is a concave function of the volume and (ii) the IS pressure, PIS(V), exhibits a van der Waals-like loop. In addition, the curvature of the PEL at the IS is anomalous, a nonmonotonic function of V. In agreement with previous studies, our work suggests that conditions (i) and (ii) are necessary (but not sufficient) signatures of the PEL for the LDA-HDA transformation to be reminiscent of a first-order phase transition. We also find that one can identify two different regions of the PEL, one associated with LDA and another with HDA. Our computer simulations are performed using a wide range of compression/decompression and cooling rates. In particular, our slowest cooling rate (0.01 K/ns) is within the experimental rates employed in hyperquenching experiments to produce LDA. Interestingly, the LDA-HDA transformation pressure that we obtain at T = 80 K and at different rates extrapolates remarkably well to the corresponding experimental pressure.
Collapse
Affiliation(s)
- Philip H Handle
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020 Innsbruck, Austria
| | - Francesco Sciortino
- Department of Physics, Sapienza-University of Rome, Piazzale Aldo Moro 5, I-00185 Roma, Italy
| | - Nicolas Giovambattista
- Department of Physics, Brooklyn College of the City University of New York, New York, New York 10016, USA
| |
Collapse
|
45
|
Verde AR, Montes de Oca JM, Accordino SR, Alarcón LM, Appignanesi GA. Comparing the performance of two structural indicators for different water models while seeking for connections between structure and dynamics in the glassy regime. J Chem Phys 2019; 150:244504. [DOI: 10.1063/1.5108796] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Alejandro R. Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Joan Manuel Montes de Oca
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Sebastián R. Accordino
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Laureano M. Alarcón
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A. Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
46
|
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
|
47
|
Camisasca G, Galamba N, Wikfeldt KT, Pettersson LGM. Translational and rotational dynamics of high and low density TIP4P/2005 water. J Chem Phys 2019; 150:224507. [PMID: 31202216 DOI: 10.1063/1.5079956] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We use molecular dynamics simulations using TIP4P/2005 to investigate the self- and distinct-van Hove functions for different local environments of water, classified using the local structure index as an order parameter. The orientational dynamics were studied through the calculation of the time-correlation functions of different-order Legendre polynomials in the OH-bond unit vector. We found that the translational and orientational dynamics are slower for molecules in a low-density local environment and correspondingly the mobility is enhanced upon increasing the local density, consistent with some previous works, but opposite to a recent study on the van Hove function. From the analysis of the distinct dynamics, we find that the second and fourth peaks of the radial distribution function, previously identified as low density-like arrangements, show long persistence in time. The analysis of the time-dependent interparticle distance between the central molecule and the first coordination shell shows that particle identity persists longer than distinct van Hove correlations. The motion of two first-nearest-neighbor molecules thus remains coupled even when this correlation function has been completely decayed. With respect to the orientational dynamics, we show that correlation functions of molecules in a low-density environment decay exponentially, while molecules in a local high-density environment exhibit bi-exponential decay, indicating that dynamic heterogeneity of water is associated with the heterogeneity among high-density and between high-density and low-density species. This bi-exponential behavior is associated with the existence of interstitial waters and the collapse of the second coordination sphere in high-density arrangements, but not with H-bond strength.
Collapse
Affiliation(s)
- Gaia Camisasca
- Department of Physics, Stockholm University, 106 91 Stockholm, Sweden
| | - Nuno Galamba
- Centre of Chemistry and Biochemistry and Biosystems and Integrative Sciences Institute, Faculty of Sciences of the University of Lisbon, C8, Campo Grande, 1749-016 Lisbon, Portugal
| | | | | |
Collapse
|
48
|
Pathak H, Späh A, Kim KH, Tsironi I, Mariedahl D, Blanco M, Huotari S, Honkimäki V, Nilsson A. Intermediate range O-O correlations in supercooled water down to 235 K. J Chem Phys 2019; 150:224506. [PMID: 31202250 DOI: 10.1063/1.5100811] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Wide angle x-ray scattering of supercooled water down to 234.8 K was studied using high energy x rays at the European Synchrotron Radiation Facility. The oxygen-oxygen pair distribution function (PDF) was calculated from the scattering pattern out to the 5th peak at an intermolecular distance, r ≈ 11 Å. We observe that the 4th peak and the 5th peak in the PDF increase in height upon supercooling. We also observe that the 4th peak position (r4) shifts to shorter distances upon supercooling consistent with previous studies, but we see a more rapid change at the lowest temperature. The running oxygen-oxygen coordination number is calculated for 5 different temperatures, and an isosbestic point at riso = 3.31 ± 0.05 Å was found corresponding to a coordination number of 4.39 ± 0.15. The comparison of the PDF of the coldest water with that of amorphous ice shows distinct differences. We propose that there are 5-member pentamer rings in low density liquid-like structures giving rise to the sharp correlations at r ≈ 9 Å and r ≈ 11 Å.
Collapse
Affiliation(s)
- Harshad Pathak
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Alexander Späh
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Kyung Hwan Kim
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Ifigeneia Tsironi
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Daniel Mariedahl
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| | - Maria Blanco
- European Synchrotron Radiation Facility, 71 Avenue des Martyrs, 38000 Grenoble, France
| | - Simo Huotari
- Department of Physics, University of Helsinki, FI-00014 Helsinki, Finland
| | - Veijo Honkimäki
- European Synchrotron Radiation Facility, 71 Avenue des Martyrs, 38000 Grenoble, France
| | - Anders Nilsson
- Department of Physics, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
49
|
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
|
50
|
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]
|