Accretion-Driven Evolution of Compact-Object Populations in Gas-Rich Environments and the Origin of Massive Gravitational-Wave Sources
Pith reviewed 2026-06-29 02:36 UTC · model grok-4.3
The pith
Gas accretion acts as mass-space transport that broadens compact-object distributions when the accretion rate scales steeper than linearly with mass.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using a continuity-equation framework, we demonstrate that population evolution is governed primarily by the mass dependence of the accretion rate, ˙m ∝ m^β. Accretion laws with β>1 naturally produce divergent evolution and generate extended high-mass tails, whereas β<1 leads to convergent evolution and compresses the population toward a narrower range of masses. We show that sustained gas accretion can substantially broaden compact-object mass distributions, populate the high-mass end of gravitational-wave catalogs, and alter the mass-ratio distribution of compact-object binaries. In particular, collective accretion within compact binaries drives their mass ratios toward unity.
What carries the argument
Continuity-equation framework in which accretion rate scaling ˙m ∝ m^β determines whether populations diverge into high-mass tails or converge to a narrow mass range.
If this is right
- Accretion laws with β>1 generate extended high-mass tails in compact-object populations.
- Collective accretion in binaries drives mass ratios toward unity.
- Massive gravitational-wave events such as GW231123 become reachable through sustained accretion.
- Gravitational-wave catalogs exhibit broader mass distributions when gas-rich environments are common.
Where Pith is reading between the lines
- The same transport process could reshape mass distributions in other accreting systems where β exceeds one.
- Environmental dependence of mass distributions could be tested by comparing field versus cluster populations.
- If β is measured from individual accretion events, the model predicts observable shifts in binary mass-ratio histograms over time.
Load-bearing premise
The population evolution is governed primarily by the mass dependence of the accretion rate without other processes such as dynamical interactions or mergers dominating in gas-rich environments.
What would settle it
A survey of compact-object masses in gas-rich galaxies that shows neither extended high-mass tails nor a preference for equal-mass binaries would falsify the predicted transport effect.
Figures
read the original abstract
The origin of the most massive gravitational-wave sources remains elusive. We show that gas accretion can be understood as a transport process in mass space, causing compact objects to migrate through a population at rates determined by the underlying growth law. Using a continuity-equation framework, we demonstrate that population evolution is governed primarily by the mass dependence of the accretion rate, $\dot m \propto m^\beta$. Accretion laws with $\beta>1$ naturally produce divergent evolution and generate extended high-mass tails, whereas $\beta<1$ leads to convergent evolution and compresses the population toward a narrower range of masses. We apply this framework to physically motivated accretion regimes and explore their consequences using analytical calculations and Monte Carlo population models. We show that sustained gas accretion can substantially broaden compact-object mass distributions, populate the high-mass end of gravitational-wave catalogs, and alter the mass-ratio distribution of compact-object binaries. In particular, collective accretion within compact binaries drives their mass ratios toward unity. Our results suggest that gaseous environments act as transport media that continuously reshape compact-object populations, providing a natural pathway toward the formation of massive mergers such as GW231123 and the high-mass tails increasingly revealed by gravitational-wave observations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript applies a continuity-equation framework to compact-object population evolution under gas accretion, with accretion rate scaling as \dot{m} o m^eta. It claims that eta > 1 produces divergent evolution and high-mass tails while eta < 1 produces convergent evolution; sustained accretion thereby broadens mass distributions, populates the high-mass end of gravitational-wave catalogs, and drives binary mass ratios toward unity via collective accretion within binaries.
Significance. If the central results hold, the work supplies a general transport-process description of how gaseous environments can reshape compact-object mass and mass-ratio distributions, offering a pathway to the high-mass gravitational-wave events (e.g., GW231123) that is complementary to purely dynamical channels.
major comments (2)
- [Abstract / framework description] Abstract and framework description: the claim that population evolution is governed primarily by the mass dependence of the accretion rate requires that accretion timescales are shorter than those of dynamical interactions and mergers, yet no quantitative comparison of these timescales is supplied for the gas-rich environments considered.
- [Abstract] Abstract: conclusions are stated to follow from analytical calculations and Monte Carlo population models, but the manuscript supplies neither explicit derivations of the continuity-equation solutions nor error analysis or direct comparison against observed catalogs, preventing verification that the math supports the stated claims.
minor comments (1)
- The exponent eta is introduced as a free parameter that controls the evolutionary outcome; a clearer mapping from specific physical accretion regimes (e.g., Bondi, Eddington-limited) to the adopted values of eta would strengthen the physical motivation.
Simulated Author's Rebuttal
We thank the referee for these constructive comments, which identify areas where the manuscript can be strengthened for clarity and verifiability. We respond to each major comment below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Abstract / framework description] Abstract and framework description: the claim that population evolution is governed primarily by the mass dependence of the accretion rate requires that accretion timescales are shorter than those of dynamical interactions and mergers, yet no quantitative comparison of these timescales is supplied for the gas-rich environments considered.
Authors: We agree that a quantitative timescale comparison is necessary to support the assumption that accretion dominates the evolution. In the revised manuscript we will add a new subsection (in the framework section) with order-of-magnitude estimates for accretion timescales in representative gas-rich environments (AGN disks, nuclear star clusters) and direct comparisons to literature values for dynamical friction, binary hardening, and merger timescales. This will delineate the parameter regime in which the continuity-equation description applies. revision: yes
-
Referee: [Abstract] Abstract: conclusions are stated to follow from analytical calculations and Monte Carlo population models, but the manuscript supplies neither explicit derivations of the continuity-equation solutions nor error analysis or direct comparison against observed catalogs, preventing verification that the math supports the stated claims.
Authors: The current text describes the continuity-equation approach and Monte Carlo results but does not provide step-by-step derivations or catalog comparisons. We will add an appendix containing the explicit analytical solutions to the continuity equation under power-law accretion, include sensitivity/error analysis for the Monte Carlo runs (e.g., variation with initial mass function and accretion duration), and insert a figure/table comparing the predicted high-mass tail and mass-ratio distribution to events in GWTC-3, with particular reference to GW231123. revision: yes
Circularity Check
No circularity: continuity-equation outcomes follow directly from assumed accretion law without reduction to fits or self-citations
full rationale
The paper applies a standard continuity-equation transport model in mass space to accretion with power-law dependence ˙m ∝ m^β. Different β values produce divergent or convergent evolution by direct integration of the continuity equation; these are explored via analytics and Monte Carlo sampling. No quoted step shows a parameter fitted to the target distribution and then relabeled as a prediction, nor any load-bearing self-citation or ansatz smuggled from prior work. The framework is self-contained against external benchmarks once the growth law and dominance assumption are stated; outcomes are model consequences rather than tautological redefinitions.
Axiom & Free-Parameter Ledger
free parameters (1)
- β (accretion rate exponent)
axioms (1)
- domain assumption Compact-object population evolution obeys a continuity equation in mass space driven primarily by the accretion rate's mass dependence
Reference graph
Works this paper leans on
-
[1]
Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, PhRvL, 116, 061102, doi: 10.1103/PhysRevLett.116.061102
-
[2]
1988, ApJ, 332, 646, doi: 10.1086/166683
Szuszkiewicz, E. 1988, ApJ, 332, 646, doi: 10.1086/166683
-
[3]
2019, ApJ, 884, 22, doi: 10.3847/1538-4357/ab3466
Antoni, A., MacLeod, M., & Ramirez-Ruiz, E. 2019, ApJ, 884, 22, doi: 10.3847/1538-4357/ab3466
-
[4]
Antonini, F., Romero-Shaw, I. M., & Callister, T. 2025, Physical Review Letters, 134, 011401, doi: 10.1103/PhysRevLett.134.011401
-
[5]
Bardeen, J. M. 1970, Nature, 226, 64, doi: 10.1038/226064a0
-
[6]
2025, arXiv e-prints, arXiv:2508.08558, doi: 10.48550/arXiv.2508.08558
Bartos, I., & Haiman, Z. 2025, arXiv e-prints, arXiv:2508.08558, doi: 10.48550/arXiv.2508.08558
-
[7]
2017, ApJ, 835, 165, doi: 10.3847/1538-4357/835/2/165
Bartos, I., Kocsis, B., Haiman, Z., & M´ arka, S. 2017, ApJ, 835, 165, doi: 10.3847/1538-4357/835/2/165
-
[8]
2018, ARA&A, 56, 83, doi: 10.1146/annurev-astro-081817-051839
Bastian, N., & Lardo, C. 2018, ARA&A, 56, 83, doi: 10.1146/annurev-astro-081817-051839
-
[9]
Blandford, R. D., & Znajek, R. L. 1977, Monthly Notices of the Royal Astronomical Society, 179, 433, doi: 10.1093/mnras/179.3.433
-
[10]
1952, MNRAS, 112, 195, doi: 10.1093/mnras/112.2.195
Bondi, H. 1952, MNRAS, 112, 195, doi: 10.1093/mnras/112.2.195
-
[11]
2001, MNRAS, 322, 231, doi: 10.1046/j.1365-8711.2001.04022.x
Bonnell, I. A., Bate, M. R., Clarke, C. J., & Pringle, J. E. 2001a, MNRAS, 323, 785, doi: 10.1046/j.1365-8711.2001.04270.x
-
[12]
2001, MNRAS, 322, 231, doi: 10.1046/j.1365-8711.2001.04022.x
Bonnell, I. A., Clarke, C. J., Bate, M. R., & Pringle, J. E. 2001b, MNRAS, 324, 573, doi: 10.1046/j.1365-8711.2001.04311.x
-
[13]
2019, MNRAS, 490, 5196, doi: 10.1093/mnras/stz2977
Rosotti, G. 2019, MNRAS, 490, 5196, doi: 10.1093/mnras/stz2977
-
[14]
Fishbach, M., Holz, D. E., & Farr, B. 2017, ApJL, 840, L24, doi: 10.3847/2041-8213/aa7045
-
[15]
2015, MNRAS, 451, 3941, doi: 10.1093/mnras/stv1214
Gerosa, D., Veronesi, B., Lodato, G., & Rosotti, G. 2015, MNRAS, 451, 3941, doi: 10.1093/mnras/stv1214
-
[16]
Second-Generation Mass Peak in the Gravitational-Wave Population as a Probe of Globular Clusters
Ginat, Y. B., Antonini, F., Flanagan, E., & Gieles, M. 2026, arXiv e-prints, arXiv:2604.07456, doi: 10.48550/arXiv.2604.07456 Gr¨ obner, M., Ishibashi, W., Tiwari, S., Haney, M., &
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.07456 2026
-
[17]
2020, A&A, 638, A119, doi: 10.1051/0004-6361/202037681
Jetzer, P. 2020, A&A, 638, A119, doi: 10.1051/0004-6361/202037681
-
[18]
Hoyle, F., & Lyttleton, R. A. 1941, MNRAS, 101, 227, doi: 10.1093/mnras/101.4.227
-
[19]
2019, ApJ, 876, 142, doi: 10.3847/1538-4357/ab158b
Kaaz, N., Antoni, A., & Ramirez-Ruiz, E. 2019, ApJ, 876, 142, doi: 10.3847/1538-4357/ab158b
-
[20]
2023, ApJ, 944, 44, doi: 10.3847/1538-4357/aca967 Kıro˘ glu, F., Kremer, K., & Rasio, F
Ramirez-Ruiz, E. 2023, ApJ, 944, 44, doi: 10.3847/1538-4357/aca967 Kıro˘ glu, F., Kremer, K., & Rasio, F. A. 2025, arXiv e-prints, arXiv:2509.05415, doi: 10.48550/arXiv.2509.05415
-
[21]
Lin, D. N. C., & Murray, S. D. 2007, ApJ, 661, 779, doi: 10.1086/515387
-
[22]
2023, MNRAS, 526, 4291, doi: 10.1093/mnras/stad2993
Natarajan, P. 2023, MNRAS, 526, 4291, doi: 10.1093/mnras/stad2993
-
[23]
McKernan, B., Ford, K. E. S., Lyra, W., & Perets, H. B. 2012, MNRAS, 425, 460, doi: 10.1111/j.1365-2966.2012.21486.x
-
[24]
McKernan, B., Ford, K. E. S., Bellovary, J., et al. 2018, ApJ, 866, 66, doi: 10.3847/1538-4357/aadae5 13
-
[25]
Curd, B. 2022, Monthly Notices of the Royal Astronomical Society, 511, 3795, doi: 10.1093/mnras/stac285 ¨Opik, E. 1924, Publications of the Tartu Astrofizica Observatory, 25, 1
-
[26]
A., Sukhbold, T., & Eldridge, J
Patton, R. A., Sukhbold, T., & Eldridge, J. J. 2022, MNRAS, 511, 903, doi: 10.1093/mnras/stab3797
-
[27]
Peters, P. C. 1964, Physical Review, 136, 1224, doi: 10.1103/PhysRev.136.B1224
-
[28]
N., Covarrubias, R., & Ramirez-Ruiz, E
Quinteros, K. N., Covarrubias, R., & Ramirez-Ruiz, E. 2025, Nature Astronomy, 9, 1770, doi: 10.1038/s41550-025-02736-y
-
[29]
Rakavy, G., & Shaviv, G. 1967, The Astrophysical Journal, 148, 803, doi: 10.1086/149204
-
[30]
2026, ApJ, 997, 13, doi: 10.3847/1538-4357/ae28c6
Rosselli-Calderon, A., Stewart, J., Shen, S., et al. 2026, ApJ, 997, 13, doi: 10.3847/1538-4357/ae28c6
-
[31]
2025, A&A, 702, A208, doi: 10.1051/0004-6361/202556434
Roupas, Z. 2025, A&A, 702, A208, doi: 10.1051/0004-6361/202556434
-
[32]
Rozner, M., & Perets, H. B. 2022, ApJ, 931, 149, doi: 10.3847/1538-4357/ac6d55
-
[33]
Rozner, M., & Perets, H. B. 2024, ApJ, 968, 80, doi: 10.3847/1538-4357/ad4bdd
-
[34]
Salpeter, E. E. 1955, ApJ, 121, 161, doi: 10.1086/145971
-
[35]
Stone, N. C., Metzger, B. D., & Haiman, Z. 2017, MNRAS, 464, 946, doi: 10.1093/mnras/stw2260
-
[36]
2025, MNRAS, 543, 1864, doi: 10.1093/mnras/staf1592
Su, Y., Rowan, C., & Rozner, M. 2025, MNRAS, 543, 1864, doi: 10.1093/mnras/staf1592
-
[37]
Janka, H. T. 2016, ApJ, 821, 38, doi: 10.3847/0004-637X/821/1/38
work page internal anchor Pith review doi:10.3847/0004-637x/821/1/38 2016
-
[38]
2020, ApJ, 898, 25, doi: 10.3847/1538-4357/ab9b8c
Tagawa, H., Haiman, Z., & Kocsis, B. 2020, ApJ, 898, 25, doi: 10.3847/1538-4357/ab9b8c
-
[39]
Properties of black hole mergers in disks of active galactic nuclei
Tagawa, H., Haiman, Z., & Kocsis, B. 2026, arXiv e-prints, arXiv:2604.25994, doi: 10.48550/arXiv.2604.25994
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.25994 2026
-
[40]
2021, ApJ, 908, 194, doi: 10.3847/1538-4357/abd555
Tagawa, H., Kocsis, B., Haiman, Z., et al. 2021, ApJ, 908, 194, doi: 10.3847/1538-4357/abd555
-
[41]
Tchekhovskoy, A., McKinney, J. C., & Narayan, R. 2012, Journal of Physics: Conference Series, 372, 012040, doi: 10.1088/1742-6596/372/1/012040
-
[42]
Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, Monthly Notices of the Royal Astronomical Society: Letters, 418, L79, doi: 10.1111/j.1745-3933.2011.01147.x The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, et al. 2025, arXiv e-prints, arXiv:2507.08219, doi: 10.48550/arXiv.2507.08219
-
[43]
Thorne, K. S. 1974, Astrophysical Journal, 191, 507, doi: 10.1086/152991
-
[44]
Woosley, S. E., & Heger, A. 2021, ApJL, 912, L31, doi: 10.3847/2041-8213/abf2c4
-
[45]
2019a, ApJ, 876, 122, doi: 10.3847/1538-4357/ab16e3
Yang, Y., Bartos, I., Haiman, Z., et al. 2019a, ApJ, 876, 122, doi: 10.3847/1538-4357/ab16e3
-
[46]
2019b, PhRvL, 123, 181101, doi: 10.1103/PhysRevLett.123.181101
Yang, Y., Bartos, I., Gayathri, V., et al. 2019b, PhRvL, 123, 181101, doi: 10.1103/PhysRevLett.123.181101
-
[47]
Young, M. D., & Clarke, C. J. 2015, MNRAS, 452, 3085, doi: 10.1093/mnras/stv1512
-
[48]
1982, Annals of the New York Academy of Sciences, 395, 226, doi: 10.1111/j.1749-6632.1982.tb43399.x
Zinnecker, H. 1982, Annals of the New York Academy of Sciences, 395, 226, doi: 10.1111/j.1749-6632.1982.tb43399.x
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.