Law of Corresponding States

Photograph the liquid–gas coexistence dome of argon in its own pressure–temperature axes, then do the same for nitrogen, methane, carbon dioxide. The four pictures look different — the substances boil, condense, and reach their critical points at wildly different pressures and temperatures. Now relabel every axis in units of that substance's own critical value. The four photographs snap into one. This is the law of corresponding states: measured in reduced variables Pr=P/PcP_r=P/P_c, Vr=V/VcV_r=V/V_c, Tr=T/TcT_r=T/T_c, the equation of state of every simple fluid collapses onto a single universal surface. It is one of the oldest and most useful unifying ideas in the thermodynamics of matter, first stated by van der Waals in 1880, and this article owns both its derivation and the honest accounting of where it succeeds and fails. Two routes reach the law. The first is almost trivial and purely algebraic: any two-parameter equation of state (van der Waals, Redlich–Kwong, Berthelot) obeys a corresponding-states law automatically, because its two molecular constants can be traded for the two independent critical coordinates, leaving a relation among reduced variables with no substance-specific information left in it. The van der Waals equation reduces exactly to (Pr+3/Vr2)(3Vr−1)=8Tr(P_r+3/V_r^2)(3V_r-1)=8T_r, a parameter-free curve. The second route is deeper and statistical-mechanical: if the intermolecular pair potential is a single energy scale times a universal shape, u(r)=ε ϕ(r/σ)u(r)=\varepsilon\,\phi(r/\sigma), then the classical configuration integral — and therefore every thermodynamic function — depends only on the two reduced groups T∗=kBT/εT^*=k_BT/\varepsilon and ρ∗=ρσ3\rho^*=\rho\sigma^3. All fluids sharing the shape ϕ\phi are thermodynamically congruent, not by coincidence but by a rescaling that removes their only two dimensional differences. The sharpest single-number test is the critical compressibility factor Zc=PcVc/(NkBTc)Z_c=P_cV_c/(Nk_BT_c). Van der Waals predicts the universal value Zc=3/8=0.375Z_c=3/8=0.375. Experiment confirms that ZcZ_c is indeed nearly constant across simple fluids — but smaller: argon 0.2910.291, nitrogen 0.2900.290, methane 0.2880.288, carbon dioxide 0.2740.274, water 0.2290.229. The near-constancy vindicates the principle; the systematic shortfall below 3/83/8 convicts the van der Waals potential of being too crude. That double verdict — right idea, wrong number — is the whole story of mean-field corresponding states, and its failure near the true critical point (where fluctuations enforce Ising, not mean-field, exponents) is exactly the doorway to universality and the renormalization group later in this track.

Key takeaways
  • When expressed in the reduced variables Pr=P/Pc,  Vr=V/Vc,  Tr=T/TcP_r=P/P_c,\;V_r=V/V_c,\;T_r=T/T_c, the equations of state of simple fluids approximately collapse onto a common universal form, f(Pr,Vr,Tr)=0f(P_r,V_r,T_r)=0.
  • It follows exactly from any two-parameter equation of state; van der Waals reduces to the parameter-free (Pr+3/Vr2)(3Vr−1)=8Tr(P_r+3/V_r^2)(3V_r-1)=8T_r.
  • More deeply it holds whenever the pair potential is a single scale times a universal shape, u(r)=ε ϕ(r/σ)u(r)=\varepsilon\,\phi(r/\sigma), so classical thermodynamics depends only on T∗=kBT/εT^*=k_BT/\varepsilon and ρ∗=ρσ3\rho^*=\rho\sigma^3.
  • The critical compressibility factor is the sharpest prediction: van der Waals gives Zc=3/8=0.375Z_c=3/8=0.375; real simple fluids cluster near 0.290.29 (Ar 0.2910.291, N₂ 0.2900.290, CH₄ 0.2880.288), below the mean-field value.
  • The near-constancy confirms the principle; the shortfall below 3/83/8 shows the van der Waals form is qualitatively right but quantitatively crude — Pitzer's acentric factor ω\omega restores accuracy.
  • Key idea: strip out the two molecular scales — an energy ε\varepsilon and a length σ\sigma (equivalently a,ba,b, equivalently Pc,TcP_c,T_c) — and every simple fluid's equation of state becomes the same dimensionless surface; what remains is universal.
The core insight

Strip out the two molecular scales — an energy ε\varepsilon and a length σ\sigma (equivalently aa and bb, equivalently PcP_c and TcT_c) — and the equation of state of every simple fluid becomes the same dimensionless surface; what remains is universal.

Where this lives today

The law of corresponding states is the daily workhorse of chemical and petroleum engineering. Given only a substance's critical constants Pc,TcP_c,T_c, generalized compressibility charts (the Nelson–Obert charts) constructed from reference-fluid data predict thermodynamic properties over wide ranges of temperature and pressure — a single universal chart replacing full PP–VV–TT measurements for every new chemical. Pitzer's three-parameter extension, adding the acentric factor ω\omega to capture departures from spherical symmetry, underlies the cubic equations of state (Soave–Redlich–Kwong, Peng–Robinson) that run process simulations worldwide. The reduced second virial coefficient collapse gives a low-density, assumption-light test that helium and argon are, at heart, the same fluid at different scales. Its deepest legacy is conceptual. Corresponding states is the classical, mean-field ancestor of universality. Guggenheim's 1945 coexistence collapse already carried the clue: his data followed an exponent near 1/31/3, not the mean-field 1/21/2, an early empirical signal that the true critical universality is not van der Waals but Ising. The modern statement is that all simple fluids flow, under the renormalization group, to the same 3D-Ising fixed point, and so share exact critical exponents and scaling functions — a far more profound "corresponding states" than the analytic van der Waals collapse, and one that only makes sense once fluctuations, not mean fields, are in charge. The law of corresponding states thus sits exactly at an important conceptual bridge between classical mean-field descriptions and modern theories of critical phenomena.

The physical picture

Each fluid carries its own natural scales—the critical pressure, temperature, and volume—which serve as reference units for reduced variables. Argon's critical temperature (151151 K) and pressure (48.648.6 atm) are set by the depth and range of the argon–argon interaction; methane's are set by methane's. But if the two interactions have the same shape and differ only in overall depth and size, then measuring each fluid in units of its own critical point divides out precisely those two differences. The residual, dimensionless equation of state has nothing substance-specific left to distinguish the fluids — so they coincide. Corresponding states is the macroscopic shadow of a microscopic sameness. At the molecular level a fluid "feels" temperature only through the ratio T∗=kBT/εT^*=k_BT/\varepsilon — thermal energy against the depth of the attractive well — and density only through ρ∗=ρσ3\rho^*=\rho\sigma^3 — how many molecular volumes are packed together. Two fluids with the same well shape ϕ\phi, differing only in depth ε\varepsilon and diameter σ\sigma, are running the identical movie at different playback speeds and zoom levels. Reduced variables synchronize the playback frame for frame. The first figure illustrates the approximate collapse of coexistence curves for several simple fluids when plotted in reduced variables.

Figure
Coexistence-curve collapse: Ar, N₂, CH₄, CO₂ plotted as TrT_r vs reduced density fall on one universal dome (Guggenheim 1945); the van der Waals dome sits nearby but its apex Z_c = 0.375 overshoots the measured ≈ 0.29
Reduced coexistence dome (data collapse) ρ/ρ_c → T_r 1.0 universal data dome (Ar,N₂,CH₄,CO₂) van der Waals dome (Z_c = 0.375, apex too high) measured apex, Z_c ≈ 0.29

The number ZcZ_c is a compact litmus test. It measures how non-ideal a fluid is exactly at its critical point: Zc=1Z_c=1 would be an ideal gas, and the actual value tells you how far from ideal the fluid sits at criticality. A universal value of ZcZ_c would imply that all simple fluids have the same reduced equation of state at the critical point. The measured values are nearly constant near 0.290.29 — the principle is real — but they are not the van der Waals 3/8=0.3753/8=0.375, which is the fingerprint of a too-crude potential. The second figure lays the prediction against the data directly.

Figure
The ZcZ_c verdict: van der Waals predicts the universal Z_c = 0.375 (amber line), but measured simple fluids cluster near 0.29 (Ar 0.291, N₂ 0.290, CH₄ 0.288, CO₂ 0.274) — near-constant, confirming the principle, but systematically below 3/8
Critical compressibility factor Z_c 0.40 0.25 vdW prediction Z_c = 0.375 Ar0.291 N₂0.290 CH₄0.288 CO₂0.274 H₂O0.229 measured cluster ≈ 0.29 (near-universal, below 3/8)

The core equations

Begin from the van der Waals equation in its extensive form, the two-parameter equation whose reduction demonstrates the law most cleanly. It is owned by aq-sm-van-der-waals-equation; here it is the vehicle.

(P+aN2V2)(V−Nb)=NkBT(1)\left(P+\frac{aN^2}{V^2}\right)(V-Nb)=Nk_BT \qquad (1)

where: aa — attraction parameter (energy × volume, ∝εσ3\propto\varepsilon\sigma^3); bb — excluded volume per molecule (∝σ3\propto\sigma^3); NN — molecule number; V,P,TV,P,T — volume, pressure, temperature. The two molecular constants aa and bb are the only substance-specific inputs. Imposing the horizontal-inflection conditions on the critical isotherm locates the critical point, which will serve as the substance's built-in unit system.

(∂P∂V)T=(∂2P∂V2)T=0 ⇒ Pc=a27b2,Vc=3Nb,kBTc=8a27b(2)\left(\frac{\partial P}{\partial V}\right)_T=\left(\frac{\partial^2 P}{\partial V^2}\right)_T=0\ \Rightarrow\ P_c=\frac{a}{27b^2},\quad V_c=3Nb,\quad k_BT_c=\frac{8a}{27b} \qquad (2)

where: Pc,Vc,TcP_c,V_c,T_c — the critical pressure, volume, temperature (derived at length in aq-sm-van-der-waals-equation). Note the three critical coordinates depend on the two constants a,ba,b, so they are not independent — a fact that itself forces a universal ZcZ_c. Trading a,ba,b for the critical coordinates and rescaling every variable by its critical value collapses the equation to a form with no molecular constants left.

(Pr+3Vr2)(3Vr−1)=8Tr(3)\left(P_r+\frac{3}{V_r^2}\right)(3V_r-1)=8T_r \qquad (3)

where: Pr=P/PcP_r=P/P_c, Vr=V/VcV_r=V/V_c, Tr=T/TcT_r=T/T_c — reduced pressure, volume, temperature. This is the reduced van der Waals equation of state. After scaling by the critical constants, every van der Waals fluid obeys the same dimensionless equation, independent of aa and bb. Evaluating the compressibility factor at the critical point gives the law's single most testable number.

Zc=PcVcNkBTc=(a/27b2)(3Nb)N(8a/27b)=38=0.375(4)Z_c=\frac{P_cV_c}{Nk_BT_c}=\frac{(a/27b^2)(3Nb)}{N(8a/27b)}=\frac{3}{8}=0.375 \qquad (4)

where: ZcZ_c — critical compressibility factor. The molecular constants cancel identically, leaving a pure number — predicted universal, and this is exactly the claim experiment scrutinizes. More generally, the corresponding-states principle follows when the intermolecular potential contains only one characteristic energy scale and one characteristic length scale.

u(r)=ε ϕ ⁣(rσ) ⇒ PVNkBT=g ⁣(kBTε,Nσ3V)≡g(T∗,ρ∗)(5)u(r)=\varepsilon\,\phi\!\left(\frac{r}{\sigma}\right)\ \Rightarrow\ \frac{PV}{Nk_BT}=g\!\left(\frac{k_BT}{\varepsilon},\frac{N\sigma^3}{V}\right)\equiv g(T^*,\rho^*) \qquad (5)

where: ε\varepsilon — well depth; σ\sigma — molecular diameter; ϕ\phi — universal dimensionless shape; T∗=kBT/εT^*=k_BT/\varepsilon; ρ∗=Nσ3/V\rho^*=N\sigma^3/V. Consequently, fluids described by the same dimensionless potential shape ϕ\phi have identical reduced thermodynamic behaviour — the law of corresponding states, now derived from the configuration integral rather than assumed from a two-parameter fit. A low-density, mean-field-free test comes from the second virial coefficient, which in reduced form is a universal function of T∗T^* alone.

B2(T)σ3=−2π ⁣∫0∞ ⁣(e−ϕ(x)/T∗−1)x2 dx≡b2∗(T∗)(6)\frac{B_2(T)}{\sigma^3}=-2\pi\!\int_0^\infty\!\left(e^{-\phi(x)/T^*}-1\right)x^2\,\mathrm{d}x\equiv b_2^*(T^*) \qquad (6)

where: B2B_2 — second virial coefficient (owned by aq-sm-second-virial-coefficient); x=r/σx=r/\sigma — reduced separation; b2∗b_2^* — universal dimensionless function of T∗T^*. Its zero defines a universal reduced Boyle temperature. Both the algebraic route (through the van der Waals critical point) and the microscopic route (through the single-scale configuration integral) are proved below and machine-checked.

Derivation AThe van der Waals critical point and the universal Z_c = 3/8

The goal is to derive the critical coordinates of the van der Waals fluid from the inflection conditions and assemble them into the compressibility factor ZcZ_c, the van der Waals model's characteristic prediction for the critical compressibility factor. The strategy is to solve the equation for the pressure, impose that the critical isotherm has a horizontal inflection (first and second volume derivatives both vanish), and solve the two simultaneous conditions for VcV_c and TcT_c, then back out PcP_c. The decisive point is that the three critical coordinates are built from only two molecular constants, so their particular combination Zc=PcVc/NkBTcZ_c=P_cV_c/Nk_BT_c must be a pure number independent of the substance — that inevitability, not the specific value, is the content of corresponding states.

Step 1. Pressure form. Solve equation (1) for pressure: P=NkBTV−Nb−aN2V2P=\dfrac{Nk_BT}{V-Nb}-\dfrac{aN^2}{V^2}. Why: the critical point is defined by the shape of the P(V)P(V) isotherm, so we need PP as an explicit function of VV at fixed TT.

Step 2. Horizontal inflection. At the critical point the isotherm is flat and inflecting: (∂P∂V)T=−NkBT(V−Nb)2+2aN2V3=0\left(\dfrac{\partial P}{\partial V}\right)_T=-\dfrac{Nk_BT}{(V-Nb)^2}+\dfrac{2aN^2}{V^3}=0 and (∂2P∂V2)T=2NkBT(V−Nb)3−6aN2V4=0\left(\dfrac{\partial^2P}{\partial V^2}\right)_T=\dfrac{2Nk_BT}{(V-Nb)^3}-\dfrac{6aN^2}{V^4}=0. Why: below TcT_c the isotherm loops (two extrema), above TcT_c it is monotonic; exactly at TcT_c the two extrema merge into a single inflection where both derivatives vanish.

Step 3. Ratio eliminates TT and aa. Dividing the second-derivative equation by the first gives 2V−Nb=3V\dfrac{2}{V-Nb}=\dfrac{3}{V}, hence Vc=3NbV_c=3Nb. Why: the ratio cancels the common NkBTNk_BT and the attraction aa, isolating a pure geometric condition on the critical volume.

Step 4. Determine TcT_c and PcP_c. Substituting Vc=3NbV_c=3Nb into the first-derivative equation gives kBTc=8a27bk_BT_c=\dfrac{8a}{27b}, and inserting both into the pressure gives Pc=a27b2P_c=\dfrac{a}{27b^2}. Why: with the critical volume fixed, one derivative condition fixes the temperature and the equation of state itself then fixes the pressure.

Step 5. Assemble ZcZ_c. Zc=PcVcNkBTc=(a/27b2)(3Nb)N(8a/27b)=38=0.375Z_c=\dfrac{P_cV_c}{Nk_BT_c}=\dfrac{(a/27b^2)(3Nb)}{N(8a/27b)}=\dfrac{3}{8}=0.375. Substituting the critical constants into Zc=PcVc/(NkBTc)Z_c=P_cV_c/(Nk_BT_c) causes all dependence on aa, bb, and NN to cancel. Why: the cancellation is the whole point — ZcZ_c is a pure number, so van der Waals predicts the same critical compressibility for every fluid; the measured spread near 0.290.29 (below 3/83/8) then measures how far real potentials depart from the van der Waals caricature. *The inflection algebra gives Vc=3NbV_c=3Nb, kBTc=8a/27bk_BT_c=8a/27b, Pc=a/27b2P_c=a/27b^2; substituting into the van der Waals equation yields the reduced EOS (Pr+3/Vr2)(3Vr−1)=8Tr(P_r+3/V_r^2)(3V_r-1)=8T_r with all constants cancelling, and Zc=3/8Z_c=3/8. Experimentally, simple fluids typically have Zc≃0.28 ⁣− ⁣0.30Z_c\simeq0.28\!-\!0.30, systematically below the van der Waals prediction 3/83/8.*

Go deeper · Derivation B · 2 minFrom a single-scale pair potential to universal reduced thermodynamics
Derivation BFrom a single-scale pair potential to universal reduced thermodynamics

The goal is to prove — with no mean-field approximation anywhere — that a pair potential of the form u(r)=ε ϕ(r/σ)u(r)=\varepsilon\,\phi(r/\sigma) forces every reduced thermodynamic function to be universal, so that corresponding states is a rigorous property of the configuration integral rather than an artifact of a two-parameter fit. The strategy is to write the classical partition function, isolate the configuration integral where all the interaction physics lives, and rescale every coordinate by the length σ\sigma. Once that is done the only dimensional combination left in the integrand is ε/kBT=1/T∗\varepsilon/k_BT=1/T^*, and the volume enters only as V/σ3V/\sigma^3; every thermodynamic function built from the integral therefore depends on T∗T^* and ρ∗\rho^* alone. The one thing to watch is the classical assumption: the thermal wavelength must be small compared to σ\sigma, or a third quantum scale survives the rescaling.

Step 1. Factor the partition function. ZN=1N! λT3NQNZ_N=\dfrac{1}{N!\,\lambda_T^{3N}}Q_N with λT=h/2πmkBT\lambda_T=h/\sqrt{2\pi mk_BT} and QN=∫∏id3ri exp⁡ ⁣[−εkBT∑i<jϕ(rij/σ)]Q_N=\int\prod_i\mathrm{d}^3r_i\,\exp\!\big[-\tfrac{\varepsilon}{k_BT}\sum_{i<j}\phi(r_{ij}/\sigma)\big]. Why: the kinetic part contributes only the thermal wavelength; all substance-specific interaction physics is confined to the configuration integral QNQ_N, so that is where the rescaling must act.

Step 2. Nondimensionalize the coordinates. Substitute ri=σ r~ir_i=\sigma\,\tilde r_i; each of the NN three-dimensional position integrals contributes a factor σ3\sigma^3, so QN=σ3N IN(T∗,V/σ3)Q_N=\sigma^{3N}\,I_N(T^*,V/\sigma^3) with IN=∫∏id3r~i exp⁡ ⁣[−1T∗∑i<jϕ(r~ij)]I_N=\int\prod_i\mathrm{d}^3\tilde r_i\,\exp\!\big[-\tfrac{1}{T^*}\sum_{i<j}\phi(\tilde r_{ij})\big]. Why: after rescaling, the integrand contains no dimensional constants — only the pure shape ϕ\phi and the ratio 1/T∗1/T^* — and the integration domain is the reduced box V/σ3V/\sigma^3; so INI_N is a pure number fixed by T∗T^*, V/σ3V/\sigma^3, and NN.

Step 3. Isolate the universal excess free energy. F−FidealNkBT=−1Nln⁡IN(V/σ3)N≡ψ(T∗,ρ∗)\dfrac{F-F_{\text{ideal}}}{Nk_BT}=-\dfrac1N\ln\dfrac{I_N}{(V/\sigma^3)^N}\equiv\psi(T^*,\rho^*). Why: subtracting the ideal-gas free energy removes the kinetic (λT\lambda_T) part entirely, provided λT≪σ\lambda_T\ll\sigma (the classical limit); what remains is a universal function of the two reduced groups, identical for every fluid sharing the shape ϕ\phi.

Step 4. Differentiate to the equation of state. P=−(∂F∂V)T,N⇒PVNkBT=g(T∗,ρ∗)P=-\left(\dfrac{\partial F}{\partial V}\right)_{T,N}\Rightarrow\dfrac{PV}{Nk_BT}=g(T^*,\rho^*). Why: a universal ψ\psi yields a universal gg; two fluids with the same ϕ\phi obey the identical reduced equation of state exactly, which is the law of corresponding states derived from first principles rather than from a two-parameter equation. The Lennard-Jones ϕ\phi is the canonical shape for the noble gases. *The σ3N\sigma^{3N} Jacobian from the coordinate rescaling and the T∗=kBT/εT^*=k_BT/\varepsilon dependence of the configuration integral confirm the universal reduced equation of state. For the Lennard-Jones shape, as for any given dimensionless potential ϕ\phi, the reduced second virial coefficient B2(T)/σ3=b2∗(T∗)B_2(T)/\sigma^3=b_2^*(T^*) is a universal function of the reduced temperature T∗T^*, and it vanishes at the reduced Boyle temperature TB∗≈3.42T^*_B\approx3.42.*

Limiting cases

The law holds precisely to the extent that a fluid is a classical, single-scale, spherically symmetric collection of particles; each way of failing that description is diagnostic.

CaseResultMeaning
Noble gases (Ar, Kr, Xe)Zc≈0.29Z_c\approx0.29, ω≈0\omega\approx0Nearly perfect Lennard-Jones spheres; the law works best here — Guggenheim's 1945 collapse used exactly these.
Quantum fluids (He, H₂, Ne)de Boer Λ∗=h/(σmε)↛0\Lambda^*=h/(\sigma\sqrt{m\varepsilon})\not\to0Light mass makes the thermal wavelength comparable to σ\sigma; a third quantum parameter is needed, and even ³He and ⁴He differ.
Polar / associating (H₂O, NH₃)ω≫0\omega\gg0, Zc≈0.23Z_c\approx0.23Directional hydrogen bonding breaks the single-scale spherical assumption; large deviations from the noble-gas surface.
Critical region, ∣T−Tc∣/Tc→0\lvert T-T_c\rvert/T_c\to0Mean-field collapse failsFluctuations dominate; the true universal behaviour is 3D-Ising (β≈0.326\beta\approx0.326, not 12\tfrac12; δ≈4.79\delta\approx4.79, not 33).

Worked examples

1 · The critical compressibility factor from the van der Waals constants (a real number)

We compute ZcZ_c for the van der Waals fluid directly from aa and bb and confirm it is substance-independent, then compare with argon.

Step 1. From the critical point, Pc=a/27b2P_c=a/27b^2, Vc=3NbV_c=3Nb, kBTc=8a/27bk_BT_c=8a/27b. Assemble Zc=PcVcNkBTc=(a/27b2)(3Nb)N(8a/27b)Z_c=\dfrac{P_cV_c}{Nk_BT_c}=\dfrac{(a/27b^2)(3Nb)}{N(8a/27b)}.

Step 2. Cancel: the numerator is 3aNb/27b2=aN/9b3aNb/27b^2=aN/9b; the denominator is 8aN/27b8aN/27b. The ratio is aN/9b8aN/27b=279⋅8=38\dfrac{aN/9b}{8aN/27b}=\dfrac{27}{9\cdot8}=\dfrac{3}{8}.

Step 3. So Zc=3/8=0.375\mathbf{Z_c=3/8=0.375}, with every factor of aa, bb, NN gone — the van der Waals prediction is the same number for every fluid.

Step 4. Argon measures Zc=0.291Z_c=0.291. The prediction is 29%29\% too high. The near constancy of ZcZ_c among many simple fluids supports the corresponding-states principle, while the systematic 0.290.29-vs-0.3750.375 gap convicts the van der Waals potential of being quantitatively crude — the right idea with the wrong well shape.

2 · The reduced Boyle temperature of the Lennard-Jones fluid (universal number)

We compute the temperature at which the second virial coefficient vanishes, in reduced units, and confirm that it is a universal reduced quantity for fluids described by the Lennard-Jones potential.

Step 1. For the Lennard-Jones shape ϕ(x)=4(x−12−x−6)\phi(x)=4(x^{-12}-x^{-6}) (with x=r/σx=r/\sigma), the reduced second virial coefficient is b2∗(T∗)=−2π∫0∞(e−ϕ(x)/T∗−1)x2 dxb_2^*(T^*)=-2\pi\int_0^\infty(e^{-\phi(x)/T^*}-1)x^2\,\mathrm{d}x, a function of the reduced temperature T∗=kBT/εT^*=k_BT/\varepsilon alone.

Step 2. At low T∗T^* the attractive well dominates: e−ϕ/T∗e^{-\phi/T^*} is large where ϕ<0\phi<0, so the integrand is positive there and b2∗<0b_2^*<0 (net attraction). At high T∗T^* the repulsive core dominates: b2∗>0b_2^*>0. Somewhere between, b2∗b_2^* crosses zero — the reduced Boyle temperature.

Step 3. Numerically integrating and root-finding gives TB∗=kBTB/ε≈3.42\mathbf{T^*_B=k_BT_B/\varepsilon\approx3.42}, a pure number set by the Lennard-Jones shape and nothing else.

Step 4. Because the reduced Boyle temperature TB∗T^*_B is the same for all Lennard-Jones fluids, their Boyle temperatures satisfy TB≈3.42 ε/kBT_B\approx3.42\,\varepsilon/k_B, scaling only with the interaction energy ε\varepsilon: argon (ε/kB≈120\varepsilon/k_B\approx120 K) gives TB≈410T_B\approx410 K, nitrogen its own value, and so on — but all collapse to the same TB∗T^*_B when plotted against T∗T^*. A single dimensionless integral fixes the Boyle temperature of every noble-gas-like fluid at once, a low-density, mean-field-free instance of corresponding states.

Origins · 1 minboyle temperature and how this was discovered
BT
boyle temperature
the temperature where the second virial coefficient vanishes and a real gas is ideal to first order in density (Boyle's 1662 law; corresponding-states reading, van der Waals 1880) · a universal reduced temperature for single-scale fluids

The Boyle temperature TBT_B is where the second virial coefficient B2(T)B_2(T) passes through zero, so that a real gas obeys PV=NkBTPV=Nk_BT to first order in density — the attractive and repulsive contributions to B2B_2 exactly cancel. It carries Robert Boyle's name because at that temperature the gas most closely obeys the pressure–volume law he published in 1662, though Boyle himself knew nothing of virial coefficients. In the corresponding-states framework the Boyle temperature is not an isolated fact but a universal marker: for any single-scale potential u=εϕ(r/σ)u=\varepsilon\phi(r/\sigma) the reduced coefficient B2/σ3B_2/\sigma^3 is a universal function of T∗=kBT/εT^*=k_BT/\varepsilon, so its zero occurs at a universal reduced temperature TB∗T^*_B — about 3.423.42 for the Lennard-Jones shape. Thus every noble-gas-like fluid has its Boyle temperature at the same TB∗T^*_B, differing only through the substance's energy scale ε\varepsilon. The Boyle temperature is therefore one of the cleanest, lowest-density tests of the law of corresponding states: it requires no mean-field approximation, only the assumption that all these fluids share one interaction shape.

History
1662
— Boyle publishes the pressure–volume law that, at the temperature later named for him, a real gas obeys to first order.
1873
— van der Waals's equation of state introduces the two molecular constants a and b.
1880
— van der Waals states the law of corresponding states, noting the collapse of reduced equations of state.
1901
— Kamerlingh Onnes casts the equation of state as the virial expansion, framing corresponding states in reduced coefficients.
1945
— Guggenheim demonstrates the coexistence-curve collapse of eight simple fluids, with an exponent near 1/3, not the mean-field 1/2.
1955
— Pitzer introduces the acentric factor ω, extending corresponding states to non-spherical molecules.

Going further

Texts. Pathria & Beale, Statistical Mechanics, 3rd ed., §12.2 (van der Waals and corresponding states). Kardar, Statistical Physics of Particles, Ch. 5 (interacting gas, van der Waals). Reif, Fundamentals of Statistical and Thermal Physics, §10.4. Sethna, Statistical Mechanics: Entropy, Order Parameters, and Complexity, Ch. 12 (continuous phase transitions, universality). Landau & Lifshitz, Statistical Physics Part 1, §84. Guggenheim, J. Chem. Phys. 13, 253 (1945). Pitzer et al., J. Am. Chem. Soc. 77, 3433 (1955). The van der Waals critical point gives the reduced EOS with all a,b,Na,b,N cancelling, Zc=3/8Z_c=3/8, and measured ZcZ_c values (Ar 0.2910.291, N₂ 0.2900.290, CH₄ 0.2880.288, CO₂ 0.2740.274, H₂O 0.2290.229, all below 3/83/8). The single-scale rescaling yields the σ3N\sigma^{3N} Jacobian and T∗T^* dependence, and the universal Lennard-Jones reduced Boyle temperature TB∗≈3.42T^*_B\approx3.42 is obtained by numerical integration. Note that a reduced Boyle temperature Tr≈2.5T_r\approx2.5 (in units of TcT_c) is sometimes quoted for noble gases — this is a different reduced group from the microscopic TB∗=kBTB/ε≈3.42T^*_B=k_BT_B/\varepsilon\approx3.42 used here; both are legitimate and consistent, as they use different reference scales, and this article states the microscopic one explicitly to avoid confusion.

Prerequisites. Van der Waals Equation of State · Maxwell Construction and Coexistence · Configuration Integral and the Interacting Partition Function · Virial Expansion of the Equation of State.

Leads to. Critical Exponents and Universality · Renormalization Group · Second Virial Coefficient.

Why do only two molecular parameters, and not three or four, appear?

Because the simplest realistic model of a simple fluid has exactly two scales: an energy (how deep the attractive well is, ε\varepsilon) and a length (how big a molecule is, σ\sigma). Everything else about the interaction is captured by the dimensionless shape ϕ(r/σ)\phi(r/\sigma), which is assumed common to the whole class. Two scales mean two dimensional constants, which can be traded for two of the three critical coordinates; the third critical coordinate is then predicted, which is why ZcZ_c comes out universal. When a fluid needs a third parameter — a quantum length for helium, or a shape/polarity parameter (Pitzer's ω\omega) for water — it is telling you the two-scale picture is incomplete for that substance.

What exactly is the difference between T^* and T_r?

They are two different reduced temperatures using two different reference scales, and both are legitimate. Tr=T/TcT_r=T/T_c measures temperature against the critical temperature — a macroscopic, directly measurable scale. T∗=kBT/εT^*=k_BT/\varepsilon measures it against the microscopic well depth ε\varepsilon — a molecular scale that comes from the pair potential. For a van der Waals fluid Tc=8a/27bkBT_c=8a/27bk_B links the two, so a given fluid has a fixed conversion factor between TrT_r and T∗T^*; for the Lennard-Jones fluid the critical temperature is kBTc≈1.3 εk_BT_c\approx1.3\,\varepsilon, so the reduced Boyle temperature TB∗≈3.4T^*_B\approx3.4 corresponds to Tr≈2.6T_r\approx2.6. This article computes and states the microscopic TB∗≈3.42T^*_B\approx3.42 to be unambiguous, since "reduced Boyle temperature" quoted without its reference scale is a common source of confusion.

Why does helium violate the classical law so badly?

Helium is light, so its thermal de Broglie wavelength λT=h/2πmkBT\lambda_T=h/\sqrt{2\pi mk_BT} is comparable to the molecular size σ\sigma near the temperatures of interest. The derivation of corresponding states assumed λT≪σ\lambda_T\ll\sigma so that quantum effects drop out of the excess free energy, leaving only the two reduced groups T∗T^* and ρ∗\rho^*. When λT∼σ\lambda_T\sim\sigma that assumption fails: a third dimensionless group, the de Boer quantum parameter Λ∗=h/(σmε)\Lambda^*=h/(\sigma\sqrt{m\varepsilon}), survives the rescaling. Helium, hydrogen, and neon all have appreciable Λ∗\Lambda^* and need a quantum law of corresponding states with Λ∗\Lambda^* as an extra axis; the difference is so real that ³He and ⁴He, chemically identical but different in mass, do not correspond to each other.

Is corresponding states an exact theorem or an empirical observation?

It is exact under a specific hypothesis and approximate in practice. If every fluid in a class truly had a pair potential u=εϕ(r/σ)u=\varepsilon\phi(r/\sigma) with one shared shape ϕ\phi, in the classical limit, then corresponding states would be an exact theorem — the rescaling argument in the second derivation proves it with no approximation. Real potentials only approximately share a shape (argon and methane are close; water is not), so in practice the law is a good approximation for simple spherical molecules and a poor one for polar or quantum fluids. The two-parameter-equation route (van der Waals) is exact for that equation but inherits the equation's own crudeness, which is why it predicts Zc=3/8Z_c=3/8 instead of the measured 0.290.29.

How is this the ancestor of universality and the renormalization group?

Corresponding states says: fluids that differ only by two scales are thermodynamically the same once rescaled. Universality says something stronger and stranger: near a critical point, fluids that differ in every microscopic detail share the same critical exponents and scaling functions, because they all flow to the same renormalization-group fixed point. The classical law is the mean-field, analytic shadow of this — it gets the collapse right away from criticality but predicts the wrong (mean-field) exponents at criticality. Guggenheim's 1945 data already showed a coexistence exponent near 1/31/3, not the mean-field 1/21/2, hinting that the true universality was different. The renormalization group, later in this track, explains why all simple fluids share the 3D-Ising exponents: corresponding states is where the idea of collapse onto a universal curve first appeared, and universality is its exact, fluctuation-driven completion.

Why is the reduced second virial coefficient a cleaner test than Z_c?

Because B2B_2 is a low-density, exact quantity — it follows directly from the pair potential with no mean-field or two-parameter-equation approximation. The reduced coefficient B2/σ3B_2/\sigma^3 is a universal function of T∗T^* for any single-scale potential, so plotting B2B_2 data for argon, nitrogen, and methane against T∗T^* (or TrT_r) collapses them onto one curve, and the Boyle temperature falls at a universal reduced value. ZcZ_c, by contrast, is a critical-point quantity where mean-field theory is at its worst, which is exactly why the van der Waals 3/83/8 is quantitatively off. So the virial test probes corresponding states in the regime where the derivation is rigorous, while ZcZ_c probes it where the derivation is weakest — both confirm the principle, but the virial test does so more cleanly.

Before this

Where this leads