In statistical mechanics, FermiDirac statistics determines the statistical distribution of fermions over the energy states for a system in thermal equilibrium. Fermions are particles which are indistinguishable and obey the Pauli exclusion principle, i.e., that no two particles may occupy the same state at the same time. Statistical thermodynamics is used to describe the behaviour of large numbers of particles. A collection of noninteracting fermions is called a Fermi gas.
FermiDirac (or FD) statistics are closely related to MaxwellBoltzmann statistics and BoseEinstein statistics. While FD statistics holds for fermions, BE statistics plays the same role for bosons – the other type of particle found in nature. MB statistics describes the velocity distribution of particles in a classical gas and represents the classical (hightemperature) limit of both FD and BE statistics. MB statistics are particularly useful for studying gases, and BE statistics are particularly useful when dealing with photons and other bosons. FD statistics are most often used for the study of electrons in solids. As such, they form the basis of semiconductor device theory and electronics. The invention of quantum mechanics, when applied through FD statistics, has made advances such as the transistor possible. For this reason, FD statistics are wellknown not only to physicists, but also to electrical engineers.
FD statistics was introduced in 1926 by Enrico Fermi and Paul Dirac and applied in 1927 by Arnold Sommerfeld to electrons in metals.
The expected number of particles in an energy state i for FD statistics is
where:
 n_{i} is the number of particles in state i

g_{i} is the degeneracy of state i
 ε_{i} is the energy of state i

μ is the chemical potential. Sometimes the Fermi energy E_{F} is used instead, as a lowtemperature approximation.

k is Boltzmann's constant

T is absolute temperature

exp is the exponential function
Derivation of the FermiDirac distribution
Say there are two fermions placed in a system with four energy levels. There are six possible arrangements of such a system, which are shown in the diagram below.
ε_{1} ε_{2} ε_{3} ε_{4}
A * *
B * *
C * *
D * *
E * *
F * *
Each of these arrangements is called a microstate of the system. The ergodic hypothesis states that at thermal equilibrium, each of these microstates will be equally likely, subject to the constraints that there be a fixed total energy and a fixed number of particles.
Depending on the values of the energy for each state, it may be that total energy for some of these six combinations is the same as others. Indeed, if we assume that the energies are multiples of some fixed value ε, the energies of each of the microstates become:
A: 3ε
B: 4ε
C: 5ε
D: 5ε
E: 6ε
F: 7ε
So if we know that the system has an energy of 5ε, we can conclude that it will be equally likely that it is in state C or state D. Note that if the particles were distinguishable (the classical case), there would be twelve microstates altogether, rather than six.
Now suppose we have a number of energy levels, labelled by index i , each level having energy ε_{i} and containing a total of n_{i} particles. Suppose each level contains g_{i} distinct sublevels, all of which have the same energy, and which are distinguishable. For example, two particles may have different momenta, in which case they are distinguishable from each other, yet they can still have the same energy. The value of g_{i} associated with level i is called the "degeneracy" of that energy level. The Pauli exclusion principle states that only one fermion can occupy any such sublevel.
Let w(n,g) be the number of ways of distributing n particles among the g sublevels of an energy level. Its clear that there are g ways of putting one particle into a level with g sublevels, so that w(1,g)=g which we will write as:
We can distribute 2 particles in g sublevels by putting one in the first sublevel and then distributing the remaining n1 particles in the remaining g1 sublevels, or we could put one in the second sublevel and then distribute the remaining n1 particles in the remaining g2 sublevels, etc. so that w(2,g)=w(1,g1)+w(1,g2)+...+w(1,1) or
where we have used the following theorem involving binomial coefficients:
Continuing this process, we can see that w(n,g) is just a binomial coefficient
The number of ways that a set of occupation numbers n_i can be realized is the product of the ways that each individual energy level can be populated:
Following the same procedure used in deriving the MaxwellBoltzmann distribution, we wish to find the set of n_{i} for which W is maximised, subject to the constraint that there be a fixed number of particles, and a fixed energy. We constrain our solution using Lagrange multipliers forming the function:
Taking the derivative with respect to n_{i} and setting the result to zero and solving for n_{i} yields the FermiDirac population numbers:
It can be shown thermodynamically that β=1/kT where k is Boltzmann's constant and T is the temperature. The term containing α is variously written:
where μ is the chemical potential and z is the absolute activity.
See Also
Maxwell Boltzmann statistics (derivation)
BoseEinstein statistics
parastatistics