For a single-particle state of energy , put . The equilibrium mean occupancies areThe Maxwell-Boltzmann distribution is the dilute classical limit. The Fermi-Dirac distribution applies to identical fermions with at most one particle in each complete single-particle state by the Pauli exclusion principle. The Bose-Einstein distribution applies to identical bosons with no such upper occupation bound. The chemical potential fixes total particle number; these are mean state occupancies, not three separately normalized single-particle probability densities. For a number-conserving ideal Bose gas, cannot exceed the lowest single-particle energy.
Take to be a wavevector and count periodic-box states with cell volume in momentum space, parameterized by . Write for internal degeneracy. The number of states in a spherical shell is . Since ,Thus the one-particle canonical partition function isHere , for one nondegenerate branch, and for the spin-half case below. The integral uses . For classical indistinguishable particles, , givingThe factor three in the energy arises from the linear dispersion relation and three-dimensional state counting; the usual nonrelativistic quadratic-dispersion value does not apply.
For spin-half fermions at , all states up to the Fermi energy are occupied. With ,Integrating the occupied energy states givesThese are the corresponding ultrarelativistic ideal Fermi gas formulae with general slope . The nonzero ground-state internal energy and degeneracy pressure reflect the Pauli exclusion principle: particles cannot all occupy the lowest-momentum state. Formally extrapolating the classical equations gives , but Maxwell–Boltzmann statistics is not valid at fixed positive as .
At large , the number equation gives , so expanding the Fermi-Dirac distribution recovers the Maxwell-Boltzmann distribution. A convenient crossover scale is the Fermi temperature,For , both gases have . At fixed , the Fermi curve begins at , has zero initial slope, lies above the classical line, and approaches it at high temperature. For example its low-temperature expansion is . The plot evaluates the finite-temperature occupation integrals at fixed number density.
If the source's symbol is interpreted literally as momentum rather than a wavevector, the state-space cell is . In that convention replace in the wavevector-based coefficients and energy scales above by . The dimensionless thermodynamics and the relation are unchanged.
Release-expansion thermometry 2026-10-18
A spatial image after a chosen dark interval is a different readout from continuously recording arrivals at a fixed probe. Here is an original idealised fluorescence imaging model: release a dilute cloud, let it move freely, then apply a short uniform probe pulse. Neglect collisions, residual forces, attenuation and motion during the pulse. Use each image's cloud centre to remove mean drift. Uniform gravity shifts that centre without changing the centred variance.
For each atom let along a horizontal imaging axis. Expanding the square about the expected values gives exactlyNo Gaussian distribution is needed for this variance identity. A jointly Gaussian distribution of initial positions and velocities does make the later density Gaussian. For a separable three-dimensional cloud viewed along , integration along the viewing direction gives the column densityHere is atom number. Integrating this image along leaves a one-dimensional Gaussian distribution with standard deviation . This is a projected density measurement, not a three-dimensional radius or a direct velocity histogram. Under the uniform, optically thin probe assumptions, expected fluorescence is proportional to this column density. A known, independent point spread function with finite variance adds its own variance to the measured width; probe duration and intensity gradients need separate treatment.
For a Maxwell-Boltzmann distribution, , where is atomic mass and the Boltzmann constant. For a nonthermal distribution this defines a directional second-moment temperature, which need not describe every feature of the distribution. A covariance between position and velocity can initially narrow the cloud even though is positive. The physical covariance matrix requires .
Choose , and , giving correlation coefficient . With in milliseconds, the numerical variance in square millimetres is . Its minimum occurs at , with standard deviation . These chosen values satisfy the covariance matrix constraint; early contraction is compatible with free motion.
Generate thirteen synthetic variance observations at equally spaced times from to , adding independent Gaussian noise of standard deviation . Linear least squares fits all three coefficients in . The example yields , and . With the chosen approximate sodium mass , the true is and the fitted value is . This is a fit to generated width data, not historical observations.
If the same data are fitted to , forcing , the inferred falls to . The omitted term is partly absorbed into both fitted coefficients. In noiseless linear least squares with a free intercept, its precise effect on the quadratic coefficient iswhere the subscript denotes equal-weight moments over the chosen observation times, using one consistent unit system. Since the ratio is positive for distinct nonnegative times, negative biases this estimate downwards. An incorrect release-time origin can also generate a linear term, so fitting it does not by itself identify a physical initial covariance. Timing calibration, spatial resolution, residual magnetic fields and probe-induced recoil energy remain separate controls. Neither a cold-looking image nor a single long expansion time identifies all three coefficients.
Sub-Doppler cooling discovery thermometry 2026-10-18
The 1988 experiment loaded sodium optical molasses from a slowed atomic beam. After stopping loading, a shutter extinguished the six cooling beams. Released atoms crossed a lower probe; a photomultiplier recorded fluorescence versus arrival time. Calculations integrated the measured initial cloud distribution and an assumed three-dimensional Maxwell-Boltzmann distribution, including gravity and detection geometry. The minimum was , versus the Doppler cooling limit.[ref]
Observation of Atoms Laser Cooled below the Doppler Limit by Paul D. Lett and others. pp. 169–170, preparation, TOF apparatus, Fig. 1 and uncertainty
Probe separation, probe intensity, cloud density and shutter timing were varied. Fountain-height, horizontal-spreading and release/recapture measurements also gave below-limit results. Geometry remained an uncertainty; agreement among complementary thermometers made a detector artifact less plausible. These measured released atoms, rather than an in-situ velocity distribution.[ref]
Observation of Atoms Laser Cooled below the Doppler Limit by Paul D. Lett and others. pp. 170–171, controls, four thermometers and release limitation
Steven Chu's account connects the explanation to optical pumping, Stark effect shifts and spatially changing electromagnetic polarization; Claude Cohen-Tannoudji's lecture separates this mechanism from subrecoil schemes. Their lectures and William D. Phillips's describe complementary experimental and theoretical routes.[ref]
The manipulation of neutral particles by Steven Chu. pp. 693–695, Fig. 6
[ref]Nobel Lecture: Manipulating atoms with photons by Claude Cohen-Tannoudji. pp. 710–713, §§II–III
Thermal bead drag and heat-bath power 2026-10-24
If the relative bead velocity is independent of its separation and has a centered Maxwell-Boltzmann distribution with , averaging the drag dissipation of a bead-spring pair gives the displayed gross loss. The cross term vanishes by the zero conditional mean. At zero velocity gradient, this remains , although mechanical stress power is zero. A thermal bath replenishes microscopic kinetic energy at equilibrium. Equating this gross loss to polymer stress power omits heat-bath input and cannot follow from the independently specified Maxwellian velocity distribution. The overdamped free-energy balance of a Hookean dumbbell uses a probability current velocity instead.
Thermal speed 2026-10-24
A characteristic random particle speed. Here temperature is measured in energy units, and the factor-two convention makes a Maxwell-Boltzmann distribution proportional to . Other conventions, such as the one-dimensional standard deviation , differ by fixed factors.

