virial theorem | Physics Library

We start with the moment of inertia about the origin for the system of particles, which is defined as

∑ 2 I = miri i

(1)

Differentiate using the chain rule. Note that a vector dotted into itself yields its magnitude square.

⃗ri ⋅ ⃗ri = r2i

(2)

This lets us make the connection that

∑ ∑ I = miri2 = mi(⃗ri ⋅ ⃗ri) i i

(3)

So after differentiating we get

dI ∑ d⃗ri ⃗dri ---= mi ---⋅ ⃗ri + ⃗ri ⋅-- dt i dt dt

(4)

Differentiating again yields

d2I- ∑ d2⃗ri ⃗dri ⃗dri ⃗dri ⃗dri d2⃗ri dt2 = mi dt2 ⋅ ⃗ri + dt ⋅ dt + dt ⋅ dt + r⃗i ⋅dt2 i

(5)

In short form

2 ∑ d-I-= 2 (mi⃗r˙i ⋅ ˙⃗ri + mi ⃗ri ⋅ ¨⃗ri) dt2 i

(6)

When dealing with a system of particles we found that the kinetic energy associated with a system of particles was

1∑ T = -- mi ˙⃗ri ⋅ ˙⃗ri 2 i

(7)

Plugging in T into (6) gives us

2 ∑ d-I-= 4T + 2 m ⃗r ⋅ ¨⃗r dt2 ii i i

(8)

Next we need to tackle the 2 ∑ _im_ir_i ⋅r_i term. We first bring in the potential energy through its connection with force. This part is named the Virial of Claussius.

f⃗i = mi ¨⃗ri = ∇iU

(9)

This gives us the equality

∑ ∑ 2 mi ⃗ri ⋅⃗r¨i = ⃗ri ⋅ ∇iU i i

(10)

Now due to Newton’s 3rd law that states for every action there is an opposite and equal reaction we have the forces on the ith particle in our system given by

⃗ ∑ ⃗ fi = Fij wherei ⁄= j j

(11)

So when we go to sum up all the forces we notice ’force pairing’ such that

∑ ⃗ ∑ ∑ ⃗ ⃗ fi = Fij + Fji i i j>i

(12)

Plugging this into the virial of Claussius yields

∑ ∑ ∑ f⃗i ⋅ ⃗ri = (F⃗ij ⋅ ⃗ri + F⃗ji ⋅ ⃗rj) i i j>i

(13)

Now we need to take the gradient of the potential energy to get the force. This is a tedious calculation which can be found here (insert link)

Gmimj Gmimj (⃗ri − ⃗rj) F⃗ij = − ∇i --r-----= − -------r3------- ij ij

(14)

Inserting this into (13) gives

∑ ∑ ∑ Gm m (⃗r − ⃗r) ⋅ ⃗r Gm m (⃗r − ⃗r) ⋅ ⃗r ⃗fi ⋅ ⃗ri = (− ----i-j--i----j---i− ----i-j--j----i---j i i j>i r3ij r3ij

(15)

The two keys in understanding the above equation is to note that

rij = |⃗ri − ⃗rj|

(16)

and that

|⃗ri − ⃗rj|2 = (⃗ri − ⃗rj) ⋅ (⃗ri − r⃗j)

(17)

So now we add the numerators of (15) to get

(⃗ri − ⃗rj) ⋅ ⃗ri − (⃗rj − ⃗ri) ⋅ ⃗rj = (r⃗i ⋅ ⃗ri − ⃗ri ⋅ ⃗rj + ⃗rj ⋅ ⃗rj − ⃗ri ⋅ ⃗rj)

(18)

Next we expand (17) to see that it is equal to (18) and this cancels a power of 2 in the denominator of (15) to finally yeild the expression for potential energy

∑ ∑ Gmimj − ---r----= − U i j>i ij

(19)

So going back to (8) we see

d2I --2-= 4T − 2U dt

(20)

From the energy equation of the system, E = T - U and the key to finish up the virial theorem is to note that the momemt of inertia does not change on average with time (”After one dynamical timescale, the time derivative of I is constant so the second derivative is zero

d2I --2-= 0 dt

(21)

This leads to

1- T = 2 U

(22)

plugging this into the energy equation gives us the result of the virial theorem which states that the total energy of a stationary system (no significant dynamical evolution) is one-half the potential energy of the system.

E = 1-U − U = − 1U 2 2

(23)

"virial theorem" is owned by bloftin.
(view preamble)