Login

Concepts

Read through concept notes with worked explanations and diagrams, chapter by chapter.

Concept 06 / 06Kinetic Theory of Gases

Real Gases, the van der Waals Equation and Deviations from Ideal Behaviour

Why the ideal gas law fails, how the a and b corrections repair it, the derivation of the critical constants, and why deviations are largest near the critical point.

Figure 1
Figure 2
Every equation in this chapter so far has rested on two idealisations: molecules are points, and they do not attract one another. Both are excellent at low density and both fail as you compress a gas. Where they fail is not a nuisance to be apologised for; it is where liquefaction, critical phenomena and refrigeration live. The van der Waals equation is the simplest correction that captures all of it. Measuring the failure: the compressibility factor Define Z=PVnRTZ = \dfrac{PV}{nRT} For an ideal gas Z=1Z = 1 at every pressure and temperature, a horizontal line. For real gases, as the graph above shows, ZZ wanders. At moderate pressures most gases dip below 11, some of them dramatically. At high pressures every gas rises above 11 and keeps climbing. Hydrogen and helium at room temperature never dip at all. Any successful theory of real gases must explain all three of those facts, and van der Waals does. The volume correction Molecules have size. If each molecule is a hard sphere, the centre of another molecule can never come closer than one diameter, so a portion of the container is permanently unavailable. Replace the free volume VV by V−nbV - nb, where bb is the excluded volume per mole. The size of bb is worth deriving because it is a neat counting argument. Two spheres of diameter dd cannot have their centres closer than dd, so each pair excludes a sphere of radius dd, of volume 43πd3\frac{4}{3}\pi d^3. Since that volume is shared between two molecules, the excluded volume per molecule is 23πd3\frac{2}{3}\pi d^3. Writing the actual volume of one molecule as v0=43π(d/2)3=16πd3v_0 = \frac{4}{3}\pi (d/2)^3 = \frac{1}{6}\pi d^3, we get 23πd3=4v0\frac{2}{3}\pi d^3 = 4v_0. So b=4NAv0b = 4 N_A v_0 The excluded volume is four times the true molecular volume, not one times. Units of bb: m3 mol−1\text{m}^3\,\text{mol}^{-1}. The pressure correction Molecules attract one another, weakly and at short range. A molecule deep in the bulk is pulled equally in all directions and feels nothing. A molecule about to strike the wall, however, has neighbours only on the inward side, so it is pulled back and arrives with less momentum than it otherwise would. The measured pressure is therefore lower than the pressure the molecules would exert if they did not attract. How large is the deficit? It is proportional to the number of molecules being pulled back near the wall, and also to the number of molecules doing the pulling. Both are proportional to density, so the deficit goes as (n/V)2(n/V)^2: Pideal=Pmeasured+a(nV)2P_{ideal} = P_{measured} + a\left(\dfrac{n}{V}\right)^2 with aa a constant measuring the strength of the attraction. Units of aa: Pa m6 mol−2\text{Pa}\,\text{m}^6\,\text{mol}^{-2}. The van der Waals equation Putting both corrections into PV=nRTPV = nRT: (P+an2V2)(V−nb)=nRT\left(P + \dfrac{a n^2}{V^2}\right)\left(V - nb\right) = nRT or for one mole, with VmV_m the molar volume, (P+aVm2)(Vm−b)=RT\left(P + \dfrac{a}{V_m^2}\right)\left(V_m - b\right) = RT Notice the physical opposition built into it. The bb term reduces the available volume, which raises the pressure relative to ideal; the aa term reflects attraction, which lowers it. Which one wins decides whether ZZ is above or below 11. Why ZZ dips and then rises Expand for one mole at moderate density. Multiply out and drop the small term ab/Vm2ab/V_m^2: PVm−Pb+aVm=RTPV_m - Pb + \dfrac{a}{V_m} = RT At low pressure Vm≈RT/PV_m \approx RT/P, so aVm≈aPRT\dfrac{a}{V_m} \approx \dfrac{aP}{RT}, giving Z=PVmRT=1+PRT(b−aRT)Z = \dfrac{PV_m}{RT} = 1 + \dfrac{P}{RT}\left(b - \dfrac{a}{RT}\right) Everything follows from the bracket. If aRT>b\dfrac{a}{RT} > b, that is at low temperature, the bracket is negative and Z<1Z < 1: attraction dominates and the gas is easier to compress than an ideal gas. If aRT<b\dfrac{a}{RT} < b, at high temperature, the bracket is positive and Z>1Z > 1 from the start: molecular size dominates and the gas resists compression. At the temperature where the bracket vanishes, TB=aRbT_B = \dfrac{a}{Rb}, the first-order deviation disappears entirely and the gas behaves almost ideally over a wide pressure range. TBT_B is called the Boyle temperature. For hydrogen TBT_B is only about 110 K110 \ \text{K}, far below room temperature, which is exactly why the hydrogen curve in the graph above never dips below 11. At very high pressure the approximation above fails and the bb term always takes over, because VmV_m approaches bb and Vm−bV_m - b approaches zero, forcing PP up without limit. That is why every gas eventually has Z>1Z > 1. Critical constants from the equation Look at the isotherm diagram above. At high temperature the van der Waals isotherms are smooth and nearly hyperbolic. Below a certain temperature they develop an S-shaped wiggle. Exactly at the boundary, the maximum and minimum of the wiggle merge into a single point of inflection with a horizontal tangent. That is the critical point, and it is found by demanding (∂P∂Vm)T=0\left(\dfrac{\partial P}{\partial V_m}\right)_T = 0 and (∂2P∂Vm2)T=0\left(\dfrac{\partial^2 P}{\partial V_m^2}\right)_T = 0 Write P=RTVm−b−aVm2P = \dfrac{RT}{V_m - b} - \dfrac{a}{V_m^2} and differentiate: ∂P∂Vm=−RT(Vm−b)2+2aVm3=0⇒RT(Vm−b)2=2aVm3\dfrac{\partial P}{\partial V_m} = -\dfrac{RT}{(V_m-b)^2} + \dfrac{2a}{V_m^3} = 0 \quad \Rightarrow \quad \dfrac{RT}{(V_m-b)^2} = \dfrac{2a}{V_m^3} ∂2P∂Vm2=2RT(Vm−b)3−6aVm4=0⇒RT(Vm−b)3=3aVm4\dfrac{\partial^2 P}{\partial V_m^2} = \dfrac{2RT}{(V_m-b)^3} - \dfrac{6a}{V_m^4} = 0 \quad \Rightarrow \quad \dfrac{RT}{(V_m-b)^3} = \dfrac{3a}{V_m^4} Dividing the first by the second kills RR, TT and aa at once: Vm−b=2Vm3⇒Vc=3bV_m - b = \dfrac{2V_m}{3} \quad \Rightarrow \quad V_c = 3b Substituting back into the first condition: RTc=2a(Vc−b)2Vc3=2a(2b)227b3=8a27b⇒Tc=8a27RbRT_c = \dfrac{2a(V_c-b)^2}{V_c^3} = \dfrac{2a(2b)^2}{27b^3} = \dfrac{8a}{27b} \quad \Rightarrow \quad T_c = \dfrac{8a}{27Rb} And then from the equation of state itself: Pc=RTcVc−b−aVc2=8a/(27b)2b−a9b2=4a27b2−3a27b2=a27b2P_c = \dfrac{RT_c}{V_c - b} - \dfrac{a}{V_c^2} = \dfrac{8a/(27b)}{2b} - \dfrac{a}{9b^2} = \dfrac{4a}{27b^2} - \dfrac{3a}{27b^2} = \dfrac{a}{27b^2} A striking prediction follows. The critical compressibility factor Zc=PcVcRTc=(a/27b2)(3b)R(8a/27Rb)=38=0.375Z_c = \dfrac{P_c V_c}{R T_c} = \dfrac{(a/27b^2)(3b)}{R\left(8a/27Rb\right)} = \dfrac{3}{8} = 0.375 contains no aa and no bb at all, so van der Waals predicts the same ZcZ_c for every substance. Real gases give 0.270.27 to 0.300.30, so this is both a triumph (a universal number does emerge, which is the origin of the law of corresponding states) and a clear demonstration that the equation is only approximate. Why deviations are worst near the critical point At the critical point (∂P/∂V)T=0(\partial P/\partial V)_T = 0, so the isothermal compressibility κT=−1V(∂V∂P)T\kappa_T = -\dfrac{1}{V}\left(\dfrac{\partial V}{\partial P}\right)_T diverges. An infinitesimal pressure change produces a huge volume change; the gas has essentially no resistance to being squeezed. Microscopically this means enormous density fluctuations: regions of gas spontaneously condense and re-evaporate on all length scales, including scales comparable to the wavelength of visible light. Those fluctuations scatter light strongly, and a normally transparent fluid turns milky. This is critical opalescence, and it is a phenomenon no equation built on "molecules are independent" could ever produce. It is exactly where the attraction term aa and the size term bb are of comparable importance and neither can be neglected, which is why both the ideal gas law and the van der Waals equation are at their worst there. Reading the isotherm diagram Above TcT_c the isotherms fall monotonically and no amount of squeezing produces a liquid. Andrews established this for carbon dioxide in 1869: above 31∘31^\circC, CO2\text{CO}_2 cannot be liquefied at any pressure. This is why oxygen and nitrogen, whose critical temperatures are 154 K154 \ \text{K} and 126 K126 \ \text{K}, resisted liquefaction for decades and were called permanent gases until people learned to cool before compressing. Below TcT_c, the van der Waals isotherm develops a loop whose middle branch has ∂P∂V>0\dfrac{\partial P}{\partial V} > 0, meaning that squeezing the gas would lower its pressure. That is mechanically unstable and cannot exist. What nature does instead is split into two phases at a single pressure, giving the flat horizontal segment drawn in the diagram. The correct pressure is fixed by Maxwell's equal-area rule: choose the horizontal line that cuts off equal areas above and below the loop. The horizontal segment is the liquid-vapour coexistence region, and the locus of its endpoints is the dome, which closes at the critical point. Key Takeaways and Formulae (P+an2V2)(V−nb)=nRT\left(P + \dfrac{an^2}{V^2}\right)(V - nb) = nRT, and Z=PVnRTZ = \dfrac{PV}{nRT}. b=4NAv0b = 4N_A v_0, four times the actual molecular volume; aa measures attraction; bb in m3mol−1\text{m}^3\text{mol}^{-1}, aa in Pa m6mol−2\text{Pa m}^6\text{mol}^{-2}. Low-pressure behaviour: Z≈1+PRT(b−aRT)Z \approx 1 + \dfrac{P}{RT}\left(b - \dfrac{a}{RT}\right). Boyle temperature TB=aRbT_B = \dfrac{a}{Rb}. Critical constants: Vc=3bV_c = 3b, Tc=8a27Rb\qquad T_c = \dfrac{8a}{27Rb}, Pc=a27b2\qquad P_c = \dfrac{a}{27b^2}, Zc=38\qquad Z_c = \dfrac{3}{8}, TBTc=278\qquad \dfrac{T_B}{T_c} = \dfrac{27}{8}. Common mistakes. The aa correction is added to PP and the bb correction is subtracted from VV, never the other way round. Both use n2n^2 and nn respectively, so watch the mole count. Z<1Z < 1 means more compressible than ideal, not less. A gas above its critical temperature cannot be liquefied by pressure alone, no matter how large the pressure. And bb is not the molecular volume; it is four times it. Brain Teaser 1 For carbon dioxide, a=3.59 atm L2mol−2a = 3.59 \ \text{atm L}^2\text{mol}^{-2} and b=0.0427 L mol−1b = 0.0427 \ \text{L mol}^{-1}. Using R=0.0821 L atm mol−1K−1R = 0.0821 \ \text{L atm mol}^{-1}\text{K}^{-1}, predict TcT_c, PcP_c and VcV_c. Compare with the measured values 304.2 K304.2 \ \text{K}, 72.9 atm72.9 \ \text{atm} and 0.094 L mol−10.094 \ \text{L mol}^{-1}, and comment on where the model does well and where it does not. Solution. Tc=8a27Rb=8(3.59)27(0.0821)(0.0427)=28.720.0947=303 KT_c = \dfrac{8a}{27Rb} = \dfrac{8(3.59)}{27(0.0821)(0.0427)} = \dfrac{28.72}{0.0947} = 303 \ \text{K} Pc=a27b2=3.5927(0.0427)2=3.590.0492=72.9 atmP_c = \dfrac{a}{27b^2} = \dfrac{3.59}{27(0.0427)^2} = \dfrac{3.59}{0.0492} = 72.9 \ \text{atm} Vc=3b=3(0.0427)=0.128 L mol−1V_c = 3b = 3(0.0427) = 0.128 \ \text{L mol}^{-1} The temperature is accurate to within 0.40.4 percent and the pressure is essentially exact. The critical volume, however, is overestimated by about 3636 percent. That pattern of errors is highly informative. Check it with ZcZ_c: Zcpredicted=38=0.375Z_c^{predicted} = \dfrac{3}{8} = 0.375, while Zcmeasured=(72.9)(0.094)(0.0821)(304.2)=6.8524.97=0.274Z_c^{measured} = \dfrac{(72.9)(0.094)}{(0.0821)(304.2)} = \dfrac{6.85}{24.97} = 0.274 Since TcT_c and PcP_c both come out right, essentially all of the error in ZcZ_c sits in VcV_c. The reason is that the hard-sphere excluded-volume model, which produced Vc=3bV_c = 3b, is the crudest part of the theory: real molecules are soft, they interpenetrate somewhat under pressure, and the effective excluded volume shrinks as the gas is compressed, so the true critical volume is smaller than 3b3b. Meanwhile TcT_c and PcP_c are governed mainly by the balance between attraction and thermal energy, which the a/V2a/V^2 term captures rather well. The practical moral: van der Waals is excellent for predicting whether a gas can be liquefied and at roughly what pressure, and poor for predicting liquid densities. As a bonus, the Boyle temperature: TB=aRb=3.59(0.0821)(0.0427)=1024 KT_B = \dfrac{a}{Rb} = \dfrac{3.59}{(0.0821)(0.0427)} = 1024 \ \text{K}, and indeed TBTc=1024303=3.375=278\dfrac{T_B}{T_c} = \dfrac{1024}{303} = 3.375 = \dfrac{27}{8}, exactly as the formulae demand. Because TBT_B for CO2\text{CO}_2 is so high, CO2\text{CO}_2 at room temperature is far below its Boyle temperature and shows a deep dip in ZZ, which is precisely what the graph above displays. Brain Teaser 2 A hypothetical gas obeys P(Vm−b)=RTP(V_m - b) = RT exactly, at all pressures. (a) Show that it can never be liquefied. (b) Find ZZ and show it always exceeds 11. (c) Does such a gas have a Boyle temperature? Solution. This gas has the size correction but no attraction, so a=0a = 0. That single fact drives everything. (a) Write P=RTVm−bP = \dfrac{RT}{V_m - b} and differentiate at constant TT: (∂P∂Vm)T=−RT(Vm−b)2\left(\dfrac{\partial P}{\partial V_m}\right)_T = -\dfrac{RT}{(V_m - b)^2} This is strictly negative for every Vm>bV_m > b and never zero. A critical point requires a point of inflection with a horizontal tangent, which requires (∂P/∂Vm)T=0(\partial P/\partial V_m)_T = 0 somewhere. Since that never happens, there is no critical point, no S-shaped loop, no coexistence dome, and therefore no liquid phase at any temperature or pressure. You can compress this gas indefinitely and it will simply get denser and stiffer, never condensing. This is the deep point of the question. Liquefaction is caused by attraction, not by compression. The bb term alone can make a gas dense; only the aa term can make it condense. (b) Z=PVmRT=VmVm−b=1+bVm−bZ = \dfrac{PV_m}{RT} = \dfrac{V_m}{V_m - b} = 1 + \dfrac{b}{V_m - b} Since Vm>bV_m > b always (the molar volume cannot fall below the excluded volume), the second term is strictly positive and Z>1Z > 1 at every pressure. Such a gas is always harder to compress than an ideal gas. In the language of the graph above, its curve starts at Z=1Z = 1 as P→0P \to 0 and rises monotonically thereafter, exactly like hydrogen at room temperature, which is a real gas so far above its Boyle temperature that the attraction term is negligible. (c) No. The Boyle temperature is where the leading deviation from ideality vanishes, TB=a/(Rb)T_B = a/(Rb), and with a=0a = 0 this gives TB=0T_B = 0. There is no finite temperature at which the gas mimics an ideal gas over a range of pressures; the positive deviation is there at every temperature, because it comes from geometry rather than from a competition between energy scales. Brain Teaser 3 Two identical rigid cylinders are filled with the same number of moles of two different gases at the same temperature, well above the critical temperatures of both. Gas A has a=0.14a = 0.14 and b=3.9×10−5b = 3.9\times 10^{-5} in SI units; gas B has a=0.36a = 0.36 and b=4.3×10−5b = 4.3\times 10^{-5}. Without computing any pressure, predict which gas has the lower pressure at moderate density and which has the higher pressure at extreme compression. Justify from the equation. Solution. Write the van der Waals pressure for one mole explicitly: P=RTVm−b−aVm2P = \dfrac{RT}{V_m - b} - \dfrac{a}{V_m^2} The two terms have opposite effects and different dependence on VmV_m, so they dominate in different regimes. The key is how fast each grows as VmV_m shrinks. At moderate density, VmV_m is still much larger than bb, so the first term is only slightly enhanced above RT/VmRT/V_m and the difference between the two gases in bb (about 1010 percent) hardly matters. The subtracted term a/Vm2a/V_m^2, however, differs by a factor of about 2.62.6. So gas B, with the much larger aa, loses more pressure to attraction and therefore has the lower pressure. Its ZZ dips further below 11, which matches the graph above, where the gas with the largest aa shows the deepest dip. At extreme compression, VmV_m approaches bb. Now (Vm−b)(V_m - b) approaches zero and the first term diverges like 1/(Vm−b)1/(V_m - b), while the attraction term a/Vm2a/V_m^2 merely approaches the finite value a/b2a/b^2. A diverging term always beats a finite one, so the repulsive size term wins and both gases end up with Z>1Z > 1. Between the two, the one whose Vm−bV_m - b vanishes first at a given VmV_m, that is the one with the larger bb, blows up sooner. Gas B has the slightly larger bb, so at extreme compression gas B reaches the higher pressure. So the ordering reverses as the gas is compressed. This crossover is the entire content of the ZZ versus PP graph: attraction is a low-and-moderate-density effect because it scales as density squared, while excluded volume is a high-density effect because it produces a divergence. Every real gas shows both, and the only question is which temperature and pressure range you happen to be looking at.