Enhanced small-angle scattering connected to the Widom line in simulations of supercooled water

Liquid-Liquid Crossover in Water Model: Local Structure vs Kinetics of Hydrogen Bonds

In equilibrium and supercooled liquids, polymorphism is manifested by thermodynamic regions defined in the phase diagram, which are predominantly of different short- and medium-range order (local structure). It is found that on the phase diagram of the water model, the thermodynamic region corresponding to the equilibrium liquid phase is divided by a line of the smooth liquid-liquid crossover. In the case of the water model TIP4P/2005, this crossover is revealed by various local order parameters and corresponds to pressures on the order of 3150 ± 350 atm at ambient temperature. In the vicinity of the crossover, the dynamics of water molecules change significantly, which is reflected, in particular, in the fact that the self-diffusion coefficient reaches its maximum values. In addition, changes in the structure also manifest themselves in changes in the kinetics of hydrogen bonding, which are captured by values of such quantities as the average lifetime of hydrogen bonding, the average lifetimes of different local coordination numbers, and the frequencies of changes in different local coordination numbers. An interpretation of the hydrogen bond kinetics in terms of the free energy landscape concept in the space of possible coordination numbers is proposed.

Charge Gradients around Dendritic Voids Cause Nanoscale Inhomogeneities in Liquid Water

Water is the matrix of life and serves as a solvent for numerous physical and chemical processes. The origins of the nature of inhomogeneities that exist in liquid water and the time scales over which they occur remains an open question. Here, we report femtosecond elastic second harmonic scattering (fs-ESHS) of liquid water in comparison to an isotropic liquid (CCl₄) and show that water is indeed a nonuniform liquid. The coherent fs-ESHS intensity was interpreted, using molecular dynamics simulations, as arising from charge density fluctuations with enhanced nanoscale polarizabilities around transient voids having an average lifetime of 300 fs. Although voids were also present in CCl₄, they were not characterized by hydrogen bond defects and did not show strong polarizability fluctuations, leading to fs-ESHS of an isotropic liquid. The voids increased in number at higher temperatures above room temperature, in agreement with the fs-ESHS results.

Anomalous temperature dependence of the experimental x-ray structure factor of supercooled water

The structural changes of water upon deep supercooling were studied through wide-angle x-ray scattering at SwissFEL. The experimental setup had a momentum transfer range of 4.5 Å^-1, which covered the principal doublet of the x-ray structure factor of water. The oxygen-oxygen structure factor was obtained for temperatures down to 228.5 ± 0.6 K. Similar to previous studies, the second diffraction peak increased strongly in amplitude as the structural change accelerated toward a local tetrahedral structure upon deep supercooling. We also observed an anomalous trend for the second peak position of the oxygen-oxygen structure factor (q₂). We found that q₂ exhibits an unprecedented positive partial derivative with respect to temperature for temperatures below 236 K. Based on Fourier inversion of our experimental data combined with reference data, we propose that the anomalous q₂ shift originates from that a repeat spacing in the tetrahedral network, associated with all peaks in the oxygen-oxygen pair-correlation function, gives rise to a less dense local ordering that resembles that of low-density amorphous ice. The findings are consistent with that liquid water consists of a pentamer-based hydrogen-bonded network with low density upon deep supercooling.

Molecular Perspective on Water Vapor Accommodation into Ice and Its Dependence on Temperature

Accommodation of vapor-phase water molecules into ice crystal surfaces is a fundamental process controlling atmospheric ice crystal growth. Experimental studies investigating the accommodation process with various techniques report widely spread values of the water accommodation coefficient on ice, α_ice, and the results on its potential temperature dependence are inconclusive. We run molecular dynamics simulations of molecules condensing onto the basal plane of ice I_h using the TIP4P/Ice empirical force field and characterize the accommodated state from this molecular perspective, utilizing the interaction energy, the tetrahedrality order parameter, and the distance below the instantaneous interface as criteria. Changes of the order parameter turn out to be a suitable measure to distinguish between the surface and bulk states of a molecule condensing onto the disordered interface. In light of the findings from the molecular dynamics, we discuss and re-analyze a recent experimental data set on α_ice obtained with an environmental molecular beam (EMB) setup [KongX.; J. Phys. Chem. A2014, 118 ( (22), ), 3973−397924814567] using kinetic molecular flux modeling, aiming at a more comprehensive picture of the accommodation process from a molecular perspective. These results indicate that the experimental observations indeed cannot be explained by evaporation alone. At the same time, our results raise the issue of rapidly growing relaxation times upon decreasing temperature, challenging future experimental efforts to cover relevant time scales. Finally, we discuss the relevance of the water accommodation coefficient on ice in the context of atmospheric cloud particle growth processes.

The nano-structural inhomogeneity of dynamic hydrogen bond network of TIP4P/2005 water

A method for studying the time dependence of the short-range molecular order of water has been proposed. In the present study, water is considered as a dynamic network between molecules at distances not exceeding 3.2 Å. The instantaneous configurations obtained with the molecular dynamics method have been sequentially analyzed. The mutual orientation of each molecule with its neighboring molecules has been studied and the interaction energy of each pair of neighbor molecules has been calculated. The majority of mutual orientation angles between molecules lie in the interval [0°; 20°]. More than 85% of the molecular pairs in each instantaneous configuration form H-bonds and the H-bond network includes all water molecules in the temperature range 233-293 K. The number of H-bonds fluctuates near the mean value and increases with decreasing temperature, and the energy of the vast majority of such bonds is much higher than the thermal energy. The interaction energy of 80% of the H-bonding molecular pairs lies in the interval [-7; -4] kcal/mol. The interaction energy of pairs that do not satisfy the H-bond angle criterion lies in the interval [-5; 4] kcal/mol; the number of such bonds does not exceed 15% and decreases with decreasing temperature. For the first time it has been found that in each instantaneous configuration the H-bond network contains built-in nanometric structural heterogeneities formed by shorter H-bonds. The fraction of molecules involved in the structural heterogeneities increases from 40% to 60% with a temperature decrease from 293 K to 233 K. Each heterogeneity has a finite lifetime and changeable structure, but they are constantly present during the entire simulation time.

A proposal for the structure of high- and low-density fluctuations in liquid water

Based on recent experimental data that can be interpreted as indicating the presence of specific structures in liquid water, we build and optimize two structural models which we compare with the available experimental data. To represent the proposed high-density liquid structures, we use a model consisting of chains of water molecules, and for low-density liquid, we investigate fused dodecahedra as templates for tetrahedral fluctuations. The computed infrared spectra of the models are in very good agreement with the extracted experimental spectra for the two components, while the extracted structures from molecular dynamics (MD) simulations give spectra that are intermediate between the experimentally derived spectra. Computed x-ray absorption and emission spectra as well as the O-O radial distribution functions of the proposed structures are not contradicted by experiment. The stability of the proposed dodecahedral template structures is investigated in MD simulations by seeding the starting structure, and remnants found to persist on an ∼30 ps time scale. We discuss the possible significance of such seeds in simulations and whether they can be viable candidates as templates for structural fluctuations below the compressibility minimum of liquid water.

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.

High and low density patches in simulated liquid water

We present insights into the nature of structural heterogeneities in liquid water by characterizing the empty space within the hydrogen bond network. Using molecular dynamics simulations, we show that density fluctuations create regions of empty space characterized by a diverse morphology - from spherical to fractal-like voids. These voids allow for the identification of low and high density patches of the liquid, encompassing short (0.3-0.5 nm) as well as long (1-2 nm) length-scales. In addition, we show that the formation of these patches is coupled to collective fluctuations involving the topology of hydrogen-bonded rings of water molecules. In particular, water molecules in the high density patches tend to be slightly more tetrahedral - which is consistent with the predictions of the hydrophobic effect.

Coherent X-rays reveal the influence of cage effects on ultrafast water dynamics

The dynamics of liquid water feature a variety of time scales, ranging from extremely fast ballistic-like thermal motion, to slower molecular diffusion and hydrogen-bond rearrangements. Here, we utilize coherent X-ray pulses to investigate the sub-100 fs equilibrium dynamics of water from ambient conditions down to supercooled temperatures. This novel approach utilizes the inherent capability of X-ray speckle visibility spectroscopy to measure equilibrium intermolecular dynamics with lengthscale selectivity, by measuring oxygen motion in momentum space. The observed decay of the speckle contrast at the first diffraction peak, which reflects tetrahedral coordination, is attributed to motion on a molecular scale within the first 120 fs. Through comparison with molecular dynamics simulations, we conclude that the slowing down upon cooling from 328 K down to 253 K is not due to simple thermal ballistic-like motion, but that cage effects play an important role even on timescales over 25 fs due to hydrogen-bonding.

The temperature dependence of intermediate range oxygen-oxygen correlations in liquid water

We analyze the recent temperature dependent oxygen-oxygen pair-distribution functions from experimental high-precision x-ray diffraction data of bulk water by Skinner et al. [J. Chem. Phys. 141, 214507 (2014)] with particular focus on the intermediate range where small, but significant, correlations are found out to 17 Å. The second peak in the pair-distribution function at 4.5 Å is connected to tetrahedral coordination and was shown by Skinner et al. to change behavior with temperature below the temperature of minimum isothermal compressibility. Here we show that this is associated also with a peak growing at 11 Å which strongly indicates a collective character of fluctuations leading to the enhanced compressibility at lower temperatures. We note that the peak at ∼13.2 Å exhibits a temperature dependence similar to that of the density with a maximum close to 277 K or 4 °C. We analyze simulations of the TIP4P/2005 water model in the same manner and find excellent agreement between simulations and experiment albeit with a temperature shift of ∼20 K.

Probing the triplet correlation function in liquid water by experiments and molecular simulations

Three-body information of liquid water is extracted using X-ray diffraction experiment as well as in molecular simulations via isothermal pressure derivative of structure factor term.

Water: A Tale of Two Liquids

Water is the most abundant liquid on earth and also the substance with the largest number of anomalies in its properties. It is a prerequisite for life and as such a most important subject of current research in chemical physics and physical chemistry. In spite of its simplicity as a liquid, it has an enormously rich phase diagram where different types of ices, amorphous phases, and anomalies disclose a path that points to unique thermodynamics of its supercooled liquid state that still hides many unraveled secrets. In this review we describe the behavior of water in the regime from ambient conditions to the deeply supercooled region. The review describes simulations and experiments on this anomalous liquid. Several scenarios have been proposed to explain the anomalous properties that become strongly enhanced in the supercooled region. Among those, the second critical-point scenario has been investigated extensively, and at present most experimental evidence point to this scenario. Starting from very low temperatures, a coexistence line between a high-density amorphous phase and a low-density amorphous phase would continue in a coexistence line between a high-density and a low-density liquid phase terminating in a liquid–liquid critical point, LLCP. On approaching this LLCP from the one-phase region, a crossover in thermodynamics and dynamics can be found. This is discussed based on a picture of a temperature-dependent balance between a high-density liquid and a low-density liquid favored by, respectively, entropy and enthalpy, leading to a consistent picture of the thermodynamics of bulk water. Ice nucleation is also discussed, since this is what severely impedes experimental investigation of the vicinity of the proposed LLCP. Experimental investigation of stretched water, i.e., water at negative pressure, gives access to a different regime of the complex water diagram. Different ways to inhibit crystallization through confinement and aqueous solutions are discussed through results from experiments and simulations using the most sophisticated and advanced techniques. These findings represent tiles of a global picture that still needs to be completed. Some of the possible experimental lines of research that are essential to complete this picture are explored.

Dynamical properties of water-methanol solutions

We study the relaxation times tα in the water-methanol system. We examine new data and data from the literature in the large temperature range 163 < T < 335 K obtained using different experimental techniques and focus on how tα affects the hydrogen bond structure of the system and the hydrophobicity of the alcohol methyl group. We examine the relaxation times at a fixed temperature as a function of the water molar fraction XW and observe two opposite behaviors in their curvature when the system moves from high to low T regimes. This behavior differs from that of an ideal solution in that it has excess values located at different molar fractions (XW = 0.5 for high T and 0.75 in the deep supercooled regime). We analyze the data and find that above a crossover temperature T ∼ 223 K, hydrophobicity plays a significant role and below it the water tetrahedral network dominates. This temperature is coincident with the fragile-to-strong dynamical crossover observed in confined water and supports the liquid-liquid phase transition hypothesis. At the same time, the reported data suggest that this crossover temperature (identified as the Widom line temperature) also depends on the alcohol concentration.

The structural origin of anomalous properties of liquid water

Water is unique in its number of unusual, often called anomalous, properties. When hot it is a normal simple liquid; however, close to ambient temperatures properties, such as the compressibility, begin to deviate and do so increasingly on further cooling. Clearly, these emerging properties are connected to its ability to form up to four well-defined hydrogen bonds allowing for different local structural arrangements. A wealth of new data from various experiments and simulations has recently become available. When taken together they point to a heterogeneous picture with fluctuations between two classes of local structural environments developing on temperature-dependent length scales.

Continuous and Discontinuous Dynamic Crossover in Supercooled Water in Computer Simulations

The dynamic crossover behavior of supercooled water as described by the first-principle based WAIL potential was investigated. Below the second liquid-liquid critical point, the viscosity shows a discontinuous jump consistent with a first-order phase transition between the high density liquid and the low density liquid. Above the critical point, a continuous transition occurs with only the first derivative of viscosity being discontinuous, and the dynamic crossover temperature is about 8 K below the thermodynamic switchover temperature. The 8 K shift can be explained by a delay in dynamic crossover, which does not occur until the more viscous liquid starts to dominate the population and jams the flow. On the basis of finite-size effects observed in our simulations, we believe that dynamic discontinuity may be observable above the critical point in confined water when the confinement is on a length scale shorter than the spatial correlation.

Thermodynamic properties of bulk and confined water

The thermodynamic response functions of water display anomalous behaviors. We study these anomalous behaviors in bulk and confined water. We use nuclear magnetic resonance (NMR) to examine the configurational specific heat and the transport parameters in both the thermal stable and the metastable supercooled phases. The data we obtain suggest that there is a behavior common to both phases: that the dynamics of water exhibit two singular temperatures belonging to the supercooled and the stable phase, respectively. One is the dynamic fragile-to-strong crossover temperature (T(L) ≃ 225 K). The second, T* ∼ 315 ± 5 K, is a special locus of the isothermal compressibility K(T)(T, P) and the thermal expansion coefficient α(P)(T, P) in the P-T plane. In the case of water confined inside a protein, we observe that these two temperatures mark, respectively, the onset of protein flexibility from its low temperature glass state (T(L)) and the onset of the unfolding process (T*).

Triplet correlation functions in liquid water

Triplet correlations have been shown to play a crucial role in the transformation of simple liquids to anomalous tetrahedral fluids [M. Singh, D. Dhabal, A. H. Nguyen, V. Molinero, and C. Chakravarty, Phys. Rev. Lett. 112, 147801 (2014)]. Here we examine triplet correlation functions for water, arguably the most important tetrahedral liquid, under ambient conditions, using configurational ensembles derived from molecular dynamics (MD) simulations and reverse Monte Carlo (RMC) datasets fitted to experimental scattering data. Four different RMC data sets with widely varying hydrogen-bond topologies fitted to neutron and x-ray scattering data are considered [K. T. Wikfeldt, M. Leetmaa, M. P. Ljungberg, A. Nilsson, and L. G. M. Pettersson, J. Phys. Chem. B 113, 6246 (2009)]. Molecular dynamics simulations are performed for two rigid-body effective pair potentials (SPC/E and TIP4P/2005) and the monatomic water (mW) model. Triplet correlation functions are compared with other structural measures for tetrahedrality, such as the O-O-O angular distribution function and the local tetrahedral order distributions. In contrast to the pair correlation functions, which are identical for all the RMC ensembles, the O-O-O triplet correlation function can discriminate between ensembles with different degrees of tetrahedral network formation with the maximally symmetric, tetrahedral SYM dataset displaying distinct signatures of tetrahedrality similar to those obtained from atomistic simulations of the SPC/E model. Triplet correlations from the RMC datasets conform closely to the Kirkwood superposition approximation, while those from MD simulations show deviations within the first two neighbour shells. The possibilities for experimental estimation of triplet correlations of water and other tetrahedral liquids are discussed.

Thermal Conductivity of Supercooled Water: An Equilibrium Molecular Dynamics Exploration

The thermal conductivity of both supercooled and ambient-temperature water at atmospheric pressure has been computed over the 140-270 K temperature range for three popular water models via equilibrium molecular dynamics in the Green-Kubo setting. No strong temperature dependence of thermal conductivity was observed. The underlying phonon modes contributing to thermal conduction processes have been examined in the present work, and it has been established that (translational) acoustic modes dominate in supercooled water.

Structural Evolution of Supercritical CO2 across the Frenkel Line

Here, we study structural properties of the supercritical carbon dioxide and discover the existence of persistent medium-range order correlations, which make supercritical carbon dioxide nonuniform and heterogeneous on an intermediate length scale. We report on the CO2 heterogeneity shell structure where, in the first shell, both carbon and oxygen atoms experience gas-like-type interactions with short-range order correlations while within the second shell, oxygen atoms essentially exhibit a liquid-like type of interactions due to localization of transverse-like phonon packets. Importantly, we highlight a catalytic role of atoms inside of the nearest-neighbor heterogeneity shell in providing a mechanism for diffusion and proving the existence of an additional thermodynamic boundary in the supercritical carbon dioxide on an intermediate length scale. Finally, we discuss important implications for answering the intriguing question whether Venus may have had CO2 oceans and urge for an experimental detection of this persistent local-order heterogeneity.

An analysis of fluctuations in supercooled TIP4P/2005 water

Large-scale, long-time molecular dynamics simulations are used to investigate fluctuations in the TIP4P/2005 water model in the supercooled region (240-190 K). Particular attention is focused in the vicinity of a previously reported liquid-liquid critical point [J. L. F. Abascal and C. Vega, J. Chem. Phys. 133, 234502 (2010)]. Water is viewed as an equimolar binary mixture with "species" defined based on a local tetrahedral order parameter. A Bhatia-Thornton fluctuation analysis is used to show that species concentration fluctuations couple to density fluctuations and completely account for the anomalous increase in the structure factor at small wave number observed under supercooled conditions. Although we find that both concentration and density fluctuations increase with decreasing temperature along the proposed critical isochore, we cannot confirm the existence of a liquid-liquid critical point. Our simulations suggest that the parameters previously reported are not a true liquid-liquid critical point and we find no evidence of two-phase coexistence in its vicinity. It is shown that very long simulations (on the order of 8 μs for 500 molecules) are necessary to obtain well converged density distributions for deeply supercooled water and this is especially important if one is seeking direct evidence of a two-phase region.

The Boson peak in supercooled water

We perform extensive molecular dynamics simulations of the TIP4P/2005 model of water to investigate the origin of the Boson peak reported in experiments on supercooled water in nanoconfined pores, and in hydration water around proteins. We find that the onset of the Boson peak in supercooled bulk water coincides with the crossover to a predominantly low-density-like liquid below the Widom line TW. The frequency and onset temperature of the Boson peak in our simulations of bulk water agree well with the results from experiments on nanoconfined water. Our results suggest that the Boson peak in water is not an exclusive effect of confinement. We further find that, similar to other glass-forming liquids, the vibrational modes corresponding to the Boson peak are spatially extended and are related to transverse phonons found in the parent crystal, here ice Ih.

Microscopic structure of water at elevated pressures and temperatures

We report on the microscopic structure of water at sub- and supercritical conditions studied using X-ray Raman spectroscopy, ab initio molecular dynamics simulations, and density functional theory. Systematic changes in the X-ray Raman spectra with increasing pressure and temperature are observed. Throughout the studied thermodynamic range, the experimental spectra can be interpreted with a structural model obtained from the molecular dynamics simulations. A spatial statistical analysis using Ripley's K-function shows that this model is homogeneous on the nanometer length scale. According to the simulations, distortions of the hydrogen-bond network increase dramatically when temperature and pressure increase to the supercritical regime. In particular, the average number of hydrogen bonds per molecule decreases to ≈ 0.6 at 600 °C and p = 134 MPa.

Transport coefficients of the TIP4P-2005 water model

A detailed understanding of the dynamics of liquid water at molecular level is of fundamental importance as well as have applications in many branches of science and technology. In this work, the diffusion of the TIP4P-2005 model of water is systematically investigated in liquid phase in the temperature range 210-310 K. The translational and rotational diffusions, as well as correlations between them, are examined. The effects of system size and shape are also probed in this study. The results suggest the presence of a temperature of dynamical arrest of molecular translations in the range of 150-180 K and of molecular rotations in the range of 80-130 K, depending on specific direction. A substantial change in the preferred directions of translations and rotations relative to the molecular coordinate system is observed slightly below (≈15 K) the melting temperature of the model. It is shown that there is a correlation between translational and rotational molecular motions essential for diffusion in the liquid. The presence of hydrodynamic size effects is confirmed and quantified; it is also shown that using a non-cubic simulation box for a liquid system leads to an anisotropic splitting in the diffusion tensor. The findings of this study enhance our general understanding of models of water, specifically the TIP4P-2005 model, as well as provide evidences of the direct connection between thermodynamics of liquid water and dynamics of its molecules.

Widom line and noise-power spectral analysis of a supercritical fluid

We have performed extensive molecular dynamics simulations to study noise-power spectra of density and potential energy fluctuations of a Lennard-Jones model of a fluid in the supercritical region. Emanating from the liquid-vapor critical point, there is a locus of isobaric specific heat maxima, called the Widom line, which is often regarded as an extension of the liquid-vapor coexistence line. Our simulation results show that the noise-power spectrum of the density fluctuations on the Widom line of the liquid-vapor transition exhibits three distinct 1/f^{γ} behaviors with exponents γ=0, 1.2, and 2, depending on the frequency f. We find that the intermediate frequency region with an exponent γ∼ 1 appears as the temperature approaches the Widom temperature from above or below. On the other hand, we do not find three distinct regions of 1/f^{γ} in the power spectrum of the potential energy fluctuations on the Widom line. Furthermore, we find that the power spectra of both the density and potential energy fluctuations at low frequency have a maximum on the Widom line, suggesting that the noise power can provide an alternative signature of the Widom line.