# Calculus

## Applications of the Derivative # Van der Waals Equation

## Conventional Derivation of the Van der Waals Equation

The state of a given amount of any substance can be described by three parameters: pressure p, volume V, and temperature T. These parameters are related to each other. Their relationship is described by the equation of state, which in the general case has the form:

$F\left( {p,V,T} \right) = 0.$

The specific form of the equation depends on the substance. For example, a rarefied gas at a sufficiently high temperature is well described by the ideal gas model. Its equation of state is the well-known ideal gas law stated by Emile Clapeyron (1799−1864) in 1834:

$pV = \frac{m}{M}RT.$

Here m is the mass of the gas, M is the molar mass (i.e. the mass of one mole of the gas), R is the universal gas constant. For one mole of gas, this equation takes the following form:

$pV = RT.$

Subsequent experiments revealed deviations in the behavior of real gases from the ideal gas law. These results were summarized by the Dutch physicist Johannes Diderik van der Waals (1837−1923) who in 1873 proposed a more accurate equation of state for a real gas.

It is called the Van der Waals equation. For one mole of a gas, it can be written as

$\left( {p + \frac{a}{{{V^2}}}} \right)\left( {V - b} \right) = RT.$

This equation takes into account the attractive and repulsive forces between molecules. The attractive forces are taken into account through the near-wall effect.

Indeed, for the particles located in the inner region, the attractive forces from other molecules are compensated on the average. However, for the particles near the walls of the container, there is an uncompensated attractive force $$f$$ directed inside the container. This force, on the one hand, is proportional to the concentration of particles $$n$$ in the container. On the other hand, it is proportional to the concentration of particles in the boundary layer. As a result, we obtain:

$f \sim {n^2} \sim \frac{1}{{{V^2}}},$

where $$n$$ is the concentration of molecules in the container, $$V$$ is the volume of one mole of gas.

The effect of attraction of the molecules of the near-wall layer reduces the pressure on the walls of the container. In the formal transition from the ideal gas law to the Van der Waals equation, this corresponds to the replacement

$p \to p + \frac{a}{{{V^2}}},$

where $$a$$ is the coefficient depending on the particular gas and the size of the container.

The repulsive forces between the molecules in the Van der Waals model are taken into account as follows: The molecules are assumed to have a spherical shape of radius $$r$$ and their centers can not come closer to each other than $$2r.$$ We can suppose that there exists a "forbidden" (excluded) volume around one of the two molecules (Figure $$2$$), equal to

$\frac{4}{3}\pi {\left( {2r} \right)^3} = 8 \cdot \frac{4}{3}\pi {r^3}.$

Hence, the excluded volume per one molecule is given by

${b_0} = 4 \cdot \frac{4}{3}\pi {r^3} = 4{V_0},$

where $${V_0}$$ is the volume of one molecule.

As a result, if for an ideal gas, the available void space where the molecules can move was equal to $$V,$$ now it becomes equal to

$V - {N_A}{b_0} = V - b,$

where $${N_A}$$ is the Avogadro's number (the number of molecules in one mole of any substance), $$b$$ is the excluded volume caused by repulsion of the particles.

Assuming that the attractive and repulsive forces between the molecules act together and replacing $$p$$ and $$V$$ with the new expressions, we can convert the ideal gas law into Van der Waals equation:

$\left( {p + \frac{a}{{{V^2}}}} \right)\left( {V - b} \right) = RT.$

## Van der Waals Isotherms

At a fixed temperature, the Van der Waals equation describes the dependence $$p\left( V \right).$$ In the $$pV$$-plane, this dependence is represented as a family of isotherms, each of which corresponds to a certain temperature. To investigate this dependency in more detail, we transform the Van der Waals equation to the following form:

$\left. {\left( {p + \frac{a}{{{V^2}}}} \right)\left( {V - b} \right) = RT\;} \right| \cdot {V^2},\;\; \Rightarrow \left( {p{V^2} + a} \right)\left( {V - b} \right) = RT{V^2},\;\; \Rightarrow p{V^3} + aV - pb{V^2} - ab - RT{V^2} = 0,\;\; \Rightarrow p{V^3} - \left( {pb + RT} \right){V^2} + aV - ab = 0\;|:p,\;\; \Rightarrow {V^3} - \left( {b + \frac{{RT}}{p}} \right){V^2} + \frac{a}{p}V - \frac{{ab}}{p} = 0.$

For a fixed value of $$p$$, the resulting equation is a third degree equation with respect to the variable $$V.$$ It is known that a cubic equation can have $$1$$ or $$3$$ real roots. The first case occurs at high temperatures $$T$$ (the green isotherm $$AB$$ in Figure $$3$$).

With lowering the temperature, an undulating region appears on the isotherm. In this case, there are three roots (the blue isotherm $$EFGHIJL$$). The transition between the two types of isotherms occurs at a certain temperature $${T_K},$$ which is called the critical temperature.

## Critical Point

At the critical point, the gas is characterized by the critical values of $${T_K}$$, $${p_K}$$ and $${V_K},$$ which are determined only by the gas properties. From the algebraic point of view, the isotherm at the critical point has three real roots, which are equal to each other. This fact makes it easy to calculate the values of $${T_K}$$, $${p_K}$$ and $${V_K}.$$ Indeed, in this case the equation of state of the gas

${V^3} - \left( {b + \frac{{RT}}{p}} \right){V^2} + \frac{a}{p}V - \frac{{ab}}{p} = 0$

should be written as

${\left( {V - {V_K}} \right)^3} = 0.$

Expanding the cube of the difference and equating the coefficients of the terms with equal powers of $$V,$$ we find the expressions for the critical parameters:

$\left( {V - {V_K}} \right)^3 = 0,\;\; \Rightarrow {V^3} - 3{V^2}{V_K} + 3V{V_K} - V_K^3 = 0.$

Take into account that at the critical point $$p = {p_K}$$, $$T = {T_K}.$$ Consequently, we obtain:

$\left\{ \begin{array}{l} - \left( {b + \frac{{RT}}{p}} \right) = - 3{V_K}\\ \frac{a}{{{p_K}}} = 3V_K^2\\ - \frac{{ab}}{{{p_K}}} = - V_K^3 \end{array} \right.,\;\; \Rightarrow \left\{ \begin{array}{l} b + \frac{{RT}}{p} = 3{V_K}\\ \frac{a}{{{p_K}}} = 3V_K^2\\ \frac{{ab}}{{{p_K}}} = V_K^3 \end{array} \right..$

Divide the third equation by the second:

$\frac{{\cancel{a}b \cdot \cancel{p_K}}}{{\cancel{p_K} \cdot \cancel{a}}} = \frac{{V_K^{\cancel{3}}}}{{3\cancel{V_K^2}}},\;\; \Rightarrow {V_K} = 3b.$

Determine $${p_K}$$ from the second equation:

${p_K} = \frac{a}{{3V_K^2}} = \frac{a}{{3 \cdot {{\left( {3b} \right)}^2}}} = \frac{a}{{27{b^2}}}.$

Finally, we find the critical temperature $${T_K}$$ from the first equation:

${T_K} = \frac{{\left( {3{V_K} - b} \right){p_K}}}{R} = \frac{{\left( {3 \cdot 3b - b} \right) \cdot \frac{a}{{27{b^2}}}}}{R} = \frac{{8a}}{{27bR}}.$

Thus, the critical parameters for a Van der Waals gas depend only on the values of $$a, b$$ and are given by the formulas

${V_K} = 3b,\;\;\;{p_K} = \frac{a}{{27{b^2}}},\;\;\;{T_K} = \frac{{8a}}{{27bR}}.$

Note that we have found the critical values without the use of the derivative. It can be shown that the first and second derivatives of the function $$p\left( V \right)$$ at a critical point are equal to zero (see Example $$1$$).

## Reduced Equation of State

Using the critical parameters $${V_K}$$, $${p_K}$$ and $${T_K}$$, it is convenient to go to the dimensionless variables

$\varphi = \frac{V}{{{V_K}}},\;\;\;\pi = \frac{p}{{{p_K}}},\;\;\;\tau = \frac{T}{{{T_K}}}$

and rewrite the Van der Waals equation in reduced form:

$V = \varphi {V_K},\;p = \pi {p_K},\;T = \tau {T_K},\;\; \Rightarrow \left( {\pi {p_K} + \frac{a}{{{{\left( {\varphi {V_K}} \right)}^2}}}} \right)\left( {\varphi {V_K} - b} \right) = R\tau {T_K},\;\; \Rightarrow \left( {\frac{{\pi a}}{{27{b^2}}} + \frac{a}{{{{\left( {3b} \right)}^2}{\varphi ^2}}}} \right)\left( {3b\varphi - b} \right) = R\tau \cdot \frac{{8a}}{{27bR}},\;\; \Rightarrow \left( {\pi + \frac{3}{{{\varphi ^2}}}} \right)\left( {3\varphi - 1} \right) = 8\tau .$

This equation is more versatile than the original Van der Waals equation. The isotherms of various substances built in accordance with the above equation will coincide at the same values of $$\tau.$$ The states of different substances, which are described by the same reduced parameters $$\pi,$$ $$\varphi,$$ $$\tau$$ are called the corresponding states. If two bodies or substances have equal values of any two of the three reduced quantities, then they have the same third quantity as well. This universal property is called the law of corresponding states.

Consider another interesting invariant. Since the critical quantitites $${V_K}$$, $${p_K}$$, $${T_K}$$ in the Van der Waals model are expressed only in terms of the two parameters $$a$$ and $$b,$$ we can construct a combination of these quantities, which is independent of $$a$$ and $$b.$$ One can write the following relationship:

$\frac{{{p_K}{V_K}}}{{{T_K}}} = \frac{{3b \cdot \frac{a}{{27{b^2}}}}}{{\frac{{8a}}{{27bR}}}} = \frac{{3\cancel{a}\cancel{b} \cdot \cancel{27}\cancel{b}R}}{{\cancel{27}\cancel{b^2} \cdot 8\cancel{a}}} = \frac{{3R}}{8}\;\;\;\text{or}\;\;\;\frac{{{p_K}{V_K}}}{{{T_K}R}} = \frac{3}{8} = 0,375.$

This ratio is called the universal compressibility factor and does not depend on any quantities characterizing the substance, that is this ratio is an invariant.

## Liquid-Gas Phase Transition and Metastable States

Consider again the Van der Waals equation at temperatures above the critical point: $$T \lt {T_K}.$$ In this case, the isotherms have an undulating shape. A portion of such isotherm between the minimum and maximum is characterized by a positive derivative $${\frac{{dp}}{{dV}}} \gt 0,$$ which corresponds to an unstable state of the substance. Any small positive perturbation of the volume $$dV \gt 0$$ causes an increase in the gas pressure (since $$dp \gt 0$$), which will lead to an explosive expansion of the gas. In the case of a negative perturbation $$dV \lt 0,$$ on the contrary, we would observe a quick collapse of the gas.

Actually, the liquid-gas phase transition occurs in this area of the $$PV$$-diagram. This transition is accompanied by a significant change in the volume $$V$$ at constant pressure $$p$$ and constant temperature $$T.$$ Such a transition is represented by a straight horizontal line in the $$PV$$-diagram (the segment line $$FJ$$ in Figure $$4$$).

The two phases coexist in the container at this moment. As the volume of the gas decreases when moving to the left along the segment line $$JF,$$ the volume of the liquid correspondingly increases.
The location of the horizontal section $$FJ$$ is determined from thermodynamic considerations − on the basis of the so-called Maxwell equal area rule, according to which the areas of the curved shapes $$FGH$$ and $$HIJ$$ should be equal (see Figure $$3$$ above).
The set of start and end points of the horizontal segment lines at different temperatures $$T$$ (where $$T \le {T_K}$$) forms a curve, which is called the binodal curve. The curve passing through the minimum and maximum points of all isotherms at $$T \le {T_K}$$ is called the spinodal curve.
The unstable states (for which $${\frac{{dp}}{{dV}}} \gt 0$$) are within the region bounded by the spinodal curve. A part of the diagram between the spinodal and binodal curves in principle satisfies the stability criterion $${\frac{{dp}}{{dV}}} \lt 0.$$ Therefore, the states in this region are called metastable states. The left half of the indicated region represents the superheated liquid, and the right half corresponds to the supercooled vapor.
It is surprising that such a simple model of a real gas, which takes into account the intermolecular forces only very roughly is capable to reflect the actual set of different states of gas and liquid. Besides the regular gas and liquid, the "liquid-gas" phase transition, and the metastable states considered above, the Van der Waals model also captures such an exotic phenomenon as stretched liquid (the curve $$MNS$$ in Figure $$4\text{).}$$ This state of the liquid is unstable.