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.

Pingback: Thermal equilibrium from entropy plots | Physics pages

Pingback: Entropy of a star | Physics pages

SamYour worked out problems from Schroeder’s Thermal Physics is much appreciated. The open-ended, approximation based questions can make it tough to figure out if I’m on the right track for a problem, but your website has helped a ton.

Thanks for doing all of these problems.

Juan Ramón González ÁlvarezI am sorry to say that negative heat capacities don’t exist. As showed by several authors including myself, those negative heat capacities follow from several mistakes. Check the below article blog article and the paper of Michaelian and Santamaría-Holek for a for mistakes made on gravitational systems

http://www.juanrga.com/2017/02/there-are-no-negative-heat-capacities.html