(1/2) d 2I /dt 2 =
2K + Ω [ = 0 in hydrostatic equilibrium].
I = total moment of inertia = ∫ r 2dm
K = total kinetic energy
Ω = total gravitational potential energy.
Derived from ideal gas of individual particles:
assume a collisionless perfect gas, then for individual particle i of
mass mi, F = ma is
mid 2ri /dt 2 = − Gmi∑ j mj (ri − rj) / |ri − rj| 3 = ∑ j Fij .
Start with the time derivative:
d ∑ i mi
vi . ri /dt =
∑ i (d mi
vi /dt) . ri +
∑ i mi
vi . d ri /dt .
The left-hand side is
(1/2) d 2I /dt 2.
The second term on the right-hand side is 2K.
The first term on the right-hand side is
∑ i Fi .
ri = ∑ pairs
Fij . ri +
Fji . rj =
∑ pairs
Fij .
(ri − rj)
which from above is
− ∑ pairs Gmimj /
|ri − rj| = Ω.
2.2. Pressure
The pressure of an isotropic fluid is
P = (1/3) ∫ n(p)
pv d 3p
where n(p) is the number distribution of particles
having momentum p per unit geometrical volume per unit momentum volume.
Homework: show that 2K = 3 ∫ P dV = 3 ∫ (P /ρ) dm
For ideal gasses, integrated internal energy U = K, and integrated total energy W = U + Ω = Ω/2 .
For gasses described by a γ-law, we write
P = (γ −1)ρE,
where E is the local specific internal energy per gram of material.
Ideal gas: γ = 5/3.
Radiation: γ = 4/3.
Relativistic Fermi degeneracy: γ = 4/3.
Then 2K = 3(γ −1) ∫ E dm =
3(γ −1)U, and the virial theorem becomes
(1/2) d 2I /dt 2 = 3(γ −1)U + Ω .
Homework: write U and Ω in terms of M, R, and f for our two-layer star. Assume all the material is at the same temperature. Pretty easy, since everything can be derived from K. Now, find the Kelvin-Helmholtz time scale tKH = R / (dR /dt) by assuming L = − dW /dt = constant. It should, of course, be tKH = constant × GM 2/LR.