Fact-checked by Grok 2 weeks ago

Gravitational collapse

Gravitational collapse is the process by which a massive astronomical object, such as a star, gas cloud, or protoplanetary disk, contracts under the influence of its own gravity when internal supporting forces like thermal pressure, radiation, or degeneracy pressure fail to counteract gravitational attraction, leading to increased density, temperature, and potentially the formation of compact objects.^[1] In the context of star formation, gravitational collapse initiates when a dense region within a molecular cloud becomes Jeans unstable—exceeding the balance between self-gravity and turbulent or thermal support—causing the cloud fragment to contract, heat up, and eventually form a protostar that accretes material until nuclear fusion begins, halting further collapse. This process typically occurs over timescales of thousands to millions of years and is crucial for the birth of stars across all masses, from low-mass red dwarfs to massive O-type stars.^[2] For evolved massive stars (greater than about 8 solar masses), gravitational collapse marks the endpoint of nuclear burning when the iron core, unable to sustain fusion, loses pressure support and implodes in seconds, triggering a core-collapse supernova that ejects outer layers while the remnant core compresses into a neutron star (if mass is 1.4–3 solar masses) or a black hole (if exceeding roughly 3 solar masses).^[2]^[1] This violent phase releases immense energy, including gravitational waves detectable by instruments like LIGO, and can produce relativistic jets or gamma-ray bursts under certain conditions.^[1]^[3] Theoretically, gravitational collapse in general relativity often culminates in spacetime singularities—regions of infinite curvature—either hidden behind an event horizon (forming black holes, as in the Oppenheimer-Snyder model of dust collapse) or, in some idealized scenarios, as naked singularities if initial conditions allow radial null geodesics to escape. These outcomes underpin the cosmic censorship hypothesis, which posits that generic collapses avoid observable naked singularities to preserve predictability in physics, though ongoing research explores quantum effects and critical phenomena near the black hole threshold.

Fundamentals

Definition and Mechanism

Gravitational collapse refers to the irreversible contraction of a massive body under its own gravitational attraction, occurring when internal pressures—such as thermal, radiation, or degeneracy pressure—fail to provide adequate support against self-gravity.^[1] This process is fundamental to the dynamics of self-gravitating systems, where the inward pull of gravity dominates over outward forces, leading to a runaway compression of the material. The underlying mechanism arises from the imbalance between the system's gravitational potential energy and its internal energy. In equilibrium, self-gravitating systems obey the virial theorem, which states that for a stable configuration, twice the total kinetic energy K plus the gravitational potential energy W (which is negative) equals zero:

$2K + W = 0

^[4] During collapse, as supportive pressures diminish, the release of gravitational potential energy is partially converted to kinetic energy, with approximately half manifesting as thermal energy that heats the collapsing material.^[4] This heating arises directly from the virial relation, as the system's contraction increases kinetic contributions from random motions, driving temperatures upward.^[4] The concept of gravitational instability in gaseous spheres was first analyzed by James Jeans in his 1902 paper 'The Stability of a Spherical Nebula',^[5] and later articulated in the context of star formation in his 1929 work Astronomy and Cosmogony,^[6] where he further developed the analysis, though it has since been extended to describe collapse in any self-gravitating system.

Physical Conditions and Instabilities

Gravitational collapse in self-gravitating systems requires the mass to exceed a critical threshold where gravitational attraction overcomes supporting pressure gradients, leading to net inward acceleration. This condition arises when the system's potential energy dominates kinetic and thermal energies, typically in regions of sufficient density where virial equilibrium cannot be maintained.^[5] The viability of collapse further depends on density \rho, temperature T, and composition, as these parameters determine the sound speed c_s = \sqrt{\gamma k_B T / (\mu m_H)}, where \gamma is the adiabatic index, k_B is Boltzmann's constant, \mu is the mean molecular weight influenced by composition (e.g., higher for ionized gas), and m_H is the hydrogen mass. Higher densities shorten the critical scale for instability, while elevated temperatures or lighter compositions increase c_s, enhancing resistance to collapse.^[5] The key mechanism triggering collapse is the Jeans instability, identified through linear perturbation analysis of a uniform, self-gravitating medium in hydrostatic equilibrium. Small density perturbations grow exponentially if their spatial scale exceeds the Jeans length, as gravitational amplification outpaces pressure restoration waves; shorter scales result in acoustic oscillations rather than runaway collapse. This criterion delineates stable from unstable configurations in dilute media like interstellar gas.^[5] The Jeans length is expressed as

\lambda_J = \sqrt{\frac{\pi c_s^2}{G \rho}},

where G is the gravitational constant; perturbations with wavelength \lambda > \lambda_J are unstable, corresponding to masses above the Jeans mass M_J \approx (4\pi/3) \rho \lambda_J^3. This sets the minimal size for viable collapsing fragments in a given medium.^[5] Following the onset of instability, the dynamical evolution proceeds on the free-fall timescale, the characteristic duration for a pressureless assembly to contract under self-gravity from initial conditions to singularity. For a uniform-density sphere, this is derived by integrating the radial equation of motion, yielding

t_{ff} = \sqrt{\frac{3\pi}{32 G \rho}}.

This timescale provides a benchmark for collapse duration, scaling inversely with the square root of density and marking the point where gravitational contraction becomes rapid and irreversible absent external support. Additional factors modulate these conditions by altering effective support or perturbation growth. Turbulence introduces supersonic motions that can suppress large-scale collapse through enhanced dispersion or, conversely, foster local overdensities that exceed the Jeans threshold and initiate fragmentation.^[7] Magnetic fields contribute stabilizing Lorentz forces via magnetic pressure and tension, raising the effective Jeans length and preventing collapse in weakly ionized plasmas unless ambipolar diffusion allows field decoupling.^[8] Rotation imparts centrifugal support, modifying the effective gravity and potentially halting radial infall or leading to flattened, disk-like structures if angular momentum is conserved during contraction.^[9]

Applications in Astrophysics

Star Formation from Gas Clouds

Gravitational collapse serves as the fundamental mechanism driving the formation of stars from interstellar gas, primarily through the fragmentation of giant molecular clouds (GMCs) into smaller, denser cores. These GMCs, with masses ranging from $10^5 to $10^6 solar masses, are turbulent structures where regions exceeding the critical density for gravitational instability begin to contract. In the global hierarchical collapse (GHC) model, the process starts with slow contraction of mostly atomic gas, accelerating as the cloud cools and becomes molecular, leading to a cascade of fragmentation into filaments, clumps, and eventually prestellar cores. This hierarchical fragmentation ensures efficient distribution of mass, with turbulence aiding in creating overdensities that initiate collapse while preventing premature dispersal.^[10] The key stages of this collapse process unfold as follows: initial cloud contraction occurs when the gravitational potential overcomes thermal and turbulent support, compressing the gas to densities around $10^{-20} g/cm³. As angular momentum is conserved, the collapsing material flattens into a rotationally supported protostellar disk, typically on scales of 100 AU, where magnetic fields and turbulence regulate the dynamics. Accretion then proceeds onto the central protostar, with material spiraling inward through the disk at rates of $10^{-5} to $10^{-4} solar masses per year, building the star's mass over the embedded phase. These stages are interconnected, with the disk serving as a reservoir that modulates the episodic nature of accretion.^[11] Central to core formation is the Jeans mass, which sets the characteristic scale for gravitational instability in the collapsing gas. Defined as M_J = \frac{4\pi}{3} \rho \left( \frac{\lambda_J}{2} \right)^3, where \rho is the density and \lambda_J is the Jeans length, this mass decreases with increasing density, allowing smaller fragments to form as collapse proceeds. In molecular clouds at temperatures of 10 K, the Jeans mass for typical densities yields cores of 0.1 to 10 solar masses, aligning with observed prestellar core masses. Observational evidence for this collapse comes from infrared surveys of Bok globules—compact, isolated dark clouds that are ideal laboratories for low-mass star formation. Near-infrared imaging reveals infalling envelopes and outflow signatures in globules like B335, indicating dynamic collapse with velocities consistent with free-fall. These observations, often using telescopes like Spitzer, show typical core masses of 1-5 solar masses and collapse timescales of $10^4 to $10^6 years, matching theoretical free-fall times for dense cores. More recent observations with the James Webb Space Telescope (JWST), as of 2025, have provided higher-resolution views of infalling protostellar envelopes and the influence of strong magnetic fields on collapse dynamics.^[12] A key theoretical description of the early collapse phase is provided by the Larson-Penston solution, which models the self-similar infall of an isothermal sphere toward a central singularity. This solution features a uniform expansion wave propagating outward at the sound speed, with supersonic infall velocities approaching four times the sound speed near the center, offering a dynamic counterpart to the static singular isothermal sphere. It has been influential in interpreting the initial conditions for protostellar collapse in numerical simulations.

Core Collapse in Stellar Evolution

In massive stars with initial masses exceeding 8 solar masses (M⊙), gravitational core collapse occurs during the final stages of stellar evolution, after the exhaustion of hydrogen and helium fuels in the core through successive nuclear burning phases. These stars progress through carbon, oxygen, neon, and silicon burning, ultimately forming an iron-nickel core that cannot sustain further fusion, as iron-group elements absorb rather than release energy upon fusion. The iron core grows via silicon shell burning at temperatures around 3–4 × 10^9 K over timescales of days, reaching a mass of approximately 1.2–1.6 M⊙, supported initially by electron degeneracy pressure. The collapse is triggered when the iron core mass surpasses the effective Chandrasekhar limit, the maximum mass sustainable by relativistic electron degeneracy pressure against gravity, leading to the failure of this pressure and initiating rapid infall. The Chandrasekhar mass is given by

M_\mathrm{Ch} = \left( \frac{\hbar c}{G} \right)^{3/2} \frac{1}{(\mu_e m_p)^2} \times K \approx 1.44 \, M_\odot,

where \hbar is the reduced Planck constant, c is the speed of light, G is the gravitational constant, \mu_e \approx 2 is the mean molecular weight per electron, m_p is the proton mass, and K \approx 3.1 is a structural constant from polytropic models. In the stellar core context, the effective limit is lower, around 1.1–1.3 M⊙, due to neutronization processes—electron captures on protons and nuclei that reduce the electron fraction Y_e from ~0.46 to ~0.30–0.40 during infall, increasing \mu_e and thus decreasing the support. This instability causes the core to contract homologously at speeds reaching 0.3–0.4c, compressing to nuclear densities (~2.8 × 10^{14} g cm^{-3}) in less than 1 second, primarily milliseconds for the inner core. The collapse proceeds in distinct stages: initial homologous infall of the inner core (r < 100–200 km), followed by a dynamical phase where the core bounces upon reaching nuclear density, halted by neutron degeneracy pressure and the equation of state of nuclear matter. This bounce, occurring at ~0.6–0.8 M⊙ core mass, generates an outgoing shock wave that dissociates infalling material, initially stalling but potentially reviving through neutrino heating to drive the explosion. The entire collapse-to-bounce sequence unfolds on dynamical timescales of ~10–100 milliseconds. This process culminates in Type II core-collapse supernovae, explosive events that eject the star's envelope at velocities of ~10,000 km s^{-1}, releasing ~10^{51} erg of kinetic energy. The remnant is typically a neutron star, though details of its formation are addressed elsewhere.^[13]

Theoretical Frameworks

Newtonian Models

Newtonian models of gravitational collapse describe the dynamics of self-gravitating systems using classical mechanics and Poisson's equation for the gravitational potential, assuming spherical symmetry and neglecting relativistic effects. These models typically consider a fluid or dust distribution where the gravitational attraction overcomes internal pressures, leading to contraction. A key assumption is the use of Lagrangian or Eulerian coordinates to track the motion of mass elements, with the gravitational field determined by the enclosed mass within a given radius.^[14] Homologous collapse refers to scenarios where all parts of the collapsing body contract uniformly, maintaining proportionality in their radial positions relative to the center. This occurs in uniform density models, such as pressureless dust spheres, where the entire structure scales self-similarly during infall. In the Newtonian limit, this approximates the dust collapse dynamics akin to the relativistic Oppenheimer-Snyder model, but without spacetime curvature; the collapse proceeds homologously until a central singularity forms in finite time. Seminal work by Goldreich and Weber demonstrated that such homologous solutions exist for polytropic gases with adiabatic index γ = 4/3, representing marginally stable stellar cores where slight pressure reductions trigger uniform contraction.^[15] Self-similar solutions extend these models by assuming the system's evolution depends only on a scaled variable, such as ξ = r / (a t), where r is the radial coordinate, t is time, and a is a constant. For pressureless collapse, Hunter-type solutions describe families of smooth, spherically symmetric inflows that develop a central singularity while preserving self-similarity. These solutions, numerically identified by Hunter, reveal discrete branches of collapse profiles that transition from expansion to contraction, highlighting the inevitability of singularity formation under pure gravitational dominance. Unlike static equilibria, they capture dynamic phases where density profiles steepen toward the center. The core dynamics in these models are governed by the equation of motion for a fluid element at radius r, incorporating the gravitational acceleration from the enclosed mass M(r):

\frac{d^2 r}{dt^2} = -\frac{G M(r)}{r^2}

This equation, derived from Newton's law of gravitation applied to spherical shells, predicts a finite-time singularity as r → 0, where acceleration diverges. Integration shows that initial uniform density distributions collapse to a point in time t_coll ≈ (3π / (32 G ρ_0))^{1/2}, with ρ_0 the initial density, establishing the scale for collapse duration. A specific application is the pressureless (cold dust) approximation, ideal for modeling primordial collapse in the early universe where thermal pressures are negligible compared to gravity. In this limit, matter behaves as collisionless particles following geodesics in the Newtonian potential, mimicking the Friedmann-Lemaître-Robertson-Walker dust-dominated cosmology on small scales. Overdense regions, triggered by instabilities like Jeans modes, contract homologously, forming compact structures that seed larger cosmic filaments without radiation or magnetic support.^[16] These Newtonian models have limitations, as they ignore pressure gradients that could halt collapse and relativistic effects that become crucial near singularities. They are valid primarily for initial phases of star formation or cosmological structure growth, where velocities remain sub-relativistic and densities are not extreme.^[17]

General Relativistic Models

In general relativity, gravitational collapse is described by solutions to the Einstein field equations that incorporate the curvature of spacetime, leading to phenomena such as event horizons and singularities absent in Newtonian treatments. These models reveal how sufficiently massive configurations inevitably undergo dynamical collapse, with spacetime geometry playing a crucial role in determining the evolution and endpoint.^[18] A foundational example is the Oppenheimer-Snyder model, which describes the spherically symmetric collapse of an incompressible, pressureless dust cloud in general relativity. In this solution, a homogeneous sphere initially at rest contracts homologously, matching interior Friedmann-Lemaître-Robertson-Walker geometry to the exterior Schwarzschild vacuum solution. As the collapse proceeds, an event horizon forms at the Schwarzschild radius, enclosing the collapsing matter and demonstrating the birth of a black hole from realistic initial conditions without pressure support.^[18] For configurations with pressure, the Tolman-Oppenheimer-Volkoff (TOV) equation governs hydrostatic equilibrium in spherically symmetric stars under general relativity, generalizing the Newtonian Lane-Emden equation by including relativistic corrections. The equation is given by

\frac{dP}{dr} = -\frac{G M(r) \rho}{r^2} \left(1 + \frac{P}{\rho c^2}\right)\left(1 + \frac{4\pi r^3 P}{M(r) c^2}\right)\left(1 - \frac{2 G M(r)}{r c^2}\right)^{-1},

where P is pressure, \rho is energy density, M(r) is the enclosed mass, G is the gravitational constant, and c is the speed of light. At high densities, the relativistic factors amplify gravitational forces, causing the equation to predict instability and collapse when central pressure diverges, as no stable equilibrium exists beyond a maximum mass.^[19] Singularity theorems, developed by Penrose and Hawking, rigorously prove that gravitational collapse leads to spacetime singularities under physically reasonable conditions. Penrose's 1965 theorem shows that if a trapped surface forms during collapse—as occurs when matter crosses its own event horizon—then geodesics are incomplete, implying a singularity where curvature becomes infinite. Hawking extended this in collaboration with Penrose, demonstrating that singularities are inevitable in collapsing spacetimes satisfying the null convergence condition (energy density non-negative) and possessing a Cauchy hypersurface with non-compact regular boundary, applicable to realistic stellar cores.^[20] Numerical studies of critical phenomena near the threshold of black hole formation reveal universal behavior in type II collapse, as discovered by Choptuik. In simulations of scalar field collapse, varying the initial amplitude p around a critical value p^* yields a black hole mass scaling as M \sim (p - p^*)^\gamma, where the exponent \gamma \approx 0.37 is universal across families of initial data, independent of the matter model, indicating self-similar critical dynamics echoing periodically.^[21] While naked singularities—unclothed by event horizons—are possible in certain idealized models violating energy conditions, the cosmic censorship hypothesis posits that generic gravitational collapse in asymptotically flat spacetimes produces hidden singularities, preserving predictability; observational cosmology and critical collapse simulations support this by favoring horizon formation over naked exposure.^[22]

Outcomes and Remnants

Neutron Stars

Neutron stars form as the remnants of gravitational collapse in the cores of massive stars with initial masses between approximately 8 and 20 solar masses (M⊙), where the iron core, typically weighing 1.4 to 3 M⊙, implodes during the final stages of stellar evolution, expelling the outer envelope in a core-collapse supernova.^[23] This process occurs when electron degeneracy pressure can no longer support the core against gravity, leading to the rapid compression of protons and electrons into neutrons, forming a proto-neutron star that cools and contracts over seconds to minutes.^[23] The structure of a neutron star is maintained by neutron degeneracy pressure, arising from the Pauli exclusion principle applied to a dense Fermi gas of neutrons, which provides the primary support against further gravitational collapse.^[24] The equation of state (EOS) governing this pressure is derived from nuclear physics, describing the relationship between pressure and density in the extreme conditions inside the star, where matter reaches nuclear densities and may include exotic phases like hyperons or quark matter, though the precise EOS remains uncertain and is constrained by observations.^[25] The Tolman-Oppenheimer-Volkoff (TOV) equation, the relativistic generalization of hydrostatic equilibrium, determines the maximum stable mass for a neutron star, known as the Oppenheimer-Volkoff limit, estimated at approximately 2 to 3 M⊙ depending on the EOS; beyond this mass, the star collapses further.^[26] Typical neutron stars have masses around 1.4 M⊙, radii of 10 to 15 kilometers, and average densities on the order of 10¹⁷ kg/m³, comparable to the density of atomic nuclei, making them among the densest objects in the universe short of black holes.^[27] Many neutron stars rotate rapidly, with periods ranging from milliseconds to seconds, manifesting as pulsars when their strong magnetic fields (∼10¹² gauss) sweep beams of radiation across our line of sight, producing observable pulses.^[28] The first neutron stars were discovered indirectly through pulsars in 1967 by Jocelyn Bell Burnell and Antony Hewish using a radio telescope at the Mullard Radio Astronomy Observatory.^[29] In the 2020s, the Neutron Star Interior Composition Explorer (NICER) mission has refined mass-radius relations through X-ray observations of rotation-powered pulsars, providing tighter constraints on the EOS and confirming radii around 12-13 km for ∼1.4 M⊙ stars.^[30]

Black Holes

Black holes represent the endpoint of complete gravitational collapse when the core mass exceeds approximately 3 solar masses, surpassing the Tolman-Oppenheimer-Volkoff limit that supports neutron stars against further implosion. In such cases, the core of a massive star, typically from progenitors exceeding 20 solar masses, undergoes unrelenting collapse following the exhaustion of nuclear fuel and supernova ejection of outer layers.^[31] Alternatively, direct collapse scenarios can form black hole seeds in the early universe through the implosion of pristine, massive gas clouds under specific conditions of low metallicity and intense radiation, bypassing star formation altogether. The boundary of a non-rotating black hole is defined by the event horizon, a surface beyond which escape is impossible, located at the Schwarzschild radius given by

R_s = \frac{2 G M}{c^2},

where M is the mass, G is the gravitational constant, and c is the speed of light.^[32] General relativistic models of collapse, such as those extending the Oppenheimer-Snyder dust solution, predict the formation of this horizon as the collapsing matter crosses R_s. For stable stars, post-Newtonian stability analysis imposes a theoretical minimum radius of about \frac{9}{8} R_s; configurations denser than this inevitably collapse due to the Buchdahl bound on spherically symmetric fluid spheres. Black holes are classified by mass into stellar-mass (roughly 3 to 100 solar masses, formed from stellar cores), intermediate-mass (hundreds to tens of thousands of solar masses, possibly from cluster mergers or direct collapse), and supermassive (millions to billions of solar masses, residing in galactic centers).^[31] According to the no-hair theorem, isolated black holes are fully characterized by only three parameters—mass, electric charge (typically negligible), and angular momentum—with no other distinguishing "hair" such as multipole moments from their progenitor structure. Observational confirmation includes the first direct image of a supermassive black hole's shadow in the galaxy M87, captured by the Event Horizon Telescope in 2019, revealing a dark central region consistent with the predicted event horizon size.^[33] By 2025, the LIGO-Virgo-KAGRA collaboration has detected over 300 gravitational wave signals from black hole mergers, including record-breaking events like the most massive binary merger observed in July 2025, affirming the existence and dynamics of these collapsed objects.^[34]^[35]

Other Compact Objects

White dwarfs represent the endpoint of gravitational collapse for low- and intermediate-mass stars, typically those with initial masses between 0.08 and 8 solar masses, where further contraction is resisted by electron degeneracy pressure rather than ongoing nuclear fusion. These stars form when the progenitor exhausts its nuclear fuel, leading to the ejection of outer layers in a planetary nebula, leaving behind a hot, dense core composed primarily of carbon and oxygen.^[36] This process prevents total collapse into more extreme remnants for masses below the critical threshold, maintaining hydrostatic equilibrium through quantum mechanical effects that limit electron velocities to near the speed of light.^[37] The stability of white dwarfs is bounded by the Chandrasekhar limit, the maximum mass at which electron degeneracy pressure can support the star against gravity, calculated as M_{\rm Ch} \approx 1.44 \, M_\odot.^[37] Above this limit, the relativistic nature of the degenerate electrons causes the pressure to increase insufficiently with density, resulting in the ignition of thermonuclear fusion and a Type Ia supernova explosion.^[36] Typical white dwarfs possess masses around 0.6 M_\odot, with radii comparable to Earth's—approximately 10,000 km—and central densities on the order of $10^9 kg/m³, making them among the densest known objects that avoid forming event horizons. Over timescales of billions of years, these remnants cool by radiating residual thermal energy, transitioning from luminous white to faint black dwarfs without further structural change.^[36] The first white dwarf, Sirius B, was predicted in 1844 by Friedrich Bessel through astrometric observations of perturbations in Sirius A's orbit and later directly observed in 1862, providing early evidence for such compact objects. Modern surveys, including data from the Gaia mission in the 2010s and 2020s, have cataloged hundreds of thousands of white dwarfs, deriving mass distributions that peak near 0.6 M_\odot and confirm the prevalence of these remnants in the solar neighborhood.^[38]^[39] Beyond established remnants, exotic compact objects like quark stars have been theoretically proposed as hypothetical alternatives where extreme densities convert neutrons into a degenerate soup of up, down, and strange quarks, potentially stable beyond neutron star limits but lacking observational confirmation. Similarly, preon stars, envisioned as aggregates of hypothetical preons—the substructure of quarks and leptons—could resist collapse at even higher densities, offering a speculative bridge to understanding fundamental particles, though their existence remains unverified.^[40]

Observational Evidence

Supernovae and Gamma-Ray Bursts

Core-collapse supernovae of Types II, Ib, and Ic arise from the gravitational collapse of massive stars with initial masses exceeding approximately 8 solar masses, where the iron core exceeds the Chandrasekhar limit and implodes, triggering an explosive ejection of the star's outer layers.^[41] These events release kinetic energy on the order of 10^51 ergs (approximately 10^44 joules) in the form of expanding ejecta, with the explosion preceded by a intense burst of neutrinos carrying away about 99% of the gravitational binding energy released during core collapse.^[42] Type II supernovae exhibit hydrogen lines in their spectra due to the retention of the progenitor's hydrogen envelope, while Types Ib and Ic lack prominent hydrogen (Ib) or both hydrogen and helium (Ic), reflecting stripped-envelope progenitors such as Wolf-Rayet stars.^[43] The explosion mechanism involves the initial implosion forming a proto-neutron star, followed by a rebound that generates an outward-propagating shock wave; however, this shock stalls near 150-200 km from the center due to energy losses from photodissociation and neutrino emission.^[13] Revival occurs through the delayed neutrino heating mechanism, where neutrinos emitted from the hot proto-neutron star deposit energy in the gain region behind the stalled shock, heating the material and driving convection that reinvigorates the shock to propagate outward and unbind the stellar envelope.^[44] This process typically unfolds over milliseconds to seconds post-bounce, with simulations confirming that sufficient heating—around 10% of the neutrino luminosity—can power the explosion.^[41] Observational confirmation of this neutrino-driven process came from Supernova 1987A (SN 1987A), the collapse of a massive star in the Large Magellanic Cloud, where detectors including Kamiokande II and the Irvine-Michigan-Brookhaven (IMB) experiment recorded a burst of 11 and 8 electron antineutrino events, respectively, about 3 hours before the optical detection on February 23, 1987, providing the first direct evidence of neutrinos from a core-collapse event.^[45]^[46] The total energy released in neutrinos was approximately 3 × 10^53 ergs, consistent with the gravitational binding energy of a forming neutron star, estimated as roughly GM^2/R ≈ 3 × 10^53 ergs for a 1.4 solar mass core with radius ~10 km.^[47] Gamma-ray bursts (GRBs), particularly long-duration GRBs (lasting >2 seconds), are linked to core collapse in rapidly rotating massive stars via the collapsar model, where the collapsing core forms a black hole surrounded by a centrifugally supported accretion disk, and relativistic jets powered by neutrino annihilation or magnetohydrodynamic processes break out to produce the burst.^[48] In this scenario, proposed by Woosley in 1993 and refined through simulations, the rapid rotation prevents fallback of the core and enables jet formation, often in low-metallicity environments that preserve angular momentum.^[49] A prominent example is GRB 080319B, detected on March 19, 2008, at redshift z=0.937, which was the intrinsically brightest GRB observed, with its optical afterglow visible to the naked eye from Earth despite being 7.5 billion light-years away, highlighting the extreme luminosity of collapsar-driven events.^[50]

Gravitational Waves from Mergers

Gravitational collapse in massive stars often results in the formation of neutron stars or black holes, which can later pair up in binary systems and merge, emitting detectable gravitational waves. These mergers include binary neutron star (NS-NS), binary black hole (BH-BH), and neutron star-black hole (NS-BH) systems, where the compact objects are remnants of prior collapse events. The gravitational radiation from such mergers provides a direct probe of the end products of stellar evolution, revealing properties like masses and spins that trace back to the collapse process. The gravitational wave signal from a binary merger evolves through distinct phases: the inspiral, where the objects spiral inward due to energy loss via gravitational radiation; the merger, a brief violent coalescence; and the ringdown, where the final remnant settles into a stable configuration, emitting waves at characteristic quasi-normal mode frequencies dictated by general relativity. During the inspiral, the waveform exhibits a "chirp" signature, with increasing frequency and amplitude as the orbital separation decreases. This phase is well-approximated by the leading-order post-Newtonian formula for the characteristic strain amplitude:

h \approx \frac{4}{D} \left( \frac{G \mathcal{M}_c}{c^2} \right)^{5/3} \left( \pi f \frac{G}{c^3} \right)^{2/3} \frac{1}{c},

where \mathcal{M}_c is the chirp mass \mathcal{M}_c = \frac{(m_1 m_2)^{3/5}}{(m_1 + m_2)^{1/5}}, f is the gravitational wave frequency (twice the orbital frequency), D is the luminosity distance to the source, and G and c are the gravitational constant and speed of light, respectively. The ringdown phase features damped oscillations at quasi-normal mode frequencies, which for non-spinning black holes are \omega_{lmn} \approx 0.5 (1 - 0.63 (1 - a)^{3/10}) / M_f for the fundamental (l=2, m=2) mode, with a the final spin and M_f the final mass, allowing tests of general relativity in the strong-field regime. Since the first direct detection in 2015 by the Advanced LIGO observatories, the LIGO-Virgo-KAGRA collaboration has identified gravitational waves from numerous compact binary mergers. The inaugural event, GW150914, was a BH-BH merger at a distance of about 410 Mpc, involving black holes of approximately 36 and 29 solar masses, confirming the existence of stellar-mass black holes formed via gravitational collapse. A landmark NS-NS merger, GW170817, detected in 2017 at 40 Mpc, produced a chirp signal followed by electromagnetic counterparts including gamma rays and a kilonova, establishing multimessenger astronomy and constraining the equation of state of neutron star matter. As of November 19, 2025, following the conclusion of the fourth observing run (O4, 2023–November 18, 2025), approximately 350 gravitational wave events have been cataloged, with over 280 being BH-BH mergers, several NS-BH mergers (such as GW200105 and GW200115), and a few confirmed NS-NS events (including GW170817). These detections, detailed in GWTC-4 and subsequent releases, inform the population statistics of collapse remnants, revealing a range of black hole masses from about 5 to over 100 solar masses and suggesting formation channels like isolated binary evolution or cluster dynamics.^[51]

References

[1]
Gravitational Collapse - an overview | ScienceDirect Topics
Gravitational collapse refers to the process where a massive object, under the influence of its own gravity, contracts and can lead to the formation of ...<|control11|><|separator|>
[2]
The Life Cycles of Stars - Imagine the Universe! - NASA
May 30, 2025 · In less than a second, the star begins the final phase of its gravitational collapse. The core temperature rises to over 100 billion degrees as ...
[3]
Gravitational Contraction and Star Formation - HyperPhysics
The source of the energy for star formation is gravitational collapse - this collapse must provide enough energy to heat the gas of the protostar to the ...Missing: definition | Show results with:definition
[4]
Astronomy And Cosmogony : The University Press, Cambridge
Jan 21, 2017 · Astronomy And Cosmogony. by: The University Press, Cambridge ... PDF WITH TEXT download · download 1 file · SINGLE PAGE PROCESSED JP2 ZIP ...
[5]
I. The stability of a spherical nebula | Philosophical Transactions of ...
(2022) Kinetic theory of Jeans' gravitational instability in millicharged dark matter system, Chinese Physics B, 10.1088/1674-1056/ac5239, 31:7, (070401), ...
[6]
The influence of turbulence during magnetized core collapse and its ...
We study the collapse of turbulent, magnetized prestellar cores, focusing on the effects of magnetic diffusion, and misalignment between rotation axis and ...
[7]
impact of turbulence and magnetic field orientation on star-forming ...
Only magnetic fields parallel to the major axis can stabilize the filament against radial collapse and reduce fragmentation. We will also show that turbulence ...
[8]
[PDF] gravitational waves from axisymmetric, rotating stellar core collapse
Jul 26, 2003 · When rotation is included, centrifugal forces counteracting gravity's pull change the dynamics of collapse. Depending on the total amount of ...
[9]
[PDF] Formation and Evolution of the Protoplanetary Disk
The disk begins to form at the growth stage of the stellar-like core and expands during accretion to the present dimensions of the solar system. Accretion at ...
[10]
The mechanism(s) of core-collapse supernovae - Journals
Sep 18, 2017 · Core-collapse supernovae (CCSNe) are the explosions that attend the deaths of massive stars. Despite decades of research, several aspects of the mechanism that ...Introduction · Modern view of stellar core... · Neutrino-driven explosions in...
[11]
Newtonian spherical gravitational collapse - IOPscience
The spherical collapse of gas clouds is reduced to quadrature by specifying the mass function n(r, t). Exact solutions for both shearing and shear-free flow are ...Missing: homologous | Show results with:homologous
[12]
Cold collapse and bounce of an FLRW cloud - Oxford Academic
Proceeding with the Newtonian code and abandoning the approximation of zero pressure, we simulate the collapse of FLRW clouds of masses 5, 20, and 1000 M ⊙ with ...
[13]
Gravitational Waves from Gravitational Collapse - PMC
In rotating stars, gravitational radiation reaction drives the r-modes toward unstable growth [5, 77]. In hot, rapidly rotating neutron stars, this instability ...
[14]
On Continued Gravitational Contraction | Phys. Rev.
Abstract. When all thermonuclear sources of energy are exhausted a sufficiently heavy star will collapse. Unless fission due to rotation, the radiation of mass, ...
[15]
On Massive Neutron Cores | Phys. Rev.
In this paper we study the gravitational equilibrium of masses of neutrons, using the equation of state for a cold Fermi gas, and general relativity.
[16]
Gravitational Collapse and Space-Time Singularities | Phys. Rev. Lett.
Oct 6, 2020 · Three scientists were recognized for proving that gravitational collapse can lead to a black hole and for observing the supermassive black hole at the center ...
[17]
Universality and scaling in gravitational collapse of a massless ...
Jan 4, 1993 · I summarize results from a numerical study of spherically symmetric collapse of a massless scalar field.
[18]
[PDF] Gravitational Collapse: the Role of General Relativity - UMD Physics
a « cosmic censor » who forbids the appearance of naked singularities, clothing ... RIVISTA DEL NUOVO CIMENTO. VOLUME I. [32] I. ROBINSON and A. TRAUTMAN: Proc ...
[19]
[1806.07267] Neutron stars formation and Core Collapse Supernovae
Jun 18, 2018 · In this chapter, we review the current status of SNe observations and theoretical modelling, the connection with their progenitor stars, and the properties of ...
[20]
[PDF] arXiv:astro-ph/0506417v1 17 Jun 2005
Jun 17, 2005 · Similarly, a neutron star is stabilised by the degeneracy pressure also. The difference is that for neutron stars the degeneracy pressure ...
[21]
EQUATION OF STATE AND NEUTRON STAR PROPERTIES ...
We discuss in detail the allowed equations of state and the impact of our results on the structure of neutron stars, the crust-core transition density, and the ...
[22]
[1210.5667] Equations of State and Maximum Mass of Neutron Stars ...
Oct 20, 2012 · By combining an equation of state (EOS) with the Tolmann-Oppenheimer-Volkoff ... maximum mass of a neutron star for a given equation of state.
[23]
[PDF] Neutron Star Structure and Measurements - N3AS
Jul 20, 2024 · Neutron stars have high density, pressure, and gravity, are mostly made of neutrons, have radii around 12km, and are formed from massive stars. ...
[24]
Chapter 6 Pulsars
These numbers motivate the definition of the canonical neutron star as being a uniform density sphere with mass M≈1.4M⊙, radius R≈10 km, and moment of inertia I ...
[25]
Observation of a Rapidly Pulsating Radio Source - Nature
Hewish et al. associated them with unusually stable oscillations in compact stars. They are now understood to be rapidly rotating, magnetized neutron stars, or ...
[26]
[PDF] Constraints on the Dense Matter Equation of State and Neutron Star ...
Sep 8, 2021 · The study uses NICER's mass-radius of PSR J0740+6620 to constrain the dense matter equation of state, using two models and providing tight ...
[27]
Black Hole Types - NASA Science
Oct 22, 2024 · If it had around 20 times the Sun's mass or more, the star's core collapses into a stellar-mass black hole. Stellar mass simulation. This ...
[28]
[physics/9905030] On the gravitational field of a mass point ... - arXiv
May 12, 1999 · Plain TeX, 7 pages, English translation of the original paper by K. Schwarzschild. Subjects: History and Philosophy of Physics (physics.hist ...
[29]
Press Release (April 10, 2019): Astronomers Capture First Image of ...
Apr 10, 2019 · The image reveals the black hole at the center of Messier 87 [1], a massive galaxy in the nearby Virgo galaxy cluster. This black hole resides ...
[30]
LIGO Detects Most Massive Black Hole Merger to Date - Caltech
Jul 14, 2025 · The LIGO-Virgo-KAGRA (LVK) Collaboration has detected the merger of the most massive black holes ever observed with gravitational waves.
[31]
Ten years later, LIGO is a black-hole hunting machine | MIT News
Sep 10, 2025 · Today, the collaborating LIGO, Virgo, and Kagra detectors have captured about 300 black hole mergers, including a new event that provides ...
[32]
White Dwarf Stars - Imagine the Universe! - NASA
Near the end of its nuclear burning stage, this type of star expels most of its outer material, creating a planetary nebula. Only the hot core of the star ...
[33]
https://eventhorizontelescope.org/press-release-april-10-2019-astronomers-capture-first-image-black-hole
[34]
Friedrich Bessel: Discoverer of White Dwarf Sirius B | AMNH
Discovering Sirius B long before technology allowed us to see it, or quantum mechanics explained the nature of white dwarfs.
[35]
Gaia white dwarf discoveries - ESA Science & Technology
Prior to the second Gaia data release, made public in 2018, only about 30 000 white dwarfs had been discovered. Using date from Gaia, 486 641 white dwarf ...Missing: mass distribution 2020s
[36]
Preon stars: a new class of cosmic compact objects - ResearchGate
Jan 13, 2025 · We show that a degenerate gas of interacting fermionic preons does allow for stable compact stars, with densities far beyond that in neutron stars and quark ...
[37]
[PDF] the physics of core-collapse supernovae - arXiv
It may well be that the range of masses that actually explodes by forming a neutron star is limited (23). The probable existence of stellar mass black holes in ...
[38]
[PDF] Neutron Stars - SLAC National Accelerator Laboratory
pared to the binding energy of the massive star that forms its progenitor, the binding energy of a neutron star is huge, about GM2/R ∼ 3 × 1053 ergs for M = 1.4 ...
[39]
[1501.04645] Investigations of supernovae and supernova remnants ...
Jan 19, 2015 · Two main physical mechanisms are used to explain supernova explosions: thermonuclear explosion of a white dwarf(Type Ia) and core collapse of a ...
[40]
Conditions for shock revival by neutrino heating in core-collapse ...
Energy deposition by neutrinos can rejuvenate the stalled bounce shock and can provide the energy for the supernova explosion of a massive star.
[41]
Observation of a neutrino burst from the supernova SN1987A
Apr 6, 1987 · A neutrino burst was observed in the Kamiokande II detector on 23 February, 7:35:35 UT (±1 min) during a time interval of 13 sec.
[42]
Neutrinos from SN1987a in the IMB detector - ScienceDirect.com
Abstract. A burst of eight neutrino events preceding the optical detection of the supernova in the Large Magellanic Cloud has been observed by the IMB detector.
[43]
[PDF] arXiv:astro-ph/0612440v1 15 Dec 2006
Dec 15, 2006 · Nearly all (∼ 99%) of the neutron star's binding energy of about 300 bethes is eventually radiated away in neutrinos of all flavors . The.
[44]
LONG GAMMA-RAY TRANSIENTS FROM COLLAPSARS
The collapsar model (Woosley 1993) was originally proposed as an explanation for common, “long-soft” gamma-ray bursts (GRBs) with durations near 20 s. This ...INTRODUCTION · THE PRODUCTION OF... · EXAMPLES OF TYPE 3...
[45]
Gamma-Ray Bursts and Explosions in "Failed Supernovae" - arXiv
Oct 19, 1998 · Collapsars - Gamma-Ray Bursts and Explosions in "Failed Supernovae". Authors:A. MacFadyen, S. E. Woosley.
[46]
Broadband observations of the naked-eye γ-ray burst GRB 080319B
Sep 11, 2008 · 7), GRB 080319B was the brightest optical burst ever observed. An observer in a dark location could have seen the prompt optical emission with ...
[47]
GWTC-4.0: Updated Gravitational-Wave Catalog Released | LIGO Lab
Aug 26, 2025 · This release marks the first in a sequence of planned public data releases from the LVK's ongoing fourth observing run (O4), currently scheduled ...