2.3 The Maxwell–Boltzmann distribution

The pressure and equipartition results used the averages vx2\langle v_x^2\rangle and ε\langle\varepsilon\rangle without ever needing the full distribution of velocities. That distribution — how many molecules move at each speed — is fixed not by mechanics but by statistical equilibrium, and it is the Maxwell–Boltzmann distribution.

Fixing the form

In equilibrium the joint distribution of the three velocity components (vx,vy,vz)(v_x, v_y, v_z) must satisfy three conditions: it is factorisable across components (the components are independent), it is isotropic (it depends only on the speed v|\mathbf v|), and it is consistent with ε=32kBT\langle\varepsilon\rangle = \tfrac32 k_B T from 2.2. Only one functional form meets all three.

Why the component distribution is Gaussian Derivation

Independence means the joint density factorises, F(vx,vy,vz)=f(vx)f(vy)f(vz)F(v_x, v_y, v_z) = f(v_x)\,f(v_y)\,f(v_z). Isotropy means it depends only on v2=vx2+vy2+vz2v^2 = v_x^2 + v_y^2 + v_z^2. A function of a sum that is also a product of functions of the separate terms must be an exponential: lnF=lnf(vx)+lnf(vy)+lnf(vz)\ln F = \ln f(v_x) + \ln f(v_y) + \ln f(v_z) depends only on vx2+vy2+vz2v_x^2 + v_y^2 + v_z^2, which forces lnf(vx)=a+bvx2\ln f(v_x) = a + b\,v_x^2. The distribution must be normalisable, so b<0b<0; writing b=m/(2kBT)b = -m/(2k_B T) and fixing aa by normalisation,

f(vx)  =  (m2πkBT)1/2emvx2/(2kBT).f(v_x) \;=\; \left(\frac{m}{2\pi k_B T}\right)^{1/2} e^{-m v_x^2/(2k_B T)}.

The choice b=m/(2kBT)b = -m/(2k_B T) is exactly what makes 12mvx2=12kBT\langle\tfrac12 m v_x^2\rangle = \tfrac12 k_B T, the equipartition condition. ✓

To get the distribution of speeds rather than components, multiply the component density by the volume of the spherical shell of radius vv in velocity space, 4πv2dv4\pi v^2\,dv:

  f(v)  =  4π(m2πkBT)3/2v2emv2/(2kBT).  \boxed{\;f(v) \;=\; 4\pi \left(\frac{m}{2\pi k_B T}\right)^{3/2} v^2\, e^{-m v^2/(2k_B T)}.\;}

The v2v^2 from the shell pushes the peak away from zero; the exponential cuts off the high-speed tail. The distribution is therefore skewed, with a long tail toward high speeds.

03266539791306speed v (m/s)f(v)v_p⟨v⟩v_rmsc
vp (most probable)410 m/s
⟨v⟩ (mean)463 m/s
vrms502 m/s
c (speed of sound)343 m/s

The most-probable, mean, and RMS speeds always stand in the ratio √2 : √(8/π) : √3 ≈ 1.41 : 1.60 : 1.73. The speed of sound is the smallest of the four: it equals √(γ/3) ≈ 0.68 of vrms for diatomic gas. Drop the temperature and the whole distribution contracts toward zero like √T.

Slide the temperature: the distribution shifts and broadens as T\sqrt{T}. Switch gases: at the same temperature, lighter molecules (H₂, He) reach far higher thermal speeds, since ε\langle\varepsilon\rangle is fixed but mm is smaller.

Three characteristic speeds

The skew means there is no single “the” speed; three are useful, each the right one for a different question:

They sit in the fixed ratio 2:8/π:31.41:1.60:1.73\sqrt2 : \sqrt{8/\pi} : \sqrt3 \approx 1.41 : 1.60 : 1.73, independent of gas or temperature — a direct fingerprint of the distribution’s shape.