Thermodynamics of supercooled water

Density isobar of water and melting temperature of ice: Assessing common density functionals

We investigate the density isobar of water and the melting temperature of ice using six different density functionals. Machine-learning potentials are employed to ensure computational affordability. Our findings reveal significant discrepancies between various base functionals. Notably, even the choice of damping can result in substantial differences. Overall, the outcomes obtained through density functional theory are not entirely satisfactory across most utilized functionals. All functionals exhibit significant deviations either in the melting temperature or equilibrium volume, with most of them even predicting an incorrect volume difference between ice and water. Our heuristic analysis indicates that a hybrid functional with 25% exact exchange and van der Waals damping averaged between zero and Becke-Johnson dampings yields the closest agreement with experimental data. This study underscores the necessity for further enhancements in the treatment of van der Waals interactions and, more broadly, density functional theory to enable accurate quantitative predictions for molecular liquids.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

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

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

Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study

Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.

Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks

In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests.

Prediction of water anomalous properties by introducing the two-state theory in SAFT

Water is one of the most abundant substances on earth, but it is still not entirely understood. It shows unusual behavior, and its properties present characteristic extrema unlike any other fluid. This unusual behavior has been linked to the two-state theory of water, which proposes that water forms different clusters, one with a high density and one with a low density, which may even form two distinct phases at low temperatures. Models incorporating the two-state theory manage to capture the unusual extrema of water, unlike traditional equations of state, which fail. In this work, we have derived the framework to incorporate the two-state theory of water into the Statistical-Associating-Fluid-Theory (SAFT). More specifically, we have assumed that water is an ideal solution of high density water molecules and low density water molecules that are in chemical equilibrium. Using this assumption, we have generalized the association term SAFT to allow for the simultaneous existence of the two water types, which have the same physical parameters but different association properties. We have incorporated the newly derived association term in the context of the Perturbed Chain-SAFT (PC-SAFT). The new model is referred to as PC-SAFT-Two-State (PC-SAFT-TS). Using PC-SAFT-TS, we have succeeded in predicting the characteristic extrema of water, such as its density and speed of sound maximum, etc., without loss of accuracy compared to the original PC-SAFT. This new framework is readily extended to mixtures, and PC-SAFT-TS manages to capture the solubility minimum of hydrocarbons in water in a straightforward manner.

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.

Thermodynamic anomalies, polyamorphism and all that

The appearance and evolution of thermodynamics anomalies, and related properties, are studied for two classes of system, modelling those dominated by covalent and ionic interactions, respectively. Such anomalies are most familiar in the density but are also present in other thermodynamics variables such as the compressibility and heat capacity. By systematically varying key model parameters the emergence and evolution of these anomalies can be tracked across the phase space. The interaction of the anomalies can often be rationalized by thermodynamics 'rules'. The emergence of these anomalies may also be correlated with the appearance of polyamorphism, the existence of multiple amorphous states which differ in density and entropy. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

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.

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

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

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

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

Capabilities and limits of autoencoders for extracting collective variables in atomistic materials science

Free energy calculations in materials science are routinely hindered by the need to provide reaction coordinates that can meaningfully partition atomic configuration space, a prerequisite for most enhanced sampling approaches. Recent studies on molecular systems have highlighted the possibility of constructing appropriate collective variables directly from atomic motions through deep learning techniques. Here we extend this class of approaches to condensed matter problems, for which we encode the finite temperature collective variable by an iterative procedure starting from 0 K features of the energy landscape i.e. activation events or migration mechanisms given by a minimum - saddle point - minimum sequence. We employ the autoencoder neural networks in order to build a scalar collective variable for use with the adaptive biasing force method. Particular attention is given to design choices required for application to crystalline systems with defects, including the filtering of thermal motions which otherwise dominate the autoencoder input. The machine-learning workflow is tested on body-centered cubic iron and its common defects, such as small vacancy or self-interstitial clusters and screw dislocations. For localized defects, excellent collective variables as well as derivatives, necessary for free energy sampling, are systematically obtained. However, the approach has a limited accuracy when dealing with reaction coordinates that include atomic displacements of a magnitude comparable to thermal motions, e.g. the ones produced by the long-range elastic field of dislocations. We then combine the extraction of collective variables by autoencoders with an adaptive biasing force free energy method based on Bayesian inference. Using a vacancy migration as an example, we demonstrate the performance of coupling these two approaches for simultaneous discovery of reaction coordinates and free energy sampling in systems with localized defects.

Liquid-Liquid Criticality in the WAIL Water Model

The hypothesis that the anomalous behavior of liquid water is related to the existence of a second critical point in deeply supercooled states has long been the subject of intense debate. Recent, sophisticated experiments designed to observe the transformation between the two subcritical liquids on nano- and microsecond time scales, along with demanding numerical simulations based on classical (rigid) models parametrized to reproduce thermodynamic properties of water, have provided support to this hypothesis. A stronger numerical proof requires demonstrating that the critical point, which occurs at temperatures and pressures far from those at which the models were optimized, is robust with respect to model parameterization, specifically with respect to incorporating additional physical effects. Here we show that a liquid-liquid critical point can be rigorously located also in the WAIL model of water [J. Chem. Phys. 137, 014510 (2012)], a model parameterized using ab-initio calculations only. The model incorporates two features not present in many previously-studied water models: it is both flexible and polarizable, properties which can significantly influence the phase behavior of water. The observation of the critical point in a model in which the water-water interaction is estimated using only quantum ab-initio calculations provides strong support to the viewpoint according to which the existence of two distinct liquids is a robust feature in the free energy landscape of supercooled water.

Evidence of a liquid-liquid phase transition in H[Formula: see text]O and D[Formula: see text]O from path-integral molecular dynamics simulations

We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H[Formula: see text]O and D[Formula: see text]O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density [Formula: see text], isothermal compressibility [Formula: see text], and self-diffusion coefficients D(T) of H[Formula: see text]O and D[Formula: see text]O are in excellent agreement with available experimental data; the isobaric heat capacity [Formula: see text] obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H[Formula: see text]O and D[Formula: see text]O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H[Formula: see text]O and D[Formula: see text]O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H[Formula: see text]O, from PIMD simulations, is located at [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D[Formula: see text]O is estimated to be [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of [Formula: see text] for D[Formula: see text]O and, particularly, H[Formula: see text]O suggest that improved water models are needed for the study of supercooled water.

Temperature-Controlled Dual-Beam Optical Trap for Single Particle Studies of Organic Aerosol

An optical trapping cell that is capable of suspending particles using two counter-propagating beams in a temperature-controlled environment is reported here. With this dual-beam optical trap, we are able to hold single micron-sized droplets at temperatures down to 253 K (-20 °C) for hours at a time and in metastable (supercooled) states. As particles are trapped at the shared focal points of two intense beams, strong cavity-enhanced Raman scattering (CERS) is observed and allows for high precision measurements of physical properties. Here, the evaporation of highly oxygenated organic systems was monitored using CERS and was used to determine temperature-dependent vapor pressures and enthalpies of vaporization. The wavelength- and temperature-dependent optical properties were also simultaneously retrieved using CERS.

Identification of an Intermediate Species along the Nitrile Hydratase Reaction Pathway by EPR Spectroscopy

A new method to trap catalytic intermediate species was employed with Fe-type nitrile hydratase from Rhodococcus equi TG328-2 (ReNHase). ReNHase was incubated with substrates in a 23% (w/w) NaCl/H₂O eutectic system that remained liquid at -20 °C, thereby permitting the observation of transient species that were present at electron paramagnetic resonance (EPR)-detectable levels in samples frozen while in the steady state. Fe^III-EPR signals from the resting enzyme were unaffected by the presence of 23% NaCl, and the catalytic activity was ∼55% that in the absence of NaCl at the optimum pH of 7.5. The reaction of ReNHase in the eutectic system at -20 °C with the substrates acetonitrile or benzonitrile induced significant changes in the EPR spectra. A previously unobserved signal with highly rhombic g-values (g₁ = 2.31) was observed during the steady state but did not persist beyond the exhaustion of the substrate, indicating that it arises from a catalytically competent intermediate. Distinct signals due to product complexes provide a detailed mechanism for product release, the rate-limiting step of the reaction. Assignment of the observed EPR signals was facilitated by density functional theory calculations, which provided candidate structures and g-values for various proposed ReNHase intermediates. Collectively, these results provide new insights into the catalytic mechanism of NHase and offer a new approach for isolating and characterizing EPR-active intermediates in metalloenzymes.

Phase classification using neural networks: application to supercooled, polymorphic core-softened mixtures

Characterization of phases of soft matter systems is a challenge faced in many physical chemical problems. For polymorphic fluids it is an even greater challenge. Specifically, glass forming fluids, as water, can have, besides solid polymorphism, more than one liquid and glassy phases, and even a liquid-liquid critical point. In this sense, we apply a neural network algorithm to analyze the phase behavior of a mixture of core-softened fluids that interact through the continuous-shouldered well (CSW) potential, which have liquid polymorphism and liquid-liquid critical points, similar to water. We also apply the neural network to mixtures of CSW fluids and core-softened alcohols models. We combine and expand methods based on bond-orientational order parameters to study mixtures, applied to mixtures of hardcore fluids and to supercooled water, to include longer range coordination shells. With this, the trained neural network was able to properly predict the crystalline solid phases, the fluid phases and the amorphous phase for the pure CSW and CSW-alcohols mixtures with high efficiency. More than this, information about the phase populations, obtained from the network approach, can help verify if the phase transition is continuous or discontinuous, and also to interpret how the metastable amorphous region spreads along the stable high density fluid phase. These findings help to understand the behavior of supercooled polymorphic fluids and extend the comprehension of how amphiphilic solutes affect the phases behavior.

Thermodynamic anomalies in silicon and the relationship to the phase diagram

The evolution of thermodynamic anomalies are investigated in the pressure-temperature (pT) plane for silicon using the well-established Stillinger-Weber potential. Anomalies are observed in the density, compressibility and heat capacity. The relationships between them and with the liquid stability limit are investigated and related to the known thermodynamic constraints. The investigations are extended into the deeply supercooled regime using replica exchange techniques. Thermodynamic arguments are presented to justify the extension to low temperature, although a region of phase space is found to remain inaccessible due to unsuppressible crystallisation. The locus corresponding to the temperature of minimum compressibility is shown to display a characteristic 'S'-shape in thepTprojection which appears correlated with the underlying crystalline phase diagram. The progression of the anomalies is compared to the known underlying phase diagrams for both the crystal/liquid and amorphous/liquid states. The locations of the anomalies are also compared to those obtained from previous simulation work and (limited) experimental observations.

Core-softened water-alcohol mixtures: the solute-size effects

Water is the most anomalous material on Earth, with a long list of thermodynamic, dynamic and structural behaviors that deviate from what is expected. Recent studies have indicated that these anomalies may be related to a competition between two liquids, which means that water has a potential liquid-liquid phase transition (LLPT) that ends at a liquid-liquid critical point (LLCP). In a recent study [J. Mol. Liq., 2020, 320, 114420], using molecular dynamics simulations and a core-softened potential approach, we have shown that adding a simple solute such as methanol can "kill" the density-anomalous behavior as the LLCP is suppressed by spontaneous crystallization in a hexagonal close packing (HCP) crystal near the LLPT. Now, we extend this work to realize how longer-chain alcohols will affect the complex behavior of water-alcohol mixtures in the supercooled regime. Besides core-softened (CS) methanol, ethanol and 1-propanol were added to a system of identical particles that interact through the continuous shouldered well (CSW) potential. We observed that the density anomaly gradually decreases its extension in phase diagrams until it disappears with the growth of the non-polar chain and the alcohol concentration, different from the liquid-liquid phase transition (and the LLCP), which remained present in all analyzed mixtures, according to Nature, 2001, 409, 692. For our model, the longer non-polar chains and higher concentrations gradually impact the competition between the scales in the CS potential, leading to a gradual disappearing of the anomalies until the TMD total disappearance is observed when the first coordination shell structure is also affected: short-range ordering is favored, leading to less competition between short- and long-range ordering and, consequently, to the extinction of anomalies. Also, the non-polar chain size and concentration have an effect on the solid phases, favoring the hexagonal close packed (HCP) solid and the amorphous solid phase over the body-centered cubic (BCC) crystal. These findings help to elucidate the behavior of water solutions in the supercooled regime and indicate that the LLCP can be observed in systems without density-anomalous behavior.

On the existence of soliton-like collective modes in liquid water at the viscoelastic crossover

The problem of large-density variations in supercooled and ambient water has been widely discussed in the past years. Recent studies have indicated the possibility of nanometer-sized density variations on the subpicosecond and picosecond time scales. The nature of fluctuating density heterogeneities remains a highly debated issue. In the present work, we address the problem of possible association of such density variations with the dynamics of terahertz longitudinal acoustic-like modes in liquid water. Our study is based on the fact that the subpicosecond dynamics of liquid water are essentially governed by the structural relaxation. Using a mode coupling theory approach, we found that for typical values of parameters of liquid water, the dynamic mechanism coming from the combination of the structural relaxation process and the finiteness of the amplitude of terahertz longitudinal acoustic-like mode gives rise to a soliton-like collective mode on a temperature-dependent nanometer length scale. The characteristics of this mode are consistent with the estimates of the amplitudes and temperature-dependent correlation lengths of density fluctuations in liquid water obtained in experiments and simulations. Thus, the fully dynamic mechanism could contribute to the formation and dynamics of fluctuating density heterogeneities. The soliton-like collective excitations suggested by our analysis may be relevant to different phenomena connected with supercooled water and can be expected to be associated with some ultrafast biological processes.

Long-Range Orientational Organization of Dipolar and Steric Liquids

Long-range orientational correlations in liquids have received recent renewed interest, in particular for the neat water case. These long-range orientational correlations, exceeding several tens of nanometers, originate from the presence of the strong permanent water dipolar moment. However, the exact dependence with the dipolar moment and the role of other local forces like steric hindrance has never been addressed. In this work, we experimentally measure long-range correlations for a set of liquids differing by their molecular weight and dipolar moment, in order to reveal the origin of their long-range organization. Hence, we show that the dipolar moment of a solvent molecule is not the unique feature determining the orientational correlation. Steric hindrance significantly helps to structure the liquids as well. In order to quantify these long-range correlations, we also derive theoretically the polarization resolved second harmonic scattering intensity as a function of the rotational invariants describing the dipolar and octupolar interaction.

Could Chemical Reaction at the Molecular Level Show Distinction between Two Liquid-Water States? Study of the Excited-State Water-Catalyzed Proton Transfer Reaction Provides a Clue

The two liquid-water states, which lead to some anomalies when temperature crosses over 50 ± 10 °C at the atmospheric pressure, have been continuously catching popular attention. In this study, using the excited-state proton transfer (ESPT) catalyzed by water molecules as a prototypical reaction, we demonstrate that the kinetics of ESPT indeed is influenced by the two liquid-water states. In the water-catalyzed ESPT of 3-cyano-7-azaindole (3CAI), a repetitive and comprehensive temperature-dependent study of ESPT in H₂O from 0 to 90 °C shows anomalous behavior. The plot of the logarithm of ESPT rate constant as a function of inverse of absolute temperature deviates from a straight line. The convex-Arrhenius behavior manifests the activation free energy for water-assisted ESPT being dependent on temperature and hence the liquid water structure. To simplify the discussion, the plot is well fitted by using two straight lines that are crossed over in the vicinity of 40 °C. The free energy difference between water-solvated 3CAI and the 1:1 H₂O:3CAI complex is deduced to be 2.29 ± 0.04 and 1.96 ± 0.04 kcal·mol^-1 in the regions of 0-40 and 40-90 °C water, respectively, which also results in different frequency factors, i.e., the proton transfer/tunneling rates of (5.83 ± 0.36) × 10¹⁰ and (3.48 ± 0.27) × 10¹⁰ s^-1, respectively. In a qualitative manner, the results are then rationalized by the different types of H-bonding configuration as proposed for two liquid-water phases, rendering experimental evidence to support the different water phases in ambient temperatures at 1 bar.

Phase Change Behavior of NaCl-H2O Binary Solution Under the Control of AC Electric Field

Cryopreservation, which refers to preservation of cells or tissues at subzero temperatures, inevitably involves the problem of cryoinjury caused by ice crystals. The application of an external electric field during the freezing process has been shown to be a promising approach to produce miniature ice grains and decrease the fraction of ice crystallization at a slow cooling rate. Thus, the dielectric and thermodynamic properties of NaCl-H₂O binary solutions at subzero temperatures were tremendously important for understanding the mechanism of ice formation under the manipulation of an AC electric field in biopreservation. However, there was still a lack of relevant information in the literature. The first objective of this study was to systematically measure the dielectric spectrum of 0.9% NaCl-H₂O binary solutions at temperatures ranging from -100°C to 0°C with a cooling/heating rate of 2°C/min. We further measured the thermodynamic properties of a 0.9% NaCl-H₂O binary solution while applying a series of electric fields near its dielectric relaxation frequency. The effect of the electric field on the crystal morphology was studied last. Pure water was selected as the control group. The results showed that an AC electric field can alter the thermodynamic process and thus the phase transition and ice crystal structure could be manipulated. It was concluded that the AC electric-assistant preservation method will be a promising technology in cryopreservation.

Rheology of water in small nanotubes

The properties of water in confinement are very different from those under bulk conditions. In some cases the melting point of ice may be shifted and one may find either ice, icelike water, or a state in which freezing is completely inhibited. Understanding the dynamics and rheology of water in confined media, such as small nanotubes, is of fundamental importance to the biological properties of micro-organisms at low temperatures, to the development of new devices for preserving DNA samples, and for other biological materials and fluids, lubrication, and development of nanostructured materials. We study rheology and dynamics of water in small nanotubes using extensive equilibrium and nonequilibrium molecular dynamics simulations. The results demonstrate that in strong confinement in nanotubes at temperatures significantly below and above bulk freezing temperature water behaves as a shear-thinning fluid at shear rates smaller than the inverse of the relaxation time in the confined medium. In addition, our results indicate the presence of regions in which the local density of water varies significantly over the same range of temperature in the nanotube. These findings may also have important implications for the design of nanofluidic systems.

Analytical Determination of the Nucleation-Prone, Low-Density Fraction of Subcooled Water

Subcooled water is the primordial matrix for ice embryo formation by homogeneous and heterogeneous nucleation. The knowledge of the specific Gibbs free energy and other thermodynamic quantities of subcooled water is one of the basic prerequisites of the theoretical analysis of ice crystallization in terms of classical nucleation theory. The most advanced equation of state of subcooled water is the IAPWS G12-15 formulation. The determination of the thermodynamic quantities of subcooled water on the basis of this equation of state requires the iterative determination of the fraction of low-density water in the two-state mixture of low-density and high-density subcooled water from a transcendental equation. For applications such as microscopic nucleation simulation models requiring highly frequent calls of the IAPWS G12-15 calculus, a new two-step predictor-corrector method for the approximative determination of the low-density water fraction has been developed. The new solution method allows a sufficiently accurate determination of the specific Gibbs energy and of all other thermodynamic quantities of subcooled water at given pressure and temperature, such as specific volume and mass density, specific entropy, isothermal compressibility, thermal expansion coefficient, specific isobaric and isochoric heat capacities, and speed of sound. The misfit of this new approximate analytical solution against the exact numerical solution was demonstrated to be smaller than or equal to the misprediction of the original IAPWS G12-15 formulation with respect to experimental values.

Experimental observation of density fluctuations in liquid metals associated with liquid-liquid, liquid-gas and metal-nonmetal transitions

We have developed a special technique and succeeded to carry out small-angle x-ray scattering measurements for some liquid metal systems. The purpose is to investigate effects of transitions such as liquid-liquid (LLT), liquid-gas (LGT) and metal-nonmetal (MNMT) transitions on mesoscopic density fluctuations in liquids. In liquid Te systems (Se-Te and Ge-Te mixtures), which show continuous LLT accompanying MNMT, parameters of density fluctuations show maxima almost in the middle of the transition, both in strength and spatial size. This work (and Kajihara et al 2012 Phys. Rev. B86 214202) was the first direct observation that density fluctuations exhibit maximum corresponding to LLT. However in this study, we could not clearly separate the effects of LLT and MNMT on the observed density fluctuations. Thus, we also investigated fluid Hg under high pressure and high temperature conditions, which shows MNMT near a critical point of LGT, to investigate how MNMT affects them. We observed distinct density fluctuations; a strength and a correlation length of them show maxima at around a critical isochore of LGT, and the former is basically consistent with a phase diagram (compressibility) of LGT; they do not show any peaks at MNMT region. Precise analysis revealed that MNMT only affects a shift of another parameter, a short-range correlation length. These results in fluid Hg indicate that the density fluctuations are mainly derived from a critical phenomena of LGT and MNMT does not play any critical role on them. We believe that the latter conclusion also holds true for liquid Te systems; MNMT plays no important role on the density fluctuations in liquid Te systems and LLT is the main origin of them.

Ice-Crystal Nucleation in Water: Thermodynamic Driving Force and Surface Tension. Part I: Theoretical Foundation

A recently developed thermodynamic theory for the determination of the driving force of crystallization and the crystal–melt surface tension is applied to the ice-water system employing the new Thermodynamic Equation of Seawater TEOS-10. The deviations of approximative formulations of the driving force and the surface tension from the exact reference properties are quantified, showing that the proposed simplifications are applicable for low to moderate undercooling and pressure differences to the respective equilibrium state of water. The TEOS-10-based predictions of the ice crystallization rate revealed pressure-induced deceleration of ice nucleation with an increasing pressure, and acceleration of ice nucleation by pressure decrease. This result is in, at least, qualitative agreement with laboratory experiments and computer simulations. Both the temperature and pressure dependencies of the ice-water surface tension were found to be in line with the le Chatelier–Braun principle, in that the surface tension decreases upon increasing degree of metastability of water (by decreasing temperature and pressure), which favors nucleation to move the system back to a stable state. The reason for this behavior is discussed. Finally, the Kauzmann temperature of the ice-water system was found to amount TK=116K, which is far below the temperature of homogeneous freezing. The Kauzmann pressure was found to amount to pK=−212MPa, suggesting favor of homogeneous freezing on exerting a negative pressure on the liquid. In terms of thermodynamic properties entering the theory, the reason for the negative Kauzmann pressure is the higher mass density of water in comparison to ice at the melting point.

Understanding the Origin of the Breakdown of the Stokes-Einstein Relation in Supercooled Water at Different Temperature-Pressure Conditions

A recent experiment has measured the viscosity of water down to approximately 244 K and up to 300 MPa. The correct viscosity and translational diffusivity data at various temperature-pressure (T-P) state points allowed for checking the validity of the Stokes-Einstein (SE) relation, which accounts for the coupling between translational self-diffusion and medium viscosity. The diffusion-viscosity decoupling increases with decreasing temperature, but the increasing pressure reduces the extent of the decoupling. Earlier simulation studies explained the breakdown of the SE relation in terms of the location of the Widom line, emanating from the liquid-liquid critical point (LLCP). Although these studies made a significant contribution to the current understanding of the above phenomena, a detailed molecular picture is still lacking. Recently, our group has explained the diffusion-viscosity decoupling from a jump-diffusion perspective. The jump-diffusion coefficient, emanating from the jump translation of water molecules, is calculated using a quantitative approach for different temperatures at ambient pressure. It has been observed that jump-diffusion is the key factor for diffusion-viscosity decoupling in supercooled water. The same method is adopted in the present work to estimate the jump-diffusion coefficient for different T-P state points and, thereby, explains the role of jump-diffusion for the different extents of the SE relation breakdown at different pressures. The residual diffusion coefficient, the other component of the total diffusion that originates from small step displacement and that is calculated by subtracting the jump-diffusion coefficient from the total diffusion, is seen to be fairly coupled to the viscosity at the entire range of temperature and pressure. Furthermore, we have calculated the average number of H-bonds per water molecule and the tetrahedral order for different T-P state points and investigated an approximate correlation between the average local structure and the contribution of the jump-diffusion to the total diffusion of water. This study, therefore, puts forward a new perspective for explaining the SE relation breakdown in supercooled water under different pressure conditions.

Thermodynamics of supercooled and stretched water: Unifying two-structure description and liquid-vapor spinodal

We have applied a two-structure approach to the description of the thermodynamic properties of supercooled and stretched water, metastable toward vapor, ice, or both, by incorporating the stability limit of liquid with respect to vapor at negative pressures. In addition to the properties of water considered in previous studies, we include new data recently obtained in deeply supercooled and stretched regions. Our model reproduces the experimentally observed anomalies in metastable water up to 400 MPa and down to -140 MPa, and can provide a physically based extrapolation in regions where no measurements are available yet. Moreover, we are able to elucidate the thermodynamic nature of the alternative "states" of liquid water, namely, high-temperature denser water (state A) and "mother-of-ice" lighter water (state B). Based on the internal consistency of the described anomalies and new data on the isothermal compressibility, we exclude the critical-point-free scenario in which the first-order liquid-liquid transition line would continue into the stretched liquid state (doubly metastable) crossing the vapor-liquid spinodal. A "singularity-free" scenario remains an option for explaining supercooled water's anomalies within the framework of two-state thermodynamics; however, the extreme case of the singularity-free scenario, ideal mixing of A and B, seems improbable. We have also clarified the concept of fast interconversion of alternative states in supercooled water as a phenomenological representation of distribution of short-ranged local structures.

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.

A polymerization scenario of the liquid-liquid transition of GeI4

SnI₄ and GeI₄ have been shown to exhibit similar polyamorphic nature. We examine the microscopic nature of the liquid-liquid transition in GeI₄ by conducting an isothermal-isobaric molecular dynamics simulation for the system composed of rigid tetrahedral molecules. The model allows us to semiquantitatively discuss the structural properties of liquid GeI₄ below 1 GPa. We define a physical bond between the nearest intermolecular iodine sites satisfying the conditions of forming the metallic I₂ bond. We then focus on the formation of molecular clusters in dynamic networks of the bonds. The clusters are mainly formed by molecules whose nearest pairs are in edge-to-edge, face-to-edge, and vertex-to-edge orientations. The clusters grow as pressure increases, and the onset of percolation is observed below 1 GPa. The finite-size scaling analysis for the percolation probability identifies that the threshold pressure is [Formula: see text] GPa for the present model, which is near the extension of the boundary between the two liquid phases. We thus speculate that the liquid-liquid transition of GeI₄ is attained by polymerization, i.e. percolation of molecular networks. The same percolation scenario with a slight modification such as introducing a bootstrap mechanism is expected to apply to the transition in liquid SnI₄.

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.

Slow dynamics coupled with cluster formation in ultrasoft-potential glasses

We numerically investigate the slow dynamics of a binary mixture of ultrasoft particles interacting with the generalized Hertzian potential. If the softness parameter, α, is small, the particles at high densities start penetrating each other, form clusters, and eventually undergo the glass transition. We find multiple cluster-glass phases characterized by a different number of particles per cluster, whose boundary lines are sharply separated by the cluster size. Anomalous logarithmic slow relaxation of the density correlation functions is observed in the vicinity of these glass-glass phase boundaries, which hints the existence of the higher-order dynamical singularities predicted by the mode-coupling theory. Deeply in the cluster glass phases, it is found that the dynamics of a single particle is decoupled from that of the collective fluctuations.

Second inflection point of water surface tension in the deeply supercooled regime revealed by entropy anomaly and surface structure using molecular dynamics simulations

The surface tension of supercooled water is of fundamental importance in physical chemistry and materials and atmospheric sciences. Controversy, however, exists over its temperature dependence in the supercooled regime, especially on the existence of the "second inflection point (SIP)". Here, we use molecular dynamics simulations of the SPC/E water model to study the surface tension of water (σw) as a function of temperature down to 198.15 K, and find a minimum point of surface excess entropy per unit area around ∼240-250 K. Additional simulations with the TIP4P/2005 water model also show consistent results. Hence, we predict an SIP of σw roughly in this region, at the boundary where the "no man's land" happens. The increase of surface entropy with decreasing temperature in the region below the inflection point is clearly an anomalous behavior, unknown for simple liquids. Furthermore, we find that σw has a near-linear correlation with the interfacial width, which can be well explained by the capillary wave theory. Deep in the supercooled regime, a compact water layer at the interface is detected in our simulations, which may be a key component that contributes to the deviation of surface tension from the International Association for the Properties of Water and Steam relationship. Our findings may advance the understanding of the origin of the anomalous properties of liquid water in the supercooled regime.

Water Is a Cagey Liquid

Liquid water is considered poorly understood. How are water's physical properties encoded in its molecular structure? We introduce a statistical mechanical model (CageWater) of water's hydrogen-bonding (HB) and Lennard-Jones (LJ) interactions. It predicts the energetic and volumetric and anomalous properties accurately. Yet, because the model is analytical, it is essentially instantaneous to compute. This model advances our understanding beyond current molecular simulations and experiments. Water has long been regarded as a "2-density liquid": a dense LJ liquid and a looser HB one. Instead, we find here a different antagonism underlying water structure-property relations: HBs in water-water pairs drive density, while HBs in cooperative cages drive openness. The balance shifts strongly with temperature and pressure. This model interprets the molecular structures underlying the liquid-liquid phase transition in supercooled water. It may have value in geophysics, biomolecular modeling, and engineering of materials for water purification and green chemistry.