Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random | Template Test |  
Login
create new user
Username:
Password:
forget your password?
Main Menu
Sections

Meta

Talkback

Downloads

Information
virial theorem (Definition)

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

$\displaystyle I = \sum_i m_i {r_i}^2$ (1)

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

$\displaystyle \vec{r_i} \cdot \vec{r_i} = r_i^2$ (2)

This lets us make the connection that

$\displaystyle I = \sum_i m_i {r_i}^2 = \sum_i m_i (\vec{r_i} \cdot \vec{r_i})$ (3)

So after differentiating we get

$\displaystyle \frac{dI}{ dt} = \sum_i m_i \frac{\vec{dr_i}}{dt} \cdot \vec{r_i} + \vec{r_i} \cdot \frac{\vec{dr_i}}{dt}$ (4)

Differentiating again yields

$\displaystyle \frac{d^2I}{dt^2} = \sum_i m_i \frac{d^2\vec{r_i}}{dt^2} \cdot \v... ...i}}{dt} \cdot \frac{\vec{dr_i}}{dt} + \vec{r_i} \cdot \frac{d^2\vec{r_i}}{dt^2}$ (5)

In short form

$\displaystyle \frac{d^2I}{dt^2} = 2\, \sum_i (m_i \dot{\vec{r_i}} \cdot \dot{\vec{r_i}} + m_i\vec{r_i} \cdot \ddot{\vec{r_i}})$ (6)

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

$\displaystyle T = \frac{1}{2} \sum_i m_i \dot{\vec{r_i}} \cdot \dot{\vec{r_i}}$ (7)

Plugging in T into (6) gives us

$\displaystyle \frac{d^2I}{dt^2} = 4T + 2 \, \sum_i m_i\vec{r_i} \cdot \ddot{\vec{r_i}}$ (8)

Next we need to tackle the $2 \, \sum_i m_i\vec{r_i} \cdot \ddot{\vec{r_i}}$ term. We first bring in the potential energy through its connection with force. This part is named the Virial of Claussius.

$\displaystyle \vec{f_i} = m_i \ddot{\vec{r_i}} = \nabla_i U$ (9)

This gives us the equality

$\displaystyle 2 \, \sum_i m_i\vec{r_i} \cdot \ddot{\vec{r_i}} = \sum_i \vec{r_i} \cdot \nabla_i U$ (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

$\displaystyle \vec{f_i} = \sum_j \vec{F_ij} \: \: \: where i \ne j$ (11)

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

$\displaystyle \sum_i \vec{f_i} = \sum_i \sum_{j>i} \vec{F_{ij}} + \vec{F_{ji}}$ (12)

Plugging this into the virial of Claussius yields

$\displaystyle \sum_i \vec{f_i} \cdot \vec{r_i} = \sum_i \sum_{j>i} (\vec{F_{ij}} \cdot \vec{r_i} + \vec{F_{ji}} \cdot \vec{r_j})$ (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)

$\displaystyle \vec{F_{ij}} = -\nabla_i \frac{G m_i m_j}{r_{ij}} = - \frac{G m_i m_j(\vec{r_i} - \vec{r_j})}{r_{ij}^3}$ (14)

Inserting this into (13) gives

$\displaystyle \sum_i \vec{f_i} \cdot \vec{r_i} = \sum_i \sum_{j>i} (- \frac{G m... ...}{r_{ij}^3} - \frac{G m_i m_j(\vec{r_j} - \vec{r_i}) \cdot \vec{r_j}}{r_{ij}^3}$ (15)

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

$\displaystyle r_{ij} = \vert \vec{r_i} - \vec{r_j} \vert$ (16)

and that

$\displaystyle \vert \vec{r_i} - \vec{r_j} \vert^2 = (\vec{r_i} - \vec{r_j}) \cdot (\vec{r_i} - \vec{r_j})$ (17)

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

$\displaystyle (\vec{r_i} - \vec{r_j}) \cdot \vec{r_i} - (\vec{r_j} - \vec{r_i})... ...c{r_i} \cdot \vec{r_j} + \vec{r_j} \cdot \vec{r_j} - \vec{r_i} \cdot \vec{r_j})$ (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

$\displaystyle \sum_i \sum_{j>i} - \frac{G m_i m_j}{r_{ij}} = -U$ (19)

So going back to (8) we see

$\displaystyle \frac{d^2I}{dt^2} = 4T - 2U$ (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

$\displaystyle \frac{d^2I}{dt^2} = 0$ (21)

This leads to

$\displaystyle T = \frac{1}{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.

$\displaystyle E = \frac{1}{2}U - U = - \frac{1}{2}U$ (23)



"virial theorem" is owned by bloftin.
View style:

Cross-references: power, gradient, Newton's 3rd law, energy, kinetic energy, square, magnitude, vector, system, moment of inertia

This is version 9 of virial theorem, born on 2004-11-19, modified 2005-08-13.
Object id is 19, canonical name is VirialTheorem.
Accessed 979 times total.

Classification:
Physics Classification: 45.40.-f (Dynamics and kinematics of rigid bodies)
Pending Errata and Addenda
None.
Discussion
Style: Expand: Order:

No messages.

Testing some escape charachters for html category with a generator has an injective cogenerator" now escape ” with "