Boltzmann Statistics, Fermi-Dirac Statistics, and Bose-Einstein Statistics
Simplify the partition function by writing the total energy of a system as a sum of individual energies. The final equations depend on whether the individual particles of the system are fermions or bosons. They obey two different types of statistics, Fermi-Dirac or Bose-Einstein. Both under certain conditions can be reduced to Boltzmann statistics.
4-1 The Special Case of Boltzmann Statistics
All known particles fall into two clasees: those whose wave function must be symmetric under the operation of the interchange of two identical particles, and those whose wave function must be antisymmetric.
<p>
</p>
There are problems in which the Hamiltonian can be written as a sum of simpler Hamiltonians. The most obvious example is the case of a dilute gas. Another example is the decomposition of the Hamiltonian of a polyatomic molecule into various degrees of freedom.
<center>
<br>
<tex>H \approx H_
+ H_
+ H_
+ H_
</tex>
<br>
</center>
Canonical Partition Function
Denote individual energy states by \left
, where the superscript denotes the particle and the subscript denotes the state.
<center>
<br>
<tex>Q(N, V, T) = \sum_j e^{-E_j / k T}</tex>
<br>
<tex>Q(N, V, T) = \sum_
e^
</tex>
<br>
<tex>Q(N, V, T) = \sum_
e^
\sum_
e^
\sum_
e^
</tex>
<br>
<tex>Q(N, V, T) = q_a q_b q_c</tex>
<br>
<tex>q(V, T) = \sum_i e^{- \epsilon_i / kT}</tex>
<br>
<tex> \mbox
</tex>
<br>
<tex>Q(N, V, T) = q(V, T)^N</tex>
</center>
The last expression requires knowledge only of the energy values of an individual particle or quasi-particle. It is called a molecular partition function.
<p>
</p>
In a perfect crystal each atom is confined to one and only one lattice point. The particles are distinguishable.
Molecular Hamiltonian
Reduce an <tex>N</tex>-body problem further into individual degrees of freedom.
<center>
<br>
<tex>q_
= q_
q_
q_
q_
...</tex>
<br>
<tex>q_
= \sum_i e^{- \epsilon_i^
/ kT }</tex>
<br>
</center>
An excellent approximation in most cases is obtained by carrying out the summation below in an unrestricted manner and dividing by <tex>N!</tex>. The number of molecular states available to a molecule at room temperature, say, is much greater than the number of molecules in the system for all but the most extreme densities. The indistinguishability of the particles has been included by dividing by <tex>N!</tex>. The condition is favored by large mass, high temperature, and low density. The results are valid for polyatomic molecules as well, since the translational energy states account for almost all of the nergy states available to any molecule.
<center>
<br>
<tex>Q(N, V, T) = \sum_
e^{-(\epsilon_i + \epsilon_j + \epsilon_k + \epsilon_l + ...)/kT}</tex>
<br>
<tex> \mbox
</tex>
<br>
<tex>Q(N, V, T) = \frac
</tex>
<br>
<tex>q(V, T) = \sum_j e^{-\epsilon_j / kT } </tex>
<br>
</center>
Requirement of Boltzmann Statistics
The number of available molecular states is much greater than the number of particles in the system. The approximation becomes increasingly better at higher temperatures.
Probability of Molecules in _j_th Energy State
The total energy of the <tex>N</tex>-body system is below.
<center>
<br>
<tex>E = N \bar \epsilon</tex>
<br>
<tex>E = kT^2 \left ( \frac
\right )_
</tex>
<br>
<tex>E = N \sum_j \epsilon_j \frac{ e^{- \epsilon_j / k T}}
</tex>
<br>
<tex>\bar \epsilon = \sum_j \epsilon_j \frac{ e^{- \epsilon_j / k T}}
</tex>
<br>
<tex>\mbox
</tex>
<br>
<tex>\pi_j = \frac{e^{- \epsilon_j / k T}}{\sum_j e^{- \epsilon_j / k T}}</tex>
<br>
<tex>\pi_j = \frac{e^{- \epsilon_j / k T}}
</tex>
<br>
</center>
The last expression is the probability that a molecule is in the _j_th vibrational state irrespective of the other degrees of freedom.
Consider the effect of symmetry requirements of N-body wave functions on the sum over states.
4-2 Fermi-Dirac and Bose-Einstein Statistics
The distribution function in the case of fermions is called Fermi-Dirac statistics, and the the distribution function in the case of bosons is called Bose-Einstein statistics. These are the only exact distributions. In the case of high temperature and/or low density, both of these go over into Boltzmann or classical distribution.
<p>
</p>
Treat the general case by means of the grand canonical ensemble. Let <tex>E_j(N, V)</tex> be the energy states available to a system containing <tex>N</tex> molecules. Let <tex>\epsilon_k</tex> be the molecular quantum states. Let <tex>n_k = n_k (E_j)</tex> be the number of molecules in the <tex>k</tex>th molecular state when the system itself is in the quantum state with energy <tex>E_j</tex>. A quantum state of the entire system is specified by the set <tex>
</tex>. The energy of the system is expressed below.
<center>
<br>
<tex>E_j = \sum_k \epsilon_k n_k</tex>
<br>
<tex>N = \sum_k n_k</tex>
<br>
<tex>Q(N, V, T) = \sum_j e^{- \beta E_j}</tex>
<br>
<tex>Q(N, V, T) = \sum_n_k
e^{- \beta \sum_i* \epsilon_i n_i}</tex>
<br>
</center>
Avoid the restriction that the sum of number of particles is equal to the total number of particles by using the grand canonical partition function.
<center>
<br>
<tex>\Theta (V, T, \mu) = \sum_
^
e^
Q(N, V, T)</tex>
<br>
<tex>\Theta (V, T, \mu) = \sum_
^
\lambda^N \sum_n_k
* e^{- \beta \sum_i* \epsilon_i n_i}</tex>
<br>
<tex>\Theta (V, T, \mu) = \sum_
^
\sum_n_k
* \lambda^
e^{- \beta \sum_i* \epsilon_j n_j}</tex>
<br>
<tex>\Theta (V, T, \mu) = \sum_
^
\sum_n_k
* \Pi_k \left ( \lambda e^{- \beta \epsilon_k} \right )^
</tex>
<br>
</center>
The next step is important and is allowed since we are summing over all values of <tex>N</tex> or since we are using the grand canonical partition function.
<center>
<br>
<tex>\Theta (V, T, \mu) = \sum_
^{n_1 \mbox{max}} \sum_
{n_2 \mbox{max}} ... \Pi_k \left ( \lambda e{-\beta \epsilon_k} \right )^
</tex>
<br>
<tex> \Theta (V, T, \mu) = \Pi_k \sum_
^{n_k \mbox
\left ( \lambda e^{-\beta \epsilon_k} \right )^
</tex>
<br>
<tex>\mbox
</tex>
<br>
<tex>\Theta_
= \Pi_k \left (1 + \lambda e^{- \beta \epsilon_k} \right )</tex>
<br>
<tex>\mbox
</tex>
<br>
<tex>\Theta_
= \Pi_k \sum_
^
\left (\lambda e^{- \beta \epsilon_k} \right )^
</tex>
<br>
<tex>\Theta_
= \Pi_k \sum_
(1 - \lambda e^{- \beta \epsilon_k} \right )^{-1}</tex>
<br>
</center>
Below is a derivation of the average number of particles in the <tex>k</tex>th quantum state.
<center>
<br>
<tex>\bar N = N</tex>
<br>
<tex>N = \sum_l \bar n_k</tex>
<br>
<tex>N = kT \left ( \frac
\right )_
</tex>
<br>
<tex>N = \lambda \left( \frac
\right )_
</tex>
<br>
<tex>N = \sum_k \frac{\lambda e^{\beta \epsilon_k}}{1 \pm \lambda e^{ \beta \epsilon_k}}</tex>
<br>
<tex>\bar n_k= \frac{ \lambda e^{\beta \epsilon_k}}{1 \pm \lambda e^{ \beta \epsilon_k}}</tex>
<br>
<tex>\bar E = N \bar \epsilon</tex>
<br>
<tex>\bar E = \sum_k \bar n_k \epsilon_k</tex>
<br>
<tex>\bar E = \sum_k \frac{\lambda \epsilon_k e^{- \beta \epsilon_k}}{1 \pm \lambda e^{- \beta \epsilon_k}}</tex>
<br>
</center>
Lastly, consider the following.
<center>
<br>
<tex>pV = kT \ln \Theta (V, T, \mu)</tex>
<br>
<tex>pV = \pm kT \sum_k \ln \left [1 \pm \lambda e^{-\beta \epsilon_k} \right]</tex>
<br>
</center>
The individual particles of the system are not independent becase of the symmetry requirements of the wave function.
<p>
</p>
Both statistics go over into Boltzmann
Both kinds of statistics should go over into Boltzmann or classical statistics in the limit of high temperature or low density, where the number of available quantum states is much greater than the number of particles. The average number of molecules in any state is very small.
<center>
<br>
<tex>\lambda \mbox
</tex>
<br>
<tex>\bar n_k = \lambda e^{- \beta \epsilon_k}</tex>
<br>
<tex>\frac
= \frac{e^-
}
</tex>
<br>
<tex>q = \sum_j e^
</tex>
<br>
<tex>\bar E \right \sum_j \lambda \epsilon_j e^
</tex>
<br>
<tex>\bar \epsilon = \frac
\right \frac{\sum_j \epsilon_j e^{\beta \epsilon_j}}{\sum_j e^{\beta \epsilon_j}}</tex>
<br>
<tex>pV \right (\pm kT) \left (\pm \lambda \sum_j e^{- \beta \epsilon_j} \right )</tex>
<br>
<tex>pV = \lambda kT \sum_j e^{- \beta \epsilon_j}</tex>
<br>
<tex>pV = \lambda kT q</tex>
<br>
<tex>\beta p V = \ln \Theta</tex>
<br>
<tex>\beta p V = \lambda q</tex>
<br>
<tex>\lambda q = N</tex>
<br>
</center>
The penultimate equation above is the perfect gas law, and thus the formulas of Fermi-Dirac and Bose-Einstein statistics reduce to those of Boltzmann statistics in the classical limits.
<center>
<br>
<tex>\Theta = e^
</tex>
<br>
<tex>\Theta = \sum_
^
\frac
</tex>
<br>
<tex>Q(N, V, T) = \frac
</tex>
<br>
</center>
Apply the limit of Boltzmann statistics to the simplest system, namely a monoatomic ideal gas.