advanced30 min studyLesson 35 in Physics

Bose and Fermi statistics

Derive ideal quantum occupation laws and connect exclusion, enhancement, and the dilute classical limit to mode counting.

On this page

What you’ll understand

  • Derive ideal Bose and Fermi mean occupations from their allowed mode occupations.
  • Compute single-mode number fluctuations and identify the common dilute limit.
  • State what finite size, interactions, and fixed total number change about condensation arguments.

Identical particles do not merely add a small correction to classical counting when wave packets overlap strongly. Bosons and fermions permit different occupation patterns, and those patterns reorganize equilibrium even without interactions. The distinction is easiest to derive one mode at a time. We will compare the two statistics at the same temperature, chemical potential, and mode energy, then distinguish that comparison from one made at the same total particle density. Those are different physical constraints.

The allowed occupation determines the sum

For a noninteracting single-particle mode of energy ϵ\epsilon, define r=e−β(ϵ−μ)r=e^{-\beta(\epsilon-\mu)}. A spinless fermionic mode can be empty or singly occupied, so its grand partition factor is 1+r1+r. A bosonic mode admits every nonnegative occupation, giving 1+r+r2+⋯=1/(1−r)1+r+r^2+\cdots=1/(1-r) when r<1r<1. Differentiating the logarithm with respect to the occupation weight produces the mean occupation. The sign difference in the final distributions comes directly from these allowed state counts.

A mode includes all quantum numbers needed to specify one single-particle state. Two fermions with opposite spins may share a spatial orbital because they occupy distinct spin modes; exclusion does not forbid that. The factorization into independent mode sums assumes the Hamiltonian is a sum of mode energies times occupations. Interactions generally correlate occupations and invalidate this simple product, even though the underlying bosonic or fermionic exchange symmetry remains exact.

nˉB(ϵ)=1eβ(ϵ−μ)−1,nˉF(ϵ)=1eβ(ϵ−μ)+1,r=e−β(ϵ−μ)\bar n_B(\epsilon)=\frac1{e^{\beta(\epsilon-\mu)}-1},\qquad \bar n_F(\epsilon)=\frac1{e^{\beta(\epsilon-\mu)}+1},\qquad r=e^{-\beta(\epsilon-\mu)}
These are ideal-gas grand canonical means, not definite occupation numbers in individual measurements. For ideal bosons at positive temperature, μ lies below the lowest single-particle energy in a finite normalizable grand ensemble.↔ Scroll wide equations horizontally.

Fluctuations distinguish enhancement and exclusion

A fermionic occupation is Bernoulli, with variance nˉF(1−nˉF)\bar n_F(1-\bar n_F). It vanishes as the mode becomes certainly empty or certainly full. A bosonic geometric occupation has variance nˉB(1+nˉB)\bar n_B(1+\bar n_B), which exceeds its mean. The extra term expresses enhanced occupation fluctuations in this ideal thermal mode. It is not a claim that every bosonic state has bunching: a fixed-number Fock state has zero number variance and is a different preparation.

When r≪1r\ll1, both means approach rr, and their variances approach their means. Exchange effects then become small, recovering Maxwell–Boltzmann counting. For a gas, the appropriate dilute condition involves the density of occupied states, often summarized by nλT3≪1n\lambda_T^3\ll1 with spin degeneracy handled consistently. Comparing formulas at fixed μ\mu does not keep density fixed. To compare two gases at the same total density, determine the chemical potential separately from the number constraint in each statistics.

Var⁡(nB)=nˉB(1+nˉB),Var⁡(nF)=nˉF(1−nˉF),r≪1⟹nˉB≃nˉF≃r\operatorname{Var}(n_B)=\bar n_B(1+\bar n_B),\qquad \operatorname{Var}(n_F)=\bar n_F(1-\bar n_F),\qquad r\ll1\Longrightarrow\bar n_B\simeq\bar n_F\simeq r
The variance formulas apply to independent ideal grand canonical modes. Fixed-total-number constraints and interactions introduce correlations or different fluctuations.↔ Scroll wide equations horizontally.

A gas requires a density of states

Mean occupations are not yet total particle number or energy. Sum them over the actual single-particle spectrum, or integrate against a density of states only when the level spacing is sufficiently fine. At low temperature, ideal fermions fill states up to a Fermi energy at fixed density, so appreciable kinetic energy survives even without thermal excitation. For ideal bosons in suitable dimensions and confinement, excited states can saturate, allowing a macroscopically occupied ground mode in the thermodynamic limit.

Dimensionality, dispersion, and trapping geometry determine whether that saturation argument works. A uniform ideal gas in one or two dimensions does not have ordinary finite-temperature Bose condensation in its standard thermodynamic limit. A finite box always retains discrete levels and smooth finite partition functions. Interactions can support other collective phenomena and alter fluctuations. A divergent grand canonical ground-mode mean as μ\mu approaches its energy signals the need to enforce the physical number constraint; it is not a prediction of infinitely many particles in a finite fixed-number experiment.

The same occupation weight, different statistics

A single mode satisfies β(ϵ−μ)=ln⁡2\beta(\epsilon-\mu)=\ln2. Compare its bosonic and fermionic mean occupations and variances at those same reservoir parameters.

  1. The occupation weight is r=e−ln⁡2=1/2r=e^{-\ln2}=1/2. For bosons, nˉB=r/(1−r)=1\bar n_B=r/(1-r)=1. For fermions, nˉF=r/(1+r)=1/3\bar n_F=r/(1+r)=1/3.

  2. The bosonic variance is 1(1+1)=21(1+1)=2. The fermionic variance is (1/3)(2/3)=2/9(1/3)(2/3)=2/9. These are number variances, so they are dimensionless.

  3. The bosonic distribution is PB(n)=2−(n+1)P_B(n)=2^{-(n+1)} for every n≥0n\geq0. The fermionic probabilities are PF(0)=2/3P_F(0)=2/3 and PF(1)=1/3P_F(1)=1/3, with every higher occupation forbidden for this one mode.

The result

The same reservoir weight produces different averages because the allowed many-particle states differ. This comparison does not hold total density fixed between two complete gases.

Explore the model

Interactive lab

Interactive lab · probability

Exclusion as a collection of Bernoulli modes

Treat each trial as an independent equal-energy fermionic mode with occupation probability p. Vary p toward zero and one, and observe the exact binomial mean and variance for n modes.

Exclusion as a collection of Bernoulli modesThese are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Current readings: Mean · np: 5; Variance · np(1−p): 2.5; Total probability: 1. All parameters are adjustable with the labeled controls below.03581000.070.150.220.3Successes kP(X = k)k = 0 · P = 9.77e-4
Exact probabilityHighlighted outcomeMean np
Mean · np
5
Variance · np(1−p)
2.5
Total probability
1
t = 0 s

These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThe lab illustrates fermionic occupation suppression when independent modes become nearly full. It does not show bosonic geometric distributions, unequal-energy Fermi seas, or interacting quantum gases.

Check your understanding

What must be done before comparing ideal Bose and Fermi gases at the same temperature and total particle density?

Try it yourself

For a finite ideal bosonic ground mode of energy zero, let μ=−kBTln⁡(1+1/M)\mu=-k_BT\ln(1+1/M) with M>0M>0. Find its grand canonical mean and variance, and explain why a strictly fixed total particle number changes the interpretation.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Substitute the chosen μ into the single-mode Bose formula, then use its variance identity.

Compare with the solution

Here e−βμ=1+1/Me^{-\beta\mu}=1+1/M, so nˉ0=M\bar n_0=M and Var⁡(n0)=M(1+M)\operatorname{Var}(n_0)=M(1+M). The large variance belongs to an open geometric mode. If the entire finite gas has fixed number NN, the ground mode cannot exceed NN and is correlated with excited occupations. One cannot transfer this unbounded single-mode fluctuation distribution unchanged to that constrained preparation.

Take it with you

  • Quantum occupation laws follow from allowed mode counts and the grand canonical weight.
  • Fermionic exclusion suppresses and ideal thermal bosonic statistics enhance single-mode number fluctuations.
  • Density constraints, spectral dimension, finite size, and interactions control whole-gas conclusions.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.