van der Waal’s equation of state

Johannes van der Waals, a Dutch physicist, introduced an equation in 1873 called the van der Waals equation to mathematically describe the physical state of a real gas.

\left [ p+a\left ( \frac{n}{V} \right )^2 \right ]\left ( V-nb \right )=nRT\; \; \; \; \; \; \; \; 3

where a and b are called the van der Waals coefficients, which are temperature-independent constants that are specific to each gas.

The equation looks similar to the ideal gas equation except for the pressure and volume factors of a\left ( \frac{n}{V} \right )^2 and nb respectively. To understand how van der Waals’ equation is derived, we need to know how the ideal gas law is conceptualised.

The ideal gas equation, pV = nRT, is derived from the combined experiments of Robert Boyle and Jacques Charles and a principle developed by Amedeo Avogadro. This means that, for a fixed amount of gas at a relatively high temperature, the product of the experimentally observed values of p and V in the equation is a constant. The equation works reasonably well at low pressures where gases approach ideality, i.e. gas particles can be assumed to be point masses with zero volume and therefore have no intermolecular forces of attraction and repulsion between them. In other words, the equation only works when the observed variables of p and V are ‘ideal’:

p_{ideal}V_{ideal}=nRT\; \; \; \; \; \; \; \; 4

In reality, gases have finite volumes and repel one another, especially at high pressures. The experimentally observed values of pressures and volumes are now preal and Vreal. For eq4 to work under such conditions, we have to modify it to include the factors: preal and Vreal.

Due to the presence of repulsive forces between real gas molecules at high pressures, the observed volume for a real gas is larger than the volume of its ideal counterpart, i.e. Vreal > Videal. The diagram below shows the volume expansion (from yellow to blue) if the point masses have finite volumes.

Letting the difference in molar volume of a real gas and an ideal gas be b, we have \frac{V_{real}}{n}-\frac{V_{ideal}}{n}=b or

V_{ideal}=V_{real}-nb\; \; \; \; \; \; \; \; 5

Substituting eq5 in eq4 and rearranging

V_{real}=\frac{nRT}{p_{ideal}}+nb\; \; \; \; \; \; \; \; 6

Eq6 predicts that the molar volume of a gas at zero Kelvin is b. Experiments indeed reveal that b has approximately the same value as the molar volume of the solid or liquid form of the gas. On the other hand, the ideal gas law, V_{ideal}=\frac{nRT}{p_{ideal}}, wrongly predicts that the volume of a gas at 0 K is zero.

Next, the pressure that a real gas exerts on the wall of a container is lower than the pressure of its ideal counterpart (pideal > preal), as the molecules of a real gas experience intermolecular forces of attraction that decelerates their motion when they are about to impact the wall of the container. This means that we need to add a pressure shortfall, Δpreal, to the observed real pressure to reflect the ideal pressure that fits eq4:

p_{ideal}=p_{real}+\Delta p_{real}\; \; \; \; \; \; \; \; 7

If we perceive the intermolecular force of attraction (i.e. the force that one molecule of gas experiences near the wall of the container) as the attraction between that molecule and the bulk of the gas, the pressure shortfall due to that particular molecule is proportional to the density of the gas:

\Delta p_{real,\, one\, molecule}=k\left ( \frac{n}{V_{real}} \right )

Furthermore, the total number of molecules N near the wall of the container is again proportional to the density of the gas:

N=k'\left ( \frac{n}{V_{real}} \right )

The total pressure shortfall is therefore:

\Delta p_{real}=N\Delta p_{real,\, one\, molecule}=k'\left ( \frac{n}{V_{real}} \right )k\left ( \frac{n}{V_{real}} \right )=a\left ( \frac{n}{V_{real}} \right )^2

where a = k’k is a proportionality constant that is specific to each gas. Eq7 becomes

p_{ideal}=p_{real}+a\left ( \frac{n}{V_{real}} \right )^2\; \; \; \; \; \; \; \; 8

Substituting eq5 and eq8 in eq4, we get the explicit form of eq3:

\left [ p_{real}+a\left ( \frac{n}{V_{real}} \right )^2 \right ]\left ( V_{real}-nb \right )=nRT\; \; \; \; \; \; \; \; 9

The pressure and volumes in eq9 now represent experimentally observed real pressure and experimentally observed real volumes respectively. We have essentially converted an equation that only works for ideal gases into one that is applicable for real gases. For simplicity, we can omit the subscripts, which gives eq3.

Using eq3, the van der Waals equation can also be expressed in the molar form by substituting it with the molar volume, V_m=\frac{V}{n}, to give:

\left ( p+\frac{a}{V_m^{\: 2}} \right )\left ( V_m-b \right )=RT

which rearranges to:

p=\frac{RT}{V_m-b}-\frac{a}{V_m^{\: 2}}\; \; \; \; \; \; \; \; 10

At high molar volumes (i.e. low pressures), VmbVm and the second term on the RHS of eq10 becomes very small, with eq10 reducing to the ideal gas equation: p\approx \frac{RT}{V_m}.

The van der Waals equation is an improvement over eq2 as it is able to describe the physical state of a real gas in terms of the van der Waals coefficients a and b, which relate to the strength of the forces of attraction between molecules and the size of the molecules respectively.

Since the critical isotherm of a gas has an inflexion point at the critical point, we can determine the van der Waals coefficients for a gas by finding the 1st and 2nd derivatives of eq10, letting these be equal to zero and solving for the critical constants:

\frac{dp}{dV_m}=-\frac{RT}{\left ( V_m-b \right )^2}+\frac{2a}{V_m^{\: 3}}=0\; \; \; \; \; \; \; \; 11

\frac{dp}{dV_m}=\frac{2RT}{\left ( V_m-b \right )^3}-\frac{6a}{V_m^{\:4}}=0\; \; \; \; \; \; \; \; 12

From eq11 and eq12, a=\frac{RTV_c^{\: 3}}{2\left ( V_c-b \right )^2} and a=\frac{RTV_c^{\: 4}}{3\left ( V_c-b \right )^3} respectively. Combining both values of a,

V_c=3b\; \; \; \; \; \; \; \; 13

Substitute eq13 back in eq11 where Vm is now Vc,

T_c=\frac{8a}{27Rb}\; \; \; \; \; \; \; \; 14

Substitute eq13 and eq14 back in eq10,

p_c=\frac{a}{27b^2}\; \; \; \; \; \; \; \; 15

Hence, by measuring the pressure, volume and temperature at the critical point of a gas (i.e. the critical constants), we can calculate the van der Waals coefficients for that gas. The table below shows the van der Waals coefficients for some common gases:

a / atm dm6 mol-2

b / 10-2 dm3 mol-1

H2

0.2420 2.65

N2

1.352

3.87

O2

1.364

3.19

CO2 3.610

4.29

 

Next article: Virial equation of state
Previous article: Critical constants
Content page of real gases
Content page of intermediate chemistry
Main content page

Leave a Reply

Your email address will not be published. Required fields are marked *