Thermodynamic Fit Functions of the Two-Phase Fluid and Critical Exponents

Two-phase fluid properties such as entropy, internal energy, and heat capacity are given by thermodynamically defined fit functions. Each fit function is expressed as a temperature function in terms of a power series expansion about the critical point. The leading term with the critical exponent dominates the temperature variation between the critical and triple points. With β being introduced as the critical exponent for the difference between liquid and vapor densities, it is shown that the critical exponent of each fit function depends (if at all) on β . In particular, the critical exponent of the reciprocal heat capacity 1 c is 1 2 α β = − and those of the entropy s and internal energy u are 2β , while that of the reciprocal isothermal compressibility 1 T κ − is 1 γ = . It is thus found that in the case of the two-phase fluid the Rushbrooke equation conjectured 2 2 α β γ + + = combines the scaling laws resulting from the two relations d d c u T = and dln d T p κ ρ = . In the context with c , the second temperature derivatives of the chemical potential μ and vapor pressure p are investigated. As the critical point is approached, 2 T 2 d d − μ diverges as c , while 2 p T 2 d d converges to a finite limit. This is explicitly pointed out for the twophase fluid, water (with 0.3155 β = ). The positive and almost vanishing internal energy of the one-phase fluid at temperatures above and close to the critical point causes conditions for large long-wavelength density fluctuations, which are observed as critical opalescence. For negative values of the internal energy, i.e. the two-phase fluid below the critical point, there are only microscopic density fluctuations. Similar critical phenomena occur when cooling a dilute gas to its Bose-Einstein condensate.


T
κ − is 1 γ = .It is thus found that in the case of the two-phase fluid the Rushbrooke equation conjectured

Introduction
An essential property of matter is its structure, i.e. the distribution of its constituents in space and time as governed by inter-particle forces [1].We are concerned here with the electrically and magnetically neutral single-component gas under steady-state conditions, which are thermodynamically defined in the immediate vicinity of the critical point and below it.
Andrews's discovery of critical opalescence in carbon dioxide in 1869 stimulated numerous investigations of critical phenomena.The experimental observations on fluids show that a colorless gas in a narrow temperature range ( ) T T T − − ≈ above its critical temperature c T suddenly becomes opaque and changes color and, at c T T ≤ , exists as a two-phase fluid of different densities in volumes that are sharply separated by an interface surface a few molecular layers thick.The endeavor to explain theoretically the observations constituted a huge challenge.The course of a century then saw the development of the familiar phenomenological theories of a van der Waals gas, of the stable and unstable thermodynamic equilibrium formulated by Gibbs, of the correlation of fluctuations, of the scaling laws, including the hierarchical reference theory (renormalization group techniques), and of the Monte Carlo computer methods (see [2]).
An insight into the nature of a fluid in the critical region is afforded by < in the sub-volume v V .This gas in thermodynamic equilibrium existing in two phases is called a saturated fluid.As thermodynamic theory teaches, as the only independent variable of the saturated fluid that can be chosen is the saturation temperature T , since the other field variables possible, viz.vapor pressure p and chemical potential µ are unique functions of T .
Every thermodynamic quantity of the saturated fluid,

{ }
free energy, entropy, internal energy, heat capacity X F S U C = , can thus be represented as a function of T .The absolute values both of the fluid, X , and of the fluid phases, vapor, v X , and condensate (liquid, solid), l X , are proportional to the mass in the volume considered.As extensive quantities they have additive properties, i.e. they satisfy the equations , , , , , , .
The mass-specific quantities x X M = and , , , v l v l v l x X M = contain the complete thermodynamic information on the fluid state [3].The difference ( ) gives the difference of the thermodynamic properties of x in the volumes v V and l V and is called an order parameter.The quantity x can be expressed by the quantities , v l x , , v l v , and v : .
The quotient ( ) ( ) ( ) is a function of the vapor pressure ( ) p T , and the quotient a function of the chemical potential ( ) The vapor pressure is a positive, convexly curved, monotonically increasing function of the temperature and the chemical potential a negative, concavely curved, monotonically decreasing function of the temperature.Open questions on the thermodynamic properties of p and µ , in particular in the critical region, are dealt with in Sections 3 and 7.
Equation (5) yield, for example, ( ) Relations (1)-( 5) represent the Gibbs equations for calculating x , v x , and l x of the saturated fluid.One might think that the job of setting up generally valid fit formulae, i.e. applicable to every two-phase fluid has already been done.But this is not so.The literature yields, e.g. for water, only formulae for industrial use ([4]- [6]) which do not correctly give the physical picture of the two-phase fluid in the critical region.
The objective of this study is to represent the quantities v x , l x , and x in the region between the triple point and the critical point as thermodynamic fit functions dependent on the temperature (Sections 4-9).The representation of a fit function by an order parameter expanded around the critical point is based on the knowledge of the behavior of the fluid in the critical region.

Thermodynamics of Critical Phenomena
The thermodynamical physics of critical phenomena above and below the critical point is extensively treated in the literature (e.g.[2] [7]- [10]).Critical phenomena occur under the natural boundary condition of the vanishing value of the internal energy, 0 U ≈ [11].In the immediate vicinity of the critical point the one-phase fluid is in unstable equilibrium on transition to the two-phase fluid, which is then in stable equilibrium.
Thermodynamics describes the macroscopic state of the fluid by means of the quantities M , V , S , and U and thus cannot delve into the microscopic processes actually occurring in the particle interactions taking place in the fluid.The effect of attractive and repulsive forces among interacting particles on the internal energy U is that U has negative sign for fluid temperatures below c T and positive sign above.U can be treated as the sum of two energy contributions, viz. the potential energy pot U , whose gradient yields the attractive forces, and the thermal energy th U , which is assigned to the sum of kinetic, vibrational and rotational particle energies.As the result, Figure 2 shows for saturated water (under the same conditions as in Figure 1) the internal energy U and the estimates pot U and th U in liquid and vapor as functions of the mean particle separation, represented by the normalized density values ( ) ,525, 647 , , 1.0002,3.1056,48.362 cm . For an estimation it is taken for granted that the thermal particle energy increases proportionally as the saturation temperature T , constant The empirical value constant is given by ( ) M .This yields in the liquid phase an initially expected increase constant and on reaching the maximum value 668 [J] at about ( ) 0.902 ( ) A. Elsner

U uM =
(with 1 ) can be calculated according to Equation ( 5), it is possible to estimate the potential pot th U U U = − numerically.pot U is always negative.The resulting repulsion and attraction forces between the particles are equal and opposite at the critical point, which is expressed in Figure 2 by the fact that the curves th U , pot U , and U all are continuous there.Qualitative information about mean-field strength of forces in liquid and vapor can be obtained from A fluid state of 0 U > is characterized by an ensemble of particles freely moving in a structureless homogeneous phase.In contrast, a fluid state of 0 U < is characterized by an ensemble of particles bound in a more or less structured form as a result of particle self-organization under certain constraints, e.g.liquid, solid, and Bose-Einstein condensate (BEC) [12], each with its specific thermodynamic property.As an example, Figure 4 shows the variations in internal energy of vapor and condensates as functions of their phase-specific volumes and the saturation temperature, respectively, for the water mass of 1 [g].
The discontinuities of d d T U S = , represented as circles in Figure 5, indicate the phase transitions of bulks of different structures.States of different aggregation exhibit qualitatively different properties.Adding energy dU to the fluid at fixed temperature distributes a surplus of the one bulk phase at the expense of the other [1].Local variations of internal energy couplings between particles change the bulk structures.A structural change is thermodynamically described by an increase in internal energy dU as a result of increasing system entropy dS , i.e. by d d U S .Changes in the bulk structure occur close to absolute zero, at the triple point and the critical point.The circles in Figure 1 and Figure 5 give the phase transitions for the functions     , as in the case of a Bose-Einstein condensate, quantum mechanically by a condensate wave function.As already mentioned, correlation functions are a measure of the number of scattering centres for light in the fluid dielectric and hence of a mean value of structural density fluctuations.The strong increase and subsequent decrease of long-wavelength fluctuations in critical fluid regions cause the observed sharp increase and decrease of scattered light intensity (e.g.[13]) and are thus experimental proof of the thermodynamic zero of the internal energy of, on the one hand, a dense gas in the critical temperature region and, on the other, a dilute-gas in the 7 10 − K region.Evidently, nature associates the problem of changing the sign of U at c T with the ability of self- organization of particles interacting in ensembles.For is to compute the order-parameter probability distribution functions by means of Monte Carlo computer methods [10].Computer simulations under different thermodynamic conditions afford quite a good picture of the phase transition in the critical region, and in conjunction with the renormalization group techniques allow calculation of critical exponents.For 0 U < the Gibbs applied theory, on the other hand, yields with Equations (1)-( 5) universal orderparameter relations for calculating material-dependent critical exponents.
Van der Waals was the first to show that the normalized difference in the phase densities ( ) ρ ρ ρ − is empirically very well described by a power law of the form ( ) T T T − appears as a variable and the so-called critical exponent β characterizes the decrease of ( ) 1).Applied thermodynamics has to determine the particular critical exponent for every order parameter, ( ) Here the value of the critical exponent depends on the choice of the variable of the order parameter.It is customary to choose the temperature as variable.
When thermodynamic relations between various thermodynamic functions x are known, these will reappear in corresponding relations between their critical exponents.A well-known exponent equation for the saturated fluid (e.g.[7]) is ( ) Consequently, the task here is to repeat the calculation of α and γ and additionally of the exponent of every function, which according to Relations (3) to ( 5) is connected with the heat capacity.Since ( ) these are the functions c , f , s , u , p , µ , , v l c , and , v l v .Obviously, c is an important measurable quantity that yields information on the phase transition at c T T < .The knowledge obtained about the temperature dependence of the functions mentioned then allows relations between critical exponents to be studied, e.g. between α , β , and γ (Sections 4-9).It is found that only a single independent critical exponent is needed to characterize all order parameters, e.g.β and the others can be expressed by it.Data for the fluid selected, saturated water, are given in the Appendix.
It remains to consider some essential properties of the above-mentioned free interface surface.
If the condition for forming a free surface between the liquid and vapor phases is given, then there is an interface particle layer, which represents a new equilibrium state described by a minimum internal energy U and simultaneously a maximum entropy S .Hence formation of the free interface surface lowers the free energy F of the fluid.The relative energy contribution of an interface quantity to the respective system quantity depends on the ratio of the numbers of interacting particles in the interface volume (of the surface area A times the layer thickness

∫
) and system volume ( ) V and is therefore extremely small.Despite the smallness of the order 8 10 − and less, surface effects play a role in nature and technology, e.g. the minimization of the free interface surface.The smallness of an interface quantity shows, on the other hand, that ignoring it when studying volume properties of the fluid is completely justified.As the existence of a surface A does not change the mass M and volume V , the property of U , S , and F being extensive quantities is maintained.If the area-specific difference of the internal energy (surface energy), then the coefficient of the surface entropy increased at constant temperature, , is assigned to the function d dT σ and that of the surface free energy decreased, Equations ( 9) give U , S , and F as functions of the system variables M , V , and A and their conjugate temperature variables µ , p , and σ .In order to establish a relation between the unknown function σ and the measurable surface tension γ , it is posssible simply to identify the negative values ( ) with the free energy, i.e. ( ) ( ) ( ) It is found that this choice gives the acceptable result of interface- specific coefficients, viz.( ) ( ) , and concave curvature ( ) 6 for water.The physical significance of the negativity of σ and ( ) ( ) and the positivity of d dT σ is readily apparent since these functions multiplied by A represent area-contributions to the negative internal and free energies and positive entropy of the two-phase fluid.

Vapor Pressure
The experimental finding that the heat capacity diverges at the critical point calls for a statement on the critical behavior of the vapor pressure and chemical potential, since . What was known about the properties of these two quantities at the time when systematic investigations of the critical behavior of fluids and magnets were initiated [14]- [17] is summarized by Stanley [7] in his book (1971), Introduction to phase transitions and critical phenomena, as follows: If c is divergent, then or both will be divergent.The exponent p ϑ is introduced as a measure of the degree of divergence (if any) of the curvature of the vapor pressure curve, i.e.
( ) . The lattice-gas model gives p ϑ α = .However, for the real gas the curvatures of p and µ might both indicate divergence, so that p ϑ might differ from α .In particular, the divergence of the heat capacity of helium-4 ( 4 He) appears to be dominated by ).This attitude, however, is not accepted by all.In general, the literature provides no uniform statement on the temperature dependence of 2 2 d d p T in the critical region.The findings range from the absence of divergence, e.g. in the case of helium-3, to explicit specification of the exponent, e.g.
We shall take up the problem and show that 0 If, on the other hand, the exponent µ ϑ is introduced in order to describe by ( ) ( In the one-phase critical region ( ) c T T > a distinction must be made between conditions at constant volume, pressure, and chemical potential since ( ) ( ) . For example, the difference of the heat capacities is positive, ( ) ( )  ( ) ( ) , one gets the areaspecific internal energy ( ) and the area-specific entropy ( ) With respect to the vapor pressure, the necessary thermodynamic proof of the finiteness of 2 2 d d p T in the critical region enlists the possibility of estimating the second temperature derivative by an expression containing as highest temperature derivative the first derivative d d p T , which itself is finite at the critical point.An estimate formula is obtained as follows.The vapor pressure can be expressed by the following relations: The function φ is defined by the ratio of the evaporation energy ( ) at the phase transition.This ratio, which is always greater than 1, decreases monotonically as the temperature rises because the binding energy difference ( )  ( ) ( )

d p T p T T p T p T p p T T
The maximum value of ( ) ( )    ; this value can be calculated according to the scaling laws ( )  , further differentiation of ( )( ) The  -expansion (where ( ) The positive functions p and d d in the  -expansion are found by fitting the given vapor pressure data ( ) A fit formula of conceptually different form, based on the expression (10) for the vapor pressure, is p T (see Figure 7).This linear function is ( ) If one introduces the dimensionless variable ( ) , y is represented as a function of  as follows: and [5], dashed lines) and Equations ( 13), ( 14), and ( 19) (solid lines).The line lnp represents published vapor pressure data from [4] and [5], which are equally well reproduced by both the dashed and solid lines.The lines ln d d p T , calculated from the respective equations, also give results that are in reasonable agreement at every temperature between the triple and critical points.In contrast, the dashed and solid lines of

( ) ( )
The essential property of the function ( )  and remain finite for 0 1 ≤ <  . This satisfies the above-stated thermodynamic requirements, viz.that d d n n y  may nowhere diverge in [ ) 0,1 .Fine fitting with measured vapor pressure data p calls for a further function ( ) z  , which, of course, is likewise arbitrarily often differentiable and nowhere diverges, and which together with ( ) y  as product function y z ⋅ fits the values lnp .As a function with m fit constants k a , the finite power series ( ) can perform the task required.
The first three derivatives of ( ) ( )( ) If the n -th derivatives of y and z are denoted for short by ( ) n y and ( ) n z , the n -th derivative of ( ) With the binomial coefficients nk n B k The following fit formula for the vapor pressure is conceived such, with due allowance for Equation (10), that it does not yield any divergent higher-order temperature derivative: Each constant k a of the series z (see Equation ( 17)) is multiplied by the prefactor 16)), yielding the constant k a , which is again denoted by k a in the fit formula (19).Equation ( 19) can be used for describing the vapor pressure of every fluid.In the literature, however, one finds fit equations (e.g. in [4] [5] [22]) that are thermodynamically incorrect, because the corresponding function z contains terms with non-integer exponents k (e.g.terms such as 1 3    or 1 2   ) that lead to divergences of Formula (19) with the ten fit constants listed in the Appendix, Equation (A3), reproduces the measured data p of water with the same accuracy as that given in [4]- [6].The calculations of 2 2 d d p T in this study and, on the other hand, according to the equations published in [4] and [5] show that the differences expected occur exactly in the critical region, as seen in Figure 7. Results of calculating the water vapor pressure p according to Equations ( 13), (14), and (19) are given in the Appendix (fit functions (1)-( 3)).

Coexistence Curve
Along the coexistence curve in the critical region, the scaling laws must of course obey each two-phase equilibrium relation, e.g.0 , and at the critical temperature c T T = the limits [11]: In accordance with the defining Equation (7) for the exponent β one obtains the scaling laws of the phase- specific volumes, internal energies and entropies in relation to their critical values as follows: , , , , , .
In these formulae the plus sign refers to the vapor phase and the minus sign to the liquid phase.According to van der Waals the temperature dependences of the volumes v v and l v can be represented as series expansions about the critical value v , where the temperature expansion variable λ gives the distance to the critical point and at the same time the decrease of the difference ( ) − , which both tend to zero on approaching the critical point: ( ) The expansion variable λ is then expressed by the power function .
The temperature dependence of λ is described by ( )

Critical Exponents of the Phase-Specific Quantities v x and l x
According to Equation (3) a quotient ( ) ( ) It thus follows that
According to Landau every quantity ( ) x x − between the stable phase-limiting values v x and l x , whose difference, as described in Equation ( 24), tends to zero on approaching the critical state c x , is based on an order parameter [28].In the case of the saturated fluid the natural order parameter is the difference of the condensation and vapor masses in relation to the fluid mass, ( ) ( ) The ratios l M M and v M M can be expressed by the quotients ( ) ( ) + , which is justified in the critical region (see curves 1 and [1] in Figure 8), then ( ) The variable λ thus represents the order parameter of the mass splitting, which as natural order parameter of the two-phase fluid is also the basis of all other order parameters: For values 0 λ → , the phase-specific volumes and internal energies on the coexistence curve conform to the scaling functions, ( The values of η and β can be determined empirically by varying them until the curves of the published data .One obtains ( ) ( ) It can now be stated that the phase-specific quantities v x and l x in relation to their critical value c x scale as ( ) In contrast, the temperature derivatives  ( ) ( ) . Published data v v and l v from Ref. [5].

T p T T T T u u p T v p T T T T T u u p T v p T T T T T
The relation ( ) = , claimed by Ref. [11], is wrong and must be replaced by Equation (29).The divergence ( ) To get temperature-dependent quantities in the entire temperature range of the liquid [ ] The sum contains m constants n a which are specific to the quantity to be fitted and are calculated by the mathematical method of conjugate gradients by fitting the given data.Quantities such as the vapor functions ( ) h and the corresponding liquid functions as well as the differences ( ) x x − are expressed in terms of the fit function (30).It is found that 10 m = constants is sufficient to generate data with the numerical exactness usual in the literature.Results are stated in the Appendix, see fit functions ( 5)-( 26).

Critical Exponent of the Fluid Energy
First it is shown that the internal energy is an order parameter: Each of the following functions has the property of an order parameter since it holds that and ( ) To obtain the fluid quantity x as a function of λ , Relation ( 2) is transformed as follows: The fluid quantity x in relation to its critical value ( ) (which takes either the minimum or maximum of ( ) In the case x u = one obtains with ( ) Below c T , u is a negative function which tends to zero in the critical region as 2 λ or ( ) u T is a positive function and scales as ( ) Quantities such as the fluid functions ( ) and ( ) See Appendix, fit functions ( 27)-(33).In particular, from ( ) ( ) ( ) ( )


described by Equation ( 34) and the scaling of 2 2 d dT µ is calculated from the difference ( ) The (positive) functions ( ) , where i a are constants.See Appendix, fit functions (38)   and (39).

Critical Exponent of the Reciprocal Isothermal Compressibility of the Two-Phase Fluid
The isothermal compressibility T κ is measured as the relative change of the fluid volume when the pressure is increased, i.e. the reciprocal isothermal compressibility 1 T κ − is defined by the relations: In the one-phase region the density > , which results in a temperature increase by dT , this in turn leading to the pressure increase dp .The isothermal compressibility is therefore infinite except on the coexistence boundary when the fluid becomes homogeneous on the liquid side.Representation of the connection between vapor pressure and phase densities along the coexistence curve allows the difference ( ) , is an isothermal variation of the entropy and volume, ( ) T s v ∂ ∂ , which according to Equation ( 3) is expressed by the isothermal phase transition ( ) ( , is expressed by ( ) is the result of the temperature increase dT or the vapor pressure increase dp since ( ) ( ) Hence the effect of the isothermal and isobaric infinitesimal changes on the two-phase fluid of constant density is ( Calculating the scaling of and yields the exponent γ of the two-phase fluid: ( ) 10 for the case of water; in the critical temperature region it goes asymptotically to zero with the proportionality factor ( ) κ under one-phase conditions for water liquid and gas about the critical region are performed by Ref. [29], yielding 1.20 0.05 γ = ± .

Fit Functions for the Saturated Fluid
To describe important thermodynamic properties of the saturated fluid, fit functions of quantities are set up.A fit function presents an appropriate power series expanded about the critical point and affords the possibility of obtaining fairly reliable information concerning the thermodynamic function.Of course, the exponent of the series is defined by the thermodynamic function.Calculation of the fit constants is done by the least squares method, published and estimated data taken as a basis.Fit functions are evaluated for water data (see Appendix).The tables each (arbitrarily) list 10 fit constants, with the aid of which each temperature value of the thermodynamic quantity in the region [ ] T T can then be calculated.As an example, some evaluated functions are represented in Figure 11, viz.,

Results and Discussion
The Gibbs theory describes the macroscopic state of matter by means of the quantities M , V , S , and U  ) ( ) ( ) ( ) T according to Ref. [30].and thus cannot delve into the microscopic processes of interacting particles which give the structure of matter.Nevertheless, the theory treats the thermodynamic equilibrium matter correctly, using interrelations between the entropy, internal energy, chemical potential, pressure, and temperature.If an additional parameter, β , is introduced for characterizing the behavior of any thermodynamic quantity under the boundary condition of a reversible isobaric and isothermal transition from the condensate structure to vapor, it follows that this β reveals in all thermodynamic relations relevant to these changes.For the two-phase fluid, the established Gibbs interrelations are summarized by Equations ( 1)- (5).Equations (8) explicitly show relations between the heat capacity, free energy, entropy, internal energy, vapor pressure plus chemical potential, and phase-specific heats and volumes as functions of the temperature and, consequently, the parameter β .Then a specific quantity is formulated, according to van der Waals, in terms of a power series expanded about the critical point as a function of T and β , which is called a thermodynamic fit function.It is valid between the critical and triple points.Thermodynamic fit functions y for 40 two-phases quantities are listed in the Appendix.They are given as functions of the temperature variable ( ) allow calculation of all thermodynamic functions with absolute figures and the unique transformation of one plotted projection plane into another.
Regarding the scaling on approach to the critical point, it is found that the phase-specific volumes, entropies, and energies are dominated by the critical exponent β , and the fluid entropy and energy by 2β , while the second temperature derivative of the chemical potential by ( ) T is described by different free energy functions [14].There are large density fluctuations in the immediate vicinity of the critical point under the condition 0 U ≥ .These unstable fluctuations do not occur under the condition 0 U < and are not the subject of the Gibbs theory, which deduces the thermodynamic stable equilibrium of the macroscopic fluid from the statistics of microscopic fluctuations of its constituents at given total energy.
The subject of quantum mechanics is to investigate in detail interactions between particles.There are statistical arguments about the relative likelihood of attractive and repulsive interactions.In particular, one has to consider the possible resonant scattering of particles with different internal degrees of freedom and spinexchange and dipolar processes [12].This is all very hard to do, but it may finally provide an explanation for large density and energy fluctuations in a fluid in the vicinity of 0 U ≥ .On the other hand, atto-and femtosecond laser techniques have made it possible to follow the movements of atoms and study subpicosecond processes in solids and liquids as well as in gases (in and out of the direct light interaction region) [31].
Almost all significant results produced in the years 1960-1980 were obtained at a time when non-classical thermodynamic behavior of the free energy of the one-phase fluid near the critical point is assumed.The efforts to obtain an equation of state that describes the singularities of the heat capacity and isothermal compressibility in the critical region are to be found in, for example, Refs.[14]- [16], and [29].As is shown in this study, however, the scaling laws of the two-phase fluid can be treated on the basis of the classical thermodynamic free energy function if the equilibrium relations between the phase-specific volumes, vapor pressure, chemical potential, and heat capacity are considered.A significant result is thus confirmation of the Rushbrooke equation V and v V , respectively.Conditions for macroscopic fluctuations and hence observation of critical opalescence exist near to the critical point where Changes in the fluid structure take place at phase transitions where the gradient ( ) , 5).This is the case at the critical point of the Bose-Einstein condensate when the BEC changes to the dilute gas, at the triple point when the solid continuously melts to liquid, and at the critical point of the fluid when the two phases disappear.The first two phase transitions mentioned are associated with heat input required to alter the fluid structure and are called phase transitions of the first kind, whereas the phase transition at the critical point occurs without additional energy and is called a phase transition of the second kind.Every isothermal transition at c T T < from the condensate to the vapor state, which requires the energy input ( ) 0 v l U U − > , is a phase transition of the first kind.Adding energy to the fluid increases its entropy and thus the mobility of particles, both being equivalent to a loss of matter structure.This feature can be described by the concept of an order parameter.The "natural" order parameter of the two-phase fluid is ( ) ( ) Area-specific order parameters are the surface energy σ and surface tension γ .They describe the effect of particle interactions in the highly inhomogeneous interface layer at the free area A between the homogeneous bulks, condensate and vapor.Since A σ ⋅ is part of the fluid internal energy U , which is negative below the critical point, σ must be negative and the temperature values can be determined according to Equation (9), if ( ) ( ) context with c , the second temperature derivatives of the chemical potential µ and vapor pressure p are investigated.As the critical point is approached, positive and almost vanishing internal energy of the one-phase fluid at temperatures above and close to the critical point causes conditions for large long-wavelength density fluctuations, which are observed as critical opalescence.For negative values of the internal energy, i.e. the two-phase fluid below the critical point, there are only microscopic density fluctuations.Similar critical phenomena occur when cooling a dilute gas to its Bose-Einstein condensate.Internal Energy, Free Energy, Heat Capacity

Figure 1 ,
which for water of mass 1 M = [g] and critical density c ρ in the volume 3.1056 c V M ρ = = [cm 3 ] shows the different fluid states as a function of the temperature T .Below c T , M is distributed as condensed mass

Figure 1 .
Figure 1.Water mass 1 M = [g] of critical density c ρ in the volume 3.1056 c V M ρ = = [cm 3 ].Below the critical temperature 647.3 c T = [K] the

Figure 2 .
Figure 2. Internal energy U , estimated thermal energy th U and potential energy pot U of water mass [ ] 1 g in liquid and vapor as functions of the normalized particle separation ( ) 1 3 positive and negative regions of the fluid internal energy U are shown in Figure 3 for water in the pressure vs volume diagram.They are separated by the isotherm ( ) , c T v p above the critical pressure c p (dashed line) and the vapor pressure ( ) p T below c p and c v v > (solid line).Along the dashed line there is a continuous change in the density passing through positive and negative regions of U .The solid lines represent the vapor pressure at temperature T and are the loci of the first-order phase transition due to the jump between low-density vapor and high-density condensate.The jump is combined with a different fluid structure in each phase.
It is seen that the phase transition near absolute zero at the critical temperature of the BEC, ,BEC

Figure 3 .
Figure 3. Volume-pressure diagram of water and the regions of positive and negative internal fluid energy U .The dashed line, i.e. the isotherm ( ) , c T v p for

Figure 4 .
Figure 4. Saturated water of mass M = 1 [g]: temperature T and internal energies U , vapor U , liquid U , and ice U versus

Figure 5 .
Figure 5. Internal energy U of saturated water of mass [ ] 1 g and critical density as a function of entropy S .It holds that d d U S T = and in the twophase regime that l v U U U = + .The mass-specific quantities (dashed lines) satisfy the relations 0 l v u u u ≤ ≤ ≤ .The curves U and u are identical since M = 1 [g].The data below the critical point are obtained from fit functions given in the appendix and those above are deduced from the T, s-Diagram of Ref. [4].

2 U U −  is equal to 2 U
energy variable U  from its average value U is itself a fluctuating variable and the mean square deviation ( ) measure of the magnitude of the fluctuations[3].The energy fluctuations become enormously close to the critical point since ( )  .Together with density fluctuations, this represents an unstable equilibrium of the fluid state.One of the best ways of finding solutions for realistic particle interactions under the condition 0 U ≥

(
Rushbrooke equation), which describes the numerical relation between the exponents of the reciprocal heat capacity ( ) α , the difference in the phase-specific volumes ( ) β , and the reciprocal isothermal compressibility ( ) γ .The temperature dependences of these functions in the immediate vicinity of the critical point are defined by have become resigned to not making any statement on µ and letting g.[18]-[22] the curvature of the two-phase chemical potential, then µ ϑ α = is valid, as shown in Section 7.
while it vanishes in the case of the two-phase fluid ( ) c T T < [11].To describe phase transitions under various conditions at c T T ≥ it isnecessary to introduce some different independent critical exponents that are inter-related[10] [23]-[27].

Figure 6 .
Figure 6.Analyse of free-surface quantities.From fitted published surface tension data ( ) γ of water [4] [6] and setting γ equal to the negative area- specific free energy, i.e.
always exceeds that of the free energy difference ( ) has a welldefined value which does not essentially change in the limits ( ) 0 v l u u − → and ( ) 0 v l v v − → and can thus already be determined at a great distance from the critical point.It is finite and almost constant in the entire critical region.From constant φ =in the critical region it follows that following etimate is generally valid: n n p T increase monotonically as the temperature up to their finite critical values c p and d d n n c p T , i.e. it holds that )r r p T and r T are reference values, e.g. the boiling temperature no T [K] at atmospheric pressure 0.101325no p = [MPa].The usual representation of measured vapor pressure data ( ) p T in the form lnp versus 1 T shows that the data can be described in first approximation by the straight line through the boundary points ( )

Figure 7 .
Figure 7. Vapor pressure of water according to the IAPWS equation (Refs.[4]and[5], dashed lines) and Equations (13),(14), and (19) (solid lines).The line lnp represents published vapor pressure data from[4] and[5], which are equally well reproduced by both the dashed and solid lines.The lines ln d d p T , calculated from the respective equations, also give results that are in reasonable agreement at every temperature between the triple and critical points.In contrast, the dashed and solid lines of

p T < ∞ for 2 n
for temperatures in the vicinity of the critical point.Vapor pressure data should obey the condition d d n n ≥ .These conditions are not satisfied in the critical region by the IAPWS equation [4] [5].
the figures in the Pascal triangle, one can calculate all ( ) 1 m − derivatives of the function ( ) y z ⋅ and they are all finite in [ ) 0,1 .

.
and ζ is the exponent defined by the function under consideration, for example, The definition range [ ) 0,1 for λ and  restricts the range of values of η and ζ .For saturated fluids one has 2 η ≈ and ζ can assume a positive or negative value.The smaller a positive exponent ζ is, the deeper the approach of the function to zero.A negative ζ corresponds to a function which diverges to infinity at the critical point.An exponent 1 ζ = leads to a series expansion with no anomalous behavior as, for example, the vapor-pressure difference ( ) c p p − in Equation (13).The properties stated are also exhibited according to Equation (21) by the functions ,

(curves 2 and 3 )
and the scaling laws (straight lines[2] [3]) satisfactorily converge in the critical region as shown for water in Figure8.From Figure9it is seen that relations(26) are valid for ( )


; they diverge at the critical point where 0

Figure 8 .
Figure 8. Scaling law of the vapor and liquid order parameters for water.Critical-point data: 647.096 c T = K, 3.101 c v = cm 3 /g, 0.3155 β =

Figure 9 .
Figure 9.Comparison of vapor data v c v v (curve 1) with scaling data ( ) 1 1 λ − (curve [1]) and liquid data l c v v (curve 2) with scaling data ( ) 1 1 λ + (curve [2]) for water.The data of each of the curves (1 and [1]) and (2 and [2]) agree in the range ( ) 1 0.01 c T T − < , thus verifying the independent variable with values between 0 and 1, and λ is replaced by β τ .A fit function expanded about the critical point with the exponent β is then


to 0. The exponent of c is denoted in accordance with Equation (7) by ( the heat capacity function c is the temperature derivative of the internal energy function u , whose exponent is 2β .From the scaling

=
14].In the two-phase region, on the other hand, the vapor pressure remain constant at constant temperature, and one has to consider how a mechanical quantity such as the isothermal compressibility is to be interpreted here.Any attempt to compress the fluid of constant density 1 v ρ − = is accompanied by condensation of the vapor mass v M and release of energy as a density function ρ , whose pressure dependence d dp ρ is a well-defined negative quantity with the reciprocal value according to the Le Ch a ˆtelier-Braun princple[3].The change of state the two infinitesimal changes of states, viz.

Figure 11 .
Figure 11.Properties of water: Curves 1: 0 v l u u ≥ ≥ , curves 2: ( ) v l u u ± − , figures proposed here are specified for saturated water ( )0.3155 β =.The critical point is thermodynamically defined by the vanishing value of the internal energy, and this ( ) ( ) , , 0 c U S T T M V = = and Nernst's thus determining the divergence of the heat capacity.The second temperature derivative of the vapor pressure converges to a finite value.The difference of phase-specific quantities likewise scales as β and the sum scales as 2β .The exponents of the reciprocal heat capacity, calculated for the two-phase fluid are not in conformity with those given in the literature for the one-phase fluid.The corresponding exponents differ since the thermodynamic behavior above and below c for the two-phase fluid.It is satisfied by the van der Waals gas because 0 with β having a value smaller than 1 2 .Since the exponent of the internal energy u is 2β , the Rushbrooke equation can be interpreted as the sum of the exponents of the reciprocal heat capacity ( ) relation to those given for the real one-phase fluid in the literature, e.g.Ref.[23], surprising indeed.Such eminently different results can be explained, as mentioned above, by the statements that the scaling hypothesis, in particular the renormalization group techniques, is valid for one-phase fluids where 0 U ≥ and p V C C > , but is not suitable to represent the behavior of two-phase fluids where 0 Nature reveals the thermodynamic energy zero of a gas at the critical point through the sudden appearance of macroscopic fluctuations in a narrow temperature range ( ) c T .The critical behavior depends on particle interactions, thus the equilibrium conditions of fluid stability and instability and the size of fluctuations in the vicinity of 0 U = have to be regarded.The condition 0 U = means that the positive internal energy contribution of the thermal energy of all fluid particles, th U , balances the negative internal energy contribution of the binding energy of all particles, pot U .Under the condition pot th U U = − the interparticle forces of attraction and repulsion cancel.This leads to an unstable equilibrium of one-phase particles and hence local density fluctuations, () .If, however, the mass M in V is split up into the condensed mass l M of overbalanced bound particles with pot 0 there is no eventuality of macroscopic fluctuations, i.e. the masses l M and v M are in stable equilibrium within the subvolumes l state).In the BEC state where macroscopic density fluctuations are suppressed and alone microscopic fluctuations occur, existence of coherent matter waves.
)), but specific order parameters may be used to specify fit functions, e.g.() entropic, energetic or volumetric first-kind phase transitions (see Appendix).

⋅Figure 5
set equal to the negative temperature values of γ .As part of the fluid free energy F , lowers the free system energy when the fluid proceeds to a new equilibrium state with an interface layer.Interesting enough, the zeros of the surface energy and surface tension, viz.0 σ = and 0 γ = at c T on approaching the critical point from below, evidently display the disappearance of the free interface M V = .The thermodynamic internal energy zero of the low-density fluid is the vanishing value v U at the critical temperature of the BEC, i.e. .The internal gas energy v U is positive above ,BEC c T and negative below, when all gas particles suddenly condense into a quantum state (the quantum state of zero momentum, zero entropy and zero temperature[3].Appendix: Fit Functions for Thermodynamic Properties of a Saturated Fluid in Liquid-Vapor EquilibriumComputer program for evaluating fit functionsTemperature values of a specific thermodynamic quantity are calculated by calling the IDL-FUNCTION evaluation in an IDL-source program:The IDL-FUNCTION evaluation is specified by the label, e.g.label k = , that is relevant to the selected fit function y , the fit arrays temperature variable x of y .The general structure of an IDL evaluation function is given below.