Van der Waals equation

Gas equation of state which accounts for non-ideal gas behavior

The van der Waals equation is a mathematical formula that describes the behavior of real gases. It is an equation of state that relates the pressure, volume, number of molecules, and temperature in a fluid. The equation modifies the ideal gas law,  p V = N k T {\displaystyle pV=NkT} in which  k {\displaystyle k} is the Boltzmann constant, in two ways: first, it considers particles to have a finite diameter (whereas an ideal gas consists of point particles); second, its particles interact with each other (unlike an ideal gas, whose particles move as though alone in the volume).

The equation is named after Dutch physicist Johannes Diderik van der Waals, who first derived it in 1873 as part of his doctoral thesis. Van der Waals based the equation on the idea that fluids are composed of discrete particles, which few scientists believed existed. However, the equation accurately predicted the behavior of a fluid around its critical point, which had been discovered a few years earlier. Its qualitative and quantitative agreement with experiments ultimately cemented its acceptance in the scientific community. These accomplishments won van der Waals the 1910 Nobel Prize in Physics.[1] Today the equation is recognized as an important model of phase change processes.[2]

## Description

[edit]

Figure A: Van der Waals equation on a wall in Leiden

One explicit way to write the van der Waals equation is:[3][4]

p = R T v − b − a v 2 , {\displaystyle \qquad p={\frac {RT}{v-b}}-{\frac {a}{v^{2}}},} |  | 1a

|
|

where  p {\displaystyle p} is pressure,  T {\displaystyle T} is temperature, and  v = V / n = N A V / N {\displaystyle v=V/n=N_{\text{A}}V/N} is molar volume, the ratio of volume,  V {\displaystyle V} , to quantity of matter,  n {\displaystyle n} ( N A {\displaystyle N_{\text{A}}} is the Avogadro constant and  N {\displaystyle N} the number of molecules). Also  a {\displaystyle a} and  b {\displaystyle b} are experimentally determinable, substance-specific constants, and  R = k N A {\displaystyle R=kN_{\text{A}}} is the universal gas constant. This form is useful for plotting isotherms (constant temperature curves).

Van der Waals wrote it in an equivalent, explicit in temperature, form in his Thesis[5][6] (although he did not, and could not, write absolute temperature in its modern form in 1873)

( p + a v 2 ) ( v − b ) = R T . {\displaystyle \qquad \left(p+{\frac {a}{v^{2}}}\right)(v-b)=RT.} |  | 1b

|
|

This form is useful for plotting isobars (constant pressure curves). With  v = V / n {\displaystyle v=V/n} , it becomes the form that appears in Figure A.[7]

When van der Waals created his equation, few scientists believed that fluids were composed of rapidly moving particles. Moreover, those who thought so did not know the atomic/molecular structure. The simplest conception of a particle, and the easiest to model mathematically was a hard sphere of volume  V 0 {\displaystyle V_{0}} ; this is what van der Waals used, and he found the total excluded volume was ⁠ B = 4 N V 0 {\displaystyle B=4NV_{0}} ⁠, namely 4 times the volume of all the particles.[8][9] The constant ⁠ b = B N A / N {\displaystyle b=BN_{\text{A}}/N} ⁠, has the dimension of molar volume, [v]. The constant  a {\displaystyle a} expresses the strength of the hypothesized inter-particle attraction. Van der Waals only had Newton's law of gravitation, in which two particles are attracted in proportion to the product of their masses, as a model. Thus he argued that, in his case, the attractive pressure was proportional to the density squared.[10] The proportionality constant, ⁠ a {\displaystyle a} ⁠, when written in the form used above, has the dimension [pv2] (pressure times molar volume squared).

Figure B: The Sutherland potential (orange) represents two hard spheres that attract according to an inverse power law, and the Lennard-Jones potential (black) represents the induced-dipole–induced-dipole interaction of two non-polar molecules. Both are simple realistic molecular models.[11]

The force magnitude between two spherically symmetric molecules is written as  F = − d φ / d r {\displaystyle F=-d\varphi /dr} , where  φ ( r ) {\displaystyle \varphi (r)} is the pair potential function, and the force direction is along the line connecting the two mass centers. The specific functional relation is most simply characterized by a single length, ⁠ σ {\displaystyle \sigma } ⁠, and a minimum energy,  − ε {\displaystyle -\varepsilon } (with ⁠ ε ≥ 0 {\displaystyle \varepsilon \geq 0} ⁠). Two of the many such functions that have been suggested are shown in Fig. B.[11]

A modern theory based on statistical mechanics produces the same result for  b = 4 N A [ ( 4 π / 3 ) ( σ / 2 ) 3 ] {\displaystyle b=4N_{\text{A}}[(4\pi /3)(\sigma /2)^{3}]} obtained by van der Waals and his contemporaries. It also produces a constant value for  a / N A ε b {\displaystyle a/N_{\text{A}}\varepsilon b} when  ε / k T {\displaystyle \varepsilon /kT} is small enough.[12][13]

Once the constants  a {\displaystyle a} and  b {\displaystyle b} are known for a given substance, the van der Waals equation can be used to predict attributes like the boiling point at any given pressure, and the critical point. These predictions are accurate for only a few substances. Most simple fluids are only a valuable approximation.

### Relationship to the ideal gas law

[edit]

The ideal gas law follows from the van der Waals equation whenever the molar volume  v {\displaystyle v} is sufficiently large (when ⁠ v ≫ b {\displaystyle v\gg b} ⁠, so ⁠ v − b ≈ v {\displaystyle v-b\approx v} ⁠), or equivalently whenever the molar density, ⁠ ρ = 1 / v {\displaystyle \rho =1/v} ⁠, is sufficiently small (when ⁠ v ≫ ( a / p ) 1 / 2 {\displaystyle v\gg (a/p)^{1/2}} ⁠, so ⁠ p + a / v 2 ≈ p {\displaystyle p+a/v^{2}\approx p} ⁠).[14]

When  v {\displaystyle v} is large enough that both inequalities are satisfied, these two approximations reduce the van der Waals equation to ⁠ p = R T / v {\displaystyle p=RT/v} ⁠, or  p v = R T {\displaystyle pv=RT} . With  R = N A k {\displaystyle R=N_{\text{A}}k} , where  k {\displaystyle k} is the Boltzmann constant, and using the definition  v = V / n {\displaystyle v=V/n} given after Eq (1a), this becomes ⁠ p V = N k T {\displaystyle pV=NkT} ⁠; either of these forms expresses the ideal gas law.[14] This is unsurprising since the van der Waals equation was constructed from the ideal gas equation to obtain an equation valid beyond the low density limit of ideal gas behavior.

What is truly remarkable is the extent to which van der Waals succeeded. Indeed, Epstein in his classic thermodynamics textbook began his discussion of the van der Waals equation by writing, "Despite its simplicity, it comprehends both the gaseous and the liquid state and brings out, in a most remarkable way, all the phenomena pertaining to the continuity of these two states".[14] Also, in Volume 5 of his Lectures on Theoretical Physics, Sommerfeld, in addition to noting that "Boltzmann[15] described van der Waals as the Newton of real gases",[6] also wrote "It is very remarkable that the theory due to van der Waals is in a position to predict, at least qualitatively, the unstable [referring to superheated liquid, and subcooled vapor, now called metastable] states" that are associated with the phase change process.[16]

## History

[edit]

In 1857 Rudolf Clausius published The Nature of the Motion which We Call Heat. In it he derived the relation  p = ( N / V ) m c 2 ¯ / 3 {\displaystyle p=(N/V)m{\overline {c^{2}}}/3} for the pressure  p {\displaystyle p} in a gas, composed of particles in motion, with number density ⁠ N / V {\displaystyle N/V} ⁠, mass ⁠ m {\displaystyle m} ⁠, and mean square speed ⁠ c 2 ¯ {\displaystyle {\overline {c^{2}}}} ⁠. He then noted that using the classical laws of Boyle and Charles, one could write  m c 2 ¯ / 3 = k T {\displaystyle m{\overline {c^{2}}}/3=kT} with a constant of proportionality ⁠ k {\displaystyle k} ⁠. Hence temperature was proportional to the average kinetic energy of the particles.[17] This article inspired further work based on the twin ideas that substances are composed of indivisible particles, and that heat is a consequence of the particle motion; movement that evolves according to Newton's laws. The work, known as the kinetic theory of gases, was done principally by Clausius, James Clerk Maxwell, and Ludwig Boltzmann. At about the same time, Josiah Willard Gibbs advanced the work by converting it into statistical mechanics.[18][19]

This environment influenced Johannes Diderik van der Waals. After initially pursuing a teaching credential, he was accepted for doctoral studies at the University of Leiden under Pieter Rijke. This led, in 1873, to a dissertation that provided a simple, particle-based equation that described the gas–liquid change of state, the origin of a critical temperature, and the concept of corresponding states.[20][21] The equation is based on two premises: first, that fluids are composed of particles with non-zero volumes, and second, that at a large enough distance each particle exerts an attractive force on all other particles in its vicinity. Boltzmann called these forces van der Waals cohesive forces.[22]

In 1869 Irish professor of chemistry Thomas Andrews at Queen's University Belfast, in a paper entitled On the Continuity of the Gaseous and Liquid States of Matter,[23] displayed an experimentally obtained set of isotherms of carbonic acid, H2CO3, that showed at low temperatures a jump in density at a certain pressure, while at higher temperatures there was no abrupt change (the figure can be seen here). Andrews called the isotherm at which the jump just disappears the critical point. Given the similarity of the titles of this paper and van der Waals' subsequent thesis, one might think that van der Waals set out to develop a theoretical explanation of Andrews' experiments; however, this is not what happened. Van der Waals began work by trying to determine a molecular attraction that appeared in Laplace's theory of capillarity, and only after establishing his equation he tested it using Andrews' results.[24][25]

By 1877 sprays of both liquid oxygen and liquid nitrogen had been produced, and a new field of research, low temperature physics, had been opened. The van der Waals equation played a part in all this, especially for the liquefaction of hydrogen and helium which was finally achieved in 1908.[26] From measurements of  p 1 , T 1 {\displaystyle p_{1},T_{1}} and  p 2 , T 2 {\displaystyle p_{2},T_{2}} in two states with the same density, the van der Waals equation produces the values[27]

b = v − R ( T 2 − T 1 ) p 2 − p 1 and a = v 2 p 2 T 1 − p 1 T 2 T 2 − T 1 . {\displaystyle b=v-{\frac {R(T_{2}-T_{1})}{p_{2}-p_{1}}}\qquad {\text{and}}\qquad a=v^{2}{\frac {p_{2}T_{1}-p_{1}T_{2}}{T_{2}-T_{1}}}.}

Thus from two such measurements of pressure and temperature, one could determine  a {\displaystyle a} and  b {\displaystyle b} , and from these values calculate the expected critical pressure, temperature, and molar volume. Goodstein summarized this contribution of the van der Waals equation as follows:[28]

> All this labor required considerable faith in the belief that gas–liquid systems were all basically the same, even if no one had ever seen the liquid phase. This faith arose out of the repeated success of the van der Waals theory, which is essentially a universal equation of state, independent of the details of any particular substance once it has been properly scaled. [...] As a result, not only was it possible to believe that hydrogen could be liquefied, but it was even possible to predict the necessary temperature and pressure.

Van der Waals was awarded the Nobel Prize in 1910, in recognition of the contribution of his formulation of this "equation of state for gases and liquids".

## Use

[edit]

The van der Waals equation has been, and remains, useful because:[29]

* Its coefficient of thermal expansion has a simple analytic expression
* It explains the existence of the critical point, and establishes the theorem of corresponding states
* Its internal energy and entropy have simple analytic expressions
* Its specific heat at constant volume  c v {\displaystyle c_{v}} is a function of  T {\displaystyle T} only
* Its specific heat at constant pressure,  c p {\displaystyle c_{p}} has a simple relationship with  c v {\displaystyle c_{v}}
* Its Joule–Thomson coefficient and associated inversion curve, which are instrumental in the commercial liquefaction of gases, have simple analytic expressions
* Together with the Maxwell construction it explains the existence of the liquid–vapor phase transition, including the observed metastable states

In addition[30]

* Its enthalpy and free energies all have simple analytic expressions
* Its isothermal compressibility has a simple analytic expression
* Its saturation curve has a simple analytic parametric solution
* It is an intermediate mathematical model that is useful as a pedagogical tool when teaching physics, chemistry, and engineering

and

* It plays an important role in the modern theory of phase transitions[31]
* It is the completely accurate equation of state for substances whose intermolecular potential matches the Sutherland potential[32]

## Critical point and corresponding states

[edit]

Figure 1: Four isotherms of the van der Waals equation along with the spinodal curve (black dash-dot curve) and the coexistence (saturation) curve (red dash-dot curve), plotted using reduced (dimensionless) variables. The critical point lies at the inflection point on the orange isotherm.

Figure 1 shows four isotherms of the van der Waals equation (abbreviated as vdW) on a  p , v {\displaystyle p,v} (pressure, molar volume) plane. The essential character of these curves is that they come in three forms:

1. At some critical temperature  T c {\displaystyle T_{\text{c}}} (orange isotherm), the slope is negative everywhere except at a single inflection point: the critical point ⁠ ( p c , v c ) {\displaystyle (p_{\text{c}},v_{\text{c}})} ⁠, where both the slope and curvature are zero, ⁠ ∂ p ∂ v | T = ∂ 2 p ∂ v 2 | T = 0 {\displaystyle \left.{\frac {\partial p}{\partial v}}\right\vert _{T}=\left.{\frac {\partial ^{2}p}{\partial v^{2}}}\right\vert _{T}=0} ⁠.
2. At higher temperatures (red isotherm), the isotherm's slope is negative everywhere. (This corresponds to values of  p , T {\displaystyle p,T} for which the vdW equation has one real root for ⁠ v {\displaystyle v} ⁠).
3. At lower temperatures (green and blue isotherms), all isotherms have two points with zero slope. (This corresponds to values of ⁠ p {\displaystyle p} ⁠,  T {\displaystyle T} for which the vdW equation has three real roots for ⁠ v {\displaystyle v} ⁠).

The critical point can be analytically determined by equating the two partial derivatives of the vdW equation, created by differentiating Eq (1a), to zero. This produces the critical values  v c = 3 b {\displaystyle v_{\text{c}}=3b} and ⁠ T c = 8 a / ( 27 R b ) {\displaystyle T_{\text{c}}=8a/(27Rb)} ⁠. Finally, using these values in Eq (1a) gives ⁠ p c = a / 27 b 2 {\displaystyle p_{\text{c}}=a/27b^{2}} ⁠.[33] These results can also be obtained algebraically by noting that at the critical point the three roots are equal. Hence, Eqs (1) can be written as either  v 3 − ( b + R T c / p c ) v 2 + ( a / p c ) v − a b / p c = 0 {\displaystyle v^{3}-(b+RT_{\text{c}}/p_{\text{c}})v^{2}+(a/p_{\text{c}})v-ab/p_{\text{c}}=0} , or  ( v − v c ) 3 = 0 {\displaystyle (v-v_{\text{c}})^{3}=0} ; two forms with the same coefficients.[34][35]

### Course of the isotherms

[edit]

Above the critical temperature  T c {\displaystyle T_{\text{c}}} , van der Waals isotherms satisfy the stability criterion that  ∂ p / ∂ v | T < 0 {\displaystyle \partial p/\partial v|_{T}<0} . Below the critical temperature, each isotherm contains an interval where this condition is violated. This unstable region is the genesis of the phase change; there is a range  v m i n ≤ v ≤ v m a x {\displaystyle v_{\rm {min}}\leq v\leq v_{\rm {max}}} , for which no observable states exist. The states for  v < v m i n {\displaystyle v<v_{\rm {min}}} are liquid, and those for  v > v m a x {\displaystyle v>v_{\rm {max}}} are vapor; the denser liquid separates and lies below the vapor due to gravity. The transition points, states with zero slope, are called spinodal points.[36] Their locus is the spinodal curve, a boundary that separates the regions of the plane for which liquid, vapor, and gas exist from a region where no observable homogeneous states exist. This spinodal curve is obtained here from the vdW equation by differentiation (or equivalently from  κ T = ∞ {\displaystyle \kappa _{T}=\infty } ) as  T s p = 2 a ( v − b ) 2 R v 3 p s p = a ( v − 2 b ) v 3 {\displaystyle T_{\rm {sp}}=2a{\frac {(v-b)^{2}}{Rv^{3}}}\qquad p_{\rm {sp}}={\frac {a(v-2b)}{v^{3}}}}

A projection of the spinodal curve is plotted in Figure 1 as the black dash-dot curve. It passes through the critical point, which is also a spinodal point.

### Principle of corresponding states

[edit]

Using the critical values to define reduced (dimensionless) variables ⁠ p r = p / p c {\displaystyle p_{r}=p/p_{\text{c}}} ⁠, ⁠ T r = T / T c {\displaystyle T_{r}=T/T_{\text{c}}} ⁠, and  v r = v / v c {\displaystyle v_{r}=v/v_{\text{c}}} renders the vdW equation in the dimensionless form (used to construct Fig. 1):

p r = 8 T r 3 v r − 1 − 3 v r 2 {\displaystyle p_{r}={\frac {8T_{r}}{3v_{r}-1}}-{\frac {3}{v_{r}^{2}}}}

This dimensionless form is a similarity relation; it indicates that all vdW fluids at the same  T r {\displaystyle T_{r}} will plot on the same curve. It expresses the law of corresponding states which Boltzmann described as follows:[37]

> All the constants characterizing the gas have dropped out of this equation. If one bases measurements on the van der Waals units [Boltzmann's name for the reduced quantities here], then he obtains the same equation of state for all gases. [...] Only the values of the critical volume, pressure, and temperature depend on the nature of the particular substance; the numbers that express the actual volume, pressure, and temperature as multiples of the critical values satisfy the same equation for all substances. In other words, the same equation relates the reduced volume, reduced pressure, and reduced temperature for all substances.

> Obviously such a broad general relation is unlikely to be correct; nevertheless, the fact that one can obtain from it an essentially correct description of actual phenomena is very remarkable.

This "law" is just a special case of dimensional analysis in which an equation containing 6 dimensional quantities, ⁠ p , v , T , a , b , R {\displaystyle p,v,T,a,b,R} ⁠, and 3 independent dimensions, [p], [v], [T], must be expressible in terms of 6 − 3 = 3 dimensionless groups.[38] Here  v ∗ = b {\displaystyle v^{*}=b} is a characteristic molar volume,  p ∗ = a / b 2 {\displaystyle p^{*}=a/b^{2}} a characteristic pressure, and  T ∗ = a / ( R b ) {\displaystyle T^{*}=a/(Rb)} a characteristic temperature, and the 3 dimensionless groups are ⁠ p / p ∗ , v / v ∗ , T / T ∗ {\displaystyle p/p^{*},v/v^{*},T/T^{*}} ⁠. According to dimensional analysis the equation must then have the form ⁠ p / p ∗ = Φ ( v / v ∗ , T / T ∗ ) {\displaystyle p/p^{*}=\Phi (v/v^{*},T/T^{*})} ⁠, a general similarity relation. In his discussion of the vdW equation, Sommerfeld also mentioned this point.[39] The reduced properties defined previously are ⁠ p r = 27 ( p / p ∗ ) {\displaystyle p_{r}=27(p/p^{*})} ⁠, ⁠ v r = ( 1 / 3 ) ( v / v ∗ ) {\displaystyle v_{r}=(1/3)(v/v^{*})} ⁠, and ⁠ T r = ( 27 / 8 ) ( T / T ∗ ) {\displaystyle T_{r}=(27/8)(T/T^{*})} ⁠. Recent research has suggested that there is a family of equations of state that depend on an additional dimensionless group, and this provides a more exact correlation of properties.[40] Nevertheless, as Boltzmann observed, the van der Waals equation provides an essentially correct description.

The vdW equation produces the critical compressibility factor ⁠ Z c = p c v c / ( R T c ) = 3 / 8 = 0.375 {\displaystyle Z_{\text{c}}=p_{\text{c}}v_{\text{c}}/(RT_{\text{c}})=3/8=0.375} ⁠, while for most real fluids ⁠ 0.23 < Z c < 0.31 {\displaystyle 0.23<Z_{\text{c}}<0.31} ⁠.[41] Thus most real fluids do not satisfy this condition, and consequently their behavior is only described qualitatively by the vdW equation. However, the vdW equation of state is a member of a family of state equations based on the Pitzer (acentric) factor, ⁠ ω {\displaystyle \omega } ⁠, and the liquid metals (mercury and cesium) are well approximated by it.[42][32]

## Thermodynamic properties

[edit]

The properties molar internal energy,  u {\displaystyle u} , and entropy,  s {\displaystyle s} , are defined by the first and second laws of thermodynamics. From these laws, they, and all other thermodynamic properties of a simple compressible substance, can be specified, up to a constant of integration, by two measurable functions. These are a mechanical equation of state,  p = p ( v , T ) {\displaystyle p=p(v,T)} , and a constant volume specific heat,  c v ( v , T ) {\displaystyle c_{v}(v,T)} .[43][44]

When  u ( v , T ) {\displaystyle u(v,T)} represents a continuous surface, it must be a continuous function with continuous partial derivatives, and its second mixed partial derivatives must be equal,  ∂ v ∂ T u = ∂ T ∂ v u {\displaystyle \partial _{v}\partial _{T}u=\partial _{T}\partial _{v}u} . Then with  c v = ∂ T u {\displaystyle c_{v}=\partial _{T}u} this condition can be written simply as  ∂ v c ( v , T ) = ∂ T [ T 2 ∂ T ( p / T ) ] {\displaystyle \partial _{v}c(v,T)=\partial _{T}[T^{2}\partial _{T}(p/T)]} . Differentiating  p / T {\displaystyle p/T} for the vdW equation gives  T 2 ∂ T ( p / T ) ] = a / v 2 {\displaystyle T^{2}\partial _{T}(p/T)]=a/v^{2}} , so  ∂ v c v = 0 {\displaystyle \partial _{v}c_{v}=0} . Consequently  c v = c v ( T ) {\displaystyle c_{v}=c_{v}(T)} for a vdW fluid exactly as it is for an ideal gas.[45] To keep things simple, it is regarded as a constant in the following,  c v = c R {\displaystyle c_{v}=cR} , with  c {\displaystyle c} a number.

### Internal energy, and entropy

[edit]

The energetic equation of state gives the internal energy, and the entropic equation of state gives the entropy as[46][44]

u − C u = ∫ c v ( v , T ) d T + ∫ [ T 2 ∂ ( p / T ) ∂ T ] d v {\displaystyle u-C_{u}=\int \,c_{v}(v,T)\,dT+\int \,\left[T^{2}{\frac {\partial (p/T)}{\partial T}}\right]\,dv} s − C s = ∫ c v ( T ) d T T + ∫ ∂ p ∂ T d v {\displaystyle s-C_{s}=\int \,c_{v}(T){\frac {dT}{T}}+\int \,{\frac {\partial p}{\partial T}}\,dv}

where  C u , C s {\displaystyle C_{u},C_{s}} are arbitrary constants of integration.

Both integrals for  u {\displaystyle u} can be easily evaluated and the result is, [47][48]

u − C u = c R T − a / v {\displaystyle \qquad u-C_{u}=cRT-a/v} |  | 2

|
|

Likewise both integrals for  s {\displaystyle s} can be evaluated with the result,[41][49]

s − C s = R ln ⁡ [ T c ( v − b ) ] {\displaystyle \qquad s-C_{s}=R\ln[T^{c}(v-b)]} |  | 3

|
|

### Free energies, and enthalpy

[edit]

The Helmholtz free energy is  f = u − T s {\displaystyle f=u-Ts} . Subtracting  T {\displaystyle T} times Eq (3) from Eq (2) gives  f {\displaystyle f} as[50]

f = C u + c T − a / v − T { C s + R ln ⁡ [ T c ( v − b ) ] } {\displaystyle \qquad f=C_{u}+cT-a/v-T\\{C_{s}+R\ln[T^{c}(v-b)]\\}} |  | 4a

|
|

The enthalpy is  h = u + p v {\displaystyle h=u+pv} , and the product  p v {\displaystyle pv} is, using Eq (1a),  p v = R T v / ( v − b ) − a / v {\displaystyle pv=RTv/(v-b)-a/v} . Adding Eq (2) gives  h {\displaystyle h} as [47][48] h − C u = R T [ c + v / ( v − b ) ] − 2 a / v {\displaystyle h-C_{u}=RT[c+v/(v-b)]-2a/v}

The Gibbs free energy is  g = h − T s {\displaystyle g=h-Ts} so subtracting  T {\displaystyle T} times Eq (3) from  h {\displaystyle h} produces  g {\displaystyle g} as[51]

g − C u = ( c + v v − b − C s − ln ⁡ [ T c ( v − b ) ] ) R T − 2 a v {\displaystyle \qquad g-C_{u}=\left(c+{v \over v-b}-C_{\text{s}}-\ln[T^{c}(v-b)]\right)RT-{2a \over v}} |  | 4b

|
|

All these results can be rendered in reduced form by using the characteristic energy  R T c {\displaystyle RT_{\text{c}}} .

### Derivatives: α, κT and cp

[edit]

Any derivative of any thermodynamic property can be expressed in terms of any three of them.[52]. A standard set is composed of  α , κ T , c v {\displaystyle \alpha ,\kappa _{T},c_{v}} . For a vdW fluid  c v ( T ) {\displaystyle c_{v}(T)} is a known function, and the other two are obtained from the first partial derivatives of the vdW equation as,  ∂ p ∂ T ) v = R v − b = α κ T and ∂ p ∂ v ) T = − R T ( v − b ) 2 + 2 a v 3 = − 1 v κ T {\displaystyle \left.{\frac {\partial p}{\partial T}}\right)_{v}={\frac {R}{v-b}}={\frac {\alpha }{\kappa _{T}}}\quad {\mbox{and}}\quad \left.{\frac {\partial p}{\partial v}}\right)_{T}=-{\frac {RT}{(v-b)^{2}}}+{\frac {2a}{v^{3}}}=-{\frac {1}{v\kappa _{T}}}}

Here  κ T = − v − 1 ∂ p v {\displaystyle \kappa _{T}=-v^{-1}\partial _{p}v} , is the isothermal compressibility, and  α = v − 1 ∂ T v p {\displaystyle \alpha =v^{-1}\partial _{T}v_{p}} , is the coefficient of thermal expansion.[53][54] Therefore,[55][56]

κ T = v 2 ( v − b ) 2 R T v 3 − 2 a ( v − b ) 2 and α = R v 2 ( v − b ) R T v 3 − 2 a ( v − b ) 2 . {\displaystyle \kappa _{T}={\frac {v^{2}(v-b)^{2}}{RTv^{3}-2a(v-b)^{2}}}\quad {\text{and}}\quad \alpha ={\frac {Rv^{2}(v-b)}{RTv^{3}-2a(v-b)^{2}}}.} |  | 5

|
|

In the limit  v -> ∞ {\displaystyle v\rightarrow \infty } ,  α = 1 / T {\displaystyle \alpha =1/T} and  κ T = v / ( R T ) {\displaystyle \kappa _{T}=v/(RT)} .[56][57] Since the vdW equation in this limit becomes  p = R T / v {\displaystyle p=RT/v} , finally  κ T = 1 / p {\displaystyle \kappa _{T}=1/p} . Both of these are the ideal gas values.

The specific heat at constant pressure,  c p {\displaystyle c_{p}} is defined as the partial derivative  c p = ∂ T h | p {\displaystyle c_{p}=\partial _{T}h|_{p}} . It is related to  c v {\displaystyle c_{v}} by the Mayer equation,  c p − c v = − T ( ∂ T p ) 2 / ∂ v p = T v α 2 / κ T {\displaystyle c_{p}-c_{v}=-T(\partial _{T}p)^{2}/\partial _{v}p=Tv\alpha ^{2}/\kappa _{T}} .[58][59] Then the two partials of the vdW equation can be used to express  c p {\displaystyle c_{p}} as,[51]

c p ( v , T ) − c v ( T ) = R 2 T v 3 R T v 3 − 2 a ( v − b ) 2 ≥ R {\displaystyle c_{p}(v,T)-c_{v}(T)={\frac {R^{2}Tv^{3}}{RTv^{3}-2a(v-b)^{2}}}\geq R} |  | 6

|
|

Here in the limit  v -> ∞ {\displaystyle v\rightarrow \infty } ,  c p − c v = R {\displaystyle c_{p}-c_{v}=R} , which is also the ideal gas result;[51] however the limit  v -> b {\displaystyle v\rightarrow b} gives the same result, which does not agree with experiments on liquids.[57]

Finally  c p , α {\displaystyle c_{p},\alpha } , and  κ T {\displaystyle \kappa _{T}} are all infinite on the curve  T = 2 a ( v − b ) 2 / ( R v 3 ) = T c ( 3 v r − 1 ) 2 / ( 4 v r 3 ) {\displaystyle T=2a(v-b)^{2}/(Rv^{3})=T_{\text{c}}(3v_{r}-1)^{2}/(4v_{r}^{3})} .[51] This is the spinodal curve defined by  κ T − 1 = 0 {\displaystyle \kappa _{T}^{-1}=0} ,[60] that was discussed in the subsection The course of the isotherms.

## Saturation

[edit]

Although the gap in  v {\displaystyle v} delimited by the two spinodal points on an isotherm (e.g.  T r = 7 / 8 {\displaystyle T_{\text{r}}=7/8} in Fig. 1) is the origin of the phase change, the change occurs at some intermediate value of pressure. This can be understood by considering that the saturated liquid and vapor states can coexist in equilibrium, with the same pressure and temperature.[61] However, the minimum and maximum spinodal points are not at the same pressure. Therefore, at a temperature  T s {\displaystyle T_{\text{s}}} , the phase change is characterized by the pressure  p s {\displaystyle p_{\text{s}}} , which lies within the range of  p {\displaystyle p} set by the spinodal points ( p min < p s < p max {\displaystyle p_{\text{min}}<p_{\text{s}}<p_{\text{max}}} ), and by the molar volume of liquid  v f {\displaystyle v_{\text{f}}} and vapor  v g {\displaystyle v_{\text{g}}} , which lie outside the range of  v {\displaystyle v} set by the spinodal points ( v f < v min {\displaystyle v_{\text{f}}<v_{\text{min}}} and  v g > v max {\displaystyle v_{\text{g}}>v_{\text{max}}} ).

Applying Eq (1a) to the saturated liquid and saturated vapor states gives:

p s = R T s v f − b − a v f 2 and p s = R T s v g − b − a v g 2 . {\displaystyle \qquad p_{\text{s}}={\frac {RT_{\text{s}}}{v_{\text{f}}-b}}-{\frac {a}{{v_{\text{f}}}^{2}}}\quad {\text{and}}\quad p_{\text{s}}={\frac {RT_{\text{s}}}{v_{\text{g}}-b}}-{\frac {a}{{v_{\text{g}}}^{2}}}.} |  | 7

|
|

Equations (7) contain four variables  p s , T s , v f , v g {\displaystyle p_{\text{s}},T_{\text{s}},v_{\text{f}},v_{\text{g}}} ), so a third equation is required to uniquely specify three of these variables in terms of the fourth. In this case of a single substance, the equation is provided by the condition of equal Gibbs free energy,[61]

g g = g f {\displaystyle g_{\text{g}}=g_{\text{f}}}

Using Eq (4b) applied to each state in this equation produces

R T s [ v g v g − b − v f v f − b − ln ⁡ ( v g − b v f − b ) ] − 2 a ( 1 v g − 1 v f ) = 0 {\displaystyle \qquad RT_{\text{s}}\left[{\frac {v_{\text{g}}}{v_{\text{g}}-b}}-{\frac {v_{\text{f}}}{v_{\text{f}}-b}}-\ln \left({\frac {v_{\text{g}}-b}{v_{\text{f}}-b}}\right)\right]-2a\left({\frac {1}{v_{\text{g}}}}-{\frac {1}{v_{\text{f}}}}\right)=0} |  | 8

|
|

This is a third equation that, along with Eqs. 7 can be solved numerically. This has been done given a value for either  T s {\displaystyle T_{\text{s}}} or  p s {\displaystyle p_{\text{s}}} , and tabular results presented;[62][63] however, the equations also admit an analytic parametric solution obtained by Lekner.[64] Details of this solution may be found in the Maxwell construction, and the dimensionless results are:

T rs ( y ) = ( 27 8 ) 2 f ( y ) [ cosh ⁡ y + f ( y ) ] g ( y ) 2 p rs = 27 f ( y ) 2 [ 1 − f ( y ) 2 ] g ( y ) 2 {\displaystyle T_{\text{rs}}(y)=\left({\frac {27}{8}}\right){\frac {2f(y)[\cosh y+f(y)]}{g(y)^{2}}}\quad p_{\text{rs}}=27{\frac {f(y)^{2}[1-f(y)^{2}]}{g(y)^{2}}}} v rf = ( 1 3 ) 1 + f ( y ) e y f ( y ) e y v rg = ( 1 3 ) 1 + f ( y ) e − y f ( y ) e − y {\displaystyle v_{\text{rf}}=\left({\frac {1}{3}}\right){\frac {1+f(y)e^{y}}{f(y)e^{y}}}\qquad \qquad \qquad \quad v_{\text{rg}}=\left({\frac {1}{3}}\right){\frac {1+f(y)e^{-y}}{f(y)e^{-y}}}} where  f ( y ) = y cosh ⁡ y − sinh ⁡ y sinh ⁡ y cosh ⁡ y − y g ( y ) = 1 + 2 f ( y ) cosh ⁡ y + f ( y ) 2 {\displaystyle f(y)={\frac {y\cosh y-\sinh y}{\sinh y\cosh y-y}}\qquad \qquad g(y)=1+2f(y)\cosh y+f(y)^{2}}

The parameter  0 ≤ y < ∞ {\displaystyle 0\leq y<\infty } is given physically by  y = ( s g − s f ) / ( 2 R ) {\displaystyle y=(s_{\text{g}}-s_{\text{f}})/(2R)} . This solution also produces values of all other property discontinuities across the saturation curve.[65] These functions define the coexistence curve (or saturation curve), which is the locus of the saturated liquid and saturated vapor states of the vdW fluid. Projections of this saturation curve are plotted in Figures 1 and 2.

Figure 2: The black dash-dot curve is the spinodal curve (stability limit) and the blue dash-dot curve is the coexistence (saturation) curve, plotted in the  p r , T r {\displaystyle p_{\text{r}},T_{\text{r}}} plane.

Referring back to Figure 1, the isotherms for  T r < 1 {\displaystyle T_{\text{r}}<1} are discontinuous. For example, the  T r = 7 / 8 {\displaystyle T_{\text{r}}=7/8} (green) isotherm consists of two separate segments. The solid green lines are composed of stable states. They terminate at dots representing the saturated liquid and vapor states forming the phase change. The dashed green lines represent metastable states (superheated liquid and subcooled vapor). They are created in the phase transition, have a finite lifetime, and then devolve into their lower energy stable alternative.

At every point in the region between the two curves in Figure 2, there are two states: one stable and one metastable. The coexistence of these states can be seen in Figure 1—for discontinuous isotherms, there are values of  p r {\displaystyle p_{\text{r}}} which correspond to two points on the isotherm: one on a solid line (the stable state) and one on a dashed region (the metastable state).

In his treatise of 1898, in which he described the van der Waals equation in great detail, Boltzmann discussed these metastable states in a section titled "Undercooling, Delayed evaporation".[66] (Today, these states are now denoted "subcooled vapor" and "superheated liquid".) Moreover, it has now become clear that these metastable states occur regularly in the phase transition process. In particular, processes that involve very high heat fluxes create large numbers of these states, and transition to their stable alternative with a corresponding release of energy that can be dangerous. Consequently, there is a pressing need to study their thermal properties.[67]

In the same section, Boltzmann also addressed and explained the negative pressures which some liquid metastable states exhibit (for example, the blue isotherm  T r = 4 / 5 {\displaystyle T_{\text{r}}=4/5} in Fig. 1). He concluded that such liquid states of tensile stresses were real, as did Tien and Lienhard many years later who wrote "The van der Waals equation predicts that at low temperatures liquids sustain enormous tension [...] In recent years measurements have been made that reveal this to be entirely correct."[68]

Even though the phase change produces a mathematical discontinuity in the homogeneous fluid properties (for example  v {\displaystyle v} ), there is no physical discontinuity.[69] As the liquid begins to vaporize, the fluid becomes a heterogeneous mixture of liquid and vapor whose molar volume varies continuously from  v f {\displaystyle v_{\text{f}}} to  v g {\displaystyle v_{\text{g}}} according to the equation of state  v = v f + x ( v g − v f ) {\textstyle v=v_{\text{f}}+x(v_{\text{g}}-v_{\text{f}})} where  x = N g / ( N f + N g ) {\textstyle x=N_{\text{g}}/(N_{\text{f}}+N_{\text{g}})} and  0 ≤ x ≤ 1 {\displaystyle 0\leq x\leq 1} is the mole fraction of the vapor. This equation is called the lever rule and applies to other properties as well.[16][69] The states it represents form a horizontal line bridging the discontinuous region of an isotherm (not shown in Fig. 1 because it is a different equation from the vdW equation).

### Extended corresponding states

[edit]

Figure 3: The family of saturation curves, showing the vdW curve as a member (blue curve). The blue dots are calculated from Lekner's solution. The orange dots are calculated from data in the ASME Steam Tables Compact Edition, 2006.

The idea of corresponding states originated when van der Waals cast his equation in the dimensionless form,  p r = p ( v r , T r ) {\displaystyle p_{\text{r}}=p(v_{\text{r}},T_{\text{r}})} . However, as Boltzmann noted, such a simple representation could not correctly describe all substances. Indeed, the saturation analysis of this form produces  p rs = p s ( T r ) {\displaystyle p_{\text{rs}}=p_{\text{s}}(T_{\text{r}})} ; namely, that all substances have the same dimensionless coexistence curve, which is not true.[70] To avoid this paradox, an extended principle of corresponding states has been suggested in which  p r = p ( v r , T r , ϕ ) {\displaystyle p_{\text{r}}=p(v_{\text{r}},T_{\text{r}},\phi )} where  ϕ {\displaystyle \phi } is a substance-dependent dimensionless parameter related to the only physical feature associated with an individual substance: its critical point.

Figure 4: A plot of the correlation including data from various substances.

One candidate for  ϕ {\displaystyle \phi } is the critical compressibility factor  Z c = p c v c / ( R T c ) {\displaystyle Z_{\text{c}}=p_{\text{c}}v_{\text{c}}/(RT_{\text{c}})} ; however, because  v c {\displaystyle v_{\text{c}}} is difficult to measure accurately, the acentric factor developed by Kenneth Pitzer,[71] ω = − log 10 ⁡ [ p r ( T r = 0.7 ) ] − 1 {\displaystyle \omega =-\log _{10}[p_{\text{r}}(T_{\text{r}}=0.7)]-1} , is more useful. The saturation pressure in this situation is represented by a one-parameter family of curves:  p rs = p s ( T r , ω ) {\displaystyle p_{\text{rs}}=p_{\text{s}}(T_{\text{r}},\omega )} . Several investigators have produced correlations of saturation data for several substances; Dong and Lienhard give[32] ln ⁡ p rs = 5.37270 ( 1 − 1 / T r ) + ω ( 7.49408 − 11.181777 T r 3 + 3.68769 T r 6 + 17.92998 ln ⁡ T r ) {\displaystyle {\begin{aligned}\ln p_{\text{rs}}=\ &5.37270(1-1/T_{\text{r}})+\omega (7.49408-\\\&11.181777\ {T_{\text{r}}}^{3}+3.68769\ {T_{\text{r}}}^{6}+17.92998\,\ln T_{\text{r}})\end{aligned}}} which has an RMS error of  ± 0.42 {\displaystyle \pm 0.42} over the range  1 ≤ T r ≤ 0.3 {\displaystyle 1\leq T_{\text{r}}\leq 0.3} .

Figure 3 is a plot of  p rs {\displaystyle p_{\text{rs}}} vs.  T r {\displaystyle T_{\text{r}}} for various values of the Pitzer factor  ω {\displaystyle \omega } as given by this equation. The vertical axis is logarithmic to show the behavior at pressures closer to zero, where differences among the various substances (indicated by varying values of  ω {\displaystyle \omega } ) are more pronounced.

Figure 4 is another plot of the same equation showing  T r {\displaystyle T_{\text{r}}} as a function of  ω {\displaystyle \omega } for various values of  p rs {\displaystyle p_{\text{rs}}} . It includes data from 51 substances, including the vdW fluid, over the range  − 0.4 < ω < 0.9 {\displaystyle -0.4<\omega <0.9} . This plot shows that the vdW fluid ( ω = − 0.302 {\displaystyle \omega =-0.302} ) is a member of the class of real fluids; indeed, the vdW fluid can quantitatively approximate the behavior of the liquid metals cesium ( ω = − 0.267 {\displaystyle \omega =-0.267} ) and mercury ( ω = − 0.21 {\displaystyle \omega =-0.21} ), which share similar values of  ω {\displaystyle \omega } . However, in general it can describe the behavior of fluids of various  ω {\displaystyle \omega } only qualitatively.

## Joule–Thomson coefficient

[edit]

The Joule–Thomson coefficient,  μ JT = ∂ p T | h {\displaystyle \mu _{\text{JT}}=\partial _{p}T|_{h}} , is of practical importance because the two end states of a throttling process ( h 2 = h 1 {\displaystyle h_{2}=h_{1}} ) lie on a constant enthalpy curve. Although ideal gases, for which  h = h ( T ) {\displaystyle h=h(T)} , do not change temperature in such a process, real gases do, and it is important in applications to know whether they heat up or cool down.[72]

This coefficient can be found in terms of the previously derived  α {\displaystyle \alpha } and  c p {\displaystyle c_{p}} as[73] μ JT = v ( α T − 1 ) c p . {\displaystyle \mu _{\text{JT}}={\frac {v(\alpha T-1)}{c_{p}}}.}

When  μ JT {\displaystyle \mu _{\text{JT}}} is positive, the gas temperature decreases as it passes through a throttling process, and when it is negative, the temperature increases. Therefore, the condition  μ JT = 0 {\displaystyle \mu _{\text{JT}}=0} defines a curve that separates the region of the  T , p {\displaystyle T,p} plane where  μ JT > 0 {\displaystyle \mu _{\text{JT}}>0} from the region where  μ JT < 0 {\displaystyle \mu _{\text{JT}}<0} . This curve is called the inversion curve, and its equation is  α T − 1 = 0 {\displaystyle \alpha T-1=0} . Evaluating this using the expression for  α {\displaystyle \alpha } derived in Eq. 5 produces,[74]

Figure 5: Curves of constant enthalpy in this plane have negative slope above this (green) inversion curve, positive slope below it, and zero slope on it; they are S-shaped. A gas entering a throttle at a state corresponding to a point on this curve to the right of its maximum will cool if the final state is below the curve. Right: a close-up of the region between zero and the critical point (1,1), showing the overlap between the inversion curve (green) and saturation curve (dashed purple).

2 a ( v − b ) 2 − R T v 2 b = 0 {\displaystyle 2a(v-b)^{2}-RTv^{2}b=0}

Note that for  v ≫ b {\displaystyle v\gg b} there will be cooling for  2 a > R T b {\displaystyle 2a>RTb} (or, in terms of the critical temperature,  T < ( 27 / 4 ) T c {\displaystyle T<(27/4)\ T_{\text{c}}} ). As Sommerfeld noted, "This is the case with air and with most other gases. Air can be cooled at will by repeated expansion and can finally be liquified."[74]

Solving for  b / v > 0 {\displaystyle b/v>0} , and using this to eliminate  v {\displaystyle v} from Eq (1a) gives the inversion curve as

p p ∗ = − 1 + 4 ( T 2 T ∗ ) 1 2 − 3 ( T 2 T ∗ ) {\displaystyle {\frac {p}{p^{*}}}=-1+4\left({\frac {T}{2T^{*}}}\right)^{\frac {1}{2}}-3\left({\frac {T}{2T^{*}}}\right)}

where, for simplicity,  a , b , R {\displaystyle a,b,R} have been replaced by  p ∗ , T ∗ {\displaystyle p^{*},T^{*}} .

A plot of the curve, in reduced variables, is shown in green in Figure 5\. Sommerfeld also displays this plot,[75] together with a curve drawn using experimental data from H2. The two curves agree qualitatively, but not quantitatively. Figure 5 shows an overlap between the saturation curve and the inversion curve plotted in the same region. This crossover means a van der Waals gas can be liquified by passing it through a throttling process under the proper conditions; real gases are liquified in this way.

## Compressibility factor

[edit]

Figure 6: Isotherms plotted against the compressibility factor  Z {\displaystyle Z} . The horizontal axis is  v ∗ ρ = ρ r / 3 {\displaystyle v^{*}\rho =\rho _{r}/3} . The isotherms, spinodal curve, and coexistence curve here are the same as in Fig. 1. Additionally plotted are the isotherm  T r = 27 / 8 {\displaystyle T_{r}=27/8} , which has zero slope at the origin, and the isotherm  T r = ∞ {\displaystyle T_{r}=\infty } . Figure 7: Generalized compressibility chart for a van der Waals gas.

Real gases are characterized by their difference from ideal gases by writing  p v = Z R T {\displaystyle pv=ZRT} , where  Z {\displaystyle Z} , called the compressibility factor. It is expressed either as  Z ( p , T ) {\displaystyle Z(p,T)} or  Z ( ρ , T ) {\displaystyle Z(\rho ,T)} , because in either case (pressure,  p {\displaystyle p} , or density,  ρ {\displaystyle \rho } ) the limit as  p {\displaystyle p} or  ρ {\displaystyle \rho } approaches zero is 1, and  Z {\displaystyle Z} takes the ideal gas value. In the second case  Z ( ρ , T ) = p ( ρ , T ) / ρ R T {\displaystyle Z(\rho ,T)=p(\rho ,T)/\rho RT} ,[76] so for a van der Waals fluid from Eq (‘’’1’’’) the compressibility factor is

Z = 1 1 − b ρ − a ρ R T {\displaystyle \qquad Z={1 \over 1-b\rho }-{a\rho \over RT}} |  | 9

|
|

or in terms of reduced variables  Z = 3 3 − ρ r − 9 ρ r 8 T r {\displaystyle Z={\frac {3}{3-\rho _{r}}}-{\frac {9\rho _{r}}{8T_{r}}}} where  0 ≤ ρ r = 1 / v r ≤ 3 {\displaystyle 0\leq \rho _{r}=1/v_{r}\leq 3} . At the critical point,  T r = ρ r = 1 {\displaystyle T_{r}=\rho _{r}=1} and  Z = Z c = 3 / 2 − 9 / 8 = 3 / 8 {\displaystyle Z=Z_{\text{c}}=3/2-9/8=3/8} .

In the limit  ρ -> 0 {\displaystyle \rho \rightarrow 0} ,  Z = 1 {\displaystyle Z=1} ; the fluid behaves like an ideal gas, as mentioned before. The derivative  ∂ ρ Z | T = b ( ( 1 − b ρ ) − 2 − a b R T ) {\displaystyle \partial _{\rho }Z|_{T}=b\left((1-b\rho )^{-2}-{a \over bRT}\right)} is never negative when  a b R T = T ∗ / T ≤ 1 {\displaystyle {a \over bRT}=T^{*}/T\leq 1} ; that is, when  T / T ∗ ≥ 1 {\displaystyle T/T^{*}\geq 1} (which corresponds to  T r ≥ 27 / 8 {\displaystyle T_{r}\geq 27/8} ). Alternatively, the initial slope is negative when  T / T ∗ < 1 {\displaystyle T/T^{*}<1} , is zero at  b ρ = 1 − ( T / T ∗ ) 1 / 2 {\displaystyle b\rho =1-(T/T^{*})^{1/2}} , and is positive for larger  b ρ ≤ 1 {\displaystyle b\rho \leq 1} (see Fig. 6). In this case, the value of  Z {\displaystyle Z} passes through  1 {\displaystyle 1} when  b ρ B = 1 − T B / T ∗ {\displaystyle b\rho _{B}=1-T_{B}/T^{*}} . Here  T B = ( 27 T c / 8 ) ( 1 − b ρ B ) {\displaystyle T_{B}=(27T_{\text{c}}/8)(1-b\rho _{B})} is called the Boyle temperature. It ranges between  0 ≤ T B ≤ 27 T c / 8 {\displaystyle 0\leq T_{B}\leq 27T_{\text{c}}/8} , and denotes a point in  T , ρ {\displaystyle T,\rho } space where the equation of state reduces to the ideal gas law. However, the fluid does not behave like an ideal gas there, because neither its derivatives  ( α , κ T ) {\displaystyle (\alpha ,\kappa _{T})} nor  c p {\displaystyle c_{p}} reduce to their ideal gas values, other than where  b ρ B ≪ 1 , T B ∼ 27 T c / 8 {\displaystyle b\rho _{B}\ll 1,\,T_{B}\sim 27T_{\text{c}}/8} the actual ideal gas region.[77]

Figure 6 plots various isotherms of  Z ( ρ , T r ) {\displaystyle Z(\rho ,T_{r})} vs  ρ r {\displaystyle \rho _{r}} . Also shown are the spinodal and coexistence curves described previously. The subcritical isotherm consists of stable, metastable, and unstable segments (identified in the same way as in Fig. 1). Also included are the zero initial slope isotherm and the one corresponding to infinite temperature.

Figure 7 shows a generalized compressibility chart for a vdW gas. Like all other vdW properties, this is not quantitatively correct for most gases, but it has the correct qualitative features.[78][79] Note the caustic generated by the crossing isotherms.

### Virial expansion

[edit]

Kammerlingh Onnes first suggested the virial expansion as an empirical alternative to the vdW equation. Subsequently, it was proven to result from Statistical mechanics,[80][81] in the form  Z ( ρ , T ) = 1 + ∑ k = 2 ∞ B k ( T ) ( ρ ) k − 1 {\displaystyle Z(\rho ,T)=1+\sum _{k=2}^{\infty }\,B_{k}(T)(\rho )^{k-1}} where  Z = p / ( ρ R T ) {\displaystyle Z=p/(\rho RT)} and the functions  B k ( T ) {\displaystyle B_{k}(T)} are the virial coefficients. The  k {\displaystyle k} th term represents a  k {\displaystyle k} -particle interaction.

Expanding the term  ( 1 − b ρ ) − 1 {\displaystyle (1-b\rho )^{-1}} in the definition of  Z {\displaystyle Z} , Eq (9), into an infinite series, absolutely convergent for  b ρ < 1 {\displaystyle b\rho <1} , produces  Z ( ρ , T ) = 1 + ( 1 − a b R T ) b ρ + ∑ k = 3 ∞ ( b ρ ) k − 1 . {\displaystyle Z(\rho ,T)=1+\left(1-{a \over bRT}\right)b\rho +\sum _{k=3}^{\infty }(b\rho )^{k-1}.}

The second virial coefficient is the slope of  Z ( ρ , T ) {\displaystyle Z(\rho ,T)} at  ρ = 0 {\displaystyle \rho =0} . It is positive when  T / T ∗ > 1 {\displaystyle T/T^{*}>1} and negative when  T / T ∗ < 1 {\displaystyle T/T^{*}<1} ( T r = T / T c > or < 27 / 8 {\displaystyle T_{\text{r}}=T/T_{\text{c}}>{\mbox{or}}<27/8} ), in agreement with the result found by differentiation. Its vdW value,  B 2 = b − a / R T {\displaystyle B_{2}=b-a/RT} agrees with a statistical mechanical calculation; however, the higher order coefficients are in error. This means that the vdW virial expansion, hence the vdW equation itself, is equivalent to a two term asymptotic approximation to the virial equation.[13] [82]

For molecules modeled as non-attracting hard spheres,  a = 0 {\displaystyle a=0} , and the vdW virial expansion becomes  Z ( ρ ) = ( 1 − b ρ ) − 1 = 1 + ∑ k = 2 ∞ ( b ρ ) k − 1 , {\displaystyle Z(\rho )=(1-b\rho )^{-1}=1+\sum _{k=2}^{\infty }(b\rho )^{k-1},} which illustrates the effect of the excluded volume alone. It was recognized early on that this was in error beginning with the term  ( b ρ ) 2 {\displaystyle (b\rho )^{2}} . Boltzmann calculated its correct value as  5 8 ( b ρ ) 2 {\textstyle {\frac {5}{8}}(b\rho )^{2}} , and used the result to propose an enhanced version of the vdW equation:  ( p + a v 2 ) ( v − b 3 ) = R T ( 1 + 2 b 3 v + 7 b 2 24 v 2 ) . {\displaystyle \left(p+{a \over v^{2}}\right)\left(v-{b \over 3}\right)=RT\left(1+{2b \over 3v}+{{7b^{2}} \over {24v^{2}}}\right).}

On expanding  ( v − b / 3 ) − 1 {\displaystyle (v-b/3)^{-1}} , this produced the correct coefficients through  ( b / v ) 2 {\displaystyle (b/v)^{2}} and also gave infinite pressure at  b / 3 {\displaystyle b/3} , which is approximately the close-packing distance for hard spheres.[83] This was one of the first of many equations of state proposed over the years that attempted to make quantitative improvements to the remarkably accurate explanations of real gas behavior produced by the vdW equation.[84]

## Mixtures

[edit]

In 1890 van der Waals published an article that initiated the study of fluid mixtures. It was subsequently included as Part III of a later published version of his thesis.[85] His essential idea was that in a binary mixture of vdW fluids described by the equations  p 1 = R T v − b 11 − a 11 v 2 and p 2 = R T v − b 22 − a 22 v 2 {\displaystyle p_{1}={\frac {RT}{v-b_{11}}}-{\frac {a_{11}}{v^{2}}}\quad {\text{and}}\quad p_{2}={\frac {RT}{v-b_{22}}}-{\frac {a_{22}}{v^{2}}}} the mixture is also a vdW fluid given by  p = R T v − b x − a x v 2 {\displaystyle p={\frac {RT}{v-b_{x}}}-{\frac {a_{x}}{v^{2}}}} where  a x = a 11 x 1 2 + 2 a 12 x 1 x 2 + a 22 x 2 2 and b x = b 11 x 1 2 + 2 b 12 x 1 x 2 + b 22 x 2 2 . {\displaystyle a_{x}=a_{11}x_{1}^{2}+2a_{12}x_{1}x_{2}+a_{22}x_{2}^{2}\quad {\text{and}}\quad b_{x}=b_{11}x_{1}^{2}+2b_{12}x_{1}x_{2}+b_{22}x_{2}^{2}.}

Here  x 1 = N 1 / N {\displaystyle x_{1}=N_{1}/N} and  x 2 = N 2 / N {\displaystyle x_{2}=N_{2}/N} , with  N = N 1 + N 2 {\displaystyle N=N_{1}+N_{2}} (so that  x 1 + x 2 = 1 {\displaystyle x_{1}+x_{2}=1} ), are the mole fractions of the two fluid substances. Adding the equations for the two fluids shows that  p ≠ p 1 + p 2 {\displaystyle p\neq p_{1}+p_{2}} , although for  v {\displaystyle v} sufficiently large  p ≈ p 1 + p 2 {\displaystyle p\approx p_{1}+p_{2}} with equality holding in the ideal gas limit. The quadratic forms for  a x {\displaystyle a_{x}} and  b x {\displaystyle b_{x}} are a consequence of the forces between molecules. This was first shown by Lorentz,[86] and was credited to him by van der Waals. The quantities  a 11 , a 22 {\displaystyle a_{11},\,a_{22}} and  b 11 , b 22 {\displaystyle b_{11},\,b_{22}} in these expressions characterize collisions between two molecules of the same fluid component, while  a 12 = a 21 {\displaystyle a_{12}=a_{21}} and  b 12 = b 21 {\displaystyle b_{12}=b_{21}} represent collisions between one molecule of each of the two different fluid components. This idea of van der Waals' was later called a one fluid model of mixture behavior.[87]

Assuming that  b 12 {\displaystyle b_{12}} is the arithmetic mean of  b 11 {\displaystyle b_{11}} and  b 22 {\displaystyle b_{22}} ,  b 12 = ( b 11 + b 22 ) / 2 {\displaystyle b_{12}=(b_{11}+b_{22})/2} , substituting into the quadratic form and noting that  x 1 + x 2 = 1 {\displaystyle x_{1}+x_{2}=1} produces  b = b 11 x 1 + b 22 x 2 {\displaystyle b=b_{11}x_{1}+b_{22}x_{2}}

Van der Waals wrote this relation, but did not make use of it initially.[88] However, it has been used frequently in subsequent studies, and its use is said to produce good agreement with experimental results at high pressure.[89]

### Common tangent construction

[edit]

In this article, van der Waals used the Helmholtz potential minimum principle to establish stability conditions. This principle states that in a system in diathermal contact with a heat reservoir  T = T R {\displaystyle T=T_{R}} ,  D F = 0 {\displaystyle DF=0} , and  D 2 F > 0 {\displaystyle D^{2}F>0} , namely at equilibrium, the Helmholtz potential is a minimum.[90] Since, like  g ( p , T ) {\displaystyle g(p,T)} , the molar Helmholtz function  f ( v , T ) {\displaystyle f(v,T)} is also a potential function whose differential is  d f = ∂ v f | T d v + ∂ T f | v d T = − p d v − s d T , {\displaystyle df=\partial _{v}f|_{T}\,dv+\partial _{T}f|_{v}\,dT=-p\,dv-s\,dT,} this minimum principle leads to the stability condition  ∂ 2 f / ∂ v 2 | T = − ∂ p / ∂ v | T > 0 {\displaystyle \partial ^{2}f/\partial v^{2}|_{T}=-\partial p/\partial v|_{T}>0} . This condition means that the function,  f ( v , T ) {\displaystyle f(v,T)} , is convex at all stable states of the system. Moreover, for those states the previous stability condition for the pressure is also necessarily satisfied.[91]

#### Single fluid

[edit]

For a single substance, the definition of the molar Gibbs free energy can be written in the form  f = g − p v {\displaystyle f=g-pv} . Thus when  p {\displaystyle p} and  g {\displaystyle g} are constant, the function  f ( v ) {\displaystyle f(v)} is a straight line with slope  − p {\displaystyle -p} , and intercept  g {\displaystyle g} . Since the curve  f ( T R , v ) {\displaystyle f(T_{R},v)} has positive curvature everywhere when  T R ≥ T c {\displaystyle T_{R}\geq T_{\text{c}}} , the curve and the straight line will have a single tangent. However, for a subcritical  T R , f ( T R , v ) {\displaystyle T_{R},\,f(T_{R},v)} is not everywhere convex. With  p = p s ( T R ) {\displaystyle p=p_{\text{s}}(T_{R})} and a suitable value of  g {\displaystyle g} , the line will be tangent to  f ( T R , v ) {\displaystyle f(T_{R},v)} at the molar volume of each coexisting phase: saturated liquid  v f ( T R ) {\displaystyle v_{f}(T_{R})} and saturated vapor  v g ( T R ) {\displaystyle v_{g}(T_{R})} ; there will be a double tangent. Furthermore, each of these points is characterized by the same values of  g {\displaystyle g} ,  p {\displaystyle p} , and  T R . {\displaystyle T_{R}.} These are the same three specifications for coexistence that were used previously.[91]

Figure 8: The tangent line (black) to the curve  f r ( 0.875 , v r ) {\displaystyle f_{\text{r}}(0.875,v_{\text{r}})} (green) at the two points  v r f = 0.576 {\displaystyle v_{rf}=0.576} and  v r g = 2.71 {\displaystyle v_{rg}=2.71} . The line's slope, given by  ∂ v r f r = − p c v c / ( R T c ) p r s {\displaystyle \partial _{v_{\text{r}}}f_{\text{r}}=-p_{\text{c}}v_{\text{c}}/(RT_{\text{c}})p_{rs}} , is  − 0.215 {\displaystyle -0.215} corresponding to  p r s = 0.5730 {\displaystyle p_{rs}=0.5730} . This is the same data as isotherm  T r = 7 / 8 {\displaystyle T_{\text{r}}=7/8} in Fig. 1. The intercept on the line is  g {\displaystyle g} , but its numerical value is arbitrary due to a constant of integration.

Figure 8 depicts an evaluation of  f ( T R , v ) {\displaystyle f(T_{R},v)} as a green curve, with  v f {\displaystyle v_{f}} and  v g {\displaystyle v_{g}} marked by the left and right green circles, respectively. The region on the green curve for  v ≤ v f {\displaystyle v\leq v_{f}} corresponds to the liquid state. As  v {\displaystyle v} increases past  v f {\displaystyle v_{f}} , the curvature of  f {\displaystyle f} (proportional to  ∂ v ∂ v f = − ∂ v p {\displaystyle \partial _{v}\partial _{v}f=-\partial _{v}p} ) continually decreases. The inflection point, characterized by zero curvature, is a spinodal point; between  v f {\displaystyle v_{f}} and this point is the metastable superheated liquid. For further increases in  v {\displaystyle v} the curvature decreases to a minimum then increases to another (zero curvature) spinodal point; between these two spinodal points is the unstable region in which the fluid cannot exist in a homogeneous equilibrium state (represented by the dotted grey curve). With a further increase in  v {\displaystyle v} the curvature increases to a maximum at  v g {\displaystyle v_{g}} , where the slope is  p s {\displaystyle p_{\text{s}}} ; the region between this point and the second spinodal point is the metastable subcooled vapor. Finally, the region  v ≥ v g {\displaystyle v\geq v_{g}} is the vapor. In this region the curvature continually decreases until it is zero at infinitely large  v {\displaystyle v} . The double tangent line (solid black) that runs between  v f {\displaystyle v_{f}} and  v g {\displaystyle v_{g}} represents states that are stable but heterogeneous, not homogeneous solutions of the vdW equation.[91] The states above this line (with larger Helmholtz free energy) are either metastable or unstable.[91] The combined solid green-black curve in Figure 8 is the convex envelope of  f ( T R , v ) {\displaystyle f(T_{R},v)} , which is defined as the largest convex curve that is less than or equal to the function.[92]

For a vdW fluid, the molar Helmholtz potential is given by Eq (4a). This is, in reduced form,

f r = f / R T c = C u + T r ( c − C s − ln ⁡ [ T r c ( 3 v r − 1 ) ] ) − 9 / ( 8 v r ) {\displaystyle f_{r}=f/RT_{\text{c}}=C_{u}+T_{\text{r}}(c-C_{\text{s}}-\ln[T_{\text{r}}^{c}(3v_{\text{r}}-1)])-9/(8v_{\text{r}})}

with derivative  ∂ v r f r = − 3 T r / ( 3 v r − 1 ) + ( 9 / ( 8 v r ) 2 = − p r {\displaystyle \partial _{v_{\text{r}}}f_{\text{r}}=-3T_{\text{r}}}/({3v_{\text{r}}-1)+(9/(8v_{\text{r}})^{2}=-p_{\text{r}}} . A plot of this function  f r {\displaystyle f_{\text{r}}} , whose slope at each point is given by  p r {\displaystyle p_{\text{r}}} of the vdW equation, for the subcritical isotherm  T r = 7 / 8 {\displaystyle T_{\text{r}}=7/8} is shown in Figure 8 along with the line tangent to it at its two coexisting saturation points. The data illustrated in Figure 8 is the same as that shown in Figure 1 for this isotherm.[91]

This double tangent construction thus provides a graphical alternative to the Maxwell construction to establish the saturated liquid and vapor points on an isotherm.[91]

#### Binary fluid

[edit]

Van der Waals used the Helmholtz function because its properties could be easily extended to the binary fluid situation. In a binary mixture of vdW fluids, the Helmholtz potential is a function of two variables,  f ( T R , v , x ) {\displaystyle f(T_{R},v,x)} , where  x {\displaystyle x} is a composition variable (for example  x = x 2 {\displaystyle x=x_{2}} so  x 1 = 1 − x {\displaystyle x_{1}=1-x} ). In this case, there are three stability conditions:[93]

∂ 2 f ∂ v 2 > 0 ∂ 2 f ∂ x 2 > 0 ∂ 2 f ∂ v 2 ∂ 2 f ∂ x 2 − ( ∂ 2 f ∂ x ∂ v ) 2 > 0 {\displaystyle {\frac {\partial ^{2}f}{\partial v^{2}}}>0\qquad {\frac {\partial ^{2}f}{\partial x^{2}}}>0\qquad {\frac {\partial ^{2}f}{\partial v^{2}}}{\frac {\partial ^{2}f}{\partial x^{2}}}-\left({\frac {\partial ^{2}f}{\partial x\partial v}}\right)^{2}>0} and the Helmholtz potential is a surface (of physical interest in the region  0 ≤ x ≤ 1 {\displaystyle 0\leq x\leq 1} ). The first two stability conditions show that the curvature in each of the directions  v {\displaystyle v} and  x {\displaystyle x} are both non-negative for stable states, while the third condition indicates that stable states correspond to elliptic points on this surface.[94] Moreover, its limit  ∂ 2 f ∂ v 2 ∂ 2 f ∂ x 2 − ∂ 2 f ∂ x ∂ v = 0 {\displaystyle {\frac {\partial ^{2}f}{\partial v^{2}}}{\frac {\partial ^{2}f}{\partial x^{2}}}-{\frac {\partial ^{2}f}{\partial x\partial v}}=0} specifies the spinodal curves on the surface.

For a binary mixture, the Euler equation[95] can be written in the form  f = − p v + μ 1 x 1 + μ 2 x 2 = − p v + ( μ 2 − μ 1 ) x + μ 1 {\displaystyle {\begin{aligned}f&=-pv+\mu _{1}x_{1}+\mu _{2}x_{2}\\\&=-pv+(\mu _{2}-\mu _{1})x+\mu _{1}\end{aligned}}} where  μ j = ∂ x j f {\displaystyle \mu _{j}=\partial _{x_{j}}f} are the molar chemical potentials of each substance,  j = 1 , 2 {\displaystyle j=1,2} . For constant values of  p {\displaystyle p} ,  μ 1 {\displaystyle \mu _{1}} , and  μ 2 {\displaystyle \mu _{2}} , this equation is a plane with slopes  − p {\displaystyle -p} in the  v {\displaystyle v} direction,  μ 2 − μ 1 {\displaystyle \mu _{2}-\mu _{1}} in the  x {\displaystyle x} direction, and intercept  μ 1 {\displaystyle \mu _{1}} . As in the case of a single substance, here the plane and the surface can have a double tangent, and the locus of the coexisting phase points forms a curve on each surface. The coexistence conditions are that the two phases have the same  T {\displaystyle T} ,  p {\displaystyle p} ,  μ 2 − μ 1 {\displaystyle \mu _{2}-\mu _{1}} , and  μ 1 {\displaystyle \mu _{1}} ; the last two are equivalent to having the same  μ 1 {\displaystyle \mu _{1}} and  μ 2 {\displaystyle \mu _{2}} individually, which are just the Gibbs conditions for material equilibrium in this situation. The two methods of producing the coexistence surface are equivalent.[93]

Although this case is similar to that of a single fluid, here the geometry can be much more complex. The surface can develop a wave (called a plait or fold) in the  x {\displaystyle x} direction as well as the one in the  v {\displaystyle v} direction. Therefore, there can be two liquid phases that can be either miscible, or wholly or partially immiscible, as well as a vapor phase.[96][97] Despite a great deal of both theoretical and experimental work on this problem by van der Waals and his successors—work which produced much useful knowledge about the various types of phase equilibria that are possible in fluid mixtures[98]—complete solutions to the problem were only obtained after 1967, when the availability of modern computers made calculations of mathematical problems of this complexity feasible for the first time.[99] The results obtained were, in Rowlinson's words,[100]

> a spectacular vindication of the essential physical correctness of the ideas behind the van der Waals equation, for almost every kind of critical behavior found in practice can be reproduced by the calculations, and the range of parameters that correlate with the different kinds of behavior are intelligible in terms of the expected effects of size and energy.

### Mixing rules

[edit]

To obtain these numerical results, the values of the constants of the individual component fluids  a 11 , a 22 , b 11 , b 22 {\displaystyle a_{11},a_{22},b_{11},b_{22}} must be known. In addition, the effect of collisions between molecules of the different components, given by  a 12 {\displaystyle a_{12}} and  b 12 {\displaystyle b_{12}} , must also be specified. In the absence of experimental data, or computer modeling results to estimate their value, the empirical combining rules, geometric and algebraic means can be used, respectively:[101] a 12 = ( a 11 a 22 ) 1 / 2 and b 12 1 / 3 = ( b 11 1 / 3 + b 22 1 / 3 ) / 2. {\displaystyle a_{12}=(a_{11}a_{22})^{1/2}\qquad {\text{and}}\qquad b_{12}^{1/3}=(b_{11}^{1/3}+b_{22}^{1/3})/2.}

These relations correspond to the empirical combining rules for the intermolecular force constants,  ϵ 12 = ( ϵ 11 ϵ 22 ) 1 / 2 and σ 12 = ( σ 11 + σ 22 ) / 2 , {\displaystyle \epsilon _{12}=(\epsilon _{11}\epsilon _{22})^{1/2}\qquad {\text{and}}\qquad \sigma _{12}=(\sigma _{11}+\sigma _{22})/2,} the first of which follows from a simple interpretation of the dispersion forces in terms of polarizabilities of the individual molecules, while the second is exact for rigid molecules.[102] Using these empirical combining rules to generalize for  n {\displaystyle n} fluid components, the quadratic mixing rules for the material constants are:[89] a x = ∑ i = 1 n ∑ j = 1 n ( a i i a j j ) 1 / 2 x i x j which can be written as a x = ( ∑ i = 1 n a i i 1 / 2 x i ) 2 {\displaystyle a_{x}=\sum _{i=1}^{n}\sum _{j=1}^{n}(a_{ii}a_{jj})^{1/2}x_{i}x_{j}\quad {\text{which can be written as}}\quad a_{x}={\Big (}\sum _{i=1}^{n}a_{ii}^{1/2}x_{i}{\Big )}^{2}} b x = ( 1 / 8 ) ∑ i = 1 n ∑ j = 1 n ( b i i 1 / 3 + b j j 1 / 3 ) 3 x i x j {\displaystyle b_{x}=(1/8)\sum _{i=1}^{n}\sum _{j=1}^{n}(b_{ii}^{1/3}+b_{jj}^{1/3})^{3}x_{i}x_{j}}

These expressions come into use when mixing gases in proportion, such as when producing tanks of air for diving[103] and managing the behavior of fluid mixtures in engineering applications. However, more sophisticated mixing rules are often necessary, to obtain satisfactory agreement with reality over the wide variety of mixtures encountered in practice.[104][105]

Another method of specifying the vdW constants, pioneered by W.B. Kay and known as Kay's rule,[106] specifies the effective critical temperature and pressure of the fluid mixture by  T c x = ∑ i = 1 n T c i x i and p c x = ∑ i = 1 n p c i x i . {\displaystyle T_{{\text{c}}x}=\sum _{i=1}^{n}T_{{\text{c}}i}x_{i}\qquad {\text{and}}\qquad p_{{\text{c}}x}=\sum _{i=1}^{n}\,p_{{\text{c}}i}x_{i}.}

In terms of these quantities, the vdW mixture constants are  a x = ( 3 4 ) 3 ( R T c x ) 2 p c x b x = ( 1 2 ) 3 R T c x p c x {\displaystyle a_{x}=\left({\frac {3}{4}}\right)^{3}{\frac {(RT_{{\text{c}}x})^{2}}{p_{{\text{c}}x}}}\qquad \qquad b_{x}=\left({\frac {1}{2}}\right)^{3}{\frac {RT_{{\text{c}}x}}{p_{{\text{c}}x}}}} which Kay used as the basis for calculations of the thermodynamic properties of mixtures. Kay's idea was adopted by T. W. Leland, who applied it to the molecular parameters  ϵ , σ {\displaystyle \epsilon ,\sigma } , which are related to  a , b {\displaystyle a,b} through  T c , p c {\displaystyle T_{\text{c}},p_{\text{c}}} by  a ∝ ϵ σ 3 {\displaystyle a\propto \epsilon \sigma ^{3}} and  b ∝ σ 3 {\displaystyle b\propto \sigma ^{3}} . Using these together with the quadratic mixing rules for  a , b {\displaystyle a,b} produces  σ x 3 = ∑ i = i n ∑ j = 1 n σ i j 3 x i x j and ϵ x = [ ∑ i = 1 n ∑ j = 1 n ϵ i j σ i j 3 x i x j ] [ ∑ i = i n ∑ j = 1 n σ i j 3 x i x j ] − 1 {\displaystyle \sigma _{x}^{3}=\sum _{i=i}^{n}\sum _{j=1}^{n}\,\sigma _{ij}^{3}x_{i}x_{j}\qquad {\text{and}}\qquad \epsilon _{x}=\left[\sum _{i=1}^{n}\sum _{j=1}^{n}\epsilon _{ij}\sigma _{ij}^{3}x_{i}x_{j}\right]\left[\sum _{i=i}^{n}\sum _{j=1}^{n}\,\sigma _{ij}^{3}x_{i}x_{j}\right]^{-1}} which is the van der Waals approximation expressed in terms of the intermolecular constants.[107][108] This approximation, when compared with computer simulations for mixtures, are in good agreement over the range  1 / 2 < ( σ 11 / σ 22 ) 3 < 2 {\displaystyle 1/2<(\sigma _{11}/\sigma _{22})^{3}<2} , namely for molecules of similar diameters. Rowlinson said of this approximation, "It was, and indeed still is, hard to improve on the original van der Waals recipe when expressed in [this] form".[109]

## Validity of the equation

[edit]

Since van der Waals presented his thesis, "[m]any derivations, pseudo-derivations, and plausibility arguments have been given" for it.[110] However, no mathematically rigorous derivation of the equation over its entire range of molar volume that begins from a statistical mechanical principle exists. Indeed, such a proof is not possible, even for hard spheres.[111][112][113][114][115] Goodstein writes, "Obviously the value of the van der Waals equation rests principally on its empirical behavior rather than its theoretical foundation."[7]

Although the use of the vdW equation is not justified mathematically, it has empirical validity. Its various applications in this region that attest to this, both qualitative and quantitative, have been described previously in this article. This point was also made by Alder, et al. who, at a conference marking the 100th anniversary of van der Waals' thesis, noted that:[116]

> It is doubtful whether we would celebrate the centennial of the Van der Waals equation if it were applicable only under circumstances where it has been proven to be rigorously valid. It is empirically well established that many systems whose molecules have attractive potentials that are neither long-range nor weak conform nearly quantitatively to the Van der Waals model. An example is the theoretically much studied system of Argon, where the attractive potential has only a range half as large as the repulsive core.

They continued by saying that this model has "validity down to temperatures below the critical temperature, where the attractive potential is not weak at all but, in fact, comparable to the thermal energy." They also described its application to mixtures "where the Van der Waals model has also been applied with great success. In fact, its success has been so great that not a single other model of the many proposed since, has equalled its quantitative predictions,[117] let alone its simplicity."[118]

Engineers have made extensive use of this empirical validity, modifying the equation in numerous ways (by one account there have been some 400 cubic equations of state produced)[119] to manage the liquids,[120] and gases of pure substances and mixtures,[121] that they encounter in practice.

This situation has been aptly described by Boltzmann:[122]

> ... van der Waals has given us such a valuable tool that it would cost us much trouble to obtain by the subtlest deliberations a formula that would really be more useful than the one that van der Waals found by inspiration, as it were.

## Notes

[edit]

1. ^ van der Waals 1967, p. 254
2. ^ Goodstein 1985, pp. 443–463
3. ^ Epstein 1937, p. 9
4. ^ Boltzmann 1995, p. 231
5. ^ van der Waals 1873, p. 56
6. ^ a b Sommerfeld 1956, p. 55
7. ^ a b Goodstein 1985, p. 446
8. ^ van der Waals 2004, pp. 168–172
9. ^ Boltzmann 1995, p. 221–224
10. ^ van der Waals 2004, p. 173
11. ^ a b Hirschfelder, Curtis & Bird 1964, pp. 31–34
12. ^ Goodstein 1985, pp. 250, 263
13. ^ a b Tien & Lienhard 1979, pp. 250, 251
14. ^ a b c Epstein 1937, p. 10
15. ^ Boltzmann, L. Enzykl. der Mathem. Wiss., V, (1), 550
16. ^ a b Sommerfeld 1956, p. 66
17. ^ Weinberg 2021, pp. 4–5
18. ^ Weinberg 2021, p. 33
19. ^ Gibbs 1948, pp. vii–xii
20. ^ van der Waals 1873
21. ^ van der Waals 2004, pp. 121–240
22. ^ Boltzmann 1995, p. 218
23. ^ Andrews 1869
24. ^ Klein 1974, p. 31
25. ^ van der Waals 2004, pp. 125, 191–194
26. ^ Goodstein 1985, pp. 450–451
27. ^ Boltzmann 1995, pp. 232–233
28. ^ Goodstein 1985, p. 452
29. ^ Sommerfeld 1956, pp. 55–68
30. ^ Johnston 2014, pp. 5, 7, 12, 18, 25
31. ^ Goodstein 1985, pp. 443–452
32. ^ a b c Dong & Lienhard 1986, pp. 158–159
33. ^ Sommerfeld 1956, pp. 56–57
34. ^ Goodstein 1985, p. 449
35. ^ Boltzmann 1995, pp. 237–238
36. ^ Lienhard, Shamsundar & Biney 1986, pp. 297–298
37. ^ Boltzmann 1995, pp. 239–240
38. ^ Barenblatt 1979, pp. 13–23
39. ^ Sommerfeld 1956, p. 57
40. ^ Peck 1982, p. 446
41. ^ a b Johnston 2014, p. 6
42. ^ Lienhard 1986, pp. 172–173
43. ^ Whitman 2023, pp. 155, 202–204
44. ^ a b Moran and Shapiro 2000, pp. 574 sfn error: no target: CITEREFMoran_and_Shapiro2000 (help)
45. ^ Sommerfeld, pp. 58–59 sfn error: no target: CITEREFSommerfeld (help)
46. ^ Whitman 2023, pp. 155, 204
47. ^ a b Johnston 2014, p. 7
48. ^ a b Whitman 2023, p. 192
49. ^ Whitman 2023, p. 257
50. ^ Johnston 2014, p. 5.
51. ^ a b c d Johnston 2014, p. 12
52. ^ Callen 1962, pp. 117–130 sfn error: no target: CITEREFCallen1962 (help)
53. ^ Moran and Shapiro 2000, p. 578 sfn error: no target: CITEREFMoran_and_Shapiro2000 (help)
54. ^ Whitman 2014, pp. 43–45, 98 sfn error: no target: CITEREFWhitman2014 (help)
55. ^ Sommerfeld 1956, p. 56
56. ^ a b Johnston 2014, p. 11
57. ^ a b Moran and Shapiro 2000, p. 580 sfn error: no target: CITEREFMoran_and_Shapiro2000 (help)
58. ^ Whitman 2023, p. 204
59. ^ Johnston 2014, p. 3
60. ^ Lienhard, Shamsundahr and Biney 1986, p. 299 sfn error: no target: CITEREFLienhard,_Shamsundahr_and_Biney1986 (help)
61. ^ a b Callen 1960, pp. 37–44
62. ^ Shamsundar & Lienhard 1983, pp. 878, 879
63. ^ Barrufet & Eubank 1989, p. 170
64. ^ Lekner 1982, pp. 161–162
65. ^ Johnston 2014, pp. 16–18
66. ^ Boltzmann 1995, pp. 248–250
67. ^ Lienhard, Shamsundar & Biney 1986, p. 297
68. ^ Tien & Lienhard 1979, p. 254
69. ^ a b Callen 1960, pp. 146–156
70. ^ Rowlinson (ed) in van der Waals 2004, p. 22 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals2004 (help)
71. ^ Pitzer et al. 1955, p. 3433
72. ^ Sommerfeld 1956, pp. 61–63
73. ^ Sommerfeld 1956, pp. 60–62
74. ^ a b Sommerfeld 1956, p. 61
75. ^ Sommerfeld 1956, p. 62, Fig.8
76. ^ Van Wylen & Sonntag 1973, p. 49
77. ^ Johnston 2014, p. 10
78. ^ Su 1946, p. 803
79. ^ Moran & Shapiro 2000, p. 113
80. ^ Hirschfelder, Curtis, and Bird 1964, pp. 137–145 sfn error: no target: CITEREFHirschfelder,_Curtis,_and_Bird1964 (help)
81. ^ Hill 1986, pp. 261–266
82. ^ Hirschfelder, Curtis, and Bird 1964, pp. 156–157, 158–159 sfn error: no target: CITEREFHirschfelder,_Curtis,_and_Bird1964 (help)
83. ^ Boltzmann 1995, pp. 353–356
84. ^ Rowlinson (ed) in van der Waals 2004, pp. 20–22 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals2004 (help)
85. ^ van der Waals 2004, pp. 243–282
86. ^ Lorentz 1881, pp. 127, 134, 600 sfn error: no target: CITEREFLorentz1881 (help)
87. ^ Rowlinson (ed) in van der Waals 2004, p. 68 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals2004 (help)
88. ^ van der Waals 2004, p. 244
89. ^ a b Redlich & Kwong 1949
90. ^ Callen 1986, p. 105 sfn error: no target: CITEREFCallen1986 (help)
91. ^ a b c d e f van der Waals 2004, pp. 245–247
92. ^ Lebowitz 1974, p. 52
93. ^ a b van der Waals 2004, pp. 247–250
94. ^ Kreyszig 1959, pp. 124–128
95. ^ Callen 1960, pp. 47–48
96. ^ Rowlinson (ed) in van der Waals 2004, pp. 23–27 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals2004 (help)
97. ^ van der Waals 2004, pp. 253–258
98. ^ DeBoer 1974, pp. 7–16
99. ^ Rowlinson (ed) in van der Waals 2004, pp. 23–27, 64–66 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals2004 (help)
100. ^ Rowlinson (ed) in van der Waals 1984, p. 66 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals1984 (help)
101. ^ Hirschfelder, Curtis & Bird 1964, pp. 252–253
102. ^ Hirschfelder, Curtis & Bird 1964, pp. 168–169
103. ^ Hewitt
104. ^ Valderrama 2003, pp. 1308–1312
105. ^ Kontogeorgis, Privat & Jaubert 2019, pp. 4626–4633
106. ^ Niemeyer
107. ^ Leland, Rowlinson & Sather 1968, p. 1447
108. ^ Rowlinson (ed) in van der Waals 1984, pp. 69–70 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals1984 (help)
109. ^ Rowlinson (ed) in van der Waals 1984, p. 70 sfn error: no target: CITEREFRowlinson_(ed)_in_van_der_Waals1984 (help)
110. ^ Goodstein 1985, p. 443
111. ^ Korteweg 1891b, p. 277
112. ^ Tonks 1936, pp. 962–963
113. ^ van Hove 1949, pp. 957–959
114. ^ Kac, Uhlenbeck & Hemmer 1963, p. 224
115. ^ Lebowitz 1974, pp. 50, 52–57
116. ^ Alder, Alley & Rigby 1974, p. 143
117. ^ Singer, J.V.R., and Singer, K., Mol. Phys.(1972), 24, 357; McDonald, J.R., (1972), 24, 391
118. ^ Alder, Alley & Rigby 1974, p. 144
119. ^ Valderrama 2003, p. 1606
120. ^ Vera & Prausnitz 1972, p. 7-10
121. ^ Kontogeorgis, Privat & Jaubert 2019, pp. 4626–4629
122. ^ Boltzmann 1995, p. 356

## References

[edit]

* Alder, B. J.; Alley, W. E.; Rigby, M. (1974). "Correction to the van der Waals model for mixtures and for the diffusion coefficient". Physica. 74 (1): 143–155\. Bibcode:1974Phy....73..143A. doi:10.1016/0031-8914(74)90231-6.143-155&rft.date=1974&rft_id=info:doi/10.1016/0031-8914(74)90231-6&rft_id=info:bibcode/1974Phy....73..143A&rft.aulast=Alder&rft.aufirst=B.+J.&rft.au=Alley,+W.+E.&rft.au=Rigby,+M.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Andrews, T. (1869). "On the Continuity of the Gaseous and Liquid States of Matter". Philosophical Transactions of the Royal Society of London. 159: 575–590.575-590&rft.date=1869&rft.aulast=Andrews&rft.aufirst=T.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Barenblatt, G.I. (1979) [1978], Similarity, Self-Similarity, and Intermediate Asymptotics, translated by Stein, Norman, NY and London: Milton
* Barrufet, M.A.; Eubank, P.T. (1989). "Generalized Saturation Properties of Pure Fluids Via Cubic Equations of State". Chemical Engineering Education. 23 (3): 168–175.168-175&rft.date=1989&rft.aulast=Barrufet&rft.aufirst=M.A.&rft.au=Eubank,+P.T.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Boltzmann, L. (1995) [1896 (Part I), 1898 (Part II) in German]. Lectures on Gas Theory. Translated by Brush, S.G. NY: Dover.
* Callen, H.B. (1960). Thermodynamics. NY: John Wiley ans Sons.
* DeBoer, J. (1974). "Van der Waals in his time and the present revival opening address". Physica. 73 (1): 1–27\. Bibcode:1974Phy....73....1D. doi:10.1016/0031-8914(74)90223-7.1-27&rft.date=1974&rft_id=info:doi/10.1016/0031-8914(74)90223-7&rft_id=info:bibcode/1974Phy....73....1D&rft.aulast=DeBoer&rft.aufirst=J.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Dong, W.G.; Lienhard, J.H. (1986). "Corresponding States Correlation of Saturated and Metastable Properties". Can J Chem Eng. 64: 158–161\. doi:10.1002/cjce.5450640123.158-161&rft.date=1986&rft_id=info:doi/10.1002/cjce.5450640123&rft.aulast=Dong&rft.aufirst=W.G.&rft.au=Lienhard,+J.H.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Epstein, P.S. (1937). Textbook of Thermodynamics. NY: John Wiley and Sons.
* Gibbs, J.W. (1948) [1901]. The Collected Works of J. Willard Gibbs Volume II Part One Elementary Principles in Statistical Mechanics. New Haven: Yale University Press.
* Goodstein, D.L. (1985) [1975]. States of Matter. NY: Dover.
* Hewitt, Nigel. "Who was Van der Waals anyway and what has he to do with my Nitrox fill?". Maths for Divers. Archived from the original on 11 March 2020. Retrieved 1 February 2019.
* Hill, Terrell L. (1986). Statistical Thermodynamics. NY: Dover.
* Hirschfelder, J. O.; Curtis, C. F.; Bird, R. B. (1964). Mollecular Theory of Gases and Liquids, corrected printing. NY: John Wiley and Sons, Inc.
* Johnston, D.C. (2014). Advances in Thermodynamics of the van der Waals Fluid. arXiv:1402.1205. Bibcode:2014atvd.book.....J. doi:10.1088/978-1-627-05532-1. ISBN 978-1-627-05532-1.
* Kac, M.; Uhlenbeck, G.E.; Hemmer, P.C. (1963). "On the van der Waals Theory of the Vapor-Liquid Equilibrium. 1. Discussion of a One-Dimensional Model". J. Math. Phys. 4 (2): 216–228\. Bibcode:1963JMP.....4..216K. doi:10.1063/1.1703946.216-228&rft.date=1963&rft_id=info:doi/10.1063/1.1703946&rft_id=info:bibcode/1963JMP.....4..216K&rft.aulast=Kac&rft.aufirst=M.&rft.au=Uhlenbeck,+G.E.&rft.au=Hemmer,+P.C.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Kreyszig, E. (1959). Differential Geometry. Toronto: University of Toronto Press.
* Klein, M. J. (1974). "The Historical Origins of the Van der Waals Equation". Physica. 73 (1): 28–47\. Bibcode:1974Phy....73...28K. doi:10.1016/0031-8914(74)90224-9.28-47&rft.date=1974&rft_id=info:doi/10.1016/0031-8914(74)90224-9&rft_id=info:bibcode/1974Phy....73...28K&rft.aulast=Klein&rft.aufirst=M.+J.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Kontogeorgis, G.M.; Privat, R.; Jaubert, J-N.J. (2019). "Taking Another Look at the van der Waals Equation of State – Almost 150 Years Later". J. Chem. Eng. Data. 64 (11): 4619–4637\. doi:10.1021/acs.jced.9b00264.4619-4637&rft.date=2019&rft_id=info:doi/10.1021/acs.jced.9b00264&rft.aulast=Kontogeorgis&rft.aufirst=G.M.&rft.au=Privat,+R.&rft.au=Jaubert,+J-N.J.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Korteweg, D.T. (1891a). "On Van Der Waals Isothermal Equation". Nature. 45 (1155): 152–154\. Bibcode:1891Natur..45..152K. doi:10.1038/045152a0.152-154&rft.date=1891&rft_id=info:doi/10.1038/045152a0&rft_id=info:bibcode/1891Natur..45..152K&rft.aulast=Korteweg&rft.aufirst=D.T.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Korteweg, D.T. (1891b). "On Van Der Waals Isothermal Equation". Nature. 45 (1160): 277. doi:10.1038/045277a0.
* Lebowitz, J.L. (1974). "Exact Derivation of the Van Der Waals Equation". Physica. 73 (1): 48–60\. Bibcode:1974Phy....73...48L. doi:10.1016/0031-8914(74)90225-0.48-60&rft.date=1974&rft_id=info:doi/10.1016/0031-8914(74)90225-0&rft_id=info:bibcode/1974Phy....73...48L&rft.aulast=Lebowitz&rft.aufirst=J.L.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Lebowitz, J.L.; Penrose, O. (1966). "Rigorous Treatment of the Van der Waals–Maxwell Theory of the Liquid–Vapor Transition". Jour Math Phys. 7 (1): 98–113\. Bibcode:1966JMP.....7...98L. doi:10.1063/1.1704821.98-113&rft.date=1966&rft_id=info:doi/10.1063/1.1704821&rft_id=info:bibcode/1966JMP.....7...98L&rft.aulast=Lebowitz&rft.aufirst=J.L.&rft.au=Penrose,+O.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Lekner, J. (1982). "Parametric solution of the van der Waals liquid–vapor coexistence curve". Am. J. Phys. 50 (2): 161–163\. Bibcode:1982AmJPh..50..161L. doi:10.1119/1.12877.161-163&rft.date=1982&rft_id=info:doi/10.1119/1.12877&rft_id=info:bibcode/1982AmJPh..50..161L&rft.aulast=Lekner&rft.aufirst=J.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Leland, T. W.; Rowlinson, J.S.; Sather, G.A. (1968). "Statistical thermodynamics of mixtures of molecules of different sizes". Trans. Faraday Soc. 64: 1447–1460.1447-1460&rft.date=1968&rft.aulast=Leland&rft.aufirst=T.+W.&rft.au=Rowlinson,+J.S.&rft.au=Sather,+G.A.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Lienhard, J.H. (1986). "The Properties and Behavior of Superheated Liquids". Lat. Am. J. Heat and Mass Transfer. 10: 169–187.169-187&rft.date=1986&rft.aulast=Lienhard&rft.aufirst=J.H.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Lienhard, J.H; Shamsundar, N.; Biney, P.O. (1986). "Spinodal Lines and Equations of State: A Review". Nuclear Engineering and Design. 95: 297–314\. Bibcode:1986NuEnD..95..297L. doi:10.1016/0029-5493(86)90056-7.297-314&rft.date=1986&rft_id=info:doi/10.1016/0029-5493(86)90056-7&rft_id=info:bibcode/1986NuEnD..95..297L&rft.aulast=Lienhard&rft.aufirst=J.H&rft.au=Shamsundar,+N.&rft.au=Biney,+P.O.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Maxwell, J.C. (1875). "On the Dynamical Evidence of the Molecular Constitution of Bodies". Nature. 11 (279): 357–359\. Bibcode:1875Natur..11..357C. doi:10.1038/011357a0.357-359&rft.date=1875&rft_id=info:doi/10.1038/011357a0&rft_id=info:bibcode/1875Natur..11..357C&rft.aulast=Maxwell&rft.aufirst=J.C.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Moran, M.J.; Shapiro, H.N. (2000). Fundamentals of Engineering Thermodynamics 4th Edition. NY: McGraw-Hill.
* Niemeyer, Kyle. "Mixture properties". Computational Thermodynamics. Archived from the original on 2 April 2024. Retrieved 2 April 2024.
* Peck, R.E. (1982). "The Assimilation of van der Waals Equation in the Corresponding States Family". Can. J. Chem. Eng. 60: 446–449\. doi:10.1002/cjce.5450600319.446-449&rft.date=1982&rft_id=info:doi/10.1002/cjce.5450600319&rft.aulast=Peck&rft.aufirst=R.E.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Pitzer, K.S.; Lippman, D.Z.; Curl, R.F.; Huggins, C.M.; Peterson, D.E. (1955). "The Volumetric and Thermodynamic Properties of Fluids. II. Compressibility Factor, Vapor Pressure and Entropy of Vaporization". J. Am. Chem. Soc. 77 (13): 3433–3440\. Bibcode:1955JAChS..77.3433P. doi:10.1021/ja01618a002.3433-3440&rft.date=1955&rft_id=info:doi/10.1021/ja01618a002&rft_id=info:bibcode/1955JAChS..77.3433P&rft.aulast=Pitzer&rft.aufirst=K.S.&rft.au=Lippman,+D.Z.&rft.au=Curl,+R.F.&rft.au=Huggins,+C.M.&rft.au=Peterson,+D.E.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Redlich, O.; Kwong, J. N. S. (1949). "On the Thermodynamics of Solutions. V. An Equation of State. Fugacities of Gaseous Solutions" (PDF). Chemical Reviews. 44 (1): 233–244\. doi:10.1021/cr60137a013. PMID 18125401. Retrieved 2 April 2024.233-244&rft.date=1949&rft_id=info:doi/10.1021/cr60137a013&rft_id=info:pmid/18125401&rft.aulast=Redlich&rft.aufirst=O.&rft.au=Kwong,+J.+N.+S.&rft_id=https://www3.nd.edu/~powers/ame.20231/redlichkwong1949.pdf&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Shamsundar, N.; Lienhard, J.H. (1983). "Saturation and Metastable Properties of the van der Waals Fluid". Can J Chem Eng. 61 (6): 876–880\. doi:10.1002/cjce.5450610617.876-880&rft.date=1983&rft_id=info:doi/10.1002/cjce.5450610617&rft.aulast=Shamsundar&rft.aufirst=N.&rft.au=Lienhard,+J.H.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Sommerfeld, A. (1956). Bopp, F.; Meixner, J. (eds.). Thermodynamics and Statistical Mechanics – Lectures on Theoretical Physics Volume V. Translated by Kestin, J. NY: Academic Press. ISBN 978-0-323-13773-7.
* Strutt, J.W., 3rd Baron Rayleigh (1891). "On the Virial of a System of Hard Colliding Bodies". Nature. 45 (1152): 80–82\. Bibcode:1891Natur..45...80R. doi:10.1038/045080a0.80-82&rft.date=1891&rft_id=info:doi/10.1038/045080a0&rft_id=info:bibcode/1891Natur..45...80R&rft.aulast=Strutt&rft.aufirst=J.W.,+3rd+Baron+Rayleigh&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">`{{cite journal}}`: CS1 maint: multiple names: authors list (link) CS1 maint: numeric names: authors list (link)
* Su, G.J. (1946). "Modified Law of Corresponding States for Real Gases". Ind. Eng. Chem. 38 (8): 803–806\. doi:10.1021/ie50440a018.803-806&rft.date=1946&rft_id=info:doi/10.1021/ie50440a018&rft.aulast=Su&rft.aufirst=G.J.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Tien, C.L.; Lienhard, J.H. (1979). Statistical Thermodynamics Revised Printing. NY: Hemisphere Publishing. Bibcode:1979wdch.book.....T.
* Tonks, L. (1936). "The Complete Equation of State of One, Two, and Three-Dimensional Gases of Hard Elastic Spheres". Phys. Rev. 50 (10): 955–963\. Bibcode:1936PhRv...50..955T. doi:10.1103/PhysRev.50.955.955-963&rft.date=1936&rft_id=info:doi/10.1103/PhysRev.50.955&rft_id=info:bibcode/1936PhRv...50..955T&rft.aulast=Tonks&rft.aufirst=L.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Valderrama, J.O. (2003). "The State of the Cubic Equations of State". Ind. Chem. Eng. Res. 42 (8): 1603–1618\. doi:10.1021/ie020447b.1603-1618&rft.date=2003&rft_id=info:doi/10.1021/ie020447b&rft.aulast=Valderrama&rft.aufirst=J.O.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* van der Waals, J.D. (1873). Over de Continuïteit van den Gas en Vloeistoftoestand (Ph.D. thesis). Leiden Univ.
* van der Waals, J.D. (2004) [1984]. Rowlinson, J.S. (ed.). On the Continuity of the Gaseous and Liquid States, edited and with an Introduction by J.S. Rowlinson. NY: Dover Phoenix Editions.
* van der Waals, Johannes D. (1967). "The Equation of State of Gases and Liquids". in Nobel Lectures, Physics 1901–1921. Amsterdam: Elsevier. pp. 254–265.254-265&rft.pub=Elsevier&rft.date=1967&rft.aulast=van+der+Waals&rft.aufirst=Johannes+D.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* van Hove, L. (1949). "Quelques Proprieties Generales De L'Integrale De Configuration D'Un Systeme De Particules Avec Interaction". Physica. 15 (11–12): 951–961\. Bibcode:1949Phy....15..951V. doi:10.1016/0031-8914(49)90059-2.11–12&rft.pages=951-961&rft.date=1949&rft_id=info:doi/10.1016/0031-8914(49)90059-2&rft_id=info:bibcode/1949Phy....15..951V&rft.aulast=van+Hove&rft.aufirst=L.&rft_id=https://dipot.ulb.ac.be/dspace/bitstream/2013/182044/4/c315aefe-4bb4-4bb6-9ac7-a013ab2fad89.txt&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Van Wylen, G.J.; Sonntag, R.E. (1973). Fundamentals of Classical Thermodynamics Second Edition. NY: John Wiley ans Sons.
* Vera, J.H.; Prausnitz, J.M. (1972). "Generalized van der Waals Theory for Dense Fluids". Chem. Eng. Jour. 3: 1–13\. doi:10.1016/0300-9467(72)85001-9.1-13&rft.date=1972&rft_id=info:doi/10.1016/0300-9467(72)85001-9&rft.aulast=Vera&rft.aufirst=J.H.&rft.au=Prausnitz,+J.M.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Weinberg, S. (2021). Foundations of Modern Physics. Cambridge: Cambridge University Press.
* Whitman, A.M. (2023). Thermodynamics: Basic Principles and Engineering Applications 2nd Edition. NY: Springer.

## See also

[edit]

* Gas laws
* Ideal gas
* Inversion temperature
* Iteration
* Maxwell construction
* Real gas
* Theorem of corresponding states
* Van der Waals constants (data page)
* Redlich–Kwong equation of state

## Further reading

[edit]

* Chandler, David (1987). Introduction to Modern Statistical Mechanics. Oxford: Oxford University Press. pp. 287–295\. ISBN 0195042778.287-295&rft.pub=Oxford+University+Press&rft.date=1987&rft.isbn=0195042778&rft.aulast=Chandler&rft.aufirst=David&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Cross, Michael (2004), "Lecture 3: First Order Phase Transitions" (PDF), Physics 127: Statistical Physics, Second Term, Pasadena, California: Division of Physics, Mathematics, and Astronomy, California Institute of Technology.
* Dalgarno, A.; Davison, W.D. (1966). "The Calculation of Van Der Waals Interactions". Advances in Atomic and Molecular Physics. 2: 1–32\. Bibcode:1966AdAMP...2....1D. doi:10.1016/S0065-2199(08)60216-X. ISBN 9780120038022.1-32&rft.date=1966&rft_id=info:doi/10.1016/S0065-2199(08)60216-X&rft_id=info:bibcode/1966AdAMP...2....1D&rft.isbn=9780120038022&rft.aulast=Dalgarno&rft.aufirst=A.&rft.au=Davison,+W.D.&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">
* Kittel, Charles; Kroemer, Herbert (1980). Thermal Physics (Revised ed.). New York: Macmillan. pp. 287–295\. ISBN 0716710889.287-295&rft.edition=Revised&rft.pub=Macmillan&rft.date=1980&rft.isbn=0716710889&rft.aulast=Kittel&rft.aufirst=Charles&rft.au=Kroemer,+Herbert&rfr_id=info:sid/en.wikipedia.org:Van+der+Waals+equation" class="Z3988">