Equipartition theorem

From Wikipedia, the free encyclopedia

Figure 1. Thermal motion of an α-helical peptide.  The jittery motion is random and complex, and the energy of any particular atom can fluctuate wildly.  Nevertheless, the equipartition theorem allows the average kinetic energy of each atom to be computed, as well as the average potential energies of many vibrational modes.  The grey, red and blue spheres represent atoms of carbon, oxygen and nitrogen, respectively; the smaller white spheres represent atoms of hydrogen.
Figure 1. Thermal motion of an α-helical peptide. The jittery motion is random and complex, and the energy of any particular atom can fluctuate wildly. Nevertheless, the equipartition theorem allows the average kinetic energy of each atom to be computed, as well as the average potential energies of many vibrational modes. The grey, red and blue spheres represent atoms of carbon, oxygen and nitrogen, respectively; the smaller white spheres represent atoms of hydrogen.

In classical statistical mechanics, the equipartition theorem is a general formula that relates the temperature of a system with its average energies. The equipartition theorem is also known as the law of equipartition, equipartition of energy, or simply equipartition. The original idea of equipartition was that, in thermal equilibrium, energy is shared equally among its various forms; for example, the average kinetic energy in the translational motion of a molecule should equal the average kinetic energy in its rotational motion.

The equipartition theorem makes quantitative predictions. Like the virial theorem, it gives the total average kinetic and potential energies for a system at a given temperature, from which the system's heat capacity can be computed. However, equipartition also gives the average values of individual components of the energy, such as the kinetic energy of a particular particle or the potential energy of a single spring. For example, it predicts that every molecule in an ideal gas has an average kinetic energy of (3/2)kBT in thermal equilibrium, where kB is the Boltzmann constant and T is the temperature. More generally, it can be applied to any classical system in thermal equilibrium, no matter how complicated. The equipartition theorem can be used to derive the classical ideal gas law, and the Dulong–Petit law for the specific heat capacities of solids. It can also be used to predict the properties of stars, even white dwarfs and neutron stars, since it holds even when relativistic effects are considered.

Although the equipartition theorem makes very accurate predictions in certain conditions, it becomes inaccurate when quantum effects are significant, namely at low enough temperatures. When the thermal energy kBT is smaller than the quantum energy spacing in a particular degree of freedom, the average energy and heat capacity of this degree of freedom are less than the values predicted by equipartition. Such a degree of freedom is said to be "frozen out" when the thermal energy is much smaller than this spacing. For example, the specific heat of a solid decreases at low temperatures as various types of motion become frozen out, rather than remaining constant as predicted by equipartition. Such decreases in specific heat were the first sign to physicists of the 19th century that classical physics was incorrect and that new physics was needed. Along with other evidence, equipartition's failure for electromagnetic radiation — also known as the ultraviolet catastrophe — led Albert Einstein to suggest that light itself was quantized into photons, a revolutionary hypothesis that spurred the development of quantum mechanics and quantum field theory.

Contents

[edit] Basic concept and simple examples

See also: Kinetic energy and Heat capacity
Figure 2. Probability density functions of the molecular speed for four noble gases at a temperature of 298.15 K (25 °C).  The four gases are helium (4He), neon (20Ne), argon (40Ar) and xenon (132Xe); the superscripts indicate their mass numbers. These probability density functions have dimensions of probability times inverse speed; since probability is dimensionless, they can be expressed in units of seconds per meter.
Figure 2. Probability density functions of the molecular speed for four noble gases at a temperature of 298.15 K (25 °C). The four gases are helium (4He), neon (20Ne), argon (40Ar) and xenon (132Xe); the superscripts indicate their mass numbers. These probability density functions have dimensions of probability times inverse speed; since probability is dimensionless, they can be expressed in units of seconds per meter.

The name "equipartition" means "share and share alike". The original concept of equipartition was that the total kinetic energy of a system is shared equally among all of its independent parts, on the average, once the system has reached thermal equilibrium. Equipartition also makes quantitative predictions for these energies. For example, it predicts that every atom of a noble gas, in thermal equilibrium at temperature T, has an average translational kinetic energy of (3/2)kBT, where kB is the Boltzmann constant. As a consequence, the heavier atoms of xenon have a lower average speed than do the lighter atoms of helium at the same temperature. Figure 2 shows the Maxwell–Boltzmann distribution for the speeds of the atoms in four noble gases.

In this example, the key point is that the kinetic energy is quadratic in the velocity. The equipartition theorem shows that in thermal equilibrium, any degree of freedom (such as a component of the position or velocity of a particle) which appears only quadratically in the energy has an average energy of ½kBT and therefore contributes ½kB to the system's heat capacity. This has many applications.

[edit] Translational energy and ideal gases

See also: Ideal gas

The (Newtonian) kinetic energy of a particle of mass m, velocity v is given by


H^{\mathrm{kin}} = \tfrac12 m |\mathbf{v}|^2 = \tfrac{1}{2} m\left( v_{x}^{2} + v_{y}^{2} + v_{z}^{2} \right),

where vx, vy and vz are the cartesian components of the velocity v. Here, H is short for Hamiltonian, and used henceforth as a symbol for energy because the Hamiltonian formalism plays a central role in the most general form of the equipartition theorem.

Since the kinetic energy is quadratic in the components of the velocity, by equipartition these three components each contribute ½kBT to the average kinetic energy in thermal equilibrium. Thus the average kinetic energy of the particle is (3/2)kBT, as in the example of noble gases above.

More generally, in an ideal gas, the total energy consists purely of (translational) kinetic energy: by assumption, the particles have no internal degrees of freedom and move independently of one another. Equipartition therefore predicts that the average total energy of an ideal gas of N particles is (3/2) N kBT.

It follows that the heat capacity of the gas is (3/2) N kB and hence, in particular, the heat capacity of a mole of a such gas particles is (3/2)NAkB=(3/2)R, where NA is Avogadro's number and R is the gas constant. Since R ≈ 2 cal/(mol·K), equipartition predicts that the molar heat capacity of an ideal gas is roughly 3 cal/(mol·K). This prediction is confirmed by experiment.[1]

The mean kinetic energy also allows the root mean square speed vrms of the gas particles to be calculated:


v_{\mathrm{rms}} = \sqrt{\langle v^{2} \rangle} = \sqrt{\frac{3 k_{B} T}{m}} = \sqrt{\frac{3 R T}{M}},

where M = NAm is the mass of a mole of gas particles. This result is useful for many applications such as Graham's law of effusion, which provides a method for enriching uranium.[2]

[edit] Rotational energy and molecular tumbling in solution

See also: Angular velocity and Rotational diffusion

A similar example is provided by a rotating molecule with principal moments of inertia I1, I2 and I3. The rotational energy of such a molecule is given by


H^{\mathrm{rot}} = \tfrac{1}{2} ( I_{1} \omega_{1}^{2} + I_{2} \omega_{2}^{2} + I_{3} \omega_{3}^{2} ),

where ω1, ω2, and ω3 are the principal components of the angular velocity. By exactly the same reasoning as in the translational case, equipartition implies that in thermal equilibrium the average rotational energy of each particle is (3/2)kBT. Similarly, the equipartition theorem allows the average (more precisely, the root mean square) angular speed of the molecules to be calculated.[3]

The tumbling of rigid molecules — that is, the random rotations of molecules in solution — plays a key role in the relaxations observed by nuclear magnetic resonance, particularly protein NMR and residual dipolar couplings.[4] Rotational diffusion can also be observed by other biophysical probes such as fluorescence anisotropy, flow birefringence and dielectric spectroscopy.[5]

[edit] Potential energy and harmonic oscillators

Equipartition applies to potential energies as well as kinetic energies: important examples include harmonic oscillators such as a spring, which has a quadratic potential energy


H^{\mathrm{pot}} = \tfrac 12 a q^{2},\,

where the constant a describes the stiffness of the spring and q is the deviation from equilibrium. If such a one dimensional system has mass m, then its kinetic energy Hkin is ½mv2 = p2/2m, where v and p = mv denote the velocity and momentum of the oscillator. Combining these terms yields the total energy[6]


H = H^{\mathrm{kin}} + H^{\mathrm{pot}} = \frac{p^{2}}{2m} + \frac{1}{2} a q^{2}.

Equipartition therefore implies that in thermal equilibrium, the oscillator has average energy


\langle H \rangle = 
\langle H^{\mathrm{kin}} \rangle + \langle H^{\mathrm{pot}} \rangle = 
\tfrac{1}{2} k_{B} T + \tfrac{1}{2} k_{B} T = k_{B} T,

where the angular brackets \left\langle \ldots \right\rangle denote the average of the enclosed quantity,[7]

This result is valid for any type of harmonic oscillator, such as a pendulum, a vibrating molecule or a passive electronic oscillator. Systems of such oscillators arise in many situations; by equipartition, each such oscillator receives an average total energy kBT and hence contributes kB to the system's heat capacity. This can be used to derive the formula for Johnson–Nyquist noise[8] and the Dulong–Petit law of solid molar heat capacities. The latter application was particularly significant in the history of equipartition.

Figure 3. Atoms in a crystal can vibrate about their equilibrium positions in the lattice.  Such vibrations account largely for the heat capacity of crystalline dielectrics; with metals, electrons also contribute to the heat capacity.
Figure 3. Atoms in a crystal can vibrate about their equilibrium positions in the lattice. Such vibrations account largely for the heat capacity of crystalline dielectrics; with metals, electrons also contribute to the heat capacity.

[edit] Specific heat capacity of solids

For more details on the molar specific heat of solids, see Einstein solid and Debye model.

An important application of the equipartition theorem is to the specific heat capacity of a crystalline solid. Each atom in such a solid can oscillate in three independent directions, so the solid can be viewed as a system of 3N independent simple harmonic oscillators, where N denotes the number of atoms in the lattice. Since each harmonic oscillator has average energy kBT, the average total energy of the solid is 3NkBT, and its heat capacity is 3NkB.

By taking N to be Avogadro's number NA, and using the relation R = NAkB between the gas constant R and the Boltzmann constant kB, this provides an explanation for the Dulong–Petit law of molar heat capacities of solids, which states that the heat capacity per mole of atoms in the lattice is 3R ≈ 6 cal/(mol·K).

However, this law is inaccurate at lower temperatures, due to quantum effects; it is also inconsistent with the experimentally derived third law of thermodynamics, according to which the molar heat capacity of any substance must go to zero as the temperature goes to absolute zero.[8] A more accurate theory, incorporating quantum effects, was developed by Albert Einstein (1907) and Peter Debye (1911).[9]

Many other physical systems can be modeled as sets of coupled oscillators. The motions of such oscillators can be decomposed into normal modes, like the vibrational modes of a piano string or the resonances of an organ pipe. On the other hand, equipartition often breaks down for such systems, because there is no exchange of energy between the normal modes. In an extreme situation, the modes are independent and so their energies are independently conserved. This shows that some sort of mixing of energies, formally called ergodicity, is important for the law of equipartition to hold.

[edit] Sedimentation of particles

See also: Sedimentation, Mason–Weaver equation, and Brewing

Potential energies are not always quadratic in the position. However, the equipartition theorem also shows that if a degree of freedom x contributes only a multiple of xs (for a fixed real number s) to the energy, then in thermal equilibrium the average energy of that part is kBT/s.

There is a simple application of this extension to the sedimentation of particles under gravity.[10] For example, the haze sometimes seen in beer can be caused by clumps of proteins that scatter light.[11] Over time, these clumps settle downwards under the influence of gravity, causing more haze near the bottom of a bottle than near its top. However, in a process working in the opposite direction, the particles also diffuse back up towards the top of the bottle. Once equilibrium has been reached, the equipartition theorem may be used to determine the average position of a particular clump of buoyant mass mb. For an infinitely tall bottle of beer, the gravitational potential energy is given by


H^{\mathrm{grav}} = m_{b} g z\,,

where z is the height of the protein clump in the bottle and g is the acceleration caused by gravity. Since s=1, the average potential energy of a protein clump equals kBT. Hence, a protein clump with a buoyant mass of 10 MDa (roughly the size of a virus) would produce a haze with an average height of about 2 cm at equilibrium. The process of such sedimentation to equilibrium is described by the Mason–Weaver equation.[12]

[edit] History

This article uses the non-SI unit of cal/(mol·K) for molar specific heat, because it offers greater accuracy for single digits.
For an approximate conversion to the corresponding SI unit of
J/(mol·K), such values should be multiplied by 4.2 J/cal.

The equipartition of kinetic energy was proposed initially in 1843, and more correctly in 1845, by John James Waterston.[13] In 1859, James Clerk Maxwell argued that the kinetic heat energy of a gas is equally divided between linear and rotational energy.[14] In 1876, Ludwig Boltzmann expanded on this principle by showing that the average energy was divided equally among all the independent components of motion in a system.[15][16] Boltzmann applied the equipartition theorem to provide a theoretical explanation of the Dulong–Petit law for the specific heat capacities of solids.

Figure 4. Idealized plot of the molar specific heat of a diatomic gas against temperature. It agrees with the value (7/2)R predicted by equipartition at high temperatures (where R is the gas constant), but decreases to (5/2)R and then (3/2)R at lower temperatures, as the vibrational and rotational modes of motion are "frozen out". The failure of the equipartition theorem led to a paradox that was only resolved by quantum mechanics.  For most molecules, the transitional temperature Trot is much less than room temperature, whereas Tvib can be ten times larger or more.  A typical example is carbon monoxide, CO, for which Trot ≈ 2.8 K and Tvib ≈ 3103 K.  For molecules with very large or weakly bound atoms, Tvib can be close to room temperature (about 300 K); for example, Tvib ≈ 308 K for iodine gas, I2.
Figure 4. Idealized plot of the molar specific heat of a diatomic gas against temperature. It agrees with the value (7/2)R predicted by equipartition at high temperatures (where R is the gas constant), but decreases to (5/2)R and then (3/2)R at lower temperatures, as the vibrational and rotational modes of motion are "frozen out". The failure of the equipartition theorem led to a paradox that was only resolved by quantum mechanics. For most molecules, the transitional temperature Trot is much less than room temperature, whereas Tvib can be ten times larger or more. A typical example is carbon monoxide, CO, for which Trot ≈ 2.8 K and Tvib ≈ 3103 K. For molecules with very large or weakly bound atoms, Tvib can be close to room temperature (about 300 K); for example, Tvib ≈ 308 K for iodine gas, I2.[17]

The history of the equipartition theorem is intertwined with that of molar heat capacity, both of which were studied in the 19th century. In 1819, the French physicists Pierre Louis Dulong and Alexis Thérèse Petit discovered that the molar-specific heats of solids at room temperature were almost all identical, roughly 6 cal/(mol·K).[18] Their law was used for many years as a technique for measuring atomic weights.[9] However, subsequent studies by James Dewar and Heinrich Friedrich Weber showed that this Dulong–Petit law holds only at high temperatures;[19] at lower temperatures, or for exceptionally hard solids such as diamond, the specific heat was lower.[20]

Experimental observations of the specific heat of gases also raised concerns about the validity of the equipartition theorem. The theorem predicts that the molar heat capacity of simple monatomic gases should be roughly 3 cal/(mole·K), whereas that of diatomic gases should be roughly 7 cal/(mol·K). Experiments confirmed the former prediction,[1] but found that molar heat capacities of diatomic gases were typically about 5 cal/(mol·K),[21] and fell to about 3 cal/(mol·K) at very low temperatures.[22] Maxwell noted in 1875 that the disagreement between experiment and the equipartition theorem was much worse than even these numbers suggest;[23] since atoms have internal parts, heat energy should go into the motion of these internal parts, making the predicted specific heats of monatomic and diatomic gases much higher than 3 cal/(mol·K) and 7 cal/(mol·K), respectively.

A third discrepancy concerned the specific heat of metals.[24] According to the classical Drude model, metallic electrons act as a nearly ideal gas, and so they should contribute (3/2) Ne kB to the heat capacity by the equipartition theorem, where Ne is the number of electrons. Experimentally, however, electrons contribute little to the heat capacity: the molar heat capacities of many conductors and insulators are nearly the same.[24]

Several explanations of equipartition's failure to account for molar heat capacities were proposed. Boltzmann defended the derivation of his equipartition theorem as correct, but suggested that gases might not be in thermal equilibrium because of their interactions with the aether.[25] Lord Kelvin suggested that the derivation of the equipartition theorem must be incorrect, since it disagreed with experiment, but was unable to show how.[26] Lord Rayleigh instead put forward a more radical view that the equipartition theorem and the experimental assumption of thermal equilibrium were both correct; to reconcile them, he noted the need for a new principle that would provide an "escape from the destructive simplicity" of the equipartition theorem.[27] Albert Einstein provided that escape, by showing in 1907 that these anomalies in the specific heat were due to quantum effects, specifically the quantization of energy in the elastic modes of the solid.[28] Einstein used the failure of equipartition to argue for the need of a new quantum theory of matter.[9] Nernst's 1910 measurements of specific heats at low temperatures[29] supported Einstein's theory, and led to the widespread acceptance of quantum theory among physicists.[30]

[edit] General formulation of the equipartition theorem

See also: Generalized coordinates, Hamiltonian mechanics, Microcanonical ensemble, and Canonical ensemble

The most general form of the equipartition theorem[10][7][3] states that under suitable assumptions (discussed below), for a physical system with Hamiltonian energy function H and degrees of freedom xn, the following equipartition formula holds in thermal equilibrium for all indices m and n:

\!
\Bigl\langle x_{m} \frac{\partial H}{\partial x_{n}} \Bigr\rangle = \delta_{mn} k_{B} T.

Here δmn is the Kronecker delta, which is equal to one if m=n and is zero otherwise. The averaging brackets \left\langle \ldots \right\rangle may refer either to the long time average of a single system, or, more commonly, the ensemble average over phase space. The ergodicity assumptions implicit in the theorem imply that these two averages agree, and both have been used to estimate internal energies of complex physical systems.

The general equipartition theorem holds in both the microcanonical ensemble,[7] when the total energy of the system is constant, and also in the canonical ensemble,[31][3] when the system is coupled to a heat bath with which it can exchange energy. Derivations of the general formula are given later in the article.

The general formula is equivalent to the following two.

  1. \Bigl\langle x_{n} \frac{\partial H}{\partial x_{n}} \Bigr\rangle = k_{B} T    for all n.
  2. \Bigl\langle x_{m} \frac{\partial H}{\partial x_{n}} \Bigr\rangle = 0    for all mn.

If a degree of freedom xn appears only as a quadratic term anxn2 in the Hamiltonian H, then the first of these formulae implies that


k_{B} T = \Bigl\langle x_{n} \frac{\partial H}{\partial x_{n}}\Bigr\rangle = 2\langle a_n x_n^2 \rangle,

which is twice the contribution that this degree of freedom makes to the average energy \langle H\rangle. Thus the equipartition theorem for systems with quadratic energies follows easily from the general formula. A similar argument, with 2 replaced by s, applies to energies of the form anxns.

The degrees of freedom xn are coordinates on the phase space of the system and are therefore commonly subdivided into generalized position coordinates qk and generalized momentum coordinates pk, where pk is the conjugate momentum to qk. In this situation, formula 1 means that for all k,


\Bigl\langle p_{k} \frac{\partial H}{\partial p_{k}} \Bigr\rangle =  \Bigl\langle q_{k} \frac{\partial H}{\partial q_{k}} \Bigr\rangle = k_{B} T.

Using the equations of Hamiltonian mechanics,[6] these formulae may also be written


\Bigl\langle p_{k} \frac{dq_{k}}{dt} \Bigr\rangle = -\Bigl\langle q_{k} \frac{dp_{k}}{dt} \Bigr\rangle = k_{B} T.

Formula 2 additionally states that the averages


\Bigl\langle q_{j} \frac{\partial H}{\partial q_{k}} \Bigr\rangle, \quad
\Bigl\langle q_{j} \frac{\partial H}{\partial p_{k}} \Bigr\rangle, \quad
\Bigl\langle p_{j} \frac{\partial H}{\partial p_{k}} \Bigr\rangle, \quad
\Bigl\langle p_{j} \frac{\partial H}{\partial q_{k}} \Bigr\rangle, \quad
\Bigl\langle q_{k} \frac{\partial H}{\partial p_{k}} \Bigr\rangle,   and   
\Bigl\langle p_{k} \frac{\partial H}{\partial q_{k}} \Bigr\rangle

are all zero for j≠k.

[edit] Relation to the virial theorem

See also: Virial theorem, Generalized coordinates, and Hamiltonian mechanics

The general equipartition theorem is an extension of the virial theorem (proposed in 1870[32]), which states that


\Bigl\langle \sum_{k} q_{k} \frac{\partial H}{\partial q_{k}} \Bigr\rangle = 
\Bigl\langle \sum_{k} p_{k} \frac{\partial H}{\partial p_{k}} \Bigr\rangle = 
\Bigl\langle \sum_{k} p_{k} \frac{dq_{k}}{dt} \Bigr\rangle = -\Bigl\langle \sum_{k} q_{k} \frac{dp_{k}}{dt} \Bigr\rangle,

where t denotes time.[6] Two key differences are that the virial theorem relates summed rather than individual averages to each other, and it does not connect them to the temperature T. Another difference is that traditional derivations of the virial theorem use averages over time, whereas those of the equipartition theorem use averages over phase space.

[edit] Applications

[edit] Ideal gas law

See also: Ideal gas and Ideal gas law

Ideal gases provide an important application of the equipartition theorem. As well as providing the formula


\begin{align}
\langle H^{\mathrm{kin}} \rangle &= \frac{1}{2m} \langle p_{x}^{2} + p_{y}^{2} + p_{z}^{2} \rangle\\
&=  \frac{1}{2} \biggl(
\Bigl\langle p_{x} \frac{\partial H^{\mathrm{kin}}}{\partial p_{x}} \Bigr\rangle + 
\Bigl\langle p_{y} \frac{\partial H^{\mathrm{kin}}}{\partial p_{y}} \Bigr\rangle +
\Bigl\langle p_{z} \frac{\partial H^{\mathrm{kin}}}{\partial p_{z}} \Bigr\rangle \biggr) = 
\frac{3}{2} k_{B} T
\end{align}

for the average kinetic energy per particle, the equipartition theorem can be used to derive the ideal gas law from classical mechanics.[3] If q = (qx, qy, qz) and p = (px, py, pz) denote the position vector and momentum of a particle in the gas, and F is the net force on that particle, then


\begin{align}
\langle \mathbf{q} \cdot \mathbf{F} \rangle &= \Bigl\langle q_{x} \frac{dp_{x}}{dt} \Bigr\rangle + 
\Bigl\langle q_{y} \frac{dp_{y}}{dt} \Bigr\rangle + 
\Bigl\langle q_{z} \frac{dp_{z}}{dt} \Bigr\rangle\\
&=-\Bigl\langle q_{x} \frac{\partial H}{\partial q_x} \Bigr\rangle -
\Bigl\langle q_{y} \frac{\partial H}{\partial q_y} \Bigr\rangle - 
\Bigl\langle q_{z} \frac{\partial H}{\partial q_z} \Bigr\rangle = -3k_{B} T,
\end{align}

where the first equality is Newton's second law, and the second line uses Hamilton's equations and the equipartition formula. Summing over a system of N particles yields


3Nk_{B} T = - \biggl\langle \sum_{k=1}^{N} \mathbf{q}_{k} \cdot \mathbf{F}_{k} \biggr\rangle.
Figure 5. The kinetic energy of a particular molecule can fluctuate wildly, but the equipartition theorem allows its average energy to be calculated at any temperature.  Equipartition also provides a derivation of the ideal gas law, an equation that relates the pressure, volume and temperature of the gas. (In this diagram five of the molecules have been colored red to track their motion; this coloration has no other significance.)
Figure 5. The kinetic energy of a particular molecule can fluctuate wildly, but the equipartition theorem allows its average energy to be calculated at any temperature. Equipartition also provides a derivation of the ideal gas law, an equation that relates the pressure, volume and temperature of the gas. (In this diagram five of the molecules have been colored red to track their motion; this coloration has no other significance.)

By Newton's third law and the ideal gas assumption, the net force on the system is the force applied by the walls of their container, and this force is given by the pressure P of the gas. Hence


-\biggl\langle\sum_{k=1}^{N} \mathbf{q}_{k} \cdot \mathbf{F}_{k}\biggr\rangle = P \oint_{\mathrm{surface}} \mathbf{q} \cdot \mathbf{dS},

where dS is the infinitesimal area element along the walls of the container. Since the divergence of the position vector q is


\boldsymbol\nabla \cdot \mathbf{q} =
\frac{\partial q_{x}}{\partial q_{x}} + 
\frac{\partial q_{y}}{\partial q_{y}} + 
\frac{\partial q_{z}}{\partial q_{z}} = 3,

the divergence theorem implies that


P \oint_{\mathrm{surface}} \mathbf{q} \cdot \mathbf{dS} = P \int_{\mathrm{volume}} \left( \boldsymbol\nabla \cdot \mathbf{q} \right) dV = 3PV,

where dV is an infinitesimal volume within the container and V is the total volume of the container.

Putting these equalities together yields


3Nk_{B} T = -\biggl\langle \sum_{k=1}^{N} \mathbf{q}_{k} \cdot \mathbf{F}_{k} \biggr\rangle = 3PV,

which immediately implies the ideal gas law for N particles:


PV = Nk_{B} T = nRT,\,

where n=N/NA is the number of moles of gas and R=NAkB is the gas constant. Although equipartition provides a simple derivation of the ideal-gas law and the internal energy, the same results can be obtained by an alternative method using the partition function[33].

[edit] Diatomic gases

See also: Two-body problem, Rigid rotor, and Harmonic oscillator

A diatomic gas can be modelled as two masses, m1 and m2, joined by a spring of stiffness a, which is called the rigid rotor-harmonic oscillator approximation.[17] The classical energy of this system is


H = 
\frac{\left| \mathbf{p}_{1} \right|^{2}}{2m_{1}} + 
\frac{\left| \mathbf{p}_{2} \right|^{2}}{2m_{2}} + 
\frac{1}{2} a q^{2},

where p1 and p2 are the momenta of the two atoms, and q is the deviation of the inter-atomic separation from its equilibrium value. Every degree of freedom in the energy is quadratic and, thus, should contribute ½kBT to the total average energy, and ½kB to the heat capacity. Therefore, the heat capacity of a gas of N diatomic molecules is predicted to be 7N · ½kB: the momenta p1 and p2 contribute three degrees of freedom each, and the extension q contributes the seventh. It follows that the heat capacity of a mole of diatomic molecules with no other degrees of freedom should be (7/2)NAkB=(7/2)R and, thus, the predicted molar heat capacity should be roughly 7 cal/(mol·K). However, the experimental values for molar heat capacities of diatomic gases are typically about 5 cal/(mol·K)[21] and fall to 3 cal/(mol·K) at very low temperatures.[22] This disagreement between the equipartition prediction and the experimental value of the molar heat capacity cannot be explained by using a more complex model of the molecule, since adding more degrees of freedom can only increase the predicted specific heat, not decrease it.[23] This discrepancy was a key piece of evidence showing the need for a quantum theory of matter.

Figure 6. A combined X-ray and optical image of the Crab Nebula. At the heart of this nebula there is a rapidly rotating neutron star which has about one and a half times the mass of the Sun but is only 25 km across (roughly the size of Madrid).  The equipartition theorem is useful in predicting the properties of such neutron stars.
Figure 6. A combined X-ray and optical image of the Crab Nebula. At the heart of this nebula there is a rapidly rotating neutron star which has about one and a half times the mass of the Sun but is only 25 km across (roughly the size of Madrid). The equipartition theorem is useful in predicting the properties of such neutron stars.

[edit] Extreme relativistic ideal gases

See also: Special relativity, White dwarf, and Neutron star

Equipartition was used above to derive the classical ideal gas law from Newtonian mechanics. However, relativistic effects become dominant in some systems, such as white dwarfs and neutron stars,[7] and the ideal gas equations must be modified. The equipartition theorem provides a convenient way to derive the corresponding laws for an extreme relativistic ideal gas.[3] In such cases, the kinetic energy of a single particle is given by the formula


H^{\mathrm{kin}} \approx cp = c \sqrt{p_{x}^{2} + p_{y}^{2} + p_{z}^{2}}.

Taking the derivative of H with respect to the px momentum component gives the formula


p_{x} \frac{\partial H^{\mathrm{kin}}}{\partial p_{x}}  = c \frac{p_{x}^{2}}{\sqrt{p_{x}^{2} + p_{y}^{2} + p_{z}^{2}}}

and similarly for the py and pz components. Adding the three components together gives


\begin{align}
\langle H^{\mathrm{kin}} \rangle
&= \biggl\langle c \frac{p_{x}^{2} + p_{y}^{2} + p_{z}^{2}}{\sqrt{p_{x}^{2} + p_{y}^{2} + p_{z}^{2}}}  \biggr\rangle\\
&= \Bigl\langle p_{x} \frac{\partial H^{\mathrm{kin}}}{\partial p_{x}} \Bigr\rangle + 
\Bigl\langle p_{y} \frac{\partial H^{\mathrm{kin}}}{\partial p_{y}} \Bigr\rangle + 
\Bigl\langle p_{z} \frac{\partial H^{\mathrm{kin}}}{\partial p_{z}} \Bigr\rangle\\
&= 3 k_{B} T
\end{align}

where the last equality follows from the equipartition formula. Thus, the average total energy of an extreme relativistic gas is twice that of the non-relativistic case: for N particles, it is 3 N kBT.

[edit] Non-ideal gases

See also: Virial expansion and Virial coefficient

In an ideal gas the particles are assumed to interact only through collisions. The equipartition theorem may also be used to derive the energy and pressure of "non-ideal gases" in which the particles also interact with one another through conservative forces whose potential U(r) depends only on the distance r between the particles.[3] This situation can be described by first restricting attention to a single gas particle, and approximating the rest of the gas by a spherically symmetric distribution. It is then customary to introduce a radial distribution function g(r) such that the probability density of finding another particle at a distance r from the given particle is equal to 4πr2ρ g(r), where ρ=N/V is the mean density of the gas.[34] It follows that the mean potential energy associated to the interaction of the given particle with the rest of the gas is


\langle h^{\mathrm{pot}} \rangle = \int_{0}^{\infty} 4\pi r^{2} \rho U(r) g(r)\, dr.

The total mean potential energy of the gas is therefore  \langle H^{pot} \rangle = \tfrac12 N \langle h^{\mathrm{pot}} \rangle , where N is the number of particles in the gas, and the factor ½ is needed because summation over all the particles counts each interaction twice. Adding kinetic and potential energies, then applying equipartition, yields the energy equation


H = 
\langle H^{\mathrm{kin}} \rangle + \langle H^{\mathrm{pot}} \rangle = 
\frac{3}{2} Nk_{B}T + 2\pi N \rho \int_{0}^{\infty} r^{2} U(r) g(r) \, dr.

A similar argument,[3] can be used to derive the pressure equation


3Nk_{B}T = 3PV + 2\pi N \rho \int_{0}^{\infty} r^{3} U'(r) g(r)\, dr.

[edit] Anharmonic oscillators

See also: Anharmonic oscillator

An anharmonic oscillator (in contrast to a simple harmonic oscillator) is one in which the potential energy is not quadratic in the extension q (the generalized position which measures the deviation of the system from equilibrium). Such oscillators provide a complementary point of view on the equipartition theorem.[35][36] Simple examples are provided by potential energy functions of the form


H^{\mathrm{pot}} = C q^{s},\,

where C and s are arbitrary real constants. In these cases, the law of equipartition predicts that


k_{B} T = \Bigl\langle q \frac{\partial H^{\mathrm{pot}}}{\partial q} \Bigr\rangle = 
\langle q \cdot s C q^{s-1} \rangle = \langle s C q^{s} \rangle = s \langle H^{\mathrm{pot}} \rangle.

Thus, the average potential energy equals kBT/s, not kBT/2 as for the quadratic harmonic oscillator (where s=2).

More generally, a typical energy function of a one-dimensional system has a Taylor expansion in the extension q:


H^{\mathrm{pot}} = \sum_{n=2}^{\infty} C_{n} q^{n}

for non-negative integers n. There is no n=1 term, because at the equilibrium point, there is no net force and so the first derivative of the energy is zero. The n=0 term need not be included, since the energy at the equilibrium position may be set to zero by convention. In this case, the law of equipartition predicts that[35]


k_{B} T = \Bigl\langle q \frac{\partial H^{\mathrm{pot}}}{\partial q} \Bigr\rangle = 
\sum_{n=2}^{\infty} \langle q \cdot n C_{n} q^{n-1} \rangle = 
\sum_{n=2}^{\infty} n C_{n} \langle q^{n} \rangle.

In contrast to the other examples cited here, the equipartition formula


\langle H^{\mathrm{pot}} \rangle = \frac{1}{2} k_{B} T - 
\sum_{n=3}^{\infty} \left( \frac{n - 2}{2} \right) C_{n} \langle q^{n} \rangle

does not allow the average potential energy to be written in terms of known constants.

[edit] Brownian motion

Figure 7. Typical Brownian motion of a particle in three dimensions.
Figure 7. Typical Brownian motion of a particle in three dimensions.
See also: Brownian motion

The equipartition theorem can be used to derive the Brownian motion of a particle from the Langevin equation.[3] According to that equation, the motion of a particle of mass m with velocity v is governed by Newton's second law


\frac{d\mathbf{v}}{dt} = \frac{1}{m} \mathbf{F} = -\frac{\mathbf{v}}{\tau} + \frac{1}{m} \mathbf{F}^{\mathrm{rnd}},

where Frnd is a random force representing the random collisions of the particle and the surrounding molecules, and where the time constant τ reflects the drag force that opposes the particle's motion through the solution. The drag force is often written Fdrag = - γv; therefore, the time constant τ equals m/γ.

The dot product of this equation with the position vector r, after averaging, yields the equation


\Bigl\langle \mathbf{r} \cdot \frac{d\mathbf{v}}{dt} \Bigr\rangle + 
\frac{1}{\tau} \langle \mathbf{r} \cdot \mathbf{v} \rangle = 0

for Brownian motion (since the random force Frnd is uncorrelated with the position r). Using the mathematical identities


\frac{d}{dt} \left( \mathbf{r} \cdot \mathbf{r} \right) = 
\frac{d}{dt} \left( r^{2} \right) = 2 \left( \mathbf{r} \cdot \mathbf{v} \right)

and


\frac{d}{dt} \left( \mathbf{r} \cdot \mathbf{v} \right) = v^{2} + \mathbf{r} \cdot \frac{d\mathbf{v}}{dt},

the basic equation for Brownian motion can be transformed into


\frac{d^{2}}{dt^{2}} \langle r^{2} \rangle + \frac{1}{\tau} \frac{d}{dt} \langle r^{2} \rangle = 
2 \langle v^{2} \rangle = \frac{6}{m} k_{B} T,

where the last equality follows from the equipartition theorem for translational kinetic energy:


\langle H^{\mathrm{kin}} \rangle = \Bigl\langle \frac{p^{2}}{2m} \Bigr\rangle = \langle \tfrac{1}{2} m v^{2} \rangle = \tfrac{3}{2} k_{B} T.

The above differential equation for \langle r^2\rangle (with suitable initial conditions) may be solved exactly:


\langle r^{2} \rangle = \frac{6k_{B} T \tau^{2}}{m} \left( e^{-t/\tau} - 1 + \frac{t}{\tau} \right).

On small time scales, with t << τ, the particle acts as a freely moving particle: by the Taylor series of the exponential function, the squared distance grows approximately quadratically:


\langle r^{2} \rangle \approx \frac{3k_{B} T}{m} t^{2} = \langle v^{2} \rangle t^{2}.

However, on long time scales, with t >> τ, the exponential and constant terms are negligible, and the squared distance grows only linearly:


\langle r^{2} \rangle \approx \frac{6k_{B} T\tau}{m} t = 6\gamma k_{B} T t.

This describes the diffusion of the particle over time. An analogous equation for the rotational diffusion of a rigid molecule can be derived in a similar way.

Figure 8. Equipartition gives an accurate estimate of the temperature at the Sun's core.
Figure 8. Equipartition gives an accurate estimate of the temperature at the Sun's core.

[edit] Stellar physics

See also: Astrophysics and Stellar structure

The equipartition theorem and the related virial theorem have long been used as a tool in astrophysics.[37] As examples, the virial theorem may be used to estimate stellar temperatures or the Chandrasekhar limit on the mass of white dwarf stars.[38][39]

The average temperature of a star can be estimated from the equipartition theorem.[40] Since most stars are spherically symmetric, the total gravitational potential energy can be estimated by integration


H^{\mathrm{grav}}_{\mathrm{tot}} = -\int_{0}^{R} \frac{4\pi r^{2} G}{r} M(r)\, \rho(r)\, dr,

where M(r) is the mass within a radius r and ρ(r) is the stellar density at radius r; G represents the gravitational constant and R the total radius of the star. Assuming a constant density throughout the star, this integration yields the formula


H^{\mathrm{grav}}_{\mathrm{tot}} = - \frac{3G M^{2}}{5R},

where M is the star's total mass. Hence, the average potential energy of a single particle is


\langle H^{\mathrm{grav}} \rangle = \frac{H^{\mathrm{grav}}_{\mathrm{tot}}}{N} = - \frac{3G M^{2}}{5RN},

where N is the number of particles in the star. Since most stars are composed mainly of ionized hydrogen, N equals roughly (M/mp), where mp is the mass of one proton. Application of the equipartition theorem gives an estimate of the star's temperature


\Bigl\langle r \frac{\partial H^{\mathrm{grav}}}{\partial r} \Bigr\rangle = \langle -H^{\mathrm{grav}} \rangle = 
k_{B} T = \frac{3G M^{2}}{5RN}.

Substitution of the mass and radius of the Sun yields an estimated solar temperature of T = 14 million kelvins, very close to its core temperature of 15 million kelvins. However, the Sun is much more complex than assumed by this model — both its temperature and density vary strongly with radius — and such excellent agreement (≈7% relative error) is partly fortuitous.[41]

[edit] Star formation

The same formulae may be applied to determining the conditions for star formation in giant molecular clouds.[42] A local fluctuation in the density of such a cloud can lead to a runaway condition in which the cloud collapses inwards under its own gravity. Such a collapse occurs when the equipartition theorem — or, equivalently, the virial theorem — is no longer valid, i.e., when the gravitational potential energy exceeds twice the kinetic energy


\frac{3G M^{2}}{5R} > 3 N k_{B} T

Assuming a constant density ρ for the cloud


M = \frac{4}{3} \pi R^{3} \rho

yields a minimum mass for stellar contraction, the Jeans mass MJ


M_{J}^{2} = \left( \frac{5k_{B}T}{G m_{p}} \right)^{3} \left( \frac{3}{4\pi \rho} \right)

Substituting the values typically observed in such clouds (T=150 K, ρ = 2×10-16 g/cm3) gives an estimated minimum mass of 17 solar masses, which is consistent with observed star formation. This effect is also known as the Jeans instability, after the British physicist James Hopwood Jeans who published it in 1902.[43]

[edit] Derivations

[edit] Kinetic energies and the Maxwell–Boltzmann distribution

The original formulation of the equipartition theorem states that, in any physical system in thermal equilibrium, every particle has exactly the same average kinetic energy, (3/2)kBT.[44] This may be shown using the Maxwell–Boltzmann distribution (see Figure 2), which is the probability distribution


f (v) = 4 \pi 
\left( \frac{m}{2 \pi k_B T}\right)^{3/2}\!\!v^2
\exp \Bigl(
\frac{-mv^2}{2k_B T}
\Bigr)

for the speed of a particle of mass m in the system, where the speed v is the magnitude \sqrt{v_x^2 + v_y^2 + v_z^2} of the velocity vector \mathbf{v} = (v_x,v_y,v_z).

The Maxwell–Boltzmann distribution applies to any system composed of atoms, and assumes only a canonical ensemble, specifically, that the kinetic energies are distributed according to their Boltzmann factor at a temperature T.[44] The average kinetic energy for a particle of mass m is then given by the integral formula


\langle H^{\mathrm{kin}} \rangle = 
\langle \tfrac{1}{2} m v^{2} \rangle = 
\int _{0}^{\infty} \tfrac{1}{2} m v^{2}\  f(v)\  dv = \tfrac{3}{2} k_{B} T,

as stated by the equipartition theorem. The same result can also be obtained by averaging the particle energy using the probability of finding the particle in certain quantum energy state[33].

[edit] Quadratic energies and the partition function

More generally, the equipartition theorem states that any degree of freedom x which appears in the total energy H only as a simple quadratic term Ax2, where A is a constant, has an average energy of ½kBT in thermal equilibrium. In this case the equipartition theorem may be derived from the partition function Z(β), where β=1/(kBT) is the canonical inverse temperature.[45] Integration over the variable x yields a factor


Z_{x} = \int_{-\infty}^{\infty} dx \ e^{-\beta A x^{2}} = \sqrt{\frac{\pi}{\beta A}},

in the formula for Z. The mean energy associated with this factor is given by


\langle H_{x} \rangle = - \frac{\partial \log Z_{x}}{\partial \beta} = \frac{1}{2\beta} = \frac{1}{2} k_{B} T

as stated by the equipartition theorem.

[edit] General proofs

General derivations of the equipartition theorem can be found in many statistical mechanics textbooks, both for the microcanonical ensemble[7][3] and for the canonical ensemble.[31][3] They involve taking averages over the phase space of the system, which is a symplectic manifold.

To explain these derivations, the following notation is introduced. First, the phase space is described in terms of generalized position coordinates qj together with their conjugate momenta pj. The quantities qj completely describe the configuration of the system, while the quantities (qj,pj) together completely describe its state.

Secondly, the infinitesimal volume

dΓ = dqidpi
i

of the phase space is introduced and used to define the volume Γ(E, ΔE) of the portion of phase space where the energy H of the system lies between two limits, E and E+ΔE:


\Gamma (E, \Delta E) = \int_{H \in \left[E, E+\Delta E \right]} d\Gamma .

In this expression, ΔE is assumed to be very small, ΔE<<E. Similarly, Σ(E) is defined to be the total volume of phase space where the energy is less than E:

Σ(E) = dΓ.
H < E

Since ΔE is very small, the following integrations are equivalent


\int_{H \in \left[ E, E+\Delta E \right]} \ldots d\Gamma  = \Delta E \frac{\partial}{\partial E} \int_{H < E} \ldots d\Gamma,

where the ellipses represent the integrand. From this, it follows that Γ is proportional to ΔE


\Gamma = \Delta E \ \frac{\partial \Sigma}{\partial E} = \Delta E \ \rho(E),

where ρ(E) is the density of states. By the usual definitions of statistical mechanics, the entropy S equals kB log Σ(E), and the temperature T is defined by


\frac{1}{T} = \frac{\partial S}{\partial E} = k_{b} \frac{\partial \log \Sigma}{\partial E} = k_{b} \frac{1}{\Sigma}\,\frac{\partial \Sigma}{\partial E} .

[edit] The canonical ensemble

In the canonical ensemble, the system is in thermal equilibrium with an infinite heat bath at temperature T (in Kelvin).[31][3] The probability of each state in phase space is given by its Boltzmann factor times a normalization factor \mathcal{N}, which is chosen so that the probabilities sum to one


\mathcal{N} \int e^{-\beta H(p, q)} d\Gamma = 1,

where β = 1/kBT. Integration by parts for a phase-space variable xk (which could be either qk or pk) between two limits a and b yields the equation


\mathcal{N} \int  \left[ e^{-\beta H(p, q)} x_{k} \right]_{x_{k}=a}^{x_{k}=b} d\Gamma_{k}+ 
\mathcal{N} \int  e^{-\beta H(p, q)} x_{k} \beta \frac{\partial H}{\partial x_{k}} d\Gamma = 1,

where k = dΓ/dxk, i.e., the first integration is not carried out over xk. The first term is usually zero, either because xk is zero at the limits, or because the energy goes to infinity at those limits. In that case, the equipartition theorem for the canonical ensemble follows immediately


\mathcal{N} \int e^{-\beta H(p, q)} x_{k} \frac{\partial H}{\partial x_{k}} \,d\Gamma = 
\Bigl\langle x_{k} \frac{\partial H}{\partial x_{k}} \Bigr\rangle = \frac{1}{\beta} = k_{B} T.

Here, the averaging symbolized by \langle \ldots \rangle is the ensemble average taken over the canonical ensemble.

[edit] The microcanonical ensemble

In the microcanonical ensemble, the system is isolated from the rest of the world, or at least very weakly coupled to it.[7] Hence, its total energy is effectively constant; to be definite, we say that the total energy H is confined between E and E+ΔE. For a given energy E and spread ΔE, there is a region of phase space Γ in which the system has that energy, and the probability of each state in that region of phase space is equal, by the definition of the microcanonical ensemble. Given these definitions, the equipartition average of phase-space variables xm (which could be either qkor pk) and xn is given by


\begin{align}
\Bigl\langle x_{m} \frac{\partial H}{\partial x_{n}} \Bigr \rangle &=
\frac{1}{\Gamma}   \, \int_{H \in \left[ E, E+\Delta E \right]}  x_{m} \frac{\partial H}{\partial x_{n}} \,d\Gamma\\
&=\frac{\Delta E}{\Gamma}\, \frac{\partial}{\partial E} \int_{H < E}  x_{m} \frac{\partial H}{\partial x_{n}} \,d\Gamma\\
&= \frac{1}{\rho} \,\frac{\partial}{\partial E} \int_{H < E}  x_{m} \frac{\partial \left( H - E \right)}{\partial x_{n}} \,d\Gamma,
\end{align}

where the last equality follows because E is a constant that does not depend on xn. Integrating by parts yields the relation


\begin{align}
\int_{H < E}  x_{m} \frac{\partial ( H - E )}{\partial x_{n}} \,d\Gamma &= 
\int_{H < E}  \frac{\partial}{\partial x_{n}} \bigl( x_{m} ( H - E ) \bigr) \,d\Gamma - 
\int_{H < E}  \delta_{mn} ( H - E ) d\Gamma\\
&=  \delta_{mn} \int_{H < E} ( E - H ) \,d\Gamma,
\end{align}

since the first term on the right hand side of the first line is zero (it can be rewritten as an integral of H - E on the hypersurface where H = E).

Substitution of this result into the previous equation yields


\Bigl\langle x_{m} \frac{\partial H}{\partial x_{n}} \Bigr\rangle = 
\delta_{mn} \frac{1}{\rho} \, \frac{\partial}{\partial E} \int_{H < E}\left( E - H \right)\,d\Gamma  = 
\delta_{mn}  \frac{1}{\rho}  \, \int_{H < E} \,d\Gamma = 
\delta_{mn}  \frac{\Sigma}{\rho}.

Since  \rho = \frac{\partial \Sigma}{\partial E} the equipartition theorem follows:


\Bigl\langle x_{m} \frac{\partial H}{\partial x_{n}} \Bigr\rangle = 
\delta_{mn} \Bigl(\frac{1}{\Sigma} \frac{\partial \Sigma}{\partial E}\Bigr)^{-1}  = 
\delta_{mn} \Bigl(\frac{\partial \log \Sigma} {\partial E}\Bigr)^{-1} = \delta_{mn} k_{B} T.

Thus, we have derived the general formulation of the equipartition theorem

\!
\Bigl\langle x_{m} \frac{\partial H}{\partial x_{n}} \Bigr\rangle = \delta_{mn} k_{B} T,

which was so useful in the applications described above.

[edit] Limitations

Figure 9. Energy is not shared among the various normal modes in an isolated system of ideal coupled oscillators; the energy in each mode is constant and independent of the energy in the other modes.  Hence, the equipartition theorem does not hold for such a system in the microcanonical ensemble (when isolated), although it does hold in the canonical ensemble (when coupled to a heat bath).  However, by adding a sufficiently strong nonlinear coupling between the modes, energy will be shared and equipartition holds in both ensembles.
Figure 9. Energy is not shared among the various normal modes in an isolated system of ideal coupled oscillators; the energy in each mode is constant and independent of the energy in the other modes. Hence, the equipartition theorem does not hold for such a system in the microcanonical ensemble (when isolated), although it does hold in the canonical ensemble (when coupled to a heat bath). However, by adding a sufficiently strong nonlinear coupling between the modes, energy will be shared and equipartition holds in both ensembles.

[edit] Requirement of ergodicity

See also: Ergodicity, Chaos theory, and Kolmogorov–Arnold–Moser theorem

The law of equipartition holds only for ergodic systems in thermal equilibrium, which implies that all states with the same energy must be equally likely to be populated.[7] Consequently, it must be possible to exchange energy among all its various forms within the system, or with an external heat bath in the canonical ensemble. The number of physical systems that have been rigorously proven to be ergodic is small; a famous example is the hard-sphere system of Yakov Sinai.[46] The requirements for isolated systems to ensure ergodicity — and, thus equipartition — have been studied, and provided motivation for the modern chaos theory of dynamical systems. A chaotic Hamiltonian system need not be ergodic, although that is usually a good assumption.[47]

A commonly cited counter-example where energy is not shared among its various forms and where equipartition does not hold in the microcanonical ensemble is a system of coupled harmonic oscillators.[47] If the system is isolated from the rest of the world, the energy in each normal mode is constant; energy is not transferred from one mode to another. Hence, equipartition does not hold for such a system; the amount of energy in each normal mode is fixed at its initial value. If sufficiently strong nonlinear terms are present in the energy function, energy may be transferred between the normal modes, leading to ergodicity and rendering the law of equipartition valid. However, the Kolmogorov–Arnold–Moser theorem states that energy will not be exchanged unless the nonlinear perturbations are strong enough; if they are too small, the energy will remain trapped in at least some of the modes.

[edit] Failure due to quantum effects

See also: Ultraviolet catastrophe, History of quantum mechanics, and Identical particles

The law of equipartition breaks down when the thermal energy kBT is significantly smaller than the spacing between energy levels. Equipartition no longer holds because it is a poor approximation to assume that the energy levels form a smooth continuum, which is required in the derivations of the equipartition theorem above.[7][3] Historically, the failures of the classical equipartition theorem to explain specific heats and blackbody radiation were critical in showing the need for a new theory of matter and radiation, namely, quantum mechanics and quantum field theory.[9]

Figure 10.  Log-log plot of the average energy of a quantum mechanical oscillator (shown in red) as a function of temperature.  For comparison, the value predicted by the equipartition theorem is shown in black.  At high temperatures, the two agree nearly perfectly, but at low temperatures when kBT << hν, the quantum mechanical value decreases much more rapidly.  This resolves the problem of the ultraviolet catastrophe: for a given temperature, the energy in the high-frequency modes (where hν >> kBT)  is almost zero.
Figure 10. Log-log plot of the average energy of a quantum mechanical oscillator (shown in red) as a function of temperature. For comparison, the value predicted by the equipartition theorem is shown in black. At high temperatures, the two agree nearly perfectly, but at low temperatures when kBT << hν, the quantum mechanical value decreases much more rapidly. This resolves the problem of the ultraviolet catastrophe: for a given temperature, the energy in the high-frequency modes (where hν >> kBT) is almost zero.

To illustrate the breakdown of equipartition, consider the average energy in a single (quantum) harmonic oscillator, which was discussed above for the classical case. Its quantum energy levels are given by En = nhν, where h is Planck's constant, ν is the fundamental frequency of the oscillator, and n is an integer. The probability of a given energy level being populated in the canonical ensemble is given by its Boltzmann factor


P(E_{n}) = \frac{e^{-n\beta h\nu}}{Z},

where β = 1/kBT and the denominator Z is the partition function, here a geometric series


Z = \sum_{n=0}^{\infty} e^{-n\beta h\nu} = \frac{1}{1 - e^{-\beta h\nu}}.

Its average energy is given by


\langle H \rangle = \sum_{n=0}^{\infty} E_{n} P(E_{n}) =
\frac{1}{Z} \sum_{n=0}^{\infty} nh\nu \ e^{-n\beta h\nu} = 
-\frac{1}{Z} \frac{\partial Z}{\partial \beta} = 
-\frac{\partial \log Z}{\partial \beta}.

Substituting the formula for Z gives the final result[7]


\langle H \rangle = h\nu \frac{e^{-\beta h\nu}}{1 - e^{-\beta h\nu}}.

At high temperatures, when the thermal energy kBT is much greater than the spacing between energy levels, the exponential argument βhν is much less than one and the average energy becomes kBT, in agreement with the equipartition theorem (Figure 10). However, at low temperatures, when >> kBT, the average energy goes to zero — the higher-frequency energy levels are "frozen out" (Figure 10). As another example, the internal excited electronic states of a hydrogen atom do not contribute to its specific heat as a gas at room temperature, since the thermal energy kBT (roughly 0.025 eV) is much smaller than the spacing between the lowest and next higher electronic energy levels (roughly 10 eV).

Similar considerations apply whenever the energy level spacing is much larger than the thermal energy. For example, this reasoning was used by Albert Einstein to resolve the ultraviolet catastrophe of blackbody radiation.[48] The paradox arises because there are an infinite number of independent modes of the electromagnetic field in a closed container, each of which may be treated as a harmonic oscillator. If each electromagnetic mode were to have an average energy kBT, there would be an infinite amount of energy in the container.[48][49] However, by the reasoning above, the average energy in the higher-ω modes goes to zero as ω goes to infinity; moreover, Planck's law of black body radiation, which describes the experimental distribution of energy in the modes, follows from the same reasoning.[48]

Other, more subtle quantum effects can lead to corrections to equipartition, such as identical particles and continuous symmetries. The effects of identical particles can be dominant at very high densities and low temperatures. For example, the valence electrons in a metal can have a mean kinetic energy of a few electronvolts, which would normally correspond to a temperature of tens of thousands of kelvins. Such a state, in which the density is high enough that the Pauli exclusion principle invalidates the classical approach, is called a degenerate fermion gas. Such gases are important for the structure of white dwarf and nuetron stars. At low temperatures, a fermionic analogue of the Bose–Einstein condensate (in which a large number of identical particles occupy the lowest-energy state) can form; such superfluid electrons are responsible for superconductivity.

[edit] See also

[edit] Notes and references

  1. ^ a b Kundt, A; Warburg E (1876). "Über die specifische Wärme des Quecksilbergases (On the specific heat of mercury gases)". Annalen der Physik 157: 353–369.  (German)
  2. ^ Fact Sheet on Uranium Enrichment U.S. Nuclear Regulatory Commission. Accessed 30th April 2007
  3. ^ a b c d e f g h i j k l Pathria, RK (1972). Statistical Mechanics. Pergamon Press, pp. 43–48, 73–74. ISBN 0-08-016747-0. 
  4. ^ Cavanagh, J; Fairbrother WJ, Palmer AG III, Skelton NJ, Rance M (2006). Protein NMR Spectroscopy: Principles and Practice, 2nd ed., Academic Press. ISBN 978-0121644918. 
  5. ^ Cantor, CR; Schimmel PR (1980). Biophysical Chemistry. Part II. Techniques for the study of biological structure and function. W. H. Freeman. ISBN 978-0716711896. 
  6. ^ a b c Goldstein, H (1980). Classical Mechanics, 2nd. ed, Addison-Wesley. ISBN 0-201-02918-9. 
  7. ^ a b c d e f g h i Huang, K (1987). Statistical Mechanics, 2nd ed., John Wiley and Sons, pp. 136–138. ISBN 0-471-81518-7. 
  8. ^ a b Mandl, F (1971). Statistical Physics. John Wiley and Sons, pp. 213–219. ISBN 0-471-56658-6. 
  9. ^ a b c d Pais, A (1982). Subtle is the Lord. Oxford University Press. ISBN 0-19-853907-X. 
  10. ^ a b Tolman, RC (1918). "A General Theory of Energy Partition with Applications to Quantum Theory". Physical Review 11: 261–275. 
  11. ^ Miedl M, Garcia M, Bamforth C (2005). "Haze formation in model beer systems". J. Agric. Food Chem. 53 (26): 10161-5. doi:10.1021/jf0506941. PMID 16366710. 
  12. ^ Mason, M; Weaver W (1924). "The Settling of Small Particles in a Fluid". Physical Review 23: 412–426. 
  13. ^ Brush, SG (1976). The Kind of Motion We Call Heat, Volume 1. Amsterdam: North Holland, 134–159. ISBN 978-0444870094. 
    Brush, SG (1976). The Kind of Motion We Call Heat, Volume 2. Amsterdam: North Holland, 336–339. ISBN 978-0444870094. 
    Waterston, JJ (1846/1893). "On the physics of media that are composed of free and elastic molecules in a state of motion". Roy. Soc. Proc. 5: 604. (abstract only). Not published in full until "" (1893). Philos. Trans. R. Soc. London A183: 1–79.  Reprinted (1928) in J.S. Haldane: The collected scientific papers of John James Waterston. Edinburgh: Oliver & Boyd. 
    Waterston, JJ (1843). Thoughts on the Mental Functions.  (reprinted in his Papers, 3, 167, 183.)
    Waterston, JJ (1851). "". British Association Reports 21: 6.  Waterston's key paper was written and submitted in 1845 to the Royal Society. After refusing to publish his work, the Society also refused to return his manuscript and stored it among its files. The manuscript was discovered in 1891 by Lord Rayleigh, who criticized the original reviewer for failing to recognize the significance of Waterston's work. Waterston managed to publish his ideas in 1851, and therefore has priority over Maxwell for enunciating the first version of the equipartition theorem.
  14. ^ Maxwell, JC (2003). "Illustrations of the Dynamical Theory of Gases", in WD Niven: The Scientific Papers of James Clerk Maxwell. New York: Dover, Vol.1, pp. 377–409. ISBN 978-0486495606.  Read by Prof. Maxwell at a Meeting of the British Association at Aberdeen on 21 September 1859.
  15. ^ Boltzmann, L (1871). "Einige allgemeine Sätze über Wärmegleichgewicht (Some general statements on thermal equilibrium)". Wiener Berichte 63: 679–711.  (German) In this preliminary work, Boltzmann showed that the average total kinetic energy equals the average total potential energy when a system is acted upon by external harmonic forces.
  16. ^ Boltzmann, L (1876). "Über die Natur der Gasmoleküle (On the nature of gas molecules)". Wiener Berichte 74: 553–560.  (German)
  17. ^ a b McQuarrie, DA (2000). Statistical Mechanics, revised 2nd ed., University Science Books, pp. 91–128. ISBN 978-1891389153. 
  18. ^ Petit, AT; Dulong PL (1819). "Recherches sur quelques points importants de la théorie de la chaleur (Studies on key points in the theory of heat)". Annales de Chimie et de Physique 10: 395–413.  (French)
  19. ^ Dewar, J (1872). "The Specific Heat of Carbon at High Temperatures". Philosophical Magazine 44: 461. 
    Weber, HF (1872). "Die specifische Wärme des Kohlenstoffs (The specific heat of carbon)". Annalen der Physik 147: 311–319.  (German)
    Weber, HF (1875). "Die specifische Wärmen der Elemente Kohlenstoff, Bor und Silicium (The specific heats of elemental carbon, boron, and silicon)". Annalen der Physik 154: 367–423, 553–582.  (German)
  20. ^ de la Rive, A; Marcet F (1840). "Quelques recherches sur la chaleur spécifique (Some research on specific heat)". Annales de Chimie et de Physique 75: 113–144.  (French)
    Regnault, HV (1841). "Recherches sur la chaleur spécifique des corps simples et des corps composés (deuxième Mémoire) (Studies of the specific heats of simple and composite bodies)". Annales de Chimie et de Physique 1 (3me Série): 129–207.  (French) Read at l'Académie des Sciences on 11 January 1841.
    Wigand, A (1907). "Über Temperaturabhängigkeit der spezifischen Wärme fester Elemente (On the temperature dependence of the specific heats of solids)". Annalen der Physik 22: 99–106.  (German)
  21. ^ a b Wüller, A (1896). Lehrbuch der Experimentalphysik (Textbook of Experimental Physics). Leipzig: Teubner, Vol. 2, 507ff.  (German)
  22. ^ a b Eucken, A (1912). "Die Molekularwärme des Wasserstoffs bei tiefen Temperaturen (The molecular specific heat of hydrogen at low temperatures)". Sitzungsberichte der königlichen Preussischen Akademie der Wissenschaften 1912: 141–151.  (German)
  23. ^ a b Maxwell, JC (1890). "On the Dynamical Evidence of the Molecular Constitution of Bodies", in WD Niven: The Scientific Papers of James Clerk Maxwell. Cambridge: At the University Press, Vol.2, pp.418–438. ASIN B000GW7DXY.  A lecture delivered by Prof. Maxwell at the Chemical Society on 18 February 1875.
  24. ^ a b Kittel, C (1996). Introduction to Solid State Physics. New York: John Wiley and Sons, 151–156. ISBN 978-0471111818. 
  25. ^ Boltzmann, L (1895). "On certain Questions of the Theory of Gases". Nature 51: 413–415. 
  26. ^ Thomson, W (1904). Baltimore Lectures. Baltimore: Johns Hopkins University Press, Sec. 27.  Re-issued in 1987 by MIT Press as Kelvin's Baltimore Lectures and Modern Theoretical Physics: Historical and Philosophical Perspectives (Robert Kargon and Peter Achinstein, editors). ISBN 978-0262111171
  27. ^ Rayleigh, JWS (1900). "The Law of Partition of Kinetic Energy". Philosophical Magazine 49: 98–118. 
  28. ^ Einstein, A (1907). "Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme (The Planck theory of radiation and the theory of specific heat)". Annalen der Physik 22: 180–190. (German)
    Einstein, A (1907). "Berichtigung zu meiner Arbeit: 'Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme' (Correction to previous article)". Annalen der Physik 22: 800.  (German)
    Einstein, A (1911). "Eine Beziehung zwischen dem elastischen Verhalten and der spezifischen Wärme bei festen Körpern mit einatomigem Molekül (A connection between the elastic behavior and the specific heat of solids with single-atom molecules)". Annalen der Physik 34: 170–174.  (German)
    Einstein, A (1911). "Bemerkung zu meiner Arbeit: 'Eine Beziehung zwischen dem elastischen Verhalten and der spezifischen Wärme bei festen Körpern mit einatomigem Molekül' (Comment on previous article)". Annalen der Physik 34: 590.  (German)
    Einstein, A (1911). "Elementare Betrachtungen über die thermische Molekularbewegung in festen Körpern (Elementary observations on the thermal movements of molecules in solids)". Annalen der Physik 35: 679–694.  (German)
  29. ^ Nernst, W (1910). "Untersuchungen über die spezifische Wärme bei tiefen Temperaturen. II. (Investigations into the specific heat at low temperatures)". Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften 1910: 262–282.  (German)
  30. ^ Hermann, Armin (1971). The Genesis of Quantum Theory (1899–1913), original title: Frühgeschichte der Quantentheorie (1899–1913), translated by Claude W. Nash, Cambridge, MA: The MIT Press, pp. 124–145. ISBN 0-262-08047-8, LCCN 73-151106. 
  31. ^ a b c Tolman, RC (1938). The Principles of Statistical Mechanics. New York: Dover Publications, pp. 93–98. ISBN 0-486-63896-0. 
  32. ^ Clausius, R (1870). "Ueber einen auf die Wärme anwendbaren mechanischen Satz". Annalen der Physik 141: 124–130.  (German)
    Clausius, RJE (1870). "On a Mechanical Theorem Applicable to Heat". Philosophical Magazine, Ser. 4 40: 122–127. 
  33. ^ a b L. Vu-Quoc, Configuration integral (statistical mechanics), 2008.
  34. ^ McQuarrie, DA (2000). Statistical Mechanics, revised 2nd ed., University Science Books, pp. 254–264. ISBN 978-1891389153. 
  35. ^ a b Tolman, RC (1927). Statistical Mechanics, with Applications to Physics and Chemistry. Chemical Catalog Company, pp. 76–77. 
  36. ^ Terletskii, YP (1971). Statistical Physics, translated: N. Fröman, Amsterdam: North-Holland, pp. 83–84. ISBN 0-7204-0221-2, LCCN 70-157006. 
  37. ^ Collins, GW (1978). The Virial Theorem in Stellar Astrophysics. Pachart Press. 
  38. ^ Chandrasekhar, S (1939). An Introduction to the Study of Stellar Structure. Chicago: University of Chicago Press, pp. 49–53. 
  39. ^ Kourganoff, V (1980). Introduction to Advanced Astrophysics. Dordrecht, Holland: D. Reidel, pp. 59–60, 134–140, 181–184. 
  40. ^ Chiu, H-Y (1968). Stellar Physics, volume I. Waltham, MA: Blaisdell Publishing. LCCN 67-17990. 
  41. ^ Noyes, RW (1982). The Sun, Our Star. Cambridge, MA: Harvard University Press. ISBN 0-674-85435-7. 
  42. ^ Ostlie, DA; Carroll BW (1996). An Introduction to Modern Stellar Astrophysics. Reading, MA: Addison-Wesley. ISBN 0-201-59880-9. 
  43. ^ Jeans, JH (1902). "The Stability of a Spherical Nebula". Phil.Trans. A 199: 1–53. 
  44. ^ a b McQuarrie, DA (2000). Statistical Mechanics, revised 2nd ed., University Science Books, pp. 121–128. ISBN 978-1891389153. 
  45. ^ Callen, HB (1985). Thermodynamics and an Introduction to Thermostatistics. New York: John Wiley and Sons, pp. 375–377. ISBN 0-471-86256-8. 
  46. ^ Arnold, VI; Avez A (1967). Théorie ergodique des systèms dynamiques. Gauthier-Villars, Paris. (French) (English edition: Benjamin-Cummings, Reading, Mass. 1968). 
  47. ^ a b Reichl, LE (1998). A Modern Course in Statistical Physics, 2nd ed., Wiley Interscience, 326–333. ISBN 978-0471595205. 
  48. ^ a b c Einstein, A (1905). "Über einen die Erzeugung und Verwandlung des Lichtes betreffenden heuristischen Gesichtspunkt (A Heuristic Model of the Creation and Transformation of Light)". Annalen der Physik 17: 132–148.  (German). An English translation is available from Wikisource.
  49. ^ Rayleigh, JWS (1900). "Remarks upon the Law of Complete Radiation". Philosophical Magazine 49: 539–540. 

[edit] Further reading

  • Mandl, F (1971). Statistical Physics. John Wiley and Sons, pp. 213–219. ISBN 0-471-56658-6. 
  • Mohling, F (1982). Statistical Mechanics: Methods and Applications. John Wiley and Sons, pp. 137–139, 270–273, 280, 285–292. ISBN 0-470-27340-2. 
  • Pathria, RK (1972). Statistical Mechanics. Pergamon Press, pp. 43–48, 73–74. ISBN 0-08-016747-0. 
  • Tolman, RC (1927). Statistical Mechanics, with Applications to Physics and Chemistry. Chemical Catalog Company, pp. 72–81.  ASIN B00085D6OO

[edit] External links