⚛️ general relativity

Self-Similar Mass Spectra of Hierarchical Black Hole Mergers

This paper models hierarchical black hole mergers as a coagulation problem to derive closed scaling laws for mass spectra, revealing that most formation channels lead to runaway growth (gelation) while only the early three-body channel allows for slow self-similar growth, with numerical results showing that energy radiated during mergers significantly alters the high-mass cutoff and total population mass.

Original authors: Liam Blum, Stefano Profumo, Ryan Schantz, Quoc Ha Tran

Published 2026-09-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Liam Blum, Stefano Profumo, Ryan Schantz, Quoc Ha Tran

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Black holes are not just the silent, invisible anchors of galaxies; they are also the ultimate recyclers of the cosmos. When two black holes collide, they do not simply stick together like clay. Instead, they spiral into one another, merging into a single, larger object while blasting a tremendous amount of energy away in the form of ripples in space-time known as gravitational waves. This process means that the new black hole is actually lighter than the two that created it, having lost a few percent of its total mass to the universe. Over billions of years, these collisions can happen again and again, building up a population of black holes that are far more massive than any single star could ever produce. Astronomers have recently begun to detect these heavyweights, including some that seem too large to have formed from a dying star, raising a fundamental question: how does the mass of these objects change as they merge repeatedly, and what does the final distribution of their sizes look like?

To answer this, a team of researchers at the University of California, Santa Cruz, treated the history of black hole mergers not as a series of unique cosmic events, but as a statistical problem of how particles clump together over time. They applied a mathematical framework originally designed to describe how tiny particles in a fluid stick to one another, adapting it to the violent, high-speed collisions of black holes. The researchers focused on a specific type of growth called hierarchical merging, where a black hole formed from a previous collision goes on to merge again. They asked whether these populations settle into a predictable pattern, a steady state where the number of small, medium, and large black holes follows a specific rule, or if the process spirals out of control, with the heaviest objects consuming everything else.

The team discovered that the answer depends entirely on the environment where the black holes live and how the merger rate changes with the size of the colliding objects. In some scenarios, particularly those involving primordial black holes—hypothetical objects formed in the very first moments of the universe—the rules of the game change drastically. The researchers found that for most of the standard ways these black holes are thought to form, the process favors the heavyweights so strongly that a stable, balanced population cannot exist. In these cases, the system is prone to "runaway growth," where the largest black holes keep growing faster and faster, eventually dominating the entire population and preventing a steady mass distribution from ever forming. This suggests that if we see a lot of heavy black holes, it is not because the population has settled into a calm equilibrium, but because the growth process is still accelerating.

However, the story is different for one specific formation channel: the early three-body interaction. In this scenario, three black holes interact in the early universe, and two of them pair up while the third is ejected. The researchers found that this specific pathway leads to a slow, steady, and predictable growth. Here, the population does approach a stable, self-similar shape where the distribution of masses follows a clear pattern. This is the only channel among the four they studied that allows for a calm, mass-conserving evolution where the total weight of the black holes remains roughly constant over time, aside from the small losses due to gravitational waves.

A critical part of their work involved correcting a common misunderstanding in how scientists translate observed merger rates into these mathematical models. Previous studies often treated the rate of collisions as a direct measure of how likely two specific black holes are to merge, but the researchers showed that this approach misses a crucial factor: the way the rate is reported in scientific literature includes a hidden multiplication by the masses of the objects involved. When they corrected for this, the mathematical description of the merger process shifted significantly. This correction was the key that reclassified the behavior of the different formation channels, moving three of them from a stable category into the runaway growth category.

The team also investigated the effect of the mass lost during each collision. Because black holes radiate away energy, the total mass of the population decreases with every merger. The researchers simulated this loss and found that it measurably changes the shape of the mass distribution. It does not simply shrink the population uniformly; instead, it shifts the cutoff point where the distribution drops off and alters the number of very massive black holes compared to the lighter ones. Their simulations showed that accounting for this energy loss is essential for accurately predicting what astronomers should see in the sky. Without including this loss, the models would predict a different shape for the mass spectrum than what actually occurs in a non-conservative system.

Ultimately, the paper provides a new, clearer lens through which to view the cosmic history of black holes. It establishes that the universe does not treat all black hole mergers the same way. While some environments allow for a slow, steady buildup of mass that follows a predictable mathematical curve, others drive the system toward a chaotic state where the heaviest objects take over. By distinguishing between these regimes and correcting how the underlying physics is calculated, the researchers have provided a more accurate foundation for interpreting the growing catalog of gravitational wave detections. Their work suggests that the heavy black holes we are seeing are likely the result of a dynamic, accelerating process rather than a static, balanced one, offering a new perspective on how the most extreme objects in the universe come to be.

Technical Summary: Self-Similar Mass Spectra of Hierarchical Black Hole Mergers

Problem Statement
The hierarchical merger of black holes (BHs) is a critical mechanism for explaining observed gravitational-wave events, such as GW190521, which lie within the pair-instability supernova (PISN) mass gap. While numerical simulations (e.g., COMPAS, CMC) and semi-analytic models have elucidated specific environmental dependencies, a systematic analytic treatment of how the BH mass distribution evolves under repeated mergers remains underdeveloped. Specifically, there is a lack of explicit scaling laws and self-similar solutions that can classify merger dynamics into universality classes and map kernel parameters directly to observable spectral features (e.g., power-law slopes, high-mass tails). Furthermore, standard coagulation theory assumes strict mass conservation, whereas BH mergers radiate a significant fraction of mass ( $\sim$ 5%) as gravitational waves, breaking the conservation assumption and altering the scaling relations.

Methodology
The authors treat hierarchical merging as a coagulation problem governed by Smoluchowski's coagulation equation. They adopt a controlled homogeneous kernel class:
$K(m, m', t) = (m + m')^{-\alpha} (mm')^{-\beta} t^{-\delta}$
where $\alpha, \beta, \delta$ encode dependence on total mass, component-mass bias, and cosmic time, respectively. The degree of homogeneity is defined as $\lambda = -(\alpha + 2\beta)$ .

The methodology proceeds in three stages:

Analytic Derivation: The authors derive the self-similar scaling form $c(m, t) = s(t)^{-2} \Phi(m/s(t))$ for the conservative case ( $L=1$ , where $L$ is the retained mass fraction). They establish the relationship between the characteristic mass growth exponent $z$ and the kernel parameters, noting that for $\delta < 1$ , $z = (1-\delta)/(1+\alpha+2\beta)$ . They also derive the exact similarity equation for the scaling function $\Phi$ and provide an exact benchmark solution for the constant kernel ( $K=1$ ).
Kernel Mapping: A critical technical contribution is the rigorous conversion of published primordial black hole (PBH) merger rates (typically quoted as rate densities in terms of total mass $M$ and symmetric mass ratio $\eta$ ) into the coagulation kernel formalism. The authors demonstrate that treating rate densities as direct kernels is incorrect; one must account for the conversion factor $m_1 m_2 = \eta M^2$ . This shifts the homogeneity degree by two units, fundamentally altering the classification of PBH formation channels.
Numerical Simulation: Using kinetic Monte Carlo methods, the authors simulate the discrete coagulation dynamics for 13 distinct kernel configurations spanning suppressive regimes ( $\lambda \le 0$ ). They compare "elastic" (conservative, $L=1$ ) and "inelastic" ( $L=0.95$ ) scenarios to quantify the impact of gravitational-wave mass loss. Simulations are stopped at a fixed surviving number of objects to ensure a fair comparison of profile shapes at equivalent stages of coarsening.

Key Contributions and Results

Reclassification of PBH Channels: The paper corrects the mapping of PBH merger channels to kernel exponents.
- Early Two-Body (E2), Late Two-Body (L2), and Late Three-Body (L3): These channels map to superlinear kernels ( $\lambda > 1$ ). Specifically, E2 yields $\lambda \approx 1.135$ , L2 yields $\lambda = 2$ , and L3 yields $\lambda > 1$ . These are identified as candidates for "gelation" or runaway growth, where a macroscopic object absorbs a finite fraction of the mass, and no steady mass-conserving spectrum exists.
- Early Three-Body (E3): This is the only channel mapping to the nongelling regime ( $\lambda < 1$ ), admitting a slow self-similar growth with $z \approx 0.0305$ (for $\gamma=1$ ) to $0.0716$ (for $\gamma=2$ ).
Non-Universality of Profile Shape: Numerical results demonstrate that the shape of the mass spectrum is not fixed by the homogeneity degree $\lambda$ alone. Kernels with identical $\lambda$ but different distributions of suppression (e.g., total-mass suppression vs. product-mass suppression) yield measurably different spectral shapes (different algebraic prefactors $p$ and tail exponents $q$ ).
Impact of Gravitational-Wave Mass Loss: The inclusion of a 5% mass loss ( $L=0.95$ $L = 0.95$ ) systematically alters the solution:
- The tail exponent $q$ increases at every sampled kernel (by $\sim 0.11$ along the $\alpha$ axis and up to $0.14$ along the $\beta$ axis), indicating a steeper high-mass cutoff.
- The characteristic scale $\xi_0$ increases.
- The conservative relation between the amplitude exponent $\theta$ and growth exponent $z$ ( $\theta = 2z$ ) breaks down. In the non-conservative case, $\theta$ and $z$ must be measured independently as the total mass decays.
Exact Benchmark: The paper provides an exact solution for the constant kernel ( $K=1$ ), showing that the discrete spectrum approaches a continuous exponential profile $\Phi(\xi) = M_1 e^{-\xi}$ , which serves as a validation point for the numerical methods.

Significance and Claims
The paper claims to provide a "common similarity framework" for hierarchical astrophysical and primordial black hole mergers. Its primary significance lies in:

Analytic Clarity: Establishing that the standard self-similar theory relies on mass conservation, which is violated in BH mergers, necessitating a distinct treatment for the radiating case.
Corrected Classification: Demonstrating that previous interpretations of PBH merger rates may have misclassified channels as nongelling when they are actually superlinear ( $\lambda > 1$ ) and prone to runaway growth. The authors emphasize that while $\lambda > 1$ is a diagnostic for potential gelation, physical gelation depends on the full kernel behavior at unequal masses and the finite coagulation history, which are not fully settled here.
Numerical Insight: Providing the first systematic quantification of how gravitational-wave mass loss reshapes the mass spectrum, showing that it cannot be treated as a simple rescaling but fundamentally changes the spectral shape and cutoffs.

The authors remain modest regarding the "runaway" regime, noting that while their classification identifies E2, L2, and L3 as candidates for gelation, establishing physical gelation requires the full population-dependent kernel and finite coagulation history, which are beyond the scope of the current analytic and numerical snapshots. The work serves as a benchmark for numerical codes and a guide for interpreting future gravitational-wave population data from LVK and next-generation detectors.

More like this