Water at the nanoscale: From filling or dewetting hydrophobic pores and carbon nanotubes to "sliding" on graphene

In this work, we study the effect of nanoconfinement on the hydration properties of model hydrophobic pores and carbon nanotubes, determining their wetting propensity and the conditions for geometrically induced dehydration. By employing a recently introduced water structural index, we aim at two main goals: (1) to accurately quantify the local hydrophobicity and predict the drying transitions in such systems, and (2) to provide a molecular rationalization of the wetting process. In this sense, we will further discuss the number and strength of the interactions required by the water molecules to promote wetting. In the case of graphene-like surfaces, an explanation for their unexpectedly significant hydrophilicity will also be provided. On the one hand, the structural index will show that the net attraction to the dense carbon network that a water molecule experiences through several simultaneous weak interactions is sufficient to give rise to hydrophilic behavior. On the other hand, we will show that an additional effect is also at play: the hydrating water molecule is retained on the surface by a smooth exchange of such simultaneous weak interactions, as if "sliding" on graphene.

TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional

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.

In-layer inhomogeneity of molecular dynamics in quasi-liquid layers of ice

Quasi-liquid layers (QLLs) are present on the surface of ice and play a significant role in its distinctive chemical and physical properties. These layers exhibit considerable heterogeneity across different scales ranging from nanometers to millimeters. Although the formation of partially ice-like structures has been proposed, the molecular-level understanding of this heterogeneity remains unclear. Here, we examined the heterogeneity of molecular dynamics on QLLs based on molecular dynamics simulations and machine learning analysis of the simulation data. We demonstrated that the molecular dynamics of QLLs do not comprise a mixture of solid- and liquid water molecules. Rather, molecules having similar behaviors form dynamical domains that are associated with the dynamical heterogeneity of supercooled water. Nonetheless, molecules in the domains frequently switch their dynamical state. Furthermore, while there is no observable characteristic domain size, the long-range ordering strongly depends on the temperature and crystal face. Instead of a mixture of static solid- and liquid-like regions, our results indicate the presence of heterogeneous molecular dynamics in QLLs, which offers molecular-level insights into the surface properties of ice.

Existence of density inhomogeneity of liquid Te associated with liquid-liquid phase transition

We performed small-angle x-ray scattering measurements of liquid Te using a synchrotron radiation facility and observed maximum scattering intensity near 620 K in the supercooled region (melting temperature 723 K). This indicates that density inhomogeneity exists in liquid Te, and the fact that this temperature coincides with the temperature at which the specific heat, sound velocity, and thermal expansion coefficient reach their maxima means that this density inhomogeneity is the cause of these thermodynamic anomalies. The thermodynamic anomalies in liquid Te had already been shown in the 1980s to be comprehensively explained by the inhomogeneity associated with the continuous liquid-liquid phase transition (LLT), but direct experimental evidence for the existence of the inhomogeneity had not been obtained. The present results, together with those already obtained for mixture systems (Te-Se, Te-Ge), indicate the existence of inhomogeneity associated with LLT in liquid Te systems, and strongly support the model. Recently, similar maximum scattering intensity has also been observed in supercooled liquid water, which exhibits thermodynamic anomalies similar to those of Te, indicating the universality of the inhomogeneous model or LLT scenario to explain the thermodynamics of such 'anomalous liquids'. Further development of the LLT scenario is expected in near future.

Isotope effects in supercooled H2O and D2O and a corresponding-states-like rescaling of the temperature and pressure

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.

Unraveling the dynamic slowdown in supercooled water: The role of dynamic disorder in jump motions

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.

Liquid-liquid transition and ice crystallization in a machine-learned coarse-grained water model

Mounting experimental evidence supports the existence of a liquid-liquid transition (LLT) in high-pressure supercooled water. However, fast crystallization of supercooled water has impeded identification of the LLT line TLL(p) in experiments. While the most accurate all-atom (AA) water models display a LLT, their computational cost limits investigations of its interplay with ice formation. Coarse-grained (CG) models provide over 100-fold computational efficiency gain over AA models, enabling the study of water crystallization, but have not yet shown to have a LLT. Here, we demonstrate that the CG machine-learned water model Machine-Learned Bond-Order Potential (ML-BOP) has a LLT that ends in a critical point at pc = 170 ± 10 MPa and Tc = 181 ± 3 K. The TLL(p) of ML-BOP is almost identical to the one of TIP4P/2005, adding to the similarity in the equation of state of liquid water in both models. Cooling simulations reveal that ice crystallization is fastest at the LLT and its supercritical continuation of maximum heat capacity, supporting a mechanistic relationship between the structural transformation of water to a low-density liquid (LDL) and ice formation. We find no signature of liquid-liquid criticality in the ice crystallization temperatures. ML-BOP replicates the competition between formation of LDL and ice observed in ultrafast experiments of decompression of the high-density liquid (HDL) into the region of stability of LDL. The simulations reveal that crystallization occurs prior to the coarsening of the HDL and LDL domains, obscuring the distinction between the highly metastable first-order LLT and pronounced structural fluctuations along its supercritical continuation.

Potential energy landscape of a flexible water model: Equation of state, configurational entropy, and Adam-Gibbs relationship

The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid-liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam-Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully).

A water structure indicator suitable for generic contexts: Two-liquid behavior at hydration and nanoconfinement conditions and a molecular approach to hydrophobicity and wetting

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.

Tuning the low-temperature phase behavior of aqueous ionic liquids

Water's anomalous behavior is often explained using a two-liquid model, where two types of water, high-density liquid (HDL) and low-density liquid (LDL), can be separated via a liquid-liquid phase transition (LLPT) at low temperature. Mixtures of water and the ionic liquid hydrazinium trifluoroacetate were suggested to also show an LLPT but with the advantage that there is no rapid ice crystallization hampering its observation. It remains controversial whether these solutions exhibit an LLPT or are instead associated with complex phase separation phenomena. We here show detailed low-temperature calorimetry and diffraction experiments on aqueous solutions containing hydrazinium trifluoroacetate and other similar ionic liquids, all at a solute mole fraction of x = 0.175. Hydrazinium trifluoroacetate, ammonium trifluoroacetate, ethylammonium trifluoroacetate and hydrazinium pentafluoropropionate all boast exothermic transitions unrelated to crystallization as well as remarkable structural changes upon cooling into the glassy state. We propose a model inspired by micelle formation and decomposition in surfactant solutions, which is complemented by MD simulations and allows rationalizing the rich phase behavior of our mixtures during cooling. The fundamental aspect of the model is the hydrophobic nature of fluorinated anions that enables aggregation, which is reversed upon cooling and culminates in the remarkable exothermic first-order transition observed at low temperature. That is, we assign the first-order transition not to an LLPT but to phase-separations similar to the ones when falling below the Krafft temperature. All other solutions merely show simple vitrification behavior. Still, they exhibit distinct differences in liquid fragility, which is decreased continuously with decreasing hydrophobicity of the anions. This might enable the systematic tuning of ionic liquids with the goal of designing aqueous solutions of specific fragility.

Deciphering Density Fluctuations in the Hydration Water of Brownian Nanoparticles via Upconversion Thermometry

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.

Identification of local structures in water from supercooled to ambient conditions

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.

The structure of water: A historical perspective

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.

An entropy trap model of thermodynamic anomalies for dual-amorphous water undergoing liquid-liquid phase transition

Water displays numerous anomalously thermodynamic behaviors. However, the working principles behind these anomalies are not well understood, and the liquid-liquid phase transition (LLPT) is often regarded as the potential reason. In this study, we developed an entropy trap model to characterize the thermodynamic LLPT in dual-amorphous water, i.e. having both low-density and high-density liquid water. From the Adam-Gibbs model and free-volume theory, thermodynamic behaviors of water have been described using the proposed model, in which the constitutive relationships among density, heat capacity, thermal expansivity and glass transition temperature have been formulated. Moreover, the glass transition and its connection to thermodynamic behaviors were also investigated for dual-amorphous water. Finally, experimental data reported in the literature were used to verify effectiveness of the proposed model. This study is expected to provide a physical insight into the anomalous thermodynamics of dual-amorphous water undergoing the LLPT.

Intrinsic Dynamics of Amorphous Ice Revealed by a Heterodyne Signal in X-ray Photon Correlation Spectroscopy Experiments

Unraveling the mechanism of water's glass transition and the interconnection between amorphous ices and liquid water plays an important role in our overall understanding of water. X-ray photon correlation spectroscopy (XPCS) experiments were conducted to study the dynamics and the complex interplay between the hypothesized glass transition in high-density amorphous ice (HDA) and the subsequent transition to low-density amorphous ice (LDA). Our XPCS experiments demonstrate that a heterodyne signal appears in the correlation function. Such a signal is known to originate from the interplay of a static component and a dynamic component. Quantitative analysis was performed on this heterodyne signal to extract the intrinsic dynamics of amorphous ice during the HDA-LDA transition. An angular dependence indicates non-isotropic, heterogeneous dynamics in the sample. Using the Stokes-Einstein relation to extract diffusion coefficients, the data are consistent with the scenario of static LDA islands floating within a diffusive matrix of high-density liquid water.

Chilling alcohol on the computer: isothermal compressibility and the formation of hydrogen-bond clusters in liquid propan-1-ol

Molecular dynamics simulations have been performed to compute the isothermal compressibility [Formula: see text] of liquid propan-1-ol in the temperature range [Formula: see text] K. A change in behaviour, from normal (high T) to anomalous (low T), has been identified for [Formula: see text] at [Formula: see text] K. The average number of hydrogen bonds (H-bond) per molecule turns to saturation in the same temperature interval, suggesting the formation of a relatively rigid network. Indeed, simulation results show a strong tendency to form H-bond clusters with distinct boundaries, with the average largest size and width of the size distribution growing upon decreasing temperature, in agreement with previous theoretical and experimental studies. These results also emphasise a connection between the behaviour of [Formula: see text] and the formation of nanometric structures.

Surface premelting of ice far below the triple point

Premelting of ice, a quasi-liquid layer (QLL) at the surface below the melting temperature, was first postulated by Michael Faraday 160 y ago. Since then, it has been extensively studied theoretically and experimentally through many techniques. Existing work has been performed predominantly on hexagonal ice, at conditions close to the triple point. Whether the same phenomenon can persist at much lower pressure and temperature, where stacking disordered ice sublimates directly into water vapor, remains unclear. Herein, we report direct observations of surface premelting on ice nanocrystals below the sublimation temperature using transmission electron microscopy (TEM). Similar to what has been reported on hexagonal ice, a QLL is found at the solid-vapor interface. It preferentially decorates certain facets, and its thickness increases as the phase transition temperature is approached. In situ TEM reveals strong diffusion of the QLL, while electron energy loss spectroscopy confirms its amorphous nature. More significantly, the premelting observed in this work is thought to be related to the metastable low-density ultraviscous water, instead of ambient liquid water as in the case of hexagonal ice. This opens a route to understand premelting and grassy liquid state, far away from the normal water triple point.

Calibrating TDDFT Calculations of the X-ray Emission Spectrum of Liquid Water: The Effects of Hartree-Fock Exchange

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.

Anisotropy in spinodal-like dynamics of unknown water at ice V-water interface

Experimentally demonstrating the existence of waters with local structures unlike that of common water is critical for understanding both the origin of the mysterious properties of water and liquid polymorphism in single component liquids. At the interfaces between water and ices Ih, III, and VI grown/melted under pressure, we previously discovered low- and high-density unknown waters, that are immiscible with the surrounding water. Here, we show, by in-situ optical microscopy, that an unknown water appears at the ice V-water interface via spinodal-like dynamics. The dewetting dynamics of the unknown water indicate that its characteristic velocity is ~ 90 m/s. The time evolution of the characteristic length of the spinodal-like undulation suggests that the dynamics may be described by a common model for spinodal decomposition of an immiscible liquid mixture. Spinodal-like dewetting dynamics of the unknown water transiently showed anisotropy, implying the property of a liquid crystal.

RNA structures and dynamics with Å resolution revealed by x-ray free-electron lasers

RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free-electron laser sources to reveal the formation and ready identification of angstrom-scale features in structured and unstructured RNAs. Previously unrecognized structural signatures of RNA secondary and tertiary structures are identified through wide-angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base-paired intermediate to assume a triple-helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. This method may help to rapidly characterize and identify structural elements in nucleic acids in both equilibrium and time-resolved experiments.

Supercooled solutions of sodium perchlorate in TIP4P/2005 water: The effect of martian solutes on thermodynamics and structure

We study the thermodynamic behavior of sodium perchlorate solutions in supercooled water through molecular dynamics numerical simulations. These solutions are of special interest because of the recent experimental results that led to hypothesize the presence of liquid water in perchlorate solutions beneath the Martian soil. We model water using the TIP4P/2005 potential. The results we obtain for solutions with concentrations 1.63 and 15.4 wt% are in agreement with those of a system undergoing a liquid-liquid phase transition where the liquid-liquid critical point shifts to slightly higher temperatures and lower pressures. The structure of the system is also analyzed, and we come to the conclusion that, even at the highest concentration considered, water retains its anomalous behavior.

Strengthening Aqueous Electrolytes without Strengthening Water

Aqueous electrolytes typically suffer from poor electrochemical stability; however, eutectic aqueous solutions-25 wt.% LiCl and 62 wt.% H3 PO4 -cooled to -78 °C exhibit a significantly widened stability window. Integrated experimental and simulation results reveal that, upon cooling, Li+ ions become less hydrated and pair up with Cl- , ice-like water clusters form, and H⋅⋅⋅Cl- bonding strengthens. Surprisingly, this low-temperature solvation structure does not strengthen water molecules' O-H bond, bucking the conventional wisdom that increasing water's stability requires stiffening the O-H covalent bond. We propose a more general mechanism for water's low temperature inertness in the electrolyte: less favorable solvation of OH- and H+ , the byproducts of hydrogen and oxygen evolution reactions. To showcase this stability, we demonstrate an aqueous Li-ion battery using LiMn2 O4 cathode and CuSe anode with a high energy density of 109 Wh/kg. These results highlight the potential of aqueous batteries for polar and extraterrestrial missions.

Turning an energy-based defect detector into a multi-molecule structural indicator for water

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.

Realistic phase diagram of water from "first principles" data-driven quantum simulations

Since the experimental characterization of the low-pressure region of water's phase diagram in the early 1900s, scientists have been on a quest to understand the thermodynamic stability of ice polymorphs on the molecular level. In this study, we demonstrate that combining the MB-pol data-driven many-body potential for water, which was rigorously derived from "first principles" and exhibits chemical accuracy, with advanced enhanced-sampling algorithms, which correctly describe the quantum nature of molecular motion and thermodynamic equilibria, enables computer simulations of water's phase diagram with an unprecedented level of realism. Besides providing fundamental insights into how enthalpic, entropic, and nuclear quantum effects shape the free-energy landscape of water, we demonstrate that recent progress in "first principles" data-driven simulations, which rigorously encode many-body molecular interactions, has opened the door to realistic computational studies of complex molecular systems, bridging the gap between experiments and simulations.

RNA structures and dynamics with Å resolution revealed by x-ray free electron lasers

RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free electron laser sources to reveal the formation and ready identification of Å scale features in structured and unstructured RNAs. New structural signatures of RNA secondary and tertiary structures are identified through wide angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base paired intermediate to assume a triple helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. In addition to understanding how RNA triplexes form and thereby function as dynamic signaling elements, this new method can vastly increase the rate of structure determination for these biologically essential, but mostly uncharacterized macromolecules.

Coherent X-ray Scattering Reveals Nanoscale Fluctuations in Hydrated Proteins

Hydrated proteins undergo a transition in the deeply supercooled regime, which is attributed to rapid changes in hydration water and protein structural dynamics. Here, we investigate the nanoscale stress-relaxation in hydrated lysozyme proteins stimulated and probed by X-ray Photon Correlation Spectroscopy (XPCS). This approach allows us to access the nanoscale dynamics in the deeply supercooled regime (T = 180 K), which is typically not accessible through equilibrium methods. The observed stimulated dynamic response is attributed to collective stress-relaxation as the system transitions from a jammed granular state to an elastically driven regime. The relaxation time constants exhibit Arrhenius temperature dependence upon cooling with a minimum in the Kohlrausch-Williams-Watts exponent at T = 227 K. The observed minimum is attributed to an increase in dynamical heterogeneity, which coincides with enhanced fluctuations observed in the two-time correlation functions and a maximum in the dynamic susceptibility quantified by the normalized variance χ_T. The amplification of fluctuations is consistent with previous studies of hydrated proteins, which indicate the key role of density and enthalpy fluctuations in hydration water. Our study provides new insights into X-ray stimulated stress-relaxation and the underlying mechanisms behind spatiotemporal fluctuations in biological granular materials.

Electron diffraction of deeply supercooled water in no man's land

A generally accepted understanding of the anomalous properties of water will only emerge if it becomes possible to systematically characterize water in the deeply supercooled regime, from where the anomalies appear to emanate. This has largely remained elusive because water crystallizes rapidly between 160 K and 232 K. Here, we present an experimental approach to rapidly prepare deeply supercooled water at a well-defined temperature and probe it with electron diffraction before crystallization occurs. We show that as water is cooled from room temperature to cryogenic temperature, its structure evolves smoothly, approaching that of amorphous ice just below 200 K. Our experiments narrow down the range of possible explanations for the origin of the water anomalies and open up new avenues for studying supercooled water.

Phase Diagram of Aqueous Solutions of LiCl: a Study of Concentration Effects on the Anomalies of Water

We perform molecular dynamics simulations in order to study thermodynamics and the structure of supercooled aqueous solutions of lithium chloride (LiCl) at concentrations c = 0.678 and 2.034 mol/kg. We model the solvent using the TIP4P/2005 potential and the ions using the Madrid-2019 force field, a force field particularly suited for studying this solution. We find that, for c = 0.678 mol/kg, the behavior of the equation of state, studied in the P-T plane, indicates the presence of a liquid-liquid phase transition, similar to what was previously found for bulk water. We estimate the position of the liquid-liquid critical point to be at T_c ≈ 174 K, P_c ≈ 1775 bar, and ρ_c ≈ 1.065 g/cm³. When the concentration is tripled to c = 2.034 mol/kg, no critical point is observed, indicating its possible disappearance at this concentration. We also study the water-water and water-ions structure in the two solutions, and we find that at the concentrations examined the effect of ions on the water-water structure is not strong, and all the features found in bulk water are preserved. We also calculate the hydration number of the Li and Cl ions, and in line with experiments, we find the value of 4 for Li⁺ and between 5.5 and 6 for Cl^-, confirming the good performances of the Madrid-2019 force field.

New insights on the catalytic center of proteins from peptidylprolyl isomerase group based on the FOD-M model

Generating the structure of the hydrophobic core is based on the orientation of hydrophobic residues towards the central part of the protein molecule with the simultaneous exposure of polar residues. Such a course of the protein folding process takes place with the active participation of the polar water environment. While the self-assembly process leading to the formation of micelles concerns freely moving bi-polar molecules, bipolar amino acids in polypeptide chain have limited mobility due to the covalent bonds. Therefore, proteins form a more or less perfect micelle-like structure. The criterion is the hydrophobicity distribution, which to a greater or lesser extent reproduces the distribution expressed by the 3D Gaussian function on the protein body. The vast majority of proteins must ensure solubility, so a certain part of it-as it is expected-should reproduce the structuring of micelles. The biological activity of proteins is encoded in the part that does not reproduce the micelle-like system. The location and quantitative assessment of the contribution of orderliness to disorder is of critical importance for the determination of biological activity. The form of maladjustment to the 3D Gauss function may be varied-hence the obtained high diversity of specific interactions with strictly defined molecules: ligands or substrates. The correctness of this interpretation was verified on the basis of the group of enzymes Peptidylprolyl isomerase-E.C.5.2.1.8. In proteins representing this class of enzymes, zones responsible for solubility-micelle-like hydrophobicity system-the location and specificity of the incompatible part in which the specific activity of the enzyme is located and coded were identified. The present study showed that the enzymes of the discussed group show two different schemes of the structure of catalytic center (taking into account the status as defined by the fuzzy oil drop model).

Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

Liquid-Liquid Criticality in TIP4P/2005 and Three-State Models of Water

Molecular dynamics simulations leading to the isothermal compressibility, the isobaric thermal expansivity, and the isobaric heat capacity of TIP4P/2005 water are found to be consistent with the coordinates of its second, liquid-liquid critical point reported recently by Debenedetti et al. [ Science 2020, 369, 289-292]. In accord with the theory of critical phenomena, we encounter that the rise in the magnitude of these response functions as temperature is lowered is especially marked along the critical isochore. Furthermore, response-function ratios provide a test for thermodynamic consistency at the critical point and manifest nonuniversal features sharply distinguishing liquid-liquid from standard gas-liquid criticality. The whole pattern of behavior revealed by simulations is qualitatively the same as the one of a three-state Ising model of water exhibiting a low-temperature liquid-liquid critical point. Exact solutions for the two-state components of such a three-state model are also provided.

Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born-Oppenheimer potential energy surface. In particular, we focus on ab initio-based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through (a) parameterization of distributed multipoles, (b) explicit fitting of the MBE, (c) inclusion of many-atom features in a neural network, and (d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

Kinetics and Mechanisms of Pressure-Induced Ice Amorphization and Polyamorphic Transitions in a Machine-Learned Coarse-Grained Water Model

Water glasses have attracted considerable attention due to their potential connection to a liquid-liquid transition in supercooled water. Here we use molecular simulations to investigate the formation and phase behavior of water glasses using the machine-learned bond-order parameter (ML-BOP) water model. We produce glasses through hyperquenching of water, pressure-induced amorphization (PIA) of ice, and pressure-induced polyamorphic transformations. We find that PIA of polycrystalline ice occurs at a lower pressure than that of monocrystalline ice and through a different mechanism. The temperature dependence of the amorphization pressure of polycrystalline ice for ML-BOP agrees with that in experiments. We also find that ML-BOP accurately reproduces the density, coordination number, and structural features of low-density (LDA), high-density (HDA), and very high-density (VHDA) amorphous water glasses. ML-BOP accurately reproduces the experimental radial distribution function of LDA but overpredicts the minimum between the first two shells in high-density glasses. We examine the kinetics and mechanism of the transformation between low-density and high-density glasses and find that the sharp nature of these transitions in ML-BOP is similar to that in experiments and all-atom water models with a liquid-liquid transition. Transitions between ML-BOP glasses occur through a spinodal-like mechanism, similar to ice crystallization from LDA. Both glass-to-glass and glass-to-ice transformations have Avrami-Kolmogorov kinetics with exponent n = 1.5 ± 0.2 in experiments and simulations. Importantly, ML-BOP reproduces the competition between crystallization and HDA→LDA transition above the glass transition temperature T_g, and separation of their time scales below T_g, observed also in experiments. These findings demonstrate the ability of ML-BOP to accurately reproduce water properties across various regimes, making it a promising model for addressing the competition between polyamorphic transitions and crystallization in water and solutions.

A "short blanket" dilemma for a state-of-the-art neural network potential for water: Reproducing experimental properties or the physics of the underlying many-body interactions?

Deep neural network (DNN) potentials have recently gained popularity in computer simulations of a wide range of molecular systems, from liquids to materials. In this study, we explore the possibility of combining the computational efficiency of the DeePMD framework and the demonstrated accuracy of the MB-pol data-driven, many-body potential to train a DNN potential for large-scale simulations of water across its phase diagram. We find that the DNN potential is able to reliably reproduce the MB-pol results for liquid water, but provides a less accurate description of the vapor-liquid equilibrium properties. This shortcoming is traced back to the inability of the DNN potential to correctly represent many-body interactions. An attempt to explicitly include information about many-body effects results in a new DNN potential that exhibits the opposite performance, being able to correctly reproduce the MB-pol vapor-liquid equilibrium properties, but losing accuracy in the description of the liquid properties. These results suggest that DeePMD-based DNN potentials are not able to correctly "learn" and, consequently, represent many-body interactions, which implies that DNN potentials may have limited ability to predict the properties for state points that are not explicitly included in the training process. The computational efficiency of the DeePMD framework can still be exploited to train DNN potentials on data-driven many-body potentials, which can thus enable large-scale, "chemically accurate" simulations of various molecular systems, with the caveat that the target state points must have been adequately sampled by the reference data-driven many-body potential in order to guarantee a faithful representation of the associated properties.

Equation of State of Liquid Water Written by Simple Experimental Polynomials and the Liquid-Liquid Critical Point

Experimental and theoretical evidence has been accumulating to support the liquid-liquid critical point (LLCP) hypothesis for water. However, no agreement has yet been reached on the pressure and temperature of LLCP. Here we made simple experimental equations that reproduced the volume of liquid water measured over a wide pressure-temperature range. They were polynomials that calculate pressure using volume and temperature as variables, and coefficients were determined by the method of least-squares. We analyzed the polynomials by changing the volume data and the number of terms in the polynomials and extrapolated them slightly to low temperatures. Consequently, the available experimental volume of liquid water indicated (but did not prove) the existence of LLCP at low temperature. Representative polynomials suggested that LLCP locates around 105 ± ∼ 9 MPa, 207 ± ∼ 5 K, and 0.993 ± ∼ 0.009 cm³/g.

Liquid Water: A Single Approach to Its Two Continuous Phase Transitions

In this paper, we show that both continuous phase transitions of liquid water, the liquid-gas and the liquid-liquid, can be articulated within a single thermodynamic analytical formalism. This result follows from a combination of the two-liquid model (TLM), recently confirmed for water, with the idea of a thermal-dependent excluded volume, v_e, concept introduced by van der Waals, in his famous state equation. Starting from the fundamentals of thermodynamics, it will be shown that the TLM naturally leads to the idea of an extensive thermal-dependent v_e that acts as a parameter of the sample thermodynamic potentials. This procedure effectively separates the thermodynamics of the system into two parts: the first concerns the clusters' thermodynamics, taken as wandering particles, and the second concerns the thermal behavior of its internal structure (geometry and number of particles). From this result, we demonstrate that the condition of mechanical instability leads to not one but two critical points, each happening in one of the above-described parts of the system.

Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice

Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale < 100 ns and widely separated from ice formation, which begins at times >1 μs.

Single and multi-pulse based X-ray photon correlation spectroscopy

The ability of pulsed nature of synchrotron radiation opens up the possibility of studying microsecond dynamics in complex materials via speckle-based techniques. Here, we present the study of measuring the dynamics of a colloidal system by combining single and multiple X-ray pulses of a storage ring. In addition, we apply speckle correlation techniques at various pulse patterns to collect correlation functions from nanoseconds to milliseconds. The obtained sample dynamics from all correlation techniques at different pulse patterns are in very good agreement with the expected dynamics of Brownian motions of silica nanoparticles in water. Our study will pave the way for future pulsed X-ray investigations at various synchrotron X-ray sources using individual X-ray pulse patterns.

Refractive Index of Supercooled Water Down to 230.3 K in the Wavelength Range between 534 and 675 nm

Knowledge of the refractive index of water in the deeply supercooled metastable liquid state is important, for example, for an accurate description of optical reflection and refraction processes occurring in clouds. However, a measurement of both the temperature and wavelength dependence of the refractive index under such extreme conditions is challenging. Here, we employ Raman spectroscopy in combination with microscopic water jets in vacuum to obtain the refractive index of supercooled water to a lowest temperature of 230.3 K. While our approach is based on the analysis of Mie resonances in Raman spectra measured by using a single excitation wavelength at 532 nm, it allows us to obtain the refractive index in a wide visible wavelength range from 534 to 675 nm. Because of a direct link between the refractive index and density of water, our results provide a promising approach to help improve our understanding of water's anomalous behavior.

Correlated Fluctuations of Structural Indicators Close to the Liquid-Liquid Transition in Supercooled Water

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.

Liquid-Liquid Transition in Water from First Principles

A long-standing question in water research is the possibility that supercooled liquid water can undergo a liquid-liquid phase transition (LLT) into high- and low-density liquids. We used several complementary molecular simulation techniques to evaluate the possibility of an LLT in an ab initio neural network model of water trained on density functional theory calculations with the SCAN exchange correlation functional. We conclusively show the existence of a first-order LLT and an associated critical point in the SCAN description of water, representing the first definitive computational evidence for an LLT in water from first principles.

Stability and Metastability of Liquid Water in a Machine-Learned Coarse-Grained Model with Short-Range Interactions

Coarse-grained water models are ∼100 times more efficient than all-atom models, enabling simulations of supercooled water and crystallization. The machine-learned monatomic model ML-BOP reproduces the experimental equation of state (EOS) and ice-liquid thermodynamics at 0.1 MPa on par with the all-atom TIP4P/2005 and TIP4P/Ice models. These all-atom models were parametrized using high-pressure experimental data and are either accurate for water's EOS (TIP4P/2005) or ice-liquid equilibrium (TIP4P/Ice). ML-BOP was parametrized from temperature-dependent ice and liquid experimental densities and melting data at 0.1 MPa; its only pressure training is from compression of TIP4P/2005 ice at 0 K. Here we investigate whether ML-BOP replicates the experimental EOS and ice-water thermodynamics along all pressures of ice I. We find that ML-BOP reproduces the temperature, enthalpy, entropy, and volume of melting of hexagonal ice up to 400 MPa and the EOS of water along the melting line with an accuracy that rivals that of both TIP4P/2005 and TIP4P/Ice. We interpret that the accuracy of ML-BOP originates from its ability to capture the shift between compact and open local structures to changes in pressure and temperature. ML-BOP reproduces the sharpening of the tetrahedral peak of the pair distribution function of water upon supercooling, and its pressure dependence. We characterize the region of metastability of liquid ML-BOP with respect to crystallization and cavitation. The accessibility of ice crystallization to simulations of ML-BOP, together with its accurate representation of the thermodynamics of water, makes it promising for investigating the interplay between anomalies, glass transition, and crystallization under conditions challenging to access through experiments.

Using coherent X-rays to follow dynamics in amorphous ices

Amorphous solid water plays an important role in our overall understanding of water's phase diagram. X-ray scattering is an important tool for characterising the different states of water, and modern storage ring and XFEL facilities have opened up new pathways to simultaneously study structure and dynamics. Here, X-ray photon correlation spectroscopy (XPCS) was used to study the dynamics of high-density amorphous (HDA) ice upon heating. We follow the structural transition from HDA to low-density amorphous (LDA) ice, by using wide-angle X-ray scattering (WAXS), for different heating rates. We used a new type of sample preparation, which allowed us to study μm-sized ice layers rather than powdered bulk samples. The study focuses on the non-equilibrium dynamics during fast heating, spontaneous transformation and crystallization. Performing the XPCS study at ultra-small angle (USAXS) geometry allows us to characterize the transition dynamics at length scales ranging from 60 nm-800 nm. For the HDA-LDA transition we observe a clear separation in three dynamical regimes, which show different dynamical crossovers at different length scales. The crystallization from LDA, instead, is observed to appear homogenously throughout the studied length scales.

Multiple insights call for revision of modern thermodynamic models to account for structural fluctuations in water

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.

Revealing the three-component structure of water with principal component analysis (PCA) of X-ray spectra

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.

Quasi-water layer sandwiched between hexagonal ice and wall and its influences on the ice tensile stress

The presence of a quasi-water/premelting layer at the interface between wall and ice when the temperature (T) is below the melting point was extensively observed in experiments. In this work, molecular dynamics simulations are performed to analyze the underlying physics of the quasi-water layer and the effects of the layer on the ice tensile stress. The results indicate that each molecule and its four nearest neighbours in the quasi-water layer representing an equilibrium structure gradually form a tetrahedral ice-like arrangement from an unstructured liquid-like structure along the direction away from the wall. The average density of the quasi-water layer is equal to or higher than the bulk density of water at T ≥ 240 K or T ≤ 240 K respectively, and reaches 1.155 g cm^-3 at T = 210 K, suggesting a structural correlation with the high-density liquid phase of water. Depending on the temperature and wall wettability, the thickness of the quasi-water layer (H_q) ranges from ∼2 Å to ∼25 Å. For prescribed hydrophilic walls, H_q increases monotonically with temperature, and is almost proportional to(T_m - T)^-1/3, where T_m is the melting temperature of ice. H_q keeps an almost constant value (2 Å) as the temperature increases and rises sharply after passing a threshold temperature of T ≈ 250 K. In the joint effects of the wall wettability and quasi-water layer's thickness, the ice tensile stress decreasing monotonically at a larger temperature shows an upward trend and then reduces to almost a constant value as the wall changes from a hydrophobic to a hydrophilic one. The results reveal the potential development of anti-icing/de-icing techniques by heating the wall or modifying its surface to increase H_q.

Water's Unusual Thermodynamics in the Realm of Physical Chemistry

While it is known since the early work by Edsall, Frank and Evans, Kauzmann, and others that the thermodynamics of solvation of nonpolar solutes in water is unusual and has implications for the thermodynamics of protein folding, only recently have its connections with the unusual temperature dependence of the density of solvent water been illuminated. Such density behavior is, in turn, one of the manifestations of a nonstandard thermodynamic pattern contemplating a second, liquid-liquid critical point at conditions of temperature and pressure at which water exists as a deeply supercooled liquid. Recent experimental and computational work unambiguously points toward the existence of such a critical point, thereby providing concrete answers to the questions posed by the 1976 pioneering experiments by Speedy and Angell and the associated "liquid-liquid transition hypothesis" posited in 1992 by Stanley and co-workers. Challenges of this phenomenology to the branch of Statistical Mechanics remain.