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

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.

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.

Complex phase diagram and supercritical matter

The supercritical region is often described as uniform with no definite transitions. The distinct behaviors of the matter therein, e.g., as liquidlike and gaslike, however, suggest "supercritical boundaries." Here we provide a mathematical description of these phenomena by revisiting the Yang-Lee theory and introducing a complex phase diagram, specifically a four-dimensional (4D) one with complex T and p. While the traditional 2D phase diagram with real temperature T and pressure p values (the physical plane) lacks Lee-Yang (LY) zeros beyond the critical point, preventing the occurrence of criticality, the off-plane zeros in this 4D scenario still induce critical anomalies in various physical properties. This relationship is evidenced by the correlation between the Widom line and LY edges in van der Waals, 2D Ising model, and water. The diverged supercritical boundaries manifest the high-dimensional feature of the phase diagram: e.g., when LY zeros of complex T or p are projected onto the physical plane, boundaries defined by isobaric heat capacity C_{p} or isothermal compression coefficient K_{T} emanates. These results demonstrate the incipient phase transition nature of the supercritical matter.

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.

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.

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.

Minimum in the pressure dependence of the interfacial free energy between ice Ih and water

Despite the importance of ice nucleation, this process has been barely explored at negative pressures. Here, we study homogeneous ice nucleation in stretched water by means of molecular dynamics seeding simulations using the TIP4P/Ice model. We observe that the critical nucleus size, interfacial free energy, free energy barrier, and nucleation rate barely change between isobars from -2600 to 500 bars when they are represented as a function of supercooling. This allows us to identify universal empirical expressions for homogeneous ice nucleation in the pressure range from -2600 to 500 bars. We show that this universal behavior arises from the pressure dependence of the interfacial free energy, which we compute by means of the mold integration technique, finding a shallow minimum around -2000 bars. Likewise, we show that the change in the interfacial free energy with pressure is proportional to the excess entropy and the slope of the melting line, exhibiting in the latter a reentrant behavior also at the same negative pressure. Finally, we estimate the excess internal energy and the excess entropy of the ice Ih-water interface.

Molecular dynamics analysis of elastic properties and new phase formation during amorphous ices transformations

Unlike conventional first-order phase transitions, the kinetics of amorphous-amorphous transitions has been much less studied. The ultrasonic experiments on the transformations between low-density and high-density amorphous ice induced by pressure or heating provided the pressure and temperature dependencies of elastic moduli. In this article, we make an attempt to build a microscopic picture of these experimentally studied transformations using the molecular dynamics method with the TIP4P/Ice water model. We study carefully the dependence of the results of elastic constants calculations on the deformation rates. The system size effects are considered as well. The comparison with the experimental data enriches our understanding of the transitions observed. Our modeling gives new information about the formation mechanisms of new phase clusters during the transition between low-density and high-density amorphous ices. We analyse the applicability of the term “nucleation” for these processes.

Modeling Solution Behavior of Poly(N-isopropylacrylamide): A Comparison between Water Models

Water is known to play a fundamental role in determining the structure and functionality of macromolecules. The same crucial contribution is also found in the in silico description of polymer aqueous solutions. In this work, we exploit the widely investigated synthetic polymer poly(N-isopropylacrylamide) (PNIPAM) to understand the effect of the adopted water model on its solution behavior and to refine the computational setup. By means of atomistic molecular dynamics simulations, we perform a comparative study of PNIPAM aqueous solution using two advanced water models: TIP4P/2005 and TIP4P/Ice. The conformation and hydration features of an atactic 30-mer at infinite dilution are probed at a range of temperature and pressure suitable to detect the coil-to-globule transition and to map the P–T phase diagram. Although both water models can reproduce the temperature-induced coil-to-globule transition at atmospheric pressure and the polymer hydration enhancement that occurs with increasing pressure, the PNIPAM–TIP4P/Ice solution shows better agreement with experimental findings. This result can be attributed to a stronger interaction of TIP4P/Ice water with both hydrophilic and hydrophobic groups of PNIPAM, as well as to a less favorable contribution of the solvent entropy to the coil-to-globule transition.

Is It Possible to Follow the Structural Evolution of Water in "No-Man's Land" Using a Pulsed-Heating Procedure?

The anomalous increase in compressibility and heat capacity of supercooled water has been attributed to its structural transformation of into a four-coordinated liquid. Experiments revealed that κ_T and C_p peak at T_W^thermo ≈ 229 K [Kim et al. Science 2017, 358, 1589; Pathak et al. Proc. Natl. Acad. Sci. 2021, 118, e2018379118]. Recently, a pulsed heating procedure (PHP) was employed to interrogate the structure of water, reporting a steep increase in tetrahedrality around T_W^PHP = 210 ± 3 K [Kringle et al. Science 2020, 369, 1490]. This discrepancy questions whether water structure and thermodynamics are decoupled, or if the shift in T_W is an artifact of PHP. Here we implement PHP in molecular simulations. We find that the stationary states captured at the bottom of the pulse are not representative of the thermalized liquid or its inherent structure. Our analysis reveals a temperature-dependent distortion that shifts T_W^PHP to ∼20 K below T_W^thermo. We conclude that 2 orders of magnitude faster rates are required to sample water's inherent structure with PHP.

Thermodynamics and kinetics of crystallization in deeply supercooled Stillinger-Weber silicon

We study the kinetics of crystallization in deeply supercooled liquid silicon employing computer simulations and the Stillinger-Weber three-body potential. The free energy barriers to crystallization are computed using umbrella sampling Monte Carlo simulations and from unconstrained molecular dynamics simulations using a mean first passage time formulation. We focus on state points that have been described in earlier work [S. Sastry and C. A. Angell, Nat. Mater. 2, 739 (2003)] as straddling a liquid-liquid phase transition (LLPT) between two metastable liquid states. It was argued subsequently [Ricci et al., Mol. Phys. 117, 3254 (2019)] that the apparent transition is due to the loss of metastability of the liquid state with respect to the crystalline state. The presence of a barrier to crystallization for these state points is therefore of importance to ascertain, which we investigate, with due attention to ambiguities that may arise from the choice of order parameters. We find a well-defined free energy barrier to crystallization and demonstrate that both umbrella sampling and mean first passage time methods yield results that agree quantitatively. Our results thus provide strong evidence against the possibility that the liquids at state points close to the reported LLPT exhibit slow, spontaneous crystallization, but they do not address the existence of a LLPT (or lack thereof). We also compute the free energy barriers to crystallization at other state points over a broad range of temperatures and pressures and discuss the effect of changes in the microscopic structure of the metastable liquid on the free energy barrier heights.