New Considerations Around the Singular Water Temperature Explaining the Anomalous Behavior of Liquid Phase

TL;DR

The paper explores how water's unique hydrogen bonding leads to unusual behavior at a specific temperature around 315 K, linked to its structural phases.

Contribution

The study proposes a universal temperature where water's thermodynamic properties change, explained by the coexistence of two liquid phases.

Findings

A singular temperature T* ≃ 315±5 K is identified where water's isothermal compressibility and expansivity show anomalies.

The behavior at T* is attributed to the coexistence of low-density and high-density liquid phases of water.

The average intermolecular distance supports the existence of two distinct liquid phases in supercooled water.

Abstract

The water thermodynamics is characterized by polydispersity, which determines its structural and dynamic properties. This is due to the specifics of its characteristic bond: the hydrogen bond (HB). The isobars of the two fundamental thermodynamic functions, the isothermal compressibility (KT(P.T)) and the isobaric expansivity (αP(P,T)), show the presence of a temperature T*≃315±5 K where both have a singular behavior. In this work, by carefully considering the thermal properties of the isobars of density ρ, specific heat CP and the self-diffusion DS, we suggest the universality characteristics of this temperature. In addition, by analyzing the average intermolecular distance dOO, in the same area of the P-T phase diagram, we demonstrate that such realities are due in the supercooled liquid state to the ratio between its two characteristic phases: the low-density liquid (LDL due to HB) and the HDL (which entirely characterizes the remaining parts of the phase diagram).

Figures

Full text

1. Introduction

As is well known, water shows chemical and physical diversity when compared to other liquids, so it is not yet fully rationalized. Although it has a simple molecular structure, water has different chemical–physical characteristics than other substances [1]. For these intriguing reasons, this liquid has always been the focus of various careful studies. Furthermore, if the corresponding studies are of interest, in science and engineering, water can be explained only by the use of physico-chemical models and experiments, in particular that of the thermodynamic and of the modern statistical physics [2]. Starting from the well-known density maximum, observed under the ambient pressure (0.1 MPa) at 4 °C, unusual properties of the liquid phase of water include the minima in the T dependences of the isothermal compressibility ( $[eqn]$ ) and isobaric specific heat ( $[eqn]$ ) located at about 40 °C. Similar behaviors are also observable in the water thermodynamic functions like the self-diffusion coefficient ( $[eqn]$ ) that near the ambient temperature increases its growth in response to an applied pressure (unlike the corresponding density ( $[eqn]$ )).

Like some other liquids, water can be easily supercooled below its melting temperature ( $[eqn]$ K), where it shows unusual dynamics with characteristic divergent behaviors for which is as a sort of prototype of glass-forming systems. The limit of this supercooling region is the homogeneous nucleation temperature ( $[eqn]$ K for $[eqn]$ MPa) at which the nucleation rate suddenly becomes very large. However, such a limit, being of a kinetic nature, is not absolute and can be bypassed if the observation time is shorter than the nucleation time. At lower temperatures, a highly viscous liquid phase (also metastable) can be obtained at about 150 K ( $[eqn]$ ) by heating water to the glassy state; liquid that crystallizes as cubic ice when T is further reduced. In such a situation, extremely high viscosity causes the nucleation rate to slow down, allowing much longer observation times. So $[eqn]$ and $[eqn]$ are then the limits of a region inaccessible to experiments on the liquid in its metastable phase. These limits can be overcome by using confined water, water in man-sized pores, or emulsified water. By choosing the confinement pore size smaller than its homogeneous nucleation process, water does not crystallize and can be further cooled below $[eqn]$ .

Thus, by using confining sizes of the order of tens of micrometers, it was possible to extend the experimental studies of the bulk water thermodynamics in a large region of its P-T phase diagram. The pioneer of this idea was Angell with collaborators, who detailed the behaviors of $[eqn]$ and $[eqn]$ in the region $[eqn]$ . A goal of significant interest in the obtained results is that by showing specific power laws they proposed the presence of a critical behavior [3].

Water in the solid crystalline phase posses another relevant characteristic: the polymorphism, characterized by many ice structures ranging from the ice Ic to ice XII (all of different density) [4].

Mischima studies on amorphous water, clearly highlighting also the glass polymorphism, have provided new information on liquid properties [5,6]. IIn fact, a real polyamorphism was discovered in the system, with the transition governed by the T, between the two amorphous forms of different densities, i.e., from high (HDA) to low (LDA). The latter, with $[eqn]$ gcm^−3^ similar to that of the ice I_h, was discovered in 1935 [7], whereas HDA can be achieved, at 77 K, up to 2 GPa, by compressing ice Ih_, and it has a density of ∼1.17 gcm^−3^ [8]. In addition, if heated, at $[eqn]$ MPa, it passes into the low-density form at ∼120 K. Such transition is reversible (HDA ⟷ LDA) as is the one between LDA and HDA by pressure-cycling at about 135 K and 2 GPa) [9]. As said, the LDA, when heated, transforms into highly viscous liquid.

An additional higher-density glass form has recently been discovered [10]: the VHDA or very-high-density amorphous [11]. The extension of this “polyamorphism”, with these singular amorphous forms, in the liquid phase has allowed us to clarify the system thermodynamics. Like in glass, the liquid polymorphism is characterized by two phases different in densities: the high- and low-density liquids, HDL and LDL [12]. While the latter is more relevant from a physical point of view, for its “open” structure originated by a clustering with a tetrahedral symmetry due to the non-covalent attractive hydrogen bonding (HB) interaction, the former is instead made up of a rather homogeneous distribution of trimers, dimers, or single water molecules. These two liquids coexist and depending on the thermodynamic variables can transform into each other through a first-order transition: the liquid–liquid transition hypothesis (LLT) [13,14,15,16,17]. All these water experiments are proven realities: polyamorphism, the continuity between its glass and liquid accompanied by a discontinuity between the two liquid forms, HDL and LDL (the liquid–liquid transition, LLT) are at the bases of the liquid–liquid critical point (LLCP) hypothesis (an additional critical point in addition to the vapor–liquid one) [18]. Three interactions, one repulsive and two attractive, govern the water properties: the first is the purely electrical Coulomb force between electron lone pairs of adjacent oxygen atoms; the attractive ones are the covalent bond CB (O-H) and the hydrogen bond HB (O-H).

The single molecule (H_2_O) originates from two CBs sharing the electron pairs with a binding energy BE $[eqn]$ eV and filling up other two orbitals by means of its non-bonding lone pairs originating the intermolecular non-covalent, O:H, van der Waals bond ( $[eqn]$ eV). Three distances are representative of its structure (two intermolecular $[eqn]$ and $[eqn]$ and the intramolecular dO-H) and two angles (the internal H-O-H ( $[eqn]$ ) and intermolecular H:O-H ( $[eqn]$ ), depending on both pressure and temperature. In general, the Bernal and Fowler ice rule holds and the oxygen atom always tends to find four neighbors to form a stable tetrahedron ( $[eqn]$ , $[eqn]$ ) whereas, on the contrary, in the liquid phase, the repulsion between electron pairs on oxygen counteract the steady tetrahedron formation. In water, such a basic building block contains two equivalent molecules and four identical O:H-O bonds (with different orientations). Despite the fluctuations in $[eqn]$ and $[eqn]$ , its average structural dimensions are independent of the thermodynamic transitions that govern water, from solid to liquid phase or to glass or vapor one. Summarizing, while HB dominates the liquid in both the stable and metastable regimes (supercooled), the repulsion mainly influences the water properties from above the boiling temperature ( $[eqn]$ ) up to the sub-critical and critical region. Finally, the vapor–liquid critical point CP coordinates are $[eqn]$ K, $[eqn]$ MPa, whereas the LLCP is estimated by Molecular Dynamics simulations to be at 200 K and $[eqn]$ MPa [19].

The observed behaviors of some thermodynamic functions of water, such as their divergences, and the decrease in entropy accompanied by polymorphism (especially that of the liquid phase) are at the origin of its complexity and related anomalies. It should be noted, however, that through these behaviors, the presence of local tetrahedral order has been experimentally demonstrated. Unfortunately, in spite of the proper theoretical suggestions and the very large number of accurate computational studies accompanied by proper suggestions [19], the suggested liquid criticality in the metastable (supercooled) phase is far from being proven [20]. Today, such a result is an intriguing chimera for experimental physics.

Despite this, we still have some certainties observed in the liquid water when it is cooled, such as the growth in size and stability of the tetrabonded clusters. Conditions for which the HB lifetime strongly increases by many orders of magnitude from the some picoseconds typical of the stable liquid region ( $[eqn]$ ) [21]. In contrast to the permanent ice tetrahedral network, the liquid tetrahedrality is local and transient. An increase in pressure opposes these ordering processes.

As said, bulk liquid water cannot exist in the region between the homogeneous nucleation temperature ( $[eqn]$ ) and that in which the ultra-viscous liquid crystallizes ( $[eqn]$ ): the so-called “No Man’s Land” water. But such a constraint can be overcome by using confined water so that water can be maintained in the liquid state inside the range $[eqn]$ – $[eqn]$ and thus explored [22].

Other ways to work inside this region is to use water in solutions, that of hydration of macromolecules and in micellar systems [23], nano droplets [24] inside and outside ice, or by melting a multi-molecular thickness of an ice surface [25]. Just these approaches have led to the discovery of new and important properties that characterize it, like the density minimum [26,27], predicted by Percy W. Bridgman [28] and confirmed by proper computational studies [29,30,31]. Regarding these reported MD data, we stress that they come from different models (TIP4P or E3B3) largely used for the water system [19].

Other important findings concern the dynamics: the crossover from a fragile to a strong glass-forming material, predicted at ambient pressure [32] and observed at $[eqn]$ K [33]; in the P-T phase diagram, this is also the locus of the Stokes–Einstein relation violation (due essentially to the decoupling between translational and rotational modes) and of the Widom line that, strongly linked to the LLCP hypothesis, identifies the maximum in the $[eqn]$ and $[eqn]$ fluctuations where thermodynamic response functions reach their extremes (minimum with negative values in the $[eqn]$ and maxima in $[eqn]$ and $[eqn]$ ).

Here, by using the liquid and amorphous water density literature data, reported in Figure 1 in a wide P-T range, we focus on the HB role and effects (structural and dynamical). Previous studies (see, e.g., [20]) have detailed that the LDL is present (together with HDL) only in a limited phase diagram area, where it is responsible of the LLT and the corresponding anomalies, whereas in the remaining parts, essentially only the HDL exists. Also from the reported data, the existence of a limiting pressure is evident above which the liquid is present only as HDL and freezes only into HDA, whereas at lower pressures, where HDL and LDL phases coexist, the liquid evolves only toward the LDA. Using these $[eqn]$ data, we evaluate the T—dependences of two fundamental thermodynamic functions, the isothermal compressibility ( $[eqn]$ ) and the isobaric expansivity ( $[eqn]$ )), at different isobars discovering the presence of a temperature $[eqn]$ K in which both show a singular behavior: all the compressibility isobars have a minimum value and the corresponding $[eqn]$ functions cross each other [34]. The “singular and universal expansivity point” is $[eqn]$ K^−1^.

Furthermore, according to the results of NMR studies, it has been highlighted that this temperature defines the boundary point between two regions of different energy configurations and therefore different dynamics. Specifically, for $[eqn]$ the water properties are governed by two defined energy levels (the Arrhenius behavior) while in opposite thermal conditions ( $[eqn]$ ) the situation changes completely, becoming that which characterizes supercooled glass-forming liquid systems. Below $[eqn]$ , cooling gives rise in water to specific correlations in the time and length scale, resulting in a dynamical clustering, i.e., the onset of the HB tetrahedral network and its effects [35]. Starting from these results, in this work, we aim to demonstrate the universality of this temperature in characterizing the physics of water in its structural and dynamic properties.

2. Results and Discussion

2.1. The Density

Figure 1 illustrates the bulk water densities ( $[eqn]$ ) for 1 < P (MPa) $[eqn]$ and $[eqn]$ T (K) $[eqn]$ K [36,37,38,39,40,41,42]; the values measured in the LDA and HDA glass are also reported, as well as the data of a MD simulation [29,30] together with the data measured at ambient pressure ( $[eqn]$ MPa) for water under MCM nanotubes confinement and emulsifieds (droplet size 1–10 μm) [42].

Being P-dependent, the density maximum disappears for $[eqn]$ MPa. Such pressure also identifies the change in the curvature of the density isobars ( $[eqn]$ ) from negative to positive. The densities measured at 130 and 150 K in the HDA and LDA are also shown. Whereas the HDA data measured at very high pressures are of the order of $[eqn]$ gcm^−3^ (or even higher), the value of the corresponding LDA for $[eqn]$ MPa and 135 K is ∼0.94 gcm^−3^. The latter value agrees with the values measured in confined water (nanotubes) inside the range $[eqn]$ – $[eqn]$ , where a minimum is present at $[eqn]$ K [26].

From the reported data, it can be observed that (i) the $[eqn]$ grows with pressure and such increase is larger in the metastable supercooled phase; (ii) the data spread at a fixed temperature is lager just in the region of the two amorphous phases (e.g., at 130 K is $[eqn]$ ) rather than in the heated liquid (at 400 K is $[eqn]$ 0.22; (iii) as said, for P ≥ 200 MPa, the curvature changes from concave to convex with an evolution similar to that of a simple fluid; (iv) the pressure of 200 MPa defines two different behaviors: all the isobars below such pressure density evolves toward the LDA values, while all those above it proceed to the HDA; (v) whereas the LDA densities are almost P independent, those of the HDA phase show a marked dependence on this variable.

2.2. The Thermodynamic Density Derivatives (The Isothermal

Compressibility and the Isobaric Thermal Expansivity)

Figure 2 shows, from 200 to 450 K, the different isobars ( $[eqn]$ to 400 MPa) of the isothermal compressibility, K_T_(T,P), obtained as $[eqn]$ , a function that also represents the volume fluctuations $[eqn]$ as $[eqn]$ . It is clearly observable that all the obtained compressibility isobars show a minimum just at $[eqn]$ . In addition, an increase in the corresponding values for $[eqn]$ can be observed for all isobars compared to the second region at higher temperatures. This growth is increasingly marked in the metastability (supercooled) regions, especially for the isobars from $[eqn]$ to 200 MPa, while at higher pressures ( $[eqn]$ MPa) the resultig compressibility values are much lower and almost symmetric in temperature around the minimum in $[eqn]$ .

The same Figure proposes in the inset the isobaric expansivity $[eqn]$ , linked to the entropy–volume cross-correlations $[eqn]$ to be $[eqn]$ . This function is proposed for $[eqn]$ T(K) $[eqn]$ K and pressures from $[eqn]$ to 800 MPa.

Regarding $[eqn]$ , at ambient P in simple liquids, the $[eqn]$ and $[eqn]$ fluctuations become smaller on decreasing T and are positively correlated. In water, the situation is opposite: they are not only more pronounced, but in the region of metastability ( $[eqn]$ K) they are anti-correlated [43], a behavior occurring due to the evolution of the water local order in the supercooled region. Furthermore, the thermal evolution of both these functions, compressibility and $[eqn]$ , shows that $[eqn]$ is a border between two different behaviors. All the $[eqn]$ isobars cross, within the error, just at $[eqn]$ . Even in this case we see that above this temperature the thermodynamic behavior of water is exactly the same as a normal fluid, but the situation changes on cooling, where anomalous behaviors are evident. More precisely, with $[eqn]$ at 200, 300, and 400 MPa, a continuous evolution occurs at all temperatures. But at the lowest pressures ( $[eqn]$ MPa) for $[eqn]$ , the $[eqn]$ evolution becomes more complex, i.e., it decreases as T decreases and, as proposed by the $[eqn]$ MPa isobar, the entropy–volume anticorrelations determine an inflection point and its minimum. This minimum decreases in value as P increases, a behavior consistent with the density minimum observed in confined water. Note that in the case of confined water (MCM-41), such minimum temperature is coincident with that of the fragile-to-strong dynamic crossover and the Widom line, which at ambient pressure is $[eqn]$ K. We have to stress that the water singular temperature $[eqn]$ has a precise thermodynamical consistence lying in the relationship connecting two of the studied response functions:

$$[eqn]$$

a thermodynamic rule that confirms the $[eqn]$ role for the liquid water [34]. Just to show the uniqueness of this temperature and the marked difference in the isobar behavior above and below it, also by considering the suggestions of these two thermodynamic response functions, we normalized the data of the isobars reported in Figure 1 by subtracting for each one the corresponding density value $[eqn]$ . The obtained result ireported in Figure 3 shows another special feature: from $[eqn]$ up to $[eqn]$ , the thermal effect on for all the density isobars is the same (and linear).

2.3. The Constant Pressure-Specific Heat

Figure 4 shows, in the range of $[eqn]$ K, two sets of the constant pressure-specific heat isobars, $[eqn]$ .

The first one illustrates the experimental measurements obtained with a differential calorimeter [44,45] (at $[eqn]$ MPa) and sound velocity [46]; the second instead presents data obtained from NMR time relaxation [47] data through the use of the Adam–Gibbs relation [48]:

$$[eqn]$$

with the A constant. Data fitting in the range of bulk water gives $[eqn]$ × $[eqn]$ $[eqn]$ $[eqn]$ , whereas A is estimated to be $[eqn]$ kJ mol^−1^; these values are found to be in agreement with those of ref. [49]. The corresponding values of $[eqn]$ and $[eqn]$ are thus obtained. Hence, $[eqn]$ is obtained by means of the derivative $[eqn]$ . From the temperature behavior of all the $[eqn]$ isobars, a slope change can be observed at about $[eqn]$ . The data behavior is characterized by a maximum (e.g., for $[eqn]$ MPa at about 225 K) in the $[eqn]$ measured in confined water. A coincidence of the values can also be observed, within the experimental error, between those evaluated by the self-diffusion coefficients and those experimentally measured in liquid.

2.4. The Self-Diffusion Coefficient

We now examine the self-diffusion coefficient $[eqn]$ data in order to gain further information about $[eqn]$ . Figure 5 shows an Arrhenius plot of the bulk water $[eqn]$ measured in different experiments. Function of pressure ( $[eqn]$ P(MPa) $[eqn]$ ) is in the range $[eqn]$ T(K) $[eqn]$ .

Data are collected from Nuclear Magnetic Resonance (NMR) experiments [21,50,51,52,53,54,55,56]. For $[eqn]$ , the reported data evolve according to the Stokes–Einstein relation, whereas below we show a super-Arrhenius behavior according to the universal one corresponding to the materials liquid state [57]. In particular, for water, $[eqn]$ represents the crossover temperature marking the system dynamic crossover from simple two energetic states (Arrhenius or strong glass forming liquids [58]) to the super-Arrhenius (or fragile liquids) characterized by an energetic multi-state distribution, which in turn is characterized by dynamic heterogeneities (with a multi-relaxation in the time evolution of the density–density correlation functions) typical of the supercooled regime and directly linked with $[eqn]$ [59,60].

2.5. The Rotational Relaxation Time

The general validity of this universal-like crossover is confirmed by the rotational relaxation time data, $[eqn]$ , at different temperatures and pressures reported in an Arrhenius plot in Figure 6. Temperature data are considered from the boiling point down to 185 K for many pressures, as indicated in the legend from $[eqn]$ to 250 MPa.

In particular, data coming from an NMR experiment are illustrated, measured by means of the system longitudinal relaxation times $[eqn]$ in emulsioned water contained in a high-pressure glass capillary [47]. Hence, the $[eqn]$ values were evaluated by taking into account only the inter- and intramolecular dipolar contributions observed in the proton relaxation as $[eqn]$ according to proper procedure [61]. As can be easily observed also for $[eqn]$ , in accordance with what has been said, $[eqn]$ represents a significant crossover between two distinct behaviors. In addition to the strong–fragile glass-forming transition, a clear dependence on pressure is evident within the experimental errors before and after this temperature: (i) for higher temperatures the obtained rotational relaxation time values are P-independent, while (ii) below $[eqn]$ the dependence on P is increasingly significant, increasing both with decreasing T and increasing P.

Finally, we considered that $[eqn]$ marks a transition from a T region, where the water dynamics is governed by only two energy levels with a single relaxation process (the Arrhenius configuration), to another typical glass-forming liquid in which temperature decrease originates large intermolecular correlations. In the time and length scale, this lead to dynamic clustering (the HB tetrahedral network) with a large distribution of possible energies and thus multirelaxations (a super-Arrhenius system).

Hence, like in complex liquids, this complex energetic situation determines a behavior for all the transport parameters of water (relaxation times, viscosity, and self-diffusion) characterized by this specific dynamic crossover, a situation unanimously proposed by the liquid state theories of crossover temperature that marks the transition from normal to supercooled liquid [62,63,64].

Physical reality is observable in the dynamics of such systems, that is, both in their time and relaxation processes; in the case of water, these thermodynamic realities are correctly associated to HB interactions. These interactions give rise to specific transient structures that become increasingly stable as the material cools within the metastability region. In the case of liquid water, all of this originates at $[eqn]$ via the HB clustering process. This temperature is therefore the magic point where the water system becomes the well-known intriguing and complex material. In our opinion, only by considering the dynamical quantities reported in this work in a unified manner can the importance of this temperature be fully appreciated. Hence, as said, all the reported data measured in the large P-T phase diagram considered and the consistency proposed by the response functions $[eqn]$ and $[eqn]$ indicate that $[eqn]$ plays a primary role in the water physics and it is the locus where its anomalies originate.

For these reasons, we conclude our analysis by directly looking at the liquid water local HB structure effect on the oxygen–oxygen distance ( $[eqn]$ ), evaluated from the density data reported in Figure 1, for temperatures below 373 K up to the deep supercooling regions.

2.6. The Intermolecular Oxygen–Oxygen Distance

The water molecules possess high asymmetry with coordination numbers ranging from two to four or even greater [65]. The water’s liquid structure, as said, can be interpreted in terms of either tetrahedrally coordinated structures with thermal fluctuation or the mixed-phase of low- and high-density fragmentation. Its symmetry is the C_v2_ group; therefore, an oxygen atom tends to find four neighbors to form a stable tetrahedron, but the repulsion between electron pairs on oxygen prevents, in the liquid phase, its stability.

The corresponding packing order follows, in its entire phase diagram except for highest T and P, the Pauling’s rule [66]. The size of liquid water molecules is essentially stable despite fluctuations in both the O:H non-covalent bond lengths and bond angles. The average intermolecular distance changes only at the transition to the solid phase, when the system assumes an ordered structure and is therefore lower in density than the liquid.

Despite these thermal fluctuations, the Ice Rule extension results for both bulk liquid and ice in an ideal tetrahedron made by two equivalent molecules and between them four identical O:H-O bonds of different orientations. It is possible (see, e.g., ref. [67]) to have a direct relation between the water density and the distance $[eqn]$ by considering the molecular mass M with $[eqn]$ 1 (gcm^−3^) at 4 °C (under the atmospheric pressure). This defines unambiguously the density-dependent molecular separation, $[eqn]$ , and the next nearest-neighboring distance $[eqn]$ (in Å). Thus, we have $[eqn]$ and $[eqn]$ . By using these simple relations and the experimental density data, the intermolecular distance between two first neighbors water molecules can be evaluated.

For the different isobars considered (from $[eqn]$ to 200 MPa, together with the experimental data of the TIP4P MD [29,30] also reported), Figure 7 reports the difference of the molecular separations obtained at the different temperatures, $[eqn]$ , with respect to the corresponding values measured at $[eqn]$ . Regarding the MD data reported here and calculated via the TIP4P model using a radial distribution function, it must be specified that they represent the average distance between oxygen atoms in water molecules. What emerges from the comparison between these data and those obtained here via the density values is a coincidence in their values, at all isobars considered, for all temperatures of the stable phase, while they show gradually increasing differences within the supercooled metastable phase. This is due to the known sensitivity of the model used in regions characterized by metastability; in such cases, it is preferable to use other models, such as E3B3 [31]. However, beyond this situation, both datasets show a similar thermal evolution. The data at 135 K correspond to the values of the low-density amorphous phase (LDA), values towards which the different isobars of the liquid phase evolve. In addition, as can be observed, these quantities increase as the temperature decreases with strong differences in their rates as the pressure increases. We can summarize this by saying that it is all due to the thermodynamic effects on the system polymorphism, created by the HB in particular, the increasing ratio under supercooling between the two phases LDL and HDL, and also the pressure effects on this.

3. Materials and Methods

Our study focused on the water system, in its liquid phase, in order to highlight the specificity of its particular thermodynamic singularity ( $[eqn]$ ). It was arranged by considering essentially the isobars of a large lot of literature data: the density, the constant pressure specific heat and the transport functions such as self diffusion and rotational relaxation times. All this over a wide range of pressures and temperatures including both the stable and metastable (supercooled) phases. As mentioned, the data considered essentially concern bulk water, although the experimental accessibility to the metastability region was amplified using confined or emulsified water.

As underlined, the nuclear magnetic resonance, NMR, technique was used in order to obtain data on the transport functions ( $[eqn]$ and $[eqn]$ ) as well as those on the specific heat. In this last case $[eqn]$ was also obtained through the Adam–Gibbs relation and compared with those measured through a differential calorimeter. The isobarsi of the Isothermal Compressibility and the Isobaric Thermal Expansivity were then obtained as the thermodynamic derivatives (of the isotherms and isobars) of the density. Finally, the Intermolecular Oxygen–Oxygen Distances, $[eqn]$ , where obtained through the extension of the Ice Rule to the liquid phase.

4. Conclusions

We considered the experimental water density isobars in a large area of the phase diagram including data from an MD simulation [29,30] and the measured data in the two characteristic glasses, LDA and HDA [5,6,7,8,9,42]. The temperature range considered was from 100 to 380 K, while the pressures considered ranged from atmospheric ( $[eqn]$ MPa) to 800 MPa. Special attention was paid to the liquid phase in the supercooled region below the melting temperature $[eqn]$ and to a particular temperature $[eqn]$ K where the compressibility ( $[eqn]$ ) and expansivity ( $[eqn]$ ) show unique thermodynamic properties [34]. For the supercooled region, we carefully considered the polymorphism of the system, characterized as is known by the presence of two liquid phases: LDL (due to HB) and HDL.

As in complex liquids, in water, the interaction process originates disordered correlated regions influencing the transport functions by means of a dynamic crossover. In water, this situation is entirely due to the specificities of the HB bond, which means tetrahedral networking for the system structure. These considerations indicate that a change in the water dynamics occurs at $[eqn]$ , the temperature at which networking takes place. Change observable through a transition from a high-T region is characterized by an Arrhenius behavior; this change is to another typical glass-forming liquid system in which the temperature decrease gives rise to increasing intermolecular interactions and correlations with growing activation barriers resulting in super-Arrhenius behavior characterized by multirelaxations. Thus, this clustering process with increasingly extended structures (LDL) in the subcooled metastable phase is, unlike the disordered phase typical of high temperatures (HDL), strongly dependent on pressure.

Taking all this into account, we focused on this crossover temperature to verify how it determines the physical chemistry of water. We then carefully examined, over a large region of its phase diagram, the isobars of various thermodynamic functions such as the specific heat $[eqn]$ , the self-diffusion $[eqn]$ , and the rotational relaxation times $[eqn]$ , finding full suggestion on the universality of this specific temperature. Obviously, this is a suggestion and not a certainty, which could be supported by new studies, especially on the dynamic properties of the system in its metastable phases, both theoretical and experimental, also using the suggestions of molecular dynamics studies.

Finally, we paid attention to the distance $[eqn]$ (the intermolecular distance between two first neighbors water molecules) evaluated from the density isobars, demonstrating that these intriguing realities are entirely due to the polymorphism of the liquid in its supercooled metastable phase.