# 1.6: Another way of dealing with real gases: the virial equation

- Page ID
- 321724

The van der Waals equation is an improvement on the ideal gas law, but it still has weaknesses. More sophisticated equations (Redlich-Kwong and Peng-Robinson are the two most prominent) use the same basic underpinning of correcting the ideal gas law with terms that represent molecular volume and molecular attraction.

On the other hand, a rawer mathematical approach to describing deviations from ideality is to express those deviations as part of a polynomial series where the constants have less intentional physical meaning, but instead are simply fine-tuned corrections. These types of approaches are called *virial equations*, so named as they are expressions of the force between molecules.

The most fundamental of the virial equations multiplies the right side of the ideal gas law by a power series:

$$PV=nRT(1+\frac{nB}{V}+\frac{{n}^{2}C}{{V}^{2}}+\frac{{n}^{3}D}{{V}^{3}}+\cdots )$$

In this equation, *B, C*, and *D* (and all the values that would follow) are dubbed *virial coefficients*. In most practical applications, even a mere single correction using *B*, also known as the *second virial coefficient*, is sufficient:

$$PV=nRT(1+\frac{nB}{V})$$

Use of the third virial coefficient *C* and the fourth virial coefficient *D*, let alone further coefficients, are reserved only for specialized technical applications. (The first virial coefficient is in the first term of the parenthetical correction - simply 1 in all conditions.)

Gas |
If T = 100 K, B = |
If T = 273 K, B = |
If T = 373 K, B = |

Argon | -0.1870 L mol^{-1} |
-0.0217 L mol^{-1} |
-0.0042 L mol^{-1} |

Helium | 0.0114 L mol^{-1} |
0.0120 L mol^{-1} |
0.0113 L mol^{-1} |

Hydrogen | -0.0020 L mol^{-1} |
0.0137 L mol^{-1} |
0.0156 L mol^{-1} |

Neon | -0.0060 L mol^{-1} |
0.0104 L mol^{-1} |
0.0123 L mol^{-1} |

Nitrogen | -0.1600 L mol^{-1} |
-0.0105 L mol^{-1} |
0.0062 L mol^{-1} |

Oxygen | -0.1975 L mol^{-1} |
-0.0220 L mol^{-1} |
-0.0037 L mol^{-1} |

*Table 4: Second virial coefficients for selected gases as a function of temperature. *

*Data directly from the corresponding table in Monk's Physical Chemistry, with units converted. *

Note that while *B* has units of L mol^{-1}, it is *not* the same as the van der Waals constant *b* and does *not* represent volume of molecules; this constant is an attempt to combine all the information included in the van der Waals *a* and *b* constants into a *single* value. For that reason, the meaning behind *B* is far more dependent on temperature than the van der Waals constants. The value of *B* for oxygen at 273 K is -0.0220 L mol^{-1}; however, that value approaches negligibility (at -0.0037 L mol^{-1}) at 373 K, and becomes far more significant (-0.198 L mol^{-1}) at 100 K.

It's not listed in the table above, but the value of *B* for oxygen at 298 K is -0.0158 L mol^{-1}. If we were to determine the volume correction for our kilogram (or 31.25 mol) of oxygen gas at standard ambient temperature, a solution of the simplified virial equation would give:

$$PV=nRT(1+\frac{nB}{V})\text{}\to \text{}P=\frac{nRT}{V}(1+\frac{nB}{V})$$

$$\to \text{}P=\frac{(31.25\text{mol})(0.082058{\text{L atm mol}}^{-1}{\text{K}}^{-1})(298\text{K})}{40.0\text{L}}(1+\frac{(31.25\text{mol})(-0.0158{\text{L mol}}^{-1})}{40.0\text{L}})$$

$$\to \text{}P=19.1\text{atm}(1-0.0123)=\mathbf{\text{18.9 atm}}$$

This is a correction very similar to that given by the van der Waals equation.

Because of the dependence on temperature that the virial coefficients show, virial coefficients are not traditionally tabulated simply. Papers and entire volumes are given over to the calculation of virial coefficients at a range of specific temperatures, and those coefficients are generally derived through gas behavior or empirically determined. Reliable approximations for many applications are genuinely valuable, and worth the price of a scientific text.

It’s worth a mention, at the end of this particular discussion, that the quantity *n / V* turns up *repeatedly *in these equations, in the virial equation in particular. It’s convenient to simplify this quantity, and we usually do so by creating a definition called the *molar volume* of a gas:

$$\overline{V}=\frac{V}{n}\text{}\text{}\text{}\text{}\text{}\text{units of}\overline{V}{\text{: L mol}}^{-1}$$

This definition allows us to simplify the virial equation reasonably well:

$$PV=nRT(1+\frac{nB}{V})\text{}\to \text{}P\frac{V}{n}=RT(1+\frac{n}{V}B)\text{}\to \text{}P\overline{V}=RT(1+\frac{B}{\overline{V}})$$

We can also apply this definition to the van der Waals equation, and achieve a similar simplification:

$$(P+\frac{{n}^{2}a}{{V}^{2}})(V-nb)=nRT\text{}\to \text{}(P+\frac{{n}^{2}a}{{V}^{2}})\left(\frac{V-nb}{n}\right)=RT$$

$$\to \text{}(P+\frac{{n}^{2}}{{V}^{2}}a)(\frac{V}{n}-b)=RT\text{}\to \text{}(P+\frac{a}{{\overline{V}}^{2}})(\overline{V}-b)=RT$$

We will use this definition repeatedly to make our lives easier as we start to dig into the actual discipline of thermodynamics. Be aware that, for clarity’s sake, I make a very deliberate nomenclature choice and write the molar quantity with a bar over the symbol. This choice will take other quantities that other textbooks are less clear about and help us understand more deeply.