Constraints on primordial black holes from the first part of LIGO-Virgo-KAGRA fourth observing run
Pith reviewed 2026-05-20 17:09 UTC · model grok-4.3
The pith
LIGO-Virgo-KAGRA O4a data yields the strongest constraints on primordial black hole abundance between 0.6 and 100 solar masses.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We analyze PBH populations using state-of-the-art modeling of PBH binaries, deriving the strongest bounds on PBH abundance in the 0.6-100 M_⊙ range from LIGO-Virgo-KAGRA O4a data and demonstrating sensitivity in the 10^{-4}-10^4 M_⊙ range. The constraints are dominated by resolvable PBH mergers, while the associated gravitational wave background provides complementary but weaker limits. Allowing PBHs to account for a subset of the cataloged events slightly relaxes these bounds. However, a joint fit with astrophysical black holes shows no compelling evidence for a PBH contribution.
What carries the argument
State-of-the-art modeling of PBH binary formation, evolution, and merger rates that predicts the expected number of resolvable events and the stochastic gravitational wave background.
If this is right
- Resolvable PBH mergers dominate the constraints on abundance in the 0.6-100 solar mass range.
- The stochastic gravitational wave background supplies complementary but weaker limits.
- Attributing some cataloged events to PBHs only mildly relaxes the upper bounds.
- A joint fit with astrophysical black holes finds no compelling evidence for a PBH contribution.
Where Pith is reading between the lines
- The full O4 dataset or later runs could produce even tighter abundance limits using the same approach.
- These gravitational-wave bounds complement independent probes such as microlensing surveys.
- Refinements to binary evolution assumptions could shift the numerical values while preserving the overall method.
Load-bearing premise
The modeling of PBH binary formation, evolution, and merger rates accurately predicts the expected number of resolvable events and the stochastic background without large systematic uncertainties from early-universe conditions or binary dynamics.
What would settle it
Detecting substantially more black hole mergers in the 0.6-100 solar mass range than allowed by the maximum permitted PBH abundance, or measuring a stochastic gravitational wave background exceeding the predicted level, would challenge the derived constraints.
Figures
read the original abstract
We analyze PBH populations using state-of-the-art modeling of PBH binaries, deriving the strongest bounds on PBH abundance in the $0.6-100 M_\odot$ range from LIGO-Virgo-KAGRA O4a data and demonstrating sensitivity in the $10^{-4}-10^4 M_\odot$ range. The constraints are dominated by resolvable PBH mergers, while the associated gravitational wave background provides complementary but weaker limits. Allowing PBHs to account for a subset of the cataloged events slightly relaxes these bounds. However, a joint fit with astrophysical black holes shows no compelling evidence for a PBH contribution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes PBH populations using state-of-the-art modeling of PBH binaries and derives upper limits on the PBH abundance f_PBH from LIGO-Virgo-KAGRA O4a data. It claims the strongest bounds to date in the 0.6-100 M_⊙ range, with constraints dominated by resolvable mergers and weaker complementary limits from the stochastic gravitational-wave background. The analysis also explores allowing PBHs to account for a subset of cataloged events (which relaxes bounds slightly) and performs a joint fit with astrophysical black holes, finding no compelling evidence for a PBH contribution.
Significance. If the central modeling assumptions hold, the work supplies timely new constraints on stellar-mass PBH dark matter using the latest O4a dataset, improving on previous limits and extending sensitivity to a wider mass range. The explicit separation of resolvable-merger and background contributions, together with the joint-fit exercise, strengthens the interpretability of the results.
major comments (1)
- [§4] §4 (PBH binary formation and merger-rate modeling): the expected number of resolvable events is computed from a specific choice of initial separation distribution, clustering prescription, and dynamical hardening. No systematic variation is shown for plausible alternatives (e.g., non-Gaussian primordial fluctuations or altered small-scale power spectra), which the text acknowledges can change the merger rate by a factor of a few. Because the headline bounds are stated to be dominated by resolvable mergers, this unquantified uncertainty directly affects the quoted f_PBH limits and must be addressed before the “strongest bounds” claim can be considered robust.
minor comments (1)
- [Abstract and §5] The abstract and §5 would benefit from a one-sentence statement of the precise early-universe initial conditions adopted for the baseline rate calculation.
Simulated Author's Rebuttal
We are grateful to the referee for the positive assessment of our work and the constructive major comment. Below we provide a detailed response and indicate the changes we will implement in the revised manuscript.
read point-by-point responses
-
Referee: [§4] §4 (PBH binary formation and merger-rate modeling): the expected number of resolvable events is computed from a specific choice of initial separation distribution, clustering prescription, and dynamical hardening. No systematic variation is shown for plausible alternatives (e.g., non-Gaussian primordial fluctuations or altered small-scale power spectra), which the text acknowledges can change the merger rate by a factor of a few. Because the headline bounds are stated to be dominated by resolvable mergers, this unquantified uncertainty directly affects the quoted f_PBH limits and must be addressed before the “strongest bounds” claim can be considered robust.
Authors: We appreciate the referee's observation regarding the modeling assumptions in §4. Our analysis utilizes the most up-to-date modeling of PBH binary formation and evolution available in the literature. We explicitly acknowledge in the text that alternative scenarios, including non-Gaussian primordial fluctuations or modified small-scale power spectra, could lead to variations in the merger rate by a factor of a few. While a comprehensive systematic study of all such alternatives is beyond the scope of this paper, we agree that this uncertainty should be better highlighted. In the revised version, we will expand the discussion in §4 to include a qualitative assessment of how these variations would propagate to the f_PBH limits, thereby qualifying the 'strongest bounds' claim appropriately. This addresses the concern while maintaining the focus on the new O4a data constraints. revision: partial
Circularity Check
No circularity: bounds derived from external LIGO data comparison
full rationale
The paper models PBH binary formation and merger rates using established prescriptions, then compares the predicted event counts and stochastic background directly to independent LIGO-Virgo-KAGRA O4a observations to extract upper limits on f_PBH. This is a standard forward-modeling exercise against external counts; the reported bounds do not reduce by construction to any parameter fitted inside the O4a dataset itself, nor to a self-citation chain that replaces the data comparison. The derivation chain therefore remains self-contained against the external benchmark.
Axiom & Free-Parameter Ledger
free parameters (1)
- PBH mass function parameters
axioms (2)
- standard math General relativity accurately describes gravitational wave generation, propagation, and detection by LIGO-Virgo-KAGRA.
- domain assumption The adopted state-of-the-art PBH binary formation and merger rate model correctly represents the expected event rate for a given abundance.
Reference graph
Works this paper leans on
-
[1]
In this case, the constraint onf PBH is derived from the non-observation of any excess of events
In theastrophysics onlycase, all the observed events originate from ABH binaries and no PBH are assumed to have been observed. In this case, the constraint onf PBH is derived from the non-observation of any excess of events. The 95% CL upper bound onf PBH corresponds to the region whereN(Λ)<3
-
[2]
In theastrophysically agnosticcase, any of the ob- served events could be primordial. We therefore set limits without modeling the ABH population, using two complementary procedures: (i) taking random subsets of the catalog and assuming that they are PBHs, as done in Refs. [10, 15]); and (ii) using a sub- set–marginalized likelihood (SML), which smoothly ...
-
[3]
In theastrophysically informedcase, we model all ob- served events from a merger rate that combines con- tributions from both PBHs and ABHs. For BH bina- ries formed through astrophysical channels, a reliable estimate of the merger rate based on first principles is missing, which is why we adopt the phenomeno- logical power-law + 2 peaks (PL+2P) parameter...
-
[4]
E. Capote et al., Phys. Rev. D111, 062002 (2025), 2411.14607
-
[5]
B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. X9, 031040 (2019), 1811.12907
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[6]
R. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. X 11, 021053 (2021), 2010.14527
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[7]
R. Abbott et al. (KAGRA, VIRGO, LIGO Scientific), Phys. Rev. X13, 041039 (2023), 2111.03606
work page internal anchor Pith review Pith/arXiv arXiv 2023
- [8]
-
[9]
S. E. Woosley, A. Heger, and T. A. Weaver, Rev. Mod. Phys.74, 1015 (2002)
work page 2002
-
[10]
B. J. Carr and S. W. Hawking, Mon. Not. Roy. Astron. Soc.168, 399 (1974)
work page 1974
-
[11]
G. H¨ utsi, M. Raidal, V. Vaskonen, and H. Veerm¨ ae, JCAP03, 068 (2021), 2012.02786
- [12]
- [13]
-
[14]
G. Franciolini, V. Baibhav, V. De Luca, K. K. Y. Ng, K. W. K. Wong, E. Berti, P. Pani, A. Riotto, and S. Vi- tale, Phys. Rev. D105, 083526 (2022), 2105.03349
-
[15]
V. De Luca, G. Franciolini, P. Pani, and A. Riotto, JCAP05, 003 (2021), 2102.03809
-
[16]
M. Andr´ es-Carcasona, A. J. Iovino, V. Vaskonen, H. Veerm¨ ae, M. Mart´ ınez, O. Pujol` as, and L. M. Mir, Phys. Rev. D110, 023040 (2024), 2405.05732
-
[17]
O. Pujolas, V. Vaskonen, and H. Veerm¨ ae, Phys. Rev. D104, 083521 (2021), 2107.03379
-
[18]
S. Mukherjee and J. Silk, Mon. Not. Roy. Astron. Soc. 506, 3977 (2021), 2105.11139
-
[19]
Gravitational Waves from Primordial Black Hole Mergers
M. Raidal, V. Vaskonen, and H. Veerm¨ ae, JCAP09, 037 (2017), 1707.01480
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[20]
T. Boybeyi, S. Clesse, S. Kuroyanagi, and M. Sakellar- iadou, Phys. Rev. D112, 023551 (2025), 2412.18318
- [21]
- [22]
-
[23]
A. Romero-Rodr´ ıguez and S. Kuroyanagi (2024), 2407.00205
-
[24]
Planck 2018 results. VI. Cosmological parameters
N. Aghanim et al. (Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], 1807.06209
work page internal anchor Pith review Pith/arXiv arXiv 2020
-
[25]
F. Crescimbeni, V. Desjacques, G. Franciolini, A. Ian- niccari, A. J. Iovino, G. Perna, D. Perrone, A. Riotto, and H. Veerm¨ ae (2025), 2502.01617
-
[26]
B. Carr, M. Raidal, T. Tenkanen, V. Vaskonen, and H. Veerm¨ ae, Phys. Rev. D96, 023514 (2017), 1705.05567
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [27]
- [28]
-
[29]
See Supplemental Material for details
- [30]
-
[31]
C. T. Byrnes, M. Hindmarsh, S. Young, and M. R. S. Hawkins, JCAP08, 041 (2018), 1801.06138
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[32]
G. Franciolini, I. Musco, P. Pani, and A. Urbano, Phys. Rev. D106, 123526 (2022), 2209.05959
-
[33]
A. Escriv` a, E. Bagui, and S. Clesse, JCAP05, 004 (2023), 2209.06196
- [34]
-
[35]
G. Ferrante, G. Franciolini, A. Iovino, Junior., and 7 A. Urbano, Phys. Rev. D107, 043520 (2023), 2211.01728
-
[36]
E. Thrane and C. Talbot, Publ. Astron. Soc. Austral. 36, e010 (2019), [Erratum: Publ.Astron.Soc.Austral. 37, e036 (2020)], 1809.02293
-
[37]
Extracting distribution parameters from multiple uncertain observations with selection biases
I. Mandel, W. M. Farr, and J. R. Gair, Mon. Not. Roy. Astron. Soc.486, 1086 (2019), 1809.02063
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[38]
R. Abbott et al. (LIGO Scientific, Virgo), Astrophys. J. Lett.913, L7 (2021), 2010.14533
-
[39]
S. Mastrogiovanni, G. Pierra, S. Perri` es, D. Laghi, G. Caneva Santoro, A. Ghosh, R. Gray, C. Karathana- sis, and K. Leyde, Astron. Astrophys.682, A167 (2024), 2305.17973
-
[40]
L. S. Finn and D. F. Chernoff, Phys. Rev. D47, 2198 (1993), gr-qc/9301003
work page internal anchor Pith review Pith/arXiv arXiv 1993
-
[41]
Multiband gravitational-wave event rates and stellar physics
D. Gerosa, S. Ma, K. W. K. Wong, E. Berti, R. O’Shaughnessy, Y. Chen, and K. Belczynski, Phys. Rev. D99, 103004 (2019), 1902.00021
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[42]
A. G. Abac et al. (LIGO Scientific, KAGRA, VIRGO), Astrophys. J. Lett.995, L18 (2025), 2508.18080
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [43]
-
[44]
GW231123: A Possible Primordial Black Hole Origin
V. De Luca, G. Franciolini, and A. Riotto (2025), 2508.09965
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [45]
- [46]
- [47]
-
[48]
A. G. Abac et al. (LIGO Scientific, VIRGO, KAGRA) (2025), 2509.04348
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [49]
-
[50]
R. Abbott et al. (LIGO Scientific, VIRGO, KAGRA), Phys. Rev. Lett.129, 061104 (2022), 2109.12197
-
[51]
R. Abbott et al. (LVK), Mon. Not. Roy. Astron. Soc. 524, 5984 (2023), [Erratum: Mon.Not.Roy.Astron.Soc. 526, 6234 (2023)], 2212.01477
- [52]
- [53]
-
[54]
M. Andres-Carcasona, A. Menendez-Vazquez, M. Mar- tinez, and L. M. Mir, Phys. Rev. D107, 082003 (2023), 2212.02829
-
[55]
A. G. Abac et al. (LIGO Scientific, VIRGO, KAGRA) (2026), 2605.05444
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [56]
- [57]
- [58]
- [59]
-
[60]
M. Andr´ es-Carcasona, O. J. Piccinni, M. Mart´ ınez, and L. M. Mir, Phys. Rev. D111, 043019 (2025), 2411.04498
-
[61]
M. Andr´ es-Carcasona, O. J. Piccinni, M. Mart´ ınez, and L.-M. Mir, PoSEPS-HEP2023, 067 (2024)
work page 2024
-
[62]
G. Alestas, G. Morras, T. S. Yamamoto, J. Garcia- Bellido, S. Kuroyanagi, and S. Nesseris, Phys. Rev. D 109, 123516 (2024), 2401.02314
- [63]
- [64]
-
[65]
P. Astone, A. Colla, S. D’Antonio, S. Frasca, and C. Palomba, Phys. Rev. D90, 042002 (2014), 1407.8333
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[66]
M. Maggiore,Gravitational Waves. Vol. 1: Theory and Experiments(Oxford University Press, 2007), ISBN 978-0-19-171766-6, 978-0-19-852074-0
work page 2007
-
[67]
B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. Lett.118, 121101 (2017), [Erratum: Phys.Rev.Lett. 119, 029901 (2017)], 1612.02029
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [68]
-
[69]
Abbottet al.(KAGRA, Virgo, LIGO Scientific), Phys
R. Abbott et al. (KAGRA, Virgo, LIGO Scientific), Phys. Rev. D104, 022004 (2021), 2101.12130
- [70]
-
[71]
R. Abbott et al. (KAGRA, Virgo, LIGO Scientific), Phys. Rev. D104, 022005 (2021), 2103.08520
-
[72]
R. Abbott et al. (KAGRA, Virgo, LIGO Scientific), Phys. Rev. D105, 122001 (2022), 2110.09834
-
[73]
A. Romero-Rodriguez, M. Martinez, O. Pujol` as, M. Sakellariadou, and V. Vaskonen, Phys. Rev. Lett. 128, 051301 (2022), 2107.11660
- [74]
-
[75]
S. Wang, Y.-F. Wang, Q.-G. Huang, and T. G. F. Li, Phys. Rev. Lett.120, 191102 (2018), 1610.08725
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [76]
-
[77]
J.-Z. Zhou, Y.-T. Kuang, Z. Chang, and H. L¨ u, Astro- phys. J.979, 178 (2025), 2410.10111
- [78]
-
[79]
B. Allen and J. D. Romano, Phys. Rev. D59, 102001 (1999), gr-qc/9710117
work page internal anchor Pith review Pith/arXiv arXiv 1999
-
[80]
J. D. Romano and N. J. Cornish, Living Rev. Rel.20, 2 (2017), 1608.06889
work page internal anchor Pith review Pith/arXiv arXiv 2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.