Gravitational binding energy

From Wikipedia, the free encyclopedia
Jump to: navigation, search

For an object consisting of loose material held together by gravity alone, the gravitational binding energy is the amount of binding energy required to pull all of that material apart, to infinity. It is also the amount of energy that is liberated (usually in the form of heat) during the accretion of such an object from material falling from infinity. An object is gravitationally bound to a massive body, if it doesn't contain enough kinetic energy to escape orbit of that massive body.

The gravitational binding energy of a system is equal to the negative of the total gravitational potential energy, considering the system as a set of small particles. For a system consisting of a celestial body and a satellite, the gravitational binding energy will have a larger absolute value than the potential energy of the satellite with respect to the celestial body, because for the latter quantity, only the separation of the two components is taken into account, keeping each intact.

For a spherical mass of uniform density, the gravitational binding energy U is given by the formula[1][2]

U = \frac{3GM^2}{5R}

where G is the gravitational constant, M is the mass of the sphere, and R is its radius. This is 80% greater than the energy required to separate to infinity the two hemispheres of the spherical mass.[citation needed]

Assuming that the Earth is a uniform sphere (which is not correct, but is close enough to get an order-of-magnitude estimate) with M = 5.97 · 1024kg and r = 6.37 · 106m, U is 2.24 · 1032J. This is roughly equal to one week of the Sun's total energy output. It is 37.5 MJ/kg, 60% of the absolute value of the potential energy per kilogram at the surface.

The actual depth-dependence of density, inferred from seismic travel times (see Adams–Williamson equation), is given in the Preliminary Reference Earth Model (PREM).[3] Using this, the real gravitational binding energy of Earth can be calculated numerically to U = 2.487 · 1032 J

According to the virial theorem, the gravitational binding energy of a star is about two times its internal thermal energy.[1]

Derivation for a uniform sphere[edit]

The gravitational binding energy of a sphere with Radius R is found by imagining that it is pulled apart by successively moving spherical shells to infinity, the outermost first, and finding the total energy needed for that.

If we assume a constant density \rho then the masses of a shell and the sphere inside it are:

m_\mathrm{shell}=4\pi r^{2}\rho\,dr      and      m_\mathrm{interior}=\frac{4}{3}\pi r^3 \rho

The required energy for a shell is the negative of the gravitational potential energy:

{\it dU}=-G\frac{m_\mathrm{shell} m_\mathrm{interior}}{r}

Integrating over all shells we get:

U=-G\int_0^R {\frac{(4\pi r^2\rho)(\tfrac{4}{3}\pi r^{3}\rho)}{r}} dr = -G{\frac{16}{15}}{\pi}^2{\rho}^2 R^5

Remembering that \rho is simply equal to the mass of the whole divided by its volume for objects with uniform density we get:

\rho=\frac{M}{\frac{4}{3}\pi R^3}

And finally, plugging this into our result we get:

U=-G\frac{16}{15} \pi^2 R^5 \left(\frac{M}{\frac{4}{3}\pi R^3}\right)^2= \frac{-3GM^2}{5R}

Non-uniform spheres[edit]

Planets and stars have radial density gradients from their lower density surfaces to their much larger density compressed cores. Degenerate matter objects (white dwarfs; neutron star pulsars) have radial density gradients plus relativistic corrections.

Neutron star relativistic equations of state provided by Jim Lattimer include a graph of radius vs. mass for various models.[4] The most likely radii for a given neutron star mass are bracketed by models AP4 (smallest radius) and MS2 (largest radius). BE is the ratio of gravitational binding energy mass equivalent to observed neutron star gravitational mass of "M" kilograms with radius "R" meters,

BE = \frac{0.60\,\beta}{1 - \frac{\beta}{2}}      \beta \ = G\,M/R\,{c}^{2}

Given current values

G = 6.6742\times10^{-11}\, m^3kg^{-1}sec^{-2} [5]
c^2 = 8.98755\times10^{16}\, m^2sec^{-2}
M_{solar} = 1.98844\times10^{30}\, kg

and star masses "M" commonly reported as multiples of one solar mass,

M_x = \frac{M}{M_\odot}

then the relativistic fractional binding energy of a neutron star is

BE = \frac{885.975\,M_x}{R - 738.313\,M_x}

See also[edit]


  1. ^ a b Chandrasekhar, S. 1939, An Introduction to the Study of Stellar Structure (Chicago: U. of Chicago; reprinted in New York: Dover), section 9, eqs. 90-92, p. 51 (Dover edition)
  2. ^ Lang, K. R. 1980, Astrophysical Formulae (Berlin: Springer Verlag), p. 272
  3. ^ Dziewonski, A. M.; Anderson, D. L.. "Preliminary Reference Earth Model". Physics of the Earth and Planetary Interiors 25: 297–356. Bibcode:1981PEPI...25..297D. doi:10.1016/0031-9201(81)90046-7. 
  4. ^ Neutron Star Masses and Radii, p. 9/20, bottom
  5. ^ Measurement of Newton's Constant Using a Torsion Balance with Angular Acceleration Feedback , Phys. Rev. Lett. 85(14) 2869 (2000)