The Water Polymorphism and the Liquid–Liquid Transition from Transport Data

NMR spectroscopic literature data are used, in a wide temperature-pressure range (180–350 K and 0.1–400 MPa), to study the water polymorphism and the validity of the liquid–liquid transition (LLT) hypothesis. We have considered the self-diffusion coefficient DS and the reorientational correlation time τθ (obtained from spin-lattice T1 relaxation times), measured, respectively, in bulk and emulsion liquid water from the stable to well inside the metastable supercooled region. As an effect of the hydrogen bond (HB) networking, the isobars of both these transport functions evolve with T by changing by several orders of magnitude, whereas their pressure dependence become more and more pronounced at lower temperatures. Both these transport functions were then studied according to the Adam–Gibbs model, typical of glass forming liquids, obtaining the water configurational entropy and the corresponding specific heat contribution. The comparison of the evaluated CP,con f isobars with the experimentally measured water specific heat reveals the full consistency of this analysis. In particular, the observed CP,con f maxima and its diverging behaviors clearly reveals the presence of the LLT and with a reasonable approximation the liquid–liquid critical point (LLCP) locus in the phase diagram.


Introduction
Water, starting from the perspective of biology, has a basic role in many research fields and technological applications, regardless of whether it is in bulk or confined [1]. In chemical physics, it is of importance due to its unusual thermodynamics, compared to normal liquids, and as a prototype of supercooled liquids [2]. This is reflected in its well-known anomalies for almost all of its properties as a function of thermodynamic variables, especially below its melting temperature T m down to the homogeneous nucleation temperature (T h ). Examples are represented by the pressure (P) and temperature (T) behaviors of its density (ρ) and the thermodynamic response functions (isobaric specific heat (C P ), the compressibility (isothermal κ T and adiabatic κ S ), and the expansivity (α P )) all related to the volume (δV) or entropy (δS) local fluctuations. For common and regular liquids, these fluctuations are positively correlated and decrease as T decreases; for water, below T m , they not only grow but become anticorrelated so that an V increase brings an entropy decrease due to a growing local order accompanied by the observed diverging (critical like) behaviors in the mentioned response functions [3].
Another relevant characteristic of water which reasonably seems to be the basis of its anomalies is the polymorphism that characterizes it in all its phases including the liquid one polymorphism [4][5][6][7]. Although the polymorphism of the solid water crystalline phase has been known for a long time, i.e., the ice has many different structural forms ranging from the ice Ic to ice XII [8], that of amorphous water is a relatively recent discovery [9][10][11]. So, from the certainty of this "polyamorphism", the idea of a liquid polymorphism [12][13][14][15][16] was proposed.
Specifically, the amorphous water phases have different densities: the water highdensity amorphous phase (HDA) [9,11] and the low-density amorphous phase (LDA) have been known since 1935 [17], and finally the VHDA (very high-density amorphous phase) [18]. Of relevant interest is that these two amorphous phases can be transformed into each other, respectively, through a reversible first order transition [10,19]. Furthermore, at ambient pressure, the LDA, if heated, undergoes a glass to liquid transition (at about 130 K) into a highly viscous fluid and then crystallizes at T x = 150 K. On this basis, the liquidliquid transition hypothesis (LLT) based on the liquid polymorphism was developed [4]. A model, related to an MD study, is nowadays central in water studies, being at the base of the liquid-liquid critical hypothesis (LLCP or second critical point in distinction to the vapor-liquid one). Like in the glass, liquid water has two liquid forms of different densities (the high-and low-density liquids, respectively, HDL and LDL); they coexist and, depending on P and T, can change one into the other by means of a first order transition: the liquid-liquid transition hypothesis (LLT). Whereas the HDL has a disordered structure (made of monomers, dimers and trimers), the LDL is characterized by an "open" structure governed by a networking process with a tetrahedral symmetry due to the noncovalent attractive hydrogen bonding (HB) interaction [4].
For precision, together with the HB, the water molecules also interact by the Coulomb repulsion between electron lone pairs on adjacent oxygen atoms and two H-O covalent bonds originating from the sharing of the electron lone pairs. The first one dominates in the stable and supercooled regimes and the repulsive potentials mainly influence the water physics from above the boiling temperature (T b ) in the sub-critical and critical regions.
In particular, the LDL tetrahedral symmetry is that of ordinary ice, with four nearest neighbors around the water molecule (also acting as a H-donor to two of them and a H-acceptor for the other two). In ice, this network is permanent while the liquid water tetrahedrality is, instead, local and transient. It should be noted that a pressure increase contrasts these ordering effects, whereas a T decrease both involves its growth in size and stability; experiments show that the HB lifetime strongly increases (many orders of magnitude) from picoseconds values characteristic of the stable liquid water [20].
The tetrahedral LDL local networking, as well as the liquid polymorphism, originates the entropy decrease and the diverging behavior observed in water functions, explaining the observed water anomalies and complexity. In addition, the cited studies have allowed a more precise definition of the water phase diagram reported in Figure 1. The LLCP is estimated to be located near 200 K and at a pressure of ∼200 MPa, very far from the locus of the vapor-liquid critical point CP (T C = 647.1 K, P C = 22.064 MPa). However, the bulk liquid water, in principle, cannot exist stably in the region between the homogeneous nucleation temperature (T h ) and that in which the ultra-viscous liquid obtained from the fusion of LDA crystallizes (T x ), and the LLCP seems to be located just inside this region, called "no man's land" [2]. Although such a region (T h − T x ) is open to MD simulation studies, some expedients have been used to gain experimental insight. Examples are confined water in nanopores (smaller than the nucleation centers) [21], around or inside macromolecules, in solutions, in ice, in emulsion [22][23][24] and micellar systems [25] or by melting a multimolecular thickness of an ice surface [26]. In particular, in this last case and for water in nanopores water can be easily maintained in the liquid state also in all the range T h − T x and the LDA can be also achieved [21,26,27]. The P − T water phase diagram. All the lines characterizing the chemical physics of the system in the liquid and disordered phases (glasses) are illustrated together with the homogeneous nucleation temperature (T h ) and that of crystallization of the ultra-viscous liquid (T x ) defining the "no man's land". The Widom line (characteristic of liquid-liquid transition) is also proposed together with those of the maximum density (ρ max ) and the melting temperature (T m ), the critical point C P and the estimated (by MD) LLCP (C ) [2]. T * is the temperature of the minimum of the isothermal compressibility κ T [28]. Finally , the regions of the polymorphism (LDL + HDL) and those for T > T* composed only of HDL are shown.
In this way, many important water properties due to the polymorphism and related to the LLT (and LLCP) were discovered. Very relevant in the present context are the results obtained from MD studies [29], pointing to the existence of the so called Widom line (WL) which identifies the P − T locus of the maximum in the δV and δS fluctuations where thermodynamic response functions reach their extremes (minimum with negative values in the α P and maxima in C P and κ T ). Other experimental results concerning the water dynamics and structure related with the Widom line and the LLT are at ambient pressure: (a) the dynamic crossover from a fragile to a strong glass-forming material, originally predicted by Angell [30] and observable at T L 225 K [31] that is also the locus of the Stokes-Einstein relation violation (due to the onset of the dynamic heterogeneities and the decoupling between the translational and rotational modes); (b) the compressibility maximum [32,33]. The P − T locus of the WL was investigated by using neutron scattering [34], and the recent studies (experimental and simulation ones) on the isothermal compressibility agree with the related findings [35][36][37][38][39]. Another phenomenon observed by using confined water is the existence of a density minimum [40,41], as predicted more than a century ago by Percy W. Bridgman [42] and subsequently confirmed by computational studies [36][37][38].
Many of the main suggestions regarding the water thermodynamics come from the large number of accurate computational studies [43], but the LLCP (inside the supercooled regime) is far off being experimentally proven in a definitive way. Due to its localization inside the "no man's land", it proved to be elusive in all the experimental attempts, although the water polymorphism and the LLT have been tried [44]. Nevertheless, the liquid polymorphism, which is favored by the temperature decrease and the corresponding growth of the hydrophilic interaction represented by the hydrogen bond, has been widely proven [5].
Very recently, according to the Adam-Gibbs model (developed to clarify the cooperative relaxation processes in glass-forming liquids [45]) and by using dynamical data (self-diffusion at ambient pressure and in a very large temperature range [21,26]), it has been shown that the specific heat in the supercooled water is largely due to configurational effects, and therefore to the LLT [46,47]. In such a way, the calculated configurational specific heat C P,con f was compared with that measured experimentally, thus confirming the fact that the temperature of the corresponding maximum is coincident (see Figure 2). Such a maximum, as theoretically predicted in terms of polymorphism [48], also defines the WL. The present work deals with an extension of such an analysis to a wider range of the P − T phase diagram including transport data for pressures up to 400 MPa and inside the supercooled region. The aim, by essentially using bulk water literature data [20,22,49], is, on the one hand, to gain more information on the water thermodynamics and, on the other, to obtain more precise signals than those currently available from the LLCP and verify its localization. It must be stressed that these frame transport function data (such as self-diffusion and relaxation) are very relevant because their values provided experimental evidence of LLT, showing, by decreasing T, both the Widom line and the violation of the Stokes-Einstein relation, as well as the HB networking and LDL's dominance over HDL [44]. The figure reports the liquid water specific heat C P (T, P), measured in the temperature range 100-430 K for 0.1 MPa, and that of the Ih ice (red lines and dots [50]); for T > 250 K data at P = 50 , 100, 150, 200, 300 and 400 MPa are also reported. For ambient pressure, data come from different experiments in bulk [51,52] and confined water (in nanotubes of 2.2 nm [47]). Data for P > 50 MPa deal instead only with bulk water [53]. The inset shows the difference ∆C P = C P,liq − C P,sol (squares) and the configurational C P,con f evaluated by using the Adam-Gibbs model [47].

Results and Discussions
Our starting point is represented by Figure 1, which illustrates the liquid water specific heat C P (T, P), measured in the temperature range 100-430 K at 0.1 MPa, and that of the Ih ice (red lines and dots [50]); for T > 250, K data at P = 50, 100, 150, 200, 300 and 400 MPa are also reported. For ambient pressure data coming from different experiments in bulk [51,52] and confined water (in nanotubes of 2.2 nm [47]). Data for P > 50 MPa deal instead only with bulk water [53]. We can assume [46] that, for water, the difference between the liquid and solid specific heat represents a good estimation of the configurational contribution, so that C P,con f ∆C P = C P,liq − C P,sol . Therefore, in the figure inset, ∆C P is reported as the difference between the specific heat values measured in the liquid and in the ice Ih, respectively, and the C P,con f evaluated by using the Adam-Gibbs model (AG) theory as described in the following [47]. For the confined water, the metastable supercooled region is explored up to about 90 K, maintaining it in the liquid state; in fact, the large peaks of latent heat that accompany solidification are absent [47,54,55].
The AG model was proposed to explain the relaxation temperature dependence in glass-forming liquids. It was detailed, in molecular-kinetic terms, by accounting the thermal effects in the size of the cooperatively rearranging regions of different energetic configurations; sizes determined by configuration restrictions, and thus expressed in terms of their configurational entropies. According to the theory, these cooperative regions have a transition probability W(T) = Fexp(−z∆µ/k B T) that can be accounted for in terms of its size z and ∆µ (the potential energy hindering cooperative rearrangements). F is a frequency factor (negligibly T-dependent) and k B is the Boltzmann constant. By expressing the cooperative region "critical size" z * as a function of the molar configurational entropy S con f , the transition probability can be expressed as W(T) = Aexp(−C/TS con f ). As the system relaxation times are related to the transition probability as τ(T) ∝ W(T) −1 , the system self-diffusion can be written as: D S0 and A = z∆µ can be assumed as constant (at a given concentration). The D S0 value can be estimated from the D S (T) in the high T limit. In such a way, the configurational entropy can be obtained from the system measured diffusion data (or the transport functions) and the configurational C P,con f was evaluated as C P,con f = T(∂S con f /∂T) P . By using such an approach, the water configurational C P,con f was evaluated at ambient pressure from the bulk water diffusion data (measured and simulated in the range 373-237 K) [46,47], obtaining D S0 = 1.07 × 10 −7 m 2 s −1 and A = 31.75 kJmol −1 . The corresponding analysis was also performed by considering confined water data [21,26,27] (373-120 K) and determining the same values of D S0 and A. The results, in terms of C P,con f , by using this procedure are reported in Figure 2 [47].
The diffusion data used in the present analysis were obtained in 1988 from the H.-D. Lüdemann team [20] with an NMR experiment by using the pulsed field gradient spin echo technique [56] in bulk water at different pressures up to 400 MPa and temperatures down to 200 K. The corresponding data are illustrated in an Arrhenius plot at the different isobars in Figure 3a, together with data measured at ambient pressure in bulk and confined water used to evaluate the C P,con f proposed in Figure 2. As can be seen, the data show a marked difference with respect to those at ambient pressure only in the supercooled regime, where the corresponding values also increase with pressure, while for higher temperatures of the liquid stable phase the corresponding variations are significantly less pronounced. However, they show a pronounced non-Arrhenius temperature dependence and an apparently diverging correlation length of the supercooled water fluctuations. This behavior was accounted for by using the dynamic scaling behavior typical of modecoupling theory, D S (T) = D 0 (T − T L ) γ [3,57].
Some years before these observations, the same team proposed another NMR study on water relaxation in a meso-sized emulsion [22], hence decreasing by more than a decade the studied temperature with respect to that of the bulk water and detailing the analysis in the pressure range 1-250 MPa in steps of 25 MPa. The longitudinal proton relaxation times T 1 of the water protons have been determined, by means of the inversion recovery pulse sequence, at 100.1 MHz, at temperatures up to 186 K, for pressures higher than 200 MPa. The corresponding proton relaxation rate (R 1 = 1/T 1 ) reported in Figure 4 at the different isobars shows a non-Arrhenius increasing dependence. It may also be noted that above the melting temperature, all the measured data, such as the D S (T) values in the same T−range, are, within the error, pressure independent. In addition, the temperature T 315 K identifies, by decreasing T, a crossover from an Arrhenius to a non-Arrhenius behavior. From the observed behavior of the thermodynamic functions, this temperature is remarkable for water, being the locus of the onset of the HB tetrahedral structure [58,59]: it is in fact the place of the minimum, at all pressures, of isothermal compressibility (κ T (P, T)) and also represents the point where all the lines of the expansivity (α P (P, T)) cross each other [28].  From R 1 , according to the current theoretical models [60][61][62], a reorientational correlation time (τ θ ) can be evaluated by assuming that the proton spin system R 1 is in general mediated via magnetic dipole couplings between nuclear magnetic moments as intramolecular or intermolecular: (R 1 ) Meas = (R 1 ) Intra + (R 1 ) Inter . The first one monitors reorientational motions only, whereas the intermolecular one is sensitive to both positional and orientational rearrangements. In particular, the (R 1 ) Intra is a function of both τ θ and the Larmor frequency ω L as: (R 1 ) Intra = (3γ 4 2 /10r 6 )((τ θ /(1 + ω 2 L τ 2 θ )) + (4τ θ /(1 + 4ω 2 L τ 2 θ ))), so that two motion regimes, a fast (ω L τ θ < 1) and a slow (ω L τ θ > 1), are possible. γ is the proton gyromagnetic ratio, is the Planck's and r is the distance between the water molecule hydrogens.
However, the R 1 isobars at about 200 MPa and above run through a maximum at about 197 K (when ω L τ θ 1) because an increase in pressure causes an increase in the correlation time, a situation originally predicted for normal liquids [61]. As the compression increases, the molecular motion is reflected in the correlation time. Furthermore, according to the asymptotic behavior (for ω L τ θ 1, it is R 1 ∼ τ θ , whereas for ω L τ θ 1 we have R 1 ∼ 1/τ θ ), and τ θ was also evaluated for 200 MPa in the very slow regime [22]. According to this, Hindmann et al.'s suggestion [63,64] that the water relaxation time T-dependence, at constant pressure, is due to two contributions (exponential) was assumed to be correct. One, at low T, is related with the cooperativity effects of the HB clustering and the other (at high T) to the breaking of a single HB; additionally, in these studies a relationship between viscosity, diffusion and rotational relaxation has been discussed [64]. Figure 5 shows, in the interval 200-303 K, all the τ θ isobars corresponding to the R 1 , of Figure 4, evaluated according to the above expression for (R 1 ) Intra . For each isobar, the lowest experimentally accessible temperature increases progressively with increasing pressure, from 237 K at 5 MPa to 205 K at 175 MPa. Only the 200 MPa isobar, evaluated according to the asymptotic behavior, is reported up to 185 K.
When the relaxation rates become dependent on the Larmor frequency (ω L ), their isobars are characterized by a dispersion and its frequency dependence revealed additional properties on the molecular motions. A study conducted at different Larmor frequencies has proved the water reorientation isotropy on a nanosecond time scale and that a single reorientational mode with a strong non-Arrhenius temperature dependence appears to be adequate to reproduce the experimental R 1 (T, P, ω L ) shape [49]. The situation was related to the increasing HB lifetime at low temperatures and to random, transient HB network fast quasi-lattice vibrations. Figure 5. The reorientational correlation time τ θ isobars, corresponding to the R 1 illustrated in Figure 3, evaluated according to the current NMR theoretical models [60][61][62]. The dotted lines correspond to a data fitting in terms of the mode coupling theory scaling law for critical phenomena.
According to the same analysis made on the self-diffusion data of bulk and confined supercooled water at ambient pressure, originally proposed to observe the dynamic crossover [30,31], we have fitted the NMR reorientational correlation time isobars according to the MCT reported as continuous lines in Figure 5. The obtained exponent values are for all the isobars γ 2, suggesting, for this reorientational relaxation, universal behavior, based on the energy landscape concept [65] typical of the dynamics of glass-forming liquids [66,67]. This situation seems to confirm some relationship in the deep supercooled regime between the rotational relaxation and transport parameters such as D S . Such an analysis, even taking into account that the accessible experimental data do not reach the lowest measured in confined water, provides a reliable estimation of the crossover temperature T L . A possible crossover at ∼ 200 K appears to be observable at 200 MPa.
Having this result, we tried to compare these values of τ θ with those of the diffusion shown in Figure 3a. For this, we have considered a linear relationship similar to that used in the scattering spectroscopy between longitudinal relaxation (τ) and the diffusion for low wave-vector Q: 1/τ = D S Q 2 . For the data normalization, we used the values corresponding to high temperatures which in both cases (τ θ and D S ) appear to be pressure independent, obtaining the used factor of 1.587 × 10 20 . Figure 3b displays such a normalization showing similar behaviors between both data, thus suggesting a sort of coupling, to the same isobars, between these two transport parameters, at least for common temperatures.
The next step is the evaluation of configurational entropy and specific heats from the self-diffusion data (Figure 3) by using the AG formalism. By also considering the behavior of the data at a high temperature, we assumed for all the isobars the same value of D S0 = 1.07 10 −7 m 2 s −1 and A = 31.75 kJmol −1 used for ambient pressure. Regarding the values corresponding to the rotational relaxation time, we will consider only the isobars of 175 and 200 MPa, just for a form of comparison with the effective self-diffusion data. The obtained S con f isobars are illustrated in Figure 6 in a linear scale for the temperature range 100-300 K. For all reported isobars, the behavior of the data obviously reflects that observed in diffusion: continuously decreasing with decreasing T and increasing with pressure (in the supercooled region) at least up to 200 MPa, a continuous and decreasing behavior can be observed with decreasing T, after which they seem to stabilize. The ambient pressure data, covering the entire temperature range from the stable region, well above the melting point (300 K), dominated by the HDL to that of the pure LDA region (100 K), illustrate the way in which water polymorphisms evolve with temperature. The LDA is the vitrified LDL, and LDA configurational entropy can be assumed to be essentially the continuous and slow evolution of the very low temperature (T < 160 K) liquid component. Taking this into consideration, such a S con f isobar reflects, as correctly shown by N. J. Hestand and J. L. Skinner (HS) by using a logistic function for growth modeling to describe the corresponding water diffusion and radial distribution function [68], the relative amount of the two liquid components, LDL and HDL, according to a logistic function. According to this, the data flex point, where the LDL and HDL are in equivalent amounts, is to this isobar the temperature of the WL. The corresponding specific heat, reported in Figure 7, was evaluated by performing the isobars' derivative C P,con f = T(∂S con f /∂T) P after a fit of the entropy data by means of the same HS procedure and the same logistic function [68]. In the case of the S con f isobars, coming from the bulk water self-diffusion (50,100,150,200,250, 300 and 400 MPa), the data fitting was carried out with a third order polynomial, only for the rotational relaxation time τ θ , evaluated from the emulsioned water relaxation rates R 1 , 175 and 200 MPa, we used a polynomial of the fourth order. Figure 6. The S con f isobars calculated according to the AG theory are illustrated, at different pressures, in a linear scale for the temperature range 100-300 K. Figure 7. The C P,con f isobars calculated according to the AG procedure by using the bulk water transport data, in the P range 0.1-400 MPa, are shown. For comparison, the bulk water experimentally measured ∆C P at 50, 100, 150, 200, 300 and 400 MPa, are also reported [53].
All the C P,con f isobars calculated with this procedure are shown in Figure 7, where are also reported, for comparison, the bulk water's experimentally measured ∆C P at 50, 100, 150, 200, 300 and 400 MPa [53]. A reasonable agreement between the temperature behavior of these latter experimental data and those evaluated in accordance with the procedure used can be observed, and the symbol size represents the experimental error. The experimental values at 323 and 298 K were considered as reference values after the entropy derivative. The entire behavior of these data is representative of the thermodynamics of water liquids in this phase of metastability. In spite of the limited temperature region inside the supercooled region (these calorimetric data coming from bulk water), the thermal evolution of the corresponding different isobars fully supports the presence of a LLCP, also giving information on where it is located in the phase diagram. In fact, if we look only at the C P,con f data derived from self-diffusion a different growth rate in isobar values can be easily observed between the isobars between 0.1 and 150 MPa and those for P > 200 MPa. The former data, at 50, 100, 150 MPa, are characterized by a divergent behavior that involves all data at all temperatures explored, even if they go more and more into the metastability regime (also if the lowest temperature value of the isobar at 150 MPa is slightly lower than the C P,con f maximum of that at 0.1 MPa). In the opposite case, P > 200 MPa, not only is the configurational specific heat growth rate less pronounced but also accompanied by a moderate change in its slope and an inflection point. This latter point therefore suggests that the LLCP is located in a pressure between 150 and 200 MPa.
The overall situation becomes clearer if we consider the specific heat data, at the 175 and 200 MPa isobars, coming from the rotational relaxation time τ θ . In this case, the increase in the temperature interval within the metastable region highlights both a maximum, and thus a WL point in C P,con f at 200 MPa (∼195 K) and a very pronounced divergent behavior at 175 MPa. In addition, the maximum value in C P,con f at 200 MPa ( 44.3 JK −1 mol −1 ) is smaller than that measured at 0.1 MPa ( 66.1 JK −1 mol −1 ), whereas the highest value obtained for the 175 MPa isobar is 72.7 JK −1 mol −1 at 206 K. These further data, in agreement with the precedents due to self-diffusion, not only give confirmation of the presence of a LLT but also indicate the presence of LLCP located in the following region of the phase diagram: an isobar P C between 175 and 200 MPa and a critical temperature 195 < T C < 206 K. Such a result obtained from dynamical data inside the no man's land, also by considering the full agreement with the finding coming from the isothermal compressibility (evaluated from the density data) studies [35,44,69], undoubtedly provide crucial information regarding our understanding of the strange behavior of liquid water. We stress that this resulting and complex evolution of the configurational modes of liquid water clearly shown both by the S con f and C P,con f isobars due to its polymorphism is certainly the result of its underlying energy configuration: the energy landscape or inherent basins of energy [65].

Conclusions
Starting from the idea that liquid water is dominated by a polymorphism generated by the HB interactions and that this polymorphism, made from the LDL and HDL, dominates the behavior of thermodynamical functions, we have analyzed a lot of experimental data that enter well inside the no man's land in order to determine what they can tell us about the presence of the LLT and more generally what they reveal about the LLCP hypothesis. By considering literature dynamical data obtained by means of the NMR technique, such as the self-diffusion coefficient (D S ) and the spin-lattice (T 1 ) relaxation times, measured, respectively, in bulk and emulsion water at several different isobars in the range 0.1 to 400 MPa, we found that their temperature dependence provides further evidence for the LLCP hypothesis. The spin-lattice times were analyzed, as the proton relaxation rate (R 1 = 1/T 1 ), according to the NMR current theoretical models [60][61][62] evaluating a reorientational correlation time (τ θ ) of the order, in the stable liquid phase, of some picoseconds. Both these data evolve with T by changes of several orders of magnitude, whereas their pressure dependence, almost nothing in the stable liquid phase (T > T m ) and as an effect of the LDL HB networking of the LDL liquid, becomes more and more pronounced at lower temperatures.
Both these quantities, diffusion isobars first and then some of τ θ , have been studied according to the Adam-Gibbs model, typical of glass-forming liquids. A situation was recently studied using the ambient pressure isobar, also considering data of confined water and the molten amorphous phase, and covering a very wide temperature range from the boiling temperature to that of LDA amorphous phase. The aim was to highlight the behavior of the water configurational entropy and the corresponding specific heat contribution in order to clarify the LLT behavior on approaching LLCP.
A comparison of the thermal behaviors of all the evaluated C P,con f isobars, in particular, the observed maxima and the evolutions of the diverging behaviors, clearly reveals the presence of the LLT and, with a reasonable approximation, the LLCP locus in the phase diagram. This latter situation is fully consistent with recently studies on the water thermodynamic functions obtained from experimental data [35,44,69] and proper MD simulations [36][37][38][39]68]. Finally, we underline that the observed significant configurational evolutions due the liquid water polymorphism are related to the local potential minima, known as inherent structures (ISs), surrounded by potential energy basins [65].
Author Contributions: All authors contributed equally. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.