Evidence of the pair instability gap in the distribution of black hole masses

This study analyzes the GWTC-4 catalog to identify a pair-instability mass gap in the distribution of secondary black hole masses, providing evidence for a subpopulation of hierarchical mergers and constraining the nuclear reaction rate of $^{12}\rm{C}(\alpha,\gamma)^{16}\rm{O}$.

The Gist

Imagine the universe as a giant cosmic bakery. For decades, astrophysicists have had a very specific recipe for how stars (the bakers) turn into black holes (the bread). According to this recipe, there is a "forbidden zone" for the size of the bread.

If a star is just right, it bakes a small loaf. If it's huge, it bakes a giant loaf. But if a star is in a specific, middle-heavy range (between about 50 and 130 times the mass of our Sun), the recipe says the oven explodes. The star blows itself apart completely, leaving no bread at all. This is called the Pair-Instability Gap.

For a long time, we looked for this missing bread in the cosmic pantry, but we couldn't find it. We kept finding loaves that were too big, making us think the recipe was wrong.

The Big Discovery
In this new paper, a team of scientists (led by Hui Tong) looked at the latest batch of data from the LIGO, Virgo, and KAGRA gravitational wave detectors. Think of these detectors as ultra-sensitive microphones listening to the "thuds" of black holes crashing into each other.

They didn't just look at the biggest black holes in the crash (the "Primary" black hole); they also looked closely at the smaller partner (the "Secondary" black hole).

The "Missing Middle" Analogy
Imagine you are sorting a pile of marbles by size.

• The Primary Mass (The Big Marble): You see marbles of all sizes, from tiny to huge. There's no obvious gap.
• The Secondary Mass (The Small Marble): Suddenly, you notice something strange. There are no marbles between the size of a pea and a golf ball. There is a clear, empty space in the middle of your pile.

That empty space is the Pair-Instability Gap. The scientists found that while the "big" black holes in a pair can be any size, the "small" partner black holes strictly avoid that forbidden middle zone. It's like a cosmic bouncer who says, "You can be small, or you can be huge, but you cannot be in the middle."

Why is the gap there? The "Recycled" Black Holes
If the gap exists, why do we still see black holes in that size range sometimes? The paper suggests these are "recycled" black holes.

Think of it like a video game:

1. Level 1 (Normal Stars): A star dies and becomes a black hole. This is a "Level 1" black hole. It follows the rules and stays out of the forbidden zone.
2. Level 2 (The Mergers): Two Level 1 black holes crash together. The result is a new, bigger black hole. This "Level 2" black hole is a "child" of a merger. Because it was built from two smaller pieces, it can end up in that forbidden middle zone where normal stars can't go.

The scientists found that the black holes in the forbidden zone are almost always the "Level 2" ones—the children of previous crashes. This explains why the gap is visible in the "smaller" partner (because the "big" partner is often the recycled one that jumped the gap).

The Spin Connection
To prove this theory, the scientists looked at how fast these black holes were spinning.

• Normal black holes spin slowly, like a lazy top.
• Recycled black holes (the ones that jumped the gap) should spin very fast, like a figure skater pulling in their arms, because they inherited the spin from their parents.

The data showed exactly this: The black holes in the forbidden zone were spinning much faster than the ones below it. This confirmed that the gap is real and that the "forbidden" black holes are indeed the recycled, "Level 2" champions.

Why Does This Matter?
This discovery is like finding a missing page in the universe's instruction manual.

1. It confirms the recipe: It proves that our theories about how stars explode are correct.
2. It measures nuclear physics: By measuring exactly where the gap starts (around 45 times the mass of the Sun), the scientists can calculate how fast carbon and oxygen atoms fuse inside dying stars. It's like using the size of the missing bread to figure out the exact temperature of the oven.

In a Nutshell
The universe has a "No-Go Zone" for black holes born from normal stars. We finally found the empty space in the data. The black holes that do appear in that zone are the cosmic equivalent of "Frankenstein monsters"—built from the wreckage of previous black hole collisions. This discovery helps us understand how stars die, how black holes are born, and even how the atoms in our own bodies were forged in the hearts of stars.

Technical Summary

1. Problem Statement

Stellar evolution theory predicts a "pair-instability supernova" (PISN) gap in the black hole (BH) mass spectrum. Stars with initial masses between approximately 100 and 260 $M_\odot$ are expected to undergo pair-instability processes where photon pressure drops due to electron-positron pair production, leading to a complete disruption of the star or a pulsational mass loss that prevents the formation of a black hole in the 50–130 $M_\odot$ range.

Despite this robust theoretical prediction, observational evidence for this gap in gravitational wave (GW) data has been elusive. Previous analyses of GW catalogs (GWTC-1 to GWTC-3) failed to find a clear gap, partly because:

• More massive binary black holes (BBHs) were discovered, ruling out a sharp cutoff in the primary mass ($m_1$) distribution.
• A suspected peak at $\sim35 M_\odot$ (attributed to pulsational PISN) was found to be in tension with stellar physics and supernova rates.
• Previous studies focused on the primary mass distribution ($m_1$), potentially missing the gap if it is obscured by "hierarchical mergers" (where a black hole formed from a previous merger populates the gap).

2. Methodology

The authors analyzed data from the fourth Gravitational Wave Transient Catalog (GWTC-4), comprising roughly 153 confident binary black hole events (including the first part of the O4a observing run).

Key Analytical Steps:

• Hierarchical Bayesian Inference: Used the GWPopulation framework to infer the underlying population distribution $\pi(m_1, m_2)$ from individual event posteriors.
• Secondary Mass Focus: Unlike previous studies focusing on $m_1$, this analysis specifically modeled the secondary mass ($m_2$) distribution. The hypothesis is that while $m_1$ may be contaminated by hierarchical mergers (2G+1G), $m_2$ in a binary is less likely to be a product of a previous merger, thus preserving the PISN gap.
• Gap Model: The mass ratio distribution $\pi(q|m_1)$ was modified to include a "forbidden" interval where the probability density is zero:
  $$\pi(q|\Lambda) \propto \begin{cases} 0 & \text{if } m_g \le q m_1 \le m_g + w_g \\ q^{\beta_q} & \text{otherwise} \end{cases}$$
  Where $m_g$ is the lower edge and $w_g$ is the width of the gap.
• Spin Correlation: The study incorporated a mass-dependent spin model (from Ref. [30]) to test for a transition in the effective inspiral spin ($\chi_{eff}$) at the same mass scale as the gap. Hierarchical mergers are expected to have high spins ($\chi \approx 0.7$).
• Model Comparison: The authors compared the "gap" model against a "no-gap" model using Bayes factors. They also tested different mass distribution models (Single Power Law + Two Peaks vs. Broken Power Law + Two Peaks) and a "pairing mass" model to ensure robustness.
• Non-Parametric Verification: Used the "$\pi$-stroke" (maximum population likelihood) method, which is model-free, to verify that the gap emerges organically from the data without parametric assumptions.

3. Key Contributions

• First Robust Detection of the PISN Gap: Provided the first statistically significant evidence of the pair-instability gap, specifically located in the secondary mass ($m_2$) distribution rather than the primary.
• Hierarchical Merger Interpretation: Demonstrated that the gap in $m_2$ is consistent with a population where the primary mass ($m_1$) is often a "second-generation" (2G) black hole formed from a previous merger, while the secondary ($m_2$) remains a "first-generation" (1G) star-collapse remnant.
• Spin-Mass Correlation: Confirmed that the mass scale where the gap begins aligns with a transition in the spin distribution. Binaries with $m_1$ above the gap threshold exhibit significantly higher spins, consistent with the hierarchical merger hypothesis.
• Nuclear Physics Constraints: Used the measured location of the gap to constrain the astrophysical S-factor for the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ nuclear reaction rate, a critical uncertainty in stellar evolution models.

4. Key Results

• Gap Location:
  • Lower Edge ($m_g$): $45^{+5}_{-4} M_\odot$ (90% credibility) when the gap lower edge is constrained to be identical to the spin transition mass.
  • Upper Edge: Constrained to $116^{+9}_{-13} M_\odot$, though this is heavily dependent on the single most massive event detected, GW231123.
  • Robustness: The lower edge is a robust feature; excluding GW231123 yields a lower edge of $43^{+7}_{-6} M_\odot$.
• Statistical Significance: The gap model is preferred over the no-gap model by a Bayes factor of $\sim 10^2$.
• Spin Transition: The spin transition mass ($\tilde{m}_1$) is found to be consistent with the lower edge of the $m_2$ gap ($\sim 45 M_\odot$). Binaries with $m_1 > 45 M_\odot$ show a transition to higher effective spins, supporting the 2G+1G merger hypothesis.
• Event Classification:
  • Identified four events (e.g., GW190519, GW190602) with strong support ($\ln B > 8$) for containing a primary black hole inside the gap (likely 2G+1G).
  • GW190521: The authors argue GW190521 is likely a "straddling binary" (primary above the gap, secondary below) rather than a 2G+2G merger, as its secondary mass is consistent with being below the gap.
• Merger Rates:
  • Lower limit on 2G+1G merger rate: $2.5^{+2.2}_{-1.2} \times 10^{-1} \text{ Gpc}^{-3} \text{ yr}^{-1}$.
  • Upper limit on 2G+2G merger rate (both components in the gap): $< 7.9 \times 10^{-2} \text{ Gpc}^{-3} \text{ yr}^{-1}$.
• Nuclear Physics Constraint: The measurement constrains the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ S-factor at 300 keV to $256^{+197}_{-104}$ keV barns (90% credibility).

5. Significance

• Validation of Stellar Theory: Confirms a long-standing prediction of stellar evolution regarding the suppression of black hole formation in the 50–130 $M_\odot$ range.
• Insight into Black Hole Formation: Resolves the "missing gap" paradox in previous catalogs by showing that the gap exists in $m_2$ but is masked in $m_1$ by hierarchical mergers. This provides strong evidence for the existence of second-generation black holes in dense stellar environments.
• Cosmological Utility: The pair-instability gap is predicted to be a stable feature across cosmic time (unlike the uncertain 35 $M_\odot$ bump). This makes it a potential "standard siren" for measuring the Hubble constant and breaking distance-redshift degeneracies.
• Nuclear Astrophysics: Provides an independent, astrophysical constraint on the $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ reaction rate, which is crucial for understanding stellar nucleosynthesis and the evolution of massive stars.

In conclusion, this paper establishes the pair-instability gap as a real feature of the black hole population, located specifically in the secondary mass distribution, and leverages this finding to simultaneously constrain black hole formation channels (hierarchical mergers) and fundamental nuclear physics.