Reference: Daniel V. Schroeder, An Introduction to Thermal Physics, (Addison-Wesley, 2000) – Problem 1.55.
Suppose we have two identical masses in a circular orbit of radius about their centre of mass. By equating centripetal and gravitational forces, we have
The total kinetic and potential energies of the system are
Therefore from 1
This is a special case of the virial theorem for gravitational orbits, which gives a relation between the average potential and kinetic energies:
The total energy of a gravitational system is therefore
The total energy is negative, which indicates that it is gravitationally bound and won’t fly apart with time. However, this has a curious consequence in that, if we increase the energy of the system by an amount such that is still negative, the kinetic energy must actually decrease.
For a system such as a star that is bound by gravitational forces and contains many particles, we can apply the equipartition theorem. Because of the high temperature within a star, the atomic nuclei become dissociated from their electrons, so the only degrees of freedom available to each particle are the translational degrees of freedom, meaning that the average kinetic energy of each particle is . Therefore, for a system with particles
This means that the heat capacity is actually negative:
[Note that for a star, since even if the star’s volume changes as energy is added to it, no expansion work is done since the star expands into a vacuum where the pressure is zero. However, gravitational work is done, so I’m not sure how that will affect the formula.]
At this point, Schroeder asks us to use dimensional analysis to get a formula for the potential energy of a star. I’m not quite sure what he means, but we can get an estimate of the potential energy as follows.
Suppose the star has a uniform density and a radius . Then the potential energy of a thin shell of the star at radius is determined by the portion of the star inside this radius (the parts of the star outside have no net force on anything inside; the proof of this is similar to that in electrostatics using Gauss’s law since gravity is also an inverse square force). The portion of the star inside the radius acts as a point mass at the centre of the star. That is
We can integrate this to get :
The average kinetic energy is therefore
We can use this to get an estimate of the temperature of the Sun, whose mass is and radius is . Taking the Sun to be composed of equal numbers of bare protons and electrons (and neglecting the mass of the electron compared that of the proton), we can estimate :
The temperature estimate is
This is within currently accepted values for the interior of the sun. The core of the Sun is estimated to be around , but drops to between and as we get beyond half the Sun’s radius from the core.