View source for Van der Waals equation

{{Use American English|date=December 2024}}
{{Short description|Gas equation of state which accounts for non-ideal gas behavior}}
{{Use dmy dates|date=December 2024}}

{{Too big|date=February 2025}}

The '''van der Waals equation''' is a mathematical formula that describes the behavior of [[Real gas|real gases]]. It is an [[equation of state]] that relates the [[pressure]], [[temperature]], and [[molar volume]] in a [[fluid]]. The equation modifies the [[ideal gas law]] 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 [[Thesis|doctoral thesis]]. Van der Waals based the equation on the idea that fluids are composed of discrete [[Particle|particles]], which few scientists believed existed. However, the equation accurately predicted the behavior of a fluid around its [[Critical point (thermodynamics)|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]].{{sfn|ps=|van der Waals|1967}} Today the equation is recognized as an important model of phase change processes.{{sfn|ps=|Goodstein|1985|pp=443-463}} Van der Waals also adapted his equation so that it applied to a binary mixture of fluids. He, and others, then used the modified equation to discover a host of important facts about the phase equilibria of such fluids.{{sfn|ps=|DeBoer|1974|pp=7-16}} This application, expanded to treat multi-component mixtures, has extended the predictive ability of the equation to fluids of industrial and commercial importance. In this arena it has spawned many similar equations in a continuing attempt by engineers to improve their ability to understand and manage these fluids;{{sfn|ps=|Valderrama|2010|pp=415-420}} it remains relevant to the present.{{sfn|ps=|Kontogeorgis|Privat|Jaubert|2019|pp=4619-4637}}

== Behavior of the equation ==

One way to write the van der Waals equation is:{{sfn|ps=|van der Waals|1967|p=174}}{{sfn|ps=|Epstein|1937|p=9}}{{sfn|ps=|Boltzmann|1995|p=231}}
<math display="block">p=\frac{RT}{v-b}-\frac {a}{v^2} ,</math>
where <math>p</math> is pressure, <math>R</math> is the [[universal gas constant]], <math>T</math> is temperature, <math>v</math> is [[molar volume]], and <math>a</math> and <math>b</math> are experimentally determinable, substance-specific constants. Molar volume is given by {{tmath|1= v=VN_\text{A}/N }}, where <math>N_\text{A}</math> is the [[Avogadro constant]], <math>V</math> is the [[volume]], and <math>N</math> is the number of molecules (the ratio <math>N/N_\text{A}</math> is the [[amount of substance]], a physical quantity with the base unit [[Mole (unit)|mole]]).

When van der Waals created his equation, few scientists believed that fluids were composed of rapidly moving particles. Moreover, those who thought so had no knowledge of the atomic/molecular structure. The simplest conception of a particle, and the easiest to model mathematically, was a hard sphere; this is what van der Waals used. In that case, two particles of diameter <math>\sigma</math> would come into contact when their centers were a distance <math>\sigma</math> apart; hence the center of the one was excluded from a spherical volume equal to <math>4\pi\sigma^3/3</math> about the other. That is 8 times {{tmath|1= V_0 }}, the volume of each particle of radius {{tmath|1= \sigma/2 }}, but there are 2 particles which gives 4 times the volume per particle. The total excluded volume is then {{tmath|1= B=4NV_0 }}; that is, 4 times the volume of all the particles. Van der Waals and his contemporaries used an alternative, but equivalent, analysis based on the mean free path between molecular collisions that gave this result.{{sfn|ps=|van der Waals|1967|pp=168-172}}{{sfn|ps=|Boltzmann|1995|p=221–224}} From the fact that the volume fraction of particles, <math>f=(V-B)/V</math> must be positive, van der Waals noted that as <math>N</math> becomes larger the factor 4 must decrease (for [[close-packing of equal spheres|spheres]] there is a known minimum {{tmath|1= 0<f_\mathrm{min}<f\le 1 }}), but he was never able to determine the nature of the decrease.{{sfn|ps=|van der Waals|1967|p=172}}{{sfn|ps=|van der Waals|1967|p=256}}  

The constant {{tmath|1= b=BN_\text{A} /N }}, and has [[Dimensional Analysis|dimension]] of molar volume, [v].<!-- [  ] is indeed the correct notation for dimensional analysis. --> The constant <math>a</math> expresses the strength of the hypothesized interparticle attraction. Van der Waals only had as a model Newton's law of gravitation, in which two particles are attracted in proportion to the product of their masses. Thus he argued that in his case the attractive pressure was proportional to the square of the density.{{sfn|ps=|van der Waals|1967|p=173}} The proportionality constant, {{tmath|1= a }}, when written in the form used above, has the dimension [pv<sup>2</sup>] (pressure times molar volume squared), which is also molar energy times molar volume.  

[[file:pair potentials.png|thumb|upright=1.6|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.{{sfn|ps=|Hirschfelder|Curtis|Bird|1964|pp=31-34}}]]

The intermolecular force was later conveniently described by the negative derivative of a [[pair potential]] function. For spherically symmetric particles, this is most simply a function of separation distance with a single characteristic length, {{tmath|1= \sigma }}, and a minimum energy, <math>-\varepsilon</math> (with {{tmath|1= \varepsilon\ge 0 }}). Two of the many such functions that have been suggested are shown in the accompanying plot.{{sfn|ps=|Hirschfelder|Curtis|Bird|1964|pp=31-34}} 

A modern theory based on [[statistical mechanics]] produces the same result for <math>b=4N_\text{A}[(4\pi/3)(\sigma/2)^3]</math> obtained by van der Waals and his contemporaries. This result is valid for any pair potential for which the increase in <math>\varphi>0</math> is sufficiently rapid. This includes the hard sphere model for which the increase is infinitely rapid and the result is exact. Indeed, the Sutherland potential most accurately models van der Waals' conception of a molecule. It also includes potentials that do not represent hard sphere force interactions provided that the increase in <math>\varphi</math> for <math>r<\sigma</math> is fast enough, but then it is approximate; increasingly better the faster the increase. In that case <math>\sigma</math> is only an "effective diameter" of the molecule. This theory also produces <math>a=IN_\text{A}\varepsilon b</math> where <math>I</math> is a number that depends on the shape of the potential function, {{tmath|1= \varphi(r)/\varepsilon }}. However, this result is only valid when the potential is weak, namely, when the minimum potential energy is very much smaller than the thermal energy, <math>\varepsilon\ll kT</math>.{{sfn|ps=|Goodstein|1985|pp=250,263}}{{sfn|ps=|Tien|Lienhard|1979|pp=250,251}}

[[Ludwig Boltzmann]] wrote equations using <math>V/M</math> ([[specific volume]]) rather than <math>VN_\text{A}/N</math> (molar volume, used here);{{sfn|ps=|Boltzmann|1995|p=231}}{{sfn|ps=|van der Waals|1967|pp=168-172}} [[Josiah Willard Gibbs]] did as well, as do most engineers. Physicists use the property <math>V/N=1/\rho_N</math> (the reciprocal of [[number density]]), but there is no essential difference between equations written with any of these properties. Equations of state written using molar volume contain <math>R</math>, those using specific volume contain <math>R/\bar{m}</math> (the substance specific <math>\bar{m}=mN_\text{A}</math> is the [[molar mass]] with <math>m</math> the mass of a single particle), and those written with number density contain {{tmath|1= k }}.

Once the constants <math>a</math> and <math>b</math> are experimentally determined 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 (thermodynamics)|critical point]] (defined by pressure <math>p_\text{c}</math> and temperature <math>T_\text{c}</math> such that the substance cannot be liquefied either when <math>p>p_\text{c}</math> no matter how low the temperature, or when <math>T>T_\text{c}</math> no matter how high the pressure; <math>p_\text{c}, T_\text{c}</math> uniquely define {{tmath|1= v_\text{c} }}). These predictions are accurate for only a few substances. For most simple fluids they are only a valuable approximation. The equation also explains why [[superheating|superheated liquid]]s can exist above their boiling point and subcooled vapors can exist below their condensation point.

[[file:Vdw isobars.png|thumb|upright=1.6|Constant pressure curves of the van der Waals equation. On the green isobar ({{tmath|1= p=0.2\,p_\text{c} }}), the saturated liquid ({{tmath|1= v_f }}) and saturated vapor ({{tmath|1= v_g }}) points at temperature <math>T_\text{s}</math> are given by green circles.]]

=== Example ===

The graph on the right plots the intersection of the surface shown in Figures&nbsp;A and C and four planes of constant pressure. Each intersection produces a curve in the <math>T, v</math> plane corresponding to the value of the pressure chosen. These curves are [[Isobar (meteorology)|isobars]], since they represent all the points with the same pressure. 

On the red isobar, {{tmath|1= p=2p_\text{c} }}, the slope is positive over the entire range, <math>b\le v<\infty</math> (although the plot only shows a finite region). This describes a fluid as [[homogeneity and heterogeneity|homogeneous]] for all <math>T</math>—that is, it does not undergo a [[phase transition]] at any temperature—which is a characteristic of all supercritical isobars {{tmath|1= p>p_\text{c} }}.

The orange isobar, {{tmath|1= p=p_\text{c} }}, is the critical one that marks the boundary between homogeneity and heterogeneity. The critical point lies on this isobar. 

The green isobar, <math>p=0.2\,p_\text{c}</math>, has a region of negative slope. This region consists of states that are unstable and therefore never observed (for this reason this region is shown dotted gray). The green curve thus consists of two disconnected branches, indicated two phase states: a [[vapor]] on the right, and a denser liquid on the left.{{sfn|ps=|Truesdell|Bharatha|1977|pp=13–15}} For this pressure, at a temperature <math>T_\text{s}</math> (specified by mechanical, thermal, and material equilibrium), the boiling (saturated) liquid <math>v_f</math> and condensing (saturated) vapor <math>v_g</math> coexist, shown on the curve as the left and right green circles, respectively. The locus of these coexistent saturation points across all subcritical isobars forms the ''saturation curve'' on the <math>p, v, T</math> surface. In this situation, the denser liquid will separate and collect below the vapor due to gravity, and a [[meniscus (liquid)|meniscus]] will form between them. This [[homogeneity and heterogeneity|heterogeneous]] combination of coexisting liquid and vapor is the [[phase transition]]. Heating the liquid in this state increases the fraction of vapor in the mixture—its {{tmath|1= v< }}, an average of <math>v_f</math> and <math>v_g</math> weighted by this fraction, increases while <math>T_\text{s}</math> remains the same. This is shown as the horizontal dotted gray line, which represents not a solution of the equation but the observed behavior. The points above {{tmath|1= T_\text{s} }}, superheated liquid, and those below it, subcooled vapor, are [[metastable]]; a sufficiently strong disturbance causes them to transform to the stable alternative. These metastable regions are shown with green dashed lines. In summary, this isobar describes a fluid as a stable vapor for {{tmath|1= T>T_\text{s} }}, a stable liquid for {{tmath|1= T<T_\text{s} }}, and a mixture of liquid and vapor at {{tmath|1= T=T_\text{s} }}, that also supports metastable states of subcooled vapor and superheated liquid. This behavior is characteristic of all subcritical isobars {{tmath|1= 0<p<p_\text{c} }}, where <math>T_\text{s}</math> is a function of {{tmath|1= p }}.{{sfn|ps=|Epstein|1937|p=11}} 

The black isobar, {{tmath|1= p=0 }}, is the [[Limit (mathematics)|limit]] of positive pressures. None of its points represent stable solutions: they are either metastable (positive or zero slope) or unstable (negative slope). Interestingly, states of negative pressure (tension) exist. Their isobars lie below the black isobar, and form those parts of the surfaces seen in Figures&nbsp;A and C that lie below the zero-pressure plane. In this <math>T,v</math> plane they have a parabola-like shape, and, like the zero-pressure isobar, their states are all either metastable (positive or zero slope) or unstable (negative slope).

=== Surface plots ===
[[File:VdWsurface2.jpg|thumb|upright=1.6|Figure A: The surface generated by the van der Waals equation.]]

{{multiple image
| align             = right
| total_width       = 450
| image1            = VdW surface ideal3.png
| alt1              = 3D graph of the surface generated by the ideal gas equation
| caption1          = Figure B: The surface generated by the ideal gas equation.
| image2            = VdW surface3c.png
| alt2              = 3D graph of the surface generated by the van der Waals equation
| caption2          = Figure C: The folded surface generated by the van der Waals equation.
}}
Figure&nbsp;B shows the surface calculated from the ideal gas equation of state. This surface is universal, meaning it represents all ideal gases. Here, the surface is normalized so that the coordinate <math>( p_0, v_0, T_0 )</math> is at <math>( 1,1,1 )</math> in the 3-dimensional plot space (the black dot). This normalization makes it easier to compare this surface with the surface generated by the van der Waals equation in Figure&nbsp;C.{{sfn|ps=|Goodstein|1985|pp=303-304,316,452}}

Figures&nbsp;A and C show the [[Plot (graphics)#Surface plot|surface]] calculated from the van der Waals equation. Note that whereas the ideal gas surface is relatively uniform, the van der Waals surface has a distinctive "fold". This fold develops from a [[critical point (thermodynamics)|critical point]] defined by specific values of pressure, temperature, and molar volume. Because the surface is plotted using dimensionless variables (formed by the ratio of each property to its respective critical value), the critical point is located at the coordinates {{tmath|1= ( 1,1,1 ) }}. When drawn using these dimensionless axes, this surface is, like that of the ideal gas, also universal. Moreover, it represents all real substances to a remarkably high degree of accuracy. This principle of [[Theorem of corresponding states|corresponding states]], developed by van der Waals from his equation, has become one of the fundamental ideas in the thermodynamics of fluids.

The fold's boundary on the surface is marked, on each side of the critical point, by the ''[[spinodal]] curve'' (identified in Fig. A, and seen in Figs. A and C). This curve delimits an unstable region wherein no observable homogeneous states exist; elsewhere on the surface, states of liquid, vapor, and gas exist. The fold in the surface is what enables the equation to predict the phenomenon of [[liquid]]–[[vapor]] [[phase transition|phase change]]. This phenomenon is described by the ''saturation curve'' (or ''coexistence curve''): the [[locus (mathematics)|locus]] of saturated liquid and vapor states which, being in equilibrium with each other, can coexist. The saturation curve is ''not'' specified by the properties of the surface alone—it is substance-dependent. The saturated liquid curve and saturated vapor curve (both identified in Fig. A) together comprise the saturation curve. The inset in Figure A shows the ''mixture states'', which are a combination of the saturated liquid and vapor states that correspond to each end of the horizontal mixture line (that is, the points of intersection between the mixture line and its isotherm). However, these mixture states are not part of the surface generated by the van der Waals equation; they are not solutions of the equation.

== Relationship to the ideal gas law ==

The [[ideal gas law]] follows from the van der Waals equation whenever the molar volume <math>v</math> is sufficiently large (when {{tmath|1= v\gg b }}, so {{tmath|1= v-b\approx v }}), or correspondingly whenever the molar density, {{tmath|1= \rho=1/v }}, is sufficiently small (when {{tmath|1= v\gg (a/p)^{1/2} }}, so {{tmath|1= p+a/v^2\approx p }}).{{sfn|ps=|Epstein|1937|p=10}}

When <math>v</math> is large enough that both inequalities are satisfied, these two approximations reduce the van der Waals equation to {{tmath|1= p=RT/v }}; rearranging in terms of <math>V</math> and <math>N</math> gives {{tmath|1= pV=NkT }}, which is the ideal gas law.{{sf|ps=|Epstein|1937|p=10}} This is not surprising since the van der Waals equation was constructed from the ideal gas equation in order to obtain an equation valid beyond the limit of ideal gas behavior. 

What is truly remarkable is the extent to which van der Waals succeeded. Indeed, [[Paul Sophus Epstein|Epstein]] in his classic thermodynamics textbook began his discussion of the van der Waals equation by writing, "In spite of 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".{{sfn|ps=|Epstein|1937|p=10}} Also, in Volume 5 of his ''Lectures on Theoretical Physics'', [[Arnold Sommerfeld|Sommerfeld]], in addition to noting that "Boltzmann<ref>Boltzmann, L. ''Enzykl. der Mathem. Wiss.'', '''V''', (1), 550</ref> described van der Waals as the [[Isaac Newton|Newton]] of [[real gas]]es",{{sfn|ps=|Sommerfeld|1956|p=55}} 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.{{sfn|ps=|Sommerfeld|1956|p=66}}

== Utility of the equation ==

The van der Waals equation has been, and remains, useful because:{{sfn|ps=|Sommerfeld|1956|pp=55–68}}
* It yields simple analytic expressions for thermodynamic properties: internal energy, entropy, enthalpy, Helmholtz free energy, Gibbs free energy, and specific heat at constant pressure {{tmath|1= c_p }}.
* It yields an analytic expression of its coefficient of thermal expansion and its isothermal compressibility.
* It yields an analytic analysis of the [[Joule–Thomson effect|Joule–Thomson coefficient]] and associated inversion curve, which were instrumental in the development of the commercial liquefaction of gases.
* It shows that the [[Specific heat capacity|specific heat]] at constant volume <math>c_v</math> is a function of <math>T</math> only.
* It explains the existence of the [[critical point (thermodynamics)|critical point]] and the liquid–vapor [[phase transition]], including the observed [[Metastability|metastable]] states.
* It establishes the [[theorem of corresponding states]].
* It is an intermediate mathematical model, useful as a pedagogical tool when teaching physics, chemistry, and engineering.

In addition, its saturation curve has an analytic solution, which can depict the liquid metals (mercury and cesium) quantitatively, and describes most real fluids qualitatively.{{sfn|ps=|Lienhard|1986|pp=172-173}} As such, this solution can be regarded as one member of a family of equations of state (known as ''extended corresponding states'').{{sfn|ps=|Peck|1982}} Consequently, the equation plays an important role in the modern theory of phase transitions.{{sfn|ps=|Goodstein|1985|pp=443–452}}

== History ==

In 1857 [[Rudolf Clausius]] published ''The Nature of the Motion which We Call Heat''. In it he derived the relation <math>p=(N/V)m\overline{c^2}/3</math> for the pressure <math>p</math> in a gas, composed of particles in motion, with number density {{tmath|1= N/V }}, mass {{tmath|1= m }}, and mean square speed {{tmath|1= \overline{c^2} }}. He then noted that using the classical laws of Boyle and Charles, one could write <math>m\overline{c^2}/3=kT</math> with a constant of proportionality {{tmath|1= k }}. Hence temperature was proportional to the average kinetic energy of the particles.{{sfn|ps=|Weinberg|2021|pp=4–5}} 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 in accordance with 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]].{{sfn|ps=|Weinberg|2021|p=33}}{{sfn|ps=|Gibbs|1948|pp=vii–xii}}

[[File:VdWaalsLeiden2020.jpg|thumb|upright=1.6|Van der Waals equation on a wall in [[Leiden]]]]

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.{{sfn|ps=|van der Waals|1873}}{{sfn|ps=|van der Waals|1984|pp=121–240}} 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 force|van der Waals cohesive forces]].{{sfn|ps=|Boltzmann|1995|p=218}}

In 1869 Irish professor of chemistry [[Thomas Andrews (scientist)|Thomas Andrews]] at Queen's University Belfast, in a paper entitled ''On the Continuity of the Gaseous and Liquid States of Matter'',{{sfn|ps=|Andrews|1869}} displayed an experimentally obtained set of isotherms of carbonic acid, {{Chem2|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 [[Thomas Andrews (scientist)#Work|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.{{sfn|ps=|Klein|1974|p=31}}{{sfn|ps=|van der Waals|1967|pp=125,191–194}}

By 1877 sprays of both [[liquid oxygen]] and [[liquid nitrogen]] had been produced, and a new field of research, [[cryogenics|low temperature physics]], had been opened. The van der Waals equation played a part in all this, especially with respect to the liquefaction of hydrogen and helium which was finally achieved in 1908.{{sfn|ps=|Goodstein|1985|pp=450–451}} From measurements of <math>p_1, T_1</math> and <math>p_2, T_2</math> in two states with the same density, the van der Waals equation produces the values{{sfn|ps=|Boltzmann|1995|pp=232–233}} 
<math display=block>b=v-\frac{R(T_2-T_1)}{p_2-p_1}\qquad\text{and}\qquad a=v^2\frac{p_2T_1-p_1T_2}{T_2-T_1}.</math> 

Thus from two such measurements of pressure and temperature one could determine <math>a</math> and {{tmath|1= b }}, and from these values calculate the expected critical pressure, temperature, and molar volume. [[David Goodstein|Goodstein]] summarized this contribution of the van der Waals equation as follows:{{sfn|ps=|Goodstein|1985|p=452}}<blockquote>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.</blockquote>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".

As noted previously, modern-day studies of first-order phase changes make use of the van der Waals equation together with the Gibbs criterion, equal chemical potential of each phase, as a model of the phenomenon. This model has an analytic coexistence (saturation) curve expressed parametrically, <math>p_\text{s}=f_p(y), T_\text{s}=f_T(y)</math> (the parameter <math>y</math> is related to the entropy difference between the two phases), that was first obtained by Plank,{{sfn|ps=|van der Waals|1984|p=19}} was known to Gibbs and others, and was later derived by Lekner.{{sfn|ps=|Lekner|1982|pp=161-162}} A summary of Lekner's solution is presented in a subsequent section, and a more complete discussion in the [[Maxwell construction#application to the van der Waals equation|Maxwell construction]].

== Critical point and corresponding states ==

[[File:vdW isotherms+2log.png|thumb|upright=1.6|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&nbsp;1 shows four [[isotherms]] of the van der Waals equation (abbreviated as vdW) on a <math>p, v</math> (pressure, molar volume) plane. The essential character of these curves is that they come in three forms: 

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

The critical point can be analytically determined by evaluating the two partial derivatives of the vdW equation in (1) and equating them to zero. This produces the critical values <math>v_\text{c}=3b</math> and {{tmath|1= T_\text{c}=8a/(27Rb) }}; plugging these back into the vdW equation gives {{tmath|1= p_\text{c}=a/27b^2 }}.{{sfn|ps=|Sommerfeld|1956|pp=56–57}}

Above the critical temperature <math>T_\text{c}</math>, van der Waals isotherms satisfy the stability criterion that <math>\partial p/\partial v|_{T}<0</math>. 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 <math>v_{\rm min}\le v\le v_{\rm max}</math>, for which no observable states exist. The states for <math>v<v_{\rm min}</math> are liquid, and those for <math>v>v_{\rm max}</math> 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''.{{sfn|ps=|Lienhard|Shamsundar|Biney|1986|pp=297-298}} Their [[Locus (mathematics)|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 <math>\kappa_T=\infty</math>) as
<math display="block">T_{\rm sp}=2a\frac{(v-b)^2}{Rv^3}=T_\text{c}\frac{(3v_r-1)^2}{4v_r^3}\qquad p_{\rm sp}= \frac{a(v-2b)}{v^3}=p_\text{c}\frac{(3v_r-2)}{v_r^3}</math>

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

Using these critical values to define reduced (dimensionless) variables {{tmath|1= p_r=p/p_\text{c} }}, {{tmath|1= T_r=T/T_\text{c} }}, and <math>v_r=v/v_\text{c}</math> renders the vdW equation in the dimensionless form (used to construct Fig. 1):
<math display="block">p_r=\frac{8T_r}{3v_r-1}-\frac{3}{v_r^2}</math>

This dimensionless form is a similarity relation; it indicates that all vdW fluids at the same <math>T_r</math> will plot on the same curve. It expresses the ''law of corresponding states'' which Boltzmann described as follows:{{sfn|ps=|Boltzmann|1995|pp=239–240}}

<blockquote>
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.
</blockquote>

This "law" is just a special case of [[Dimensional analysis#Formulation|dimensional analysis]] in which an equation containing 6 dimensional quantities, {{tmath|1= p, v, T, a, b, R }}, and 3 independent dimensions, [p], [v], [T] (independent means that "none of the dimensions of these quantities can be represented as a product of powers of the dimensions of the remaining quantities",{{sfn|ps=|Barenblatt|1979|p=16}} and {{tmath|1= [R = pv/T] }}), must be expressible in terms of 6 − 3 = 3 dimensionless groups.{{sfn|ps=|Barenblatt|1979|pp=13–23}} Here <math>v^*=b</math> is a characteristic molar volume, <math>p^*=a/b^2</math> a characteristic pressure, and <math>T^*=a/(Rb)</math> a characteristic temperature, and the 3 dimensionless groups are {{tmath|1= p/p^*, v/v^*,T/T^* }}. According to dimensional analysis the equation must then have the form {{tmath|1= p/p^*=\Phi(v/v^*,T/T^*) }}, a general similarity relation. In his discussion of the vdW equation, Sommerfeld also mentioned this point.{{sfn|ps=|Sommerfeld|1956|p=57}} The reduced properties defined previously are {{tmath|1= p_r=27(p/p^*) }}, {{tmath|1= v_r=(1/3)(v/v^*) }}, and {{tmath|1= 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. Nevertheless, as Boltzmann observed, the van der Waals equation provides an essentially correct description.

The vdW equation produces the critical compressibility factor {{tmath|1= Z_\text{c}=p_\text{c}v_\text{c}/(RT_\text{c})=3/8=0.375 }}, while for most real fluids {{tmath|1= 0.23<Z_\text{c}< 0.31 }}.{{sfn|ps=|Johnston|2014|p=6}} 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|acentric]]) factor, {{tmath|1= \omega }}, and the liquid metals (mercury and cesium) are well approximated by it.{{sfn|ps=|Lienhard|1986|pp=172-173}}{{sfn|ps=|Dong|Lienhard|1986|pp=158-159}}

== Thermodynamic properties ==

The properties of molar internal energy <math>u</math> and entropy <math>S</math>—defined by the [[Thermodynamics equations#Laws of thermodynamics|first and second laws of thermodynamics]], hence all thermodynamic properties of a simple compressible substance—can be specified, up to a constant of integration, by two measurable functions: a mechanical equation of state {{tmath|1= p=p(v,T) }}, and a constant volume specific heat {{tmath|1= c_v(v,T) }}.

=== Internal energy and specific heat at constant volume ===

The internal energy is given by the energetic equation of state,{{sfn|ps=|Whitman|2023|p=155}}{{sfn|ps=|Moran|Shapiro|2000|p=574}}
<math display=block>u-C_u=\int\,c_v(v,T)\,dT+\int\,\left[T\frac{\partial p}{\partial T}-p(v,T)\right]\,dv=\int\,c_v(v,T)\,dT+\int\,\left[T^2\frac{\partial (p/T)}{\partial T} \right]\,dv
</math>
where <math>C_u</math> is an arbitrary constant of integration.

Now in order for <math>du(v,T)</math> to be an exact differential—namely that <math>u(v,T)</math> be continuous with continuous partial derivatives—its second mixed partial derivatives must also be equal, <math>\partial_{v}\partial_Tu=\partial_T\partial_v u</math>. Then with <math>c_v=\partial_Tu</math> this condition can be written as <math>\partial_vc(v,T)= \partial_T[T^2\partial_T(p/T)]</math>. Differentiating <math>p/T</math> for the vdW equation gives <math>T^2\partial_T(p/T)=a/v^2</math>, so <math>\partial_vc_v=0</math>. Consequently <math>c_v=c_v(T)</math> for a vdW fluid exactly as it is for an ideal gas. For simplicity, it is regarded as a constant here, <math>c_v=cR</math> for some constant number <math>c</math>. Then both integrals can be evaluated, resulting in

{{NumBlk||<math display=block>u-C_u=cRT-a/v</math>|{{EquationRef|2}}}}

This is the energetic equation of state for a perfect vdW fluid.

=== Enthalpy ===

The [[enthalpy]] of a system is given by <math>h=u+pv</math>. Substituting with {{EquationNote|2|Eq. 2}} and <math>pv={RTv \over (v-b)}-{a \over v}</math> (the vdW equation multiplied by <math>v</math>) gives  
<math display="block">h-C_u=RT \left( c+{v \over (v-b)} \right) -{2a \over v}.</math>  

This is the enthalpic equation of state for a perfect vdW fluid.

=== Entropy and free energies ===

The [[entropy]] <math>S</math> is given by the entropic equation of state:{{sfn|ps=|Whitman|2023|p=203}}{{sfn|ps=|Moran|Shapiro|2000|p=574}}
<math display="block">S-C_\text{s}=\int\,c_v(T)\,\frac{dT}{T}+\int\,\frac{\partial p}{\partial T}\,dv</math>

Using <math>c_v=cR</math> as before, and integrating the second term using <math>\partial_Tp=R/(v-b)</math> we obtain
<math display="block">S-C_\text{s}=R\ln[T^c(v-b)].</math>

This is the entropic equation of state for a perfect vdW fluid.

The [[Helmholtz free energy]] is <math>f=u-TS</math>, so combining the previous results gives
<math display="block">f=C_u+cT-{a \over v} -T\left(C_\text{s}+R\ln[T^c(v-b)]\right).</math>

This is the Helmholtz free energy for a perfect vdW fluid.

Likewise, the [[Gibbs free energy]] is <math>g=h-TS</math>, which evaluates to

{{NumBlk||<math display="block">g-C_u = \left(c+{v \over v-b}-C_\text{s}-\ln[T^c(v-b)]\right)RT-{2a \over v}</math>|{{EquationRef|3}}}}

This is the Gibbs free energy for a perfect vdW fluid.

=== Isothermal compressibility, thermal expansion and specific heat ===

The two first partial derivatives of the vdW equation are
<math display="block">\left.\frac{\partial p}{\partial T}\right|_v=\frac{R}{v-b}=\frac{\alpha}{\kappa_T}\quad\text{and}\quad\left.\frac{\partial p}{\partial v}\right|_T=-\frac{RT}{(v-b)^2}+
\frac{2a}{v^3}=-\frac{1}{v\kappa_T}</math>
where <math>\kappa_T=-v^{-1}\partial_pv</math> is the isothermal compressibility (a measure of the relative increase of volume from an increase of pressure, at constant temperature), and <math>\alpha=v^{-1}\partial_Tv_p</math> is the coefficient of thermal expansion (a measure of the relative increase of volume from an increase of temperature, at constant pressure). Therefore,{{sfn|ps=|Sommerfeld|1956|p=56}}{{sfn|ps=|Johnston|2014|p=11}}

{{NumBlk|:|<math>\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}.</math>|{{EquationRef|4}}}}

In the limit <math>v\rightarrow\infty</math>, the vdW equation becomes <math>p={RT \over v}</math>, and <math>\alpha={1 \over T}</math>, and <math>\kappa_T = {v \over RT} = p^{-1}</math>. Both these limits of <math>\alpha</math> and <math>\kappa_T</math> are the ideal gas values, which is consistent because, as noted earlier, a vdW fluid behaves like an ideal gas in this limit.

The specific heat at constant pressure <math>c_p</math> is defined as the partial derivative <math>c_p=\partial_Th|_p</math>. However, it is not independent of <math>c_v</math>, as they are related by the [[Julius von Mayer|Mayer]] equation, <math>c_p-c_v=-T(\partial_Tp)^2/\partial_vp=Tv\alpha^2/\kappa_T</math>.{{sfn|ps=|Whitman|2023|p=204}}{{sfn|ps=|Moran|Shapiro|2000|p=580}}{{sfn|ps=|Johnston|2014|p=3}} Then the two partials of the vdW equation can be used to express <math>c_p</math> as{{sfn|ps=|Johnston|2014|p=12}}
<math display="block">c_p(v,T)-c_v(T)=\frac{R^2Tv^3}{RTv^3-2a(v-b)^2}\ge R.</math>

Here in the limit <math>v\rightarrow\infty</math>, <math>c_p-c_v= R</math>, which is also the ideal gas result as expected;{{sfn|ps=|Johnston|2014|p=12}} however the limit <math>v\rightarrow b</math> gives the same result, which does not agree with experiments on liquids.

In this liquid limit we also find <math>\alpha=\kappa_T=0</math>, namely that the vdW liquid is incompressible. Moreover, since <math>\partial_Tp=-\partial_Tv/\partial_pv=\alpha/\kappa_T=\infty</math>, it is also mechanically incompressible, that is, <math>\kappa_T</math> approaches 0 faster than <math>\alpha</math> does.

Finally, <math>c_p</math>, <math>\alpha</math>, and <math>\kappa_T</math> are all infinite on the curve <math>T={2a(v-b)^2 \over Rv^3}=T_\text{c}{(3v_r-1)^2 \over 4v_r^3}</math>.{{sfn|ps=|Johnston|2014|p=12}} This curve, called the [[spinodal]] curve, is defined by <math>\kappa_T^{-1}=0</math>.

== Saturation ==

Although the gap in <math>v</math> delimited by the two spinodal points on an isotherm (e.g. <math>T_\text{r}=7/8</math> in Fig. 1) is the origin of the phase change, the change occurs as some intermediate value. This can be seen by considering that both saturated liquid and saturated vapor can coexist in equilibrium, at which they have the same pressure and temperature.{{sfn|ps=|Callen|1960|pp=37–44}} However, the minimum and maximum spinodal points are not at the same pressure. Therefore, at a temperature <math>T_\text{s}</math>, the phase change is characterized by the pressure <math>p_\text{s}</math>, which lies within the range of <math>p</math> set by the spinodal points (<math>p_\text{min}<p_\text{s}<p_\text{max}</math>), and by the molar volume of liquid <math>v_\text{f}</math> and vapor <math>v_\text{g}</math>, which lie outside the range of <math>v</math> set by the spinodal points (<math>v_\text{f}<v_\text{min}</math> and <math>v_\text{g}>v_\text{max}</math>).

Applying the vdW equation to the saturated liquid (fluid) and saturated vapor (gas) states gives:

{{NumBlk||<math display="block">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}.
</math>|{{EquationRef|5}}}}

These two equations contain four variables (<math>p_\text{s},T_\text{s},v_\text{f},v_\text{g}</math>), so a third equation is required in order to uniquely specify three of these variables in terms of the fourth. The following is a derivation of this third equation (the result is {{EquationNote|Eq. 6}}).

Now, the [[Latent heat#Specific latent heat|energy required to vaporize a mole]] at constant pressure is <math>Q_\text{vap}=h_\text{g}-h_\text{f}</math> (from the first law of thermodynamics) and at constant temperature is <math>Q_\text{vap}=T(S_\text{g}-S_\text{f})</math> (from the second law). Thus,
<math display="block">\begin{align}
Q_\text{vap} = h_\text{g}-h_\text{f} &= T(S_\text{g}-S_\text{f}) \\ 
h_\text{g}-TS_\text{g} &= h_\text{f}-TS_\text{f} \\ 
g_\text{g} &= g_\text{f}
\end{align}
</math>

That is, in this case, the Gibbs free energy in the saturated liquid state equals that in the saturated vapor state.{{sfn|ps=|Callen|1960|p=153}} The Gibbs free energy is one of the four thermodynamic potentials whose partial derivatives produce all other thermodynamics state properties;{{sfn|ps=|Callen|1960|pp=85–101}} its differential is <math>dg=\partial_pg\,dp+\partial_Tg\,dT =v\,dp-S\,dT</math>. Integrating this over an isotherm from <math>p_\text{s}, v_\text{f}</math> to <math>p_\text{s}, v_\text{g}</math>, noting that the pressure is the same at each endpoint, and setting the result to zero yields
<math display="block">g_\text{g}-g_\text{f}=\sum_{j=1}^3\int_{p_j}^{p_j+1}\,v_j\,dp=-\int_{v_\text{f}}^{v_\text{g}}\,p\,dv+p_\text{s}(v_\text{g}-v_\text{f})=0</math>

Here because <math>v</math> is a multivalued function, the <math>p</math> integral must be divided into 3 parts corresponding to the 3 real roots of the vdW equation in the form, <math>v(p,T)</math> (this can be visualized most easily by imagining Fig. 1 rotated <math>90^\circ</math>); the result is a special case of material equilibrium.{{sfn|ps=|Callen|1960|pp=146–156}} The last equality, which follows from integrating <math>v\,dp=d(pv)-p\,dv</math>, is the Maxwell [[Maxwell construction|equal area rule]], which requires that the upper area between the vdW curve and the horizontal through <math>p_\text{s}</math> be equal to the lower area.{{sfn|ps=|Maxwell|1875|pp=358-359}} This form means that the thermodynamic restriction that fixes <math>p_\text{s}</math> is specified by the equation of state itself, <math>p=p(v,T)</math>. Using the equation for the Gibbs free energy for the vdW equation ({{EquationNote|3|Eq. 3}}), the difference <math>g_\text{g}-g_\text{f}</math> can be evaluated as

{{NumBlk||<math display="block">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</math>|{{EquationRef|6}}}}

This is a third equation that along with the two vdW equations of {{EquationNote|5|Eqs. 5}} can be solved numerically. This has been done given a value for either <math>T_\text{s}</math> or <math>p_\text{s}</math>, and tabular results presented;{{sfn|ps=|Shamsundar|Lienhard|1983|pp=878,879}}{{sfn|ps=|Barrufet|Eubank|1989|p=170}} however, the equations also admit an analytic parametric solution obtained by Lekner.{{sfn|ps=|Lekner|1982|pp=161-162}} Details of this solution may be found in the [[Maxwell construction#Application to the van der Waals Equation|Maxwell construction]]; the results are:
<math display="block">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}</math>
<math display="block">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}}</math>
where
<math display="block">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</math>
and the parameter <math>0\le y<\infty</math> is given physically by <math>y=(S_\text{g}-S_\text{f})/(2R)</math>. The values of all other property discontinuities across the saturation curve also follow from this solution.{{sfn|ps=|Johnston|2014|pp=16-18}} 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. Various projections of this saturation curve are plotted in Figures&nbsp;1, 2a, and 2b.

[[File:Vdw5annotatedwith dropping.png|thumb|upright=1.6|Figure 2a: The coexistence (saturation) curve is the locus of the saturated liquid and vapor states, which smoothly meet at the critical point. The isotherm <math>T_\text{r}=0.9</math> is also displayed, showing the intersection.]]
[[File:Vdw stability-saturation.png|thumb|upright=1.6|Figure 2b: 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 <math>p_\text{r}, T_\text{r}</math> plane.]]
Referring back to Figure&nbsp;1, the isotherms for <math>T_\text{r}<1</math> are discontinuous. For example, the <math>T_\text{r}=7/8</math> (green) isotherm consists of two separate segments. The solid green lines represent stable states, and terminate at dots that represent the saturated liquid and vapor states that comprise the phase change. The dashed green lines represent metastable states ([[superheated liquid]] and subcooled vapor) that are created in the process of phase transition, have a short lifetime, and then devolve into their lower energy stable alternative.

At every point in the region between the two curves in Figure&nbsp;2b, there are two states: one stable and one metastable. The coexistence of these states can be seen in Figure&nbsp;1—for discontinuous isotherms, there are values of <math>p_\text{r}</math> 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".{{sfn|ps=|Boltzmann|1995|pp=248–250}} (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.{{sfn|ps=|Lienhard|Shamsundar|Biney|1986|p=297}}

In the same section, Boltzmann also addressed and explained the negative pressures which some liquid metastable states exhibit (for example, the blue isotherm <math>T_\text{r}=4/5</math> in Fig. 1). He concluded that such liquid states of tensile stresses were real, as did [[Chang-Lin Tien|Tien]] and [[John H. Lienhard|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."{{sfn|ps=|Tien|Lienhard|1979|p=254}}

Even though the phase change produces a mathematical discontinuity in the homogeneous fluid properties (for example <math>v</math>), there is no physical discontinuity.{{sfn|ps=|Callen|1960|pp=146–156}} As the liquid begins to vaporize, the fluid becomes a heterogeneous mixture of liquid and vapor whose molar volume varies continuously from <math>v_\text{f}</math> to <math>v_\text{g}</math> according to the equation of state <math display="inline">v=v_\text{f}+x(v_\text{g}-v_\text{f})</math> where <math display="inline">x=N_\text{g}/(N_\text{f}+N_\text{g})</math> and <math>0\le x\le 1</math> is the mole fraction of the vapor. This equation is called the [[lever rule]] and applies to other properties as well.{{sfn|ps=|Sommerfeld|1956|p=66}}{{sfn|ps=|Callen|1960|pp=146–156}} 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 ===

[[File:vapor pressure vs temperature1.png|thumb|upright=1.6|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, <math>p_\text{r}=p(v_\text{r},T_\text{r})</math>. However, as Boltzmann noted, such a simple representation could not correctly describe all substances. Indeed, the saturation analysis of this form produces <math>p_\text{rs}=p_\text{s}(T_\text{r})</math>; namely, that all substances have the same dimensionless coexistence curve, which is not true.{{sfn|ps=|van der Waals|1984|p=22}} To avoid this paradox, an extended principle of corresponding states has been suggested in which <math>p_\text{r}=p(v_\text{r},T_\text{r},\phi)</math> where <math>\phi</math> is a substance-dependent dimensionless parameter related to the only physical feature associated with an individual substance: its critical point.

[[File:Tr vs Pitzer factor.jpg|thumb|upright=1.6|Figure 4: A plot of the correlation including data from various substances.]]
One candidate for <math>\phi</math> is the critical compressibility factor <math>Z_\text{c}=p_\text{c}v_\text{c}/(RT_\text{c})</math>; however, because <math>v_\text{c}</math> is difficult to measure accurately, the [[Acentric factor|acentric]] factor developed by [[Kenneth Pitzer]],{{sfn|ps=|Pitzer|Lippman|Curl|Huggins|1955|p=3433}} <math>\omega=-\log_{10}[p_\text{r}(T_\text{r}=0.7)]-1</math>, is more useful. The saturation pressure in this situation is represented by a one-parameter family of curves: <math>p_\text{rs}=p_\text{s}(T_\text{r},\omega)</math>. Several investigators have produced correlations of saturation data for a number of substances; Dong and Lienhard give{{sfn|ps=|Dong|Lienhard|1986|pp=158-159}}
<math display="block">
\begin{align}
\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{align}</math>
which has an [[RMS error]] of <math>\pm 0.42</math> over the range <math>1\le T_\text{r}\le0.3</math>.

Figure&nbsp;3 is a plot of <math>p_\text{rs}</math> vs. <math>T_\text{r}</math> for various values of the Pitzer factor <math>\omega</math> as given by this equation. The vertical axis is logarithmic in order to show the behavior at pressures closer to zero, where differences among the various substances (indicated by varying values of <math>\omega</math>) are more pronounced.

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

== Joule–Thomson coefficient ==

The Joule–Thomson coefficient, <math>\mu_\text{JT}=\partial_pT|_h</math>, is of practical importance because the two end states of a throttling process (<math>h_2=h_1</math>) lie on a constant enthalpy curve. Although ideal gases, for which <math>h=h(T)</math>, 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.{{sfn|ps=|Sommerfeld|1956|pp=61–63}}

This coefficient can be found in terms of the previously derived <math>\alpha</math> and <math>c_p</math> as{{sfn|ps=|Sommerfeld|1956|pp=60-62}}
<math display="block">\mu_\text{JT}=\frac{v(\alpha T-1)}{c_p}.</math>

When <math>\mu_\text{JT}</math> is positive, the gas temperature decreases as it passes through a throttling process, and when it is negative, the temperature increases. Therefore, the condition <math>\mu_\text{JT}=0</math> defines a curve that separates the region of the <math>T,p</math> plane where <math>\mu_\text{JT}>0</math> from the region where <math>\mu_\text{JT}<0</math>. This curve is called the ''inversion curve'', and its equation is <math>\alpha T-1=0</math>. Evaluating this using the expression for <math>\alpha</math> derived in {{EquationNote|4|Eq. 4}} produces
<math display="block">\frac{2a(v-b)^2-RTv^2b}{RTv^3-2a(v-b)^2}=0\quad \Rightarrow \quad 2a(v-b)^2-RTv^2b=0</math>

Note that for <math>v\gg b</math> there will be cooling for <math>2a>RTb</math> (or, in terms of the critical temperature, <math>T<(27/4)\ T_\text{c}</math>). 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.''"{{sfn|ps=|Sommerfeld|1956|p=61}}

[[File:Vdw inversion2.png|thumb|upright=1.6|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).]]
In terms of <math>b/v</math>, the equation has a simple positive solution <math display="inline">b/v=1-\sqrt{RTb/(2a)} </math>, which for <math>b/v=0</math> produces <math>T=2a/(Rb)=27T_\text{c}/4</math>. Using this to eliminate <math>v</math> from the vdW equation then gives the inversion curve as 
<math display="block">\frac{p}{p^*}=-1+4\left(\frac{T}{2T^*}\right)^{\frac{1}{2}}-3\left(\frac{T}{2T^*}\right)</math> 
where, for simplicity, <math>a,b,R</math> have been replaced by <math>p^*,T^*</math>.

The maximum of this quadratic curve occurs with <math>z^2=T/(2T^*)</math>, for
<math display="block">D_zp/p^*=4-6(T/2T^*)^{\frac{1}{2}}=0</math>
which gives <math>(T/2T^*)^{\frac{1}{2}}_\text{max}=2/3</math>, or <math>T_\text{max}=8T^*/9</math>, and the corresponding <math>p_\text{max}=p^*/3</math>. By the quadratic formula, the zeros of the curve <math>3z^2- 4z+1=0</math> are <math>z=1/3</math> and <math>1</math> (<math>T/T^*=2/9=0.\overline{2}</math> and <math>2</math>). In terms of the dimensionless variables <math>T_\text{r}, p_\text{r}</math>, the zeros are at <math>T_\text{r}=3/4</math> and <math>27/4</math>, while the maximum is <math>p_\text{rmax}=9</math>, and occurs at <math>T_\text{rmax}=3</math>. A plot of the curve is shown in green in Figure&nbsp;5. Sommerfeld also displays this plot,{{sfn|ps=|Sommerfeld|1956|p=62|loc=Fig.8}} together with a curve drawn using experimental data from H<sub>2</sub>. The two curves agree qualitatively, but not quantitatively. For example the maximum on these two curves differ by about 40% in both magnitude and location.{{Excessive detail inline|date=January 2025}}

Figure&nbsp;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 ==

[[File:Vdw Z rho.png|thumb|upright=1.6|Figure 6: Isotherms plotted against the compressibility factor <math>Z</math>. The horizontal axis is <math>v^*\rho=\rho_r/3</math>. The isotherms, spinodal curve, and coexistence curve here are the same as in Fig. 1. Additionally plotted are the isotherm <math>T_r=27/8</math>, which has zero slope at the origin, and the isotherm <math>T_r=\infty</math>.]] [[File:Vdw Z p_r 1.png|thumb|upright=1.6|Figure 7: Generalized compressibility chart for a van der Waals gas.]] 
Real gases are characterized by their difference from ideal gases by writing <math>pv=ZRT</math>. Here <math>Z</math>, called the compressibility factor, is expressed either as <math>Z(p,T)</math> or <math>Z(\rho,T)</math>.{{Clarify|reason=Are p and rho just different notations of the same variable? If so, use just one (whichever is more common in the literature). If they are different, explain.|date=January 2025}} In either case, the limit as <math>p</math> or <math>\rho</math> approaches zero is 1, and <math>Z</math> takes the ideal gas value. In the second case <math>Z(\rho,T)=p(\rho,T)/\rho RT</math>,{{sfn|ps=|Van Wylen|Sonntag|1973|p=49}} so for a van der Waals fluid the compressibility factor is   

{{Numblk|:|<math>Z={1 \over 1-b\rho}-{a\rho \over RT}</math>|{{EquationRef|7}}}} 

or in terms of reduced variables 
<math display="block">Z=\frac{3}{3-\rho_r}-\frac{9\rho_r}{8T_r}</math>
where <math>0\le\rho_r=1/v_r\le 3</math>. At the critical point, <math>T_r=\rho_r=1</math> and <math>Z=Z_\text{c}=3/2-9/8=3/8</math>.

In the limit <math>\rho\rightarrow 0</math>, <math>Z=1</math>; the fluid behaves like an ideal gas, as mentioned before. The derivative <math>\partial_{\rho}Z|_T = b\left( (1-b\rho)^{-2}-{a \over bRT}\right)</math> is never negative when <math>{a \over bRT}=T^*/T\le 1</math>; that is, when <math>T/T^*\ge 1</math> (which corresponds to <math>T_r\ge27/8</math>). Alternatively, the initial slope is negative when <math>T/T^*<1</math>, is zero at <math>b\rho=1-(T/T^*)^{1/2}</math>, and is positive for larger <math>b\rho\le 1</math> (see Fig. 6). In this case, the value of <math>Z</math> passes through <math>1</math> when <math>b\rho_B=1-T_B/T^*</math>. Here <math>T_B=(27T_\text{c}/8)(1-b\rho_B)</math> is called the [[Boyle temperature]]. It ranges between <math>0 \le T_B \le 27T_\text{c}/8</math>, and denotes a point in <math>T,\rho</math> 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 <math>(\alpha, \kappa_T)
</math> nor <math>c_p
</math> reduce to their ideal gas values, other than where <math>b\rho_B\ll 1,\,T_B\sim 27T_\text{c}/8</math> the actual ideal gas region.{{sfn|ps=|Johnston|2014|p=10}}

Figure&nbsp;6 plots various isotherms of <math>Z(\rho,T_r)</math> vs <math>\rho_r</math>. 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&nbsp;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.{{sfn|ps=|Su|1946|p=803}}{{sfn|ps=|Moran|Shapiro|2000|p=113}} Note the [[Caustic (mathematics)|caustic]] generated by the crossing isotherms.

=== Virial expansion ===
[[Statistical mechanics]] suggests that the compressibility factor <math>Z</math> can be expressed by a power series, called a [[virial expansion]]:{{sfn|ps=|Tien|Lienhard|1979|pp=247–248}}
<math display="block">Z(\rho,T)=1+\sum_{k=2}^\infty\,B_k(T)(\rho)^{k-1}</math>
where the functions <math>B_k(T)</math> are the virial coefficients; the <math>k</math>th term represents a <math>k</math>-particle interaction.

Expanding the term <math>(1-b\rho)^{-1}</math> in the definition of <math>Z</math> ({{EquationNote|7|Eq. 7}}) into an infinite series, convergent for <math>b\rho<1</math>, produces 
<math display="block">Z(\rho,T)=1+\left(1-{a \over bRT}\right)b\rho+\sum_{k=3}^\infty\ B_k(T)(b\rho)^{k-1}.
</math> 

The corresponding expression for <math>Z(\rho_r, T_r)</math> when <math>\rho_r<3</math> is
<math display="block">Z(\rho_r,T_r)=1+\left(1-{27 \over {8T_r}}\right)^{-1}\left({\rho_r \over 3}\right)+\sum_{k=2}^\infty\left({\rho_r \over 3}\right)^{k-1}.</math>

These are the virial expansions, one dimensional and one dimensionless, for the van der Waals fluid. The second virial coefficient is the slope of <math>Z(\rho_r,T_r)</math> at <math>\rho_r=0</math>. Notice that it can be positive when <math>T_r>27/8</math> or negative when <math>T_r<27/8</math>, which agrees with the result found previously by differentiation.

For molecules modeled as non-attracting hard spheres, <math>a=0</math>, and the vdW virial expansion becomes 
<math display="block">Z(\rho)=(1-b\rho)^{-1}=1+\sum_{k=2}^\infty (b\rho)^{k-1},</math> 
which illustrates the effect of the excluded volume alone. It was recognized early on that this was in error beginning with the term <math>(b\rho)^2</math>. Boltzmann calculated its correct value as <math display="inline">\frac{5}{8} (b\rho)^2</math>, and used the result to propose an enhanced version of the vdW equation:
<math display="block">\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).</math>

On expanding <math>(v-b/3)^{-1}</math>, this produced the correct coefficients through <math>(b/v)^2</math> and also gave infinite pressure at <math>b/3</math>, which is approximately the [[Close-packing of equal spheres|close-packing]] distance for hard spheres.{{sfn|ps=|Boltzmann|1995|pp=353-356}} 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.{{sfn|ps=|van der Waals|1984|pp=20-22}}

== Mixtures ==

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.{{sfn|ps=|van der Waals|1967|pp=243-282}} His essential idea was that in a binary mixture of vdW fluids described by the equations
<math display=block>
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}</math>
the mixture is also a vdW fluid given by
<math display=block>p=\frac{RT}{v-b_x}-\frac{a_x}{v^2}</math>
where
<math display="block">a_x=a_{11}x_1^2+2a_{12}x_1x_2+a_{22}x_2^2\quad\text{and}\quad b_x=b_{11}x_1^2+2b_{12}x_1x_2+b_{22}x_2^2.</math>

Here <math>x_1=N_1/N</math> and <math>x_2=N_2/N</math>, with <math>N=N_1+N_2</math> (so that <math>x_1+x_2=1</math>), are the mole fractions of the two fluid substances. Adding the equations for the two fluids shows that <math>p\ne p_1+p_2</math>, although for <math>v</math> sufficiently large <math>p\approx p_1+p_2</math> with equality holding in the ideal gas limit. The quadratic forms for <math>a_x</math> and <math>b_x</math> are a consequence of the forces between molecules. This was first shown by Lorentz,{{sfn|ps=|Lorentz|1881|pp=127,134,600}} and was credited to him by van der Waals. The quantities <math>a_{11},\,a_{22}</math> and <math>b_{11},\,b_{22}</math> in these expressions characterize collisions between two molecules of the same fluid component, while <math>a_{12}=a_{21}</math> and <math>b_{12}=b_{21}</math> 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''.{{sfn|ps=|van der Waals|1984|p=68}}

Assuming that <math>b_{12}</math> is the arithmetic mean of <math>b_{11}</math> and <math>b_{22}</math>, <math>b_{12}=(b_{11}+b_{22})/2</math>, substituting into the quadratic form and noting that <math>x_1+x_2=1</math> produces
<math display=block>b=b_{11}x_1+b_{22}x_2</math>

Van der Waals wrote this relation, but did not make use of it initially.{{sfn|ps=|van der Waals|1967|p=244}} However, it has been used frequently in subsequent studies, and its use is said to produce good agreement with experimental results at high pressure.{{sfn|ps=|Redlich|Kwong|1949}}

=== Common tangent construction ===

In this article, van der Waals used the Helmholtz potential minimum principle to establish the conditions of stability. This principle states that in a system in [[Diathermal wall|diathermal]] contact with a heat reservoir <math>T=T_R</math>, <math>DF=0</math>, and <math>D^2F>0</math>, namely at equilibrium, the Helmholtz potential is a minimum.{{sfn|ps=|Callen|1960|p=105}} Since, like <math>g(p,T)</math>, the molar Helmholtz function <math>f(v,T)</math> is also a potential function whose differential is
<math display="block">df=\partial_vf|_T\,dv+\partial_Tf|_v\,dT=-p\,dv-S\,dT, </math> 
this minimum principle leads to the stability condition <math>\partial^2f/\partial v^2|_T=-\partial p/\partial v|_T> 0</math>. This condition means that the function, <math>f</math>, is [[convex function|convex]] at all stable states of the system. Moreover, for those states the previous stability condition for the pressure is necessarily satisfied as well.  

==== Single fluid ====
For a single substance, the definition of the molar Gibbs free energy can be written in the form <math>f=g-pv</math>. Thus when <math>p</math> and <math>g</math> are constant along with temperature, the function <math>f(T_R, v)</math> represents a straight line with slope <math>-p</math>, and intercept <math>g</math>. Since the curve <math>f(T_R, v)</math> has positive curvature everywhere when <math>T_R\ge T_\text{c}</math>, the curve and the straight line will have a single tangent. However, for a subcritical <math>T_R,\, f(T_R,v)</math> is not everywhere convex. With <math>p=p_\text{s}(T_R)</math> and a suitable value of <math>g</math>, the line will be tangent to <math>f(T_R,v)</math> at the molar volume of each coexisting phase: saturated liquid <math>v_f(T_R)</math> and saturated vapor <math>v_g(T_R)</math>; there will be a double tangent. Furthermore, each of these points is characterized by the same values of <math>g</math>, <math>p</math>, and <math>T_R.</math> These are the same three specifications for coexistence that were used [[#Saturation|previously]]. 

[[file:Vdw helmholtz1.png|thumb|upright=1.6|Figure 8: The tangent line (black) to the curve <math>f_\text{r}(0.875,v_\text{r})</math> (green) at the two points <math>v_{rf}=0.576</math> and <math>v_{rg}=2.71</math>. The line's slope, given by <math>\partial_{v_\text{r}}f_\text{r}=-p_\text{c}v_\text{c}/(RT_\text{c})p_{rs}</math>, is <math>-0.215 </math> corresponding to <math>p_{rs}=0.5730</math>. This is the same data as isotherm <math>T_\text{r}=7/8</math> in Fig. 1. The intercept on the line is <math>g</math>, but its numerical value is arbitrary due to a constant of integration.]]
Figure&nbsp;8 depicts an evaluation of <math>f(T_R,v)</math> as a green curve, with <math>v_f</math> and <math>v_g</math> marked by the left and right green circles, respectively. The region on the green curve for <math>v\le v_f</math> corresponds to the liquid state. As <math>v</math> increases past <math>v_f</math>, the curvature of <math>f</math> (proportional to <math>\partial_v\partial_vf=-\partial_vp</math>) continually decreases. The [[inflection point]], characterized by zero curvature, is a spinodal point; between <math>v_f</math> and this point is the metastable superheated liquid. For further increases in <math>v</math> 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 <math>v</math> the curvature increases to a maximum at <math>v_g</math>, where the slope is <math>p_\text{s}</math>; the region between this point and the second spinodal point is the metastable subcooled vapor. Finally, the region <math>v\ge v_g</math> is the vapor. In this region the curvature continually decreases until it is zero at infinitely large <math>v</math>. The double tangent line (solid black) that runs between <math>v_f</math> and <math>v_g</math> represents states that are stable but heterogeneous, not homogeneous solutions of the vdW equation.{{sfn|ps=|van der Waals|1967|pp=245-247}} The states above this line (with larger Helmholtz free energy) are either metastable or unstable.{{sfn|ps=|van der Waals|1967|pp=245-247}} The combined solid green-black curve in Figure&nbsp;8 is the convex envelope of <math>f(T_R,v)</math>, which is defined as the largest convex curve that is less than or equal to the function.{{sfn|ps=|Lebowitz|1974|p=52}}

For a vdW fluid, the molar Helmholtz potential is
<math display="block">\begin{align}
f_\text{r}&=u_\text{r}-T_\text{r}S_\text{r} \\ 
&=\mathrm{C}_u+T_\text{r}\left(c-\mathrm{C}_\text{s}-\ln[T_\text{r}^c(3v_\text{r}-1)]\right)-{9 \over {8v_\text{r}}}
\end{align}</math>
where <math>f_\text{r}={f \over {RT_\text{c}}}</math>. Its derivative is
<math display="block">\partial_{v_\text{r}}f_\text{r}=-{{3T_\text{r}} \over {3v_\text{r}-1}}+{9 \over {(8v_\text{r})^2}}=-p_\text{r} </math> 
which is the vdW equation, as expected. A plot of this function <math>f_\text{r}</math>, whose slope at each point is specified by the vdW equation, for the subcritical isotherm <math>T_\text{r}=7/8</math> is shown in Figure&nbsp;8 along with the line tangent to it at its two coexisting saturation points. The data illustrated in Figure&nbsp;8 is exactly the same as that shown in Figure&nbsp;1 for this isotherm. This double tangent construction thus provides a graphical alternative to the Maxwell construction to establish the saturated liquid and vapor points on an isotherm.

==== Binary fluid ====
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, <math>f(T_R,v,x)</math>, where <math>x</math> is a composition variable (for example <math>x=x_2</math> so <math>x_1=1-x</math>). In this case, there are three stability conditions:
<math display="block">\frac{\partial^2f}{\partial v^2}>0\qquad\frac{\partial^2f}{\partial x^2}>0\qquad\frac{\partial^2f}{\partial v^2}\frac{\partial^2f}{\partial x^2}-\left(\frac{\partial^2f}{\partial x\partial v}\right)^2>0</math>
and the Helmholtz potential is a surface (of physical interest in the region <math>0\le x\le 1</math>). The first two stability conditions show that the curvature in each of the directions <math>v</math> and <math>x</math> are both non-negative for stable states, while the third condition indicates that stable states correspond to elliptic points on this surface.{{sfn|ps=|Kreyszig|1959|pp=124-128}} Moreover, its limit 
<math display="block">\frac{\partial^2f}{\partial v^2}\frac{\partial^2f}{\partial x^2}
-\frac{\partial^2f}{\partial x\partial v}=0</math>
specifies the spinodal curves on the surface.

For a binary mixture, the [[Homogeneous function#Euler's theorem|Euler equation]]{{sfn|ps=|Callen|1960|pp=47-48}} can be written in the form 
<math display="block">\begin{align}
f&=-pv+\mu_1x_1+\mu_2x_2 \\ 
&=-pv+(\mu_2-\mu_1)x+\mu_1
\end{align}</math>
where <math>\mu_j=\partial_{x_j}f</math> are the molar [[chemical potential]]s of each substance, <math>j=1,2</math>. For constant values of <math>p</math>, <math>\mu_1</math>, and <math>\mu_2</math>, this equation is a plane with slopes <math>-p</math> in the <math>v</math> direction, <math>\mu_2-\mu_1</math> in the <math>x</math> direction, and intercept <math>\mu_1</math>. 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 <math>T</math>, <math>p</math>, <math>\mu_2-\mu_1</math>, and <math>\mu_1</math>; the last two are equivalent to having the same <math>\mu_1</math> and <math>\mu_2</math> individually, which are just the Gibbs conditions for material equilibrium in this situation. The two methods of producing the coexistence surface are equivalent. 

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 [[Straw plaiting#etyomology|plait]] or fold) in the <math>x</math> direction as well as the one in the <math>v</math> direction. Therefore, there can be two liquid phases that can be either [[Miscibility|miscible]], or wholly or partially immiscible, as well as a vapor phase.{{sfn|ps=|van der Waals|1984|pp=23-27}}{{sfn|ps=|van der Waals|1967|pp=253-258}} 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{{sfn|ps=|DeBoer|1974|pp=7-16}}—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.{{sfn|ps=|van der Waals|1984|pp=23-27,64-66}} The results obtained were, in Rowlinson's words,
<blockquote>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.{{sfn|ps=|van der Waals|1984|p=66}}</blockquote>

=== Mixing rules ===

In order to obtain these numerical results, the values of the constants of the individual component fluids <math>a_{11}, a_{22},b_{11}, b_{22}</math> must be known. In addition, the effect of collisions between molecules of the different components, given by <math>a_{12}</math> and <math>b_{12}</math>, 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:{{sfn|ps=|Hirschfelder|Curtis|Bird|1964|pp=252-253}} 
<math display="block">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.</math>

These relations correspond to the empirical combining rules for the intermolecular force constants,
<math display="block">\epsilon_{12}=(\epsilon_{11}\epsilon_{22})^{1/2}\qquad\text{and}\qquad\sigma_{12}=(\sigma_{11}+\sigma_{22})/2,</math>
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.{{sfn|ps=|Hirschfelder|Curtis|Bird|1964|pp=168-169}} Using these empirical combining rules to generalize for <math>n</math> fluid components, the quadradic mixing rules for the material constants are:{{sfn|ps=|Redlich|Kwong|1949}}
<math display="block">a_x = \sum_{i=1}^{n} \sum_{j=1}^{n} (a_{ii}a_{jj})^{1/2}x_ix_j\quad\text{which can be written as}\quad  a_x=\Big(\sum_{i=1}^{n}a_{ii}^{1/2}x_i \Big)^2</math>
<math display="block">b_x=(1/8)\sum_{i=1}^{n}\sum_{j=1}^{n}(b_{ii}^{1/3} +b_{jj}^{1/3})^3 x_ix_j</math>

These expressions come into use when mixing gases in proportion, such as when producing tanks of air for diving{{sfn|ps=|Hewitt}} and managing the behavior of fluid mixtures in engineering applications. However, more sophisticated mixing rules are often necessary, in order to obtain satisfactory agreement with reality over the wide variety of mixtures encountered in practice.{{sfn|ps=|Valderrama|2003|pp=1308-1312}}{{sfn|ps=|Kontogeorgis|Privat|Jaubert|2019|pp=4626-4633}}

Another method of specifying the vdW constants, pioneered by W.B. Kay and known as ''Kay's rule'',{{sfn|ps=|Niemeyer}} specifies the effective critical temperature and pressure of the fluid mixture by
<math display="block">T_{\text{c}x}=\sum_{i=1}^nT_{\text{c}i}x_i\qquad\text{and}\qquad p_{\text{c}x}=\sum_{i=1}^n\,p_{\text{c}i}x_i.</math>

In terms of these quantities, the vdW mixture constants are 
<math display="block">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}}</math>
which Kay used as the basis for calculations of the thermodynamic properties of mixtures.{{sfn|ps=|van der Waals|1984|p=69}} Kay's idea was adopted by T. W. Leland, who applied it to the molecular parameters <math>\epsilon, \sigma</math>, which are related to <math>a, b</math> through <math>T_\text{c}, p_\text{c}</math> by <math>a\propto\epsilon\sigma^3</math> and <math> b\propto\sigma^3</math>. Using these together with the quadratic mixing rules for <math>a, b</math> produces
<math display="block">\sigma_x^3=\sum_{i=i}^n\sum_{j=1}^n\, \sigma_{ij}^3x_ix_j\qquad\text{and}\qquad \epsilon_x=\left[\sum_{i=1}^n\sum_{j=1}^n\epsilon_{ij}\sigma_{ij}^3x_ix_j\right]\left[\sum_{i=i}^n\sum_{j=1}^n\, \sigma_{ij}^3x_ix_j\right]^{-1}</math> 
which is the van der Waals approximation expressed in terms of the intermolecular constants.{{sfn|ps=|Leland|Rowlinson|Sather|Watson|1968|p=1447}}{{sfn|ps=|van der Waals|1984|pp=69-70}} This approximation, when compared with computer simulations for mixtures, are in good agreement over the range <math>1/2< (\sigma_{11}/\sigma_{22})^3<2</math>, namely for molecules of similar diameters. In fact, 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".{{sfn|ps=|van der Waals|1984|p=70}}

== Mathematical and empirical validity ==

Since van der Waals presented his thesis, "[m]any derivations, pseudo-derivations, and plausibility arguments have been given" for it.{{sfn|ps=|Goodstein|1985|p=443}} However, no mathematically rigorous derivation of the equation over its entire range of molar volume that begins from a [[statistical mechanics|statistical mechanical]] principle exists. Indeed, such a proof is not possible, even for hard spheres.{{sfn|ps=|Korteweg|1891b|p=277}}{{sfn|ps=|Tonks|1936|pp=962-963}}{{sfn|ps=|Kac|Uhlenbeck|Hemmer|1963|p=224}} Goodstein writes, "Obviously the value of the van der Waals equation rests principally on its empirical behavior rather than its theoretical foundation."{{sfn|ps=|Goodstein|1985|p=446}}

=== History of derivations ===

The classical [[Canonical ensemble|canonical partition function]], <math>Q(V,T,N)</math>, of statistical mechanics for a three dimensional <math>N</math> particle macroscopic system is

{{Numblk|:|<math display="block">Q_N=\Lambda^{-3N}N!^{-1}{\cal Z}_N\qquad\text{with}\qquad{\cal Z}_N=\int_{V^{\bf N}}\,\exp[-\Phi({\bf r^N})/kT]\,d{\bf r^N} </math>|{{EquationRef|8}}}}

where <math>Q_N(V,T)\equiv Q(V,T,N)</math>, <math>\Lambda=h(2\pi mkT)^{-1/2}</math> is the [[DeBroglie wavelength|de Broglie wavelength]] (alternatively <math>\Lambda^{-3}\equiv n_Q</math> is the quantum concentration), <math>{\cal Z}_N(V,T)\equiv{\cal Z}(V,T,N)</math> is the <math>N</math> particle configuration integral, and <math>\Phi</math> is the intermolecular potential energy, which is a function of the <math>N</math> particle position vectors <math>{\bf r^N}={\bf r}_1,{\bf r}_2\ldots{\bf r}_N</math>. Lastly <math>d{\bf r^N}=d{\bf r}_1d{\bf r}_2\ldots d{\bf r}_N</math> is the volume element of <math>V^{\bf N}</math>, which is a <math>3N</math>-dimensional space.{{sfn|ps=|Goodstein|1985|pp=51,61-68}}{{sfn|ps=|Tien|Lienhard|1979|pp=241-252}}{{sfn|ps=|Hirschfelder|Curtis|Bird|1964|pp=132-141}}{{sfn|ps=|Hill|1986|pp=112-119}}  

The connection of <math>Q_N</math> with thermodynamics is made through the Helmholtz free energy, <math>F=-kT\ln Q_N</math>, from which all other properties can be found; in particular <math> p=-\partial_VF|_T=kT\partial_V\ln\cal{Z}_N</math>. For point particles that have no force interactions (<math>\Phi=0</math>), all <math>3N</math> integrals of <math>{\cal Z}</math> can be evaluated producing <math>Q_N=(n_QV)^N/N!</math>. In the thermodynamic limit, <math>N\rightarrow\infty,\,V\rightarrow\infty</math> with <math>V/N</math> finite, the Helmholtz free energy per particle (or per mole, or per unit mass) is finite; for example, per mole it is <math>N_\text{A}F/N=f=-RT[\ln(n_Qv/N_\text{A})+1]</math>. The thermodynamic state equations in this case are those of a monatomic ideal gas, specifically
<math>-\partial_vf|_T=p=RT/v.</math>{{sfn|ps=|Hirschfelder|Curtis|Bird|1964|p=133}} 

Early derivations of the vdW equation were criticized mainly on two grounds.{{sfn|ps=|Kac|Uhlenbeck|Hemmer|1963|p=223}} First, a rigorous derivation from the partition function should produce an equation that does not include unstable states for which, <math>\partial_vp|_T>0</math>. Second, the constant <math>b=4N_{\text{A}}V_0</math> in the vdW equation (here <math>V_0</math> is the volume of a single molecule) gives the maximum possible number of molecules as <math>4N_{\text{max}}V_0=V</math>, or a close packing density of 1/4=0.25, whereas the known [[close-packing of equal spheres|close-packing density]] of spheres is <math>\pi/(3\sqrt{2})\approx 0.740</math>.{{sfn|ps=|Korteweg|1891b|p=277}} Thus a single value of <math>b</math> cannot describe both gas and liquid states.

The second criticism is an indication that the vdW equation cannot be valid over the entire range of molar volume. Van der Waals was well aware of this problem; he devoted about 30% of his Nobel lecture to it, and also said that it is{{sfn|ps=|van der Waals|1967|p=256}}
<blockquote>... the weak point in the study of the equation of state. I still wonder whether there is a better way. In fact this question continually obsesses me, I can never free myself from it, it is with me even in my dreams.
</blockquote> 
In 1949 the first criticism was confirmed by [[Leon Van Hove|van Hove]] when he showed that in the thermodynamic limit, hard spheres with finite-range attractive forces have a finite Helmholtz free energy per particle. Furthermore, this free energy is a continuously decreasing function of the volume per particle (see Fig. 8 where <math>f, v</math> are molar quantities). In addition, its derivative exists and defines the pressure, which is a non-increasing function of the volume per particle.{{sfn|ps=|van Hove|1949|p=951}} Since the vdW equation has states for which the pressure increases with increasing volume per particle, this proof means it cannot be derived from the partition function, without an additional constraint that precludes those states.

In 1891 [[Diederik Korteweg|Korteweg]] used kinetic theory ideas to show{{sfn|ps=|Korteweg|1891a|p=153}} that a system of <math>N</math> hard rods of length <math>\delta</math>, constrained to move along a straight line of length <math>L>N\delta</math> and exerting only direct contact forces on one another, satisfy a vdW equation with <math>a=0</math>; [[John William Strutt, 3rd Baron Rayleigh|Rayleigh]] also knew this.{{sfn|ps=|Strutt|1891|p=81|loc=footnote 1}} [[Lewi Tonks|Tonks]], by evaluating the configuration integral,{{sfn|ps=|Tonks|1936|p=959}} later showed that the force exerted on a wall by this system is given by <math>F_W=NkT/(L-N\delta) =kT/(l-\delta)</math> with <math>l=L/N</math>. This can be put in a more recognizable, molar, form by dividing by the rod cross sectional area <math>A</math>, and defining <math>N_\text{A}AL/N=v>b=N_\text{A}A\delta</math>. This produces <math>p=RT/(v-b)</math>; there is no condensation, as <math>\partial_vp|_T<0</math> for all <math>b\le v<\infty</math>. This result is obtained because in one dimension, particles cannot pass by one another as they can in higher dimensions; their mass center coordinates <math>x_i,\, i=1, 2, \dots, N</math> satisfy the relations <math>\delta/2\le x_1\le x_2-\delta,\cdots,\le x_N\le L-\delta/2</math>. As a result, the configuration integral is <math>{\cal Z}_N=(L-N\delta)^N</math>.{{sfn|ps=|Kac|1958|p=224}} 

In 1959 this one-dimensional gas model was extended by [[Marc Kac|Kac]] to include particle pair interactions through an attractive potential, 
<math>\varphi(x)=-\varepsilon\exp(-x/\ell),\, x\ge\delta</math>. This specific form allowed evaluation of the [[Statistical Mechanics#Three thermodynamic ensembles|grand partition function]],

{{Numblk|:|<math>Q_G(T,V,\mu)=\sum_{{\cal N}\ge 0}\,Q_{\cal N}(T,V)e^{\mu{\cal N}/kT},</math>|{{EquationRef|9}}}}

in the thermodynamic limit, in terms of the eigenfunctions and eigenvalues of a homogeneous integral equation.{{sfn|ps=|Kac|1958}} Although an explicit equation of state was not obtained, it was proved that the pressure was a strictly decreasing function of the volume per particle, hence condensation did not occur. 

[[file:vdw isotherms-maxwell.png|thumb|upright=1.6|Figure 9: Isotherms of the vdW equation + the Maxwell construction. The subcritical isotherm (<math>T_r=7/8</math>) is a composite curve: the black section is composed of heterogeneous states (liquid and vapor), while the green sections contain homogeneous states.]]
Four years later, in 1963, Kac together with [[George Uhlenbeck|Uhlenbeck]] and [[Per Christian Hemmer|Hemmer]] modified the pair potential of Kac's previous work as <math>\varphi(x)=-\varepsilon(\delta/\ell)\exp(-x/\ell),\,x\ge\delta</math>, so that
<math display=block>a_N=-\int_0^\infty\,\varphi(x)\,dx =\varepsilon\delta</math> 
was independent of <math>\ell</math>.{{sfn|ps=|Kac|Uhlenbeck|Hemmer|1963|pp=216-217}} They found that a second limiting process they called the ''van der Waals limit'', <math>\ell\rightarrow\infty</math> (in which the pair potential becomes both infinitely long range and infinitely weak) and performed after the thermodynamic limit, produced the one-dimensional vdW equation (here rendered in molar form)
<math display="block">p=\frac{RT}{v-b}-\frac{a}{v^2} \quad\text{where}\quad 
\begin{cases} a=N_\text{A}^2Aa_N=N_\text{A}\varepsilon b \\ b=N_\text{A}V_0=N_\text{A}A\delta 
\end{cases} </math>
as well as the Gibbs criterion, <math>\mu_f=\mu_g</math> (equivalently the [[Maxwell construction]]). As a result, all isotherms satisfy the condition <math>\partial_vp|_T\le 0</math> as shown in Figure&nbsp;9, and hence the first criticism of the vdW equation is not as serious as originally thought.{{sfn|ps=|Kac|Uhlenbeck|Hemmer|1963|p=224}}

Then, in 1966, [[Joel Lebowitz|Lebowitz]] and [[Oliver Penrose|Penrose]] generalized what they called the ''Kac potential'', <math>\varphi</math>, to apply to a nonspecific function of <math>\nu</math> dimensions:
<math display="block">\varphi=-\varepsilon(\sigma/\ell)^\nu\bar{\varphi}(x/\ell).</math>

For <math>\nu=1</math> and <math>\bar{\varphi}=\exp(-x/\ell)</math> this reduces to the specific one-dimensional function considered by Kac, et al., and for <math>\nu=3</math> it is an arbitrary function (although subject to specific requirements) in physical three-dimensional space. In fact, the function <math>\bar{\varphi}(x/\ell)</math> must be bounded, non-negative, and one whose integral
<math display="block">a_N=-\frac{1}{2}\int_{V^\nu}\,\varphi\left(\frac{r}{\ell}\right)\,d{\bf r}^\nu=\frac{\varepsilon\sigma^\nu}{2}\int_{V^\nu}\,\bar{\varphi}(x)\,d{\bf x}^\nu</math> 
is finite, independent of <math>\ell</math>.{{sfn|ps=|Lebowitz|Penrose|1966|p=98}}{{sfn|ps=|Lebowitz|1974|pp=50-52}} By obtaining upper and lower bounds on <math>{\cal Z}(V,T,N,\ell)</math> and hence on <math>F</math>, taking the thermodynamic limit (<math>\lim N,V\rightarrow\infty</math> with <math>N/V</math> finite) to obtain upper and lower bounds on the function <math>F(V,T,N/V,\ell)/V</math>, then subsequently taking the van der Waals limit, they found that the two bounds coalesced and thereby produced a unique limit (here written in terms of the free energy per mole and the molar volume):
<math display="block"> f(v,T)= \mathrm{CE}\{f^0(v,T)-a/v\}.</math>

The abbreviation <math>\mathrm{CE}</math> stands for "convex envelope"; this is a function which is the largest [[convex function]] that is less than or equal to the original function. The function <math>f^0(v,T)</math> is the limit function when <math>\varphi=0</math>; also here <math>a=N_A^2a_N</math>. This result is illustrated in Figure&nbsp;8 by the solid green curves and black line, which is the convex envelope of <math>f(T_R,v)</math>. 

The corresponding limit for the pressure is a generalized form of the vdW equation
<math display=block>
p(v,T)=p^0(v,T)-a/v^2</math>
together with the Gibbs criterion, <math>\mu_f=\mu_g</math> (equivalently the [[Maxwell construction]]). Here <math>p^0=-\partial_vf^0|_T</math> is the pressure when attractive molecular forces are absent. 

The conclusion from all this work is that a rigorous mathematical derivation from the partition function produces a generalization of the vdW equation together with the Gibbs criterion if the attractive force is infinitely weak with an infinitely long range. In that case, <math>p^0</math>, the pressure that results from direct particle collisions (or more accurately the core repulsive forces), replaces <math>RT/(v-b)</math>. This is consistent with the second criticism, which can be stated as <math>p^0\ne RT/(v-b)\text{ for all }b\le v<\infty</math>. Consequently, the vdW equation cannot be rigorously derived from the configuration integral over the entire range of <math>v</math>. 

Nevertheless, it is possible to derive the vdW equation using the [[virial expansion]] in the limits <math>\tau, \rho b\rightarrow 0</math>. Consequently the equation produces an accurate approximation in a region defined by <math>v\gg b,\, T\gg\varepsilon/k</math> (on a molecular basis <math>\rho_N\sigma^3\ll 1</math>), which corresponds to a dilute gas. But as the density increases, the behavior of the vdW approximation and the two-term virial expansion differ markedly. Whereas the virial approximation in this instance either increases or decreases continuously, the vdW approximation together with the Maxwell construction expresses physical reality in the form of a phase change, while also indicating the existence of metastable states. This difference in behaviors was pointed out by Korteweg{{sfn|ps=|Korteweg|1891a|p=}} and Rayleigh (see Rowlinson{{sfn|ps=|van der Waals|1984|p=20}}) in the course of their dispute with [[Peter Guthrie Tait|Tait]] about the vdW equation.  

===Empirical justification===
In this extended region, use of the vdW equation is not justified mathematically; however, 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 [[Berni Alder|Alder]], et al. who, at a conference marking the 100th anniversary of van der Waals' thesis, noted that:{{sfn|ps=|Alder|Alley|Rigby|1974|p=143}} 
<blockquote>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 quantatively 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.</blockquote> 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,<ref>Singer, J.V.R., and Singer, K., Mol. Phys.(1972), '''24''', 357; McDonald, J.R., (1972), '''24''', 391</ref> let alone its simplicity."{{sfn|ps=|Alder|Alley|Rigby|1974|p=144}}

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){{sfn|ps=|Valderrama|2003|p=1606}} in order to manage the liquids,{{sfn|ps=|Vera|Prausnitz|1972|p=7-10}} and gases of pure substances and mixtures,{{sfn|ps=|Kontogeorgis|Privat|Jaubert|2019|pp=4626-4629}} that they encounter in practice. 

This situation has been aptly described by Boltzmann:{{sfn|ps=|Boltzmann|1995|p=356}}
<blockquote>  
... 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.
</blockquote>

== Notes ==
{{reflist}}

== References ==

* {{cite journal |last1=Alder|first1=B. J. |last2=Alley|first2=W. E. |last3=Rigby|first3=M.|date=1974|title=Correction to the van der Waals model for mixtures and for the diffusion coefficient|journal=Physica|volume=74|issue=1 |pages=143–155|doi=10.1016/0031-8914(74)90231-6 |bibcode=1974Phy....73..143A |author1-link=Berni Alder}}
* {{cite journal |last1=Andrews |first1=T. |date=1869 |title=On the Continuity of the Gaseous and Liquid States of Matter |journal=Philosophical Transactions of the Royal Society of London |volume=159 |pages=575–590 }}
* {{citation |last1=Barenblatt |first1=G.I. |author-link1=Grigory Barenblatt |date=1979 |orig-date=1978 |title=Similarity, Self-Similarity, and Intermediate Asymptotics |translator-last1=Stein |translator-first1=Norman |publisher=Milton |publication-place=NY and London }}
* {{cite journal|last1=Barrufet|first1= M.A.|last2= Eubank|first2= P.T.|date=1989|title=Generalized Saturation Properties of Pure Fluids Via Cubic Equations of State|journal= Chemical Engineering Education|volume=23|number=3|pages= 168–175}}
* {{cite book| last = Boltzmann	| first = L.|author-link=Ludwig Boltzmann |date = 1995| orig-date = 1896 (Part I), 1898 (Part II) in German | title = Lectures on Gas Theory| location = NY| publisher = Dover |translator-last1= Brush |translator-first1= S.G.}}
* {{cite book|last=Callen|first= H.B.|author-link=Herbert Callen|date=1960|title= Thermodynamics|location=NY|publisher= John Wiley ans Sons}}
* {{cite journal|last=DeBoer|first= J.|date=1974| title=Van der Waals in his time and the present revival opening address|journal= Physica|volume= 73|issue= 1|pages= 1–27|doi= 10.1016/0031-8914(74)90223-7|bibcode= 1974Phy....73....1D}} 
* {{cite journal|last1=Dong|first1= W.G.|last2=Lienhard|first2= J.H.|date=1986|title=Corresponding States Correlation of Saturated and Metastable Properties|journal= Can J Chem Eng|volume= 64|pages= 158–161|doi= 10.1002/cjce.5450640123}} 
* {{cite book | last = Epstein | first = P.S. |author-link=Paul Sophus Epstein| date = 1937 | title = Textbook of Thermodynamics | location = NY | publisher = John Wiley and Sons}}
* {{cite book | last = Gibbs | first= J.W.|author-link=Josiah Willard Gibbs| date=1948|orig-date=1901|title=The Collected Works of J. Willard Gibbs Volume II Part One Elementary Principles in Statistical Mechanics | location = New Haven | publisher= Yale University Press}}
* {{cite book |last = Goodstein | first = D.L.|author-link=David Goodstein | date = 1985 |orig-date=1975| title = States of Matter |location = NY | publisher = Dover}}
* {{cite web | url = http://www.nigelhewitt.co.uk/diving/maths/vdw.html | title = Who was Van der Waals anyway and what has he to do with my Nitrox fill? | first = Nigel | last = Hewitt | work = Maths for Divers | access-date = 1 February 2019 | archive-date = 11 March 2020 | archive-url = https://web.archive.org/web/20200311003630/http://www.nigelhewitt.co.uk/diving/maths/vdw.html | url-status = live }}
* {{cite book|last=Hill|first=Terrell L.|date=1986|title=Statistical Thermodynamics|location=NY|publisher=Dover}}
* {{cite book|last1=Hirschfelder|first1= J. O.|last2= Curtis|first2= C. F.|last3=Bird|first3= R. B.|author-link1=Joseph O. Hirschfelder|author-link3=Robert Byron Bird|date=1964|title=Mollecular Theory of Gases and Liquids, corrected printing|location=NY|publisher=John Wiley and Sons, Inc.}}
* {{cite book|last=Johnston|first= D.C.|title=Advances in Thermodynamics of the van der Waals Fluid |date=2014 |chapter=|doi=10.1088/978-1-627-05532-1 | arxiv=1402.1205|bibcode= 2014atvd.book.....J|isbn=978-1-627-05532-1 }}
* {{cite journal|last=Kac|first=Marc|author-link=Marc Kac|date=1958|title=On the Partition Function of a One-Dimensional Gas|journal=Phys Fluids|volume=1|pages=8–12}}
* {{cite journal|last1=Kac|first1=M.|last2=Uhlenbeck|first2=G.E.|last3=Hemmer|first3=P.C.|date=1963|title=On the van der Waals Theory of the Vapor-Liquid Equilibrium. 1. Discussion of a One-Dimensional Model|journal=J. Math. Phys.|volume=4|issue=2 |pages=216–228|doi=10.1063/1.1703946 |bibcode=1963JMP.....4..216K }}
* {{cite book|last=Kreyszig|first= E.|author-link=Erwin Kreyszig|date=1959 |title=Differential Geometry|location=Toronto|publisher= University of Toronto Press}} 
* {{cite journal|last=Klein|first= M. J.|author-link=Martin J. Klein|date= 1974|title=The Historical Origins of the Van der Waals Equation|journal= Physica|volume=73|issue= 1|pages= 28–47|doi= 10.1016/0031-8914(74)90224-9|bibcode= 1974Phy....73...28K}} 
* {{cite journal|last1=Kontogeorgis|first1=G.M.|last2=Privat|first2=R.|last3=Jaubert|first3=J-N.J.|date=2019|title=Taking Another Look at the van der Waals Equation of State – Almost 150 Years Later|journal=J. Chem. Eng. Data|volume=64|issue=11 |pages=4619–4637|doi=10.1021/acs.jced.9b00264 }}
* {{cite journal|last=Korteweg|first=D.T.|date=1891a|author-link=Diederik Korteweg|title=On Van Der Waals Isothermal Equation|journal=Nature|volume=45|number=1155|pages=152–154|doi=10.1038/045152a0 |bibcode=1891Natur..45..152K }}
* {{cite journal|last=Korteweg|first=D.T.|date=1891b|title=On Van Der Waals Isothermal Equation|journal=Nature|volume=45|number=1160|pages=277|doi=10.1038/045277a0 }}
* {{cite journal|last=Lebowitz|first=J.L.|date=1974|title=Exact Derivation of the Van Der Waals Equation|journal=Physica|volume=73|issue=1 |pages=48–60|doi=10.1016/0031-8914(74)90225-0 |bibcode=1974Phy....73...48L }}
* {{cite journal|last1=Lebowitz|first1=J.L.|last2=Penrose|first2=O.|author-link1=Joel Lebowitz|author-link2=Oliver Penrose|date=1966|title=Rigorous Treatment of the Van der Waals–Maxwell Theory of the Liquid–Vapor Transition|journal=Jour Math Phys|volume=7|issue=1 |pages=98–113|doi=10.1063/1.1704821 |bibcode=1966JMP.....7...98L }}
* {{cite journal|last=Lekner|first= J.|date=1982|title=Parametric solution of the van der Waals liquid–vapor coexistence curve|journal= Am. J. Phys.|volume= 50|issue= 2|pages= 161–163|doi= 10.1119/1.12877|bibcode= 1982AmJPh..50..161L}}
* {{cite journal |last1=Leland |first1=T. W. |last2=Rowlinson |first2=J.S. |last3=Sather |first3=G.A. |last4=Watson |first4=I.D. |journal=Trans. Faraday Soc. |volume=65 | year=1968 }}
* {{cite journal |last = Lienhard |first = J.H. |author-link=John Lienhard| date = 1986 | title =  The Properties and Behavior of Superheated Liquids | journal = Lat. Am. J. Heat and Mass Transfer | volume = 10 | pages = 169–187}}
* {{cite journal |last1=Lienhard|first1= J.H|last2= Shamsundar|first2= N.|last3= Biney|first3= P.O.|date=1986|title=Spinodal Lines and Equations of State: A Review|journal= Nuclear Engineering and Design|volume=95|pages= 297–314|doi= 10.1016/0029-5493(86)90056-7|bibcode= 1986NuEnD..95..297L}}
* {{cite journal |last1=Lorentz |first1=H. A. |date=1881 |title=<!--unknown-->|journal=Ann. Der Physik und Chemie |volume=12}}
* {{cite journal |last=Maxwell|first=J.C.|author-link=James Clerk-Maxwell|date=1875|title= On the Dynamical Evidence of the Molecular Constitution of Bodies|journal= Nature|volume=11|number=279|pages= 357–359|doi=10.1038/011357a0 |bibcode=1875Natur..11..357C }} 
* {{cite book|last1=Moran|first1=M.J.|last2= Shapiro|first2= H.N.|date= 2000|title=Fundamentals of Engineering Thermodynamics 4th Edition| location=NY|publisher=McGraw-Hill}}
* {{cite web
 | url = https://kyleniemeyer.github.io/computational-thermo/content/mixtures/mixtures.html
 | title = Mixture properties
 | last =  Niemeyer | first = Kyle
 | work = Computational Thermodynamics
 | access-date = 2 April 2024
 | archive-date = 2 April 2024
 | archive-url = https://web.archive.org/web/20240402220705/https://kyleniemeyer.github.io/computational-thermo/content/mixtures/mixtures.html
 | url-status = live
}}
* {{cite journal | last = Peck | first = R.E. | date = 1982 | title = The Assimilation of van der Waals Equation in the Corresponding States Family | journal =  Can. J. Chem. Eng. | volume = 60 | pages = 446–449| doi = 10.1002/cjce.5450600319 }}
* {{cite journal | last1=Pitzer |first1= K.S.|last2= Lippman |first2=D.Z.|last3= Curl|first3= R.F.|last4= Huggins|first4= C.M.|last5= Peterson|first5= D.E.|author-link1=Kenneth Pitzer|author-link3=Robert Curl|date=1955 |title= The Volumetric and Thermodynamic Properties of Fluids. II. Compressibility Factor, Vapor Pressure and Entropy of Vaporization |journal=  J. Am. Chem. Soc. |volume= 77 |issue= 13|pages= 3433–3440|doi= 10.1021/ja01618a002|bibcode= 1955JAChS..77.3433P}} 
* {{cite journal |last1=Redlich |first1=O. |last2=Kwong |first2=J. N. S. |year=1949 |title=On the Thermodynamics of Solutions. V. An Equation of State. Fugacities of Gaseous Solutions. |url=https://www3.nd.edu/~powers/ame.20231/redlichkwong1949.pdf |journal=Chemical Reviews |volume=44 |issue=1 |pages=233–244 |doi=10.1021/cr60137a013 |pmid=18125401 |accessdate=2 April 2024}}
* {{cite journal|last1=Shamsundar|first1= N.|last2= Lienhard|first2= J.H.|date= 1983|title=Saturation and Metastable Properties of the van der Waals Fluid|journal= Can J Chem Eng|volume=61|issue= 6|pages= 876–880|doi= 10.1002/cjce.5450610617}}
* {{cite book |last=Sommerfeld |first=A. |editor-last1=Bopp |editor-first1=F. |editor-last2=Meixner |editor-first2=J. |translator-last=Kestin |translator-first=J. |author-link=Arnold Sommerfeld |date=1956 |url=https://books.google.com/books?id=7EejC63EUAYC |title=Thermodynamics and Statistical Mechanics – Lectures on Theoretical Physics Volume V |location=NY |publisher=Academic Press|isbn=978-0-323-13773-7 }}
* {{cite journal |last=Strutt |first=J.W., 3rd Baron Rayleigh |author-link=John William Strutt, 3rd Baron Rayleigh |date=1891 |title=On the Virial of a System of Hard Colliding Bodies |journal=Nature|volume=45|issue=1152 |pages=80–82 |doi=10.1038/045080a0 |bibcode=1891Natur..45...80R }}
* {{cite journal |last1=Su |first1=G.J. |date=1946 |title=Modified Law of Corresponding States for Real Gases |journal=Ind. Eng. Chem. |volume=38 |issue=8 |pages=803–806 |doi=10.1021/ie50440a018 }}
* {{cite book | last1 = Tien | first1 = C.L. | last2 = Lienhard | first2 = J.H. |author-link1=Chang-Lin Tien|date = 1979 | title = Statistical Thermodynamics Revised Printing | location = NY | publisher = Hemisphere Publishing| bibcode = 1979wdch.book.....T }}
* {{cite journal|last=Tonks|first=L.|author-link=Lewi Tonks|date=1936|title=The Complete Equation of State of One, Two, and Three-Dimensional Gases of Hard Elastic Spheres|journal=Phys. Rev.|volume=50|issue=10 |pages=955–963|doi=10.1103/PhysRev.50.955 |bibcode=1936PhRv...50..955T }}
* {{cite book | last1 = Truesdell | first1 = C. |last2 =Bharatha |first2 = S.|author-link1=Clifford Truesdell |date = 1977 | title = Classical Thermodynamics as a Theory of Heat Engines | location = NY | publisher = Springer-Verlag}}
* {{cite journal|last=Valderrama|first=J.O.|date=2003|title=The State of the Cubic Equations of State|journal=Ind. Chem. Eng. Res.|volume=42|issue=8 |pages=1603–1618|doi=10.1021/ie020447b }}
* {{cite thesis |last=van der Waals |first=J.D. |date=1873 |title=Over de Continuïteit van den Gas en Vloeistoftoestand |degree=Ph.D. |location=Leiden Univ }}
* {{cite book | last = van der Waals |first = J.D. |author-link=Johannes Diderik van der Waals| editor-last = Rowlinson | editor-first = J.S. |editor-link=John Shipley Rowlinson| date = 1984 | title = On the Continuity of the Gaseous and Liquid States, edited and with an Introduction by J.S. Rowlinson |  location = NY | publisher = Dover Phoenix Editions}}
* {{cite book |last1=van der Waals|first1=Johannes D.|date=1967|title=in Nobel Lectures, Physics 1901–1921|chapter=The Equation of State of Gases and Liquids|pages=254–265|location=Amsterdam|publisher=Elsevier}}
* {{cite journal|last=van Hove|first=L.|author-link=Leon Van Hove|date=1949|title=Quelques Proprieties Generales De L'Integrale De Configuration D'Un Systeme De Particules Avec Interaction|journal=Physica|volume=15|issue=11–12 |pages=951–961|doi=10.1016/0031-8914(49)90059-2 |bibcode=1949Phy....15..951V |url=https://dipot.ulb.ac.be/dspace/bitstream/2013/182044/4/c315aefe-4bb4-4bb6-9ac7-a013ab2fad89.txt }}
* {{cite book|last1=Van Wylen|first1= G.J.|last2=Sonntag|first2= R.E.|date=1973|title=Fundamentals of Classical Thermodynamics Second Edition|location=NY|publisher=John Wiley ans Sons|author1-link= Gordon Van Wylen}} 
* {{cite journal|last1=Vera|first1= J.H.|last2=Prausnitz|first2=J.M.|date=1972|title=Generalized van der Waals Theory for Dense Fluids|journal=Chem. Eng. Jour.|volume=3|pages=1–13|doi= 10.1016/0300-9467(72)85001-9}}
* {{cite book |last = Weinberg | first = S.|author-link=Steven Weinberg |date=2021|title= Foundations of Modern Physics| location=Cambridge|publisher= Cambridge University Press}}
* {{cite book |last=Whitman|first= A.M.|date=2023 |title=Thermodynamics: Basic Principles and Engineering Applications 2nd Edition|location=NY|publisher= Springer}}

== See also ==
<!--For consistecy of format, moving phrase "To solve implicit equations like Van der Waals" from the "Iteration" bullet up to here, until someone can check these wikilinks, mark the ones in the article as being so, remove ones unneeded, and for ones kept but not in the article, annotate them consistently.-->
* [[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 ==
* {{cite book | first = David | last = Chandler | year = 1987 | title = Introduction to Modern Statistical Mechanics | pages = 287–295 | location = Oxford | publisher = Oxford University Press | isbn = 0195042778}}
* {{citation | first = Michael | last = Cross | year = 2004 | chapter-url = http://www.pmaweb.caltech.edu/~mcc/Ph127/b/Lecture3.pdf | chapter = Lecture 3: First Order Phase Transitions | url = http://www.pmaweb.caltech.edu/~mcc/Ph127/b/index.html | title =  Physics 127: Statistical Physics, Second Term | location = Pasadena, California | publisher = Division of Physics, Mathematics, and Astronomy, California Institute of Technology}}.
* {{cite journal | first1 = A. | last1 = Dalgarno | first2 = W.D. | last2 = Davison | year = 1966 | title = The Calculation of Van Der Waals Interactions | journal = Advances in Atomic and Molecular Physics | volume = 2 | pages = 1–32 | doi = 10.1016/S0065-2199(08)60216-X| bibcode = 1966AdAMP...2....1D | isbn = 9780120038022 }}
* {{cite book | first1 = Charles | last1 = Kittel | first2 = Herbert | last2 = Kroemer | year = 1980 | title = Thermal Physics | edition = Revised | pages = 287–295 | location = New York | publisher = Macmillan | isbn = 0716710889}}

[[Category:Eponymous equations of physics]]
[[Category:Equations of state]]
[[Category:Gas laws]]
[[Category:Engineering thermodynamics]]
[[Category:Johannes Diderik van der Waals|Equation]]