San José State University

applet-magic.com
Thayer Watkins
Silicon Valley
USA

 The Virial Theorem and the Matter of the Origin of the Position Vectors of the Particles

The Virial Theorem is an important theorem in mechanics. It is marvelously simple to prove. In the material below the variables which are vectors will be displayed in red.

Consider a system of n point particles indexed by i. Let ri, vi and pi be the position, velocity and momentum vectors, respectively, for the i-th particle. Its mass is denoted as mi. Let Fi be the net force, internal and external, impinging upon the i-th particle.

• Statement of the Virial Theorem: For the n point particles bound together into a system the time average of the kinetic energy of the particles, Σ½mivi², plus one half of the time average of ΣFi·ri, the virial of the particles, is equal to zero.

## Proof

For each particle

#### vi = dri/dtpi = mividpi/dt = Fi

Define H as Σpi·ri. (Note that the dot product of two vectors is a scalar.) The quantity H does not have a generally accepted name, although it seems it should have. It is in the nature of a moment.

Then

#### dH/dt = Σ(dpi/dt)·ri + Σ pi·(dri/dt)

Because (dpi/dt) = Fi the first term on the right in the above equation reduces to ΣFi·ri. This quantity was called the virial of the system by Clausius in 1870, a term he coined based upon the Latin term vis, meaning force or energy.

Because pi=mivi =mi(dri/dt) the second term on the right in the previous equation reduces to

#### Σ pi·(dri/dt) = Σ mi(dri/dt)² =Σmivi²

This last expression is just twice the kinetic energy K of the system; i.e., 2Σ½mivi². Thus

## The Time Average of the Variables

The time average of a variable y(t) over the interval 0 to τ is defined as

#### y = (1/τ)∫0τ y(t)dt

Time averaging the equation of dH/dt gives

#### dH/dt = ΣFi·ri + 2K.

The time average of (dH/dt) is just

#### dH/dt = (1/τ)[H(τ)−H(0)]

If the system is cyclical such that it returns to its initial state after an interval then τ can be chosen equal to the cycle period and dH/dt reduces to 0. If the system is not cyclical then for a system which is bounded the limit of dH/dt as τ increases without bound is zero.

Thus

## The Nature of Vectors

A vector is defined by three elements; its location, its direction and its magnitude. Vectors are often represented as arrows. The origin of the arrow specifies its location, the direction of the arrow is the direction of the vector and the length of the arrow is the magnitude of the vector. The subtle problem in the virial theorem is the specification of the origins for the vectors. In case of the forces, location of a force on a particle is simply the location of the particle. Forces cannot impinge upon the system of particles except through the particles.

The location of the origin for the position vectors, ri, is a different matter. The origin for the position vectors would seem to be entirely arbitrary. However a change in the origin will affect the virial by changing the magnitudes of the position vectors and also through changing the angles between the position vectors and the force vectors.

Suppose that for one specification of the origin of the position vectors the value of the H is determined

#### H = Σ pi·ri

Then suppose the origin of the position vectors is moved and let roo represent the movement of the origin. The new position vectors are then given by

#### r'i = ri + roo

The forces within a system typically depend only upon the distances between the particles and these distances are not affected by a change in the origin of the position vectors. Likewise the momenta and kinetic energies would seem to be not affected by a change in the origin of the position vectors. But the values of H and ΣFi ri would seem to be affected.

The new value of H' is then

Likewise

#### ΣFi·r'i = ΣFi·ri + ΣFi·roo

The origin of the position vectors could be located at the center of the Andromeda Galaxy, which is 2.5 million light years away, thus seemingly affecting the value of H and ΣFi·ri enormously .

If the change of origin were to change the H only by a constant amount then the dynamics of the system would be unaffected. However by the virial theorem the average kinetic energy of the system would change by the negative of half the amount of the change in the virial. This would seem to indicate that the origin of the position vectors cannot be chosen arbitrarily, even though the proof of the virial theorem does not seem to impose any restriction on that origin.

It is known that the virial theorem holds for the choice of the origin of the position vectors being at the center of mass of the system. To see that the virial theorem holds for other choices for the origin of the position vectors consider the case for a two particle system.

## Two Particle Systems

Consider a two-particle system held together by gravitation and that the orbits of the particles are circles with their centers at the center of mass of the particles. Furthermore assume the masses of the particles are both m. If s is the separation distance then the orbits have a radius of s/2. Then the gravitational attraction must be sufficient to generate the centripetal acceleration necessary to keep the particles in circular orbits. This requires that

#### mv²/(s/2) = Gm²/s²

where v is the tangential velocity and G is the graviational constant.

The kinetic energy of one particle is ½mv² so the total kinetic energy of the system, K, is mv². From the above relation

#### K = mv² = ½Gm²/s

The RHS of this equation happens to be one half of negative of the potential energy V of the system, so

#### K = −½V or, equivalently K + ½V = 0

The position vectors with respect to the center of mass have magnitudes equal to (s/2). They are parallel with the forces and point in the opposite direction. Thus H of the system is given by

#### ΣFiri = −2(Gm²/s²)(s/2) = −Gm²/s

Thus the virial theorem is exactly satisfied.

Now assume the origin of the coordinate system is at the center of one of the particles, say Particle 1. Then r1 is equal to zero and r2 is equal in magnitude to the separation distance for Particle 2.

The force between the particles is given by

#### Fi = −Gm1m2/|r2-r1|² which reduces to Fi = −Gm12/|r2|² = −Gm²/s²

The value of H is then given by

#### H = −(Gm²/s²)s = −Gm²/s

Again the virial theorem is exactly satisfied.

## The Effect of a Change in the Origin of the Coordinate System

Consider the position vectors for one origin r and another origin r'. (Here the prime (') denotes a different case rather than differentiation.) Then

#### r'(t) = r(t) + roo(t)

where roo(t) is the vector between the origin of the first coordinate system and the origin of the second coordinate system. The velocities are then related by

#### dr'/dt = dr/dt + droo/dt or, equivalentlyv'(t) = v(t) + voo(t)

If the origins of the position vectors are motionless then v'(t)=v(t), and likewise the momenta vectors, p'(t)=p(t). However this is not the case with the value of p(t)·r(t).

#### p'(t)r'(t) = p(t)·r(t) + p(t)·roo(t)

For forces depending only upon the separation of the particles F'(t)=F(t). But

#### F'(t)·r'(t) = F(t)·r(t) + F(t)·roo(t)

For summation over all the particles of the system then

#### ΣF'i(t)·r'i(t) = ΣF(t)·r(t) + (ΣFi(t))·roo(t)

For systems subject only to internal forces the net force is zero because the forces between any two particles cancel out. (Newton's third law of mechanics.) Thus the virials for the systems under the two coordinate systems are the same. Thus if (but not only if) the origins of the coordinate systems are motionless and the particles are subject only to internal forces dependent only upon the separation of the particles the virial theorem holds for arbitrary locations of the origin of the coordinate system.

## Conclusion

For systems of particles subject only to internal forces the location of the origin of the coordinate system does not affect the value of ΣFi·ri and the virial theorem holds regardless of the location of the origin of the coordinate system.

(To be continued.)