REVIEW
The Mass-Ratio Distribution of the Low-Mass Binary Black Hole Subpopulation: A Natural Outcome of Isolated Binary Evolution
High Energy Astrophysical Phenomena · Submitted 1 Sep 2026 · 184 references
Pith's verdict: A well-motivated forward model from POSYDON showing a two-component q distribution that is robust across model variations, but the observational comparison is qualitative: only a selected subset of models is shown, and no quantitative consistency measure is provided.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The mass-ratio distribution of low-mass black hole mergers emerges naturally from isolated binary evolution across a range of model assumptions, as shown by detailed population synthesis simulations.
The reading
What carries the argument
The key machinery is the POSYDON binary population synthesis framework, which couples grids of detailed 1D stellar evolution simulations (MESA) with a binary evolution module that self-consistently handles Roche-lobe overflow, common-envelope evolution, stellar winds, supernova kicks, and remnant formation. By sampling initial binary parameters (masses, separations, metallicities) from observationally motivated distributions and weighting by cosmic star formation history (IllustrisTNG-100), POSYDON produces a physically motivated population of merging BBHs. The central role of POSYDON is to replace rapid, analytic prescriptions with full stellar-structure calculations for the pre-common-envelope and mass-transfer phases, reducing the number of free parameters (e.g., the binding energy parameter $\lambda$ is computed directly from density profiles). The paper uses 15 model configurations varying supernova engines (Sukhbold16-engineN20, Sukhbold16-engineW20, PattonSukhbold20-engineW20), common-envelope efficiency ($\alpha_{\rm CE}=0.5,1.0,2.0$), and black hole accretion efficiency (Eddington-limited vs. a GRRMHD-informed prescription).
Extended reading notes
Core claim
The central discovery is that the mass-ratio distribution of merging binary black holes formed through isolated binary evolution (as simulated with the POSYDON code) consistently displays a two-component structure: a broad peak at asymmetric mass ratios ($q\simeq0.5$–$0.7$) and a secondary peak or plateau at near-equal masses ($q\simeq1$). This structure persists across a wide range of model assumptions—three different supernova remnant prescriptions, two common-envelope efficiencies ($\alpha_{\rm CE}=1.0,2.0$), and both Eddington-limited and GRRMHD-informed black hole accretion efficiencies. The asymmetric component originates from the combined contributions of standard common-envelope evolution and stable mass transfer, while the near-equal-mass component traces double-core common-envelope and contact mass-transfer systems. When compared to the $q$ distribution inferred for the low-mass ($m_1\lesssim15\,M_\odot$) subpopulation in GWTC-5.0, the theoretical predictions fall within the 90% credible intervals of both a non-parametric binned Gaussian process reconstruction and a more parametric mixture model. This agreement suggests that the observed low-mass BBH subpopulation is substantially composed of systems formed through isolated binary evolution, without the need for exotic physics or fine-tuning.
Load-bearing premise
The observational q distribution used for comparison comes from a parametric model that was deliberately constructed to describe systems with the spin and mass properties expected from isolated binary evolution, so the agreement between theory and observation may be partly built into the inference rather than being an independent test.
Editorial extensions
If this is right
- If the paper is correct, the distinct two-component $q$ structure (asymmetric peak + near-equal-mass tail) should become more sharply defined as the gravitational-wave catalog grows, providing a clean observational signature to distinguish isolated binary evolution from dynamical or primordial formation channels.
- The relative strength of the asymmetric versus near-equal-mass components will depend on common-envelope efficiency and accretion physics, meaning that future measurements of the $q$ distribution can constrain these poorly known stellar and binary parameters.
- The paper predicts that the low-mass subpopulation ($m_1\lesssim15\,M_\odot$) should show a combined $q$ distribution that is inconsistent with dynamical formation predictions (which favor near-equal-mass systems) but consistent with the isolated-binary shape, supporting multi-channel population models.
- As improved models of hydrogen-rich mass-transfer efficiency become available (e.g., Xing et al. 2026), the predicted asymmetric peak is expected to shift toward higher $q$, potentially improving agreement with observational inferences that show a plateau near $q=1$.
- The identification of double-core common-envelope and contact systems as the origin of the near-equal-mass component offers a specific prediction: systems with $q\approx1$ should have experienced a phase where both stars filled their Roche lobes simultaneously, a testable feature through detailed binary modelling.
Reading between the lines
- The strong agreement between the theoretical $q$ shape and the parametric model of Ray et al. (2026) is partially circular: that model was constructed targeting systems with properties expected from isolated evolution (aligned spins, broad $q$). A future test would be to use an observational inference that does not assume any formation channel, and to quantify the Bayesian evidence for or against
- The paper does not explore variations in the natal kick prescription beyond a single Maxwellian with $\sigma=265$ km/s. The conclusion that SMT systems produce asymmetric $q$ depends sensitively on kick magnitude; models with no kicks produce near-equal-mass SMT systems. The paper's result may therefore depend on the adopted kick distribution, which is not well constrained.
- The merger rate predicted by the models is systematically higher than the observed rate (up to a factor of 37), which the paper acknowledges but does not resolve. If the rate problem is solved by reducing common-envelope efficiency or increasing accretion feedback, the relative balance between the two $q$ components could shift, possibly breaking the current consistency with the observed $q$ shape
Original abstract
The observed merging binary black hole population is increasingly consistent with being composed of a mixture of subpopulations, each likely the result of different formation mechanisms. In particular, the low-mass subpopulation, with primary black hole masses below $\approxeq 15\,M_{\odot}$, has been attributed to mergers formed through isolated binary evolution. We use the binary population synthesis code POSYDON to study the mass ratio $(q)$ distribution of binary black hole mergers from isolated binary evolution and compare to the observed low-mass subpopulation. We explore variations in supernova remnant prescriptions, common-envelope efficiency, and black hole accretion efficiency. We find that our models have a preference for asymmetric $q$, most with a broad peak near $q\approxeq0.5-0.7$, and a near-equal-mass component whose relative strength varies across models. This resultant $q$ distribution is consistent with the $q$ distribution of the low-mass subpopulation observed with gravitational waves. We find that for the majority of models these features arise from physically distinct formation subchannels: the asymmetric peak reflects contributions from common-envelope and stable mass-transfer systems, while the near-equal-mass component traces double-core common envelope and contact systems. We conclude that the features in the $q$ distribution of the observed low-mass subpopulation emerge naturally from isolated binary evolution across a range of model assumptions. As the gravitational-wave catalog continues to grow, the relative strength of these $q$ features will provide an increasingly powerful diagnostic of isolated binary evolution and its formation subchannels.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
The observed merging binary black hole population is increasingly consistent with being composed of a mixture of subpopulations, each likely the result of different formation mechanisms. In particular, the low-mass subpopulation, with primary black hole masses below $\approxeq 15\,M_{\odot}$, has been attributed to mergers formed through isolated binary evolution. We use the binary population synthesis code POSYDON to study the mass ratio $(q)$ distribution of binary black hole mergers from isolated binary evolution and compare to the observed low-mass subpopulation. We explore variations in supernova remnant prescriptions, common-envelope efficiency, and black hole accretion efficiency. We find that our models have a preference for asymmetric $q$, most with a broad peak near $q\approxeq0.5-0.7$, and a near-equal-mass component whose relative strength varies across models. This resultant $q$ distribution is consistent with the $q$ distribution of the low-mass subpopulation observed with gravitational waves. We find that for the majority of models these features arise from physically distinct formation subchannels: the asymmetric peak reflects contributions from common-envelope and stable mass-transfer systems, while the near-equal-mass component traces double-core common envelope and contact systems. We conclude that the features in the $q$ distribution of the observed low-mass subpopulation emerge naturally from isolated binary evolution across a range of model assumptions. As the gravitational-wave catalog continues to grow, the relative strength of these $q$ features will provide an increasingly powerful diagnostic of isolated binary evolution and its formation subchannels.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | The observed merging binary black hole population is increasingly consistent with being composed of a mixture of subpopulations, each likely the result of different formation mechanisms. | background 1.00 | 0.07 | 0.03 | hedged possibility 0.93 |
| s1 | In particular, the low-mass subpopulation, with primary black hole masses below $\approxeq 15\,M_{\odot}$, has been attributed to mergers formed through isolated binary evolution. | background 1.00 | 0.05 | 0.03 | hedged possibility 0.60 |
| s2 | We use the binary population synthesis code POSYDON to study the mass ratio $(q)$ distribution of binary black hole mergers from isolated binary evolution and compare to the observed low-mass subpopulation. | method 0.63 | 0.95 | 0.02 | asserted 0.96 |
| s3 | We explore variations in supernova remnant prescriptions, common-envelope efficiency, and black hole accretion efficiency. | method 0.92 | 0.87 | 0.01 | asserted 0.60 |
| s4 | We find that our models have a preference for asymmetric $q$, most with a broad peak near $q\approxeq0.5-0.7$, and a near-equal-mass component whose relative strength varies across models. | result 1.00 | 0.97 | 0.12 | asserted 0.98 |
| s5 | This resultant $q$ distribution is consistent with the $q$ distribution of the low-mass subpopulation observed with gravitational waves. | result 0.73 | 0.72 | 0.20 | hedged possibility 0.76 |
| s6 | We find that for the majority of models these features arise from physically distinct formation subchannels: the asymmetric peak reflects contributions from common-envelope and stable mass-transfer systems, while the near-equal-mass component traces double-core common envelope and contact systems. | interpretation 0.86 | 0.95 | 0.12 | asserted 1.00 |
| s7 | We conclude that the features in the $q$ distribution of the observed low-mass subpopulation emerge naturally from isolated binary evolution across a range of model assumptions. | interpretation 0.95 | 0.78 | 0.95 | asserted 0.99 |
| s8 | As the gravitational-wave catalog continues to grow, the relative strength of these $q$ features will provide an increasingly powerful diagnostic of isolated binary evolution and its formation subchannels. | future work 0.80 | 0.12 | 0.05 | asserted 0.50 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2609.01784 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.