Real gases and real solutions

Every formula you have built so far in this unit — chemical potential, equilibrium constants, Maxwell relations — was derived assuming an ideal gas or an ideal solution. Real CO2\text{CO}_2 in a cylinder, or real chloroform mixed with acetone, does not quite obey PV=nRTPV=nRT or Raoult's law pi=xipi∗p_i=x_ip_i^*. Rather than rederive the whole thermodynamic machinery, chemists patch it: a single corrected variable — fugacity ff in place of pressure pp, activity aia_i in place of mole fraction xix_i — is substituted into the same ideal-gas-shaped formulas, and everything downstream (equilibrium constants, Gibbs energies) carries over unchanged. This article builds that patch: the compressibility factor ZZ that measures the deviation, the van der Waals equation that models it quantitatively, and fugacity/activity that repair the thermodynamics.

The core insight

Real gases and real solutions are handled not by abandoning the ideal formulas but by replacing pp with fugacity ff and xix_i with activity aia_i inside them, with the van der Waals equation supplying the quantitative correction for gases via excluded volume bb and attractive-force parameter aa.

The core idea

Compressibility factor. The deviation of a real gas from ideal behaviour is measured directly by comparing the actual PVPV product to the ideal one:

Z=PVnRTZ = \frac{PV}{nRT}

where: PP is pressure, VV is volume, nn is moles, RR is the gas constant, TT is temperature. Z=1Z=1 for an ideal gas at every condition. Real gases show Z<1Z<1 at moderate pressure, because net intermolecular attraction pulls the gas in and makes it more compressible than ideal, and Z>1Z>1 at high pressure, where the finite size of the molecules (excluded volume) dominates and resists further compression.

Van der Waals equation. The simplest model that corrects for both effects at once adds an attractive-pressure term and subtracts an excluded-volume term from the ideal-gas law:

(P+an2V2)(V−nb)=nRT\left(P + \frac{an^2}{V^2}\right)(V-nb) = nRT

where: aa is a substance-specific constant capturing the strength of intermolecular attraction (it reduces the effective pressure the container walls feel), and bb is a substance-specific constant capturing the volume molecules physically occupy (it reduces the volume actually available for the gas to move in). Both are empirical fits, not universal constants — different for CO2\text{CO}_2, N2\text{N}_2, He\text{He}, and so on. The Boyle temperature TB=a/(Rb)T_B=a/(Rb) is the special temperature at which these two corrections cancel over a range of pressures, so the gas behaves ideally there even though a,b≠0a,b\ne0.

Principle of corresponding states. Normalising pressure, volume and temperature by the critical constants (Pr=P/PcP_r=P/P_c, Vr=V/VcV_r=V/V_c, Tr=T/TcT_r=T/T_c) puts different real gases on approximately the same footing: at equal reduced variables, most gases show approximately the same ZZ. This is why a single reduced equation of state can (approximately) describe many different substances at once.

Figure
A curve of compressibility factor Z against pressure showing the ideal gas flat line at Z equals one, a real gas dipping below one at moderate pressure from attraction and rising above one at high pressure from excluded volume
Compressibility factor Z vs pressure pressure P Z = PV/nRT ideal gas, Z = 1 Z less than 1: attraction dominates Z greater than 1: excluded volume dominates P to 0, Z to 1

Fugacity. Instead of discarding the ideal-gas chemical-potential formula, real-gas non-ideality is folded into a corrected effective pressure called fugacity, f=γPf=\gamma P, where γ\gamma is the fugacity coefficient and γ→1\gamma\to1 as P→0P\to0:

μ=μ∘+RTln⁡(fp∘)\mu = \mu^\circ + RT\ln\left(\frac{f}{p^\circ}\right)

where: μ\mu is the chemical potential, μ∘\mu^\circ the standard chemical potential, RR the gas constant, TT the temperature, and p∘p^\circ the standard-state pressure. This is formally identical to the ideal-gas formula with ff standing in for pp, so every equilibrium relation built earlier (KK, ΔrG∘=−RTln⁡K\Delta_rG^\circ=-RT\ln K) survives unchanged if fugacities/activities are used consistently throughout.

Real solutions. Analogously, deviations from Raoult's law are absorbed into an activity coefficient γi\gamma_i:

ai=γixi,μi=μi∘+RTln⁡aia_i = \gamma_i x_i, \qquad \mu_i = \mu_i^\circ + RT\ln a_i

Two reference states are used: the Raoult's-law reference (γi→1\gamma_i\to1 as xi→1x_i\to1, natural for the majority solvent) and the Henry's-law reference (γi→1\gamma_i\to1 as xi→0x_i\to0, natural for a dilute solute, where pi=KHxip_i=K_Hx_i and KH≠pi∗K_H\ne p_i^*). The excess Gibbs energy GE=Greal−GidealG^E=G^{real}-G^{ideal} (and analogous HEH^E, SES^E, VEV^E) measures, introductorily, how far a real mixture departs from ideal-solution behaviour.

Derivation ASolving the van der Waals equation for pressure

The goal is to isolate PP explicitly from the van der Waals equation (P+an2V2)(V−nb)=nRT\left(P+\dfrac{an^2}{V^2}\right)(V-nb)=nRT, so that a real-gas pressure can be computed directly from nn, VV, TT and the substance's aa, bb constants, the same way P=nRT/VP=nRT/V is computed for an ideal gas — the only new work is algebraic rearrangement, but keeping track of which term corrects for attraction and which for excluded volume is what actually matters physically.

Step 1. Start from the van der Waals equation as given and divide both sides by the volume factor (V−nb)(V-nb), which isolates the bracket containing PP: P+an2V2=nRTV−nbP+\dfrac{an^2}{V^2}=\dfrac{nRT}{V-nb}. Notice V−nbV-nb, not VV, appears here — the excluded volume nbnb has already reduced the space available to the gas, so the "ideal-like" term uses the smaller effective volume.

Step 2. Move the attraction term to the right-hand side by subtracting an2/V2an^2/V^2 from both sides: P=nRTV−nb−an2V2P=\dfrac{nRT}{V-nb}-\dfrac{an^2}{V^2}. This is now an explicit formula for PP built from two competing pieces.

Step 3. Read the two terms physically. The first term, nRT/(V−nb)nRT/(V-nb), is what an ideal gas would exert if squeezed into the smaller available volume V−nbV-nb — it is always larger than the naive ideal-gas nRT/VnRT/V, because excluded volume alone would raise the pressure. The second term, an2/V2an^2/V^2, is subtracted because intermolecular attraction pulls molecules away from the walls, lowering the pressure they actually exert. The net sign of the correction — whether PvdWP_{vdW} ends up above or below the naive ideal value — depends on which term wins at the given density and temperature.

Step 4. Confirm the ideal-gas limit. As V→∞V\to\infty (or n→0n\to0) at fixed TT, both nb/V→0nb/V\to0 and an2/V2→0an^2/V^2\to0, so P→nRT/VP\to nRT/V, the ordinary ideal-gas law — the correction terms vanish exactly when molecules are far enough apart to stop interacting or colliding, as they should.

Worked example

One mole of CO2\text{CO}_2 gas at 400 K occupies 5.0 L. Using van der Waals constants a=3.640 L2atm mol−2a=3.640\ \text{L}^2\text{atm mol}^{-2}, b=0.04267 L mol−1b=0.04267\ \text{L mol}^{-1}, find the pressure predicted by the van der Waals equation, compare it with the ideal-gas prediction, and compute ZZ.

Step 1. Ideal-gas prediction first, as a baseline: Pideal=nRTV=(1)(0.08206)(400)5.0=6.565 atmP_{ideal}=\dfrac{nRT}{V}=\dfrac{(1)(0.08206)(400)}{5.0}=6.565\ \text{atm}, using R=0.08206 L atm K−1mol−1R=0.08206\ \text{L atm K}^{-1}\text{mol}^{-1} to match the given units.

Step 2. Apply the van der Waals formula for PP derived above, P=nRTV−nb−an2V2P=\dfrac{nRT}{V-nb}-\dfrac{an^2}{V^2}. The first term: (1)(0.08206)(400)5.0−(1)(0.04267)=32.8244.9573=6.621 atm\dfrac{(1)(0.08206)(400)}{5.0-(1)(0.04267)}=\dfrac{32.824}{4.9573}=6.621\ \text{atm}. The second term: (3.640)(1)2(5.0)2=3.64025.0=0.1456 atm\dfrac{(3.640)(1)^2}{(5.0)^2}=\dfrac{3.640}{25.0}=0.1456\ \text{atm}.

Step 3. Combine: PvdW=6.621−0.1456=6.475 atmP_{vdW}=6.621-0.1456=6.475\ \text{atm}. This is lower than the ideal prediction of 6.565 atm6.565\ \text{atm} — at these moderate conditions the attractive term outweighs the excluded-volume term for CO2\text{CO}_2.

Step 4. Compute ZZ using the van der Waals pressure: Z=PVnRT=(6.475)(5.0)(1)(0.08206)(400)=32.37532.824=0.9863Z=\dfrac{PV}{nRT}=\dfrac{(6.475)(5.0)}{(1)(0.08206)(400)}=\dfrac{32.375}{32.824}=0.9863. Since Z<1Z<1, real CO2\text{CO}_2 here is slightly more compressible than an ideal gas, with actual pressure a little below the ideal-gas prediction — exactly the direction expected for a gas with significant intermolecular attraction.

Common traps

  • Applying the fugacity or activity correction and then plugging the original uncorrected pp or xix_i back into an equilibrium formula — once ff or aia_i is introduced, every downstream formula must use it consistently.
  • Confusing the Raoult's-law reference (γi→1\gamma_i\to1 as xi→1x_i\to1) with the Henry's-law reference (γi→1\gamma_i\to1 as xi→0x_i\to0) — the numerical γi\gamma_i for the same real solution differs depending on which convention is used.
  • Treating van der Waals aa and bb as universal constants rather than substance-specific empirical fits, and forgetting the equation is only a first-order correction that breaks down near the critical point.
  • Reading Z≠1Z\ne1 as "the gas disobeys thermodynamics" — non-ideality is fully absorbed by using ff (or ZZ) consistently inside the ideal-gas-shaped formulas; the framework itself remains exact.

Ask the tutor: if the Raoult's-law versus Henry's-law reference-state choice still feels arbitrary, ask the AI tutor to compute γi\gamma_i for the same chloroform data under both conventions side by side.

Why do we bother with fugacity instead of just using a different equation for real gases?

Because it lets every formula already built for ideal gases keep working without rederivation. The chemical-potential formula, the equilibrium-constant expression, the Gibbs-energy relations — all were derived assuming pp behaves ideally. Swapping in fugacity f=γPf=\gamma P for pp makes those same formulas exact for real gases too, since γ→1\gamma\to1 recovers the ideal case automatically as pressure drops. It is a bookkeeping trick that saves you from re-deriving thermodynamics from scratch for every non-ideal system.

Are the van der Waals constants a and b the same for every gas?

No — they are empirical, fit separately for each substance from real PVT data. A gas with strong intermolecular attraction, like CO2\text{CO}_2 or NH3\text{NH}_3, has a large aa; a gas made of small, weakly interacting molecules, like He\text{He}, has a small aa. Likewise bb scales with molecular size. That is also why the van der Waals equation is only a first-order model: it captures the right qualitative trends but is not exact, especially near the critical point where molecular clustering gets complicated.

Why does chloroform in acetone need the Raoult's-law reference and not Henry's-law?

The choice depends on the mole fraction range and which component you are treating as "mostly present." At xchloroform=0.2x_{chloroform}=0.2, chloroform is still a fairly major component (not a trace solute), and if we want γi→1\gamma_i\to1 as chloroform approaches purity (xi→1x_i\to1), the Raoult's-law convention is the consistent choice. If instead chloroform were present only in trace amounts, dissolved in a large excess of another solvent, the Henry's-law convention (γi→1\gamma_i\to1 as xi→0x_i\to0) would be more natural. Always state which reference state a quoted γi\gamma_i uses, since the same real data give different numerical activity coefficients under the two conventions.

Before this

Where this leads