Warning: mkdir() [
function.mkdir]: Permission denied in
/home/webs/affiliatelib2/CacheManager.php on line
12
Warning: mkdir() [
function.mkdir]: No such file or directory in
/home/webs/affiliatelib2/CacheManager.php on line
12
Warning: fopen(/home/templatecore2cache//*cluesnet.com/af/af6d03e14a902294f91112655d9554f7f68f7ed3.tc2cache) [
function.fopen]: failed to open stream: No such file or directory in
/home/webs/affiliatelib2/CacheManager.php on line
130
Warning: fwrite(): supplied argument is not a valid stream resource in
/home/webs/affiliatelib2/CacheManager.php on line
131
Warning: fclose(): supplied argument is not a valid stream resource in
/home/webs/affiliatelib2/CacheManager.php on line
132
In
statistical mechanics,
Fermi-Dirac statistics is a particular case of particle statistics developed by Enrico Fermi and Paul Dirac that determines the statistical distribution of
fermions over the energy states for a system in thermal equilibrium. In other words, it is a probability of a given energy level to be occupied by a fermion.
More generally,
Fermi-Dirac statistics means that the
Slater determinant of
fermions must be
Skew-symmetric matrix under an exchange of every pair of
fermions (that is, if one exchanges any fermion with another, the
wavefunction gets an overall minus sign).
Fermions are particles which are identical particles and obey the
Pauli exclusion principle, i.e., no more than one particle may occupy the same quantum state at the same time. Statistical thermodynamics is used to describe the behaviour of large numbers of particles. A collection of non-interacting fermions is called a Fermi gas.
F-D statistics was introduced in 1926 by
Enrico Fermi and
Paul Dirac and applied in 1926 by
Ralph Fowler to describe the collapse of a star to a white dwarf and in 1927 by Arnold Sommerfeld to electrons in metals.
Pascual Jordan developed in 1925 the same statistics which he called
Pauli statistics. The problem was that his referee Max Born forgot the paper for six months before finding it again. In the meantime it was independently discovered by Enrico Fermi and Paul Dirac.http://scienceweek.com/2000/sw000519.htm
For F-D statistics, the expected number of particles in states with energy \epsilon _i is
n_i = \frac{g_i}{e^{(\epsilon_i-\mu) / k T} + 1}
where:
n_i \ is the number of particles in state
i,
\epsilon_i \ is the energy of state
i,
g_i \ is the
Degenerate energy level of state
i (the number of states with energy \epsilon_i \ ),
\mu\ is the
chemical potential (Sometimes the Fermi energy E_F \ is used instead, as a low-temperature approximation),
\ k\ is Boltzmann's constant, and
\ T\ is absolute
temperature.
In the case where \mu is the Fermi energy E_F \ and g_i = 1 \ , the function is called the
Fermi function:
F(E) = \frac{1}{e^{(\epsilon_i-E_F) / k T} + 1}
Which distribution to use
A derivation
Consider a single-particle state of a multiparticle system, whose energy is \mathbf{\epsilon}. For example, if our system is some quantum gas in a box, then a state might be a particular single-particle wave function. Recall that, for a grand canonical ensemble in general, the partition function (statistical mechanics) is
:Z \;= \sum_s e^{ -( E(s) - \mu N(s) ) / kT}
where
E(s) is the energy of a state
s,
N(s) is the number of particles possessed by the system when in the state
s,
\mu denotes the chemical potential, and
s is an index that runs through all possible microstates of the system.
In the present context, we take our system to be a fixed single-particle state (
not a particle). So our system has energy n \cdot \epsilon when the state is occupied by
n particles, and
0 if it is unoccupied. Consider the balance of single-particle states to be the
reservoir. Since the system and the reservoir occupy the same physical space, there is clearly exchange of particles between the two (indeed, this is the very phenomenon we are investigating). This is why we use the grand partition function, which, via chemical potential, takes into consideration the flow of particles between a system and its thermal reservoir.
For fermions, a state can only be either occupied by a single particle or unoccupied. Therefore our system has multiplicity two: occupied by one particle, or unoccupied, called s_1 and s_2 respectively. We see that E(s_1) = \; \epsilon, N(s_1) = \; 1, and E(s_2) = \; 0, N(s_2) = \; 0. The partition function is therefore
Z = \sum_{i = 1} ^2 e^{ -( E(s_i) - \mu N(s_i) ) / kT}
= e^{ -( \epsilon - \mu ) / kT} + 1 .
For a grand canonical ensemble, probability of a system being in the microstate s_{\alpha} is given by
P( s_{\alpha} ) = \frac{e^ {-( E(s_{\alpha}) - \mu N(s_{\alpha}) ) / kT} }{Z}.
Our state being occupied by a particle means the system is in microstate s_1, whose probability is
\bar{n} = P( s_1 ) =
\frac{ e^{ -( E(s_1) - \mu N(s_1) ) / kT} }{Z}= \frac{e^{ -( \epsilon - \mu ) / kT-->{e^{ -( \epsilon - \mu)/ kT} + 1}= \frac{1}{e^{ ( \epsilon - \mu)/ kT} + 1}.
\bar{n} is called the
Fermi-Dirac distribution. For a fixed temperature
T, \bar{n}(\epsilon) is the probability that a state with energy ε will be occupied by a fermion. Notice\bar{n} is a decreasing function in ε. This is consistent with our expectation that higher energy states are less likely to be occupied.
Note that if the energy level ε has degeneracy \; g_{\epsilon}, then we would make the simple modification:
\bar{n} = g_{\epsilon} \cdot \frac{1}{e^{ ( \epsilon - \mu)/ kT} + 1}.
This number is then the expected number of particles in the totality of the states with energy ε.
For all temperature
T, \bar{n}(\mu) = \frac{1}{2} , that is, the states whose energy is μ will always have equal probability of being occupied or unoccupied.
In the limit T \rightarrow 0, \bar{n} becomes a step function (
see graph above). All states whose energy is below the chemical potential will be occupied with probability
1 and those states with energy above μ will be unoccupied. The chemical potential at zero temperature is called Fermi energy, denoted by E _F, i.e.
E _F = \; \mu(T = 0).
It may be of interest here to note that, in general the chemical potential is temperature-dependent. However, for systems well below the
Fermi temperature T_F = \frac{ E _F }{k}, it is often sufficient to use the approximation \mathbf{\mu} ≈ \; E_F .
Another derivation
In the previous derivation, we have made use of the grand partition function (or Gibbs sum over states). Equivalently, the same result can be achieved by directly analyzing the multiplicities of the system.
Suppose 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. Assume 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, labeled by index
i , each levelhaving energy ε
i and containing a total of
ni particles. Suppose each level contains
gi 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
gi 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 whichwe will write as:
w(1,g)=\frac{g!}{1!(g-1)!}
We can distribute 2 particles in
g sublevels by putting one in the first sublevel and then distributing the remaining
n − 1 particles in the remaining
g − 1 sublevels, or we could put one in the second sublevel and then distribute the remaining
n − 1 particles in the remaining
g − 2 sublevels, etc. so that
w'(2, g
) = w
(1, g
− 1) + w
(1,g
− 2) + ... + w
(1, 1) or
w(2,g)=\sum_{k=1}^{g-1}w(1,g-k) =\sum_{k=1}^{g-1}\frac{(g-k)!}{1!(g-k-1)!}=\sum_{g-k=1}^{g-1}\frac{(g-k)!}{1!(g-k-1)!}=\frac{g!}{2!(g-2)!}
where we have used the following theorem involving
binomial coefficients:
\sum_{k=n}^g \frac{k!}{n!(k-n)!}=\frac{(g+1)!}{(n+1)!(g-n)!}
Continuing this process, we can see that
w(
n,
g) is just a binomial coefficient
w(n,g)=\frac{g!}{n!(g-n)!}
The number of ways that a set of occupation numbers
ni can be realized is the product of the ways that each individual energy level can be populated:
W = \prod_i w(n_i,g_i) = \prod_i \frac{g_i!}{n_i!(g_i-n_i)!}
Following the same procedure used in deriving the Maxwell-Boltzmann statistics,we wish to find the set of
ni for which
W is maximized, 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:
f(n_i)=\ln(W)+\alpha(N-\sum n_i)+\beta(E-\sum n_i \epsilon_i)
Again, using
Stirling's approximation for the factorials and taking the derivative with respect to
ni, and setting the result to zero and solving for
ni yields the Fermi-Dirac population numbers:
n_i = \frac{g_i}{e^{\alpha+\beta \epsilon_i}+1}
It can be shown thermodynamically that β = 1/
kT where
k is
Boltzmann's constant and
T is the
temperature, and that α = -μ/
kT where μ is the
chemical potential, so that finally:
n_i = \frac{g_i}{e^{(\epsilon_i-\mu)/kT}+1}
Note that the above formula is sometimes written:
n_i = \frac{g_i}{e^{\epsilon_i/kT}/z+1}
where z=exp(\mu/kT) is the absolute activity.
References
Charles Kittel and Herbert Kroemer,
Thermal Physics, 2nd ed. (W.H. Freeman, 1980)
See also
Fermi–Dirac statistics - Wikipedia, the free encyclopedia
In statistical mechanics, Fermi-Dirac statistics is a particular case of particle statistics developed by Enrico Fermi and Paul Dirac that determines the statistical distribution ...
Fermi-Dirac statistics definition of Fermi-Dirac statistics in the ...
Fermi-Dirac statistics, class of statistics that applies to particles called fermions. Fermions have half-integral values of the quantum mechanical property called spin and are ...
Fermi-Dirac statistics - definition of Fermi-Dirac statistics by the ...
Fermi-Dirac statistics. See under statistical mechanics. Thesaurus Legend: Synonyms Related Words Antonyms. Noun: 1. Fermi-Dirac statistics - (physics) law obeyed by a systems of ...
Derivation of Fermi-Dirac Statistics
Provides a brief technical explanation and a derivation.
Fermi-Dirac statistics
Fermi-Dirac statistics ... Photon statistics Up: Quantum statistics Previous: Formulation of the statistical
Maxwell–Boltzmann statistics - Wikipedia, the free encyclopedia
Fermi-Dirac 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 Fermi-Dirac Distribution
The Fermi-Dirac Distribution The Fermi-Dirac distribution applies to fermions, particles ... Applied statistics concepts
Fermi-Dirac Statistics
Fermi-Dirac Statistics ... Derivation of Fermi-Dirac Statistics Up: Derivation of Fermi-Dirac Statistics Previous: Derivation of Fermi-Dirac Statistics
Enrico Fermi: Fermi-Dirac Statistics
Fermi-Dirac Statistics. During his postgraduate study at the University of Gottingen, Fermi developed the distinct statistical model, based on the behavior of a perfect gas ...
Fermi Dirac Statistics And the Number Theory