Dynamic crossover in deeply cooled water confined in MCM-41 at 4 kbar and its relation to the liquid-liquid transition hypothesis

Complex dynamics of partially freezable confined water revealed by combined experimental and computational studies

We investigate water dynamics in mesoporous silica across partial crystallization by combining broadband dielectric spectroscopy (BDS), nuclear magnetic resonance (NMR), and molecular dynamics simulations (MDS). Exploiting the fact that not only BDS but also NMR field-cycling relaxometry and stimulated-echo experiments provide access to dynamical susceptibilities in broad frequency and temperature ranges, we study both the fully liquid state above the melting point Tm and the dynamics of coexisting water and ice phases below this temperature. It is found that partial crystallization leads to a change in the temperature dependence of rotational correlation times τ, which occurs in addition to previously reported dynamical crossovers of confined water and depends on the pore diameter. Furthermore, we observe that dynamical susceptibilities of water are strongly asymmetric in the fully liquid state, whereas they are much broader and nearly symmetric in the partially frozen state. Finally, water in the nonfreezable interfacial layer below Tm does not exhibit a much debated dynamical crossover at ∼220 K. We argue that its dynamics is governed by a static energy landscape, which results from the interaction with the bordering silica and ice surfaces and features a Gaussian-like barrier distribution. Consistently, our MDS analysis of the motional mechanism reveals a hopping motion of water in thin interfacial layers. The rotational correlation times of the confined ice phases follow Arrhenius laws. While the values of τ depend on the pore diameter, freezable water in various types of confinements and mixtures shows similar activation energies of Ea ≈ 0.43 eV.

Dynamical Susceptibilities of Confined Water from Room Temperature to the Glass Transition

We confine water to narrow silica pores, where crystallization is suppressed, and determine the dynamical susceptibilities of the liquid from room temperature down to the glass transition by combining broadband dielectric spectroscopy (BDS) with ¹H and ²H nuclear magnetic resonance (NMR), in particular, by establishing NMR field-cycling relaxometry. For the correlation times, derivative analysis reveals Vogel-Fulcher-Tammann and Arrhenius regimes at T ≥ 215 K and T ≤ 160 K, respectively, which are separated by a broad crossover region. The continuous transition in the temperature dependence is accompanied by a gradual change from asymmetric high-temperature shapes of the dynamical susceptibilities to symmetric low-temperature ones and by a steady decrease of the dielectric relaxation strength. In the Arrhenius regime (E_a = 0.48 eV) at T ≤ 160 K, 2D ²H NMR spectra reveal quasi-isotropic water reorientation. We rationalize these results in terms of a crossover to an interface-affected, noncooperative relaxation involving both rotational and translational motions.

Non-monotonic composition dependence of the breakdown of Stokes-Einstein relation for water in aqueous solutions of ethanol and 1-propanol: explanation using translational jump-diffusion approach

A series of experimental and simulation studies examined the validity of the Stokes-Einstein relationship (SER) of water in binary water/alcohol mixtures of different mixture compositions. These studies revealed a strong non-monotonic composition dependence of the SER with maxima at the specific alcohol mole fraction where the non-idealities of the thermodynamic and transport properties are observed. The translational jump-diffusion (TJD) approach elucidated the breakdown of the SER in pure supercooled water as caused by the jump translation of molecules. The breakdown of SER in the supercooled water/methanol binary mixture was successfully explained using the same TJD approach. To further generalize the picture, here we focus on the non-monotonic composition dependence of SER breakdown of water in two water/alcohol mixtures (water/ethanol and water/propanol) for a broad temperature range. In agreement with previous studies, maximum breakdown of SER is observed for the mixture with alcohol mole fraction x = 0.2. Diffusion of the water molecules at the maximum SER breakdown point is largely contributed by jump-diffusion. The residual-diffusion, obtained by subtracting the jump-diffusion from the total diffusion, approximately follows the SER for different compositions and temperatures. We also performed hydrogen (H-)bond dynamics and observed that the contribution of jump-diffusion is proportional to the total free energy of activation of breaking all H-bonds that exist around a molecule. This study, therefore, suggests that the more a molecule is trapped by H-bonding, the more likely it is to diffuse through the jump-diffusion mechanism, eventually leading to an increasing degree of SER breakdown.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

The anomalies and criticality of liquid water

The origin of water's anomalies has been a matter of long-standing debate. A two-state model, dating back to Röntgen, relies on the dynamical coexistence of two types of local structures-locally favored tetrahedral structure (LFTS) and disordered normal-liquid structure (DNLS)-in liquid water. Phenomenologically, this model not only explains water's thermodynamic anomalies but also can rationalize the existence of a liquid-liquid critical point (LLCP) if there is a cooperative formation of LFTS. We recently found direct evidence for the coexistence of LFTS and DNLS in the experimental structure factor of liquid water. However, the existence of the LLCP and its impact on water's properties has remained elusive, leaving the origin of water's anomalies unclear. Here we propose a unique strategy to locate the LLCP of liquid water. First, we make a comprehensive analysis of a large set of experimental structural, thermodynamic, and dynamic data based on our hierarchical two-state model. This model predicts that the two thermodynamic and dynamical fluctuation maxima lines should cross at the LLCP if it exists, which we confirm by hundred-microsecond simulations for model waters. Based on recent experimental results of the compressibility and diffusivity measurements in the no man's land, we reveal that the two lines cross around 184 K and 173 MPa for real water, suggesting the presence of the LLCP around there. Nevertheless, we find that the criticality is almost negligible in the experimentally accessible region of liquid water because it is too far from the LLCP. Our findings would provide a clue to settle the long-standing debate.

Static field gradient NMR studies of water diffusion in mesoporous silica

NMR diffusometry is used to ascertain the pore-size dependent water diffusion in MCM-41 and SBA-15 silica over broad temperature ranges. Detailed analysis of 1H and 2H NMR stimulated-echo decays reveals that fast water motion through voids between different silica particles impairs such studies in the general case. However, water diffusion inside single pores is probed in the present approach, which applies high static field gradients to enhance the spatial resolution of the experiment and uses excess water in combination with subzero temperatures to embed the silica particles in an ice matrix and, thus, to suppress interparticle water motion. It is found that the diffusion of confined water slows down by almost two orders of magnitude when the pore diameter is reduced from 5.4 nm to 2.1 nm at weak cooling. In the narrower silica pores, the temperature dependence of the self-diffusion coefficient of water is well described by an Arrhenius law with an activation energy of Ea = 0.40 eV. The Arrhenius behavior extends over a broad temperature range of at least 207-270 K, providing evidence against a fragile-to-strong crossover in response to a proposed liquid-liquid phase transition near 225 K. In the wider silica pores, partial crystallization results in a discontinuous temperature dependence. Explicitly, the diffusion coefficients drop when cooling through the pore-size dependent melting temperatures Tm of confined water. This finding can be rationalized by the fact that water can explore the whole pore volumes above Tm, but is restricted to narrow interfacial layers sandwiched between silica walls and ice crystallites below this temperature. Comparing our findings for water diffusion with previous results for water reorientation, we find significantly different temperature dependencies, indicating that the Stokes-Einstein-Debye relation is not obeyed.

The Proton Density of States in Confined Water (H2O)

The hydrogen density of states (DOS) in confined water has been probed by inelastic neutron scattering spectra in a wide range of its P–T phase diagram. The liquid–liquid transition and the dynamical crossover from the fragile (super-Arrhenius) to strong (Arrhenius) glass forming behavior have been studied, by taking into account the system polymorphism in both the liquid and amorphous solid phases. The interest is focused in the low energy region of the DOS (E<10 meV) and the data are discussed in terms of the energy landscape (local minima of the potential energy) approach. In this latest research, we consider a unit scale energy (EC) linked to the water local order governed by the hydrogen bonding (HB). All the measured spectra, scaled according to such energy, evidence a universal power law behavior with different exponents (γ) in the strong and fragile glass forming regions, respectively. In the first case, the DOS data obey the Debye squared-frequency law, whereas, in the second one, we obtain a value predicted in terms of the mode-coupling theory (MCT) (γ≃1.6).

Unravelling the contribution of local structures to the anomalies of water: The synergistic action of several factors

We investigate the microscopic origin of water's anomalies by inspecting the hydrogen bond network (HBN) and the spatial organization of low-density-liquid (LDL) like and high-density-liquid (HDL) like environments. Specifically, we simulate-via classical molecular dynamics simulations-the isobaric cooling of a sample composed of 512 water molecules from ambient to deeply undercooled conditions at three pressures, namely, 1 bar, 400 bars, and 1000 bars. In correspondence with the Widom line (WL), (i) the HDL-like dominating cluster undergoes fragmentation caused by the percolation of LDL-like aggregates following a spinodal-like kinetics; (ii) such fragmentation always occurs at a "critical" concentration of ∼20%-30% in LDL; (iii) the HBN within LDL-like environments is characterized by an equal number of pentagonal and hexagonal rings that create a state of maximal frustration between a configuration that promotes crystallization (hexagonal ring) and a configuration that hinders it (pentagonal ring); (iv) the spatial organization of HDL-like environments shows a marked variation. Moreover, the inspection of the global symmetry shows that the intermediate-range order decreases in correspondence with the WL and such a decrease becomes more pronounced upon increasing the pressure, hence supporting the hypothesis of a liquid-liquid critical point. Our results reveal and rationalize the complex microscopic origin of water's anomalies as the cooperative effect of several factors acting synergistically. Beyond implications for water, our findings may be extended to other materials displaying anomalous behaviours.

Nanoscale dynamics of water confined in ordered mesoporous carbon

The single particle dynamics of water confined in ordered mesoporous carbon matrix was investigated in the temperature range from 290 K to 170 K by quasielastic neutron scattering using three high resolution neutron spectrometers.

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.

Water-like anomalies as a function of tetrahedrality

Water is the most common and yet least understood material on Earth. Despite its simplicity, water tends to form tetrahedral order locally by directional hydrogen bonding. This structuring is known to be responsible for a vast array of unusual properties, e.g., the density maximum at 4 °C, which play a fundamental role in countless natural and technological processes, with the Earth’s climate being one of the most important examples. By systematically tuning the degree of tetrahedrality, we succeed in continuously interpolating between water-like behavior and simple liquid-like behavior. Our approach reveals what physical factors make water so anomalous and special even compared with other tetrahedral liquids.

Tetrahedral interactions describe the behavior of the most abundant and technologically important materials on Earth, such as water, silicon, carbon, germanium, and countless others. Despite their differences, these materials share unique common physical behaviors, such as liquid anomalies, open crystalline structures, and extremely poor glass-forming ability at ambient pressure. To reveal the physical origin of these anomalies and their link to the shape of the phase diagram, we systematically study the properties of the Stillinger–Weber potential as a function of the strength of the tetrahedral interaction λ. We uncover a unique transition to a reentrant spinodal line at low values of λ, accompanied with a change in the dynamical behavior, from non-Arrhenius to Arrhenius. We then show that a two-state model can provide a comprehensive understanding on how the thermodynamic and dynamic anomalies of this important class of materials depend on the strength of the tetrahedral interaction. Our work establishes a deep link between the shape of the phase diagram and the thermodynamic and dynamic properties through local structural ordering in liquids and hints at why water is so special among all substances.

Counterion-Regulated Dynamics of Water Confined in Lyotropic Liquid Crystalline Morphologies

The dynamics of confined water is of fundamental and long-standing interest. In technologically important forms of confinement, such as proton-exchange membranes, electrostatic interactions with the confining matrix and counterions play significant roles on the properties of water. There has been recent interest on the dynamics of water confined to the lyotropic liquid crystalline (LLC) morphologies of Gemini dicarboxylate surfactants. These systems are exciting because the nature of confinement, for example, size and curvature of channels and surface functionality is dictated by the chemistry of the self-assembling surfactant molecules. Quasielastic neutron scattering experiments have shown an interesting dependence of the water self-diffusion constant, D_α, on the identity (denoted α) of the counterion: at high hydration, the magnitude of the water self-diffusion constant is in the order D_TMA < D_Na < D_K, where TMA, Na, and K refer to tetramethyl ammonium, sodium, and potassium counterions, respectively. This sequence is similar to what is seen in bulk electrolyte solutions. At low hydrations, however, the order of water self-diffusion is different, that is, D_Na < D_TMA < D_K. In this work, we present molecular dynamics simulations for the dynamics of water in the LLC phases of dicarboxylate Gemini surfactants. The simulations reproduce the trends seen in experiments. From an analysis of the trajectories, we hypothesize that two competing factors play a role: the volume accessible to the water molecules and the correlations between the water and the counterion. The excluded volume effect is the largest with TMA⁺, and the electrostatic correlation is the strongest with Na⁺. The observed trend is a result of which of these two effects is dominant at a given water to surfactant ratio.

Hydration-dependent dynamics of water in calcium-silicate-hydrate: A QENS study by global model

HYPOTHESIS

In a saturated cement paste, there are three different types of water: the structural water chemically reacted with cement, the constrained water absorbed to the surface of the pores, and the free water in the center of the pores. Each type has different physicochemical state and unique relation to cement porosity. The different water types have different dynamics which can be detected using quasi-elastic neutron scattering (QENS). Since the porosity of a hardened cement paste is impacted strongly by the water to cement ratio (w/c), it should be possible to extract the hydration dependence of the pores by exploiting the dynamical parameters of the confined water.

EXPERIMENTS

Three C-S-H samples with different water levels, 8%, 17% and 30% were measured using QENS. The measurements were carried out in the scattering vector, Q, range from 0.5 Å^-1 to 1.3 Å^-1, and in the temperature interval from 230 K to 280 K. The data were analyzed using a novel global model developed for cement QENS spectra.

FINDINGS

The results show that while increasing the water content, the structural water index (SWI) decreases and the confining radius, a, increases. Both SWI and a have a linear relationship with the water content. The Arrhenius plot of the translational relaxation time shows that the constrained water dominates the non-structural water at water contents lower than 17%. The rotational activation energy is smaller for lower water content. The analysis demonstrated that our newly proposed global model is practical and useful for analyzing cement QENS data.

Slow Dynamics and Structure of Supercooled Water in Confinement

We review our simulation results on properties of supercooled confined water. We consider two situations: water confined in a hydrophilic pore that mimics an MCM-41 environment and water at interface with a protein. The behavior upon cooling of the α relaxation of water in both environments is well interpreted in terms of the Mode Coupling Theory of glassy dynamics. Moreover, we find a crossover from a fragile to a strong regime. We relate this crossover to the crossing of the Widom line emanating from the liquid-liquid critical point, and in confinement we connect this crossover also to a crossover of the two body excess entropy of water upon cooling. Hydration water exhibits a second, distinctly slower relaxation caused by its dynamical coupling with the protein. The crossover upon cooling of this long relaxation is related to the protein dynamics.

Pressure Effect on the Boson Peak in Deeply Cooled Confined Water: Evidence of a Liquid-Liquid Transition

The boson peak in deeply cooled water confined in nanopores is studied to examine the liquid-liquid transition (LLT). Below ∼180  K, the boson peaks at pressures P higher than ∼3.5  kbar are evidently distinct from those at low pressures by higher mean frequencies and lower heights. Moreover, the higher-P boson peaks can be rescaled to a master curve while the lower-P boson peaks can be rescaled to a different one. These phenomena agree with the existence of two liquid phases with different densities and local structures and the associated LLT in the measured (P, T) region. In addition, the P dependence of the librational band also agrees with the above conclusion.