Virial theorem

In statistical mechanics, the virial theorem provides a general equation that relates the average over time of the total kinetic energy of a stable system of discrete particles, bound by a conservative force (forces characterized exclusively by their work),^{[dubious – discuss]} with that of the total potential energy of the system. Mathematically, the theorem states

\left\langle T\right\rangle =-{\frac {1}{2}}\,\sum _{k=1}^{N}{\bigl \langle }\mathbf {F} _{k}\cdot \mathbf {r} _{k}{\bigr \rangle }

where

T

is the total kinetic energy of the

N

particles,

F k

represents the force on the

k

th particle, which is located at position

r k

, and angle brackets represent the average over time of the enclosed quantity. The word virial for the right-hand side of the equation derives from vis, the Latin word for "force" or "energy", and was given its technical definition by Rudolf Clausius in 1870.^[1]

The significance of the virial theorem is that it allows the average total kinetic energy to be calculated even for very complicated systems that defy an exact solution, such as those considered in statistical mechanics; this average total kinetic energy is related to the temperature of the system by the equipartition theorem. However, the virial theorem does not depend on the notion of temperature and holds even for systems that are not in thermal equilibrium. The virial theorem has been generalized in various ways, most notably to a tensor form.

If the force between any two particles of the system results from a potential energy $V (r) = αr n$ that is proportional to some power $n$ of the interparticle distance $r$ , the virial theorem takes the simple form

2\langle T\rangle =n\langle V_{\text{TOT}}\rangle .

Thus, twice the average total kinetic energy $⟨ T ⟩$ equals $n$ times the average total potential energy $⟨ V TOT ⟩$ . Whereas $V (r)$ represents the potential energy between two particles of distance $r$ , $V TOT$ represents the total potential energy of the system, i.e., the sum of the potential energy $V (r)$ over all pairs of particles in the system. A common example of such a system is a star held together by its own gravity, where $n$ equals −1.

History

In 1870, Rudolf Clausius delivered the lecture "On a Mechanical Theorem Applicable to Heat" to the Association for Natural and Medical Sciences of the Lower Rhine, following a 20-year study of thermodynamics. The lecture stated that the mean vis viva of the system is equal to its virial, or that the average kinetic energy is equal to 1/2 the average potential energy. The virial theorem can be obtained directly from Lagrange's identity^{[moved resource?]} as applied in classical gravitational dynamics, the original form of which was included in Lagrange's "Essay on the Problem of Three Bodies" published in 1772. Karl Jacobi's generalization of the identity to N bodies and to the present form of Laplace's identity closely resembles the classical virial theorem. However, the interpretations leading to the development of the equations were very different, since at the time of development, statistical dynamics had not yet unified the separate studies of thermodynamics and classical dynamics.^[2] The theorem was later utilized, popularized, generalized and further developed by James Clerk Maxwell, Lord Rayleigh, Henri Poincaré, Subrahmanyan Chandrasekhar, Enrico Fermi, Paul Ledoux, Richard Bader and Eugene Parker. Fritz Zwicky was the first to use the virial theorem to deduce the existence of unseen matter, which is now called dark matter. Richard Bader showed the charge distribution of a total system can be partitioned into its kinetic and potential energies that obey the virial theorem.^[3] As another example of its many applications, the virial theorem has been used to derive the Chandrasekhar limit for the stability of white dwarf stars.

Illustrative special case

Consider $N = 2$ particles with equal mass $m$ , acted upon by mutually attractive forces. Suppose the particles are at diametrically opposite points of a circular orbit with radius $r$ . The velocities are $v 1 (t)$ and $v 2 (t) = - v 1 (t)$ , which are normal to forces $F 1 (t)$ and $F 2 (t) = - F 1 (t)$ . The respective magnitudes are fixed at $v$ and $F$ . The average kinetic energy of the system in an interval of time from $t 1$ to $t 2$ is

\langle T\rangle ={\frac {1}{t_{2}-t_{1}}}\int _{t_{1}}^{t_{2}}\sum _{k=1}^{N}{\frac {1}{2}}m_{k}\left|\mathbf {v} _{k}(t)\right|^{2}dt={\frac {1}{t_{2}-t_{1}}}\int _{t_{1}}^{t_{2}}\left({\frac {1}{2}}m|\mathbf {v} _{1}(t)|^{2}+{\frac {1}{2}}m|\mathbf {v} _{2}(t)|^{2}\right)dt=mv^{2}.

Taking center of mass as the origin, the particles have positions

r 1 (t)

and

r 2 (t) = - r 1 (t)

with fixed magnitude

r

. The attractive forces act in opposite directions as positions, so

F 1 (t) \cdot r 1 (t) = F 2 (t) \cdot r 2 (t) = - Fr

. Applying the centripetal force formula

F = mv 2 / r

results in:

-{\frac {1}{2}}\sum _{k=1}^{N}{\bigl \langle }\mathbf {F} _{k}\cdot \mathbf {r} _{k}{\bigr \rangle }=-{\frac {1}{2}}(-Fr-Fr)=Fr={\frac {mv^{2}}{r}}\cdot r=mv^{2}=\langle T\rangle ,

as required. Note: If the origin is displaced then we'd obtain the same result. This is because the dot product of the displacement with equal and opposite forces

F 1 (t)

,

F 2 (t)

results in net cancellation.

Statement and derivation

Although the virial theorem depends on averaging the total kinetic and potential energies, the presentation here postpones the averaging to the last step.

For a collection of $N$ point particles, the scalar moment of inertia $I$ about the origin is defined by the equation

I=\sum _{k=1}^{N}m_{k}\left|\mathbf {r} _{k}\right|^{2}=\sum _{k=1}^{N}m_{k}r_{k}^{2}

where

m k

and

r k

represent the mass and position of the

k

th particle.

r k = | r k |

is the position vector magnitude. The scalar

G

is defined by the equation

G=\sum _{k=1}^{N}\mathbf {p} _{k}\cdot \mathbf {r} _{k}

where

p k

is the momentum vector of the

k

th particle.^[4] Assuming that the masses are constant,

G

is one-half the time derivative of this moment of inertia

{\begin{aligned}{\frac {1}{2}}{\frac {dI}{dt}}&={\frac {1}{2}}{\frac {d}{dt}}\sum _{k=1}^{N}m_{k}\mathbf {r} _{k}\cdot \mathbf {r} _{k}\\&=\sum _{k=1}^{N}m_{k}\,{\frac {d\mathbf {r} _{k}}{dt}}\cdot \mathbf {r} _{k}\\&=\sum _{k=1}^{N}\mathbf {p} _{k}\cdot \mathbf {r} _{k}=G\,.\end{aligned}}

[1]

[2]

[3]

[4]