Reexamining Evidence of a Pair-Instability Mass Gap in the Binary Black Hole Population

Using flexible models on the GWTC-4 catalog, this study finds no current evidence for a sharp pair-instability mass gap near 40–50 solar masses, suggesting that high-mass binary black holes are likely not dominated by hierarchical mergers and instead constraining the potential gap's lower bound to approximately 57 solar masses.

The Gist

The Big Picture: The "Missing Middle" Mystery

Imagine the universe as a giant gym where black holes are the weightlifters. For a long time, physicists believed there was a specific weight class that simply didn't exist. They thought that if a star tried to grow too heavy (between 40 and 70 times the mass of our Sun), it would explode so violently that it would shatter itself into pieces, leaving no black hole behind.

This theoretical "forbidden zone" is called the Pair-Instability Mass Gap. It's like a gym rule that says, "No one can lift between 400 and 700 pounds."

Recently, scientists started finding black holes that should have been in this forbidden zone. This led to a big debate:

1. The Old Theory: Maybe the rule is real, and these heavy black holes are actually "second-generation" lifters—created when two smaller black holes crashed into each other and merged.
2. The New Theory (This Paper): Maybe the rule doesn't exist at all, or it's just much higher up than we thought.

The Investigation: Re-weighing the Evidence

The authors of this paper, Anarya Ray and Vicky Kalogera, took a fresh look at the latest data from the GWTC-4 catalog (a list of 153 black hole collisions detected by gravitational wave sensors).

Think of previous studies as trying to find a specific type of fish in a lake by using a net with very large, rigid holes. If the fish didn't fit through the holes, the researchers assumed the fish didn't exist. They found a "gap" in the data and concluded, "Aha! The heavy black holes are missing because they are second-generation mergers!"

This paper says: "Wait a minute. Let's use a flexible net instead."

Instead of forcing the data to fit a specific shape (a sharp gap), the authors used flexible models. They let the data tell the story without assuming the "gap" was there to begin with.

The Findings: No Sharp Cliff, Just a Gentle Slope

Here is what they found when they used their flexible net:

1. No Sharp Drop-Off: The data doesn't show a sudden cliff where black holes stop existing at 40–50 solar masses. Instead, it looks more like a gentle hill. The number of heavy black holes gets smaller as they get heavier, but they don't just vanish.

   • Analogy: Imagine a staircase. Previous studies thought there was a step where the stairs just ended abruptly. This paper says, "Actually, it's just a ramp that gets steeper and steeper, but there are still people walking up it."

2. The "Second-Generation" Theory is Weak: If those heavy black holes were made by merging two smaller ones (a "2G+1G" merger), we would expect them to be very lopsided (one heavy, one light) and spinning in weird, chaotic directions.

   • The Twist: The heavy black holes the authors found are actually very balanced (the two black holes are similar in size) and have spins that don't quite match the "chaotic merger" prediction.
   • Analogy: If you found a pile of bricks that were supposed to be made by smashing two small bricks together, you'd expect a messy, uneven pile. Instead, you found a perfectly symmetrical tower. This suggests they weren't made by smashing; they might have grown naturally.

3. The Gap Might Be Higher Up: The authors aren't saying the gap doesn't exist at all. They are saying, "If there is a gap, it's probably much higher up, maybe around 60 solar masses or more."

   • Analogy: We thought the "no-fly zone" started at 10,000 feet. The data suggests the zone might actually start at 20,000 feet. We just haven't flown high enough yet to see the ceiling.

Why Does This Matter? (The Nuclear Recipe)

Why do we care about the exact weight where black holes stop forming? Because it tells us about the recipe of the universe.

The existence of this gap depends on a specific nuclear reaction inside dying stars (how Carbon and Oxygen interact).

• If the gap is low (40–50 solar masses), the nuclear reaction rate must be high.
• If the gap is high (57+ solar masses), the reaction rate is lower.

By finding that the gap is likely higher (or non-existent in the 40–50 range), the authors are telling nuclear physicists: "Your recipe needs a little less of this ingredient."

The Takeaway

The paper is essentially a "reality check" for the astrophysics community.

• Previous View: "We see a gap, so heavy black holes must be rare, second-generation monsters."
• This Paper's View: "We don't see a sharp gap. The heavy black holes are there, they are balanced, and they might just be normal stars that grew very large. We need more data to be sure, but we shouldn't force the data to fit a theory that might be wrong."

In short: The universe might not have a "forbidden zone" for black holes between 40 and 50 suns. The heavyweights are still in the gym, lifting weights, and they look more like natural athletes than the "Frankenstein monsters" we thought they were.

Technical Summary

1. Problem Statement

The paper addresses the ongoing debate regarding the existence and location of the Pair-Instability Supernova (PISN) mass gap in the binary black hole (BBH) mass spectrum.

• Theoretical Context: Stellar evolution models predict a gap in black hole masses between approximately $40-70 M_\odot$ and up to $\sim130 M_\odot$ due to the complete disruption of progenitor stars by pulsational pair-instability supernovae (PISN).
• Recent Claims: Previous studies analyzing the fourth gravitational wave transient catalog (GWTC-4), specifically H. Tong et al. (2025), claimed evidence of a sharp drop-off (gap) in the secondary black hole mass spectrum starting at $m_g \approx 45 M_\odot$. They attributed the presence of black holes above this mass to a sub-population of hierarchical mergers (specifically 2G+1G systems, where a second-generation black hole merges with a first-generation one).
• The Conflict: These claims rely heavily on specific prior assumptions regarding the mass-ratio distribution and the width of the gap. Ray and Kalogera question whether these findings are robust or if they are artifacts of "strong model assumptions" (biased priors) that force a gap-like feature into the data.

2. Methodology

The authors employ Bayesian hierarchical inference to re-analyze the GWTC-4 dataset (153 events) using flexibly parametrized models designed to avoid imposing sharp features a priori.

• Data: 153 GWTC-4 events with a false-alarm rate of $\le 1$ per year. They utilize parameter estimation (PE) posteriors and ranked simulated signals (injections) to correct for Malmquist biases.
• Primary Mass Model: They utilize a Broken Power Law with two peaks (BPL2P) for the primary mass distribution, consistent with recent literature, and test an alternative Single Power Law with two peaks (SPL2P).
• Flexible Mass-Ratio Distribution: Instead of assuming a specific gap shape, they model the mass-ratio distribution $p(q|m_1)$ as a mixture of:
  1. A truncated Gaussian ($TN_q$).
  2. A continuous broken power-law ($BPL_q$).
  • This allows the data to determine if a sharp break (gap) exists or if the distribution is smooth.
• Spin Transition Modeling: They simultaneously model transitions in effective aligned spin ($\chi_{eff}$) and effective precessing spin ($\chi_p$) at a transition mass $\tilde{m}$. They test if the sub-population above $\tilde{m}$ exhibits the specific spin characteristics predicted for hierarchical mergers (broad spin distributions, low alignment).
• Model Comparison: They perform rigorous model selection using Bayes Factors (BF) to compare:
  1. Flexible Models: No enforced gap.
  2. Gap-Enforced Models: Strong priors forcing a gap (mimicking Tong et al.).
  3. Tong et al. (2025) Models: Direct comparison with the previous study's specific parametrization.
• Cluster Simulations: They compare their inferred distributions against Cluster Monte Carlo (CMC) simulations of dense stellar clusters to test the consistency of the "2G+1G hierarchical merger" hypothesis.

3. Key Contributions

• Rejection of the Sharp Gap at 40–50 $M_\odot$: The study provides strong evidence that GWTC-4 data does not support a sharp drop-off in the secondary mass spectrum near $40-50 M_\odot$.
• Demonstration of Prior Bias: They show that the "gap" reported in previous studies is largely driven by restrictive priors (e.g., enforcing a gap width $>20 M_\odot$). When using flexible priors, the data prefers a smooth fall-off.
• Refutation of the 2G+1G Dominance: Even if a spin transition exists, the mass-ratio distribution of the high-mass sub-population is inconsistent with the 2G+1G hierarchical merger scenario. The data shows a significant preference for symmetric mass binaries ($q \in 0.6-1$), whereas 2G+1G models predict highly asymmetric systems.
• New Constraints on PISN Physics: By allowing for a gap only at higher masses, they provide updated constraints on the nuclear reaction rate $S$-factor for $^{12}C(\alpha, \gamma)^{16}O$.

4. Key Results

A. Mass Distribution

• Smooth Fall-off: The inferred secondary mass distribution shows a smooth fall-off rather than a sharp gap. The break location is constrained to $m_{2,b1} = 40.0^{+9}_{-7} M_\odot$.
• Bayes Factors: The flexible model is preferred over the gap-enforced model by a factor of $\sim3$ (BF $\approx 191$ vs. 60 for BPL2P; BF $\approx 353$ vs. 100 for SPL2P).
• Gap Location: If a PISN gap exists, its lower bound is constrained to be higher: $57^{+17}_{-10} M_\odot$. This implies the $S$-factor for the $^{12}C(\alpha, \gamma)^{16}O$ reaction is $\le 101^{+87}_{-23}$ keV barns, which overlaps significantly with independent nuclear physics measurements (unlike the marginal overlap found in gap-enforced studies).

B. Spin and Mass-Ratio Transitions

• Spin Transition: A transition in spin distributions is confirmed at $\tilde{m} = 43^{+11}_{-5} M_\odot$. Above this mass, the distribution of $\chi_{eff}$ and $\chi_p$ broadens, consistent with dynamical assembly.
• Mass-Ratio Inconsistency: Despite the spin transition, the mass-ratio distribution above $\tilde{m}$ is inconsistent with 2G+1G hierarchical mergers.
  • Finding: $52^{+18}_{-23}\%$ of the high-mass sub-population has mass ratios $q \in [0.6, 1]$.
  • Implication: 2G+1G systems are expected to be highly asymmetric ($q \ll 1$). The observed symmetry suggests that 2G+1G mergers do not dominate this population. Contributions from 2G+2G systems are insufficient to resolve this discrepancy.

C. Alternative Explanations

Since the 2G+1G scenario is disfavored, the authors propose alternative formation channels for the high-mass, symmetric, broad-spin sub-population:

1. BH-Star Collisions: Dynamical mergers involving black holes accreting from stars (1G+1G), which can produce equal-mass, high-spin systems.
2. Isolated Binary Evolution: Scenarios like super-Eddington accretion or chemically homogeneous evolution.
3. Mixed Populations: A combination of isolated and hierarchical channels.

5. Significance

• Astrophysical Interpretation: The paper challenges the consensus that the high-mass BBH population in GWTC-4 is dominated by hierarchical mergers in dense clusters. It suggests that isolated binary evolution or accretion-driven growth may play a more significant role than previously thought.
• Nuclear Physics: The revised lower bound on the PISN cutoff ($\sim57 M_\odot$) provides a constraint on the $^{12}C(\alpha, \gamma)^{16}O$ reaction rate that is fully consistent with independent laboratory measurements, whereas previous "gap-enforced" analyses were only marginally consistent.
• Methodological Impact: The study serves as a critical warning against prior-induced biases in gravitational wave population inference. It demonstrates that enforcing "physically motivated" but restrictive priors (like a specific gap width) can lead to incorrect astrophysical conclusions.
• Future Outlook: The authors emphasize that larger sample sizes (future GW catalogs) and non-parametric inference methods are necessary to definitively locate the PISN gap and disentangle the relative contributions of different formation channels.

In summary, Ray and Kalogera argue that the current evidence for a sharp PISN mass gap at $40-50 M_\odot$ is an artifact of modeling assumptions, and the data actually supports a smoother mass distribution with a potential gap onset at significantly higher masses, driven by formation channels other than simple 2G+1G hierarchical mergers.