Fact-checked by Grok 2 weeks ago

Gravitational binding energy

Gravitational binding energy is the minimum energy required to disassemble a gravitationally bound system, such as a planet, star, or galaxy cluster, into its individual components dispersed to infinite separation, overcoming the attractive forces of gravity.^[1] For a uniform-density sphere of mass M and radius R, this energy is given by E_b = \frac{3}{5} \frac{G M^2}{R}, where G is the gravitational constant, representing the magnitude of the system's gravitational self-potential energy.^[2] In astrophysics, gravitational binding energy plays a central role in understanding the stability and evolution of celestial bodies through the virial theorem, which equates twice the total kinetic energy of a self-gravitating system in equilibrium to the negative of its gravitational potential energy, implying that internal thermal energy balances roughly half the binding energy.^[3] For stars, this relation governs hydrostatic equilibrium, where contraction releases binding energy that heats the core to initiate nuclear fusion, sustaining the star against further collapse.^[4] In the context of stellar death, the formation of a neutron star core releases gravitational binding energy on the order of $10^{53} ergs, mostly carried away by neutrinos, with a small fraction powering the supernova explosion whose kinetic energy (\sim 10^{51} ergs) must exceed the binding energy of the stellar envelope to unbind and disperse it.^[5]^[6] Notable examples include Earth's binding energy of approximately $2.24 \times 10^{32} joules, equivalent to about one week of the Sun's total luminosity output, and the Sun's own binding energy of approximately $2.3 \times 10^{41} joules, which underscores the immense scales involved in cosmic structures.^[7]

Basic Concepts

Definition

Gravitational binding energy U is the negative of the gravitational self-energy of a system, representing the minimum energy required to overcome the mutual gravitational attraction and disperse all mass elements to infinite separation while leaving them with zero kinetic energy.^[8] This quantity quantifies the work needed to disassemble a gravitationally bound object, such as a planet or star, into unbound particles at rest at infinity.^[9] In Newtonian gravity, the gravitational binding energy for a discrete system of point masses m_i at positions \mathbf{r}_i is given by

U = -\frac{G}{2} \sum_{i \neq j} \frac{m_i m_j}{r_{ij}},

where G is the gravitational constant and r_{ij} = |\mathbf{r}_i - \mathbf{r}_j| is the distance between masses i and j.^[8] For a continuous mass distribution with density \rho(\mathbf{r}), it takes the form

U = -\frac{1}{2} \int \rho(\mathbf{r}) \Phi(\mathbf{r}) \, dV,

where \Phi(\mathbf{r}) is the gravitational potential due to the mass distribution itself, satisfying Poisson's equation \nabla^2 \Phi = 4\pi G \rho.^[8] The factor of \frac{1}{2} in both expressions accounts for the double-counting of pairwise gravitational interactions in the summation or integration; without it, each interaction between a pair of mass elements would be included twice.^[8] These formulations assume a Newtonian framework, where gravity is instantaneous and described by inverse-square attraction.^[8] The binding energy is typically expressed in units of joules (SI) or ergs (cgs); for example, Earth's gravitational binding energy is approximately $2.24 \times 10^{32} J.^[1]

Physical Significance

The negative value of gravitational binding energy signifies that a system is gravitationally bound, meaning its total mechanical energy is less than zero when referenced to the state of infinite dispersion where particles are at rest relative to each other at infinite separation.^[10] This negative energy reflects the work done by gravity in assembling the system, resulting in a stable configuration where escape requires additional positive energy input to reach zero total energy.^[10] In contrast, a positive total energy would indicate an unbound system prone to dispersal.^[10] The magnitude of this binding energy sets a critical threshold for system disruption: external energy input must surpass the absolute value of the binding energy to disassemble the structure into unbound components.^[11] For instance, in asteroid collisions, impact kinetic energy exceeding the gravitational binding energy leads to catastrophic fragmentation, as modeled in simulations of gravity-dominated bodies. Similarly, tidal forces from a nearby massive body can disrupt a satellite if the work done by those forces exceeds the binding energy, causing the object to break apart along Roche limits. Unlike electromagnetic or nuclear binding energies, which diminish or remain constant with increasing system size and dominate at atomic or subatomic scales, gravitational binding energy scales favorably with mass and radius—for fixed density, it grows proportionally to M^{5/3}—rendering larger bodies progressively more resistant to disruption.^[12] This inverse scaling relative to other forces underscores gravity's role in stabilizing macroscopic structures like planets and stars, where electromagnetic and nuclear contributions are negligible.^[12] The concept of gravitational binding energy was formalized in the mid-19th century through the work of William Thomson (Lord Kelvin) and Hermann von Helmholtz, who invoked gravitational contraction as the primary energy source for the Sun's luminosity in the 1860s.^[13] Helmholtz's 1856 proposal and Kelvin's 1862 refinements estimated the Sun's age at around 20-30 million years based on the release of this binding energy during slow contraction, bridging thermodynamics and astronomy until nuclear processes were identified.^[13]

Newtonian Derivations

Uniform Sphere

The gravitational binding energy of a spherically symmetric body with uniform density is derived under the assumptions of Newtonian gravity, constant mass density \rho, total mass M = \frac{4}{3} \pi R^3 \rho, and spherical symmetry.^[14]^[15] One standard approach to compute the binding energy U employs the shell method, imagining the sphere assembled from infinitesimal mass shells added sequentially from the center outward.^[14] The energy contribution dU when adding a shell of mass dM at radius r to the existing mass M(r) enclosed within r is given by the potential energy of the shell in the field of M(r):

dU = -\frac{G M(r) \, dM}{r}.

For uniform density, the enclosed mass is M(r) = \frac{4}{3} \pi r^3 \rho, and the shell mass is dM = 4 \pi r^2 \, dr \, \rho.^[15] Substituting these yields

dU = -\frac{G}{r} \left( \frac{4}{3} \pi r^3 \rho \right) \left( 4 \pi r^2 \, dr \, \rho \right) = - \frac{16 \pi^2 G \rho^2}{3} r^4 \, dr.

Integrating from r = 0 to r = R,

U = -\frac{16 \pi^2 G \rho^2}{3} \int_0^R r^4 \, dr = -\frac{16 \pi^2 G \rho^2}{3} \left[ \frac{r^5}{5} \right]_0^R = -\frac{16 \pi^2 G \rho^2 R^5}{15}.

Expressing in terms of total mass M = \frac{4}{3} \pi R^3 \rho, so \rho = \frac{3M}{4 \pi R^3} and \rho^2 R^5 = \left( \frac{3M}{4 \pi R^3} \right)^2 R^5 = \frac{9 M^2}{16 \pi^2 R}, substitution gives

U = -\frac{16 \pi^2 G}{15} \cdot \frac{9 M^2}{16 \pi^2 R} = -\frac{3}{5} \frac{G M^2}{R}.

The factor of $3/5 arises from the volume average of the gravitational potential over the sphere, which for uniform density is \langle \Phi \rangle = -\frac{3}{5} \frac{G M}{R}, since the potential inside is \Phi(r) = -\frac{G M}{2 R^3} (3 R^2 - r^2), and the self-energy is half the integral of \rho \Phi to avoid double-counting pairwise interactions.^[15]^[14] As an illustrative example, approximating Earth as a uniform sphere with mass M = 5.972 \times 10^{24} kg and mean radius R = 6.371 \times 10^6 m yields U \approx -2.24 \times 10^{32} J using the formula above.^[16] This value underestimates the actual binding energy because Earth's density increases toward the core, requiring a more general integral for non-uniform cases.^[14]

Non-Uniform Density Distributions

For bodies with non-uniform density distributions, the Newtonian gravitational binding energy, also known as the gravitational self-energy, requires an integral formulation to account for the varying mass distribution within the object. Assuming spherical symmetry, the total gravitational potential energy U is given by

U = -\int_0^R \frac{G m(r)}{r} \, dm(r),

where G is the gravitational constant, R is the outer radius, m(r) = 4\pi \int_0^r \rho(s) s^2 \, ds is the mass enclosed within radius r, \rho(r) is the density as a function of radius, and dm(r) = 4\pi r^2 \rho(r) \, dr is the infinitesimal mass of a thin spherical shell at radius r.^[17] This expression arises from the pairwise gravitational interactions but simplifies under spherical symmetry by treating the assembly of concentric mass shells. The derivation proceeds by considering the incremental assembly of the body from the center outward. Each infinitesimal shell of mass dm(r) at radius r is brought from infinity and placed in the gravitational potential created by the already assembled inner mass m(r), which generates a potential \phi(r) = -[G](/page/G) m(r)/r at that location. The work done to assemble the shell is thus dU = \phi(r) \, dm(r) = - [G m(r)/r] dm(r), and integrating over the entire mass yields the total binding energy U.^[17] This shell-by-shell approach avoids the full double integral over all mass elements while correctly capturing the negative energy required to bind the system against dispersion. A particularly useful application of this formula occurs in polytropic models of stars, where the density profile follows \rho(r) \propto \theta^n from solutions to the Lane-Emden equation, with polytropic index n characterizing the equation of state P \propto \rho^{1 + 1/n}. For such configurations, the binding energy takes the scaled form

U = -\frac{3}{5 - n} \frac{G M^2}{R},

where M is the total mass and R is the radius, valid for $0 \leq n < 5.^[18]^[17] For the uniform density case (n = 0), this recovers the familiar factor of $3/5; as n approaches 5, the structure becomes more centrally condensed with an effectively infinite radius, altering the binding efficiency. These models provide realistic approximations for stellar interiors, such as n \approx 1.5 for convective regions in main-sequence stars. This integral approach assumes spherical symmetry, which simplifies calculations for many astrophysical objects but limits applicability to non-spherical distributions, where more complex methods like multipole expansions or full pairwise summations are necessary.^[17]

Extensions and Applications

Relativistic Contexts

In general relativity, the concept of gravitational binding energy extends beyond the Newtonian framework, particularly in regimes of strong gravitational fields where spacetime curvature significantly influences the total energy of a system. The Arnowitt-Deser-Misner (ADM) mass provides a measure of the total energy at spatial infinity for asymptotically flat spacetimes, encompassing both matter contributions and the negative gravitational binding energy. In the weak-field limit, the ADM mass reduces to the Newtonian total mass, which is the integral of the rest-mass density plus the Newtonian gravitational potential energy U, confirming the compatibility of general relativity with Newtonian gravity for weakly bound systems. For compact objects in hydrostatic equilibrium, the Tolman-Oppenheimer-Volkoff (TOV) equation governs the structure, incorporating relativistic effects that modify the binding energy. The mass defect \Delta M = M - \int \rho \, dV, where M is the total gravitational mass and \int \rho \, dV is the integrated proper rest mass, arises from the negative gravitational binding energy, satisfying |\Delta M| c^2 \approx |U| in the relativistic context. This relation highlights how gravitational binding reduces the observed mass compared to the dispersed rest mass of the constituents. A prominent example occurs in neutron stars, where the fractional binding energy |U|/Mc^2 typically ranges from 0.1 to 0.2, reflecting the strong gravitational fields near these compact objects. This substantial binding fraction contributes to observable effects, such as gravitational redshift of emission lines from the surface, which encodes information about the star's compactness and equation of state.^[19] In the case of black holes, such as the Schwarzschild solution, the gravitational binding energy is finite and corresponds to the difference between the ADM mass-energy Mc^2 and the total rest mass-energy of the infalling matter, with the binding energy (approximately 0.4 Mc^2 for non-rotating black holes) radiated away primarily as gravitational waves during collapse.^[20] The first relativistic treatments of gravitational binding in collapsing stars, incorporating these effects, were developed in the 1930s by Oppenheimer and Snyder, who modeled the dust collapse leading to black hole formation.

Astrophysical Role

In self-gravitating astrophysical systems, the virial theorem provides a key connection to gravitational binding energy, stating that for a stable, bound configuration, twice the total kinetic energy K plus the gravitational potential energy \Omega (approximately equal to the binding energy U) sums to zero: $2K + \Omega = 0.^[21] This implies that the total energy E = K + U = U/2 < 0, indicating a bound state where the system's gravitational binding dominates over kinetic support, enabling long-term stability in structures like stars and galaxies.^[22] This relation underscores how gravitational binding energy governs the equilibrium and evolution of such systems against dispersion. During planetary formation, gravitational binding energy plays a crucial role in determining the efficiency of accretion processes, as it represents the energy barrier that incoming material must overcome to integrate into a growing protoplanet. In collisional accretion models, the binding energy of planetesimals and planetary embryos influences the outcomes of impacts, with higher binding energies promoting more efficient merging and growth by dissipating collision energy effectively.^[23] For instance, the Moon's relatively low gravitational binding energy of approximately $1.2 \times 10^{29} J renders it susceptible to tidal disruption risks if its orbit were to decay toward Earth's Roche limit, where tidal forces from the primary body exceed the satellite's self-gravity. In the structure of main-sequence stars, gravitational binding energy sets the scale for the energy budget required to maintain hydrostatic equilibrium, with nuclear fusion output balancing the release of binding energy over the star's lifetime. For the Sun, the magnitude of the binding energy |U| is approximately $6 \times 10^{41} J, corresponding to a main-sequence lifetime of about $10^{10} years as fusion gradually counters the virial-driven contraction.^[24]^[25] This balance ensures that the star's luminosity arises from converting a fraction of its mass into energy, preventing immediate collapse while the binding energy provides the reservoir for evolutionary changes. Gravitational binding energy also defines critical thresholds for collapse in astrophysical media, such as the Jeans mass, where the thermal kinetic energy of a gas cloud falls below the magnitude of its gravitational potential energy, triggering fragmentation and star formation. Clouds exceeding this mass scale become unstable to self-gravity, leading to protostellar core formation in molecular clouds.^[26] Similarly, for white dwarfs, the Chandrasekhar limit of approximately 1.4 solar masses marks the point where the star's gravitational binding energy overwhelms electron degeneracy pressure, causing relativistic instability and potential collapse to a neutron star.^[27]

References

[1]
Lecture 4: Matter & Gravity
Jan 4, 2006 · For example, a sphere of uniform density would have a constant of 3/5. For the Earth, the Gravitational Binding Energy is about 2x1032 Joules, ...
[2]
[PDF] Chapter 2 Stellar Dynamics in Galaxies
The total gravitational potential energy of a uniform sphere of mass M and radius R. (a standard result) is. U = −. 3. 5. GM2. R. 20. Page 7. where G is the ...Missing: binding | Show results with:binding
[3]
[PDF] 10. Timescales in Stellar Interiors
The Virial Theorem relates a system's internal kinetic energy to its gravitational binding energy, and is applied to systems in equilibrium. For stars, it ...
[4]
http://www.bartol.udel.edu/~owocki/phys633/Notes_Part_3.pdf
[5]
Neutrinos and supernova theory
### Summary of Gravitational Binding Energy from the Article
[6]
Energy in Newtonian Gravity | Foundations of Physics
Dec 12, 2022 · We have adopted the convention that the potential \Phi is positive. At some points we may set the constant G to 1. The formula using the Green ...
[7]
Gravitational Binding Energy - an overview | ScienceDirect Topics
Gravitational binding energy (BE) is defined as the energy required to disperse a mass against the gravitational forces holding it together, which increases ...
[8]
Specific properties of gravitational energy - AstroNuclPhysics
It is the potential energy Ep of bodies in the gravitational field and the associated gravitational binding energy, and during the motion of bodies their ...<|control11|><|separator|>
[9]
Gravity-dominated Collisions: A Model for the Largest Remnant ...
Their disruption threshold K* is a multiple of the gravitational binding energy of the two-body system at the point of collision, U G , mutual ; we introduce U ...<|control11|><|separator|>
[10]
Consequences of strong compression in tidal disruption events
... tidal forces do an amount of work exceeding the star's gravitational binding energy, and the star's fluid elements begin moving on roughly geodesic trajectories ...
[11]
Is the whole the sum of its parts? - Einstein-Online
The forces binding protons and neutrons together to form atomic nuclei are considerably stronger, with binding energies that are a few million or even billion ...
[12]
[PDF] The Source of Solar Energy, ca. 1840-1910: From Meteoric ... - arXiv
According to the Helmholtz or Helmholtz-Thomson theory, gravitational contraction was the source of the Sun's thermal energy. The theory predicted a solar ...
[13]
[PDF] 8.901 Lecture Notes - MIT
which is the gravitational potential (or binding energy) of a uniform sphere. All together, this is equivalent to the energy it would take to disassemble this.
[14]
Potential due to uniform sphere
The gravitational potential outside a uniform sphere of mass $M$ is the same as that generated by a point mass $M$ located at the sphere's center.Missing: binding | Show results with:binding
[15]
Facts About Earth - NASA Science
With an equatorial diameter of 7,926 miles (12,756 kilometers), Earth is the biggest of the terrestrial planets and the fifth largest planet in our solar system ...
[16]
An introduction to the study of stellar structure - ADS
An introduction to the study of stellar structure. Chandrasekhar, Subrahmanyan. Abstract. Publication: Chicago. Pub Date: 1939; Bibcode: 1939isss.book.
[17]
[PDF] POLYTROPES
GM2. 2R − n + 1. 6 s. Zc. GMrdMr r. = = −. GM2. 2R. + n + 1. 6. Ω=Ω. (poly.20). The last line of eq. (poly.20) gives the desired answer, i.e. Ω = − 3. 5−n. GM2.
[18]
The Effect of Neutron Star Gravitational Binding Energy on ...
This difference in the rate of change of the masses has a simple relation with the star's gravitational redshift. We show that this effect has the potential ...
[19]
[PDF] Poisson Equation & Virial Theorem Gravity in Astrophysical Fluids
In order for the star to remain in virial equilibrium its kinetic energy, which is proportional to temperature, has to increase; the star's energy loss results.
[20]
[PDF] GRAVITATIONAL POTENTIAL
Gravitational potential is the gravitational energy per unit mass. It is used to find force or potential field of a stellar distribution.
[21]
[PDF] Similar-Sized Collisions and the Diversity of Planets
Oct 11, 2009 · Figure 6: Accretion efficiency for colliding planetary ... The total gravitational binding energy of a uniform planet of mass M and radius R is.
[22]
[PDF] Stellar energy sources: gravitation In the first class we had fun ...
Then the gravitational binding energy is GMdM/R; note that this is positive, because it's the energy we would need to put into the layer dM to get it to ...Missing: definition | Show results with:definition
[23]
The Jeans mass
The gravitational energy $ U$ within the perturbation, is of order $ U\sim G ... So you've learnt how important cooling is for the formation of stars.Missing: binding | Show results with:binding
[24]
[PDF] arXiv:1804.05176v2 [gr-qc] 3 Sep 2018
Sep 3, 2018 · The largest allowed mass for electron degeneracy pressure to balance gravitational collapse is the Chandrasekhar limit. Here we have set G ...<|separator|>