Picture a crowded bumper-car rink. Everyone starts at roughly the same speed, but after a few minutes of crashing into each other, the speeds are all over the place — a few cars crawl, a few zoom, and most are somewhere in the middle. Nobody planned that spread; the random collisions just produced it. The molecules in any gas do exactly the same thing. They are not all moving at one "temperature speed" — they share out their energy through countless collisions until a stable pattern of speeds appears. That pattern is the Maxwell–Boltzmann distribution.
To pin it down we need three quantities: the speed of a molecule $v$, the temperature $T$ of the gas, and the mass $m$ of one molecule. The plain-language rule is: hotter or lighter means faster; colder or heavier means slower. The cleanest single number is the root-mean-square speed, a kind of energy-weighted average:
Try it for the air in this room. A nitrogen molecule has $m \approx 4.65\times10^{-26}$ kg, and at $T = 300$ K this gives $v_{\text{rms}} \approx \sqrt{3(1.38\times10^{-23})(300)/(4.65\times10^{-26})} \approx 517$ m/s — faster than a passenger jet. The molecules whizzing past you right now are doing roughly half a kilometre per second, in every direction at once.
The full curve, $f(v)\propto v^2\,e^{-mv^2/2k_BT}$, is a tug-of-war between two factors. The Boltzmann factor $e^{-mv^2/2k_BT}$ says high speeds cost a lot of energy and are rare. The geometric factor $v^2$ says there are very few ways to be almost stationary, so $v=0$ is rare too. Their product peaks in between, at the most probable speed $v_p$. That is why the distribution always has a hump rather than a single spike, and why three averages appear — $v_p < \langle v\rangle < v_{\text{rms}}$. Raising $T$ doesn't lift every molecule equally; it shifts and flattens the whole curve, while $\langle KE\rangle = \tfrac{3}{2}k_B T$ ties temperature directly to energy. (The animation runs in 2D, so its on-screen curve is $f(v)\propto v\,e^{-mv^2/2k_BT}$ and $\langle KE\rangle = k_B T$ — same physics, one fewer dimension.)
First, drag Temperature T from 300 K up to 1200 K and watch the histogram slide right and flatten — quadrupling $T$ only doubles the speeds, because $v\propto\sqrt{T}$. Second, open Different Gas Masses mode: light H₂, medium N₂ and heavy Ar share the same $T$, yet the light gas spreads far into the high-speed tail. Third, switch to Escape Velocity mode and lower $v_{\text{esc}}$ — the readout shows what fraction of molecules sit in the tail above it, which is exactly why Earth keeps its nitrogen but lets hydrogen and helium leak away.
| Symbol | Meaning | SI Unit |
|---|---|---|
| $m$ | Mass of one molecule | kg |
| $T$ | Absolute temperature | K |
| $k_B$ | Boltzmann constant ($1.381\times10^{-23}$ J/K) | J K⁻¹ |
| $v$ | Molecular speed (magnitude of velocity) | m s⁻¹ |
| $f(v)$ | Probability density of speeds | s m⁻¹ |
| $v_p$ | Most probable speed (peak of $f(v)$) | m s⁻¹ |
| $\langle v \rangle$ | Mean speed | m s⁻¹ |
| $v_{\text{rms}}$ | Root-mean-square speed | m s⁻¹ |
Slider T directly enters all three speed formulas. Slider m represents particle mass — heavier molecules ⇒ slower at the same temperature. The "Different Gas Masses" mode shows H₂ (light), N₂ (intermediate), and Ar (heavy) at the same $T$ — the histograms shift accordingly. The "Escape Velocity" mode highlights the high-speed tail: at low $v_{\text{esc}}/v_{\text{rms}}$ ratio, a significant fraction of molecules can escape (atmospheric loss). The simulation uses a 2D version with $\langle KE\rangle = k_B T$ instead of $\tfrac32 k_B T$ — physics is the same, just one fewer degree of freedom.
Atmospheric escape (light gases like H and He escape Earth's gravity faster because their high-speed tails extend past escape velocity), evaporation cooling (only the fastest molecules leave the liquid surface, lowering average kinetic energy), Doppler broadening of spectral lines (atoms moving toward/away shift emission frequency), Bose-Einstein condensation (laser cooling cuts the high-speed tail to lower $T$), and chemical reaction rates (Arrhenius factor $e^{-E_a/k_B T}$ comes directly from the high-speed tail).
A swarm of particles bouncing elastically inside a 2D box. Their initial speeds are sampled from a Maxwell-Boltzmann distribution at temperature $T$. The histogram below tracks the live speed distribution — you'll see it stay close to the theoretical curve (red line) once the system reaches equilibrium. Heat the gas (raise $T$) and the curve flattens and shifts right; cool it and it sharpens to lower speeds.
Two competing effects shape $f(v)$. The Boltzmann factor $e^{-mv^2/2k_BT}$ is largest at $v=0$ — energetically, low speeds are favoured. But the geometric factor $4\pi v^2$ (volume of a velocity-space shell) goes to zero at $v=0$ — there are very few ways to be exactly stationary. Their product peaks at $v_p$. So the peak comes from balancing energy cost against geometric availability.
They average different things. $v_p$ = location of histogram peak (most likely value). $\langle v\rangle$ = arithmetic mean (good for collision rates). $v_{\text{rms}}$ = $\sqrt{\langle v^2\rangle}$ (good for kinetic energy: $\tfrac32 k_B T = \tfrac12 m v_{\text{rms}}^2$). They are different because the distribution is asymmetric — a long high-speed tail pulls $\langle v\rangle$ above $v_p$, and pulls $v_{\text{rms}}$ even higher.
The exponent comes from energy ($\propto v^2$). The polynomial prefactor comes from dimensionality: in $D$ dimensions, the surface of a $D$-sphere of radius $v$ scales as $v^{D-1}$. So 1D has $v^0=$ constant (Gaussian), 2D has $v^1$, 3D has $v^2$. In 4 dimensions you'd get $v^3 e^{-mv^2/2k_BT}$. The $v^{D-1}$ factor is just the volume of the velocity-space shell.
Earth's gravitational escape velocity at the exobase is ~10.4 km/s. At thermospheric temperatures (~1000 K), $v_{\text{rms}}$ for hydrogen is ~5 km/s and for helium ~2.5 km/s — significant fractions of these molecules sit in the high-speed tail above escape velocity, so they leak away over geological time. For nitrogen ($v_{\text{rms}}$ ~1 km/s) the tail past 10.4 km/s is astronomically negligible. This is why all of Earth's primordial He and H₂ are gone, but mantle outgassing replenishes some He today.
Laser cooling uses photon scattering to selectively slow atoms moving toward the laser (Doppler-shifted into resonance). Each photon absorption removes a small momentum from the atom. Over millions of cycles, fast atoms get preferentially slowed, "shaving off" the high-speed tail of the distribution. Repeated stages can reach microkelvins. Removing the tail also means $\langle KE\rangle$ drops — the Maxwell-Boltzmann curve gets sharper and shifts to lower $v$. This is how Bose-Einstein condensates are made.