Population can move between the levels by electron-impact excitation and de-excitation, spontaneous emission, radiative absorption and stimulated emission. Other collision partners and radiative cascades can matter in a fuller atom. In the elementary coronal approximation, the plasma is dilute, optically thin, steady, and weakly illuminated at the transition frequency. Electron impacts predominantly excite the lower level, while spontaneous emission predominantly empties the upper level; collisional de-excitation and photon trapping are negligible. The population balance is
A key condition is , together with insignificant external pumping and competing channels. A Maxwellian electron Maxwell-Boltzmann distribution lets the tabulated rate coefficient be expressed using one temperature.
Landau damping is the collisionless decay of a plasma wave through resonant energy exchange with particles whose velocity is close to the wave's phase velocity . Particles slightly slower than the wave tend to gain energy from it; particles slightly faster can give energy to it. For the decreasing resonant part of a Maxwell-Boltzmann distribution, the net transfer is from the wave to the electrons, so the wave amplitude decays.
In the linear Vlasov equation description this is collisionless phase mixing, not a requirement for frequent particle collisions. Fine-scale structure develops in the velocity dependence of the perturbation even while the macroscopic electric field decreases. The resonant contribution to the growth rate has the sign of at a positive wave phase velocity: a negative slope produces damping, whereas a sufficiently positive slope can drive an instability. A population comparison must be local to the resonance; merely counting all particles faster and slower than the wave is not the criterion.
For the shifted beam Maxwell-Boltzmann distribution,
On the rising flank put . The positive slope is proportional to , whose derivative is . Its unique positive maximum is at , giving
This is the region where beam drive can most readily overcome the negative plasma slope when that background changes little over the narrow beam width. It is an estimate for the first resonant unstable wave, not a proof that every beam-plasma dispersion relation has its exact marginal phase velocity at this value. Both the total slope and the real dispersion relation determine the exact marginal mode.
For a single-particle state of energy , put . The equilibrium mean occupancies are
The 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 is
Here , for one nondegenerate branch, and for the spin-half case below. The integral uses . For classical indistinguishable particles, , giving
The 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 gives
These 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.
Figure 1.
Pressure of classical and spin-half Fermi gases with linear dispersion
. Original fixed-density calculation. The Fermi curve is obtained by solving the number equation for the chemical potential at each temperature; its positive zero-temperature pressure and the shared classical limit are marked. Details.
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.
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 exactly
No 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 density
Here 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.
Figure 1.
Cloud expansion with an initial position–velocity covariance
. Original synthetic variance observations and ideal column-density marginals. All parameters and noise are chosen here. The three-dimensional Gaussian is integrated along the viewing axis and then along the image's other axis. The profiles omit imaging blur. Details.
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 is
where 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.
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]
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]
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][ref]
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.