Mass functions as evidence for hierarchical black hole mergers ✛ ✎ ✦
Context
How stellar-mass black holes grow after birth is one of astrophysics' central questions. According to Einstein's general relativity, accelerating masses emit gravitational waves, and Schwarzschild first found the exact solution for a non-rotating black hole. Observations with LIGO, Virgo, and KAGRA detectors (created with key contributions from Rainer Weiss) revealed a population of merging black holes with large masses and high spins. Two competing hypotheses for their origin exist: growth through repeated black hole mergers in dense star clusters (hierarchical mergers) or gas accretion onto isolated holes. Distinguishing them means understanding stellar populations and cluster dynamics, and tracing the chain from stellar evolution to the formation of massive compact objects.
Methods
Researchers applied hierarchical Bayesian analysis to 259 black hole merger events from the GWTC-5 gravitational-wave catalog. They built a flexible mixture model of two subpopulations describing mass, spin, and redshift distributions within the framework of Einstein's equations, which predict signal characteristics. Mass functions were modeled with splines, accounting for selection effects and detector noise. Separate distributions were recovered for the first subpopulation (low spins, interpreted as products of stellar evolution and core collapse) and the second (high spins — hierarchical remnant candidates). Gravitational wave emission during merger is described by the quadrupole formula, allowing source parameter extraction.
Results
The main result: the mass function shape of the second subpopulation is nearly identical to the final mass distribution after black hole mergers of the first. This similarity is shown in the figure: up to 80 solar masses, the curves practically coincide. The Bhattacharyya coefficient, measuring overlap of normalized distributions, reaches 0.95 in the [40,80] M⊙ interval, indicating a statistically significant match. Mass peaks of the first subpopulation at ~10 M⊙ and ~35 M⊙ reappear in the second at ~17 M⊙ and ~58 M⊙, consistent with expected mass buildup from mergers. Additionally, a sharp cutoff of the low-spin population at 54.2+7.7-7.2 M⊙ is interpreted as the lower edge of the mass gap caused by pair-instability supernovae. This allowed computation of the S-factor for the 12C(α,γ)16O reaction, crucial for nucleosynthesis in massive stars: S300 = 151+30-26 keV·barn, in excellent agreement with theoretical estimates.
Implications
The discovered morphological link is direct proof that high-spin massive black holes are products of earlier black hole mergers of the first subpopulation. An alternative scenario, such as growth via accretion from an accretion disk, would require fine-tuning to reproduce the same spectral features. Thus, hierarchical assembly dominates the formation of the observed binary black hole population, and primordial black holes are not needed to explain it. The derived upper limit on the maximum mass of low-spin holes also challenges single stellar evolution models predicting slowly rotating black holes up to ~93 M⊙.
Future development
As data accumulate from new observing runs of LIGO–Virgo–KAGRA and future detectors like Einstein Telescope and Cosmic Explorer, mass functions will be measured with even greater precision. This will not only strengthen the evidence for hierarchical black hole mergers but also allow detailed study of contributions from higher generations (third and beyond) and link them to the birth of intermediate-mass black holes. Modeling of cluster dynamics and accretion flows will also gain new observational constraints.
Impact
The results will impact compact object astrophysics, stellar evolution theory, nuclear astrophysics (via refined nucleosynthesis), and cosmology, refining the black hole merger rate and star formation history.
Next steps
Next steps include incorporating time delays between merger generations and analyzing contributions from different dynamical environments (globular clusters, active galactic nucleus disks) via detailed numerical simulations and further gravitational wave observations.
Key open problems
The work directly addresses open problems: the origin of the observed black hole mass function, the nature of the mass gap from pair-instability supernovae, the rate of the 12C(α,γ)16O nuclear reaction, and the formation mechanisms of binaries in dense environments.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
The authors applied hierarchical Bayesian analysis to 259 black hole merger events from the GWTC-5 catalog, building a mixture model of two subpopulations. They recovered separate mass and spin distributions and found nearly identical shapes between the mass function of the high-spin subpopulation and the final mass distribution of low-spin black holes after merger, supporting a hierarchical origin.
Where its strength is
The morphological match of the mass functions is a direct and striking statistical argument for hierarchical mergers, independent of the details of cluster dynamics modeling. The obtained constraint on the mass gap and the S-factor of the 12C(α,γ)16O reaction links gravitational-wave observations to nuclear astrophysics.
Where it could go next
- One could compare the recovered subpopulations with results from stellar cluster dynamics simulations, where hierarchical mergers leave characteristic mass–spin correlations.based on: The Smoking Gun of Hierarchy: Black Holes Betray Their Ancestors arXiv:2607.01121v1
- It would be interesting to check how the model's predictions for the effective spin distribution align with independent estimates from individual event analyses, especially for high-mass events.based on: Twin Black Holes: Why Light Ones Are Lazy and Heavy Ones Are Tops arXiv:2605.24281
- The analysis could be extended to intermediate-mass events to trace the continuity of the hierarchical chain and possible transitions between subpopulations.based on: Four Layers of Black Holes: Gravitational Waves Reveal Merger History arXiv:2607.22011v1
Worth trying
- Test whether the mass function match persists under alternative spline parameterizations.
- Assess how detector calibration uncertainties affect the conclusion of morphological similarity.
- Compare the obtained S-factor with results from recent laboratory experiments.
Why this may matter
If the hierarchical origin of high-spin black holes is confirmed, it will provide a new tool for studying the dynamics of dense stellar clusters and allow the mass gap to be used to refine nuclear reaction rates in massive stars.
Nearby works in our archive
- Black Holes Are Sorted by Weightsimilarity 0.71
- The Smoking Gun of Hierarchy: Black Holes Betray Their Ancestorssimilarity 0.70
- A Gravitational Census of Invisible Giantssimilarity 0.69
- Twin Black Holes: Why Light Ones Are Lazy and Heavy Ones Are Topssimilarity 0.69
- Why Black Holes in a Pair Look So Alikesimilarity 0.69
- Cosmic Scales: How Black Holes Swap Rolessimilarity 0.69
- Four Layers of Black Holes: Gravitational Waves Reveal Merger Historysimilarity 0.68
The archive is crowded around this work: 9 reviews nearby, closest at 0.71. This ground is well walked — and there is plenty to compare against.
The pair-instability supernova mass gap, where black holes are almost never born, was long considered a 'forbidden zone,' but now we see how hierarchical mergers fill it with remnants of previous generations, like a cosmic nesting doll.
The idea of hierarchical mergers echoes the 'nesting-doll universe' from Liu Cixin's 'The Three-Body Problem,' where civilizations embed worlds within each other, but here nature itself builds increasingly massive black holes through a cascade of collisions.
Key numbers
- Bhattacharyya coefficient: 0.95
- maximum mass of low-spin black holes: 54.2 M☉
- S-factor for ¹²C(α,γ)¹⁶O at 300 keV: 151 keV·barn
- number of BBH events in GWTC-5: 259
- fraction of hierarchical mergers: ~10%