Measurement prospects for the pair-instability mass cutoff with gravitational waves
Pith reviewed 2026-05-16 05:33 UTC · model grok-4.3
The pith
Gravitational wave data similar to GWTC-4 cannot yet guarantee detection of the pair-instability black hole mass cutoff.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using full Bayesian parameter estimation on simulated catalogs with and without a mass cutoff, the analysis shows that for GWTC-4-like volumes confident identification is not guaranteed, but strongest constraints are compatible with real-data results when the model includes the cutoff; spurious identification is unlikely. For O4 catalogs uncertainty on the cutoff mass reduces by at least 20 percent, giving S_300 greater than or equal to 125 keV b at 90 percent credibility for a 40-50 solar mass cutoff, while relative uncertainties on the Hubble parameter from gravitational waves alone can reach 100 percent. Nonparametric PixelPop modeling of GWTC-4 shows prominent mass features but does not
What carries the argument
Full Bayesian inference on simulated catalogs using parametric population models that include or exclude a sharp mass cutoff, validated against the nonparametric PixelPop model applied to GWTC-4 data.
If this is right
- Uncertainty on the cutoff mass decreases by at least 20 percent with O4 data volumes
- Only a lower bound on the carbon-alpha-oxygen reaction rate is obtained for cutoffs at 40-50 solar masses
- Hubble parameter constraints from gravitational waves alone remain uncertain by up to 100 percent
- The nonparametric PixelPop model applied to GWTC-4 reveals mass features without requiring a sharp cutoff, consistent with simulations that include one
Where Pith is reading between the lines
- Astrophysical claims about pair-instability supernovae from gravitational-wave populations require explicit model-comparison and predictive tests to be robust.
- Higher-sensitivity detectors or longer observing runs will likely be needed to tightly locate the cutoff and derive strong constraints on stellar evolution.
- The current compatibility with GWTC-4 results leaves open the possibility that observed mass features arise from the pair-instability gap, but more events are required to distinguish this from other mass-distribution effects.
Load-bearing premise
The parametric population models and simulated selection effects accurately capture the true underlying black-hole mass distribution and detection biases.
What would settle it
A future observed catalog in which the posterior shows no support for a cutoff when the simulated truth includes one, or a strong preference for a cutoff when the truth has none.
Figures
read the original abstract
Pair-instability supernovae leave behind no compact remnants, resulting in a predicted gap in the distribution of stellar black-hole masses. Gravitational waves from binary black-hole mergers probe the relevant mass range and analyses of the LIGO-Virgo-KAGRA catalog (GWTC-4) indicate a possible mass cutoff at $40$-$50\,M_\odot$. However, the robustness of this result remains unclear. To this end, using full Bayesian parameter estimation, we simulate gravitational-wave catalogs with and without such a mass cutoff, then test whether its presence or absence is correctly inferred with parametric population models. For catalogs similar to GWTC-4, confident identification of a cutoff is not guaranteed, but the best constraints among our simulations are compatible with results from GWTC-4 when the model includes a cutoff. Conversely, spurious identification of a cutoff is unlikely. For catalogs expected by the end of the O4 observing run, uncertainty in the cutoff mass is reduced by $\gtrsim20\,\%$, but a cutoff at 40-50$M_\odot$ yields only a lower bound on the $^{12}\mathrm{C}(\alpha,\gamma)^{16}\mathrm{O}$ reaction rate, our most stringent constraints on the S-factor at $300\,\mathrm{keV}$ being $S_{300}\gtrsim125\,\mathrm{keV}\,\mathrm{b}$ at $90\,\%$ credibility. Relative uncertainties on the Hubble parameter $H_0$ from gravitational-wave data alone can still be up to $100\,\%$. We also analyze GWTC-4 with the nonparametric PixelPop population model, finding that some mass features are more prominent than in parametric models but a sharp cutoff is not required. However, the parametric model passes a likelihood-based predictive test in GWTC-4 and the PixelPop results are consistent with those from our simulated catalogs with a cutoff. Such tests are necessary to make astrophysical claims from gravitational-wave catalogs.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript simulates gravitational-wave catalogs with and without a pair-instability mass cutoff at 40-50 solar masses, then recovers the catalogs using the same parametric population models applied to real GWTC-4 data. It reports that for GWTC-4-like catalogs, confident cutoff identification is not guaranteed but the strongest recovered constraints are compatible with existing results when a cutoff is present, while spurious detections are unlikely. For O4-scale catalogs the cutoff-mass uncertainty shrinks by at least 20 percent, yet only a lower bound on the 12C(α,γ)16O S-factor at 300 keV is obtained. A nonparametric PixelPop analysis of GWTC-4 is also presented and shown to be consistent with the parametric results under a cutoff model.
Significance. If the central simulation-based validation holds, the work supplies a necessary controlled test of the inference pipeline used to claim mass cutoffs from GWTC-4. The combination of forward simulations, likelihood-based predictive checks, and a nonparametric cross-check provides a concrete assessment of detection prospects and of the robustness of current parametric conclusions. The quantitative statements on uncertainty reduction and on the limited constraining power for nuclear rates are directly useful for planning future observing runs and for interpreting any future cutoff detection.
major comments (1)
- [§4.2] §4.2 and the simulation setup: the claim that the parametric models accurately capture the true underlying mass distribution and selection effects is load-bearing for the recovery tests. If the real population contains features (e.g., additional peaks or a smoother taper) not spanned by the chosen parametric family, the reported false-positive and false-negative rates may not translate to real data. A quantitative sensitivity test to model misspecification would strengthen the result.
minor comments (3)
- [Abstract] The abstract states that the best constraints are 'compatible with results from GWTC-4'; it would help to quote the numerical GWTC-4 cutoff posterior (median and 90 % interval) for direct comparison.
- [Figures 2-4] Figure captions for the simulated catalogs should explicitly state the injected cutoff location and the number of events in each realization.
- [§5.3] Notation for the S-factor (S_300) is introduced without a reference to the nuclear-physics definition; a brief parenthetical or footnote would remove ambiguity.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the positive assessment. We address the single major comment below and will incorporate the requested sensitivity test in the revised version.
read point-by-point responses
-
Referee: [§4.2] §4.2 and the simulation setup: the claim that the parametric models accurately capture the true underlying mass distribution and selection effects is load-bearing for the recovery tests. If the real population contains features (e.g., additional peaks or a smoother taper) not spanned by the chosen parametric family, the reported false-positive and false-negative rates may not translate to real data. A quantitative sensitivity test to model misspecification would strengthen the result.
Authors: We agree that model misspecification is a valid concern for translating the reported rates to real data. Our simulations are intentionally performed within the same parametric family used for the GWTC-4 analysis to test the inference pipeline under its own assumptions, and the nonparametric PixelPop results on both simulated and real data already provide a consistency check. To directly quantify the impact of misspecification, we will add a new test in the revised manuscript: we will generate additional mock catalogs using an alternative mass distribution (a smoother logistic taper plus an optional secondary peak outside the parametric family) and recover them with the standard parametric model, reporting the resulting changes in false-positive and false-negative rates for cutoff detection. revision: yes
Circularity Check
No significant circularity
full rationale
The paper generates simulated catalogs independently of the real GWTC-4 data, injects known cutoffs or their absence, and recovers them using parametric models; it then applies both parametric and nonparametric (PixelPop) models to the actual catalog for comparison. These steps are forward simulations and separate data analyses rather than any fitted parameter being relabeled as a prediction or any result being defined in terms of itself. No self-citation chain is load-bearing for the central claims about detection prospects, and the nonparametric check serves as an internal consistency test without reducing the conclusions to the inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The simulated catalogs faithfully reproduce the selection biases and parameter-estimation uncertainties of the LIGO-Virgo-KAGRA network.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We simulate gravitational-wave catalogs with and without such a mass cutoff, then test whether its presence or absence is correctly inferred with parametric population models.
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We also analyze GWTC-4 with the nonparametric PixelPop population model
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 3 Pith papers
-
Posterior Predictive Checks for Gravitational-wave Populations: Limitations and Improvements
Maximum-likelihood-based posterior predictive checks detect model misspecification better than event-level versions for uncertain spin tilts, but current detector sensitivity limits their power; the Gaussian Component...
-
Gravitational-wave constraints on $H_0$ are robust to (putative) redshift evolution in the binary black hole mass spectrum at current sensitivity
Spectral-siren H0 constraints from GWTC-4.0 binary black holes remain robust when the mass spectrum is permitted to evolve with redshift at current detector sensitivity.
-
Emergent structure in the binary black hole mass distribution and implications for population-based cosmology
B-spline agnostic reconstruction of binary black hole masses from GWTC-4.0 reveals multiple features and a logarithmic hierarchy that impacts Hubble constant measurements, with a low-mass subpopulation isolation metho...
Reference graph
Works this paper leans on
-
[1]
need to be revised. The distributions do not de- cline sharply at low and high masses, reflecting the ag- nosticPixelPopprior [75], unlike the assumption of strict minima and maxima in the parametric models. In- stead, in Figure 1 we show the 99%, 99.9%, and 99.99% percentile masses, computed from the marginal cumula- tive distributions form 1 andm 2 from...
-
[2]
but, unlike common approaches in GW population studies, we compare observed and predicted sources self consistently (see Ref. [71] for a discussion). For each ob- served event, we select the parameter estimate with the highest likelihood1 among its posterior samples. We com- pare the observed highest-likelihood secondary masses to predictions from theSing...
-
[3]
Over the broader range 30< m 2/M⊙ <150, the GWTC-4 median is typically shallower than results from catalogs of all three simulated populations. From the Single Power Law + 2 Peaks + Cutoffcatalogs, PixelPopidentifies that them 2 slope over 30–60M ⊙ is steeper than over 10–20M ⊙ with posterior probabilities (≳92%) consistently more similar to the GWTC-4 re...
work page 2000
-
[4]
J. Aasiet al., Class. Quant. Grav.32, 074001 (2015), arXiv:1411.4547 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[5]
E. Capoteet al., Phys. Rev. D111, 062002 (2025), arXiv:2411.14607 [gr-qc]
-
[6]
(LIGO Instrument Science Collaboration) 2025 Class
S. Soniet al., Class. Quant. Grav.42, 085016 (2025), arXiv:2409.02831 [astro-ph.IM]
-
[7]
Advanced Virgo: a 2nd generation interferometric gravitational wave detector
F. Acerneseet al., Class. Quant. Grav.32, 024001 (2015), arXiv:1408.3978 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[8]
Akutsuet al.(KAGRA), PTEP2021, 05A101 (2021), arXiv:2005.05574 [physics.ins-det]
T. Akutsuet al., PTEP2021, 05A101 (2021), arXiv:2005.05574 [physics.ins-det]
-
[9]
B. P. Abbottet al., Phys. Rev. X9, 031040 (2019), arXiv:1811.12907 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[10]
R. Abbottet al., Phys. Rev. X11, 021053 (2021), arXiv:2010.14527 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[11]
R. Abbottet al., Phys. Rev. D109, 022001 (2024), arXiv:2108.01045 [gr-qc]
work page internal anchor Pith review arXiv 2024
-
[12]
R. Abbottet al., Phys. Rev. X13, 041039 (2023), arXiv:2111.03606 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[13]
A. G. Abacet al., arXiv:2508.18082 [gr-qc] (2025)
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[14]
M. Mapelli, Front. Astron. Space Sci.7, 38 (2020), arXiv:2105.12455 [astro-ph.HE]
-
[15]
Mapelli, inHandbook of Gravitational Wave Astron- omy(Springer, Singapore, 2020) pp
M. Mapelli, inHandbook of Gravitational Wave Astron- omy(Springer, Singapore, 2020) pp. 1–65
work page 2020
-
[16]
I. Mandel and F. S. Broekgaarden, Living Rev. Rel.25, 1 (2022), arXiv:2107.14239 [astro-ph.HE]
- [17]
-
[18]
R. Abbottet al., Astrophys. J. Lett.913, L7 (2021), arXiv:2010.14533 [astro-ph.HE]
-
[19]
The population of merging compact binaries inferred using gravitational waves through GWTC-3
R. Abbottet al., Phys. Rev. X13, 011048 (2023), arXiv:2111.03634 [astro-ph.HE]
work page internal anchor Pith review arXiv 2023
-
[20]
A. G. Abacet al., arXiv:2508.18083 [astro-ph.HE] (2025)
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[21]
W. A. Fowler and F. Hoyle, Astrophys. J. Suppl.9, 201 (1964)
work page 1964
- [22]
- [23]
-
[24]
G. S. Fraley, Astrophysics and Space Science2, 96 (1968)
work page 1968
-
[25]
S. E. Woosley, Astrophys. J.836, 244 (2017), arXiv:1608.08939 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[26]
The Effect of Pair-Instability Mass Loss on Black Hole Mergers
K. Belczynskiet al., Astron. Astrophys.594, A97 (2016), arXiv:1607.03116 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[27]
The Impact of Pair-instability Mass Loss on the Binary Black Hole Mass Distribution,
S. Stevenson, M. Sampson, J. Powell, A. Vigna-G´ omez, C. J. Neijssel, D. Sz´ ecsi, and I. Mandel, The Astrophys- ical Journal882, 121 (2019), arXiv:1904.02821 [astro- ph.HE]
- [28]
- [29]
- [30]
- [31]
-
[32]
F. Gabrielliet al., Mon. Not. Roy. Astron. Soc.534, 151 (2024), arXiv:2408.16823 [astro-ph.HE]
- [33]
-
[34]
Impact of the Rotation and Compactness of Progenitors on the Mass of Black Holes,
M. Mapelli, M. Spera, E. Montanari, M. Limongi, A. Chieffi, N. Giacobbo, A. Bressan, and Y. Bouffanais, Astrophys. J.888, 76 (2020), arXiv:1909.01371 [astro- ph.HE]
- [35]
-
[36]
The impact of stellar rotation on the black hole mass-gap from pair- instability supernovae,
P. Marchant and T. Moriya, Astron. Astrophys.640, L18 (2020), arXiv:2007.06220 [astro-ph.HE]
-
[37]
G. Costa, A. Bressan, M. Mapelli, P. Marigo, G. Iorio, and M. Spera, Mon. Not. Roy. Astron. Soc.501, 4514 (2021), arXiv:2010.02242 [astro-ph.SR]
-
[38]
R. J. deBoeret al., Rev. Mod. Phys.89, 035007 (2017), arXiv:1709.03144 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[39]
The Low Detection Rate of Pair Instability Supernovae and the Effect of the Core Carbon Fraction
K. Takahashi, Astrophys. J.863, 153 (2018), arXiv:1807.05373 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [40]
-
[41]
Measuring the binary black hole mass spectrum with an astrophysically motivated parameterization
C. Talbot and E. Thrane, Astrophys. J.856, 173 (2018), arXiv:1801.02699 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[42]
F. Antonini, M. Gieles, F. Dosopoulou, and D. Chat- topadhyay, Mon. Not. Roy. Astron. Soc.522, 466 (2023), arXiv:2208.01081 [astro-ph.HE]
- [43]
- [44]
-
[45]
2025, arXiv e-prints, arXiv:2511.22093, doi: 10.48550/arXiv.2511.22093
O. Sridhar, A. Ray, and V. Kalogera, arXiv:2511.22093 [astro-ph.HE] (2025)
-
[46]
Y.-Z. Wang, S.-P. Tang, Y.-F. Liang, M.-Z. Han, X. Li, Z.-P. Jin, Y.-Z. Fan, and D.-M. Wei, Astrophys. J.913, 42 (2021), arXiv:2104.02566 [astro-ph.HE]
- [47]
- [48]
-
[49]
D. Gerosa and M. Fishbach, Nature Astron.5, 749 (2021), arXiv:2105.03439 [astro-ph.HE]
-
[50]
C. Karathanasis, S. Mukherjee, and S. Mastrogio- vanni, Mon. Not. Roy. Astron. Soc.523, 4539 (2023), arXiv:2204.13495 [astro-ph.CO]
-
[51]
Y.-Z. Wang, Y.-J. Li, J. S. Vink, Y.-Z. Fan, S.-P. Tang, Y. Qin, and D.-M. Wei, Astrophys. J. Lett.941, L39 (2022), arXiv:2208.11871 [astro-ph.HE]
- [52]
- [53]
-
[54]
F. Antonini, I. M. Romero-Shaw, and T. Callister, Phys. Rev. Lett.134, 011401 (2025), arXiv:2406.19044 [astro- ph.HE]
-
[55]
F. Antonini, T. Callister, F. Dosopoulou, I. M. Romero- Shaw, and D. Chattopadhyay, Phys. Rev. D112, 063040 (2025), arXiv:2506.09154 [astro-ph.HE]
-
[56]
I. Maga˜ na Hernandez and A. Palmese, arXiv:2508.19208 [astro-ph.HE] (2025)
-
[57]
Evidence of the pair instability gap from black hole masses
H. Tonget al., Nature 10.1038/s41586-026-10359-0 (2026), arXiv:2509.04151 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1038/s41586-026-10359-0 2026
-
[58]
N. Guttman, E. Payne, P. D. Lasky, and E. Thrane, arXiv:2509.09876 [astro-ph.HE] (2025)
-
[59]
A. G. Abacet al., Astrophys. J. Lett.993, L25 (2025), arXiv:2507.08219 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[60]
Mandel, The Astrophysical Journal996, L4 (2026), arXiv:2509.05885 [astro-ph.HE]
I. Mandel, Astrophys. J. Lett.996, L4 (2026), arXiv:2509.05885 [astro-ph.HE]
- [61]
-
[62]
Gravitational-wave constraints on the pair-instability mass gap and nuclear burning in massive stars
F. Antonini, I. Romero-Shaw, T. Callister, F. Dosopoulou, D. Chattopadhyay, M. Gieles, and M. Mapelli, arXiv:2509.04637 [astro-ph.HE] (2025)
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [63]
- [64]
-
[65]
S. Banagiri, E. Thrane, and P. D. Lasky, arXiv:2509.15646 [astro-ph.HE] (2025)
- [66]
-
[67]
Signatures of a subpopulation of hierarchical mergers in the GWTC-4 gravitational-wave dataset
C. Plunkett, T. Callister, M. Zevin, and S. Vitale, arXiv:2601.07908 [gr-qc] (2026)
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[68]
A. M. Farah, A. Vijaykumar, and M. Fishbach, arXiv:2601.03456 [astro-ph.HE] (2026)
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[69]
A. Vijaykumar, A. M. Farah, and M. Fishbach, Astro- phys. J. Lett.999, L30 (2026), arXiv:2601.03457 [astro- ph.HE]
-
[70]
A new group of low-spin $50-70M_\odot$ Black Holes and the high pair-instability mass cutoff
Y.-Z. Wang, Y.-J. Li, S.-J. Gao, S.-P. Tang, and Y.-Z. Fan, arXiv:2510.22698 [astro-ph.HE] (2025)
work page internal anchor Pith review arXiv 2025
- [71]
- [72]
-
[73]
A. Toubiana, M. L. Katz, and J. R. Gair, Mon. Not. Roy. Astron. Soc.524, 5844 (2023), arXiv:2305.08909 [gr-qc]
- [74]
-
[75]
S. Biscoveanu, Astrophys. J.996, 95 (2026), arXiv:2502.04278 [astro-ph.HE]
- [76]
-
[77]
A. G. Abacet al., arXiv:2508.18079 [gr-qc] (2025)
work page internal anchor Pith review arXiv 2025
-
[78]
J. Heinzel, M. Mould, S. ´Alvarez-L´ opez, and S. Vitale, Phys. Rev. D111, 063043 (2025), arXiv:2406.16813 [astro-ph.HE]
-
[79]
J. Heinzel, M. Mould, and S. Vitale, Phys. Rev. D111, L061305 (2025), arXiv:2406.16844 [astro-ph.HE]
-
[80]
S. Alvarez-Lopez, J. Heinzel, M. Mould, and S. Vitale, arXiv:2506.20731 [astro-ph.HE] (2025)
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.