1
|
Förster M, Ukoji N, Sahle CJ, Niskanen J, Sakrowski R, Surmeier G, Weis C, Irifune T, Imoto S, Yavas H, Huotari S, Marx D, Sternemann C, Tse JS. Generating interstitial water within the persisting tetrahedral H-bond network explains density increase upon compressing liquid water. Proc Natl Acad Sci U S A 2024; 121:e2403662121. [PMID: 39284048 PMCID: PMC11441526 DOI: 10.1073/pnas.2403662121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2024] [Accepted: 07/16/2024] [Indexed: 10/02/2024] Open
Abstract
Despite its ubiquitous nature, the atomic structure of water in its liquid state is still controversially debated. We use a combination of X-ray Raman scattering spectroscopy in conjunction with ab initio and path integral molecular dynamics simulations to study the local atomic and electronic structure of water under high pressure conditions. Systematically increasing fingerprints of non-hydrogen-bonded H[Formula: see text]O molecules in the first hydration shell are identified in the experimental and computational oxygen K-edge excitation spectra. This provides evidence for a compaction mechanism in terms of a continuous collapse of the second hydration shell with increasing pressure via generation of interstitial water within locally tetrahedral hydrogen-bonding environments.
Collapse
Affiliation(s)
- Mirko Förster
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Nnanna Ukoji
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SKS7N 5E2, Canada
| | | | - Johannes Niskanen
- Department of Physics and Astronomy, University of Turku, Turun yliopistoFI-20014, Finland
| | - Robin Sakrowski
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Göran Surmeier
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Christopher Weis
- Fakultät Physik/DELTA, Technische Universität Dortmund, Dortmund44227, Germany
| | - Tetsuo Irifune
- Geodynamics Research Center, Ehime University, Matsuyama790, Japan
| | - Sho Imoto
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum44780, Germany
| | - Hasan Yavas
- Deutsches Elektronen-Synchrotron, Hamburg22607, Germany
- Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, CA94025
| | - Simo Huotari
- Department of Physics, University of Helsinki, HelsinkiFI-00014, Finland
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum44780, Germany
| | | | - John S. Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SKS7N 5E2, Canada
| |
Collapse
|
2
|
Gao Y, Wu J, Feng Y, Han J, Fang H. Effects of Hydrogen Bond Networks on Viscosity in Aqueous Solutions. J Phys Chem B 2024; 128:8984-8996. [PMID: 39236306 DOI: 10.1021/acs.jpcb.4c03856] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
In aqueous solutions, the impact of ions on hydrogen bond networks plays a crucial role in transport properties. We used molecular dynamics simulations to explain how ions affect viscosity through structural changes. We developed a quantitative model to describe the effect of ions on viscosity. The model comprises two parts: the addition of ions alters hydrogen bond networks, and changes in hydrogen bond networks exponentially lead to changes in viscosity. The influence of ions on hydrogen bond networks involves the following mechanisms: first, ions can disrupt the tetrahedral structures within the first solvation shell into three-coordinated structures through substitution; second, structural changes within the first shells affect the global hydrogen bond network through electrostatic forces and the hindrance of ionic volumes. By analyzing the mechanisms of how hydrogen bond networks determine viscosity through the decomposition of viscosity, we found that the proportion of potential viscosity in aqueous solutions primarily increases due to the enhancement of non-hydrogen bonding interactions, and the proportion of hydrogen bonding viscosity decreases accordingly. Our results demonstrate that hydrogen bond networks are crucial for describing the changes in transport phenomena affected by external factors.
Collapse
Affiliation(s)
- Yitian Gao
- State Key Laboratory of Hydro-Science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
- China Renewable Energy Engineering Institute, Beijing 100120, China
| | - Jian Wu
- Department of Physics, Tsinghua University, Beijing 100084, China
| | - Yixuan Feng
- State Key Laboratory of Hydro-Science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Jiale Han
- State Key Laboratory of Hydro-Science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Hongwei Fang
- State Key Laboratory of Hydro-Science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
3
|
Malosso C, Manko N, Izzo MG, Baroni S, Hassanali A. Evidence of ferroelectric features in low-density supercooled water from ab initio deep neural-network simulations. Proc Natl Acad Sci U S A 2024; 121:e2407295121. [PMID: 39083416 PMCID: PMC11317578 DOI: 10.1073/pnas.2407295121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2024] [Accepted: 06/30/2024] [Indexed: 08/02/2024] Open
Abstract
Over the last decade, an increasing body of evidence has emerged, supporting the existence of a metastable liquid-liquid critical point in supercooled water whereby two distinct liquid phases of different densities coexist. Analyzing long molecular dynamics simulations performed using deep neural-network force fields trained to accurate quantum mechanical data, we demonstrate that the low-density liquid phase displays a strong propensity toward spontaneous polarization, as witnessed by large and long-lived collective dipole fluctuations. Our findings suggest that the dynamical stability of the low-density phase, and hence the transition from high-density to low-density liquid, is triggered by a collective process involving an accumulation of rotational angular jumps, which could ignite large dipole fluctuations. This dynamical transition involves subtle changes in the electronic polarizability of water molecules which affects their rotational mobility within the two phases. These findings hold the potential for catalyzing activity in the search for dielectric-based probes of the putative second critical point.
Collapse
Affiliation(s)
- Cesare Malosso
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
| | - Natalia Manko
- Condensed Matter and Statistical Physics (CMSP), The Abdus Salam Centre for Theoretical Physics, Trieste34151, Italy
| | - Maria Grazia Izzo
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
| | - Stefano Baroni
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
- Consiglio Nazionale delle Ricerche-Istituto Officina dei Materiali, Scuola Internazionale Superiore di Studi Avanzati Unit, Trieste34136, Italy
| | - Ali Hassanali
- Condensed Matter and Statistical Physics (CMSP), The Abdus Salam Centre for Theoretical Physics, Trieste34151, Italy
| |
Collapse
|
4
|
Di Fonte N, Faccio C, Zanetti-Polzi L, Daidone I. Early prediction of spinodal-like relaxation events in supercooled liquid water. J Chem Phys 2024; 161:034505. [PMID: 39012812 DOI: 10.1063/5.0211031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Accepted: 07/01/2024] [Indexed: 07/18/2024] Open
Abstract
Several computational studies on different water models reported evidence of a phase transition in supercooled conditions between two liquid states of water differing in density: the high-density liquid (HDL) and the low-density liquid (LDL). Yet, conclusive experimental evidence of the existence of a phase transition between the two liquid water phases could not be obtained due to fast crystallization in the region where the phase transition should occur. For the same reason, the investigation of possible transition mechanisms between the two phases is committed to computational investigations. In this work, we simulate an out-of-equilibrium temperature-induced transition from the LDL to the HDL-like state in the TIP4P/2005 water model. To structurally characterize the system relaxation, we use the node total communicability (NTC) we recently proposed as an effective order parameter to discriminate the two liquid phases differing in density. We find that the relaxation process is compatible with a spinodal-like scenario. We observe the formation of HDL-like domains in the LDL phase and we characterize their fluctuating behavior and subsequent coarsening and stabilization. Furthermore, we find that the formation of stable HDL-like domains is favored in the regions where the early formation of small patches of highly connected HDL-like molecules (i.e., with very high NTC values) is observed. Besides characterizing the LDL- to HDL-like relaxation from a structural point of view, these results also show that the NTC order parameter can serve as an early-time predictor of the regions from which the transition process initiates.
Collapse
Affiliation(s)
- Nico Di Fonte
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| | - Chiara Faccio
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
- Department of Mathematics "Tullio Levi-Civita," University of Padova, Via Trieste 63, Padova 35121, Italy
| | - Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, Modena 41125, Italy
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| |
Collapse
|
5
|
Zhang H, Zhou Y, Zhang Z, Sun H, Pan Z, Mou M, Zhang W, Ye Q, Hou T, Li H, Hsieh CY, Zhu F. Large Language Model-Based Natural Language Encoding Could Be All You Need for Drug Biomedical Association Prediction. Anal Chem 2024. [PMID: 39011990 DOI: 10.1021/acs.analchem.4c01793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/17/2024]
Abstract
Analyzing drug-related interactions in the field of biomedicine has been a critical aspect of drug discovery and development. While various artificial intelligence (AI)-based tools have been proposed to analyze drug biomedical associations (DBAs), their feature encoding did not adequately account for crucial biomedical functions and semantic concepts, thereby still hindering their progress. Since the advent of ChatGPT by OpenAI in 2022, large language models (LLMs) have demonstrated rapid growth and significant success across various applications. Herein, LEDAP was introduced, which uniquely leveraged LLM-based biotext feature encoding for predicting drug-disease associations, drug-drug interactions, and drug-side effect associations. Benefiting from the large-scale knowledgebase pre-training, LLMs had great potential in drug development analysis owing to their holistic understanding of natural language and human topics. LEDAP illustrated its notable competitiveness in comparison with other popular DBA analysis tools. Specifically, even in simple conjunction with classical machine learning methods, LLM-based feature representations consistently enabled satisfactory performance across diverse DBA tasks like binary classification, multiclass classification, and regression. Our findings underpinned the considerable potential of LLMs in drug development research, indicating a catalyst for further progress in related fields.
Collapse
Affiliation(s)
- Hanyu Zhang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, Alibaba-Zhejiang University Joint Research Center of Future Digital Healthcare, Hangzhou 330110, China
- College of Pharmaceutical Sciences, The Second Affiliated Hospital, Zhejiang University School of Medicine, State Key Laboratory of Advanced Drug Delivery and Release Systems, Zhejiang University, Hangzhou 310058, China
| | - Yuan Zhou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Zhichao Zhang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Huaicheng Sun
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Ziqi Pan
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Minjie Mou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Wei Zhang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Qing Ye
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Honglin Li
- Innovation Center for AI and Drug Discovery, East China Normal University, Shanghai 200062, China
| | - Chang-Yu Hsieh
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | - Feng Zhu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, Alibaba-Zhejiang University Joint Research Center of Future Digital Healthcare, Hangzhou 330110, China
- College of Pharmaceutical Sciences, The Second Affiliated Hospital, Zhejiang University School of Medicine, State Key Laboratory of Advanced Drug Delivery and Release Systems, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
6
|
Sun Q, Chen YN, Liu YZ. The Effects of External Interfaces on Hydrophobic Interactions I: Smooth Surface. Molecules 2024; 29:3128. [PMID: 38999080 PMCID: PMC11243484 DOI: 10.3390/molecules29133128] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2024] [Revised: 06/25/2024] [Accepted: 06/28/2024] [Indexed: 07/14/2024] Open
Abstract
External interfaces, such as the air-water and solid-liquid interfaces, are ubiquitous in nature. Hydrophobic interactions are considered the fundamental driving force in many physical and chemical processes occurring in aqueous solutions. It is important to understand the effects of external interfaces on hydrophobic interactions. According to the structural studies on liquid water and the air-water interface, the external interface primarily affects the structure of the topmost water layer (interfacial water). Therefore, an external interface may affect hydrophobic interactions. The effects of interfaces on hydrophobicity are related not only to surface molecular polarity but also to the geometric characteristics of the external interface, such as shape and surface roughness. This study is devoted to understanding the effects of a smooth interface on hydrophobicity. Due to hydrophobic interactions, the solutes tend to accumulate at external interfaces to maximize the hydrogen bonding of water. Additionally, these can be demonstrated by the calculated potential mean forces (PMFs) using molecular dynamic (MD) simulations.
Collapse
Affiliation(s)
- Qiang Sun
- Key Laboratory of Orogenic Belts and Crustal Evolution, Ministry of Education, The School of Earth and Space Sciences, Peking University, Beijing 100871, China (Y.-Z.L.)
| | | | | |
Collapse
|
7
|
Fransson T, Pettersson LGM. TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional. J Chem Phys 2024; 160:234105. [PMID: 38884399 DOI: 10.1063/5.0209719] [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: 06/03/2024] [Indexed: 06/18/2024] Open
Abstract
We investigate the performance of time-dependent density functional theory (TDDFT) for reproducing high-level reference x-ray absorption spectra of liquid water and water clusters. For this, we apply the integrated absolute difference (IAD) metric, previously used for x-ray emission spectra of liquid water [T. Fransson and L. G. M. Pettersson, J. Chem. Theory Comput. 19, 7333-7342 (2023)], in order to investigate which exchange-correlation (xc) functionals yield TDDFT spectra in best agreement to reference, as well as to investigate the suitability of IAD for x-ray absorption spectroscopy spectrum calculations. Considering highly asymmetric and symmetric six-molecule clusters, it is seen that long-range corrected xc-functionals are required to yield good agreement with the reference coupled cluster (CC) and algebraic-diagrammatic construction spectra, with 100% asymptotic Hartree-Fock exchange resulting in the lowest IADs. The xc-functionals with best agreement to reference have been adopted for larger water clusters, yielding results in line with recently published CC theory, but which still show some discrepancies in the relative intensity of the features compared to experiment.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| | - Lars G M Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| |
Collapse
|
8
|
Liu YZ, Chen YN, Sun Q. The Dependence of Hydrophobic Interactions on the Shape of Solute Surface. Molecules 2024; 29:2601. [PMID: 38893477 PMCID: PMC11173737 DOI: 10.3390/molecules29112601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2024] [Revised: 05/26/2024] [Accepted: 05/29/2024] [Indexed: 06/21/2024] Open
Abstract
According to our recent studies on hydrophobicity, this work is aimed at understanding the dependence of hydrophobic interactions on the shape of a solute's surface. It has been observed that dissolved solutes primarily affect the structure of interfacial water, which refers to the top layer of water at the interface between the solute and water. As solutes aggregate in a solution, hydrophobic interactions become closely related to the transition of water molecules from the interfacial region to the bulk water. It is inferred that hydrophobic interactions may depend on the shape of the solute surface. To enhance the strength of hydrophobic interactions, the solutes tend to aggregate, thereby minimizing their surface area-to-volume ratio. This also suggests that hydrophobic interactions may exhibit directional characteristics. Moreover, this phenomenon can be supported by calculated potential mean forces (PMFs) using molecular dynamics (MD) simulations, where different surfaces, such as convex, flat, or concave, are associated with a sphere. Furthermore, this concept can be extended to comprehend the molecular packing parameter, commonly utilized in studying the self-assembly behavior of amphiphilic molecules in aqueous solutions.
Collapse
Affiliation(s)
| | | | - Qiang Sun
- Key Laboratory of Orogenic Belts and Crustal Evolution, Ministry of Education, The School of Earth and Space Sciences, Peking University, Beijing 100871, China; (Y.-Z.L.); (Y.-N.C.)
| |
Collapse
|
9
|
Kimmel GA. Isotope effects in supercooled H2O and D2O and a corresponding-states-like rescaling of the temperature and pressure. J Chem Phys 2024; 160:194508. [PMID: 38767262 DOI: 10.1063/5.0207719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2024] [Accepted: 04/30/2024] [Indexed: 05/22/2024] Open
Abstract
Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H2O and D2O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations-with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a "corresponding-states-like" rescaling for water properties in three classical water models that all exhibit a liquid-liquid transition and critical point [Uralcan et al., J. Chem. Phys. 150, 064503 (2019)], the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H2O and D2O is investigated here. Utilizing previously published data and equations-of-state for H2O and D2O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid-liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water's anomalies.
Collapse
Affiliation(s)
- Greg A Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| |
Collapse
|
10
|
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
|
11
|
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
|
12
|
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
|
13
|
Loubet NA, Verde AR, Appignanesi GA. A water structure indicator suitable for generic contexts: Two-liquid behavior at hydration and nanoconfinement conditions and a molecular approach to hydrophobicity and wetting. J Chem Phys 2024; 160:144502. [PMID: 38587223 DOI: 10.1063/5.0203989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2024] [Accepted: 03/25/2024] [Indexed: 04/09/2024] Open
Abstract
In a recent work, we have briefly introduced a new structural index for water that, unlike previous indicators, was devised specifically for generic contexts beyond bulk conditions, making it suitable for hydration and nanoconfinement settings. In this work, we shall study this metric in detail, demonstrating its ability to reveal the existence of a fine-tuned interplay between the local structure and energetics in liquid water. This molecular principle enables the establishment of an extended hydrogen bond network, while simultaneously allowing for the existence of network defects by compensating for uncoordinated sites. By studying different water models and different temperatures encompassing both the normal liquid and the supercooled regime, this molecular mechanism will be shown to underlie the two-state behavior of bulk water. In addition, by studying functionalized self-assembled monolayers and diverse graphene-like surfaces, we shall show that this principle is also operative at hydration and nanoconfinement conditions, thus generalizing the validity of the two-liquid scenario of water to these contexts. This approach will allow us to define conditions for wettability, providing an accurate measure of hydrophobicity and a reliable predictor of filling and drying transitions. Hence, it might open the possibility of elucidating the active role of water in the broad fields of biophysics and materials science. As a preliminary step, we shall study the hydration structure and hydrophilicity of graphene-like systems (parallel graphene sheets and carbon nanotubes) as a function of the confinement dimensionality.
Collapse
Affiliation(s)
- Nicolás A Loubet
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Alejandro R Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
14
|
Gao Y, Wu J, Feng Y, Han J, Fang H. Structural effects of water clusters on viscosity at high shear rates. J Chem Phys 2024; 160:104502. [PMID: 38456533 DOI: 10.1063/5.0187906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Accepted: 02/19/2024] [Indexed: 03/09/2024] Open
Abstract
In this study, we use molecular dynamics simulations of liquid water to investigate how shear thinning affects the viscosity of liquid water by structural changes of the hydrogen bond network. The effect of shear on viscosity can be divided into two parts: shear-induced destruction of the hydrogen bond network and the influence of the water structure on shear viscosity. First, strong shear destroys tetrahedral structures and thus reduces the connectivity of the hydrogen bond network. It is mainly because shear deformation, characterized by compression and expansion axes, respectively, triggers the destruction and formation of hydrogen bonds, resulting in anisotropic effects on water structures. At the same time, shear destroys large clusters and enhances the formation of small ones, resulting in a decrease in average cluster sizes. Second, the change of viscosity obeys a power law relationship with the change of hydrogen bond structures, highlighting a one-to-one correspondence between structure and property. Meanwhile, in order to explain why the structure affects viscosity, we define hydrogen-bond viscosity and find that the cooperative motion of the water structures can promote momentum transfer in the form of aggregations. Hydrogen-bond viscosity accounts for 5%-50% of the total viscosity. Our results elucidate that water structures are the important structural units to explain the change of water properties.
Collapse
Affiliation(s)
- Yitian Gao
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Jian Wu
- Department of Physics, Tsinghua University, Beijing 100084, China
| | - Yixuan Feng
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Jiale Han
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Hongwei Fang
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
15
|
Novak N, Liang X, Kontogeorgis GM. Prediction of water anomalous properties by introducing the two-state theory in SAFT. J Chem Phys 2024; 160:104505. [PMID: 38465683 DOI: 10.1063/5.0186752] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Accepted: 02/14/2024] [Indexed: 03/12/2024] Open
Abstract
Water is one of the most abundant substances on earth, but it is still not entirely understood. It shows unusual behavior, and its properties present characteristic extrema unlike any other fluid. This unusual behavior has been linked to the two-state theory of water, which proposes that water forms different clusters, one with a high density and one with a low density, which may even form two distinct phases at low temperatures. Models incorporating the two-state theory manage to capture the unusual extrema of water, unlike traditional equations of state, which fail. In this work, we have derived the framework to incorporate the two-state theory of water into the Statistical-Associating-Fluid-Theory (SAFT). More specifically, we have assumed that water is an ideal solution of high density water molecules and low density water molecules that are in chemical equilibrium. Using this assumption, we have generalized the association term SAFT to allow for the simultaneous existence of the two water types, which have the same physical parameters but different association properties. We have incorporated the newly derived association term in the context of the Perturbed Chain-SAFT (PC-SAFT). The new model is referred to as PC-SAFT-Two-State (PC-SAFT-TS). Using PC-SAFT-TS, we have succeeded in predicting the characteristic extrema of water, such as its density and speed of sound maximum, etc., without loss of accuracy compared to the original PC-SAFT. This new framework is readily extended to mixtures, and PC-SAFT-TS manages to capture the solubility minimum of hydrocarbons in water in a straightforward manner.
Collapse
Affiliation(s)
- Nefeli Novak
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, 2800 Kgs Lyngby, Denmark
| | - Xiaodong Liang
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, 2800 Kgs Lyngby, Denmark
| | - Georgios M Kontogeorgis
- Center for Energy Resources Engineering, Department of Chemical and Biochemical Engineering, Technical University of Denmark, 2800 Kgs Lyngby, Denmark
| |
Collapse
|
16
|
Grabowska J, Kuffel A, Zielkiewicz J. Long-range, water-mediated interaction between a moderately active antifreeze protein molecule and the surface of ice. J Chem Phys 2024; 160:095101. [PMID: 38445741 DOI: 10.1063/5.0187663] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Accepted: 02/16/2024] [Indexed: 03/07/2024] Open
Abstract
Using molecular dynamics simulations, we show that a molecule of moderately active antifreeze protein (type III AFP, QAE HPLC-12 isoform) is able to interact with ice in an indirect manner. This interaction occurs between the ice binding site (IBS) of the AFP III molecule and the surface of ice, and it is mediated by liquid water, which separates these surfaces. As a result, the AFP III molecule positions itself at a specific orientation and distance relative to the surface of ice, which enables the effective binding (via hydrogen bonds) of the molecule with the nascent ice surface. Our results show that the final adsorption of the AFP III molecule on the surface of ice is not achieved by chaotic diffusion movements, but it is preceded by a remote, water-mediated interaction between the IBS and the surface of ice. The key factor that determines the existence of this interaction is the ability of water molecules to spontaneously form large, high-volume aggregates that can be anchored to both the IBS of the AFP molecule and the surface of ice. The results presented in this work for AFP III are in full agreement with the ones obtained by us previously for hyperactive CfAFP, which indicates that the mechanism of the remote interaction of these molecules with ice remains unchanged despite significant differences in the molecular structure of their ice binding sites. For that reason, we can expect that also other types of AFPs interact with the ice surface according to an analogous mechanism.
Collapse
Affiliation(s)
- Joanna Grabowska
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| | - Anna Kuffel
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| | - Jan Zielkiewicz
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| |
Collapse
|
17
|
Maturi F, Raposo Filho RS, Brites CDS, Fan J, He R, Zhuang B, Liu X, Carlos LD. Deciphering Density Fluctuations in the Hydration Water of Brownian Nanoparticles via Upconversion Thermometry. J Phys Chem Lett 2024; 15:2606-2615. [PMID: 38420927 PMCID: PMC10926164 DOI: 10.1021/acs.jpclett.4c00044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 02/14/2024] [Accepted: 02/16/2024] [Indexed: 03/02/2024]
Abstract
We investigate the intricate relationship among temperature, pH, and Brownian velocity in a range of differently sized upconversion nanoparticles (UCNPs) dispersed in water. These UCNPs, acting as nanorulers, offer insights into assessing the relative proportion of high-density and low-density liquid in the surrounding hydration water. The study reveals a size-dependent reduction in the onset temperature of liquid-water fluctuations, indicating an augmented presence of high-density liquid domains at the nanoparticle surfaces. The observed upper-temperature threshold is consistent with a hypothetical phase diagram of water, validating the two-state model. Moreover, an increase in pH disrupts the organization of water molecules, similar to external pressure effects, allowing simulation of the effects of temperature and pressure on hydrogen bonding networks. The findings underscore the significance of the surface of suspended nanoparticles for understanding high- to low-density liquid fluctuations and water behavior at charged interfaces.
Collapse
Affiliation(s)
- Fernando
E. Maturi
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
- Institute
of Chemistry, São Paulo State University
(UNESP), 14800-060 Araraquara, SP, Brazil
| | - Ramon S. Raposo Filho
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Carlos D. S. Brites
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| | - Jingyue Fan
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Ruihua He
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Bilin Zhuang
- Harvey
Mudd College, 301 Platt
Boulevard, Claremont, California 91711, United States
| | - Xiaogang Liu
- Department
of Chemistry, National University of Singapore, Singapore 117543
| | - Luís D. Carlos
- Phantom-g,
CICECO - Aveiro Institute of Materials, Department of Physics, University of Aveiro, 3810-193 Aveiro, Portugal
| |
Collapse
|
18
|
Foffi R, Sciortino F. Identification of local structures in water from supercooled to ambient conditions. J Chem Phys 2024; 160:094504. [PMID: 38436442 DOI: 10.1063/5.0188764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 02/09/2024] [Indexed: 03/05/2024] Open
Abstract
Studies of water thermodynamics have long been tied to the identification of two distinct families of local structures, whose competition could explain the origin of the many thermodynamic anomalies and the hypothesized liquid-liquid critical point in water. Despite the many successes and insights gained, the structural indicators proposed throughout the years were not able to unequivocally identify these two families over a wide range of conditions. We show that a recently introduced indicator, Ψ, which exploits information on the hydrogen bond network connectivity, can reliably identify these two distinct local environments over a wide range of thermodynamic conditions (188-300 K and 0-13 kbar) and that close to the liquid-liquid critical point, the spatial correlations of density fluctuations are identical to those of the Ψ indicator. Our results strongly support the idea that water thermodynamic properties arise from the competition between two distinct and identifiable local environments.
Collapse
Affiliation(s)
- Riccardo Foffi
- Department of Civil, Environmental and Geomatic Engineering, Institute for Environmental Engineering, ETH Zürich, Laura-Hezner-Weg 7, 8093 Zürich, Switzerland
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Rome, Italy
| |
Collapse
|
19
|
Finney JL. The structure of water: A historical perspective. J Chem Phys 2024; 160:060901. [PMID: 38341786 DOI: 10.1063/5.0182665] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/11/2024] [Indexed: 02/13/2024] Open
Abstract
Attempts to understand the molecular structure of water were first made well over a century ago. Looking back at the various attempts, it is illuminating to see how these were conditioned by the state of knowledge of chemistry and physics at the time and the experimental and theoretical tools then available. Progress in the intervening years has been facilitated by not only conceptual and theoretical advances in physics and chemistry but also the development of experimental techniques and instrumentation. Exploitation of powerful computational methods in interpreting what at first sight may seem impenetrable experimental data has led us to the consistent and detailed picture we have today of not only the structure of liquid water itself and how it changes with temperature and pressure but also its interactions with other molecules, in particular those relevant to water's role in important chemical and biological processes. Much remains to be done in the latter areas, but the experimental and computational techniques that now enable us to do what might reasonably be termed "liquid state crystallography" have opened the door to make possible further advances. Consequently, we now have the tools to explore further the role of water in those processes that underpin life itself-the very prospect that inspired Bernal to develop his ideas on the structure of liquids in general and of water in particular.
Collapse
Affiliation(s)
- John L Finney
- Department of Physics and Astronomy and London Centre for Nanotechnology, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
20
|
Santis GD, Herman KM, Heindel JP, Xantheas SS. Descriptors of water aggregation. J Chem Phys 2024; 160:054306. [PMID: 38341703 DOI: 10.1063/5.0179815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Accepted: 01/05/2024] [Indexed: 02/13/2024] Open
Abstract
We rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (rOH) and HOH bond angle (θHOH), (ii) two-body: OO distance (rOO), OHO angle (θOHO), and HOOX dihedral angle (ϕHOOX), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θOOO), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (qm, m = 2, 3, 4, 5) and radius of gyration (Rg). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles. We mined a previously reported database of 4 948 959 water cluster minima for (H2O)n, n = 3-25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that rOH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.
Collapse
Affiliation(s)
- Garrett D Santis
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MSIN J7-10, Richland, Washington 99352, USA
- Computational and Theoretical Chemistry Institute (CTCI), Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
21
|
Fransson T, Pettersson LGM. Calibrating TDDFT Calculations of the X-ray Emission Spectrum of Liquid Water: The Effects of Hartree-Fock Exchange. J Chem Theory Comput 2023; 19:7333-7342. [PMID: 37787584 PMCID: PMC10601479 DOI: 10.1021/acs.jctc.3c00728] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Indexed: 10/04/2023]
Abstract
The structure and dynamics of liquid water continue to be debated, with insight provided by, among others, X-ray emission spectroscopy (XES), which shows a split in the high-energy 1b1 feature. This split is yet to be reproduced by theory, and it remains unclear if these difficulties are related to inaccuracies in dynamics simulations, spectrum calculations, or both. We investigate the performance of different methods for calculating XES of liquid water, focusing on the ability of time-dependent density functional theory (TDDFT) to reproduce reference spectra obtained by high-level coupled cluster and algebraic-diagrammatic construction scheme calculations. A metric for evaluating the agreement between theoretical spectra termed the integrated absolute difference (IAD), which considers the integral of shifted difference spectra, is introduced and used to investigate the performance of different exchange-correlation functionals. We find that computed spectra of symmetric and asymmetric model water structures are strongly and differently influenced by the amount of Hartree-Fock exchange, with best agreement to reference spectra for ∼40-50%. Lower percentages tend to yield high density of contributing states, resulting in too broad features. The method introduced here is useful also for other spectrum calculations, in particular where the performance for ensembles of structures are evaluated.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| |
Collapse
|
22
|
Fijan D, Wilson M. Thermodynamic anomalies, polyamorphism and all that. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220336. [PMID: 37634531 PMCID: PMC10460645 DOI: 10.1098/rsta.2022.0336] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2023] [Accepted: 05/18/2023] [Indexed: 08/29/2023]
Abstract
The appearance and evolution of thermodynamics anomalies, and related properties, are studied for two classes of system, modelling those dominated by covalent and ionic interactions, respectively. Such anomalies are most familiar in the density but are also present in other thermodynamics variables such as the compressibility and heat capacity. By systematically varying key model parameters the emergence and evolution of these anomalies can be tracked across the phase space. The interaction of the anomalies can often be rationalized by thermodynamics 'rules'. The emergence of these anomalies may also be correlated with the appearance of polyamorphism, the existence of multiple amorphous states which differ in density and entropy. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.
Collapse
Affiliation(s)
- Domagoj Fijan
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford, UK
| | - Mark Wilson
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford, UK
| |
Collapse
|
23
|
Škvára J, Nezbeda I, Urbic T. Supercooled water in two dimensions: Structure and thermodynamics of the Mercedes-Benz model. J Mol Liq 2023; 386:122445. [PMID: 37435361 PMCID: PMC10331298 DOI: 10.1016/j.molliq.2023.122445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2023]
Abstract
The two-dimensional Mercedes-Benz model of water has been studied by molecular simulations over a wide range of thermodynamic conditions as an attempt to locate the supercooled region where a liquid-liquid separation and, potentially, also other structures may occur. Both the correlation functions and a number of local structure factors have been used to identify different structural arrangements. These include, in addition to the hexatic phase, also the hexagon, pentagon, and quadruplet arrangements. All these structures result from the competition between the hydrogen bonding and Lennard-Jones interactions and their effect at different temperatures and pressures. Based on the obtained results, an attempt is made to sketch a (rather complex) phase diagram of the model.
Collapse
Affiliation(s)
- Jiří Škvára
- Institute of Chemical Process Fundamentals, Czech Acad. Sci., 16502 Prague, Czech Republic
- Faculty of Science, J. E. Purkinje University, 40096 Ústí nad Labem, Czech Republic
| | - Ivo Nezbeda
- Faculty of Science, J. E. Purkinje University, 40096 Ústí nad Labem, Czech Republic
| | - Tomaz Urbic
- University of Ljubljana, Faculty of Chemistry and Chemical Technology, Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
24
|
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
|
25
|
Zielkiewicz J. Mechanism of antifreeze protein functioning and the "anchored clathrate water" concept. J Chem Phys 2023; 159:085101. [PMID: 37622597 DOI: 10.1063/5.0158590] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Accepted: 08/04/2023] [Indexed: 08/26/2023] Open
Abstract
In liquid water, there is a natural tendency to form aggregates that consist of water molecules linked by hydrogen bonds. Such spontaneously formed aggregates are surrounded by a "sea" of disordered water molecules, with both forms remaining in equilibrium. The process of creating water aggregates also takes place in the solvation water of proteins, but in this case, the interactions of water molecules with the protein surface shift the equilibrium of the process. In this paper, we analyze the structural properties of the solvation water in antifreeze proteins (AFPs). The results of molecular dynamics analysis with the use of various parameters related to the structure of solvation water on the protein surface are presented. We found that in the vicinity of the active region responsible for the binding of AFPs to ice, the equilibrium is clearly shifted toward the formation of "ice-like aggregates," and the solvation water has a more ordered ice-like structure. We have demonstrated that a reduction in the tendency to create "ice-like aggregates" results in a significant reduction in the antifreeze activity of the protein. We conclude that shifting the equilibrium in favor of the formation of "ice-like aggregates" in the solvation water in the active region is a prerequisite for the biological functionality of AFPs, at least for AFPs having a well-defined ice binding area. In addition, our results fully confirm the validity of the "anchored clathrate water" concept, formulated by Garnham et al. [Proc. Natl. Acad. Sci. U. S. A. 108, 7363 (2011)].
Collapse
Affiliation(s)
- Jan Zielkiewicz
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| |
Collapse
|
26
|
Loubet NA, Verde AR, Lockhart JA, Appignanesi GA. Turning an energy-based defect detector into a multi-molecule structural indicator for water. J Chem Phys 2023; 159:064512. [PMID: 37578063 DOI: 10.1063/5.0159060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2023] [Accepted: 07/27/2023] [Indexed: 08/15/2023] Open
Abstract
Recent studies have provided conclusive evidence for the existence of a liquid-liquid critical point in numerical models of water. Such a scenario implies the competition between two local molecular arrangements of different densities: a high-density liquid (HDL) and a low-density liquid (LDL). Within this context, the development of accurate structural indicators to properly characterize the two interconverting local structures is demanded. In a previous study, we introduced a reliable energy-based structural descriptor that properly discriminates water molecules into tetrahedrally arranged molecules (T molecules) and distorted molecules (D molecules). The latter constitute defects in terms of hydrogen bond (HB) coordination and have been shown to represent a minority component, even at high temperatures above the melting point. In addition, the D molecules tend to form high-quality HBs with three T molecules and to be surrounded by T and D molecules at further distances. Thus, it became evident that, while the LDL state might consist of a virtually pure T state, the HDL state would comprise mixed molecular arrangements including the D molecules. Such a need to abandon the single-molecule description requires the investigation of the degree of structural information to be incorporated in order to build an appropriate multi-molecule indicator. Hence, in this work, we shall study the effect of the local structural constraints on the water molecules in order to discriminate the different molecular arrangements into two disjoint classes. This will enable us to build a multi-molecule structural indicator for water whose performance will then be investigated within the water's supercooled regime.
Collapse
Affiliation(s)
- Nicolás A Loubet
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Alejandro R Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Jano A Lockhart
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
27
|
Chen BWJ, Zhang X, Zhang J. Accelerating explicit solvent models of heterogeneous catalysts with machine learning interatomic potentials. Chem Sci 2023; 14:8338-8354. [PMID: 37564405 PMCID: PMC10411631 DOI: 10.1039/d3sc02482b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Accepted: 07/11/2023] [Indexed: 08/12/2023] Open
Abstract
Realistically modelling how solvents affect catalytic reactions is a longstanding challenge due to its prohibitive computational cost. Typically, an explicit atomistic treatment of the solvent molecules is needed together with molecular dynamics (MD) simulations and enhanced sampling methods. Here, we demonstrate the utility of machine learning interatomic potentials (MLIPs), coupled with active learning, to enable fast and accurate explicit solvent modelling of adsorption and reactions on heterogeneous catalysts. MLIPs trained on-the-fly were able to accelerate ab initio MD simulations by up to 4 orders of magnitude while reproducing with high fidelity the geometrical features of water in the bulk and at metal-water interfaces. Using these ML-accelerated simulations, we accurately predicted key catalytic quantities such as the adsorption energies of CO*, OH*, COH*, HCO*, and OCCHO* on Cu surfaces and the free energy barriers of C-H scission of ethylene glycol over Cu and Pd surfaces, as validated with ab initio calculations. We envision that such simulations will pave the way towards detailed and realistic studies of solvated catalysts at large time- and length-scales.
Collapse
Affiliation(s)
- Benjamin W J Chen
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, #16-16 Connexis Singapore 138632 Singapore
| | - Xinglong Zhang
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, #16-16 Connexis Singapore 138632 Singapore
| | - Jia Zhang
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, #16-16 Connexis Singapore 138632 Singapore
| |
Collapse
|
28
|
Eltareb A, Lopez GE, Giovambattista N. The Importance of Nuclear Quantum Effects on the Thermodynamic and Structural Properties of Low-Density Amorphous Ice: A Comparison with Hexagonal Ice. J Phys Chem B 2023; 127:4633-4645. [PMID: 37178124 DOI: 10.1021/acs.jpcb.3c01025] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We study the nuclear quantum effects (NQE) on the thermodynamic properties of low-density amorphous ice (LDA) and hexagonal ice (Ih) at P = 0.1 MPa and T ≥ 25 K. Our results are based on path-integral molecular dynamics (PIMD) and classical MD simulations of H2O and D2O using the q-TIP4P/F water model. We show that the inclusion of NQE is necessary to reproduce the experimental properties of LDA and ice Ih. While MD simulations (no NQE) predict that the density ρ(T) of LDA and ice Ih increases monotonically upon cooling, PIMD simulations indicate the presence of a density maximum in LDA and ice Ih. MD and PIMD simulations also predict a qualitatively different T-dependence for the thermal expansion coefficient αP(T) and bulk modulus B(T) of both LDA and ice Ih. Remarkably, the ρ(T), αP(T), and B(T) of LDA are practically identical to those of ice Ih. The origin of the observed NQE is due to the delocalization of the H atoms, which is identical in LDA and ice Ih. H atoms delocalize considerably (over a distance ≈ 20-25% of the OH covalent-bond length) and anisotropically (preferentially perpendicular to the OH covalent bond), leading to less linear hydrogen bonds HB (larger HOO angles and longer OO separations) than observed in classical MD simulations.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, New York 11210, United States
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Gustavo E Lopez
- Department of Chemistry, Lehman College of the City University of New York, Bronx, New York 10468, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Nicolas Giovambattista
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, New York 11210, United States
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, New York 10016, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| |
Collapse
|
29
|
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
|
30
|
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
|
31
|
Habitzreuter MA, Barbosa MC. Density anomaly in water-alcohol mixtures: Minimum model for structure makers and breakers. Phys Rev E 2023; 107:034601. [PMID: 37072948 DOI: 10.1103/physreve.107.034601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 02/21/2023] [Indexed: 04/20/2023]
Abstract
We modeled the change in the temperature of maximum density (TMD) of a waterlike solvent when small amounts of solute are added to the mixture. The solvent is modeled as a two length scales potential, which is known to exhibit waterlike characteristic anomalies, while the solute is chosen to have an attractive interaction with the solvent which is tuned from small to large attractive potential. We show that if the solute exhibits high attraction with the solvent it behaves as a structure maker and the TMD increases with the addition of solute, while if the solute shows a low attraction with the solvent the TMD decreases, with the solute behaving as a structure breaker.
Collapse
Affiliation(s)
- Marco A Habitzreuter
- Instituto de Física - Universidade Federal do Rio Grande do Sul Porto Alegre, Rio Grande do Sul 91501-970, Brazil
| | - Marcia C Barbosa
- Instituto de Física - Universidade Federal do Rio Grande do Sul Porto Alegre, Rio Grande do Sul 91501-970, Brazil
| |
Collapse
|
32
|
Gorfer A, Dellago C, Sega M. High-density liquid (HDL) adsorption at the supercooled water/vapor interface and its possible relation to the second surface tension inflection point. J Chem Phys 2023; 158:054503. [PMID: 36754827 DOI: 10.1063/5.0132985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
We investigate the properties of water along the liquid/vapor coexistence line in the supercooled regime down to the no-man's land. Extensive molecular dynamics simulations of the TIP4P/2005 liquid/vapor interface in the range 198-348 K allow us to locate the second surface tension inflection point with a high accuracy at 283 ± 5 K, close to the temperature of maximum density. This temperature also coincides with the appearance of a density anomaly at the interface known as the apophysis. We relate the emergence of the apophysis to the observation of high-density liquid (HDL) water adsorption in the proximity of the liquid/vapor interface.
Collapse
Affiliation(s)
- Alexander Gorfer
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, Wien A-1090, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, Wien A-1090, Austria
| | - Marcello Sega
- Department of Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
| |
Collapse
|
33
|
Takahashi O, Pettersson LGM. Dynamical and interference effects in X-ray emission spectroscopy of H-bonded water – origin of the split lone-pair peaks. Mol Phys 2023. [DOI: 10.1080/00268976.2023.2170686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Osamu Takahashi
- Basic Chemistry Program, Graduate School of Advanced Science and Engineering, Hiroshima University, Higashi-Hiroshima, Japan
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm, Sweden
| |
Collapse
|
34
|
Jedrecy A, Saitta AM, Pietrucci F. Free energy calculations and unbiased molecular dynamics targeting the liquid-liquid transition in water no man's land. J Chem Phys 2023; 158:014502. [PMID: 36610960 DOI: 10.1063/5.0120789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The existence of a first-order phase transition between a low-density liquid (LDL) and a high-density liquid (HDL) form of supercooled water has been a central and highly debated issue of physics and chemistry for the last three decades. We present a computational study that allows us to determine the free-energy landscapes of supercooled water over a wide range of pressure and temperature conditions using the TIP4P/2005 force field. Our approach combines topology-based structural transformation coordinates, state-of-the-art free-energy calculation methods, and extensive unbiased molecular dynamics. All our diverse simulations cannot detect any barrier within the investigated timescales and system size, for a discontinuous transition between the LDL and HDL forms throughout the so-called "no man's land," until the onset of the solid, non-diffusive amorphous forms.
Collapse
Affiliation(s)
- Alexandre Jedrecy
- Insitut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, CNRS, MNHN, UMR 7590, Paris, France
| | - A Marco Saitta
- Insitut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, CNRS, MNHN, UMR 7590, Paris, France
| | - Fabio Pietrucci
- Insitut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, CNRS, MNHN, UMR 7590, Paris, France
| |
Collapse
|
35
|
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
|
36
|
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]
|
37
|
Š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]
|
38
|
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
|
39
|
Tsochantaris E, Muthachikavil AV, Peng B, Liang X, Kontogeorgis GM. Multiple insights call for revision of modern thermodynamic models to account for structural fluctuations in water. AIChE J 2022; 68:e17891. [PMID: 36591369 PMCID: PMC9787682 DOI: 10.1002/aic.17891] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Revised: 06/29/2022] [Accepted: 08/15/2022] [Indexed: 01/04/2023]
Abstract
Modern thermodynamic models incorporate the concept of association (hydrogen bonding) and they can describe very satisfactorily many properties of water containing mixtures. They have not been successful in representing water's anomalous properties and this work provides a possible explanation. We have analyzed and interpreted recent experimental data, molecular simulation results, and two-state theory approaches and compared against the predictions from thermodynamic models. We show that the dominance of the tetrahedral structure implemented in modern thermodynamic models may be the reason for their failure for describing water systems. While this study does not prove the two-state theories for water, it indicates that a high level of tetrahedral structure of water is not in agreement with water's anomalous properties when used in thermodynamic models.
Collapse
Affiliation(s)
- Evangelos Tsochantaris
- Department of Chemical and Biochemical Engineering, Center for Energy Resources EngineeringTechnical University of DenmarkLyngbyDenmark
| | - Aswin V. Muthachikavil
- Department of Chemical and Biochemical Engineering, Center for Energy Resources EngineeringTechnical University of DenmarkLyngbyDenmark
| | - Baoliang Peng
- Research Institute of Petroleum Exploration & Development (RIPED), PetroChinaBeijingChina
| | - Xiaodong Liang
- Department of Chemical and Biochemical Engineering, Center for Energy Resources EngineeringTechnical University of DenmarkLyngbyDenmark
| | - Georgios M. Kontogeorgis
- Department of Chemical and Biochemical Engineering, Center for Energy Resources EngineeringTechnical University of DenmarkLyngbyDenmark
| |
Collapse
|
40
|
Jin Z, Zhao J, Chen G, Chen G, Luo Z, Xu L. Revealing the three-component structure of water with principal component analysis (PCA) of X-ray spectra. SOFT MATTER 2022; 18:7486-7496. [PMID: 36000526 DOI: 10.1039/d2sm00576j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Combining principal component analysis (PCA) of X-ray spectra with MD simulations, we experimentally reveal the existence of three basic components in water. These components exhibit distinct structures, densities, and temperature dependencies. Among the three, the two major components correspond to the low-density liquid (LDL) and the high-density liquid (HDL) predicted by the two-component model, and the third component exhibits a unique 5-hydrogen-bond configuration with ultra-high local density. As the temperature increases, the LDL component decreases and the HDL component increases, while the third component varies non-monotonically with a peak around 20 °C to 30 °C. The 3D structure of the third component is further illustrated as the uniform distribution of five hydrogen-bonded neighbors on a spherical surface. Our study reveals experimental evidence for water's possible three-component structure, which provides a fundamental basis for understanding water's special properties and anomalies.
Collapse
Affiliation(s)
- Zhipeng Jin
- Department of Physics, the Chinese University of Hong Kong, Hong Kong, China.
| | - Jiangtao Zhao
- National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230026, China.
| | - Gang Chen
- School of Physical Science and Technology, ShanghaiTech University, Shanghai 201210, China
| | - Guo Chen
- Chongqing Key Laboratory of Soft Condensed Matter Physics and Smart Materials, College of Physics, Chongqing University, Chongqing 400044, China.
| | - Zhenlin Luo
- National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230026, China.
| | - Lei Xu
- Department of Physics, the Chinese University of Hong Kong, Hong Kong, China.
| |
Collapse
|
41
|
Lee E, Baiz CR. How cryoprotectants work: hydrogen-bonding in low-temperature vitrified solutions. Chem Sci 2022; 13:9980-9984. [PMID: 36128234 PMCID: PMC9430440 DOI: 10.1039/d2sc03188d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Accepted: 08/05/2022] [Indexed: 11/30/2022] Open
Abstract
Dimethyl sulfoxide (DMSO) increases cell and tissue viability at low temperatures and is commonly used as a cryoprotectant for cryogenic storage of biological materials. DMSO disorders the water hydrogen-bond networks and inhibits ice-crystal growth, though the specific DMSO interactions with water are difficult to characterize. In this study, we use a combination of Fourier Transform infrared spectroscopy (FTIR), molecular dynamics simulations, and vibrational frequency maps to characterize the temperature-dependent hydrogen bonding interactions of DMSO with water from 30 °C to -80 °C. Specifically, broad peaks in O-D stretch vibrational spectra of DMSO and deuterated water (HDO) cosolvent systems show that the hydrogen bond networks become increasingly disrupted compared to pure water. Simulations demonstrated that these disrupted hydrogen bond networks remain largely localized to the first hydration shell of DMSO, which explains the high DMSO concentrations needed to prevent ice crystal formation in cryopreservation applications.
Collapse
Affiliation(s)
- Euihyun Lee
- Department of Chemistry, The University of Texas at Austin Austin TX 78712 USA
| | - Carlos R Baiz
- Department of Chemistry, The University of Texas at Austin Austin TX 78712 USA
| |
Collapse
|
42
|
Using Car-Parrinello simulations and microscopic order descriptors to reveal two locally favored structures with distinct molecular dipole moments and dynamics in ambient liquid water. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119936] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
43
|
Martelli F. Steady-like topology of the dynamical hydrogen bond network in supercooled water. PNAS NEXUS 2022; 1:pgac090. [PMID: 36741425 PMCID: PMC9896910 DOI: 10.1093/pnasnexus/pgac090] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 06/08/2022] [Indexed: 02/07/2023]
Abstract
We investigate the link between topology of the hydrogen bond network (HBN) and large-scale density fluctuations in water from ambient conditions to the glassy state. We observe a transition from a temperature-dependent topology at high temperatures, to a steady-like topology below the Widom temperature TW ∼ 220 K signaling the fragile-to-strong crossover and the maximum in structural fluctuations. As a consequence of the steady topology, the network suppresses large-scale density fluctuations much more efficiently than at higher temperatures. Below TW , the contribution of coordination defects of the kind A 2 D 1 (two acceptors and one donor) to the kinetics of the HBN becomes progressively more pronounced, suggesting that A 2 D 1 configurations may represent the main source of dynamical heterogeneities. Below the vitrification temperature, the freezing of rotational and translational degrees of freedom allow for an enhanced suppression of large-scale density fluctuations and the sample reaches the edges of nearly hyperuniformity. The formed network still hosts coordination defects, hence implying that nearly hyperuniformity goes beyond the classical continuous random network paradigm of tetrahedral networks and can emerge in scenarios much more complex than previously assumed. Our results unveil a hitherto undisclosed link between network topology and properties of water essential for better understanding water's rich and complex nature. Beyond implications for water, our findings pave the way to a better understanding of the physics of supercooled liquids and disordered hyperuniform networks at large.
Collapse
|
44
|
Tsai E, Gallage Dona HK, Tong X, Du P, Novak B, David R, Rick SW, Zhang D, Kumar R. Unraveling the Role of Charge Patterning in the Micellar Structure of Sequence-Defined Amphiphilic Peptoid Oligomers by Molecular Dynamics Simulations. Macromolecules 2022; 55:5197-5212. [PMID: 35784657 PMCID: PMC9245439 DOI: 10.1021/acs.macromol.2c00141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2022] [Revised: 05/26/2022] [Indexed: 11/28/2022]
Abstract
![]()
Electrostatic interactions
play a significant role in regulating
biological systems and have received increasing attention due to their
usefulness in designing advanced stimulus-responsive materials. Polypeptoids
are highly tunable N-substituted peptidomimetic polymers that lack
backbone hydrogen bonding and chirality. Therefore, polypeptoids are
suitable systems to study the effect of noncovalent interactions of
substituents without complications of backbone intramolecular and
intermolecular hydrogen bonding. In this study, all-atom molecular
dynamics (MD) simulations were performed on micelles formed by a series
of sequence-defined ionic polypeptoid block copolymers consisting
of a hydrophobic segment and a hydrophilic segment in an aqueous solution.
By combining the results from MD simulations and experimental small-angle
neutron scattering data, further insights were gained into the internal
structure of the formed polypeptoid micelles, which is not always
directly accessible from experiments. In addition, information was
gained into the physics of the noncovalent interactions responsible
for the self-assembly of weakly charged polypeptoids in an aqueous
solution. While the aggregation number is governed by electrostatic
repulsion of the negatively charged carboxylate (COO–) substituents on the polypeptoid chain within the micelle, MD simulations
indicate that the position of the charge on singly charged chains
mediates the shape of the micelle through the charge–dipole
interactions between the COO– substituent and the
surrounding water. Therefore, the polypeptoid micelles formed from
the single-charged series offer the possibility for tailorable micelle
shapes. In contrast, the polypeptoid micelles formed from the triple-charged
series are characterized by more pronounced electrostatic repulsion
that competes with more significant charge–sodium interactions,
making it difficult to predict the shape of the micelles. This work
has helped further develop design principles for the shape and structure
of self-assembled micelles by controlling the position of charged
moieties on the backbone of polypeptoid block copolymers.
Collapse
Affiliation(s)
- Erin Tsai
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | | | - Xinjie Tong
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Pu Du
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Brian Novak
- Department of Mechanical and Industrial Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Rolf David
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Steven W. Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, United States
| | - Donghui Zhang
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Revati Kumar
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
- Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
45
|
Herrero C, Pauletti M, Tocci G, Iannuzzi M, Joly L. Connection between water's dynamical and structural properties: Insights from ab initio simulations. Proc Natl Acad Sci U S A 2022; 119:e2121641119. [PMID: 35588447 PMCID: PMC9173753 DOI: 10.1073/pnas.2121641119] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 04/12/2022] [Indexed: 01/25/2023] Open
Abstract
SignificanceFirst-principles calculations, which explicitly account for the electronic structure of matter, can shed light on the molecular structure and dynamics of water in its supercooled state. In this work, we use density functional theory, which relies on a functional to describe electronic exchange and correlations, to evaluate which functional best describes the temperature evolution of bulk water transport coefficients. We also assess the validity of the Stokes-Einstein relation for all the functionals in the temperature range studied, and explore the link between structure and dynamics. Based on these results, we show how transport coefficients can be computed from structural descriptors, which require shorter simulation times to converge, and we point toward strategies to develop better functionals.
Collapse
Affiliation(s)
- Cecilia Herrero
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Michela Pauletti
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
- Institut Universitaire de France (IUF), 75005 Paris, France
| |
Collapse
|
46
|
Gao Y, Fang H, Ni K, Feng Y. Water clusters and density fluctuations in liquid water based on extended hierarchical clustering methods. Sci Rep 2022; 12:8036. [PMID: 35577839 PMCID: PMC9110331 DOI: 10.1038/s41598-022-11947-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2022] [Accepted: 04/25/2022] [Indexed: 11/29/2022] Open
Abstract
The microscopic structures of liquid water at ambient temperatures remain a hot debate, which relates with structural and density fluctuations in the hydrogen bond network. Here, we use molecular dynamics simulations of liquid water to study the properties of three-dimensional cage-like water clusters, which we investigate using extended graph-based hierarchical clustering methods. The water clusters can cover over 95% of hydrogen bond network, among which some clusters maximally encompass thousands of molecules extending beyond 3.0 nm. The clusters imply fractal behaviors forming percolating networks and the morphologies of small and large clusters show different scaling rules. The local favored clusters and the preferred connections between adjacent clusters correspond to lower energy and conformational entropy depending on cluster topologies. Temperature can destroy large clusters into small ones. We show further that the interior of clusters favors high-density patches. The water molecules in the small clusters, inside which are the void regarded as hydrophobic objects, have a preference for being more tetrahedral. Our results highlight the properties and changes of water clusters as the fundamental building blocks of hydrogen bond networks. In addition, the water clusters can elucidate structural and density fluctuations on different length scales in liquid water.
Collapse
|
47
|
Kirchner KA, Cassar DR, Zanotto ED, Ono M, Kim SH, Doss K, Bødker ML, Smedskjaer MM, Kohara S, Tang L, Bauchy M, Wilkinson CJ, Yang Y, Welch RS, Mancini M, Mauro JC. Beyond the Average: Spatial and Temporal Fluctuations in Oxide Glass-Forming Systems. Chem Rev 2022; 123:1774-1840. [PMID: 35511603 DOI: 10.1021/acs.chemrev.1c00974] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Atomic structure dictates the performance of all materials systems; the characteristic of disordered materials is the significance of spatial and temporal fluctuations on composition-structure-property-performance relationships. Glass has a disordered atomic arrangement, which induces localized distributions in physical properties that are conventionally defined by average values. Quantifying these statistical distributions (including variances, fluctuations, and heterogeneities) is necessary to describe the complexity of glass-forming systems. Only recently have rigorous theories been developed to predict heterogeneities to manipulate and optimize glass properties. This article provides a comprehensive review of experimental, computational, and theoretical approaches to characterize and demonstrate the effects of short-, medium-, and long-range statistical fluctuations on physical properties (e.g., thermodynamic, kinetic, mechanical, and optical) and processes (e.g., relaxation, crystallization, and phase separation), focusing primarily on commercially relevant oxide glasses. Rigorous investigations of fluctuations enable researchers to improve the fundamental understanding of the chemistry and physics governing glass-forming systems and optimize structure-property-performance relationships for next-generation technological applications of glass, including damage-resistant electronic displays, safer pharmaceutical vials to store and transport vaccines, and lower-attenuation fiber optics. We invite the reader to join us in exploring what can be discovered by going beyond the average.
Collapse
Affiliation(s)
- Katelyn A Kirchner
- Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Daniel R Cassar
- Department of Materials Engineering, Federal University of São Carlos, São Carlos, Sao Paulo 13565-905, Brazil
- Ilum School of Science, Brazilian Center for Research in Energy and Materials, Campinas, Sao Paulo 13083-970, Brazil
| | - Edgar D Zanotto
- Department of Materials Engineering, Federal University of São Carlos, São Carlos, Sao Paulo 13565-905, Brazil
| | - Madoka Ono
- Research Institute for Electronic Science, Hokkaido University, Sapporo, Hokkaido 001-0021, Japan
- Materials Integration Laboratories, AGC Incorporated, Yokohama, Kanagawa 230-0045, Japan
| | - Seong H Kim
- Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Karan Doss
- Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Mikkel L Bødker
- Department of Chemistry and Bioscience, Aalborg University, Aalborg 9220, Denmark
| | - Morten M Smedskjaer
- Department of Chemistry and Bioscience, Aalborg University, Aalborg 9220, Denmark
| | - Shinji Kohara
- Research Center for Advanced Measurement and Characterization National Institute for Materials Science, 1-2-1, Sengen, Tsukuba, Ibaraki 305-0047, Japan
| | - Longwen Tang
- Department of Civil and Environmental Engineering, University of California, Los Angeles, California 90095, United States
| | - Mathieu Bauchy
- Department of Civil and Environmental Engineering, University of California, Los Angeles, California 90095, United States
| | - Collin J Wilkinson
- Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
- Department of Research and Development, GlassWRX, Beaufort, South Carolina 29906, United States
| | - Yongjian Yang
- Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Rebecca S Welch
- Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Matthew Mancini
- Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - John C Mauro
- Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
48
|
Gartner TE, Hunter KM, Lambros E, Caruso A, Riera M, Medders GR, Panagiotopoulos AZ, Debenedetti PG, Paesani F. Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model. J Phys Chem Lett 2022; 13:3652-3658. [PMID: 35436129 DOI: 10.1021/acs.jpclett.2c00567] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
For the past 50 years, researchers have sought molecular models that can accurately reproduce water's microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ∼208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. This work suggests that MB-pol may provide predictive capability for studies of liquid water's physical properties across broad ranges of thermodynamic states, including the so-called water's "no man's land" which is difficult to probe experimentally.
Collapse
Affiliation(s)
- Thomas E Gartner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Kelly M Hunter
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Gregory R Medders
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | | | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California, San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
49
|
Kono Y, Ohara K, Kondo NM, Yamada H, Hiroi S, Noritake F, Nitta K, Sekizawa O, Higo Y, Tange Y, Yumoto H, Koyama T, Yamazaki H, Senba Y, Ohashi H, Goto S, Inoue I, Hayashi Y, Tamasaku K, Osaka T, Yamada J, Yabashi M. Experimental evidence of tetrahedral symmetry breaking in SiO 2 glass under pressure. Nat Commun 2022; 13:2292. [PMID: 35484122 PMCID: PMC9051114 DOI: 10.1038/s41467-022-30028-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Accepted: 04/11/2022] [Indexed: 11/09/2022] Open
Abstract
Bimodal behavior in the translational order of silicon's second shell in SiO2 liquid at high temperatures and high pressures has been recognized in theoretical studies, and the fraction of the S state with high tetrahedrality is considered as structural origin of the anomalous properties. However, it has not been well identified in experiment. Here we show experimental evidence of a bimodal behavior in the translational order of silicon's second shell in SiO2 glass under pressure. SiO2 glass shows tetrahedral symmetry structure with separation between the first and second shells of silicon at low pressures, which corresponds to the S state structure reported in SiO2 liquid. On the other hand, at high pressures, the silicon's second shell collapses onto the first shell, and more silicon atoms locate in the first shell. These observations indicate breaking of local tetrahedral symmetry in SiO2 glass under pressure, as well as SiO2 liquid.
Collapse
Affiliation(s)
- Yoshio Kono
- Geodynamics Research Center, Ehime University, 2-5 Bunkyo-cho, Matsuyama, 790-8577, Japan.
| | - Koji Ohara
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan
| | - Nozomi M Kondo
- Geodynamics Research Center, Ehime University, 2-5 Bunkyo-cho, Matsuyama, 790-8577, Japan
| | - Hiroki Yamada
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan
| | - Satoshi Hiroi
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan
| | - Fumiya Noritake
- Graduate Faculty of Interdisciplinary Research, University of Yamanashi, 4-3-11 Takeda, Kofu, Yamanashi, 400-8511, Japan
| | - Kiyofumi Nitta
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan
| | - Oki Sekizawa
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan
| | - Yuji Higo
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan
| | - Yoshinori Tange
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan
| | - Hirokatsu Yumoto
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Takahisa Koyama
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Hiroshi Yamazaki
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Yasunori Senba
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Haruhiko Ohashi
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Shunji Goto
- Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5198, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Ichiro Inoue
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Yujiro Hayashi
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Kenji Tamasaku
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Taito Osaka
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Jumpei Yamada
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Makina Yabashi
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| |
Collapse
|
50
|
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
|