Supercooled water and the kinetic glass transition

Evidence of ferroelectric features in low-density supercooled water from ab initio deep neural-network simulations

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.

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl2 from Tm + 85 K to Tm + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl2. We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. This result may contribute to new understandings of the structural relaxation of other network liquids.

Mode coupling behavior and fragile to strong transition of trehalose in a binary mixture with water upon supercooling

We perform molecular dynamics simulations of a binary mixture of water and trehalose with the TIP4P/Ice water model. We analyze the slow dynamics of trehalose molecules in the mildly supercooled region for concentrations of 3.66 and 18.57 wt. %. We previously studied the dynamics of water in the same mixtures. Supercooled TIP4P/Ice water solvating trehalose molecules was found to follow the Mode Coupling Theory (MCT) and to undergo a transition from a fragile to a strong behavior for both concentrations. Here, we show that also the dynamics of trehalose molecules follows the MCT and displays a fragile to strong crossover (FSC). The results show that trehalose in binary mixtures with water shares with it the dynamical behavior typical of glass forming liquids. Moreover, the FSC for trehalose structural relaxation times is found to occur at temperatures close to those previously obtained for water in the same solutions, showing that the dynamics of the solute is strongly coupled to that of the solvent. We also perform a MCT test showing that the trehalose dynamics obeys the MCT time-temperature superposition principle and that the exponents derived from the theory and the ones obtained from fitting procedure of the relaxation times are comparable, confirming that trehalose molecules in supercooled water solutions follow the MCT of glassy dynamics. Moreover, as predicted by the theory, trehalose particles have MCT parameters comparable to those of water in the same mixtures. This is an important result, given that MCT was originally formulated for monoatomic particles.

Structure and slow dynamics of protein hydration water with cryopreserving DMSO and trehalose upon cooling

We study, through molecular dynamics simulations, three aqueous solutions with one lysozyme protein and three different concentrations of trehalose and dimethyl sulfoxide (DMSO). We analyze the structural and dynamical properties of the protein hydration water upon cooling. We find that trehalose plays a major role in modifying the structure of the network of HBs between water molecules in the hydration layer of the protein. The dynamics of hydration water presents, in addition to the α-relaxation, typical of glass formers, a slower long-time relaxation process, which greatly slows down the dynamics of water, particularly in the systems with trehalose, where it becomes dominant at low temperatures. In all the solutions, we observe, from the behavior of the α-relaxation times, a shift of the Mode Coupling Theory crossover temperature and the fragile-to-strong crossover temperature toward higher values with respect to bulk water. We also observe a strong-to-strong crossover from the temperature behavior of the long-relaxation times. In the aqueous solution with only DMSO, the transition shifts to a lower temperature than in the case with only lysozyme reported in the literature. We observe that the addition of trehalose to the mixture has the opposite effect of restoring the original location of the strong-to-strong crossover. In all the solutions analyzed in this work, the observed temperature of the protein dynamical transition is slightly shifted at lower temperatures than that of the strong-to-strong crossover, but their relative order is the same, showing a correlation between the motion of the protein and that of the hydration water.

Water slowing down drives the occurrence of the low temperature dynamical transition in microgels

The protein dynamical transition marks an increase in atomic mobility and the onset of anharmonic motions at a critical temperature (T d), which is considered relevant for protein functionality. This phenomenon is ubiquitous, regardless of protein composition, structure and biological function and typically occurs at large protein content, to avoid water crystallization. Recently, a dynamical transition has also been reported in non-biological macromolecules, such as poly(N-isopropyl acrylamide) (PNIPAM) microgels, bearing many similarities to proteins. While the generality of this phenomenon is well-established, the role of water in the transition remains a subject of debate. In this study, we use atomistic molecular dynamics (MD) simulations and elastic incoherent neutron scattering (EINS) experiments with selective deuteration to investigate the microscopic origin of the dynamical transition and distinguish water and PNIPAM roles. While a standard analysis of EINS experiments would suggest that the dynamical transition occurs in PNIPAM and water at a similar temperature, simulations reveal a different perspective, also qualitatively supported by experiments. From room temperature down to about 180 K, PNIPAM exhibits only modest changes of dynamics, while water, being mainly hydration water under the probed extreme confinement, significantly slows down and undergoes a mode-coupling transition from diffusive to activated. Our findings therefore challenge the traditional view of the dynamical transition, demonstrating that it occurs in proximity of the water mode-coupling transition, shedding light on the intricate interplay between polymer and water dynamics.

Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum

In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

Multilayered noise model for transport in complex environments

Transport in complex fluidic environments often exhibits transient subdiffusive dynamics accompanied by non-Gaussian probability density profiles featuring a nonmonotonic non-Gaussian parameter. Such properties cannot be adequately explained by the original theory of Brownian motion. Based on an extension of kinetic theory, this study introduces a chain of hierarchically coupled random walks approach that effectively captures all these intriguing characteristics. If the environment consists of a series of independent white noise sources, then the problem can be expressed as a system of hierarchically coupled Ornstein-Uhlenbech equations. Due to the linearity of the system, the most essential transport properties have a closed analytical form.

Glassy dynamics of water in TIP4P/Ice aqueous solutions of trehalose in comparison with the bulk phase

We perform molecular dynamics simulations of TIP4P/Ice water in solution with trehalose for 3.65 and 18.57 wt. % concentrations and of bulk TIP4P/Ice water at ambient pressure, to characterize the structure and dynamics of water in a sugar aqueous solution in the supercooled region. We find here that TIP4P/Ice water in solution with trehalose molecules follows the Mode Coupling Theory and undergoes a fragile to strong transition up to the highest concentration investigated, similar to the bulk. Moreover, we perform a Mode Coupling Theory test, showing that the Time Temperature Superposition principle holds for both bulk TIP4P/Ice water and for TIP4P/Ice water in the solutions and we calculate the exponents of the theory. The direct comparison of the dynamical results for bulk water and water in the solutions shows upon cooling along the isobar a fastening of water dynamics for lower temperatures, T < 240 K. We found that the counter-intuitive behavior for the low temperature solutions can be explained with the diffusion anomaly of water leading us to the conclusion that the fastening observed below T = 240 K in water dynamics is only fictitious, due to the fact that the density of water molecules in the solutions is higher than the density of the bulk at the same temperature and pressure. This result should be taken into account in experimental investigations which are often carried out at constant pressure.

The influence of force fields on the structure and dynamics of water confined in ZIF-8 from atomistic simulations

The complexity of modeling flexible crystals, such as ZIF-8, mainly stems from the handling of intramolecular interactions. Numerous force fields have been proposed in the literature to describe the interactions between atoms in ZIF-8. We employ seven force fields to examine the structure and dynamic behavior of water molecules confined in ZIF-8, with the aim of investigating the impact of force fields on simulation results. Various structural characterization methods consistently indicate that the choice of different force fields has quantitative effects but no qualitative effects on the structural characteristics of confined water. Additionally, the force fields do not impact the qualitative description of the diffusion mechanism. Both mean-square displacement and van Hove autocorrelation function reveal two characteristic movements of water molecules diffusing in ZIF-8: a short-time intra-cavity hopping process and a long-time inter-cavity hopping process. However, the framework flexibility is found to play a crucial role in determining the order of spatial arrangement and local structure, self-diffusion coefficient and reorientational dynamics of confined water. Specifically, the DREIDING force field gives rise to an unrealistic stiff framework, enhancing the order of spatial arrangement and diminishing the local ordered structure of confined water. Meanwhile, it results in much slower translational and reorientational dynamics. Hence, the general DREIDING force field cannot be considered for providing a quantitative description of the water structure and dynamics.

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.

Molecular Dynamics Simulation of Droplet Impact on a Hydrophobic 3D Elastic Surface

The phenomenon of droplets impacting elastic surfaces is common in nature and in many engineering applications. It has been shown that droplet impact on an elastic surface drastically reduces droplet contact time and hinders droplet spreading. However, most of the current studies are based on experiments, and the analysis of the influence mechanism of the elastic substrate on the dynamic behavior of droplets is not complete. In addition, the simulations of droplet impact on elastic substrates are mainly focused on 2D elastic films or vibrating rigid substrates, ignoring the effect of 3D elastic substrate deformation on the droplet dynamic behavior. Therefore, in this paper, we propose to model the droplet impact on a 3D hydrophobic elastic substrate using the molecular dynamics method. We find that droplet pancake rebound can substantially reduce the droplet contact time. Moreover, we record the conditions required for the pancake rebound of the droplet. Furthermore, we investigated the effects of the elastic modulus of the substrate and the initial velocity of the droplet on the droplet contact time, contact area, and spreading factor. This study further elucidates the influence mechanism of the elastic substrate on the dynamic behavior of the droplet and provides theoretical guidance for regulating the dynamic behavior of the droplet in related fields.

Revealing the hidden dynamics of confined water in acrylate polymers: Insights from hydrogen-bond lifetime analysis

Polymers contain functional groups that participate in hydrogen bond (H-bond) with water molecules, establishing a robust H-bond network that influences bulk properties. This study utilized molecular dynamics (MD) simulations to examine the H-bonding dynamics of water molecules confined within three poly(meth)acrylates: poly(2-methoxyethyl acrylate) (PMEA), poly(2-hydroxyethyl methacrylate) (PHEMA), and poly(1-methoxymethyl acrylate) (PMC1A). Results showed that H-bonding dynamics significantly slowed as the water content decreased. Additionally, the diffusion of water molecules and its correlation with H-bond breakage were analyzed. Our findings suggest that when the H-bonds between water molecules and the methoxy oxygen of PMEA are disrupted, those water molecules persist in close proximity and do not diffuse on a picosecond time scale. In contrast, the water molecules H-bonded with the hydroxy oxygen of PHEMA and the methoxy oxygen of PMC1A diffuse concomitantly with the breakage of H-bonds. These results provide an in-depth understanding of the impact of polymer functional groups on H-bonding dynamics.

Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.

Correlations between defect propensity and dynamical heterogeneities in supercooled water

A salient feature of supercooled liquids consists in the dramatic dynamical slowdown they undergo as temperature decreases while no significant structural change is evident. These systems also present dynamical heterogeneities (DH): certain molecules, spatially arranged in clusters, relax various orders of magnitude faster than the others. However, again, no static quantity (such as structural or energetic measures) shows strong direct correlations with such fast-moving molecules. In turn, the dynamic propensity approach, an indirect measure that quantifies the tendency of the molecules to move in a given structural configuration, has revealed that dynamical constraints, indeed, originate from the initial structure. Nevertheless, this approach is not able to elicit which structural quantity is, in fact, responsible for such a behavior. In an effort to remove dynamics from its definition in favor of a static quantity, an energy-based propensity has also been developed for supercooled water, but it could only find positive correlations between the lowest-energy and the least-mobile molecules, while no correlations could be found for those more relevant mobile molecules involved in the DH clusters responsible for the system's structural relaxation. Thus, in this work, we shall define a defect propensity measure based on a recently introduced structural index that accurately characterizes water structural defects. We shall show that this defect propensity measure provides positive correlations with dynamic propensity, being also able to account for the fast-moving molecules responsible for the structural relaxation. Moreover, time dependent correlations will show that defect propensity represents an appropriate early-time predictor of the long-time dynamical heterogeneity.

Quantum uncertainty effects in the dynamics of supercooled liquids: A molecular dynamics study

Dynamics of density fluctuations in quantum supercooled liquids is analyzed using molecular dynamics simulations. In contrast to the classical case, the uncertainty in the particle position (delocalization of quantum particle in space) leads to significant differences in the dynamics of quantum liquids, both in the short- and long-time limits. The effect of uncertainty is found to be significant for length scales smaller than the uncertainty itself, and diminishes as the length scale grows. The dynamic heterogeneity of the system at short times is enhanced due to uncertainty. In the intermediate (β-relaxation) time regime, the heterogeneity tends to get suppressed due to quantum uncertainty. The probability distribution of particle displacements shows highly nonclassical behavior with double-peak structure at short timescales.

Water Diffusion Proceeds via a Hydrogen-Bond Jump Exchange Mechanism

The self-diffusion of water molecules plays a key part in a broad range of essential processes in biochemistry, medical imaging, material science, and engineering. However, its molecular mechanism and the role played by the water hydrogen-bond network rearrangements are not known. Here we combine molecular dynamics simulations and analytic modeling to determine the molecular mechanism of water diffusion. We establish a quantitative connection between the water diffusion coefficient and hydrogen-bond jump exchanges, and identify the features that determine the underlying energetic barrier. We thus provide a unified framework to understand the coupling between translational, rotational, and hydrogen-bond dynamics in liquid water. It explains why these different dynamics do not necessarily exhibit identical temperature dependences although they all result from the same hydrogen-bond exchange events. The consequences for the understanding of water diffusion in supercooled conditions and for water transport in complex aqueous systems, including ionic, biological, and confined solutions, are discussed.

Dynamical Manifestations of Supercooling in Amyloid Hydration

The effect of extreme temperature on amyloidogenic species remains sparsely explored. In a recent study (J. Phys. Chem. Lett., 2019, 10, (10)), we employed exhaustive molecular dynamics simulations to explore the cold thermal response of a putative small amyloid oligomer and to elicit the role of solvent modulation. Herein, we investigate the dynamical response of the hydration waters of the oligomer within the supercooled states. Using NMR-based formalism, we delineate the entropic response in terms of the side-chain conformational entropy that corroborates the weakening of the hydrophobic core with lowering of temperature. The translational dynamics of the protein and hydration waters reveal the coupling of protein dynamical fluctuations with solvent dynamics under supercooled conditions. Probing the translational motion as a space-time correlation indicates glassy dynamics exhibited by hydration waters in the supercooled regime. Caging of the water molecules with lowering of temperature and the resultant hopping dynamics are reflected in the longer β-relaxation timescales of translational motion. Furthermore, we utilized mode-coupling theory (MCT) and derived the ideal glass transition temperature from translational and rotational dynamics, around ∼196 and 209 K, respectively. Interestingly, rotational motion in the supercooled regime deviates from the MCT law, exhibits Arrhenius motion, and marks a fragile-to-strong crossover at 227 K. The low-frequency vibrational modes also coincide with the dynamical transition. This exposition lends dynamical insights into the hydration coupling of an amyloid aggregate under cryogenic conditions.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

Structure and dynamics of nanoconfined water and aqueous solutions

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Dynamical crossover and its connection to the Widom line in supercooled TIP4P/Ice water

We perform molecular dynamics simulations with the TIP4P/Ice water model to characterize the relationship between dynamics and thermodynamics of liquid water in the supercooled region. We calculate the relevant properties of the phase diagram, and we find that TIP4P/Ice presents a retracing line of density maxima, similar to what was previously found for atomistic water models and models of other tetrahedral liquids. For this model, a liquid-liquid critical point between a high-density liquid and a low-density liquid was recently found. We compute the lines of the maxima of isothermal compressibility and the minima of the coefficient of thermal expansion in the one phase region, and we show that these lines point to the liquid-liquid critical point while collapsing on the Widom line. This line is the line of the maxima of correlation length that emanates from a second order critical point in the one phase region. Supercooled water was found to follow mode coupling theory and to undergo a transition from a fragile to a strong behavior right at the crossing of the Widom line. We find here that this phenomenology also happens for TIP4P/Ice. Our results appear, therefore, to be a general characteristic of supercooled water, which does not depend on the interaction potential used, and they reinforce the idea that the dynamical crossover from a region where the relaxation mechanism is dominated by cage relaxation to a region where cages are frozen and hopping dominates is correlated in water to a phase transition between a high-density liquid and a low-density liquid.

Transition pathway of hydrogen bond switching in supercooled water analyzed by the Markov state model

In this work, we examine hydrogen-bond (H-bond) switching by employing the Markov State Model (MSM). During the H-bond switching, a water hydrogen initially H-bonded with water oxygen becomes H-bonded to a different water oxygen. MSM analysis was applied to trajectories generated from molecular dynamics simulations of the TIP4P/2005 model from a room-temperature state to a supercooled state. We defined four basis states to characterize the configuration between two water molecules: H-bonded ("H"), unbound ("U"), weakly H-bonded ("w"), and alternative H-bonded ("a") states. A 16 × 16 MSM matrix was constructed, describing the transition probability between states composed of three water molecules. The mean first-passage time of the H-bond switching was estimated by calculating the total flux from the HU to UH states. It is demonstrated that the temperature dependence of the mean first-passage time is in accordance with that of the H-bond lifetime determined from the H-bond correlation function. Furthermore, the flux for the H-bond switching is decomposed into individual pathways that are characterized by different forms of H-bond configurations of trimers. The dominant pathway of the H-bond switching is found to be a direct one without passing through such intermediate states as "w" and "a," the existence of which becomes evident in supercooled water. The pathway through "w" indicates a large reorientation of the donor molecule. In contrast, the pathway through "a" utilizes the tetrahedral H-bond network, which is revealed by the further decomposition based on the H-bond number of the acceptor molecule.

Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water

We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.

Water dynamics within nanostructured amphiphilic statistical copolymers from quasielastic neutron scattering

Understanding the properties of water under either soft or hard confinement has been an area of great interest, but nanostructured amphiphilic polymers that provide a secondary confinement have garnered significantly less attention. Here, a series of statistical copolymers of 2-hydroxyethyl acrylate (HEA) and 2-(N-ethylperfluorooctane sulfonamido)ethyl methacrylate (FOSM) are swollen to equilibrium in water to form nanostructured physically cross-linked hydrogels to probe the effect of soft confinement on the dynamics of water. Changing the composition of the copolymer from 10 to 21 mol. % FOSM decreases the average size of the assembled FOSM cross-link, but also the spacing between the cross-links in the hydrogels with the mean distance between the FOSM aggregates decreasing from 3.9 to 2.7 nm. The dynamics of water within the hydrogels were assessed with quasielastic neutron scattering. These hydrogels exhibit superior performance for inhibition of water crystallization on supercooling in comparison to analogous hydrogels with different hydrophilic copolymer chemistries. Despite the lower water crystallinity, the self-diffusion coefficient for these hydrogels from the copolymers of HEA and FOSM decreases precipitously below 260 K, which is a counter to the nearly temperature invariant water dynamics reported previously with an analogous hydrogel [Wiener et al., J. Phys. Chem. B 120, 5543 (2016)] that exhibits nearly temperature invariant dynamics to 220 K. These results point to chemistry dependent dynamics of water that is confined within amphiphilic hydrogels, where the interactions of water with the hydrophilic segments can qualitatively alter the temperature dependent dynamics of water in the supercooled state.

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

Insight into the Structure and Dynamics of Polymers by Neutron Scattering Combined with Atomistic Molecular Dynamics Simulations

Combining neutron scattering and fully atomistic molecular dynamics simulations allows unraveling structural and dynamical features of polymer melts at different length scales, mainly in the intermolecular and monomeric range. Here we present the methodology developed by us and the results of its application during the last years in a variety of polymers. This methodology is based on two pillars: (i) both techniques cover approximately the same length and time scales and (ii) the classical van Hove formalism allows easily calculating the magnitudes measured by neutron scattering from the simulated atomic trajectories. By direct comparison with experimental results, the simulated cell is validated. Thereafter, the information of the simulations can be exploited, calculating magnitudes that are experimentally inaccessible or extending the parameters range beyond the experimental capabilities. We show how detailed microscopic insight on structural features and dynamical processes of various kinds has been gained in polymeric systems with different degrees of complexity, and how intriguing questions as the collective behavior at intermediate length scales have been faced.

Protein-Water and Water-Water Long-Time Relaxations in Protein Hydration Water upon Cooling-A Close Look through Density Correlation Functions

We report results on the translational dynamics of the hydration water of the lysozyme protein upon cooling obtained by means of molecular dynamics simulations. The self van Hove functions and the mean square displacements of hydration water show two different temperature activated relaxation mechanisms, determining two dynamic regimes where transient trapping of the molecules is followed by hopping phenomena to allow to the structural relaxations. The two caging and hopping regimes are different in their nature. The low-temperature hopping regime has a time scale of tenths of nanoseconds and a length scale on the order of 2–3 water shells. This is connected to the nearest-neighbours cage effect and restricted to the supercooling, it is absent at high temperature and it is the mechanism to escape from the cage also present in bulk water. The second hopping regime is active at high temperatures, on the nanoseconds time scale and over distances of nanometers. This regime is connected to water displacements driven by the protein motion and it is observed very clearly at high temperatures and for temperatures higher than the protein dynamical transition. Below this temperature, the suppression of protein fluctuations largely increases the time-scale of the protein-related hopping phenomena at least over 100 ns. These protein-related hopping phenomena permit the detection of translational motions of hydration water molecules longly persistent in the hydration shell of the protein.

Dynamical properties of densely packed confined hard-sphere fluids

Numerical solutions of the mode-coupling theory (MCT) equations for a hard-sphere fluid confined between two parallel hard walls are elaborated. The governing equations feature multiple parallel relaxation channels which significantly complicate their numerical integration. We investigate the intermediate scattering functions and the susceptibility spectra close to structural arrest and compare to an asymptotic analysis of the MCT equations. We corroborate that the data converge in the β-scaling regime to two asymptotic power laws, viz. the critical decay and the von Schweidler law. The numerical results reveal a nonmonotonic dependence of the power-law exponents on the slab width and a nontrivial kink in the low-frequency susceptibility spectra. We also find qualitative agreement of these theoretical results to event-driven molecular dynamics simulations of polydisperse hard-sphere systems. In particular, the nontrivial dependence of the dynamical properties on the slab width is well reproduced.

Electron Paramagnetic Resonance Measurements of Four Nitroxide Probes in Supercooled Water Explained by Molecular Dynamics Simulations

Electron paramagnetic resonance (EPR) measurements of the rotational diffusion of small nitroxide probes have been demonstrated to be a powerful technique for experimentally investigating the properties of supercooled liquids, such as water. However, since only the rotational diffusion of the probe molecules is measured and EPR measurements are indirect, it is not clear what the relationship between the behavior of water and the probe molecule is. To address this, we have performed molecular dynamics simulations of four nitroxide probes in TIP4P-Ew and OPC water models to directly compare with EPR experiments and to determine the behavior of the water and the underlying microscopic coupling between the water and the probes. In all, 200 ns simulations were run for 23 temperatures between 253 and 283 K for all four probes with each water model for an aggregate of 36.8 μs of simulation time. Simulations for both water models systematically underestimated the rotational diffusion coefficients for both water and probes, though OPC simulations were generally in better agreement with the experiments than TIP4P-Ew simulations. Despite this, when the temperature dependence of the data was fit to a power law, fit parameters for TIP4P-Ew were generally in better agreement with the experiments than OPC. For probe molecules, the singular temperature was found to be T₀ = 226.5 ± 0.4 K from experiments, T₀ = 208 ± 2 K for OPC water, and T₀ = 215 ± 2 K for TIP4P-Ew water. While for water molecules, the singular temperature was found to be T₀ = 220.3 ± 0.2 K from experiments, T₀ = 208 ± 2 K for OPC water, and T₀ = 220 ± 1 K for TIP4P-Ew water. Systematic underestimation of the rotational diffusion coefficients was most pronounced at lower temperatures and was clearly observed in changes to the Arrhenius activation energy. Above the maximum density temperature of T_ρmax = 277 K, an activation energy of E_A ≈ 16.7 kJ/mol was observed for the probes from experiments, while OPC had E_A ≈ 15.2 kJ/mol and TIP4P-Ew had E_A ≈ 14.6 kJ/mol. Below the maximum density temperature, the activation energy jumped to E_A ≈ 32.5 kJ/mol for experiments but only E_A ≈ 23 kJ/mol for OPC and E_A ≈ 22 kJ/mol for TIP4P-Ew. In all cases, we saw good agreement between the behavior of the probe molecules and water. To understand why, we calculated the average number of hydrogen bonds between the probe molecules and water. From this, we were able to explain the rotational diffusion times for all of the probes. These results show that current molecular models are sufficient to capture physical phenomena observed with EPR and to help elucidate why the probes provide accurate insights into the behavior of supercooled water.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Structure, Dynamics, and Wettability of Water at Metal Interfaces

The water/metal interface often governs important chemophysical processes in various technologies. Therefore, from scientific and engineering perspectives, the detailed molecular-level elucidation of the water/metal interface is of high priority, but the related research is limited. In experiments, the surface-science techniques, which can provide full structural details of the surface, are not easy to directly apply to the interfacial systems under ambient conditions, and the well-defined facets cannot be entirely free from contamination at the contact with water. To answer long-standing debates regarding the wettability, structure, and dynamics of water at metal interfaces, we here develop reliable first-principles-based multiscale simulations. Using the state-of-the-art simulations, we find that the clean metal surfaces are actually superhydrophilic and yield zero contact angles. Furthermore, we disclose an inadequacy of widespread ice-like bilayer model of the water adlayers on metal surfaces from both averaged structural and dynamic points of view. Our findings on the nature of water on metal surfaces provide new molecular level perspectives on the tuning and design of water/metal interfaces that are at the heart of many energy applications.

Dynamics of TMAO and urea in the hydration shell of the protein SNase

Using all-atom molecular dynamics simulations of aqueous solutions of the globular protein SNase, the dynamic behavior of water molecules and cosolvents (trimethylamine-N-oxide (TMAO) and urea) in the hydration shell of the protein was studied for different solvent compositions. TMAO is a potent protein-stabilizing osmolyte, whereas urea is known to destabilize proteins. For molecules that are initially located in successive narrow layers at a given distance from the protein, the mean displacements and the distribution of displacements for short time intervals are calculated. For molecules that are initially located in solvation shells of a given thickness around the protein, the characteristic residence times in these shells are determined to characterize the dynamic behavior of the solvent molecules as a function of the distance to the protein. A combined consideration of these characteristics allows to reveal additional features of the dynamics of the cosolvents. It is shown that TMAO molecules leave the nearest vicinity of the protein faster than urea molecules, despite the fact that the mobility of TMAO molecules, measured by their mean displacements, is lower than that of urea. Moreover, we show that the rate of release of TMAO molecules from the hydration shell is lower in ternary (TMAO + urea + H₂O) solvent mixtures than in the binary ones. This is consistent with a recent observation that the fraction of TMAO near the protein decreases in the presence of urea. From the analysis of the decay of the number of particles initially located in the region of the first peak of the distribution function of solvent molecules around the protein, we estimated that about 20 water molecules and 6-7 urea molecules stay near the protein for more than 1000 ps.

Characterization of hydration water in supercooled water-trehalose solutions: The role of the hydrogen bonds network

The structural and dynamical properties of hydration water in aqueous solutions of trehalose are studied with molecular dynamics simulation. We simulate the systems in the supercooled region to investigate how the interaction with the trehalose molecules modifies the hydrogen bond network, the structural relaxation, and the diffusion properties of hydration water. The analysis is performed by considering the radial distribution functions, the residence time of water molecules in the hydration shell, the two body excess entropy, and the hydrogen bond water-water and water-trehalose correlations of the hydration water. The study of the two body excess entropy shows the presence of a fragile to strong crossover in supercooled hydration water also found in the relaxation time of the water-water hydrogen bond correlation function, and this is in agreement with predictions of the mode coupling theory and of previous studies of the oxygen-oxygen density correlators [A. Iorio et al., J. Mol. Liq. 282, 617 (2019); Sci. China: Phys., Mech. Astron. 62, 107011 (2019)]. The water-trehalose hydrogen bond correlation function instead evidences a strong to strong crossover in the relaxation time, and this crossover is related to a trehalose dynamical transition. This signals the role that the strong interplay between the soluted molecules and the surrounding solvent has in determining the dynamical transition common to both components of the system that happens upon cooling and that is similar to the well known protein dynamical transition. We connect our results with the cryoprotecting role of trehalose molecules.

Transport dynamics of complex fluids

Thermal motion in complex fluids is a complicated stochastic process but ubiquitously exhibits initial ballistic, intermediate subdiffusive, and long-time diffusive motion, unless interrupted. Despite its relevance to numerous dynamical processes of interest in modern science, a unified, quantitative understanding of thermal motion in complex fluids remains a challenging problem. Here, we present a transport equation and its solutions, which yield a unified quantitative explanation of the mean-square displacement (MSD), the non-Gaussian parameter (NGP), and the displacement distribution of complex fluids. In our approach, the environment-coupled diffusion kernel and its time correlation function (TCF) are the essential quantities that determine transport dynamics and characterize mobility fluctuation of complex fluids; their time profiles are directly extractable from a model-free analysis of the MSD and NGP or, with greater computational expense, from the two-point and four-point velocity autocorrelation functions. We construct a general, explicit model of the diffusion kernel, comprising one unbound-mode and multiple bound-mode components, which provides an excellent approximate description of transport dynamics of various complex fluidic systems such as supercooled water, colloidal beads diffusing on lipid tubes, and dense hard disk fluid. We also introduce the concepts of intrinsic disorder and extrinsic disorder that have distinct effects on transport dynamics and different dependencies on temperature and density. This work presents an unexplored direction for quantitative understanding of transport and transport-coupled processes in complex disordered media.

Translational and rotational dynamics of high and low density TIP4P/2005 water

We use molecular dynamics simulations using TIP4P/2005 to investigate the self- and distinct-van Hove functions for different local environments of water, classified using the local structure index as an order parameter. The orientational dynamics were studied through the calculation of the time-correlation functions of different-order Legendre polynomials in the OH-bond unit vector. We found that the translational and orientational dynamics are slower for molecules in a low-density local environment and correspondingly the mobility is enhanced upon increasing the local density, consistent with some previous works, but opposite to a recent study on the van Hove function. From the analysis of the distinct dynamics, we find that the second and fourth peaks of the radial distribution function, previously identified as low density-like arrangements, show long persistence in time. The analysis of the time-dependent interparticle distance between the central molecule and the first coordination shell shows that particle identity persists longer than distinct van Hove correlations. The motion of two first-nearest-neighbor molecules thus remains coupled even when this correlation function has been completely decayed. With respect to the orientational dynamics, we show that correlation functions of molecules in a low-density environment decay exponentially, while molecules in a local high-density environment exhibit bi-exponential decay, indicating that dynamic heterogeneity of water is associated with the heterogeneity among high-density and between high-density and low-density species. This bi-exponential behavior is associated with the existence of interstitial waters and the collapse of the second coordination sphere in high-density arrangements, but not with H-bond strength.

Water-Polymer Coupling Induces a Dynamical Transition in Microgels

The long debated protein dynamical transition was recently found also in nonbiological macromolecules, such as poly- N-isopropylacrylamide (PNIPAM) microgels. Here, by using atomistic molecular dynamics simulations, we report a description of the molecular origin of the dynamical transition in these systems. We show that PNIPAM and water dynamics below the dynamical transition temperature T _d are dominated by methyl group rotations and hydrogen bonding, respectively. By comparing with bulk water, we unambiguously identify PNIPAM-water hydrogen bonding as mainly responsible for the occurrence of the transition. The observed phenomenology thus crucially depends on the water-macromolecule coupling, being relevant to a wide class of hydrated systems, independently from the biological function.

Common microscopic structural origin for water's thermodynamic and dynamic anomalies

Water displays a vast array of unique properties, known as water's anomalies, whose origin remains subject to hot debate. Our aim in this article is to provide a unified microscopic physical picture of water's anomalies in terms of locally favored structures, encompassing both thermodynamic and dynamic anomalies, which are often attributed to different origins. We first identify locally favored structures via a microscopic structural descriptor that measures local translational order and provide direct evidence that they have a hierarchical impact on the anomalies. At each state point, the strength of thermodynamic anomalies is directly proportional to the amount of locally favored structures, while the dynamic properties of each molecule depend on the local structure surrounding both itself and its nearest neighbors. To incorporate this, we develop a novel hierarchical two-state model. We show by extensive simulations of two popular water models that both thermodynamic and kinetic anomalies can be almost perfectly explained by the temperature and pressure dependence of these local and non-local versions of the same structural descriptor, respectively. Moreover, our scenario makes three unique predictions in supercooled water, setting it apart from other scenarios: (1) Presence of an "Arrhenius-to-Arrhenius" crossover upon cooling, as the origin of the apparent "fragile-to-strong" transition; (2) maximum of dynamic heterogeneity around 20 K below the Widom line and far above the glass transition; (3) violation of the Stokes-Einstein-Debye relation at ∼2T _g, rather than 1.2T _g typical of normal glass-formers. These predictions are verified by recent measurement of water's diffusion at very low temperatures (point 1) and discoveries from our extensive simulations (points 2-3). We suggest that the same scenario may generally apply to water-like anomalies in liquids tending to form locally favored structures, including not only other important tetrahedral liquids such as silicon, germanium, and silica, but also metallic and chalcogenide liquids.

Evidence of a low-temperature dynamical transition in concentrated microgels

A low-temperature dynamical transition has been reported in several proteins. We provide the first observation of a "protein-like" dynamical transition in nonbiological aqueous environments. To this aim, we exploit the popular colloidal system of poly-N-isopropylacrylamide (PNIPAM) microgels, extending their investigation to unprecedentedly high concentrations. Owing to the heterogeneous architecture of the microgels, water crystallization is avoided in concentrated samples, allowing us to monitor atomic dynamics at low temperatures. By elastic incoherent neutron scattering and molecular dynamics simulations, we find that a dynamical transition occurs at a temperature T _d ~ 250 K, independently from PNIPAM mass fraction. However, the transition is smeared out on approaching dry conditions. The quantitative agreement between experiments and simulations provides evidence that the transition occurs simultaneously for PNIPAM and water dynamics. The similarity of these results with hydrated protein powders suggests that the dynamical transition is a generic feature in complex macromolecular systems, independently from their biological function.

On the time required to freeze water

By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avrami's expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude.