4.3 Fermi-Dirac statistics

Fermi-Dirac statistics describe how fermions (particles like electrons, protons, and neutrons) distribute themselves across energy states in quantum systems. This framework explains why electrons fill energy levels the way they do in atoms and solids, and it's built directly on the Pauli exclusion principle. Mastering this topic is essential for understanding the electronic properties of metals, semiconductors, and even astrophysical objects like white dwarf stars.

Fundamentals of Fermi-Dirac statistics

Fermi-Dirac (FD) statistics govern the behavior of fermions in quantum mechanical systems. Where classical statistics treat particles as distinguishable and place no limits on state occupancy, FD statistics account for the fact that identical fermions cannot share a quantum state. This distinction has enormous consequences across condensed matter physics and astrophysics.

Fermions and the exclusion principle

Fermions are particles with half-integer spin ($1/2, 3/2, \ldots$). Electrons, protons, and neutrons are all fermions. The defining constraint on fermions is the Pauli exclusion principle: no two identical fermions can occupy the same quantum state simultaneously.

• This forces electrons to stack into progressively higher energy levels rather than all piling into the ground state.
• It directly explains electron shell structure in atoms and band structure in solids.
• Contrast this with bosons (integer spin), which can pile into the same state, leading to completely different collective behavior like Bose-Einstein condensation.

Derivation of the Fermi-Dirac distribution

The FD distribution is derived from the grand canonical ensemble. The key step is recognizing that each single-particle state can be occupied by either 0 or 1 fermion (never more), which constrains the partition function.

1. Write the grand partition function for a single state with energy $E$: since occupancy is restricted to 0 or 1, the sum has only two terms.
2. Compute the average occupation number $\langle n \rangle$ by taking the appropriate thermodynamic derivative.
3. The result is the Fermi-Dirac distribution function:

$f(E) = \frac{1}{e^{(E - \mu)/k_BT} + 1}$

Here $E$ is the energy of the state, $\mu$ is the chemical potential, $k_B$ is Boltzmann's constant, and $T$ is temperature. The function $f(E)$ gives the probability that a state at energy $E$ is occupied at temperature $T$.

Notice the "+1" in the denominator. That's what distinguishes FD from Bose-Einstein statistics (which has "$-1$") and from the classical Maxwell-Boltzmann limit (which has no correction term at all).

Fermi energy and Fermi level

The Fermi energy $E_F$ is the energy of the highest occupied single-particle state at absolute zero ($T = 0$). At $T = 0$, the distribution function becomes a sharp step: all states below $E_F$ are filled ($f = 1$), and all states above are empty ($f = 0$).

• The Fermi level is the chemical potential evaluated at $T = 0$, and it marks the boundary between occupied and unoccupied states.
• You calculate $E_F$ by integrating the density of states up to the total number of particles in the system.
• This single energy scale controls a huge range of material properties: thermal conductivity, electrical conductivity, heat capacity, and magnetic susceptibility all depend on where $E_F$ sits relative to the band structure.

Properties of Fermi-Dirac systems

Density of states

The density of states $g(E)$ counts how many quantum states are available per unit energy interval. It depends on the dimensionality of the system:

• 3D free electron gas: $g(E) \propto E^{1/2}$
• 2D system: $g(E)$ is constant (independent of energy)
• 1D system: $g(E) \propto E^{-1/2}$

The density of states matters because physical observables are computed by integrating $g(E) \cdot f(E)$ over energy. A large density of states near $E_F$ means more electrons participate in thermal and transport processes, directly influencing specific heat and electrical conductivity. Experimentally, $g(E)$ can be measured using techniques like scanning tunneling spectroscopy.

Chemical potential at low temperatures

At temperatures well below the Fermi temperature $T_F = E_F / k_B$, the chemical potential stays close to $E_F$ but shifts slightly downward:

$\mu(T) \approx E_F\left[1 - \frac{\pi^2}{12}\left(\frac{k_BT}{E_F}\right)^2\right]$

This correction is small for typical metals at room temperature because $T_F$ is on the order of tens of thousands of kelvin (e.g., $T_F \approx 80{,}000$ K for copper). The chemical potential determines which states near the Fermi level are thermally accessible, which is why it's central to thermoelectric effects like the Seebeck effect.

Temperature dependence of Fermi energy

For most practical purposes, the Fermi energy is nearly temperature-independent. The correction follows the same form as the chemical potential:

$E_F(T) \approx E_F(0)\left[1 - \frac{\pi^2}{12}\left(\frac{k_BT}{E_F(0)}\right)^2\right]$

Since $k_BT / E_F$ is typically very small (on the order of $10^{-2}$ at room temperature for metals), this shift is negligible in most calculations. Still, it becomes relevant when computing thermodynamic quantities like the electronic specific heat, which depends sensitively on how the distribution broadens near $E_F$.

Applications of Fermi-Dirac statistics

Electron gas in metals

The simplest model of a metal treats conduction electrons as a non-interacting Fermi gas confined to the volume of the solid. Despite its simplicity, this model explains several key observations:

• Electronic specific heat scales linearly with temperature ($C_e \propto T$), in contrast to the classical prediction of a constant value. Only electrons within $\sim k_BT$ of the Fermi level contribute.
• Pauli paramagnetism gives a temperature-independent magnetic susceptibility, because only a thin shell of electrons near $E_F$ can flip their spins in response to a field.
• The model predicts a $T^2$ contribution to electrical resistivity from electron-electron scattering at low temperatures.
• It forms the starting point for Fermi liquid theory, which incorporates interactions perturbatively.

Semiconductor physics

FD statistics are indispensable in semiconductor physics. The distribution function determines how electrons populate the conduction band and how holes populate the valence band.

• The bandgap separates the valence and conduction bands; FD statistics predict how carrier concentrations in each band vary with temperature.
• In doped semiconductors, the Fermi level shifts toward the conduction band (n-type) or valence band (p-type), and its exact position is calculated using FD statistics combined with the density of states.
• Device physics for p-n junctions, transistors, and solar cells all rely on knowing the carrier distribution at equilibrium and under bias.

White dwarf stars

White dwarf stars provide a dramatic astrophysical application of FD statistics. After a low-to-intermediate mass star exhausts its nuclear fuel, the remaining core is supported against gravitational collapse by electron degeneracy pressure.

• At the extreme densities inside a white dwarf, electrons are packed so tightly that the Pauli exclusion principle generates an outward pressure even at zero temperature.
• This degeneracy pressure sets a maximum mass for stable white dwarfs: the Chandrasekhar limit ($\approx 1.4 \, M_\odot$). Beyond this mass, electron degeneracy pressure cannot resist gravity, and the star collapses further.
• FD statistics also predict the mass-radius relationship and cooling rate of white dwarfs.

Quantum effects in Fermi systems

Pauli paramagnetism

When you apply an external magnetic field to a metal, electrons with spins aligned opposite to the field can flip to align with it, lowering their energy. But only electrons near the Fermi level have empty states available to flip into.

• The resulting paramagnetic susceptibility is temperature-independent and proportional to the density of states at $E_F$.
• This is much smaller than the classical prediction (which would have all electrons contributing), precisely because the Pauli exclusion principle "freezes out" most of the electron population.
• In heavy fermion compounds, the large effective mass enhances the density of states at $E_F$, producing an anomalously large Pauli susceptibility.

Landau diamagnetism

The orbital motion of electrons in a magnetic field produces a diamagnetic response (opposing the applied field). This is distinct from the spin-based Pauli paramagnetism.

• For free electrons, Landau diamagnetism is exactly $-1/3$ the magnitude of the Pauli paramagnetic susceptibility, so the net response is paramagnetic.
• The magnitude is inversely proportional to the effective mass, making it particularly important in systems with light carriers like graphene.
• Both Pauli and Landau contributions must be summed to get the total magnetic susceptibility of a metal.

Quantum oscillations

When a strong magnetic field is applied to a metal, the continuous energy spectrum splits into discrete Landau levels. As the field strength changes, these levels sweep through the Fermi energy, causing periodic oscillations in measurable quantities.

• de Haas-van Alphen effect: oscillations in magnetization
• Shubnikov-de Haas effect: oscillations in resistivity

The oscillation frequency is directly related to the extremal cross-sectional area of the Fermi surface perpendicular to the field. This makes quantum oscillations one of the most powerful tools for mapping Fermi surface geometry, measuring effective masses, and probing the electronic structure of metals, semimetals, and topological materials.

Fermi-Dirac vs Bose-Einstein statistics

Key differences and similarities

Transition temperature

Quantum statistical effects become significant below a characteristic temperature scale.

• For fermions, this is the Fermi temperature: $T_F = E_F / k_B$. In metals, $T_F \sim 10^4$ to $10^5$ K, so electrons are deep in the quantum regime even at room temperature.
• For bosons, the critical temperature marks the onset of Bose-Einstein condensation: $T_c \propto n^{2/3}/m$, where $n$ is the number density and $m$ is the particle mass.
• Below these temperatures, classical approximations fail and the full quantum distributions must be used.

Quantum degeneracy

A system becomes quantum degenerate when the thermal de Broglie wavelength $\lambda$ becomes comparable to the average interparticle spacing. The degeneracy parameter $n\lambda^3$ quantifies this:

• $n\lambda^3 \ll 1$: classical regime, Maxwell-Boltzmann statistics apply
• $n\lambda^3 \gtrsim 1$: quantum degenerate regime

For fermions, degeneracy produces a filled Fermi sea with degeneracy pressure. For bosons, it triggers condensation into the ground state. The thermal de Broglie wavelength is $\lambda = \sqrt{2\pi\hbar^2 / mk_BT}$, so lower temperatures and lighter particles push systems toward degeneracy.

Experimental observations

Photoemission spectroscopy

Photoemission spectroscopy directly measures the energy distribution of electrons in a material. A photon ejects an electron, and measuring its kinetic energy reveals the binding energy of the state it came from.

• The occupied portion of the band structure is mapped this way, and the sharp cutoff at $E_F$ is a direct signature of the Fermi-Dirac distribution.
• Angle-resolved photoemission spectroscopy (ARPES) adds momentum resolution, allowing full mapping of the band structure $E(\mathbf{k})$.
• ARPES has been instrumental in studying superconductors, topological insulators, and strongly correlated materials where electron interactions reshape the band structure.

de Haas-van Alphen effect

This effect manifests as oscillations in the magnetic susceptibility of a metal as the applied magnetic field is varied. The physical origin is the passage of Landau levels through the Fermi energy.

• The oscillation period in $1/B$ is related to the extremal Fermi surface cross-section $A$ by: the frequency $F = (\hbar / 2\pi e) A$.
• By measuring oscillation frequencies for different field orientations, you can reconstruct the full 3D Fermi surface topology.
• The temperature dependence of the oscillation amplitude yields the effective mass, and the field dependence gives information about scattering rates.

Quantum Hall effect

In a two-dimensional electron system subjected to a strong perpendicular magnetic field, the Hall conductance becomes quantized in integer multiples of $e^2/h$.

• This occurs when Landau levels are completely filled, and the Fermi level sits in a gap between levels. The system is incompressible, and the Hall resistance shows exact plateaus.
• The integer quantum Hall effect is well explained by single-particle FD statistics combined with Landau quantization and disorder-induced localization.
• The fractional quantum Hall effect requires many-body physics beyond single-particle FD statistics, involving strongly correlated electron states with fractional charge excitations.
• These discoveries opened the field of topological phases of matter.

Advanced topics

Fermi liquid theory

Real metals have electron-electron interactions, yet many of their properties look qualitatively similar to the free Fermi gas. Fermi liquid theory (developed by Landau) explains why: interactions dress the bare electrons into quasiparticles that behave like free fermions but with renormalized parameters (effective mass, magnetic moment, etc.).

• The Fermi surface survives interactions, and low-energy excitations near it are long-lived quasiparticles.
• Thermodynamic and transport properties retain the same functional forms as the free gas (linear specific heat, $T^2$ resistivity) but with renormalized coefficients.
• The theory breaks down in strongly correlated systems where interactions are too strong for the quasiparticle picture to hold.

Strongly correlated electron systems

In some materials, electron-electron interactions are comparable to or larger than the kinetic energy, and the independent-particle picture fails entirely.

• Examples include transition metal oxides, heavy fermion compounds, and high-$T_c$ cuprate superconductors.
• Mott insulators are a striking case: band theory predicts a metal, but strong Coulomb repulsion localizes the electrons and opens a gap.
• Theoretical treatment requires advanced methods like dynamical mean-field theory (DMFT) and quantum Monte Carlo simulations.
• These systems host exotic phases including unconventional superconductivity, quantum spin liquids, and quantum critical points.

Topological insulators

Topological insulators have an insulating bulk but host robust conducting states on their surfaces or edges. These surface states are protected by time-reversal symmetry and the topology of the bulk band structure.

• The surface electrons behave as massless Dirac fermions with spin-momentum locking: the spin direction is tied to the momentum direction, suppressing backscattering.
• They arise from strong spin-orbit coupling combined with band inversion.
• The surface states exhibit a quantized Hall conductance without an external magnetic field (in certain configurations).
• Potential applications include spintronics, low-dissipation electronics, and platforms for topological quantum computation.