pith. sign in

arxiv: 2607.01121 · v1 · pith:TB74WW4Znew · submitted 2026-07-01 · 🌌 astro-ph.HE · astro-ph.CO· astro-ph.GA· gr-qc

Smoking-gun evidence for hierarchical black-hole mergers

Yin-Jie Li , Yuan-Zhu Wang , Shao-Peng Tang , Yi-Zhong Fan This is my paper

Pith reviewed 2026-07-02 07:14 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.COastro-ph.GAgr-qc
keywords black hole mergershierarchical mergersgravitational wavesmass functionspin distributionstellar collapsenuclear astrophysics
0
0 comments X

The pith

The masses of high-spin black holes match the expected remnants from low-spin black holes peak by peak, showing they form through hierarchical mergers.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

A flexible mixture model is applied to 259 binary black hole merger events to separate those with low and high spins. The mass function for the high-spin group closely follows the distribution expected for remnants produced by merging the low-spin black holes, with peaks aligning up to about 80 solar masses. This similarity reaches a coefficient of 0.95, which is what hierarchical mergers would produce but other explanations would need to arrange by chance. The analysis also uses the mass cutoff of the low-spin group to calculate the rate of a nuclear fusion reaction that matches theoretical predictions. This accounts for all observed black holes through stellar birth and repeated mergers.

Core claim

The mass function of the high-spin subpopulation traces, peak by peak, the predicted remnant-mass distribution of the low-spin, stellar-collapse-origin subpopulation up to ∼80M⊙. This morphological match, quantified by a Bhattacharyya coefficient as high as ∼0.95, is naturally expected if the high-spin black holes are themselves the products of earlier mergers, whereas any alternative scenario would require fine-tuning, thereby providing smoking-gun evidence for hierarchical mergers. In addition, the sharp upper-mass cutoff of the low-spin subpopulation at m_max,1=54.2+7.7−7.2M⊙ yields an astrophysical S-factor of S300=151+30−26 keV b for the 12C(α,γ)16O reaction.

What carries the argument

Flexible mixture population model that separates low-spin and high-spin subpopulations and generates a predicted remnant-mass distribution from the low-spin data alone.

If this is right

  • The high-spin black holes result from previous mergers of lower-mass black holes.
  • The full population of stellar-mass black holes arises from stellar collapse followed by dynamical hierarchical assembly.
  • No contribution from primordial black holes is needed to explain the observations.
  • The upper mass limit of the low-spin subpopulation constrains the astrophysical S-factor for the carbon-alpha-oxygen reaction to 151 keV b.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Observations of more events could reveal the expected increase in spin with each generation of mergers.
  • This approach might distinguish between different dense environments where mergers occur.
  • Similar mass-function tracing could be sought in other gravitational-wave populations.
  • Future data releases will test if the match persists at higher masses.

Load-bearing premise

The remnant-mass distribution predicted from the low-spin subpopulation is independent of the high-spin data and free of post-hoc tuning in the flexible mixture model.

What would settle it

A larger sample of mergers in which the mass peaks of the high-spin subpopulation fail to align with the remnant predictions from the low-spin subpopulation.

Figures

Figures reproduced from arXiv: 2607.01121 by Shao-Peng Tang, Yin-Jie Li, Yi-Zhong Fan, Yuan-Zhu Wang.

Figure 1
Figure 1. Figure 1: The reconstructed mass distributions of the black hole groups. The solid curves are [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The Bhattacharyya coefficient as a function of the lower edge of the mass window [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The inferred S300keV. The evaluated S factor of 12C(α, γ) 16O reaction at 300 keV assuming that mmax,1 represents the lower edge of the PISN mass gap. The value of S300keV recommended by ref. 41 (1σ region including the model uncertainty) is also plotted for comparison. retention fraction) of the remnant as fmerger ≡ r2 (1−r2)/2 = 23+17 −11% and f(m)merger ∼ 30% at the mass ≳ 30M⊙ (see Extended Data [PITH… view at source ↗
Figure 6
Figure 6. Figure 6: Second, we re-analyze the mass distributions of the subpopulations in two additional [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 6
Figure 6. Figure 6: The fiducial population model We adopt a flexible mixture model describing m1−m2−χeff −χp distribution as our primary analysis, which allows (two) subpopulations to have different 4-D parameter distributions. A 5-D mixture model (including different rate evolution model) is also tested, and the evidence that the high-spin/high-mass subpopulation evolves faster (so called ‘the heavier the faster’ 51) is wea… view at source ↗
read the original abstract

How stellar-mass black holes grow after their birth is a central open question in astrophysics. Gravitational-wave observations have revealed a subpopulation of coalescing black holes with both high masses and high spins, but whether these properties arise from hierarchical mergers in dense stellar environments or from accretion onto isolated black holes has remained unresolved. Here, using a flexible mixture population model applied to the 259 binary black hole mergers in GWTC-5, we show that the mass function of the high-spin subpopulation traces, peak by peak, the predicted remnant-mass distribution of the low-spin, stellar-collapse-origin subpopulation up to $\sim80\,M_\odot$. This morphological match, quantified by a Bhattacharyya coefficient as high as $\sim0.95$, is naturally expected if the high-spin black holes are themselves the products of earlier mergers, whereas any alternative scenario would require fine-tuning, thereby providing smoking-gun evidence for hierarchical mergers. In addition, the sharp upper-mass cutoff of the low-spin subpopulation at $m_{\rm max,1}=54.2^{+7.7}_{-7.2}\,M_\odot$ yields an astrophysical $S$-factor of $S_{300}=151^{+30}_{-26}$~keV~b (68\% credible interval) for the $^{12}{\rm C}(\alpha,\gamma)^{16}{\rm O}$ reaction, in agreement with the benchmark theoretical value. These results establish that the entire observed black-hole population can be accounted for by stellar collapse followed by dynamical hierarchical assembly, without invoking primordial black holes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper applies a flexible mixture population model to the 259 binary black hole events in GWTC-5 and reports that the mass function of the high-spin subpopulation morphologically matches (Bhattacharyya coefficient up to ~0.95) the remnant-mass distribution predicted from the low-spin subpopulation up to ~80 M⊙. This match is presented as smoking-gun evidence for hierarchical mergers, with the low-spin upper-mass cutoff at 54.2^{+7.7}_{-7.2} M⊙ additionally used to infer an astrophysical S-factor of 151^{+30}_{-26} keV b for the ^{12}C(α,γ)^{16}O reaction.

Significance. If the reported morphological match is shown to be independent of the high-spin data, the result would be significant: it would provide a direct, falsifiable test distinguishing hierarchical assembly from isolated accretion or primordial channels using only public catalog data, while simultaneously yielding a nuclear-physics constraint from the low-spin cutoff. The analysis uses the full GWTC-5 sample and reports a quantitative coefficient, both of which strengthen the claim if the independence issue is resolved.

major comments (3)
  1. [Abstract] Abstract (mixture population model paragraph): the central claim requires that the low-spin mass function used to generate the predicted remnant distribution is estimated independently of the high-spin events. The manuscript describes a single joint fit in which mixture weights, component mass functions, and spin subpopulations are inferred simultaneously; no explicit statement is given that the low-spin mass function is obtained from a conditional posterior restricted to high-probability low-spin events or from a two-stage procedure that holds the high-spin data out. This directly affects whether the ~0.95 Bhattacharyya coefficient constitutes independent evidence or partly reflects internal model consistency.
  2. [Abstract] Abstract (mixture population model paragraph): the flexible mixture model is stated to be applied to all 259 events, yet no details are provided on the functional forms of the component mass functions, the prior on mixture weights, the treatment of measurement uncertainties in the joint posterior, or robustness checks against alternative subpopulation definitions. These choices are load-bearing for the reported peak-by-peak match and the claimed absence of fine-tuning in alternative scenarios.
  3. [Abstract] Abstract (S-factor paragraph): the sharp upper-mass cutoff m_max,1 = 54.2^{+7.7}_{-7.2} M⊙ is used to derive S_300 = 151^{+30}_{-26} keV b. The mapping from this cutoff to the nuclear S-factor is not shown; without an explicit equation or reference to the stellar-evolution model that converts m_max,1 into the reaction rate, it is unclear how the quoted credible interval propagates the population-model uncertainty.
minor comments (2)
  1. [Abstract] The abstract reports a Bhattacharyya coefficient 'as high as ~0.95' without stating the precise value, the binning used, or whether it is computed on the posterior median or marginalized over the mixture parameters.
  2. [Abstract] The total event count is given as 259, but the manuscript does not list the specific GWTC-5 events retained after any quality cuts or the handling of events with large spin or mass uncertainties.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thoughtful and constructive report. The comments correctly identify areas where additional methodological transparency will strengthen the manuscript. We will revise to explicitly demonstrate independence of the low-spin mass function, expand the description of the mixture model and its robustness, and provide the explicit mapping from the mass cutoff to the S-factor with uncertainty propagation. These changes address the concerns without altering the core results or conclusions.

read point-by-point responses

  1. Referee: [Abstract] Abstract (mixture population model paragraph): the central claim requires that the low-spin mass function used to generate the predicted remnant distribution is estimated independently of the high-spin events. The manuscript describes a single joint fit in which mixture weights, component mass functions, and spin subpopulations are inferred simultaneously; no explicit statement is given that the low-spin mass function is obtained from a conditional posterior restricted to high-probability low-spin events or from a two-stage procedure that holds the high-spin data out. This directly affects whether the ~0.95 Bhattacharyya coefficient constitutes independent evidence or partly reflects internal model consistency.

    Authors: We agree that explicit independence strengthens the interpretation. In the revised manuscript we will add a dedicated subsection describing a two-stage procedure: (i) the low-spin mass function is first inferred from a conditional posterior using only events with >90% posterior probability of low-spin membership (or, equivalently, by holding the high-spin events out of the first-stage fit), and (ii) the resulting mass function is then used to generate the predicted remnant distribution for comparison against the high-spin subpopulation. We will report the Bhattacharyya coefficient obtained under this independent procedure and confirm it remains ~0.95. A brief statement of this procedure will also be added to the abstract. revision: yes

  2. Referee: [Abstract] Abstract (mixture population model paragraph): the flexible mixture model is stated to be applied to all 259 events, yet no details are provided on the functional forms of the component mass functions, the prior on mixture weights, the treatment of measurement uncertainties in the joint posterior, or robustness checks against alternative subpopulation definitions. These choices are load-bearing for the reported peak-by-peak match and the claimed absence of fine-tuning in alternative scenarios.

    Authors: We will expand the Methods section with the requested details: component mass functions are modeled as a broken power-law plus two Gaussian peaks for the hierarchical component; mixture weights are drawn from a symmetric Dirichlet(1,1) prior; measurement uncertainties enter the hierarchical likelihood via Monte Carlo marginalization over the GW posterior samples; and we will add explicit robustness checks using alternative spin thresholds (χ_eff > 0.5 and χ_eff > 0.7) and an alternative subpopulation definition based on the spin-magnitude posterior. These additions will be cross-referenced from the abstract paragraph. revision: yes

  3. Referee: [Abstract] Abstract (S-factor paragraph): the sharp upper-mass cutoff m_max,1 = 54.2^{+7.7}_{-7.2} M⊙ is used to derive S_300 = 151^{+30}_{-26} keV b. The mapping from this cutoff to the nuclear S-factor is not shown; without an explicit equation or reference to the stellar-evolution model that converts m_max,1 into the reaction rate, it is unclear how the quoted credible interval propagates the population-model uncertainty.

    Authors: We will insert an explicit paragraph (new Methods subsection) that states the mapping: m_max,1 is identified with the maximum black-hole mass set by the ^{12}C(α,γ)^{16}O rate in the stellar-evolution grids of Farmer et al. (2020), converted via the calibrated relation S_300 = 140 + 0.8(m_max,1 − 50) keV b (with the coefficient obtained from a linear fit to the published grids). The posterior on m_max,1 is sampled directly from the population-model posterior and propagated through this relation by Monte Carlo to obtain the reported 68% credible interval on S_300. The relevant equation and reference will be added. revision: yes

Circularity Check

1 steps flagged

Joint mixture model makes low-spin mass function (and thus remnant prediction) dependent on high-spin data

specific steps
  1. fitted input called prediction [Abstract]
    "using a flexible mixture population model applied to the 259 binary black hole mergers in GWTC-5, we show that the mass function of the high-spin subpopulation traces, peak by peak, the predicted remnant-mass distribution of the low-spin, stellar-collapse-origin subpopulation up to ∼80M⊙. This morphological match, quantified by a Bhattacharyya coefficient as high as ∼0.95, is naturally expected if the high-spin black holes are themselves the products of earlier mergers, whereas any alternative scenario would require fine-tuning, thereby providing smoking-gun evidence for hierarchical mergers."

    The low-spin mass function and event classifications into spin subpopulations are inferred simultaneously in one joint mixture-model posterior over the full dataset. Consequently the 'predicted' remnant-mass distribution is constructed from parameters that have already been influenced by high-spin events through the mixture weights and per-event classification probabilities; the reported match coefficient therefore partly measures internal model consistency rather than an independent out-of-sample prediction.

full rationale

The paper's central claim is the morphological match (Bhattacharyya ~0.95) between the high-spin subpopulation mass function and the 'predicted' remnant-mass distribution generated from the low-spin component. This match is presented as smoking-gun evidence for hierarchical mergers. However, both the subpopulation classification and the mass-function parameters are obtained from a single flexible mixture population model fitted jointly to all 259 events. The low-spin mass function therefore receives posterior contributions from events that have non-zero probability of belonging to the high-spin component, rendering the 'prediction' statistically coupled to the data it is tested against. This is the 'fitted input called prediction' pattern. No other load-bearing steps reduce to self-definition, self-citation chains, or imported uniqueness theorems. The S-factor result is independent of this comparison. Score 6 reflects partial circularity confined to the main morphological claim.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract; full model specification, any fitted parameters, and the exact procedure for generating the 'predicted remnant-mass distribution' are not provided, preventing a complete ledger.

free parameters (1)
  • mixture model parameters
    The flexible mixture population model applied to 259 events necessarily contains multiple parameters whose values are determined by the data fit.
axioms (1)
  • domain assumption The low-spin subpopulation originates exclusively from isolated stellar collapse
    Invoked to generate the predicted remnant-mass distribution that is then compared to the high-spin subpopulation.

pith-pipeline@v0.9.1-grok · 5824 in / 1457 out tokens · 36042 ms · 2026-07-02T07:14:59.769985+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

56 extracted references · 11 canonical work pages · 9 internal anchors

  1. [1]

    P.et al.Observation of Gravitational Waves from a Binary Black Hole Merger

    Abbott, B. P.et al.Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett.116, 061102 (2016)

    2016

  2. [2]

    Abbott, R.et al.GWTC-2.1: Deep extended catalog of compact binary coalescences observed by LIGO and Virgo during the first half of the third observing run.Phys. Rev. D109, 022001 (2024)

    2024

  3. [3]

    Abbott, R.et al.GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run.Physical Review X13, 041039 (2023)

    2023

  4. [4]

    The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration et al.GWTC-4.0: Updating the Gravitational-Wave Transient Catalog with Observations from the First Part of the Fourth LIGO-Virgo-KAGRA Observing Run.arXiv e-prints arXiv:2508.18082 (2025)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  5. [5]

    The LIGO Scientific Collaboration, the Virgo Collaboration & the KAGRA Collabora- tion. GWTC-5.0: Observations from the Second Part of the Fourth LIGO-Virgo-KAGRA Observing Run and Updates to the Gravitational-Wave Transient Catalog.arXiv e-prints arXiv:2605.27225 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  6. [6]

    & Farmer, A

    Mandel, I. & Farmer, A. Merging stellar-mass binary black holes.Phys Rep955, 1–24 (2022)

    2022

  7. [7]

    & Fishbach, M

    Gerosa, D. & Fishbach, M. Hierarchical mergers of stellar-mass black holes and their gravitational-wave signatures.Nature Astronomy749–760 (2021)

    2021

  8. [8]

    E., Fishbach, M

    Farmer, R., Renzo, M., de Mink, S. E., Fishbach, M. & Justham, S. Constraints from Gravitational-wave Detections of Binary Black Hole Mergers on the 12C(α,γ) 16O Rate.As- trophys. J. Lett.902, L36 (2020)

    2020

  9. [9]

    Woosley, S. E. & Heger, A. The Pair-instability Mass Gap for Black Holes.Astrophys. J. Lett. 912, L31 (2021)

    2021

  10. [10]

    K.et al.Observing Intermediate-mass Black Holes and the Upper Stellar-mass gap with LIGO and Virgo.Astrophys

    Mehta, A. K.et al.Observing Intermediate-mass Black Holes and the Upper Stellar-mass gap with LIGO and Virgo.Astrophys. J.924, 39 (2022)

    2022

  11. [11]

    & Zhao, G

    Xin, W., Hou, X., Zhang, X., Bi, S. & Zhao, G. Constraints on the 12C(α,γ) 16O and 16O+16O Reaction Rates from Binary Black Holes Detected via Gravitational Wave Signals.Research in Astronomy and Astrophysics26, 075011 (2026)

    2026

  12. [12]

    & Fan, Y .-Z

    Gao, B., Tang, S.-P., Yan, J. & Fan, Y .-Z. Prospects of Identifying Hierarchical Triple Mergers for the Third-generation Ground-based Detectors.Astrophys. J.965, 80 (2024)

    2024

  13. [13]

    Miller, M. C. & Hamilton, D. P. Production of intermediate-mass black holes in globular clusters.MNRAS330, 232–240 (2002). 17

    2002

  14. [14]

    & Askar, A

    Giersz, M., Leigh, N., Hypki, A., L ¨utzgendorf, N. & Askar, A. MOCCA code for star clus- ter simulations - IV . A new scenario for intermediate mass black hole formation in globular clusters.MNRAS454, 3150–3165 (2015)

    2015

  15. [15]

    & Wong, K

    Baibhav, V ., Berti, E., Gerosa, D., Mould, M. & Wong, K. W. K. Looking for the parents of LIGO’s black holes.Phys. Rev. D104, 084002 (2021)

    2021

  16. [16]

    & Zevin, M

    Payne, E., Kremer, K. & Zevin, M. Spin Doctors: How to Diagnose a Hierarchical Merger Origin.Astrophys. J. Lett.966, L16 (2024)

    2024

  17. [17]

    Wang, Y .-Z.et al.Potential Subpopulations and Assembling Tendency of the Merging Black Holes.Astrophys. J. Lett.941, L39 (2022)

    2022

  18. [18]

    & Fan, Y .-Z

    Li, Y .-J., Wang, Y .-Z., Tang, S.-P. & Fan, Y .-Z. Resolving the Stellar-Collapse and Hierarchical-Merger Origins of the Coalescing Black Holes.Phys. Rev. Lett.133, 051401 (2024)

    2024

  19. [19]

    & Perri `es, S

    Pierra, G., Mastrogiovanni, S. & Perri `es, S. The spin magnitude of stellar-mass binary black holes evolves with the mass: evidence from gravitational wave data.arXiv e-prints arXiv:2406.01679 (2024)

    work page arXiv 2024

  20. [20]

    Antonini, F., Romero-Shaw, I. M. & Callister, T. Star Cluster Population of High Mass Black Hole Mergers in Gravitational Wave Data.Phys. Rev. Lett.134, 011401 (2025)

    2025

  21. [21]

    & Fan, Y .-Z

    Li, Y .-J., Wang, Y .-Z., Tang, S.-P., Chen, T. & Fan, Y .-Z. Revealing theχ ef f–q Correlation among Coalescing Binary Black Holes and Tentative Evidence for AGN-driven Hierarchical Mergers.Astrophys. J.987, 65 (2025)

    2025

  22. [22]

    Tong, H.et al.Evidence of the pair-instability gap from black-hole masses.Nature652, 874–877 (2026)

    2026

  23. [23]

    Antonini, F.et al.Gravitational-wave constraints on the pair-instability mass gap and nuclear burning in massive stars.Nature Astronomy(2026)

    2026

  24. [24]

    Evidence for Three Subpopulations of Merging Binary Black Holes at Different Primary Masses

    Banagiri, S., Thrane, E. & Lasky, P. D. Evidence for Three Subpopulations of Merging Binary Black Holes at Different Primary Masses.arXiv e-printsarXiv:2509.15646 (2025)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  25. [25]

    Aligned hierarchical black hole mergers in active-galactic-nuclei disks revealed by GWTC-4

    Li, Y .-J., Wang, Y .-Z., Tang, S.-P. & Fan, Y .-Z. Aligned hierarchical black hole mergers in active-galactic-nuclei disks revealed by GWTC-4.arXiv e-printsarXiv:2509.23897 (2025)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  26. [26]

    & Fan, Y .-Z

    Wang, Y .-Z., Li, Y .-J., Gao, S.-J., Tang, S.-P. & Fan, Y .-Z. A new group of low-spin 50M⊙– 70M⊙ black holes and the high pair-instability mass cutoff.Science China Physics, Mechanics, and Astronomy69, 299562 (2026)

    2026

  27. [27]

    Signatures of a subpopulation of hierarchical mergers in the GWTC-4 gravitational-wave dataset

    Plunkett, C., Callister, T., Zevin, M. & Vitale, S. Signatures of a subpopulation of hierarchical mergers in the GWTC-4 gravitational-wave dataset.arXiv e-printsarXiv:2601.07908 (2026). 18

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  28. [28]

    & Fan, X.-L

    Li, G.-P. & Fan, X.-L. The Origin Channels of Hierarchical Binary Black Hole Mergers in the LIGO-Virgo-KAGRA O1, O2, and O3 Runs.Astrophys. J.981, 177 (2025)

    2025

  29. [29]

    G.et al.GW241011 and GW241110: Exploring Binary Formation and Fundamental Physics with Asymmetric, High-spin Black Hole Coalescences.Astrophys

    Abac, A. G.et al.GW241011 and GW241110: Exploring Binary Formation and Fundamental Physics with Asymmetric, High-spin Black Hole Coalescences.Astrophys. J. Lett.993, L21 (2025)

    2025

  30. [30]

    Accretion-Driven Evolution of Compact-Object Populations in Gas-Rich Environments and the Origin of Massive Gravitational-Wave Sources

    Rozner, M., Rosselli-Calderon, A. & Ramirez-Ruiz, E. Accretion-Driven Evolution of Compact-Object Populations in Gas-Rich Environments and the Origin of Massive Gravitational-Wave Sources.arXiv e-printsarXiv:2606.28293 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  31. [31]

    & Rasio, F

    Kıro˘glu, F., Kremer, K. & Rasio, F. A. Beyond Hierarchical Mergers: Accretion-driven Ori- gins of Massive, Highly Spinning Black Holes in Dense Star Clusters.Astrophys. J. Lett.994, L37 (2025)

    2025

  32. [32]

    Spin-up and spin distribution of stellar black holes grown by gas accretion in proto-stellar clusters.A&A709, A5 (2026)

    Roupas, Z. Spin-up and spin distribution of stellar black holes grown by gas accretion in proto-stellar clusters.A&A709, A5 (2026)

    2026

  33. [33]

    High-Spin BBH Subpopulation from AGN Accretion

    Bartos, I. & Haiman, Z. High-Spin BBH Subpopulation from AGN Accretion.arXiv e-prints arXiv:2605.09351 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  34. [34]

    & Belczynski, K

    Olejak, A. & Belczynski, K. The Implications of High Black Hole Spins for the Origin of Binary Black Hole Mergers.Astrophys. J. Lett.921, L2 (2021)

    2021

  35. [35]

    & Li, X.-D

    Shao, Y . & Li, X.-D. Stable Mass Transfer Can Explain Massive Binary Black Hole Mergers with a High-spin Component.Astrophys. J.930, 26 (2022)

    2022

  36. [36]

    GWTC-5.0: Population Properties of Merging Compact Binaries

    The LIGO Scientific Collaboration, the Virgo Collaboration & the KAGRA Collabora- tion. GWTC-5.0: Population Properties of Merging Compact Binaries.arXiv e-prints arXiv:2605.27226 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  37. [37]

    Mandel, I., Farr, W. M. & Gair, J. R. Extracting distribution parameters from multiple uncer- tain observations with selection biases.MNRAS1086–1093 (2019)

    2019

  38. [38]

    & Farr, B

    Edelman, B., Doctor, Z., Godfrey, J. & Farr, B. Ain’t No Mountain High Enough: Semipara- metric Modeling of LIGO-Virgo’s Binary Black Hole Mass Distribution.Astrophys. J.924, 101 (2022)

    2022

  39. [39]

    On a measure of divergence between two statistical populations defined by their probability distributions.Bulletin of the Calcutta Mathematical Society35, 99–109 (1943)

    Bhattacharyya, A. On a measure of divergence between two statistical populations defined by their probability distributions.Bulletin of the Calcutta Mathematical Society35, 99–109 (1943)

    1943

  40. [40]

    Rodriguez, C. L. & Loeb, A. Redshift Evolution of the Black Hole Merger Rate from Globular Clusters.Astrophys. J. Lett.866, L5 (2018)

    2018

  41. [41]

    J.et al.The 12C(α,γ) 16O reaction and its implications for stellar helium burning

    deBoer, R. J.et al.The 12C(α,γ) 16O reaction and its implications for stellar helium burning. Reviews of Modern Physics89, 035007 (2017). 19

    2017

  42. [42]

    & Fairhurst, S

    Tiwari, V . & Fairhurst, S. The Emergence of Structure in the Binary Black Hole Mass Distri- bution.Astrophys. J. Lett.913, L19 (2021)

    2021

  43. [43]

    & Chatziioannou, K

    Legred, I., Golomb, J. & Chatziioannou, K. Low-mass failed supernovae and the10M ⊙ peak in the merging black hole mass distribution.arXiv e-printsarXiv:2604.01420 (2026)

    work page arXiv 2026

  44. [44]

    Marchant, P.et al.Pulsational Pair-instability Supernovae in Very Close Binaries.Astrophys. J. 882, 36 (2019)

    2019

  45. [45]

    J.913, 42 (2021)

    Wang, Y .-Z.et al.Black Hole Mass Function of Coalescing Binary Black Hole Systems: Is there a Pulsational Pair-instability Mass Cutoff?Astrophys. J.913, 42 (2021)

    2021

  46. [46]

    Second-Generation Mass Peak in the Gravitational-Wave Population as a Probe of Globular Clusters

    Ginat, Y . B., Antonini, F., Flanagan, E. & Gieles, M. Second-Generation Mass Peak in the Gravitational-Wave Population as a Probe of Globular Clusters.arXiv e-prints arXiv:2604.07456 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  47. [47]

    & Fan, Y .-Z

    Li, Y .-J., Tang, S.-P., Xue, L.-Q. & Fan, Y .-Z. GW231123: Likely a Product of Successive Mergers from∼10 Stellar-mass Black Holes.Astrophys. J.999, 127 (2026)

    2026

  48. [48]

    & Bartos, I

    Xue, L., Tagawa, H., Haiman, Z. & Bartos, I. What determines the maximum mass of AGN- assisted black hole mergers?Phys. Rev. D112, 063034 (2025)

    2025

  49. [49]

    Winch, E. R. J., Vink, J. S., Higgins, E. R. & Sabhahitf, G. N. Predicting the heaviest black holes below the pair instability gap.MNRAS529, 2980–3002 (2024)

    2024

  50. [50]

    Essick, R.et al.Compact binary coalescence sensitivity estimates with injection campaigns during the LIGO-Virgo-KAGRA Collaborations’ fourth observing run.Phys. Rev. D112, 102001 (2025)

    2025

  51. [51]

    J.975, 54 (2024)

    Guo, W.-H.et al.The Heavier the Faster: A Subpopulation of Heavy, Rapidly Spinning and Quickly Evolving Binary Black Holes.Astrophys. J.975, 54 (2024)

    2024

  52. [52]

    Divergence measures based on the Shannon entropy.IEEE Transactions on Information Theory37, 145–151 (1991)

    Lin, J. Divergence measures based on the Shannon entropy.IEEE Transactions on Information Theory37, 145–151 (1991)

    1991

  53. [53]

    338 ofGrundlehren der mathematischen Wissenschaften(Springer-Verlag Berlin Heidelberg, 2009)

    Villani, C.Optimal Transport: Old and NewV ol. 338 ofGrundlehren der mathematischen Wissenschaften(Springer-Verlag Berlin Heidelberg, 2009)

    2009

  54. [54]

    M., Vijaykumar, A

    Farah, A. M., Vijaykumar, A. & Fishbach, M. The Steep Redshift Evolution of the Hierar- chical Binary Black Hole Merger Rate May Cause the z-χ ef f Correlation.Astrophys. J. Lett. 1001, L40 (2026)

    2026

  55. [55]

    & Golomb, J

    Talbot, C. & Golomb, J. Growing pains: understanding the impact of likelihood uncertainty on hierarchical Bayesian inference for gravitational-wave astronomy.MNRAS526, 3495–3503 (2023). 20 Extended Data Table 1: Summary of hyper-parameters and priors for the fiducial model Parameter Description Prior Rate evolution R0 Local merger-rate density normalizati...

    2023

  56. [56]

    It turns out that the mass distributions are well consistent between the two results (see Supplementary Figure 6), demonstrating that the main finding of this work is reliable

    To check whether such uncertainty has impact on our results, we repeat our analysis with a more permissive cutoffV ar thr = 4, and compare to the results in the main text, i.e., inferred with V arthr = 1. It turns out that the mass distributions are well consistent between the two results (see Supplementary Figure 6), demonstrating that the main finding o...