Approach of a System of Particles towards Thermal Equilibrium

In this simple model, a rigid spherical container contains identical spherical, perfectly elastic particles with the same radius and initially with the same energy. Colliding particles that overlap by cause a repelling force proportional to . This leaves each particle's energy constant for binary collisions, but redistributes the energy of more than two particles. This redistribution broadens the initially concentrated distribution to an exponential Boltzmann distribution. Such an energy distribution is known from statistical thermodynamics to apply to systems in thermal equilibrium. This Demonstration shows how amazingly fast this energy distribution is attained, although short-lived fluctuations may be observed. Control labels and graphics are explained in tooltips that appear when you mouseover them.

The total energy of the system is equal to the sum over the energies of the individual particles. The energy of each particle is the sum of its kinetic energy, its potential energy with respect to the wall, and half the sum of its potential energy of interaction with all the other particles. If one forms energy bins using only the kinetic energies (which most textbooks restrict in their treatment of gases), one would approach a Boltzmann distribution in a similar manner, as long as not too many particles overlap with one another or with the wall.

Snapshot 1: sufficiently many multiple-particle collisions have occurred to approach the Boltzmann energy distribution

Snapshot 2: a statistical fluctuation leads us away from the Boltzmann distribution

Snapshot 3: the small size of the particles has prevented a sufficient number of many-body collisions, so we remain far from the Boltzmann distribution

Snapshot 4: a state with maximum number of particles

The autorun feature animates the particle motions, according to the selected options for the initial state.