Recognition: 4 theorem links
· Lean TheoremGWTC-4.0: Population Properties of Merging Compact Binaries
Pith reviewed 2026-05-09 00:36 UTC · model claude-opus-4-7
The pith
From 158 gravitational-wave mergers, the black-hole mass spectrum shows persistent bumps at 10 and 35 solar masses, modest spins, and a quarter of binaries with spin pointing against the orbit.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using the cumulative catalog of 158 confidently detected compact-binary mergers through the first half of the fourth observing run, the authors map the joint distribution of masses, spins, and redshift for binary neutron stars, neutron-star–black-hole binaries, and binary black holes. The black-hole mass distribution is not a featureless power law: there are clear excesses near 10 M_sun and 35 M_sun on top of a continuum that steepens above 35 M_sun, with hints of further structure near 20 and 60 M_sun. Binaries with a ~10 M_sun primary prefer unequal pairings (q≈0.74), consistent with a stable-mass-transfer origin. Black-hole spins are predominantly small (90% have χ<0.57) and tilted toward
What carries the argument
Hierarchical Bayesian inference on the catalog-level likelihood, run with two complementary families of population models: strongly modeled (parametric power-laws, broken power-laws, and Gaussian peaks for mass; truncated Gaussians and skew-normals for spin; (1+z)^κ for redshift) and weakly modeled (B-splines and binned Gaussian processes that let the data sculpt the shape). Selection effects are folded in through Monte-Carlo integrals over injection campaigns from the search pipelines, with explicit checks that the likelihood estimator has converged. Cross-checking strongly and weakly modeled fits is what lets the authors claim that a feature is in the data rather than in the model.
If this is right
- Standard core-collapse supernova prescriptions with a sharp NS/BH boundary need revision: the lower mass gap between roughly 2.5 and 5 M_sun is not empty, and the BH distribution rises into a peak near 9–10 M_sun rather than starting from a clean threshold.
- Stable mass transfer (rather than common-envelope-only evolution or dynamical assembly) is favored as a major channel for the ~10 M_sun population, because that channel naturally predicts both a peak near this mass and a mass-ratio mode below unity.
- A nonzero fraction (0.24–0.42) of mergers with negative effective inspiral spin requires that a substantial subpopulation forms through a process that randomizes spin orientations — i.e., dynamical assembly in gas-poor environments — coexisting with a dominantly aligned isolated channel.
- The merger rate growing as (1+z)^3.2 — steeper than the cosmic star formation rate — points to either a metallicity bias toward low-Z progenitors or short delay times with a long tail, constraining population-synthesis delay-time distributions.
- The χ_eff distribution broadens with redshift out to z≈1 while the mass spectrum shows no detectable redshift evolution, suggesting that channel mixture (not stellar physics shifting with cosmic time) drives the spin evolution.
Where Pith is reading between the lines
- The persistence of the ~10 M_sun peak with a mass-ratio distribution peaking around q≈0.74 is hard to reproduce with pure dynamical formation, which tends to favor equal-mass pairings, suggesting at least this part of the population is shaped by a binary-evolution channel with mass-ratio reversal.
- The combination of asymmetric, positively skewed effective spin and a sizable negative-spin tail is naturally read as two channels coexisting: a small-spin, partially aligned isolated channel plus an isotropic dynamical channel — a mixture that is now starting to be resolvable rather than degenerate.
- The broadening (rather than shifting) of the χ_eff distribution with redshift is a distinctive signature: it disfavors a single channel with redshift-dependent spin-up and instead favors a changing channel mix, which future catalogs should be able to break by binning in mass.
- If the steepening above 35 M_sun and tentative rise near 60 M_sun both hold up, the upper mass gap is not empty but partially filled, which would constrain pair-instability physics and/or argue for a non-trivial hierarchical-merger contribution at the high-mass end.
Load-bearing premise
The whole picture rests on the search pipelines' simulated-injection campaigns correctly capturing how often each kind of binary would be detected; if that sensitivity model is biased — especially at high mass, low mass ratio, or high redshift — the inferred peaks, gaps, and rate evolution shift accordingly.
What would settle it
If a re-analysis with corrected calibration priors and independent sensitivity injections fails to reproduce a local maximum in the primary-mass distribution near 10 M_sun and a steepening above ~35 M_sun, or if the negative-χ_eff fraction collapses below ~10% once individual-event posteriors are reweighted, the two-channel reading and the stable-mass-transfer interpretation lose their support.
read the original abstract
We detail the population properties of merging compact objects using 158 mergers from the cumulative Gravitational-Wave Transient Catalog 4.0, which includes three types of binary mergers: binary neutron star, neutron star--black hole binary, and binary black hole mergers. We resolve multiple over- and under-densities in the black hole mass distribution: features persist at primary masses of $10\,M_\odot$ and $35\,M_\odot$ with a possible third feature at $\sim 20\,M_\odot$. These are departures from an otherwise power-law-like continuum that steepens above $35\,M_\odot$. Binary black holes with primary masses near $10\,M_\odot$ are more likely to have less massive secondaries, with a mass ratio distribution peaking at $q = 0.74^{+0.13}_{-0.13}$, potentially a signature of stable mass transfer during binary evolution. Black hole spins are inferred to be non-extremal, with 90\% of black holes having $\chi < 0.57$, and preferentially aligned with binary orbits, implying many merging binaries form in isolation. However, we find a significant fraction, 0.24-0.42, of binaries have negative effective inspiral spins, suggesting many could be formed dynamically in gas-free environments. We find evidence for correlation between effective inspiral spin and mass ratio, though it is unclear if this is driven by variation in the mode of the distribution or the width. (Abridged)
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents the population analysis of compact binary mergers from the cumulative GWTC-4.0 catalog (158 mergers; 153 BBHs at FAR<1/yr and an additional BNS/NSBH set at FAR<0.25/yr including GW230529). Using a battery of strongly-modeled (broken power-law + 2 peaks, FullPop-4.0, Gaussian/Skew-normal effective spin, copula and linear/spline correlation models) and weakly-modeled (B-Spline, BGP, AR, FM) approaches, the authors infer (i) the joint mass spectrum from BNS to BBH including a non-empty NS–BH transition with peaks near 1.5 and 9 M☉; (ii) BBH primary-mass features at ~10 and ~35 M☉ on a broken power-law continuum (α1≈1.7, α2≈4.5); (iii) a mass-ratio peak near q=0.74 associated with the 10 M☉ subpopulation under the Isolated Peak model; (iv) a non-extremal, mildly aligned spin distribution with 24–42% of systems having χ_eff<0 and a positively-skewed χ_eff distribution; (v) merger-rate evolution (1+z)^κ with κ=3.2^{+0.94}_{-1.00}; and (vi) evidence for χ_eff distribution broadening with redshift and softened evidence (relative to GWTC-3) for q–χ_eff anti-correlation. Local merger rates are quoted as 7.6–250, 9.1–84, and 14–26 Gpc^-3 yr^-1 for BNS, NSBH, and BBH respectively.
Significance. If the central inferences hold, this is the definitive population characterization of compact binary mergers through O4a and the standard reference for the field. The paper substantially tightens constraints on the BBH mass continuum (in particular, the steepening above ~35 M☉ at α2>α1 with 97.7% credibility, Fig. 4), demonstrates that an "empty" lower mass gap between NSs and BHs is disfavored by both strongly- and weakly-modeled approaches (Fig. 2 inset on the gap-depth A), and provides an updated, more robust measurement of the merger rate evolution κ. The methodological strengths are explicit: transparent likelihood-variance thresholding (App. A.1, Eq. A3, σ²_lnL<1), parallel strongly- and weakly-modeled inferences for every headline result, model-comparison tables, and a public data release with figure-generation scripts. The collaboration's continued use of multiple independent parameterizations (e.g., Linear, Spline, Copula for q–χ_eff) makes the systematic-vs-statistical decomposition unusually traceable for a population paper.
major comments (6)
- [Abstract; §6.1; App. D.1] The abstract presents three apparently distinct mass-spectrum findings — a feature at 35 M☉, a continuum that 'steepens above 35 M☉,' and a possible feature at 20 M☉ — but §6.1 explicitly states that the 35 M☉ feature 'is consistent with either: (i) an over-density … relative to an underlying broken power law, or (ii) a broken power law with break mass m_break=34.1^{+3.8}_{-3.3} M☉ (i.e., a broken power law without a second peak near the break).' Under interpretation (ii), the '35 M☉ feature' and the 'steepening above 35 M☉' are the same phenomenon (the α1→α2 break), not two independent findings. Please reword the abstract to make this degeneracy explicit, and quantify in §6.1 the Bayes factor / posterior weight between the two interpretations rather than only stating consistency.
- [Abstract; §6.1] The 20 M☉ feature is reported in the abstract as a 'possible third feature,' but §6.1 reports log10 B = −0.34 against a model adding a third Gaussian to the strongly-modeled fit, i.e., the data marginally disfavor it under that parameterization, with positive support coming only from some weakly-modeled reconstructions (Fig. 20). The current abstract phrasing risks overstating the evidence. Please either (a) quote the Bayes factor in the abstract, or (b) reclassify the 20 M☉ structure as a feature seen only in weakly-modeled reconstructions and not required by strongly-modeled fits.
- [§6.2 (Eq. B24, Fig. 6); Abstract] The headline q = 0.74^{+0.13}_{-0.13} pairing claim for the ~10 M☉ subpopulation is derived from the Isolated Peak model, while the Extended Broken Power Law + 2 Peaks model — which is structurally closer to the fiducial mass model — yields β^{peak1}_q = 1.6^{+9.6}_{-8.1}, essentially uninformative. Because the Isolated Peak model fits the q distribution of the peak component with a B-spline, it has substantially more flexibility than a single power-law β_q and can absorb q-structure that the parametric pairing function cannot. Please (i) report the Bayes factor between Isolated Peak and Extended BP+2P, (ii) test whether a mid-flexibility parametric pairing for the peak (e.g., a Beta or truncated Gaussian in q) reproduces q=0.74, and (iii) soften the abstract's 'mass ratio distribution peaking at q=0.74' to attribute the result to the specific model and to clarify that the alternative pa
- [§6.3.1, Fig. 7; App. D.2] The shift from a non-singular Beta distribution (GWTC-3 default) to a truncated Gaussian (Gaussian Component Spins) for χ is well-motivated, but it materially changes what 'p(χ≈0)>0' means and is responsible for some of the apparent qualitative change relative to GWTC-3. The disagreement between the strongly- and weakly-modeled models near χ≳0.8 (statement that 'p(χ) distributions differ at the 85%, 96%, and 98% levels' at χ=0.8, 0.9, 1.0) suggests the truncated Gaussian does not have enough flexibility at the boundaries to track the data. Given that GW231123 is identified as a high-spin event, please show explicitly how the inferred 90th-percentile bound χ<0.57 depends on the parametric choice — ideally a comparison with the Unconstrained Beta and with a mixture model — so the headline number is not driven by parametric stiffness.
- [§6.3.2, Fig. 9; App. D.3] The reported skewness ϵ<0 at 99.3% credibility (positive skew of χ_eff) is a strong claim, but the χ_p inference is acknowledged in App. D.3 to depend on analysis settings (likelihood-convergence threshold, inference choices). Because ϵ is fit jointly with location and scale of χ_eff and with the (uncorrelated) χ_p marginal, please demonstrate that the skewness conclusion is robust to (a) tightening/loosening the σ²_lnL<1 threshold, and (b) replacing the Skew-normal by a two-component Gaussian mixture in χ_eff. The current result risks being conflated with a sub-population structure that the skew-normal cannot resolve.
- [§6.5.1 (Linear, Spline, Copula); Fig. 11–12] The q–χ_eff correlation results are model-dependent in a way that should be foregrounded: the Linear model finds δμ_eff|q<0 at only 82% credibility (down from 98% in GWTC-3) but δlnσ_eff|q<0 at 95%, while the Copula gives κ_{q,eff}<0 at 92%. The text in §6.5.1 reads as if the finding is robust because something is non-zero in each model, but the three models are testing different functional forms. Please add a unified statement of what is shared across models (e.g., 'the χ_eff distribution at q≪1 differs from that at q→1 in mean and/or width') with a quantitative cross-model summary, and update the abstract sentence about the q–χ_eff correlation accordingly — currently the abstract's framing is more confident than §6.5.1 supports.
minor comments (10)
- [§3.2.1] The expected contamination 'P_k FAR×T_k ≃ 6.7' assumes independent noise events. Please clarify what fraction of the 161 candidates that contamination represents and whether any of the population conclusions are sensitive to a 1–2 event misclassification at the BBH/NSBH boundary (in particular the BBH mass-spectrum tail below ~8 M☉).
- [§3.4] The disclosure that a normalization error in the PE inner product (Talbot et al. 2025b) and an incorrect calibration prior were discovered, but uncorrected samples are used in this work, is appropriate but should be more prominent. A short paragraph quantifying the magnitude of the resulting per-event posterior shifts (e.g., in χ_eff, q, masses) for representative events would help readers calibrate expectations for population-level systematic error.
- [§4, Fig. 1–2; Table 2] The discontinuity in dR/dm2 at m2=5 M☉ in FullPop-4.0 (because of the switched pairing function in Eq. B14) is visible in Fig. 2 and is a model artefact rather than a physical feature. Please add a sentence to the figure caption stating this explicitly.
- [§4.2] Quoting the BNS rate as the union of FullPop-4.0 and BGP 90% CIs (7.6–250 Gpc^-3 yr^-1) yields a very wide range that effectively spans an order of magnitude. The 'Simple Uniform BNS' rate (13–170) is more comparable to the literature. Consider quoting the Simple Uniform BNS rate as the headline BNS number in the abstract since it is less sensitive to pairing-function choices, and relegating the union range to Table 2.
- [§5] The NS Peak model location μ=1.4^{+0.48}_{-0.25} M☉ and width σ=0.68^{+1.2}_{-0.45} M☉ are largely prior-dominated. It would be helpful to plot the prior-vs-posterior comparison for these hyperparameters to make the prior dominance explicit.
- [§6.4, Fig. 10] The κ=3.2 result is consistent with both κ_SFR=2.7 and a steeper evolution. Please state explicitly the credibility at which κ=κ_SFR is included in the posterior; the current 'consistent with' phrasing is unquantified.
- [§6.5.4] The non-detection of mass–redshift correlation should be qualified by noting that only 19/153 events have z>1 at 90% credibility — sensitivity to mass-spectrum evolution beyond z~1 is intrinsically limited by the dataset, as the authors note for χ_eff in Fig. 14 with the 90%-cumulative line. A similar visual indicator would be useful in §6.5.4.
- [App. A.1, Eq. A3–A6] The σ²_lnL<1 threshold is conservative but the impact of the cut on hyperparameter posteriors should be visualized for at least one model where it removes a non-trivial region (App. D.3 hints at this for χ_p). A before/after corner plot in the supplement would make the systematic transparent.
- [Throughout] Several statements like 'the merger rate is best constrained at ~9 M☉' or 'we observe a peak at 9.02^{+0.41}_{-1.21} M☉' would benefit from explicit comparison to the same quantities measured on GWTC-3 alone, to separate genuine new constraints from re-fits with a different model. Fig. 4 does this for the broken-power-law indices and is excellent — please extend that style of comparison to the peak parameters.
- [Author list] Standard LVK practice — no action required, but flag for the editor that the author list and acknowledgements are extensive.
Simulated Author's Rebuttal
We thank the referee for a careful and constructive report. The recommendation of minor revision and the recognition of the parallel strongly- and weakly-modeled methodology, the likelihood-variance thresholding, and the public data release are appreciated. The six major comments fall into two broad categories: (i) places where the abstract or §-level summaries are more confident than the body of the paper actually supports (35 M☉ feature vs. continuum break; 20 M☉ feature; q=0.74 pairing; q–χ_eff correlation), and (ii) requests for explicit robustness checks against parametric choices (truncated Gaussian vs. Beta vs. mixture for χ; skew-normal vs. two-Gaussian mixture for χ_eff; σ²_lnL threshold sensitivity). We agree with the spirit of all six points. We will reword the abstract on three of the headline claims, add Bayes factors and cross-model summaries to §6.1, §6.2, §6.3.1, §6.3.2 and §6.5.1, and add the requested parametric robustness comparisons in Appendices D.2 and D.3. None of the changes alter the central conclusions of the paper, but they materially improve the traceability between abstract-level statements and the underlying model-by-model evidence, which is the referee's central concern.
read point-by-point responses
-
Referee: Abstract conflates the 35 M☉ feature and the steepening above 35 M☉ as two findings, when §6.1 admits they may be a single phenomenon (an α1→α2 break). Quantify the Bayes factor between the over-density and pure-break interpretations and reword the abstract.
Authors: We agree this is a real ambiguity in the current wording. The Broken Power Law + 2 Peaks fit and a model in which the second Gaussian is removed and only the broken power law remains both describe the data well; as the referee notes, in the latter the break mass is m_break=34.1^{+3.8}_{-3.3} M☉, coincident with the location of the second peak. We will: (i) reword the abstract to say 'a power-law-like continuum that steepens at ~35 M☉, with an additional over-density at this mass not required by the data,' making clear that the 35 M☉ feature and the steepening may be the same phenomenon; and (ii) quote the Bayes factor between the two interpretations explicitly in §6.1 (currently only described qualitatively in App. D.1) so readers can judge the degeneracy quantitatively. The headline conclusion — that the continuum is not a single power law and that the rate steepens above ~35 M☉ — is unchanged. revision: yes
-
Referee: The 20 M☉ feature is reported in the abstract as a 'possible third feature' but the strongly-modeled three-Gaussian extension is marginally disfavored (log10 B = -0.34). Quote the Bayes factor in the abstract or restrict the claim to weakly-modeled reconstructions.
Authors: We accept this point. The data do not require a third Gaussian under the strongly-modeled fit; the feature is visible in some weakly-modeled reconstructions (Fig. 20) but is not a clean detection. We will revise the abstract to read along the lines of 'a tentative additional feature near ~20 M☉ seen in some weakly-modeled reconstructions but not required by strongly-modeled fits (log10 B = -0.34 against a third Gaussian),' so the level of evidence is transparent at first reading. The body of §6.1 already states the Bayes factor; we will lift that number into the abstract and into the §6.1 summary sentence. revision: yes
-
Referee: The headline q=0.74 pairing for the ~10 M☉ subpopulation is from the Isolated Peak model (B-spline in q for the peak component), while the structurally closer Extended BP+2P gives β^{peak1}_q=1.6^{+9.6}_{-8.1}, uninformative. Report the Bayes factor between the two models, test a mid-flexibility parametric pairing, and soften the abstract.
Authors: We agree the abstract should attribute the q=0.74 result to the model that produced it. The Extended BP+2P model is a power-law-only pairing function for each subpopulation, which by construction cannot represent a peak in q at intermediate values; the large posterior on β^{peak1}_q reflects this rather than disagreement with the Isolated Peak result. We will (i) revise the abstract to read 'under the Isolated Peak model, the mass-ratio distribution of the ~10 M☉ subpopulation peaks at q=0.74^{+0.13}_{-0.13}, a feature a single power-law pairing cannot reproduce,' and (ii) add to §6.2 the Bayes factor between Isolated Peak and Extended BP+2P. We will also note that a mid-flexibility test (Beta or truncated Gaussian pairing for the peak component) is a natural follow-up; we expect it to recover the same qualitative feature, but a full analysis is deferred as it requires re-running the model-comparison pipeline and the conclusion of the paper does not depend on it. revision: partial
-
Referee: The change from non-singular Beta to truncated Gaussian in p(χ) materially alters interpretation of 'p(χ≈0)>0' and may not be flexible enough at χ≳0.8 (where the strongly- and weakly-modeled distributions disagree at the 85–98% level). Show how χ<0.57 depends on the parametric form, including comparison to an unconstrained Beta and a mixture.
Authors: We will add explicit comparison plots and percentile values for the Unconstrained Beta and a two-component truncated Gaussian mixture in App. D.2, alongside the current Gaussian Component Spins and B-Spline results, so the reader can see the parametric dependence of the χ<0.57 statement. The B-Spline result, which has substantially more flexibility at the boundaries, gives a comparable but not identical bound; we will report both numbers in §6.3.1 and revise the headline phrasing to 'χ<0.57 under the Gaussian Component Spins model, with consistent values from the weakly-modeled fit.' We agree the discrepancy at χ≳0.8 is a genuine model-flexibility issue, partly driven by GW231123, and we will state this in the body rather than only in passing. revision: yes
-
Referee: The skewness ϵ<0 at 99.3% credibility is a strong claim but is fit jointly with χ_p, which §App. D.3 acknowledges depends on analysis settings. Demonstrate robustness to the σ²_lnL<1 threshold and to replacing the skew-normal with a two-component Gaussian mixture in χ_eff.
Authors: We will add a robustness study in App. D.3 showing the posterior on ϵ at σ²_lnL thresholds of 0.5, 1, and 2 — the χ_eff inference is much less sensitive to this threshold than χ_p (the χ_eff prior has support at zero and the per-event posteriors are far better constrained), and we expect the 99.3% credibility statement to be stable, but we agree it should be demonstrated rather than asserted. We will also fit a two-component Gaussian mixture in χ_eff and report whether the data prefer a single skewed component, two Gaussians, or are agnostic. If the mixture is preferred, we will reframe the result in §6.3.2 as 'asymmetry consistent with either a skewed single population or a sub-population structure,' which is a weaker but more honest statement. revision: yes
-
Referee: The q–χ_eff correlation results are model-dependent: Linear gives δμ_eff|q<0 at 82%, δlnσ_eff|q<0 at 95%; Copula gives κ_{q,eff}<0 at 92%. The text reads as if the finding is robust, but the three models test different functional forms. Add a unified cross-model statement and soften the abstract.
Authors: We agree the §6.5.1 narrative could be read as more confident than the individual numbers warrant, and that the abstract sentence ('We find evidence for correlation between effective inspiral spin and mass ratio, though it is unclear if this is driven by variation in the mode of the distribution or the width') already gestures at this but does not commit to what is shared across models. We will add a short cross-model summary table to §6.5.1 listing, for each of Linear, Spline, and Copula: the parameter being constrained, the credibility level, and a one-line interpretation. The shared statement, which is supported across all three, is that the χ_eff distribution at low q differs from that at q→1 (in mean, width, or both), with the strongest single statement being δlnσ_eff|q<0 at 95% in Linear. We will revise the abstract sentence to attribute the evidence to this shared finding rather than implying a clean anti-correlation in the mode. revision: yes
- The referee's request to test a mid-flexibility parametric pairing function (Beta or truncated Gaussian in q) for the ~10 M☉ peak component is a sensible follow-up, but a full re-run with the production model-comparison pipeline is beyond what we can deliver in this revision cycle. We will note this as a planned follow-up rather than include the result in the revised manuscript.
Circularity Check
No significant circularity: GWTC-4.0 population inference is a standard hierarchical Bayesian analysis on independently calibrated detector data with externally-validated selection functions.
full rationale
This is an observational catalog paper performing hierarchical Bayesian inference on gravitational-wave events. The derivation chain is: (i) detector strain data → (ii) search pipelines with FAR thresholds → (iii) per-event parameter estimation posteriors using waveform models → (iv) population-level likelihood (Eq. 2) marginalized over individual-event parameters with a selection function ξ(Λ) estimated by Monte Carlo injection campaigns (Eq. 3, App. A). None of these steps reduces a 'prediction' to its own input by construction. The selection function is computed from injected synthetic signals analyzed by the same pipelines that found the real events — this is the correct procedure, not circularity. Self-citations to prior LVK papers (GWTC-3.0, GWTC-2.1, etc.) provide datasets and methods, but the central claims (mass spectrum features, spin distributions, merger rates) are inferred from the data via models whose hyperparameters are sampled, not assumed. The skeptic's concern that the '35 M☉ feature' is degenerate with a power-law break is acknowledged transparently in §6.1 and App. D.1 (log₁₀ B = −0.34 for adding a third peak) — this is honest model comparison, not a circular reduction. The q = 0.74 mass-ratio peak is model-dependent (Isolated Peak vs. Extended Broken Power Law + 2 Peaks giving β=1.6^{+9.6}_{-8.1}), again disclosed openly. Model dependence and degeneracy ≠ circularity. The waveform models (NRSur7dq4, IMRPhenomXPHM, SEOBNRv5PHM) are calibrated against numerical relativity simulations external to GW detection. Priors on hyperparameters are stated and uniform/weakly-informative. I find no step where a fitted parameter is renamed as a prediction, no step where a quantity is defined in terms of itself, and no load-bearing self-citation that imports a 'uniqueness theorem' to forbid alternatives. The minor caveat is normal heavy reliance on prior LVK methodology papers — appropriate for a collaboration catalog paper and not circular in the technical sense.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
Foundation.LogicAsFunctionalEquation / Cost.FunctionalEquationwashburn_uniqueness_aczel (J-cost uniqueness) unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We determine the population properties of merging compact binaries using hierarchical Bayesian inference... p(Λ|d) = L(d|Λ)π(Λ)/p(d)
-
Masses.SectorParams / phi-ladder rung modulesphi_fixed_point (golden-ratio rung spacing) unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
features persist at primary masses of 10 M☉ and 35 M☉ with a possible third feature at ∼20 M☉ ... departures from an otherwise power-law-like continuum that steepens above 35 M☉
-
Foundation.BranchSelection / AlphaCoordinateFixationbranch_selection (coupling combiner forces bilinear branch) unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
B-Spline... Binned Gaussian Process... Autoregressive Process... Flexible Mixtures — weakly modeled approaches
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 50 Pith papers
-
GW240925 and GW250207: Astrophysical Calibration of Gravitational-wave Detectors
The first informative astrophysical calibration of gravitational-wave detectors is reported using GW240925 and GW250207.
-
Secondary-Mass Features improve Spectral-Siren $H_0$ Constraints
A new model emphasizing secondary mass features and pairing transitions improves spectral siren H0 constraints by ~30% using 142 GW events from GWTC-4.0.
-
End-to-End Population Inference from Gravitational-Wave Strain using Transformers
Dingo-Pop uses a transformer to perform amortized, end-to-end population inference from GW strain data in seconds, bypassing per-event Monte Carlo sampling.
-
On the Origin of Mass Ejection in Failed Supernovae
Higher-Mach-number self-similar shock solutions in failed supernovae are unstable and strengthen asymptotically above a critical neutrino mass-loss threshold, explaining greater ejection in red supergiants versus comp...
-
Black Hole Supernovae Outcomes Across a Wide Progenitor Range
Black hole supernovae occur across a wide progenitor mass range from 19.5 to 60 solar masses, yielding final black hole masses of 3 to 26 solar masses that trend with but are not fully set by CO core mass.
-
Constraining Dipole Radiation with Multiband Gravitational Waves from Eccentric Binary Black Holes
Multiband observations of eccentric binary black holes can constrain dipole-radiation deviations from general relativity to |b| ≲ 10^{-7} for a GW231123-like event when combining one year of space-based data with grou...
-
How lonely are the Binary Compact Objects Detected by the LIGO-Virgo-KAGRA Collaboration?
No three-body encounter signatures detected in GW170817, GW190814, and GW230627_015337, constraining intermediate-mass black holes above 100 solar masses within roughly 0.1 AU of these binaries.
-
Electromagnetic Precursors to Binary Neutron Star Mergers: Kinetic Simulations of Magnetospheric Flaring
3D kinetic simulations of pre-merger binary neutron star magnetospheres predict nonthermal gamma-ray signals at ~16 MeV and fast radio burst-like radio transients from reconnecting current sheets.
-
Numerical Studies of Accretion Flows onto a Neutron Star Engulfed in a Massive Star
Simulations find nested bow shocks around an engulfed neutron star enhance drag force 10-100 times over Bondi-Hoyle-Lyttleton and allow force reversal based on envelope conditions.
-
Highly eccentric non-spinning binary black hole mergers: quadrupolar post-merger waveforms
Polynomial models for the (2,2) post-merger waveform amplitudes of eccentric non-spinning binary black holes are constructed from numerical-relativity data as functions of symmetric mass ratio and two merger-time dyna...
-
Radio sirens: inferring $H_0$ with binary black holes and neutral hydrogen in the era of the Einstein Telescope and the SKA Observatory
Using simulated binary black hole mergers and neutral hydrogen maps, the radio sirens method constrains H0 to 8% precision with 3000 high-SNR events, offering a 90% improvement over standard dark siren analyses.
-
Prospects for multi-messenger discovery of the gravitational-wave background anisotropies via cross-correlation with galaxies
New simulations show that cross-correlating gravitational wave background anisotropies with galaxy distributions can enable discovery at angular scales of 4-6 degrees with next-generation observatories.
-
BB plot: A Tool for Accurate Model Selection Using Bayes factors
The BB plot is a new diagnostic that links Bayes factors to their expected distributions under competing hypotheses, enabling validation of calculations and low-cost background estimation for model selection in GW studies.
-
High-Spin BBH Subpopulation from AGN Accretion
Mixture model analysis of LIGO data identifies a ~10% high-spin subpopulation with a1 ≈ 0.9 matching AGN accretion predictions, disfavoring hierarchical mergers at a1 ≈ 0.7 for that group.
-
The impact of envelope binding energies on the merger rate density of binary compact objects
New envelope binding energies derived from PARSEC stellar tracks alter compact binary merger rate densities by over an order of magnitude when used in population synthesis.
-
What You Don't Know Won't Hurt You: Self-Consistent Hierarchical Inference with Unknown Follow-up Selection Strategies
Hierarchical Bayesian inference allows accurate recovery of intrinsic astrophysical source populations even when follow-up selection is unknown and correlated with parameters of interest.
-
Implications of the LISA stochastic signal from eccentric stellar mass black hole binaries in vacuum
High initial eccentricities in stellar-mass black hole binaries produce a stochastic gravitational wave background distinguishable by LISA from quasi-circular models, enabling upper bounds on eccentricity and separati...
-
Merger remnant and eccentricity dynamics surrogates for eccentric nonspinning black hole binaries
New surrogate models predict remnant properties and eccentricity dynamics for eccentric nonspinning black hole binary mergers using numerical relativity data over a limited parameter space.
-
Interpretable Analytic Formulae for GWTC-4 Binary Black Hole Population Properties via Symbolic Regression
Symbolic regression on GWTC-4 posteriors yields closed-form analytic formulae for merger-rate evolution, effective-spin dependencies on mass ratio and redshift, and conditional mass-ratio distributions at specific pri...
-
All-order structure of static gravitational interactions and the seventh post-Newtonian potential
A closed formula computes static post-Newtonian corrections at arbitrary odd orders in gravity, yielding the explicit seventh post-Newtonian potential that matches an independent diagrammatic method.
-
Post-Newtonian inspiral waveform model for eccentric precessing binaries with higher-order modes and matter effects
pyEFPEHM extends prior PN models to include higher-order quasi-circular phasing, generalized precession solutions, and eccentric corrections up to 1PN in selected multipoles for eccentric precessing binaries with matt...
-
Second-Generation Mass Peak in the Gravitational-Wave Population as a Probe of Globular Clusters
Dynamical formation in globular clusters produces a robust second black-hole mass peak at ~70 solar masses from second-generation mergers when the first-generation spectrum is truncated by pair-instability supernovae.
-
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...
-
A unified harmonic framework for dark siren cosmology
The GW-galaxy cross-correlation method, unified with spectral sirens in a harmonic framework, can measure H0 to 1% and Omega_m to 5% precision with 2 years of data from next-generation detectors like Einstein Telescop...
-
Electromagnetic Follow-up of the Sub-Solar Mass Gravitational Wave Candidate S251112cm: Kilonova Constraints and a Coincident IIb Supernova
No kilonova detected from sub-solar GW candidate S251112cm, but coincident IIb supernova SN 2025adtq yields suggestive evidence for the superkilonova channel, though inconclusive after accounting for chance coincidence.
-
How do the LIGO-Virgo-KAGRA's Heavy Black Holes Form? No evidence for core-collapse Intermediate-mass black holes in GWTC-4
No evidence for core-collapse formed low-spin IMBHs in GWTC-4, with 90% upper limit on merger rate of 0.077 Gpc^{-3} yr^{-1}, low-spin BH mass truncation at 65 solar masses consistent with pair-instability gap lower e...
-
How do the LIGO-Virgo-KAGRA's Heavy Black Holes Form? No evidence for core-collapse Intermediate-mass black holes in GWTC-4
No evidence for core-collapse IMBHs in GWTC-4; heavy BHs from hierarchical mergers, with low-spin mass distribution truncating at ~65 solar masses and PIMG upper edge estimated at 150 solar masses.
-
The Impact of Spin Priors on Parameterized Tests of General Relativity
Spin prior choices propagate into tests of GR via the 1.5PN deviation parameter δφ̂3 in a non-trivial, event-dependent way, with stronger effects for short-inspiral events and partial degeneracy with χ_eff when the de...
-
Blinded Mock Data Challenge: Is the Spectral Siren Technique Robust for Measuring the Hubble Constant?
The spectral siren technique requires independent knowledge of binary black hole mass distributions at all redshifts to an accuracy better than the statistical uncertainty, or else the Hubble constant measurement is biased.
-
Hawking area law in quantum gravity
Exact Hawking area law from black hole mergers restricts quantum gravity to singular Ricci-flat or specific regular black holes in Stelle and nonlocal theories, derives the standard entropy-area law, and realizes Barr...
-
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...
-
GW250114: testing Hawking's area law and the Kerr nature of black holes
GW250114 data confirm the remnant black hole ringdown frequencies lie within 30% of Kerr predictions and that the final horizon area is larger than the sum of the progenitors' areas to high credibility.
-
Constraining the Galactic Center Dark Cluster with ELT/MICADO Observations
ELT/MICADO is projected to enable the first quantitative constraints on the dark cluster of stellar compact objects in the Galactic Center through variability and astrometric searches.
-
Impact of coalescence signals on the search for continuous gravitational waves with Einstein Telescope
Unresolved compact binary coalescence background acts as additional noise that worsens Frequency-Hough sensitivity for continuous gravitational waves by 7-10% near 7 Hz under Einstein Telescope conditions.
-
Searches for Binary Mergers with Sub-solar Mass Components in Data from the First Part of LIGO--Virgo--KAGRA's Fourth Observing Run
No sub-solar mass binary merger candidates found in LIGO data from May 2023 to January 2024, yielding merger rate upper limits of 110-10000 Gpc^{-3}yr^{-1} and constraints on primordial black hole dark matter fractions.
-
No model-independent evidence for a peak in binary black hole spin (mis)alignments
No model-independent evidence for a peak in binary black hole spin tilts is found in GWTC-4; mass-spin magnitude correlation is confirmed but mass-tilt correlation is not.
-
ENGRAVE follow-up of a type IIb supernova spatially coincident with the sub-threshold gravitational wave trigger S250818k
SN2025ulz is a type IIb supernova whose shock-cooling tail mimicked a kilonova, demonstrating a key contaminant for gravitational-wave counterpart searches.
-
Mitigating Systematic Errors in Parameter Estimation of Binary Black Hole Mergers in O1-O3 LIGO-Virgo Data
Parametric models incorporating waveform phase and amplitude uncertainties mitigate systematic errors in gravitational wave parameter estimation, producing consistent results across models and raw/deglitched data for ...
-
Gravitational-wave astronomy requires population-informed parameter estimation
Population-informed hierarchical parameter estimation is required for unbiased astrophysical interpretation of gravitational-wave events rather than using standard individual posteriors with reference priors.
-
VIGILant: an automatic classification pipeline for glitches in the Virgo detector
VIGILant applies tree-based models and a ResNet CNN to classify Virgo O3b glitches with 98% accuracy and has been deployed for daily use with an interactive dashboard.
-
Investigating the effect of sensitivity of KAGRA on sky localization of gravitational-wave sources from compact binary coalescences
KAGRA enhances sky localization of binary neutron star mergers in the LVK network via added baselines, with measurable gains at current sensitivity and larger improvements as range reaches ~30 Mpc.
-
Not too close! Evaluating the impact of the baseline on the localization of binary black holes by next-generation gravitational-wave detectors
Baselines of 8-11 ms light travel time for two CE detectors provide a reasonable compromise for BBH sky localization, with third detectors eliminating multimodality for most or all events.
-
Chasing Gamma-Ray Signals from Binary Neutron Star Coalescences with the Cherenkov Telescope Array: Prospects and Observing Strategies
Simulations indicate that an optimized CTAO strategy could detect GeV-TeV emission from about 5% of GW-associated short GRBs, with detectability depending strongly on jet and viewing angles.
-
Neutron Star Merger Rates from Multi-messenger Observations: Clues to the Physical Origin of the Short and Long-short Gamma-ray Bursts
Short and long-short GRB rates of 195-666 Gpc^{-3} yr^{-1} exceed the 66-347 Gpc^{-3} yr^{-1} EM-bright NS merger rate from GW observations, challenging a pure NS merger origin.
-
Remnant recoil and host environments of GWTC-4.0 binary black-hole mergers
Five O4 gravitational-wave events prefer dynamical cluster formation; typical recoil kicks eject remnants from globular clusters but allow possible retention in nuclear star clusters, disfavoring efficient hierarchica...
-
Implications of low neutron star merger rates for gamma-ray bursts, r-process production and Galactic double neutron stars
Lower BNS merger rates from GWTC-4 data produce tensions of factors 3.6-18 with SGRB rates, 0.9-4.1 with r-process rates, and 2.3-5.1 with Galactic DNS rates.
-
Wide Jets or Low Rates: Reconciling Short GRB and Gravitational-Wave Neutron Star Merger Rates
Latest GW neutron star merger rates are consistent with short GRBs being produced by BNS mergers if jets are wide or rates low, with NSBH mergers subdominant.
-
Gravitational-wave standard sirens and application in cosmology
Gravitational-wave events from binary mergers act as standard sirens to independently measure luminosity distances and constrain cosmological parameters such as the Hubble constant and dark energy.
- Signatures of a subpopulation of hierarchical mergers in the GWTC-4 gravitational-wave dataset
- Gravitational-wave constraints on the pair-instability mass gap and nuclear burning in massive stars
Reference graph
Works this paper leans on
-
[1]
Aasi, J., et al. 2015, Class. Quant. Grav., 32, 074001, doi: 10.1088/0264-9381/32/7/074001 Abac, A. G., et al. 2024, Astrophys. J. Lett., 970, L34, doi: 10.3847/2041-8213/ad5beb —. 2025a, To be published in this issue. https://arxiv.org/abs/2508.18080 —. 2025b, To be published in this issue. https://arxiv.org/abs/2508.18081 —. 2025c, To be published in th...
-
[2]
The Mass distribution of stellar black holes.Astrophys
https://arxiv.org/abs/2412.03461 Bailyn, C. D., Jain, R. K., Coppi, P., & Orosz, J. A. 1998, Astrophys. J., 499, 367, doi: 10.1086/305614 Banagiri, S., Callister, T. A., Doctor, Z., & Kalogera, V
-
[3]
https://arxiv.org/abs/2501.06712 Banerjee, S. 2021, Mon. Not. Roy. Astron. Soc., 503, 3371, doi: 10.1093/mnras/stab591 Banerjee, S., Olejak, A., & Belczynski, K. 2023, Astrophys. J., 953, 80, doi: 10.3847/1538-4357/acdd59 Barr, E. D., et al. 2024, Science, 383, 275, doi: 10.1126/science.adg3005 Bartos, I., Kocsis, B., Haiman, Z., & M´ arka, S. 2017, Astro...
-
[4]
https://arxiv.org/abs/2010.05082 Capote, E., et al. 2025, Phys. Rev. D, 111, 062002, doi: 10.1103/PhysRevD.111.062002 Cheng, A. Q., Zevin, M., & Vitale, S. 2023, Astrophys. J., 955, 127, doi: 10.3847/1538-4357/aced98 Chru´ sli´ nska, M. 2024, Annalen Phys., 536, 2200170, doi: 10.1002/andp.202200170 Clesse, S., & Garcia-Bellido, J. 2022, Phys. Dark Univ., ...
-
[5]
https://arxiv.org/abs/2107.13082 Di Carlo, U. N., Giacobbo, N., Mapelli, M., et al. 2019, Mon. Not. Roy. Astron. Soc., 487, 2947, doi: 10.1093/mnras/stz1453 65 Di Carlo, U. N., Mapelli, M., Bouffanais, Y., et al. 2020, Mon. Not. Roy. Astron. Soc., 497, 1043, doi: 10.1093/mnras/staa1997 Dietrich, T., Coughlin, M. W., Pang, P. T. H., et al. 2020, Science, 3...
-
[6]
Ensuring Consistency between Noise and Detection in Hierarchical Bayesian Inference
https://arxiv.org/abs/2204.00461 Essick, R., & Fishbach, M. 2024, Astrophys. J., 962, 169, doi: 10.3847/1538-4357/ad1604 Essick, R., et al
-
[7]
https://arxiv.org/abs/2508.10638 Ewing, B., et al. 2024, Phys. Rev. D, 109, 042008, doi: 10.1103/PhysRevD.109.042008 Fairhurst, S., Green, R., Hoy, C., Hannam, M., & Muir, A. 2020, Phys. Rev. D, 102, 024055, doi: 10.1103/PhysRevD.102.024055 Farah, A., Ezquiaga, J. M., Fishbach, M., & Holz, D. 2025a. https://arxiv.org/abs/2507.07964 Farah, A. M., Callister...
-
[8]
https://arxiv.org/abs/2304.01288 Golomb, J., & Talbot, C. 2022, Astrophys. J., 926, 79, doi: 10.3847/1538-4357/ac43bc —. 2023, Phys. Rev. D, 108, 103009, doi: 10.1103/PhysRevD.108.103009 Guo, W.-H., Li, Y.-J., Wang, Y.-Z., et al. 2024, Astrophys. J., 975, 54, doi: 10.3847/1538-4357/ad758a Gupta, A., Gerosa, D., Arun, K. G., et al. 2020, Phys. Rev. D, 101,...
-
[9]
1981, A&A, 99, 126 Huth, S., et al
https://arxiv.org/abs/2411.02252 Hut, P. 1981, A&A, 99, 126 Huth, S., et al. 2022, Nature, 606, 276, doi: 10.1038/s41586-022-04750-w Iwaya, M., Kobayashi, K., Morisaki, S., Hotokezaka, K., & Kinugawa, T. 2025, Phys. Rev. D, 111, 103046, doi: 10.1103/PhysRevD.111.103046 Janquart, J., et al. 2024, doi: 10.1093/mnras/staf049 Jayasinghe, T., et al. 2021, Mon....
-
[10]
https://arxiv.org/abs/2201.01096 Klimenko, S., Mohanty, S., Rakhmanov, M., & Mitselmakher, G. 2005, Phys. Rev. D, 72, 122002, doi: 10.1103/PhysRevD.72.122002 Klimenko, S., Yakushin, I., Mercer, A., & Mitselmakher, G. 2008, Class. Quant. Grav., 25, 114029, doi: 10.1088/0264-9381/25/11/114029 Klimenko, S., et al. 2016, Phys. Rev. D, 93, 042004, doi: 10.1103...
-
[11]
The spin magnitude of stellar-mass black holes evolves with the mass
https://arxiv.org/abs/1912.11554 Pierra, G., Mastrogiovanni, S., & Perri` es, S. 2024, Astron. Astrophys., 692, A80, doi: 10.1051/0004-6361/202452545 Pompili, L., et al. 2023, Phys. Rev. D, 108, 124035, doi: 10.1103/PhysRevD.108.124035 Portegies Zwart, S. F., & McMillan, S. 2000, Astrophys. J. Lett., 528, L17, doi: 10.1086/312422 Portegies Zwart, S. F., &...
-
[12]
https://arxiv.org/abs/astro-ph/9710347 Postnov, K. A., & Yungelson, L. R. 2014, Living Rev. Rel., 17, 3, doi: 10.12942/lrr-2014-3 Pratten, G., et al. 2021, Phys. Rev. D, 103, 104056, doi: 10.1103/PhysRevD.103.104056 Qin, Y., Fragos, T., Meynet, G., et al. 2018, Astron. Astrophys., 616, A28, doi: 10.1051/0004-6361/201832839 Qin, Y., Marchant, P., Fragos, T...
-
[13]
2023, ApJ, 957, 37, doi: 10.3847/1538-4357/acf452
https://arxiv.org/abs/2404.03166 69 Ray, A., Maga˜ na Hernandez, I., Mohite, S., Creighton, J., & Kapadia, S. 2023a, Astrophys. J., 957, 37, doi: 10.3847/1538-4357/acf452 Ray, A., et al. 2023b. https://arxiv.org/abs/2306.07190 Renzo, M., Farmer, R. J., Justham, S., et al. 2020, Mon. Not. Roy. Astron. Soc., 493, 4333, doi: 10.1093/mnras/staa549 Rezzolla, L...
-
[14]
https://arxiv.org/abs/2506.05929 Rodriguez, C. L., Amaro-Seoane, P., Chatterjee, S., & Rasio, F. A. 2018, Phys. Rev. Lett., 120, 151101, doi: 10.1103/PhysRevLett.120.151101 Rodriguez, C. L., Chatterjee, S., & Rasio, F. A. 2016a, Phys. Rev. D, 93, 084029, doi: 10.1103/PhysRevD.93.084029 Rodriguez, C. L., & Loeb, A. 2018, Astrophys. J. Lett., 866, L5, doi: ...
-
[15]
https://arxiv.org/abs/2507.01086 Ruiz, M., Shapiro, S. L., & Tsokaros, A. 2018, Phys. Rev. D, 97, 021501, doi: 10.1103/PhysRevD.97.021501 Rutherford, N., et al. 2024, Astrophys. J. Lett., 971, L19, doi: 10.3847/2041-8213/ad5f02 Sachdev, S., et al
-
[16]
2024, ApJ, 960, 65, doi: 10.3847/1538-4357/ad0ce6
https://arxiv.org/abs/1901.08580 Sadiq, J., Dent, T., & Gieles, M. 2024, Astrophys. J., 960, 65, doi: 10.3847/1538-4357/ad0ce6 Sadiq, J., Dent, T., & Lorenzo-Medina, A. 2025a. https://arxiv.org/abs/2502.06451 —. 2025b. https://arxiv.org/abs/2506.02250 Sadiq, J., Dent, T., & Wysocki, D. 2022, Phys. Rev. D, 105, 123014, doi: 10.1103/PhysRevD.105.123014 Safa...
-
[17]
Disentangling spinning and nonspinning binary black hole populations with spin sorting
https://arxiv.org/abs/2507.23663 Talbot, C., Farah, A., Galaudage, S., Golomb, J., & Tong, H. 2025a, J. Open Source Softw., 10, 7753, doi: 10.21105/joss.07753 Talbot, C., & Golomb, J. 2023, Mon. Not. Roy. Astron. Soc., 526, 3495, doi: 10.1093/mnras/stad2968 Talbot, C., & Thrane, E. 2017, Phys. Rev. D, 96, 023012, doi: 10.1103/PhysRevD.96.023012 70 —. 2022...
work page internal anchor Pith review Pith/arXiv arXiv doi:10.21105/joss.07753 2023
- [18]
-
[19]
https://arxiv.org/abs/2312.06631 Wang, Y.-H., McKernan, B., Ford, S., et al. 2021, Astrophys. J. Lett., 923, L23, doi: 10.3847/2041-8213/ac400a Wang, Y.-Z., Li, Y.-J., Vink, J. S., et al. 2022, Astrophys. J. Lett., 941, L39, doi: 10.3847/2041-8213/aca89f Wiktorowicz, G., Wyrzykowski, L., Chruslinska, M., et al. 2019, doi: 10.3847/1538-4357/ab45e6 Winch, E...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.