_{1}

^{*}

A revised phase diagram for water shows three distinct fluid phases. There is no continuity of liquid and gas, and no “critical point” on Gibbs’ density surface as hypothesized by van der Waals. A supercritical colloidal mesophase bounded by percolation transition loci separates supercritical liquid water and gas-phase steam. The water phase is bounded by a percolation transition (PA) of available volume, whereas steam is bounded by the loci of a percolation transition (PB) at a density whereupon a bonded molecular cluster suddenly percolates large distances. At the respective percolation densities, there is no barrier to nucleation of water to steam (PA) or steam to water (PB). Below the critical temperature, the percolation loci become the metastable spinodals in the two-phase coexistence region. A critical divide is defined by the interception of PA and PB the p-T plane. Critical parameters are obtainable from slopes and intercepts of pressure-density supercritical isotherms within the mesophase. The supercritical mesophase is a fourth equilibrium state besides ice, water and steam. A thermodynamic state function rigidity (dp/dρ)_{T }defines a distinction between liquid and gas, and shows a remarkable symmetry due to an equivalence in number density fluctuations, arising from available volume and molecular clusters, in liquid and gas respectively. Following an earlier debate in these pages [“Fluid phases of argon: A debate on the absence of van der Waals’ critical point” Natural Science 5 (2) 194-206 (2013)], we here report further debate on a science of criticality applied to water and steam (APPENDIX 1).

The thermodynamic diagram of water phases is described in a recent educational review [

The conventional phase diagram of water, and indeed all other pure fluids, is presently based upon the van der Waals hypothesis [_{c},p_{c} where upon an equation-of-state for which the first two derivatives of pressure, i.e. with respect to changes in volume at constant T, go to zero. This has been the accepted description of liquid-gas criticality since van der Waals published his renowned thesis “On the continuity of liquid and gas”. The hypothesis further implies that for all supercritical temperatures the thermodynamic state functions e.g. density ρ(T,p) or Gibbs energy G(T,p), from the gas at low density, steam, to the liquid at high density, water, are continuous in all their derivatives for supercritical temperatures across the whole density range.

Another hypothesis that has been fashionable amongst theoretical physicists for about 50 years is the concept of universality [

“A critical point is a special example of a phase transition. Consider, e.g. the water-steam transition. Suppose that water and steam are placed under pressure always at the boiling temperature. At the critical point, the distinction between water and steam disappears, and the whole boiling phenomena vanishes. The principal distinction between water and steam is that they have different densities. As the pressure and temperature approach their critical values the difference in density between water and steam goes to zero”.

It is evident from modern 5-figure precision density data (_{c}, the density difference does not go to zero. Rather, we have found that a liquid-gas critical “point”, as hypothesised by van der Waals and accepted by the physics community for 140 years, does not exist a such. Moreover, there is no universality of critical phenomena that spans the dimensions. In the ρ(T,p) surface there is a dividing line of uniform

chemical potential connecting maximum and minimum coexisting gas and liquid densities, respectively, above which there exists a supercritical mesophase bounded by weak higher-order percolation transitions, and below which there is the familiar liquid-vapour coexistence region. Furthermore, and contrary to Wilson’s statement above [_{c} water and steam remain distinctly different phases, and are separated by a supercritical mesophase.

This new description of liquid-gas criticality has also been shown to describe all the experimental thermodynamic phenomenology of liquid and gaseous argon [

Percolation transitions occur when structural properties of a molecular system, for example clusters of molecules in the gas phase, or holes in the liquid phase, on reaching a certain density suddenly coalesce from being of molecular dimension to become macroscopic, i.e. system spanning. This gives rise to higher-order thermodynamic phase transitions, in which there are discontinuities in second and/or higher derivatives of Gibbs chemical potential (μ) with respect to temperature or pressure, notably: isothermal compressibility (d_{2}μ/dp^{2})_{T}, heat capacity (d_{2}μ/dT^{2})_{p} and thermal expansivity (d_{2}μ/dpdT) all of which undergo some degree of change at, or in the vicinity of, percolation transitions.

We begin therefore by taking a look at the behaviour of these three second-order thermodynamic state functions in the case of supercritical water (

The first observation from _{c}. The observation of these lines has a long history that goes back to J. D. Bernal [_{p}), which can be clearly observed in _{p}) (

The second observation that we make from _{p}, α_{p} and κ_{T}, all increase with temperature along an isobar, whereas in the gaseous state they all decrease with temperature. Thus there is a fundamental thermodynamic defining distinction between liquid and gas that extends to supercritical temperatures. Below, we find that defining property is the work required to change the density, i.e.an inverse compressibility, or rigidity (dp/dρ)_{T} where ρ is the density. The rigidity of a gas decreases with density and increases for a liquid, i.e. (d_{2}p/dρ^{2})_{T} is negative for a gas and positive for a liquid; its value is zero in the mesophase. Lines of higher-order discontinuity in the rigidity can be characterized, not just for water, but all liquids, in terms of percolation loci that bound the supercritical mesophase separating gas from liquid.

Percolation transitions in model fluids of hard-spheres and square-well molecules are well defined by a characteristic distance, which unambiguously divides configuration space volume into “sites” and “holes”. The densities at which small clusters of occupied sites in the gas phase, or unoccupied holes in the liquid phase, first span the whole system are referred to as percolation transitions, and denoted by PB and PA, referring to “bonded cluster” (sites) and “available volume” (holes) respectively [

water, however, the steep hydrogen bond is fairly sharply defined for water clusters. At a given temperature, the total volume of any equilibrium configuration of water or steam can be divided in occupied sites, not accessible to an incoming molecule, and holes which would be accessible to an additional molecule. At these percolation transitions, thermodynamic state functions can change form due to sudden changes in the state-dependence of density and/or energy fluctuations. It is the change in (dp/dρ)_{T} that defines PA and PB phenomenologically.

In order to understand “mysterious lines” [

The critical divide from _{c} is shown in _{T} data [_{c}, are obtained from the EXCEL trend line fits to the densities of water and steam in _{c} = 374˚C). The percolation transitions define the limits of existence of the supercritical steam and water phases on the Gibbs density surface.

Both percolation transitions can be observed on all the supercritical isotherms as a weak discontinuity in the slope of the p(ρ) isotherms. For temperatures up to 800˚C, all the isotherms show three distinctly different

regions of behavior for the state function p(ρ)_{T}. In the steam phase ρ < ρ_{PB} the rigidity (dp/dρ)_{T} decreases with density, in the mesophase ρ_{PB} < ρ < ρ_{PA}(dp/dρ)_{T}_{ }is constant, and for water ρ > ρ_{PA} and (dp/dρ)_{T} increases with density. The slope of the isotherm in the mesophase is zero at the critical temperature. This fundamental extensive state function, called “rigidity”, is the work required to change the density reversibly, with the dimensions of a molar energy, hence equal to the change in Gibbs energy (G) with density

It will eventually explain all the mysterious maxima in

The reason for linearity of p(ρ) in the mesophase, can be understood if the percolation loci, although lacking a definition at the molecular level, are defined phenomenologically similar to subcritical spinodal lines. The limit of existence of the gas phase (PB) corresponds to the density (ρ_{PB}) at which the first system-spanning cluster of “occupied sites” or stable liquid state “nucleite” spontaneously appears in equilibrium configurations. The chemical potential of the small cluster species in the gas phase, becomes equal to the chemical potential of a liquid “droplet” at a higher pressure. At this density the gas phase ceases to exist as such, and for densities above ρ_{PB}, i.e. in the mesophase, the fluid is essentially a homogeneous mixture of gaseous molecules, i.e. dimers, trimers and small clusters, etc. plus liquid-like macro clusters. There is a bimodal distribution of cluster sizes [

Likewise, as the liquid state density is reduced across PA, from above ρ_{PA} to below ρ_{PA}, the “unoccupied holes” spontaneously coalesce to form system-spanning “holes” of gas which can contain single molecules, or small clusters, which will have the same chemical potential as the liquid, albeit at the lower gas phase pressure. Thus, in this region there are larger fluctuations in energy, pressure and density that explains the maxima loci of C_{p}, α_{p} and κ_{T} in

The discontinuities in (dp/dρ)_{T}, and linearity of (dp/dρ)_{T} in the mesophase region, seen in _{c} and PB will resemble a liquid- in-gas “mist” with wide polydispersity of the micro-droplets, whereas in the vicinity of PA a colloid with water being the continuous phase, and steam dispersed rather like small bubbles as in a gas-in-liquid “foam”. The different refractive index and size distribution of a dispersed phase give rise to a supercritical opalescence, or white appearance, known in colloid science as Tyndall scattering [

For temperatures below T_{c}, PA and PB become the metastable spinodals of liquid to gas, and gas to liquid, spontaneous nucleation loci respectively. Both above and below the critical divide at T_{c}, PA and PB represent the limits of existence of the liquid water and gaseous steam phases respectively.

Another feature of second or higher-order phase transitions, which is consistent with thermodynamic properties, is the non-existence of supercritical metastability. Whereas first-order phase transitions on Gibbs density surfaces are characterized by a possible metastable existence of either phase on both sides of the transition point of equal chemical potential, there can be no metastability across a second-order phase transition. The reason being since, in a higher-order transition, (dp/dρ)_{T} is the same for both phases at the point of the transition it follows from equation (1) that the phase with the higher rigidity would have a lower chemical potential, and hence be the more stable, on both sides of the transition on the Gibbs density surface. Since this is a thermodynamic impossibility, PA and PB loci are minima and maxima density limits for the existence of the liquid and gas phases, respectively.

Supercritical isotherms for experimental ρ(p,T) data in the near-critical mesophase region can be used to obtain independent critical parameters without resort to experimental data in the two-phase region. In this section we perform the same analysis that resulted in the determination or revised critical properties of argon, but the present exercise is less accurate as there are fewer data points on each isotherm in the tables used [

Every supercritical isotherm, shown for example in

where R_{m} is the constant isothermal rigidity (dp/dρ)_{T} within the mesophase. It is independent of density, and as shown in _{c}) and critical pressure (p_{c}), respectively, and also to obtain the pressures and densities, of the percolation loci PB and PA. For densities both below PB, and densities above PA, on every isotherm, the departures in pressure p = (p‒p_{m})_{T,}_{ }are found to be quadratic in density. A plot of density against the modulus |(p ‒ p_{m})|^{1/2} gives a straight line that interpolates at ∆p = 0 to obtain the percolation densities at PA and PB and, hence also the pressures using Equation (2), all as given in

The present thermodynamic description of criticality has given rise to slightly different values for the critical constants than those previously reported in the case of argon [_{c}, which is then obtainable from the experimental measurements of the rigidity constants in the mesoscale region by extrapolation of the slopes of the isotherms against temperature to zero (_{c} by this procedure (374.66) is slightly higher than previous literature values (given as 373.746˚C [

Likewise, a critical pressure can then also be obtained from the supercritical isotherms by plotting the intercept of the linear plots (p_{o} in Equation (2)) against T ‒ T_{c}, and interpolating to T ‒ T_{c} = 0 (

The critical pressure we obtain by this procedure, p_{c} = 221.7 ± 1 bar lies within the uncertainty of the literature results [^{−3} and the highest vapor mass density they can observe near T_{c} is 220 kg∙m^{−3} as seen in ^{−3}. This is essentially the mean value obtained from experimental saturated vapor pressures in the coexistence region, and interpolation using the LRD, historically referred to as a “critical density” [

The coexisting gas and liquid densities at the critical temperature and pressure can also be conveniently and ac- curately obtained from the pressure-dependence of the supercritical percolation loci (

isotherm (˚C) | slope (J/kg) | intercept (bar) | density PB | pressure PB | density PA | pressure PA |
---|---|---|---|---|---|---|

800 | 4.14579 | 43.28 | 162.2 | 715.8 | 230.6 | 999.5 |

600 | 2.25628 | 132.9 | 222.6 | 635.1 | 310.0 | 799.7 |

500 | 1.21997 | 186.6 | 241.3 | 480.9 | 349.5 | 572.8 |

450 | 0.72309 | 204.2 | 239.1 | 377.1 | 364.0 | 444.8 |

420 | 0.41230 | 216.2 | 242.9 | 316.4 | 373.4 | 355.2 |

410 | 0.32650 | 218.1 | 249.9 | 299.6 | 358.7 | 335.2 |

400 | 0.22003 | 222.2 | 254.2 | 278.1 | 357.7 | 300.9 |

390 | 0.13225 | 222.0 | 247.8 | 254.8 | 356.0 | 269.1 |

380 | 0.05756 | 217.9 | 238.1 | 231.7 | 370.2 | 239.3 |

374.7 | 0 | 221.2 | 234.3 | 221.2 | 370.9 | 221.2 |

densities in the range close to T_{c} are linear functions of pressure so that the coexisting densities can be obtained by plotting the percolation densities of gas (PB) and liquid (PA) against (p-p_{c)} and interpolating to zero as seen in

The data for the percolation locican now are used to construct a preliminary phase diagram of water and steam. _{p}, α_{p} and κ_{T} lie broadly in the region of the supercritical mesophase bounded by the percolation loci.

The next observation from

predicted existence was based upon a hypothesis regarding the form of the equation of state. Now, we find that two different equations of state are required for water and steam at all temperatures. Also plotted on _{c} and continue the below T_{c} to become the spinodals. This behavior was also observed for liquid argon [

At T_{c}, PA and PB cross the critical coexistence line, to become subcritical limits of existence of the metastable compressed gas (p > p_{sat}) and expanded liquid phases (p < p_{sat}), respectively. The two percolation loci define the limits of existence of gas and liquid phases, not only above T_{c}, i.e. in the supercritical meso phase, but also below T_{c}, in the two-phase water-steam coexistence region. This is consistent with a phenomenological definition of PA and PB when liquid and gas have the same values of the rigidity (dp/dρ)_{T} on the same isotherm whereupon (d_{2}p/dρ^{2})_{T} = 0 at both PA and PB, albeit at different pressures. At this point, there is also zero surface tension between water and steam, and consequently, within the mesophase, there is no barrier to spontaneous nucleation of steam from water (at PA) or water from steam (at PB).

None of the three supercritical maxima lines in _{T} on the p, T plane is quite close to PA. A better understanding of lines of maxima in these properties can perhaps be seen from the Gibbs density surface where the data points for PA and PB are plotted in _{p}. α_{p} and κ_{T}, at constant temperature, as functions of density, are obtained from the tables [_{p} maximum, which could be Bernal’s “hypercritical line” [_{p} is associated with energy fluctuations, it could represent the colloidal mesophase inversion line when the dispersed phase changes from steam to water. Curiously, both α_{p} and κ_{T} also stem from the central mean density at T_{c}, and both cease even to exhibit a maximum at all nearby the points at which they intersect the two percolations lines PA and PB, respectively.

At T_{c}, the rigidity difference between liquid and gas is zero, hence the surface tension should go to zero at a finite density difference [

In this preliminary investigation of supercritical water and steam, we have seen that the principle difference between liquid and gas is not just the density. Moreover the density difference does not disappear at the critical temperature. Rather, the structures of liquid and gas phases are fundamentally different in their description. Steam is a spatial distribution of small clusters of water molecules. Every species will have the same chemical

potential with monomers being the most probable, then dimers, trimers, etc. with ever decreasing probability, to balance the decreasing entropy with the increased enthalpy of association. These hydrogen bonded molecular clusters are occupied “sites” distributed in one large volume-spanning “hole”. Water, on the other hand is onelarge volume-spanning cluster, i.e. a network, of hydrogen bonded molecular “sites”, with a distribution within it of available “holes”. The foregoing analysis shows that this distinction between water and steam extends to supercritical temperatures.

It has been deduced from statistical thermodynamics [_{T} i.e. the work required to increase the density at constant temperature.

The relative slope of the rigidity functions is the distinguishing thermodynamic property between gases and liquids. For steam it is seen to decrease quadratically with density, and for water it increases quadratically. This can be explained by the exact statistical mechanical relationship between the rigidity function, and fluctuations in the molecular number density (N) in a thermodynamically large fixed sub volume (V) of the equilibrium fluid. Infinite snapshots of such a dynamical system in statistical thermodynamics are a grand canonical ensemble. From the statistical properties of this ensemble, it is proven that the rigidity function is exactly related to the change in Gibbs chemical potential (μ) with number density (N) according to

and from the thermodynamic theory of fluctuations [

where R is the gas constant, T is the absolute temperature and V is a molar volume. For the ideal gas the variance in Equation (4) is 1, and it will approach zero for an amorphous ground state [

amorphous ground state. The symmetry seen in

In conclusion, the previous conjecture that the present description of liquid-gas criticality, already verified in the case of argon [

We would like to thank the Scientific Committee of the 10th Asian Thermophysical Conference in South Korea Oct.1-5 2013 for the invitation to attend and present this paper. The presentation was well-received by the conference participants, with encouraging discussion and positive feedback. This manuscript was submitted to the International Journal of Thermophysics (IJTP) website in accord with ATPC guidelines for consideration by University of Korea Professor Jeong Won Kang as IJTP Editor of ATPC 2013 Proceedings, i.e. a member of the ATPC Scientific Technical Committee who invited the paper and attended the conference. On 29th Nov. 2013 I received E-mail from IJTP administrator Dr. D. Friend (of NIST, Boulder, USA), who had not attended ATPC-2013, and who was PRL reviewer A in a previous debate [

We therefore also thank the Editors of Natural Science, who provide a valuable new forum for many areas of physical sciences, by publishing open debate with open access, where advancement of scientific ideas may have been impeded by prejudice and bigotry of the peer review processes in conventional journals, and in the IJTP case, by the absence of any appeal procedure.

Reviews of manuscript “Gibbs density surface of water and steam”:

submitted (28-9-13) International Journal of Thermophysics:

reviewing Editor: Dr. D. Friend (NIST: Boulder CO)

Reviewer 1: (Prof. J. V. Sengers, University of Maryland wrote:)

This is a revolutionary paper challenging both the famous principle of Van der Waals concerning the continuity of the gaseous and liquid state and the modern theory of critical phenomena in terms of renormalization-theory and critical-point universality as applied to the vapor-liquid critical point in fluids. Basically, this paper resurrects an old picture in terms of liquidons and gasons, reviewed by Levelt Sengers [

First, I object to using the IAPWS-97 Industrial Formulation for the Thermodynamic Properties of Water and Steam [Ref. 11 in manuscript] to draw any conclusions about the thermodynamic critical behavior of water and steam. All figures appear to be based on IAPWS-97. However, IAPWS-97 is an industrial formulation developed for high computational speeds, but this equation does not even claim to be qualitatively reliable in the critical region. For example, from

My second problem, which is a more fundamental one and which also pertains to a previous article of the author [Ref. 9 in manuscript], is that the author completely ignores the experimental evidence in support of critical-point universality and of the resulting critical behavior of thermodynamic properties in terms of critical power laws [_{V} diverges with a critical exponent α = 0.1105 ± 0.027 in excellent agreement with the theoretical value α = 0.110 ± 0.003 for Ising-like systems [_{V}. The author claims in a previous publication [Ref. 9 in manuscript] that his theory is in agreement with the experimentally observed thermodynamic behavior of argon in the critical region. However, I have not found any comparison in that paper with the experimental data for the isochoric heat capacity of argon in the critical region [

As long as the author cannot demonstrate that his theory implies a similar divergence of C_{V} as observed experimentally, I do not think that the manuscript should be accepted for publication.

This manuscript proposes a major reinterpretation of the fluid region phase diagram of water, and presumably all simple fluid systems. Significantly, there is no critical point in this new description. The vapor and liquid states extend above the critical temperature, bounding a supercritical “mesophase” with colloid-like characteristics. This alternative view of the phase diagram is supported quite convincingly by the figures presented in the manuscript, and does appear to offer a plausible alternative explanation to the usual one involving the critical point which could explain a number of unresolved experimental observations (well summarized in the introduction of the manuscript). Which one is correct still needs to be established, presumably by further experiments and maybe theory. It is important that this new perspective is presented to the scientific community so that it will spur on such investigations. I support publication of the manuscript, which is well written.

While challenges to conventional interpretations of physical phenomena are certainly worth considering, their threshold for publication must be set reasonably high in terms of an ability to interpret experimental results, predict new phenomena, or simplify our understanding of the physical system. This paper challenges about 150 years of theory and experiment, and offers very little in return.

The premise of the approach seems to stem from an interpretation of

Similarly the statement that there is “no universality” contradicts most evidence (whatever “spans the dimensions” might mean). The statement about the mesophase above T_{c} is also given without evidence.

I can see no reason to re-fit tabulated EOS results to get new liquid and vapor saturation lines. (Evidently experimental data were again not considered.) The equations (in

While “hole” models of fluids have been published in numerous papers, the lack of experimental evidence or predictive powers have limited there utility: nothing changes in the current version. As noted, there is “no such well-defined distance to define PA and PB.” The argument about a “steep [??] hydrogen bond” seems irrelevant for a model that is supposed to be common for other (all?) fluids: presumable this argument was not used for argon, and such ad hoc explanations do not replace the universality of critical phenomena. On p. 5, lines 43 - 44, it is not clear what is meant by “change in dp/drho|T” and no precise definition is given.

The conjecture of linearity in the mesophase seems to be made without experimental support. The discussion of clusters is made without experimental support. The spontaneous coalescence is described without experimental support. The critical opalescence has been explained in a conventional critical region theory―what is added here?

The calculations of

In summary, particle-hole theories can be useful in statistical models (such as lattice models), but they seem very limited for interpreting existing experimental data. This paper presents nothing that contradicts existing approaches to criticality, and provides no experimental based evidence for the model under consideration. There seems to be no reason to adjust the critical values of water, nor reject the concept of critical behavior.

Sengers: “This is a revolutionary paper challenging both the famous principle of Van der Waals concerning the continuity of the gaseous and liquid state and the modern theory of critical phenomena in terms of renormalization-theory and critical-point universality as applied to the vapor-liquid critical point in fluids.”

Author: In fact this particular manuscript on water and steam is not so “revolutionary”. There are already several (5+ pipeline) refereed and published papers of this new science of criticality with all the supporting evidence, including references [

Sengers: “Basically, this paper resurrects an old picture in terms of liquidons and gasons, reviewed by Levelt Sengers [

Author: This hitherto obscure historical review by Levelt-Sengers, is interesting, and indeed very relevant, but it is not an objective review of the science. It is politically prejudiced throughout, and highly subjective in its conclusions, especially in favour of the Dutch school and van der Waals hypothesis. The Levelt-Sengers review, interestingly, reveals the names of several early scientists expressed alternative descriptions or reservations of the van der Waals “continuity of liquid and gas” hypothesis; including Faraday, Mendeleev, Ramsay, Gouy, Mathias and Cailletet.

For Sengers to describe the new science of liquid-gas criticality and the supercritical mesophase as a “resurrection” of this nebulous early nineteenth-century idea is quite a stretch. There is no definition of “liquidons and gasons”; and didn’t even know this paper existed until he sent it to me by e-mail last September.

Sengers: “First, I object to using the IAPWS-97 Industrial Formulation for the Thermodynamic Properties of Water and Steam [Ref. 11 in manuscript] to draw any conclusions about the thermodynamic critical behavior of water and steam. All figures appear to be based on IAPWS-97. However, IAPWS-97 is an industrial formulation developed for high computational speeds, but this equation does not even claim to be qualitatively reliable in the critical region. For example, from

Author: Use of IAPWS-95 rather than-97may improve the precision, but would not change the science. However, I do agree with Sengers that both the IAPWS-97 formulation of Properties of Water and Steam (which I use), and also the IAPWS-95 formulation, cannot be accurate in the immediate vicinity of the critical temperature. What Sengers fails to recognize, however, is the reason for the inadequacy of these modern steam tables in the vicinity of T_{c}. It is precisely because the functional forms used for the parameterization of the data, in both IAPWS-95 and IAPWS-97 wrongly assume the LRD in conjunction with existence of a critical point singularity on the Gibbs density surface, that is the cause of the inadequacy. This interesting point can be included. Indeed, it is a good reason for publication since the new science explains the root of the problem!

Sengers: My second problem, which is a more fundamental one and which also pertains to a previous article of the author [Ref. 9 in manuscript], is that the author completely ignores the experimental evidence in support of critical-point universality and of the resulting critical behavior of thermodynamic properties in terms of critical power laws [

Author: Non-universality is discussed in references [

Sengers: “For instance, the isochoric heat capacity has been shown to diverge at the vapor-liquid critical point [

Author: It is rather unusual to go back 50 years to dig out these three papers. On looking through them, it is clear this cannot be “reproducible experimental evidence”. This perhaps explains the 40-year gap.

Both Sengers and Voronel appear to have been under a life-long misapprehension regarding the divergence of the state function C_{v} = (dU/dT)_{V}, where U is the internal energy. This property depends only upon energy fluctuations, and, unlike C_{p} and K_{T}_{, }which depend on density fluctuations, for gas-liquid equilibria, C_{v} does not diverge at the critical temperature. This misapprehension appears to have prejudiced the design, the measurement procedures, and the interpretation of the experimental results.

The paper by Haupt and Strabe (Sengers ref [_{v}(T) along any isochore in a liquid gas system. Even the hypothetical van der Waals “critical point”, for a real gas-liquid system, could not have such divergence of C_{v}.

Sengers: The most recent experiment concerning the behavior of the isochoric heat capacity of a fluid near the vapor-liquid critical point is the one reported by Haupt and Straub [_{V} diverges with a critical exponent α = 0.1105 ± 0.027 in excellent agreement with the theoretical value α = 0.110 ± 0.003 for Ising-like systems [

Author: The experiment referred to here is a rather obscure one-off experiment on SF_{6} aboard the space shuttle, and is not easily verifiable or reproducible. Again, I must say that I find it hard to understand how these authors “found” a divergence in a property (C_{v}) that does not, and could not, diverge at any point on its Gibbs T-density surface of liquid-gas equilibria. Moreover, C_{v} would not diverge at the critical temperature even in the van der Waals hypothesis of a singular critical density were true. It is only those properties that depend on a divergence of density fluctuations, arising from the thermodynamic condition (dp/dV)_{T} ® 0 that diverge, e.g. isothermal compressibility or isobaric heat capacity. All isochoric reversible processes on Gibbs T-density surface, e.g. (dU/dT)v (= C_{v}) are continuous without divergences, basically because there is no latent heat associated with expansion in reversible isochoric processes.

Any comparison, or coincidence, between whatever it is that these authors actually measured, and the published value of an exponent for Ising models is therefore not relevant. There is certainly no valid reason for rejection. Quite the opposite, there is now an overwhelming case to open up this debate.

Sengers: Unfortunately, the author is silent about the critical behavior of C_{V}. The author claims in a previous publication [Ref. 9 in manuscript] that his theory is in agreement with the experimentally observed thermodynamic behavior of argon in the critical region. However, I have not found any comparison in that paper with the experimental data for the isochoric heat capacity of argon in the critical region [

Author: Reason for the “silence” is that there isn’t much of interest to say about C_{v}; unlike C_{p} there is no divergence. Along any isochore that traverses within the two-phase region at T_{c}C_{v} is everywhere well-behaved. The science of liquid-gas criticality regarding C_{v} is as follows.

For all temperatures below T_{c}, the thermodynamic state function C_{v} = (dU/dT)v where U is the internal energy is the mean of heat capacities of coexisting liquid per mole of liquid and of the gas per mole of gas as given by the Lever rule of apportionment. At the critical temperature this reaches a maximum value but, unlike the heat capacity at constant pressure (C_{p}), or the isothermal compressibility, C_{v} does not diverge to infinity. For T > T_{c}, along any isochore, C_{v} decreases monotonically from its maximum value at T_{c}. There is no discontinuity in C_{v} at T_{c} with any similarity to the singularities observed in Ising lattice gases. In this respect, C_{v} is a rather uninteresting property of a thermodynamic liquid-gas system. Lattice gases are not thermodynamic systems in the Gibbs definition, i.e. of being able to change energy by reversible exchange of both heat and work, because of the configurational constraints, they can exhibit energy fluctuations that behave quite differently.

Sengers: As long as the author cannot demonstrate that his theory implies a similar divergence of C_{V} as observed experimentally, I do not think that the manuscript should be accepted for publication.

Author: Sengers doesn’t seem to understand that there is no new “theory” here. What has been proposed is an alternative phenomenological description of the thermodynamic phase limits of liquids and gases on the Gibbs density surface that is consistent with Gibbs phase rule, and not based upon any hypothesis (such as van der Waals) or various other ideas involving mean-field theories and universality.

What can be demonstrated is the invalidity of Sengers principle reason for rejection using literature experimental data for the isochoric heat capacity of water and steam (see graph below). The literature value of T_{c} for water is 373.95˚C ± 0.02˚C. On the plot below is given the C_{v} values for the near-critical/supercritical isotherm (small black circles) 374.00˚C (i.e. Kelvin T/T_{c} = 1.0001), alongside the two values for saturated steam (red circle) and saturated water (blue circle) in coexistence at the subcritical 2-phase region temperature of 373.80˚ (T/T_{c} = 0.9996). There is no divergence of C_{v}!

By contrast, the values of C_{p}, which does diverge at T_{c}, for steam and water at 373.8˚C are 51,520 J/(molK) and 37,000 J/(molK) respectively, i.e. of the Order 1000 times Cv!

Friend: “While challenges to conventional interpretations of physical phenomena are certainly worth considering, their threshold for publication must be set reasonably high in terms of an ability to interpret experimental results, predict new phenomena, or simplify our understanding of the physical system. This paper challenges about 150 years of theory and experiment, and offers very little in return.”

Author: The paper that first challenged “150 years of theory and experiment”, reference [^{2} and accepted by the physics community for 140 years ^{4,23−27}, does not exist as such. There is no universality of critical phenomena that spans the dimensions. ^{4,32} Rather, in the V(T, p) surface, we find a line of uniform chemical potential at the critical temperature and pressure, above which there exists a supercritical meso-phase bounded by the second-order percolation transitions, and below which there is the familiar liquid-vapor coexistence region”.

What this manuscript “offers in return” is a vindication of that early conclusion which was based upon percolation transitions in model hard-sphere and square-well model systems. Here we now see, the conclusion is entirely consistent with all known thermodynamic properties of water and steam.

Friend: “The premise of the approach seems to stem from an interpretation of

Author: The IAWPS-97 water steam tables are themselves the “experimental evidence” in this case. Here is some more experimental evidence, with a simple question for Friend? Take a look at all the raw experimental data in

Why are there no experimental density data for saturated gas density between ρ/ρ_{c} = 0.8 and 1.0? Why are there no liquid state density points between ρ/ρ_{c} =1.0 and 1.2? Is this not the ultimate simple experimental evidence?

Friend: In the first full paragraph of p. 3, it is stated that the “5-figure precision density data” indicates that Wilson’s critical point “does not accord with experiment.” No experimental evidence is given―only results from an EOS, which, by construction conforms to some model behavior and some sets of data.

Author: _{c} (i.e. 0.1 K). Use of original experimental data points rather than IAPWS-97 parameterization, which is not model-based) would give exactly the same graph.

Friend: Similarly the statement that there is “no universality” contradicts most evidence (whatever “spans the dimensions” might mean).

Author: The statement “no universality” indeed contradicts many theoretical physicists publications on the subject, but theory is not “evidence”. The statement does not contradict the experimental thermodynamic property investigations of liquid-gas equilibria and liquid-gas criticality going back 150 years to the first measurements on CO_{2} by Andrews.

“Spanning the dimensions” refers to a widely held misapprehension that has arisen from the exact solution for a partition function showing a similar critical singularity in 1, 2 and 3 dimensions. (see: M. Kac, G.E. Uhlenbeck and J. Hammer “On the van der Waals theory of the vapor-liquid equilibrium” (I) J., Math. Phys., 4 216 (1963) (II) ibid 4 2239 (1963) (III) ibid 5 60 (1964).) for Kac-van der Waals mean-field model (reference 32 in reference 7 of my manuscript). The Kac model is a state (temperature)-dependent effective pair potential, and as such, it does not correspond to any model molecular Hamiltonian and is therefore inapplicable to real molecules that comprise Gibbs thermodynamic systems. The Kac model effective-pair potential artifact has no real molecular counterpart, and is not applicable to real molecular liquid-gas thermodynamic equilibria.

Friend: The statement about the mesophase above T_{c} is also given without evidence.

Author: The nature of the mesophase is evidenced in references [

Friend: I can see no reason to re-fit tabulated EOS results to get new liquid and vapor saturation lines. (Evidently experimental data were again not considered.) The equations (in

Author: Yes, of course they contradict most interpretations etc…that’s precisely the point of my manuscript. I’ll repeat it again: the thermodynamic equilibrium data for water-steam system is inconsistent with previous interpretations based upon van der Waals hypothesis, and entirely consistent with two-phase coexistence up to a critical dividing line with a 20% density gap at T_{c}.

Friend: _{T} would be the same with and without the conventional critical point, so not much is added in this discussion.

Author: The behavior of (dp/dρ)_{T} would be quite similar with and without the conventional critical point, but definitely not “the same”. All we are showing here is that the data for water and steam is entirely consistent with the existence of the supercritical mesophase in which, over a wide range of density, (d_{2}p/dρ^{2})_{T} = 0 to within the 5-figure precision uncertainty.

Friend: While “hole” models of fluids have been published in numerous papers, the lack of experimental evidence or predictive powers have limited there utility: nothing changes in the current version. As noted, there is “no such well-defined distance to define PA and PB.” The argument about a “steep [??] hydrogen bond” seems irrelevant for a model that is supposed to be common for other (all?) fluids: presumable this argument was not used for argon, and such ad hoc explanations do not replace the universality of critical phenomena.

Author: This manuscript has nothing to do with the “hole models” referred to.

The “hole volume” in the liquid state referred to in this manuscript is not a “model”; it is actually a thermodynamic state function chemical potential exp (−μ/kT) = V_{hole}/V

where V_{hole}/V is a normalized probability of inserting the additional molecule anywhere which implies “holes”-but there is no “molecular model”. Presumption of “the universality of critical phenomena” is evidence of prejudice.

Friend: On p. 5, lines 43 - 44, it is not clear what is meant by “change in (dp/dρ)_{T}” and no precise definition is given.

Author: Change in rigidity (dp/dρ)_{T} with density are of course manifested in the second derivative (d_{2}p/dρ^{2})_{T}_{. }The following changes or inequalities define the percolation loci thermodynamically along any isotherm

Friend: There is no explanation of why portions of the isotherms are assumed to have constant slopes, nor why the transitions lines are assumed to be at constant density in

Author: The constant slopes is not an assumption; the data points, when plotted on an EXCEL sheet, give linear trendlines with a regression coefficient of 0.999…etc. with remarkable consistency in the mesophase for every model system or real liquid that we have looked at, including the present water-steam.PA and PB in _{c}.

Reviewer 3: The conjecture of linearity in the mesophase seems to be made without experimental support.

Author: The linearity is found in all the model square-well model fluids, from sticky spheres to mean field limit, in all Lennard-Jones fluids, in liquid argon, and now, here, in water-steam. It is not a conjecture. It is an experimental observation!

Friend: The discussion of clusters is made without experimental support. The spontaneous coalescence is described without experimental support.

Author: Whilst presently there is not yet real experimental verification of the detailed nature of the clustering at the molecular level, computer simulations in reference 8 and 9 of the manuscript, and also in another recently published paper (D. M. Heyes and L. V. Woodcock , “Critical and Supercritical Properties of Lennar-Jones Fluids” Fluid Phase Equilibria 356 351 - 308 (2013) show that the distribution of clusters in the mesophase, when a cluster is defined by a distance of maximum attractive force, is bi-modal. There are system-spanning liquid clusters or Order N, and small gas clusters of Order 1 to n (where n is a small number Order 10).

Friend: The critical opalescence has been explained in a conventional critical region theory--what is added here?

Author: van der Waaals hypothesis of the supercritical region around the mean critical density as a homogeneous fluid, the only possible origin of fluctuations that could give rise to light scattering is molecular Brownian motion as described by Einstein 100 years ago and later Smoluchowsky. What is suggested here is an alternative and more plausible explanation based on the colloidal nature of the mesophase. All otherwise transparent colloidal systems, e.g. clouds or foams, scatter white light when one phase is dispersed with droplets of a different refractive index to the continuous phase. It is called Tyndall scattering.

Friend: The calculations of

Author: Yes, I agree here. The paper could be improved if I use the equations rather than the tabulations and indeed also evaluate the derivatives if at all possible. I could do this to improve the paper, and also the precision of PB and PA in the Figures, but it will not change any fundamental conclusions.

Friend: In summary, particle-hole theories can be useful in statistical models (such as lattice models), but they seem very limited for interpreting existing experimental data.

Author: I don’t disagree with this statement. This manuscript, however, is not about particle-hole model theories.

Friend: This paper presents nothing that contradicts existing approaches to criticality, and provides no experimental based evidence for the model under consideration.

Author: There is no “model” under consideration. It is exactly the opposite, all thermodynamic data used in this manuscript is extracted from precise numerical parameterization of validated literature experimental thermodynamic properties of water and steam. This original thermodynamic data is entirely consistent with the non-existence of a “critical point” on the Gibbs density surface.

Friend: There seems to be no reason to adjust the critical values of water, nor reject the concept of critical behavior.

Author: Except for “the critical density”, which does not exist, I agree that the data used here in the water- steam supercritical mesophase needs sharpening up before an accurate (5-figure) revision may be possible. This preliminary manuscript however, now consolidates the many compelling reasons, to abandon the concept of a van der Waals critical singularity, and brings new insight into the supercritical region between gas and liquid state limits.