next up previous
Next: Harmonic Oscillators Up: Applications of Statistical Thermodynamics Previous: General Paramagnetism


Equipartition Theorem

The internal energy of a monatomic ideal gas containing $ N$ particles is $ (3/2) N
 k T$ . This means that each particle possess, on average, $ (3/2) k T$ units of energy. Monatomic particles have only three translational degrees of freedom, corresponding to their motion in three dimensions. They possess no internal rotational or vibrational degrees of freedom. Thus, the mean energy per degree of freedom in a monatomic ideal gas is $ (1/2) k T$ . In fact, this is a special case of a rather general result. Let us now try to prove this.

Suppose that the energy of a system is determined by $ f$ generalized coordinates, $ q_k$ , and $ f$ corresponding generalized momenta, $ p_k$ , so that

$\displaystyle E = E(q_1, \cdots, q_f, p_1,\cdots, p_f).$ (7.124)

Suppose further that:
  1. The total energy splits additively into the form

    $\displaystyle E = \epsilon_i(p_i) + E'(q_1,\cdots, p_f),$ (7.125)

    where $ \epsilon_i$ involves only one variable, $ p_i$ , and the remaining part, $ E'$ , does not depend on $ p_i$ .
  2. The function $ \epsilon_i$ is quadratic in $ p_i$ , so that

    $\displaystyle \epsilon_i(p_i) = b p_i^{ 2},$ (7.126)

    where $ b$ is a constant.
The most common situation in which the previous assumptions are valid is where $ p_i$ is a momentum. This is because the kinetic energy is usually a quadratic function of each momentum component, whereas the potential energy does not involve the momenta at all. However, if a coordinate, $ q_i$ , were to satisfy assumptions 1 and 2 then the theorem that we are about to establish would hold just as well.

What is the mean value of $ \epsilon_i$ in thermal equilibrium if conditions 1 and 2 are satisfied? If the system is in equilibrium at absolute temperature $ T\equiv (k \beta)^{-1}$ then it is distributed according to the canonical probability distribution. In the classical approximation, the mean value of $ \epsilon_i$ is expressed in terms of integrals over all phase-space:

$\displaystyle \overline{\epsilon_i} = \frac{ \int_{-\infty}^{\infty} \exp[-\bet...
... {\int_{-\infty}^{\infty} \exp[-\beta  E(q_1,\cdots, p_f)]  dq_1\cdots dp_f}.$ (7.127)

Condition 1 gives

$\displaystyle \overline{\epsilon}_i$ $\displaystyle =\frac{ \int_{-\infty}^{\infty} \exp[-\beta (\epsilon_i + E')] ...
...p_f} {\int_{-\infty}^{\infty} \exp[-\beta (\epsilon_i + E')] dq_1\cdots dp_f}$    
  $\displaystyle = \frac{\int_{-\infty}^{\infty} \exp(-\beta  \epsilon_i) \epsil...
...psilon_i)  dp_i  \int_{-\infty}^{\infty} \exp(-\beta  E') dq_1\cdots dp_f},$ (7.128)

where use has been made of the multiplicative property of the exponential function, and where the final integrals in both the numerator and denominator extend over all variables, $ q_k$ and $ p_k$ , except for $ p_i$ . These integrals are equal and, thus, cancel. Hence,

$\displaystyle \overline{\epsilon}_i = \frac{\int_{-\infty}^{\infty} \exp(-\beta...
... \epsilon_i  dp_i} {\int_{-\infty}^{\infty}\exp(-\beta  \epsilon_i)  dp_i}.$ (7.129)

This expression can be simplified further because

$\displaystyle \int_{-\infty}^\infty \exp(-\beta  \epsilon_i) \epsilon_i  dp_...
...ial \beta} \left[\int_{-\infty}^\infty \exp(-\beta  \epsilon_i)  dp_i\right],$ (7.130)

so

$\displaystyle \overline{\epsilon}_i = - \frac{\partial}{\partial\beta} \ln\left[\int_{-\infty}^\infty \exp(-\beta  \epsilon_i)  dp_i\right].$ (7.131)

According to condition 2,

$\displaystyle \int_{-\infty}^{\infty} \exp(-\beta  \epsilon_i)  dp_i = \int_{...
...rac{1}{\sqrt{\beta}} \int_{-\infty}^{\infty} \exp\left(- b  y^{ 2}\right) dy,$ (7.132)

where $ y = \sqrt{\beta}  p_i$ . Thus,

$\displaystyle \ln \int_{-\infty}^{\infty} \exp(-\beta  \epsilon_i)  dp_i = - ...
...}{2} \ln \beta + \ln \int_{-\infty}^{\infty} \exp\left(- b  y^{ 2}\right) dy.$ (7.133)

Note that the integral on the right-hand side is independent of $ \beta$ . It follows from Equation (7.131) that

$\displaystyle \overline{\epsilon}_i = - \frac{\partial}{\partial \beta} \left( - \frac{1}{2} \ln \beta\right) = \frac{1}{2 \beta},$ (7.134)

giving

$\displaystyle \overline{\epsilon}_i = \frac{1}{2}  k  T.$ (7.135)

This is the famous equipartition theorem of classical physics. It states that the mean value of every independent quadratic term in the energy is equal to $ (1/2) k T$ . If all terms in the energy are quadratic then the mean energy is spread equally over all degrees of freedom. (Hence, the name ``equipartition.'')


next up previous
Next: Harmonic Oscillators Up: Applications of Statistical Thermodynamics Previous: General Paramagnetism
Richard Fitzpatrick 2016-01-25