Equipartition theorem Created 2026-09-24 Updated 2026-10-20
The equipartition theorem says that each independent quadratic term in a Hamiltonian described by classical mechanics has mean energy in a canonical ensemble. Here is the Boltzmann constant and the absolute temperature. Count terms in the energy, rather than vaguely counting everything that can move: one harmonic oscillator has one position degree of freedom but TWO quadratic terms, one in momentum and one in displacement.
A precise sufficient setting is with , where is one of the canonical variables, denotes all remaining variables, and the phase space integration measure factors as . The range of is the whole real line, independently of ; there is no constraint coupling its allowed values to other variables. The partition function must be finite, and the relevant expected values must exist. The Boltzmann factor then gives
The probability density function is a normal distribution with variance . A larger coefficient narrows the fluctuations enough to cancel its greater energy cost. This is why the mean contribution is independent of mass or spring constant.[ref]
For a proof outline, the one-variable Gaussian integral is , so . More generally, integration by parts in a canonical variable gives when its boundary term vanishes and the integrals converge. For this is twice its mean energy. At a finite hard boundary the boundary term need not vanish; hard-wall correction to harmonic equipartition illustrates the resulting departure. A constraint, a nonquadratic potential energy, or a noncanonical integration measure therefore requires its own calculation. The theorem is an ensemble average, not equality of instantaneous energies, and applying it to a recorded trajectory requires equilibration and sufficient sampling.
For a classical one-dimensional harmonic oscillator, with positive mass and spring constant ,
The origin is the equilibrium position. Its heat capacity is at fixed , not . Although the two mean energies agree, their measured amplitudes do not: and . In a harmonic optical trap, the position measurement observes the potential energy contribution; it does not measure both contributions by observing alone.[ref]
For particles in a three-dimensional classical monatomic gas, neglect interactions and thermally accessible internal excitations. Its three translational momentum terms give internal energy , heat capacity at constant volume , and heat capacity at constant pressure . Per mole, replace by the molar gas constant . For example, using stipulated , the molar heat capacity at fixed volume is . These results assume a dilute, nonrelativistic ideal gas with classical translational statistics, rather than a quantum degeneracy regime.
For an ideal gas of diatomic molecules, modeled as two chemical atoms with a bond, there are three translational quadratic terms, two accessible rotational kinetic energy terms for changes of bond orientation, and one bond-stretching harmonic oscillator with two quadratic terms. Rotation about the bond axis is not an additional active mode in this linear-molecule model. The three familiar limits of heat capacity at constant volume are
Corresponding molar heat capacities are , and . Counting the stretching harmonic oscillator only once as would give the wrong last value.[ref]
The frozen regimes require quantum mechanics. In the simple molecular rotation model, with state degeneracy , , and moment of inertia . For the bond-stretching quantum harmonic oscillator, with angular frequency and from the Planck constant. When , the three limits above correspond to , , and , provided translational motion remains classical. They are smooth crossovers, not abrupt switches. Real molecules can have different allowed rotational populations; at high temperature, thermal dissociation, internal excitations and nonharmonic bond motion can invalidate the fixed-molecule model before its last limit is reached.[ref]
For a quantum harmonic oscillator with fixed , put . Its exact mean energy and heat capacity are
At , the mean energy approaches and approaches . At , thermal excitation is exponentially suppressed and tends to zero, although the zero-point energy survives. Frozen thus means little thermal response, not zero motion or zero energy.[ref]
Figure 1. . Original computed equilibrium curves, not measured gas data. Left: approaches its classical value 1 as grows. Right: for stipulated . The rotational partition function is with every nonnegative integer allowed; is the variance of divided by . No molecular nuclear-spin restrictions, thermal dissociation or nonharmonic effects are modeled. Horizontal guides mark the three limiting values, not exact plateaus at finite temperature. Details.
A motivating measurement is thermal calibration of an optical trap. A focused laser through a microscope objective confines a bead in water; imaging or a quadrant photodetector records a position time series. Gieseler and colleagues describe a polystyrene bead of radius at , a trapping laser, and photodetector acquisition at . Their recorded bead continues undergoing Brownian motion; higher trapping power produces a narrower equilibrium position distribution. Infer the centered variance, not a squared distance from an arbitrary camera origin:
The normal distribution position histogram and agreement with dynamical calibration methods support the local harmonic optical trap interpretation. For a stipulated example, and imply . This calculated scale is not an experimental reading.[ref]
Position calibration is essential: if detector output has conversion , its physical variance is . An unknown scale cannot be determined together with from that one variance. A complementary calibration demonstrated by Tolić-Nørrelykke and colleagues combines the power spectrum of thermal motion with the bead response to known sinusoidal motion of a translation stage. The driven response supplies an additional length scale, and their method does not require bead size, liquid dynamic viscosity or wall distance as inputs. It still uses a thermal dynamical model; agreement between methods is a consistency check, not a theory-free proof of equipartition theorem.[ref]
The readout itself can hide fluctuations. Wong and Halvorsen used a optical trap and adjustable-exposure CCD imaging of a bead in water. They varied exposure and trapping power, held the bead away from the nearest wall, adjusted illumination to maintain comparable images, and measured tracking noise using a stuck bead. The measured variance fell as exposure increased; corrected results agreed with their time-averaging model. For a stationary overdamped harmonic optical trap, with optical trap relaxation time , an exposure of length measures . Its ideal noise-free variance is
For , : using the blurred variance in the uncorrected formula overestimates by . At , the calculated factor is . Additive tracking noise instead raises the measured variance and biases the inferred downward if left uncorrected.[ref]
Figure 2.
Equal mean quadratic energies with unequal fluctuations, and exposure-biased readout
. Original analytic models, not bead trajectories or measured curves. Left: equal-energy contours in scaled phase space, with , and , for . Increasing stiffness narrows position fluctuations but leaves both mean quadratic contributions equal to . Right: the exact noise-free exposure reduction ; at , and naive inferred stiffness is times true stiffness. The exposure calculation assumes stationary overdamped harmonic Brownian motion and uniform averaging. Details.
Useful controls are to record over many optical trap relaxation times, compare blocks for drift, calibrate the position scale, measure fixed-bead noise, vary exposure, and check the position histogram for departures from a normal distribution. A sparse sequence of accurate instantaneous positions can estimate an equilibrium variance; a long exposure averages away motion even if many frames are collected. Flow, appreciable laser heating, escape from the harmonic region or nonconservative forces can break the assumed thermal equilibrium. Inferring by imposing equipartition theorem is a calibration, whereas testing the relation requires an independently established stiffness and position scale. The beautiful practical connection is that microscopic randomness can set a force scale: once calibrated, a small constant applied force displaces the mean by . The same quadratic-term counting predicts macroscopic heat capacity and makes microscopic fluctuations useful measurements.
Write and . In the steady-state approximation for the dimer, chlorine peroxide chlorine peroxide formation is balanced by its two competing removal channels:
The association event rate is ; the associated disappearance rate of chlorine monoxide is twice that value. Each successful photodissociation event produces two chlorine atoms after the subsequent fast fragmentation, and their return to chlorine monoxide consumes two ozone molecules. Thermal dissociation instead returns the dimer to two chlorine monoxide radicals without destroying ozone. Thus the chemical branching fraction for an ozone-destroying event is , giving
This is the sunlit chlorine monoxide dimer cycle, with net conversion . It assumes that the dimer relaxes rapidly, downstream fragmentation and chlorine–ozone reactions are sufficiently fast, and competing chlorine sinks can be neglected. Counting the formation event twice and then multiplying by two again would overestimate loss by a factor of two.
Lower winter temperatures favour formation and persistence of polar stratospheric clouds, providing particle surfaces on which reservoir chlorine can be activated, for example
These are heterogeneous surface processes; cooling changes particle abundance, adsorption, phase and reaction uptake, so a simple gas-phase Arrhenius argument is insufficient. Removal of nitric acid in sedimenting particles can also delay conversion of active chlorine back into chlorine nitrate. The relevant condensation thresholds depend on pressure and water/nitric-acid abundances, rather than a single universal temperature.
Once sunlight is present, provides active chlorine. In the chlorine monoxide dimer cycle,
The net reaction is , with chlorine regenerated. Low temperature favours formation of the ClOOCl dimer and suppresses its thermal dissociation, making photodissociation rather than thermal return of two ClO radicals more competitive. Laboratory flow experiments formed the dimer at 220–240 K and distinguished its infrared and ultraviolet signatures, supporting the chemical identity needed for this mechanism.[ref]
Colder polar conditions can prepare and sustain greater springtime ozone loss, but darkness alone does not drive the photolytic cycle. In polar night reactive precursors can accumulate; sunlight is required for rapid dimer-cycle ozone destruction. This differs from the middle-stratospheric gas-phase cooling response derived above.[ref]
The radiative family is a stylized approximation to gas-pressure-supported main sequence stars with mostly radiative interiors and moderately temperature-sensitive hydrogen burning, such as the roughly solar-mass regime. Its opacity exponents are local approximations, not a universal composition-independent opacity law. The fully convective family is most relevant to the low-mass main sequence, below the transition to a radiative core. Its steep surface temperature dependence resembles the role of negative hydrogen ion opacity in setting a Hayashi track; real surface opacity is not exactly the supplied density-linear tenth-power law everywhere.
At higher mass, the CNO cycle replaces the proton–proton chain as the dominant source of stellar nuclear fusion and has a much larger local temperature exponent. This promotes a convective core. Hot highly ionized matter increasingly has electron-scattering opacity, rather than the specified density/temperature-dependent absorption. Radiation pressure becomes important in the equation of state, changing and eventually making the luminosity relation approach an Eddington luminosity scaling. Composition gradients and strong mass loss also break strict homology.
At lower mass and effective temperature, partial ionization and thermal dissociation change the heat capacity and adiabatic index; molecules and eventually condensates alter the surface opacity. Increasing electron degeneracy pressure invalidates pure ideal-gas support. Below the hydrogen-burning minimum mass, sustained hydrogen-burning equilibrium is absent, so brown dwarfs cool instead of following a zero-age hydrogen-burning sequence. These changes explain why neither derived slope should be extrapolated indefinitely.
An atmospheric thermal inversion is an altitude interval with . Absorption of incoming stellar radiation above the usual thermal-emitting layers can heat the upper gas faster than it cools, producing an inversion. In a semi-grey irradiated atmosphere, a large shortwave-to-infrared opacity ratio favors such high-altitude energy deposition; local infrared emitters and the intrinsic flux also matter.
In the Solar system, Earth and all four giant planets have well-known stratospheric inversions. The ozone layer absorbs ultraviolet sunlight on Earth; methane and photochemical hydrocarbons absorb solar radiation in the giant planets, with aerosols contributing. The giant planets are Jupiter, Saturn, Uranus and Neptune. This refers to their stratospheric temperature rise, not to the gradient at every atmospheric level.
For hot Jupiters, influential factors include the stellar flux and spectrum; the abundances of high-altitude absorbers such as titanium monoxide and vanadium monoxide; atmospheric metallicity of a giant planet and atmospheric carbon-to-oxygen ratio; thermal dissociation, atmospheric photochemistry and condensation; a atmospheric cold trap or atmospheric condensate rainout that removes absorbers; replenishment by vertical mixing; exoplanet clouds and atmospheric hazes; and heat redistribution by circulation. The ratio of visible heating to infrared cooling, rather than a single chemical species in isolation, determines whether an inversion persists.
The inferred dayside water mixing ratio is about times below the solar-equilibrium reference. Several explanations are possible. A high atmospheric carbon-to-oxygen ratio, especially near or above unity in hot carbon monoxide-dominated gas, locks much of the oxygen in CO and leaves little water. A genuinely oxygen-poor or low-metallicity envelope also reduces water, although bulk elemental abundances must apply to both hemispheres. Atmospheric photochemistry can destroy water at low pressures. In portions of the dayside that are substantially hotter than the stated minimum, thermal dissociation can also reduce water, especially at low pressure; alone does not establish strong dissociation throughout the emitting region. Temperature/opacity degeneracies or incomplete treatment of exoplanet clouds and horizontal structure in an atmospheric retrieval can bias the inferred abundance.
At the cooler terminator, is only a factor of five below the reference. Moderately reduced oxygen abundance, an enhanced C/O ratio, and exoplanet cloud/atmospheric haze dilution of spectral features are plausible explanations. Water is comparatively stable at , so strong thermal dissociation is not the natural explanation there. A exoplanet transmission spectrum samples a slant path through the limb, and the degeneracy among exoplanet cloud height, reference pressure and gas abundance can mimic a lower mixing ratio.
The two measurements need not describe the same pressure range or longitude. Local dissociation/photochemistry on the hotter dayside and reformation on the cooler limb can produce a real spatial difference. Even at one elemental C/O ratio, hot carbon monoxide-rich chemistry can leave less oxygen for water than cooler methane-rich chemistry, an example of carbon partition and atmospheric water abundance. Transport can modify or homogenize those tendencies, depending on the reaction and advection timescales. Alternatively, inconsistent assumptions in the emission and transmission retrievals can create an apparent discrepancy. A reconciliation should use one bulk elemental inventory, separate dayside/limb temperature profiles and contribution pressures, and consistent exoplanet cloud and transport physics; it should not assign independent planetary metallicities to the two hemispheres.
For dimer photodissociation frequency and competing thermal dissociation frequency , the successful branch fraction is . With steady dimer abundance and rapid subsequent recycling, ozone destruction is . The factor two counts ozone molecules consumed per productive dimer event.