Isobaric Critical Exponents: Test of Analyticity Against NIST Reference Data

Finite systems may undergo first or second order phase transitions under not isovolumetric but isobaric condition. The “analyticity” of a finite-system partition function has been argued to imply universal values for isobaric critical exponents, α_P, β_P, and γ_P. Here we test this prediction by analyzing NIST REFPROP data for twenty major molecules, including H₂O, CO₂, O₂, etc. We report they are consistent with the prediction for temperature range, $10^{-5}<|T/T_c-1|<10^{-3}$. For each molecule, there appears to exist a characteristic natural number, n = 2, 3, 4, 5, 6, which determines all the critical exponents for T < T_c as ${\alpha }_P={\gamma }_P=\frac{n}{n+1}$ and ${\beta }_P={\delta }^{-1}=\frac{1}{n+1}$. For the opposite T > T_c, all the fluids seem to indicate the universal value of n = 2.

1. Introduction

The celebrated paradigm for critical phenomena in modern statistical physics assumes the thermo-dynamic limit [1] with the intention of achieving non-analyticity [2, 3]. In this way, correlation length diverges at the critical point, which gives rise to the scaling theory of criticality. In order to describe the scaling, renormalization group methods were developed as for a mathematical framework (e.g., [4–12]).

However, the assertion that the thermodynamic limit for non-analyticity is a necessary condition to realize phase transitions is valid only in the isovolumetric case and does not necessarily apply to the isobaric situation [13–19]. In particular, a counterexample based on the ‘analyticity’ of a finite-system partition function has been suggested [20], which predicts its own universal values for isobaric critical exponents defined in an analogous manner to the usual ones except that the critical point is approached via isobaric transformations.

In order to boil water while heating, one may better keep the pressure constant, rather than fix the density and take the thermodynamic limit. In the low energy regime, the particle number, N, is conserved, and thus fixing the density is equivalent to keeping the volume constant. In statistical physics, the partition function, Z_N(T, V), is the sum of Boltzmann factors, where the energy eigenvalues are generically quantized with respect to the volume of the system. One crucial difference between volume, V, and pressure, P, is then that the former is a fundamental variable of the partition function along with temperature, T, while the latter is a conjugate variable derivable from the partition function, which is a priori a single-valued function of the two continuous fundamental variables: P(T, V).

For finite N < < ∞, the partition function is analytic, and hence, changing the temperature alone leaving the volume fixed cannot generate any singularity nor a discrete phase transition. To do so, the conventional statistical physics paradigm suggests to take the thermodynamic limit, V → ∞ with the density, N/V, fixed. Yet, an alternative scheme has been also around in the literature, prescribing to keep not the density but the pressure constant, notably by London already in 1954 [13]: For related discussions see Imry et al. [14], Chaba and Pathria [15] (theoretical), Span & Wagner [16] (experimental), and Gross et al. [17–19] (micro-canonical ensemble). The pressure, not as the fundamental but as the conjugate variable, is to be fixed while the system undergoes any discrete phase transition.

In fact, under constant pressure condition, the analytic partition function of a finite system can feature a discrete phase transition provided the volume becomes a multi-valued “function” of T and P. This happens when the partition function features a spinodal curve, which is defined to be the collection of all the points on the (T, V) plane where the isothermal compressibility becomes infinite, or

$\begin{array}{ll}\Phi (T,V):=-\frac{\partial P(T,V)}{\partial V}=0.& (1)\end{array}$

Equivalently, the second volume derivative of the free energy vanishes on the spinodal curve [e.g., [21–26]], which thus marks out the frontier of local stability. Crucially, on the spinodal curve, the temperature derivative at fixed pressure acting on an arbitrary physical quantity which is a function of T and V, diverges since the isobaric derivative contains the inverse, Φ⁻¹,

$\begin{array}{ll}\begin{array}{c}\frac{\text{d }}{\text{d}T}{|}_P=\frac{\partial }{\partial T}{|}_V+\frac{\text{d}V}{\text{d}T}{|}_P\frac{\partial }{\partial V}{|}_T\hfill \\ =\frac{\partial }{\partial T}{|}_V+\frac{\frac{\partial P(T,V)}{\partial T}}{\Phi (T,V)}\frac{\partial }{\partial V}{|}_T.\hfill \end{array}& (2)\end{array}$

The existence of the spinodal curve for a given analytic partition function of a finite system has been verified: the standard textbook systems of ideal Bose gas feature a spinodal curve, only if the number of the identical particles is sufficiently yet finitely large [20, 27–30]. Small perturbations due to extra interactions may well deform the shape of the spinodal curve but hardly eliminate its existence. In the present paper, we test the theoretical assertion that the spinodal curve exists, or emerges [30], even for a generic interacting system with finitely many particles, or real molecules, while the corresponding partition function is analytic. Under this assumption, it follows that, the isothermal compressibility (or, equivalently, the inverse of Φ) is the only source of divergence when the isobaric derivative (2) acts on a physical quantity which is an analytic function of T and V. This essentially implies that the corresponding exponent of the singularity in such an isobaric situation is determined by the power series expansion of Φ in a simple manner: the isobaric scaling behavior near a critical point, (T_c, V_c), can be classified by a natural number which fixes the isobaric critical exponents in a universal and analytic manner, as we review below.

Specifically below, we analyze the (commercially available) experimental data from “NIST Reference Fluid Thermodynamic and Transport Properties Database (REFPROP)" [31], which has been regarded as the world-best dataset for various fluids (especially H₂O and CO₂¹). Among numerous molecules, we focus on twenty simple molecules including H₂O, CO₂, O₂, etc. and report that their critical exponents seem to agree with the theoretical prediction based on analyticity: for each molecule, there is a characteristic natural number, n = 2, 3, 4, 5, 6, which determines all the critical exponents for T < T_c as ${\alpha }_P={\gamma }_P=\frac{n}{n+1}$ and ${\beta }_P={\delta }^{-1}=\frac{1}{n+1}$. For the opposite T > T_c, all the fluids feature the universal value of n = 2.

2. Review: Prediction from Analyticity

We start with the canonical partition function, Z_N(T, V), of a generic N-particle system. For finite N < < ∞, the partition function is an analytic function of the two continuous variables, T and V. The temperature dependence originates from the Boltzmann factor, e^−βE with β = (kBT)⁻¹, and the volume dependence arises due to the fact that the energy eigenvalues are quantized with respect to the size of the system. Naturally, all the physical quantities computable from the partition function are a priori functions of T and V too. For example, the pressure is given by

$\begin{array}{ll}P(T,V)=k\text{B}T\frac{\partial ln{Z}_N(T,V)}{\partial V}.& (3)\end{array}$

It follows infinitesimally,

$\begin{array}{ll}\text{d}P=\frac{\partial P(T,V)}{\partial T}\text{d}T+\frac{\partial P(T,V)}{\partial V}\text{d}V,& (4)\end{array}$

and hence with the pressure kept fixed, i.e., on an isobar, we get

$\begin{array}{ll}\frac{\text{d}T}{\text{d}V}{|}_P=-\left(\frac{\partial P(T,V)}{\partial V}\right){\left(\frac{\partial P(T,V)}{\partial T}\right)}^{-1}.& (5)\end{array}$

The spinodal curve (1) is defined to be the collection of all the points on the (T, V) plane where the partial volume derivative of the pressure vanishes. Physically, it corresponds to the boundary between the thermodynamically stable points, Φ > 0, and the unstable points, Φ < 0. The existence of a spinodal curve may not be always guaranteed for a given (analytic) partition function. Yet, it has been shown, first numerically [20, 27] and then through analytic method [28, 29] that, standard textbook systems of ideal Bose gas feature a spinodal curve, and hence the isobar on the (T, V) plane zigzags, if the number of particles is sufficiently large: specifically N ≥ 7, 616 for the 3D cubic box or N ≥ 35, 131 for the 2D square box. Surely, these dimensionless critical numbers of particles are mathematically rather unnatural: the natural one would have been ∞, i.e., the thermodynamic limit. But, the numbers are quantum mechanically generated characteristics of the geometric shape of the box, i.e., “cube” or “square.” Further, their existence is anyhow physically natural, as no experiment deals a genuinely infinite system. Extra perturbative interaction may well change the dimensionless critical numbers and the shapes of the spinodal curves. Yet, what matters is the appearance of a spinodal curve for large enough N². Hereafter, even for a generic interacting real system with finitely many particles, we shall assume, merely that the spinodal curve exists or emerges, and further tactically that the temperature derivative of the pressure does not vanish (at least) on the spinodal curve,

$\begin{array}{ll}\frac{\partial P(T,V)}{\partial T}{|}_{\Phi =0}\ne 0.& (6)\end{array}$

It follows from (5) that $\frac{\text{d}T}{\text{d}V}{|}_P$ vanishes on the spinodal curve, such that the isothermal line (dT = 0, dV≠0) is tangent to the isobar (dP = 0). This is also consistent with the definition of the spinodal curve (1).

Crucially, on the spinodal curve, the temperature derivative at fixed pressure acting on an arbitrary physical quantity which is a function of temperature and volume diverges, since the derivative contains the inverse, Φ⁻¹, as shown in Equation (2). This realizes a discrete phase transition for a finite system under constant pressure condition. Further, if the physical quantity on which the derivative (2) acts is an analytic function of T and V, the inverse, Φ⁻¹, is the only source of a divergence: none of the partial derivatives, neither $\frac{\partial }{\partial T}{|}_V$ nor $\frac{\partial }{\partial V}{|}_T$, can produce any singularity. This implies that the corresponding exponent of the singularity is determined by the power series expansion of Φ on the isobar and therefore, the exponents can be classified to take universal and exact values, under the assumptions of analytic partition function and isobaric transformations.

The critical point, (T_c, V_c), is defined to be a particular point on the spinodal curve on which the higher order volume derivative of the pressure vanishes:

$\begin{array}{ll}\begin{array}{c}\Phi (T_c,V_c)=-\frac{\partial P(T_c,V_c)}{\partial V}=0,\\ \frac{{\partial }^{2}P(T_c,V_c)}{\partial V^2}=-\frac{\partial \Phi (T_c,V_c)}{\partial V}=0.\end{array}& (7)\end{array}$

The first condition is the definition of the spinodal curve such that, as explained above, the isothermal line is tangent to the isobar, while the second condition implies in a similar fashion that the isothermal line is also tangent to the spinodal curve. Thus, the three curves, isothermal, spinodal and isobar, are all tangent to each other at a critical point. In general, for a certain natural number, n ≥ 2, we let

$\begin{array}{ll}\begin{array}{ccc}\frac{{\partial }^{k}P(T_c,V_c)}{\partial V^k}=0& \mathrm{\text{for }}& 1\le k\le n,\end{array}& (8)\end{array}$

and

$\begin{array}{ll}\frac{{\partial }^n\Phi (T_c,V_c)}{\partial V^n}=-\frac{{\partial }^{n+1}P(T_c,V_c)}{\partial V^{n+1}}\ne 0.& (9)\end{array}$

Now, at least locally, around any spinodal point, we invert P(T, V) to express the temperature as a function of P and V. This is justified by the assumption of (6) via the implicit function theorem which also ensures that the function T(P, V) is analytic as well. We further fix the pressure to be critical, i.e., P_c = P(T_c, V_c). This gives a local expression of the temperature around the critical point as a function of the volume under the fixed critical pressure,

$\begin{array}{ll}\begin{array}{cc}T=f_c(V),& T_c=f_c(V_c),\end{array}& (10)\end{array}$

satisfying

$\begin{array}{ll}P(f_c(V),V)=P_c.& (11)\end{array}$

The analyticity of the partition function, combined with the implicit function theorem, guarantees that f_c(V) is an analytic function of V. By taking the volume-derivative of the relation (11), which should vanish, it is straightforward to show that the power series expansion of the function, f_c(V), around the critical volume starts at the (n + 1)th order,

$\begin{array}{ll}T-T_c=\frac{1}{(n+1)!}\left[\frac{{\text{d}}^{n+1}{f}_c(V_c)}{\text{d}{V}^{n+1}}\right]{(V-V_c)}^{n+1}+\mathrm{\text{higher orders}},& (12)\end{array}$

where, from our assumptions of (6), (8), and (9), the leading order coefficient is non-trivial,

$\begin{array}{ll}\frac{{\text{d}}^{n+1}{f}_c(V_c)}{\text{d}{V}^{n+1}}=\left[\frac{{\partial }^n\Phi (T_c,V_c)}{\partial V^n}\right]{\left[\frac{\partial P(T_c,V_c)}{\partial T}\right]}^{-1}\ne 0.& (13)\end{array}$

Clearly, Equation (12) leads to the following isobaric critical exponent:

$\begin{array}{ll}\begin{array}{cc}V/V_c-1~{|T/T_c-1|}^{{\beta }_P},& {\beta }_P=\frac{1}{n+1},\end{array}& (14)\end{array}$

which implies

$\begin{array}{ll}\frac{\text{d}V}{\text{d}T}{|}_{P=P_c}~{|T/T_c-1|}^{-\frac{n}{n+1}}.& (15)\end{array}$

Similarly, with (9) and (12), the expansion of Φ on the critical isobar starts from the nth order in V−V_c such that, in agreement with (15),

$\begin{array}{ll}\Phi (f(V),V)~{(V/V_c-1)}^n~{|T/T_c-1|}^{\frac{n}{n+1}}.& (16)\end{array}$

Thus, from (2), any physical quantity given by the temperature derivative along the critical isobar diverges with the universal exponent, $\frac{n}{n+1}$. This includes the critical exponent of the specific heat under constant pressure:

$\begin{array}{l}\begin{array}{ccc}{C}_P~{|T/T_c-1|}^{-{\alpha }_P},& {\alpha }_P=\frac{n}{n+1}.& (17)\end{array}\end{array}$

On the other hand, since the specific heat at constant volume or constant density is finite for any finite system, the corresponding critical exponent should be trivial (if sufficiently zoomed in) [16]³:

$\begin{array}{l}\begin{array}{cccc}{C}_V=\mathrm{\text{finite }}& \mathrm{\text{i.e.}},& {\alpha }_V=0.& (18)\end{array}\end{array}$

Further, the inverse of Φ gives the critical exponent of the isothermal compressibility on the critical isobar:

$\begin{array}{l}\begin{array}{ccc}{\kappa }_T:=-\frac{\partial lnV(P_c,T)}{\partial P}~{|T/T_c-1|}^{-{\gamma }_P},& {\gamma }_P=\frac{n}{n+1}.& (19)\end{array}\end{array}$

This isothermal compressibility exponent, ${\gamma }_P=\frac{n}{n+1}$, is strictly less than one (19), and hence differs from the typical estimations of its constant volume analog γ_V [4–11]. Nevertheless, the equality of α_P = γ_P was reported experimentally for some liquid-crystalline materials, such as nCB and $\bar{m}\text{S}5$ [32].

Finally, from (9), we obtain the isothermal critical exponent at the critical temperature, T = T_c:

$\begin{array}{ll}\begin{array}{cc}P/P_c-1~{|V/V_c-1|}^{\delta },& \delta =n+1.\end{array}& (20)\end{array}$

Note that, unlike α_P, β_P and γ_P, our last exponent, δ, is not defined under the isobaric condition, and yet there should be no ambiguity in the definition (20).

In summary, under the combined hypotheses of analyticity and existence of a spinodal curve, we obtain the following universal values for the (isobaric) critical exponents [20]:

$\begin{array}{ll}\begin{array}{cc}{\alpha }_P={\gamma }_P=\frac{n}{n+1},& {\beta }_P={\delta }^{-1}=\frac{1}{n+1}.\end{array}& (21)\end{array}$

They are all determined by a single natural number, n, which should be greater than one, n = 2, 3, 4, 5, ⋯, and corresponds to the characteristic of a critical point. Henceforth, we name the natural number, critical index. As can be easily checked, the isobaric critical exponents satisfy the “scaling laws” [4–11]:

$\begin{array}{l}\begin{array}{cc}\mathrm{\text{Rushbrooke}}:& {\alpha }_P+2{\beta }_P+{\gamma }_P=2,\\ \mathrm{\text{Widom}}:& {\gamma }_P={\beta }_P(\delta -1).& (22)\end{array}\end{array}$

It is worthwhile to note that, at a generic spinodal point with lower pressure, P<P_c, the case of n = 1 for (8) and (9) corresponds to the first-order phase transition between the liquid and the gas, which features the superheating and supercooling phenomena. Even at these points, our analysis implies the universal exponent:

$\begin{array}{ll}{\alpha }_P={\beta }_P={\gamma }_P={\delta }^{-1}=\frac{1}{2}.& (23)\end{array}$

We stress that our analysis and the resulting critical exponents (21) as well as (23) should hold universally not only for liquid-vapor transitions but also for any isobaric phase transition where its own spinodal curve can be defined. Such examples include the Nematic to Smetic-A phase transition of the liquid-crystalline experiment (without δ) [32] and a relativistic ideal Bose gas theory (see Figure 5 of [20]). Both examples feature the same critical index, “n = 2” and agree with (23) when P < P_c.

3. Van Der Waals Fluid: Toy Model of n = 2

The first historical example of mean field theory was the van der Waals fluid, which we show here to provide a simple example of the above general consideration. The van der Waals equation reads

$\begin{array}{ll}\left(P_r+\frac{3}{v_r^2}\right)\left(v_r-\frac{1}{3}\right)=\frac{8}{3}{T}_r,& (24)\end{array}$

where P_r = P/P_c, v_r = v/v_c, and T_r = T/T_c are reduced pressure, volume per particle and temperature, such that the critical point is given by P_r = v_r = T_r = 1. Around the critical point, the critical isobar, P_r ≡ 1, is given by

$\begin{array}{ll}{T}_r=\frac{1}{8}(3v_r-1)\left(1+\frac{3}{v_r^2}\right)=1+\frac{3{(v_r-1)}^3}{8v_r^2}.& (25)\end{array}$

Therefore, near the critical point, $T_r-1~{(v_r-1)}^3$, we get n = 2, and hence ${\beta }_P=\frac{1}{3}$. On the other hand, taking the v_r-derivative of (24), with fixed temperature, gives

$\begin{array}{ll}\left(\frac{\partial P_r(v_r,T_r)}{\partial v_r}-\frac{6}{v_r^3}\right)\left(v_r-\frac{1}{3}\right)+P_r+\frac{3}{v_r^2}=0,& (26)\end{array}$

and hence, combined with (25), we get ${\gamma }_P=\frac{2}{3}$, as follows

$\begin{array}{ll}\begin{array}{c}{\kappa }_T=-\frac{\partial ln{v}_r(P_r,T_r)}{\partial P_r}{|}_{P_r=1}=\frac{v_r^2(v_r-\frac{1}{3})}{(v_r+2){(v_r-1)}^2}\\ ~\frac{1}{{(v_r-1)}^2}~{(T_r-1)}^{-\frac{2}{3}}.\end{array}& (27)\end{array}$

From (24), the critical isothermal relation between pressure and temperature reads, with T_r ≡ 1,

$\begin{array}{ll}{P}_r-1=-\frac{{(v_r-1)}^3}{(v_r-\frac{1}{3})v_r^2}~{(v_r-1)}^3.& (28)\end{array}$

Thus, we note δ = 3. Lastly, in order to obtain α_P, one should first remember that the specific heat at constant volume of any fluid whose equation of state is linear in temperature is equal to that of the classical ideal gas. Thus, it is finite and constant, so the corresponding critical exponent should be trivial, i.e., α_V = 0 (18). On the other hand, the specific heat at constant pressure, C_P, is given by the thermodynamic identity (see e.g., [10]),

$\begin{array}{ll}{C}_P-C_V=-T_r{\left(\frac{\partial P_r}{\partial T_r}{|}_{v_r}\right)}^2\frac{\partial v_r}{\partial P_r}{|}_{T_r},& (29)\end{array}$

which is proportional to κ_T. Therefore, α_P, the critical exponent of the specific heat under constant pressure takes the same value as ${\gamma }_P=\frac{2}{3}$. Put all together, the critical point of van der Waals fluid belongs to the class of n = 2. Note that the textbook approach to obtain the critical exponents of van der Waals fluid is based on the concept of Maxwell's equal area construction [9, 10, 33], which, we assert, cannot be applied nor defined at the critical point itself. More precisely, the conventional derivation of the exponent, “$\beta =\frac{1}{2}$" is based on the behavior of the binodal curve of liquid-vapor coexistence near the critical point, while our isobaric critical exponent, ${\beta }_P=\frac{1}{3}$, is obtained here by simply fixing the pressure to be critical, P_r ≡ 1, in the very van der Waals equation (25).

4. Result and Discussion

We recall the theoretical prediction of the critical exponents (21), ${\alpha }_P={\gamma }_P=\frac{n}{n+1}$, ${\beta }_P={\delta }^{-1}=\frac{1}{n+1}$, which can be easily derived not by assuming the thermodynamic limit nor any singularity of the free energy, but by postulating the analyticity of the partition function and the existence of a spinodal curve. Our main result in this work is to verify that such a simple approach appears relevant to the description of the critical phenomena of the liquid-gas phase transition of ‘real fluids’: the analyticity-based theoretical prediction is consistent with the high-precision experimental data, REFPROP of NIST [31].

The NIST data are sufficiently precise even close to the critical points: the uncertainties for most fluids do not exceed 2%. For examples, the uncertainty in pressure of nitrogen in the critical region is estimated to be 0.02%. The uncertainty in pressure of water in the critical region is 0.1%. The uncertainty in heat capacities of argon is within 0.3% for the vapor and 2% for the liquid. The uncertainties in vapor pressure of helium are less than 0.02% and for the heat capacities are about 2%. The detailed information on the data can be obtained by downloading the (commercially available) REFPROP database program through its website [31]. We tend to believe that such high-precision data were not available while the conventional paradigm of taking the thermodynamic limit and achieving non-analyticity was proposed and established during the mid twentieth century. In particular, the isothermal relation between pressure and volume appears to have been rather poorly measured [e.g., [5]].

By analyzing the data for randomly chosen twenty simple molecules, we are able to classify the critical phenomena of those real fluids in terms of the critical index, n₋ = 2, 3, 4, 5, 6, for T < T_c. For the opposite range of temperature, T > T_c, the natural number seems to be universal as n₊ = 2, which agrees with van der Waals fluid (24), a relativistic ideal Bose gas (n₋ = n₊ = 2 [20]), and an experimental result on SF₆ [34]. Table 1 and Figures 1–4 show the main results.

TABLE 1

Table 1. Classification of the 20 major fluids plus the relativistic ideal Bose gas and the van der Waals fluid, according to the critical indices: n₋ for T < T_c and n₊ for T > T_c.

FIGURE 1

Figure 1. Critical exponent, ${\alpha }_P=\frac{n}{n+1}$ : n₋ = 2, 3, 4, 5, 6 for T < T_c and n₊ = 2 for T > T_c.

FIGURE 2

Figure 2. Critical exponent, ${\beta }_P=\frac{1}{n+1}$ : n₋ = 2, 3, 4, 5, 6 for T < T_c and n₊ = 2 for T > T_c.

FIGURE 3

Figure 3. Critical exponent, ${\gamma }_P=\frac{n}{n+1}$ : n₋ = 2, 3, 4, 5, 6 for T < T_c and n₊ = 2 for T > T_c.

FIGURE 4

Figure 4. Critical exponent, δ = n+1 : n₋ = 2, 3, 4, 5, 6 for V < V_c and n₊ = 2 for V > V_c.

Table 1 contains our classification of the twenty molecular fluids, plus the relativistic ideal Bose gas and the van der Waals fluid, according to the critical indices: n₋ for T < T_c and n₊ for T > T_c. Note the universal value, n₊ = 2. We have determined the critical indices from the experimental data of NIST (REFPROP) [31] especially over the reduced temperature range, $10^{-5}<|T/T_c-1|<10^{-2}$.

Figures 1–3 show the log-log plot of the scaled physical quantities (C_P, |V/V_c−1|, and κ_T) vs. the reduced temperature (|T/T_c−1|) for the twenty molecules, where data are logarithmically-binned from REFPROP dataset for our analysis. On the other hand, Figure 4 shows the log-log plot of the scaled reduced pressure (|P/P_c−1|) vs. the reduced volume (|V/V_c−1|) for the same molecules, where data are also logarithmically-binned. All the figures represent the best fit plots of experimental data obtained from the twenty molecules, which confirm the scaling behavior characterized by the critical indices, n₋ = 2, 3, 4, 5, 6 for T < T_c and n₊ = 2 for T > T_c.

Generically, except C₅H₁₀, we have obtained two distinct critical indices, n₋ and n₊, depending on T < T_c and T > T_c. This rather surprising result leads us to speculate —still assuming the analyticiy— that there may exist more than one critical points quite close to each other, such as (T_c−, V_c−:n₋) and (T_c+, V_c+:n₊) (Lee et al., in preparation).

We may regard our result as a tantalizing hint for the finiteness of the Avogadro number in real experimental laboratories, such that analyticity appears to persist at least for the close up range like $10^{-5}<|T/T_c-1|<10^{-3}$. On the other hand, slightly further away from the critical point, e.g., $10^{-3}<|T/T_c-1|<10^{-1}$, the singularity or the critical index appears to get enhanced still in a universal manner for β₋, γ₋, δ₋. This might correspond to a thermodynamic limiting behavior.

As a matter of fact, NIST REFPROP assumes certain form of extrapolation formulas for data fitting, such that it appears able to generate data for extremely short range, e.g., $10^{-9}<|T/T_c-1|<10^{-5}$. Yet, such data turn out to deviate from the analytic prediction of our interest, which should be due to the artificially designed fitting formulas. While for the range of $10^{-5}<|T/T_c-1|<10^{-3}$, we trust NIST REFPROP represents true experimental data. For this range it is a highly non-trivial verification of us that all the four critical exponents of a given real molecule are fixed by a ‘common’ natural number, i.e., the critical index.

5. Conclusion

Finite systems can undergo first or second order phase transitions under isobaric condition without resorting to the thermodynamic limit. The “analyticity” of a finite-system partition function predicts then the universal values of the isobaric critical exponents. In this work, analyzing NIST REFPROP data for 20 simple molecules, including H₂O, CO₂, O₂, etc. over the the reduced temperature range of $10^{-5}<|T/T_c-1|<10^{-3}$, we have tested and verified the prediction: for each molecule, there appears to exist a characteristic natural number, n = 2, 3, 4, 5, 6, which determines all the critical exponents for T < T_c as ${\alpha }_P={\gamma }_P=\frac{n}{n+1}$ and ${\beta }_P={\delta }^{-1}=\frac{1}{n+1}$. For the opposite T > T_c, all the fluids seem to indicate the universal value of n = 2.

In the sense of Kuhn [35], our result might be viewed as “anomaly" within the standard paradigm of taking the thermodynamic limit for critical phenomena, e.g., the (non-analytic) 3D Ising model or XY model [36]. Nevertheless, our analysis represents truthfully the NIST experimental data for the specific temperature window, $10^{-5}<|T/T_c-1|<10^{-3}$. We call for further theoretical as well as experimental investigations.

Author Contributions

J-HP proposed the research. WC contributed the NIST data handling and analysis. D-HK led the interpretation of the data. J-HP and D-HK wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We wish to thank Xavier Bekaert for numerous enlightening discussions and Eric W. Lemmon at NIST for helpful correspondences. This work was supported by Sogang University (Grant Nos. 201610033.01 and 201710066.01) and the National Research Foundation of Korea (Grant Nos. 2015K1A3A1A21000302 and 2016R1D1A1B01015196).

Footnotes

1. ^Private communications with E. W. Lemmon at NIST.

2. ^In order to identify any spinodal curve, or more generically to observe nontrivial quantum effects, one should rather keep the analyticity of the partition function and aim to perform not a somewhat hasty approximate integral but an exact discrete sum over countable quantum states, c.f. Nekrasov NA. Seiberg-Witten prepotential from instanton counting. Adv. Theor. Math. Phys. 7 831 (2003).

3. ^On the other hand, if appropriately zoomed out, it might be possible to observe pseudo-thermodynamic-limiting behavior, with finite peak.

References

1. Anderson PW. More is different. Science (1972) 177:393–6.

PubMed Abstract | Google Scholar

2. Yang CN, Lee TD. Statistical theory of equations of state and phase transitions. I. Theory of condensation. Phys Rev. (1952) 87:404.

Google Scholar

3. Yang CN, Lee TD. Statistical theory of equations of state and phase transitions. II. Lattice gas and Ising model. Phys Rev. (1952) 87:410.

Google Scholar

4. Fisher ME. The theory of equilibrium critical phenomena. Rep Prog Phys. (1967) 30:615–713.

Google Scholar

5. Heller P. Experimental investigations of critical phenomena. Rep Prog Phys. (1967) 30:731–826.

Google Scholar

6. Kadanoff LP, Götze W, Hamblen D, Hecht R, Lewis EAS, Palciauskas VV, et al. Static phenomena near critical points: theory and experiment. Rev Mod Phys. (1967) 39:395.

Google Scholar

7. Stanley HE. Introduction to Phase Transitions and Critical Phenomena. New York, NY: Oxford University Press (1971).

Google Scholar

8. Wilson KG, Kogut J. The renormalization group and the ϵ expansion. Phys Rep. (1974) 12:75–200.

Google Scholar

9. Huang K. Statistical Mechanics. 2nd ed. New York, NY: Wiley (1987).

10. Goldenfeld N. Lectures on Phase Transitions and the Renormalization Group. Vol. 85. Reading, MA: Addison-Wesley (1992).

Google Scholar

11. Kadanoff LP. Statistical Physics: Statics, Dynamics and Renormalization. Singapore: World Scientific (2000).

Google Scholar

12. Kadanoff LP. More is the same; phase transitions and mean field theories. J Stat Phys. (2009) 137:777. doi: 10.1007/s10955-009-9814-1

CrossRef Full Text | Google Scholar

13. London F. Superfluids, Vol. 2, Macroscopic Theory of Superfluid Helium. New York, NY: Dover Publications (1954).

14. Imry Y, Bergman DJ, Gunther L. Bose-Einstein condensation in two-dimensional systems. In: Chen TC, Nordheim LW, editors. Low Temperature Physics-LT 13. Vol. 1. New York, NY: Springer; Dover Publications, Inc. (1974). p. 80.

Google Scholar

15. Chaba AN, Pathria RK. Bose-Einstein condensation in a two-dimensional system at constant pressure. Phys Rev B (1975) 12:3697.

Google Scholar

16. Span R, Wagner W. A new equation of state for carbon dioxide covering the fluid region from the triple-point temperature to 1100 K at pressures up to 800 MPa. J Phys Chem Ref Data (1996) 25:1509–96.

Google Scholar

17. Gross DHE, Votyakov E. Phase transitions in ‘small’ systems. Eur Phys J B (2000)15:115. doi: 10.1063/1.555991

CrossRef Full Text | Google Scholar

18. Gross DHE. Geometric foundation of thermo-statistics, phase transitions, Second Law of Thermodynamics, but without thermodynamic limit. Phys Chem Chem Phys. (2002) 4:863–75. doi: 10.1039/B108775D

CrossRef Full Text | Google Scholar

19. Gross DHE, Kenney JF. The microcanonical thermodynamics of finite systems: the Microscopic origin of condensation and phase separations: and the conditions for heat flow from lower to higher temperatures. J Chem Phys. (2005) 122:224111. doi: 10.1063/1.1901658

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Park JH, Kim SW. Existence of a critical point in the phase diagram of ideal relativistic neutral Bose gas. New J Phys.(2011) 13:033003. doi: 10.1088/1367-2630/13/3/033003

CrossRef Full Text | Google Scholar

21. Gibbs JW. The Scientific Papers of J. Willard Gibbs New York, NY: Dover Publications, Inc. (1961).

Google Scholar

22. Skripov VP. Metastable Liquids J. Wiley (1973).

23. Chandler D. Introduction to Modern Statistical Mechanics. New York, NY: Oxford University Press (1987).

24. Chimowitz EH. Introduction to Critical Phenomena in Fluids. New York, NY: Oxford University Press (2005).

25. O'Connell JP, Haile JM. Thermodynamics: Fundamentals for Applications. Cambridge: Cambridge University Press (2005).

Google Scholar

26. Scott Shell M. Thermodynamics and Statistical Mechanics: An Integrated Approach. Cambridge: Cambridge University Press (2015).

27. Park JH, Kim SW. Thermodynamic instability and first-order phase transition in an ideal Bose gas. Phys Rev A (2010) 81:063636. doi: 10.1103/PhysRevA.81.063636

CrossRef Full Text | Google Scholar

28. Jeon I, Kim SW, Park JH. Isobar of an ideal Bose gas within the grand canonical ensemble. Phys Rev A (2011) 84:023636. doi: 10.1103/PhysRevA.84.023636

CrossRef Full Text | Google Scholar

29. Cho W, Kim SW, Park JH. Two-dimensional Bose-Einstein condensate under pressure. New J Phys. (2015) 17:013038. doi: 10.1088/1367-2630/17/1/013038

CrossRef Full Text | Google Scholar

30. Park JH. How many is different? Answer from ideal Bose gas. J Phys Conf Ser. (2014) 490:012018. doi: 10.1088/1742-6596/490/1/012018

CrossRef Full Text | Google Scholar

31. Lemmon EW, Huber ML, McLinden MO. NIST Reference Fluid Thermodynamic and Transport Properties Database (REFPROP): Version 9.1, National Institute of Standards and Technology (NIST). Gaithersburg, MD: Standard Reference Data Program (2013). Available online at: https://www.nist.gov/srd/refprop

32. Zywocinski A, Wieczorek SA. Critical exponents for thermal expansion and isothermal compressibility near the Nematic to Smectic-A phase transition. J Phys Chem B (1997) 101:6970–6. doi: 10.1021/jp971234g

CrossRef Full Text | Google Scholar

33. Levelt Sengers JMH. From Van der Waals' equation to the scaling laws. Physica (1974) 73:73–106. doi: 10.1016/0031-8914(74)90227-4

CrossRef Full Text | Google Scholar

34. Wagner W, Kurzeja N, Pieperbeck B. The thermal behaviour of pure fluid substances in the critical region - experiences from recent pρT measurements on SF₆ with a multi-cell apparatus. Fluid Phase Equilibria (1992) 79:151–74.

Google Scholar

35. Kuhn TS. The Structure of Scientific Revolutions. (Chicago, IL: University of Chicago Press) (1962) (2nd ed.: 1970; 3rd ed.: 1996; 4th ed.: 2012).

36. Garland CW, Nounesis G. Critical behavior at nematic–smectic–A phase transitions. Phys Rev E (1994) 49:2964.

PubMed Abstract | Google Scholar

37. Guggenheim EA. The principle of corresponding states. J Chem Phys. (1945) 13:253.

Google Scholar

Appendix: Comment on Guggenheim's Fitting

Around the critical point, with Φ_c ≡ Φ(T_c, V_c) = 0 and from (9), we may expand

$\begin{array}{ll}\begin{array}{c}\Phi (T,V)=\left(\frac{\partial {\Phi }_c}{\partial T}+\frac{{\partial }^2{\Phi }_c}{\partial T\partial V}\Delta V+\frac{1}{2}\frac{{\partial }^2{\Phi }_c}{\partial T^2}\Delta T\right)\Delta T\\ +\frac{1}{n!}\frac{{\partial }^n{\Phi }_c}{\partial V^n}{\Delta V}^n+\cdots ,\end{array}& (30)\end{array}$

where ΔT = T−T_c, ΔV = V−V_c. Then, near to the critical point, in comparison with the critical isobar (14), the spinodal curve can be approximated by

$\begin{array}{ll}0<T/T_c-1~{(V/V_c-1)}^n.& (31)\end{array}$

Our result implies then that, the supercooling (T < T_c) and the superheating (T > T_c) spinodal curves read separately in terms of the density, ρ_± = N/V_±,

$\begin{array}{ll}{\rho }_{\pm }/{\rho }_c-1~|T/T_c-1{|}^{1/n_{\pm }}.& (32)\end{array}$

From the empirical result reported in this work, n₊<n₋ (except C₅H₁₀), we expect close to the critical point,

$\begin{array}{ll}|T/T_c-1{|}^{1/n_{+}}<<|T/T_c-1{|}^{1/n_{-}}.& (33)\end{array}$

Thus, the difference between the densities of the supercooled gas and the superheated liquid can be approximated by

$\begin{array}{ll}{\rho }_{+}-{\rho }_{-}~|T/T_c-1{|}^{1/n_{-}}.& (34)\end{array}$

This is comparable with the renowned Guggenheim's fitting [37] with n = 3.