REVIEW 4 major objections 4 minor 4 cited by
Supermassive black hole growth in hierarchically merging nuclear star clusters
T0 review · 4 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read The paper claims that supermassive black hole seeds formed in compact nuclear star clusters at high redshift, and that growing them along galaxy merger histories reproduces observed black hole–galaxy scaling relations.
desk verdict A transparent, code-backed plausibility argument for NSC-born SMBH seeds, with the main caveat being an unverified compact-radius assumption and some tuned parameters. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The carrying object is Nuce, a semianalytic code for nuclear cluster evolution that relies on Hénon's principle, the balanced evolution by which energy generation in the cluster core drives gradual expansion. Nuce evolves cluster structure and the central black hole through four growth channels: runaway stellar collisions, Bondi gas accretion capped at the Eddington limit, repeated mergers of stellar-mass black holes in a core-collapsed subsystem, and loss-cone consumption of stars. Coupled to galaxy merger trees from the NewHorizon simulation, the code decides which nuclear star clusters coalesce after major galaxy mergers and which black hole binaries merge, using published fitting formulae for remnant mass, spin, and gravitational-wave recoil; whether the merger remnant stays in the nucleus is set by comparing the kick velocity to the cluster escape velocity.
What would settle it
A high-redshift census of star-forming galaxies at redshifts 6 to 10 that finds nuclear star clusters are rarely as compact as about one parsec would undercut the model's seed channel; the compact clusters JWST has already found would need to be representative rather than rare exceptions.
Extended reading notes
Core claim
At the paper's center is the claim that galactic nuclei assemble their supermassive black holes hierarchically from light seeds. Each protogalaxy above a critical stellar mass forms a nuclear star cluster; inside compact clusters, stellar-mass black holes formed from massive stars undergo repeated mergers, with an escape-velocity threshold of roughly 300 kilometers per second, and are boosted by gas accretion to produce seeds of intermediate mass. When the seeding model is evolved along galaxy merger trees, the local population contains about 0.13 intermediate-mass black holes per cubic megaparsec, an occupation fraction near 80 percent in dwarf galaxies below one billion solar masses, and a black hole mass versus galaxy stellar mass relation with slope around 1.1. Supermassive black holes grow further during major galaxy mergers, those that grow by coherent gas accretion emerge with near-maximal spins, and gravitational-wave recoil occasionally ejects them into off-nuclear wandering or isolated states.
Load-bearing premise
The model assumes that almost every small protogalaxy above about ten million solar masses formed a dense nuclear star cluster at high redshift with sub-parsec scale sizes, and that such compact clusters are common enough; if they are not, the seed-assembly engine loses its fuel.
Editorial extensions
If this is right
- If the scenario is right, supermassive black hole seeds can form from ordinary stellar-mass black holes, removing the need for massive Population III or primordial seeds.
- The local universe should contain roughly 0.13 intermediate-mass black holes per cubic megaparsec, about half nuclear and half satellite, concentrated in dwarf galaxies.
- Space-based gravitational-wave observatories should see a massive black hole merger rate of about 5.3 events per year out to redshift 5, dominated by intermediate-mass black hole binaries.
- Nuclear star clusters would naturally produce a high tidal disruption event rate density of 10,000 to 100,000 per year per cubic gigaparsec at redshifts below 5.
- Nuclear supermassive black holes grown by coherent gas accretion should appear highly spinning, while intermediate-mass black holes should show a bimodal spin distribution depending on formation channel.
Reading between the lines
- A test the paper does not run is to compare its predicted near-maximal spins for supermassive black holes against X-ray reflection measurements with selection effects accounted for; if accretion is chaotic rather than coherent, spins would be lower and the growth channel would need revision.
- Because the simulated volume is small and the merger trees stop at redshift 0.25, the model says little about the rare, highly luminous quasars seen at redshift above 6; applying the same seeding prescription to larger-volume trees would show whether an additional massive-seed channel is needed.
- The predicted black hole occupation fraction of about 80 percent in dwarf galaxies is effectively an upper bound under the assumption that every galaxy above the mass threshold forms a nuclear cluster; patchy cluster formation would lower the implied number densities proportionally.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a semi-analytic model, NSC-tree, that couples the Nuce nuclear star cluster evolution code to galaxy merger trees extracted from the NewHorizon simulation. The model assumes that each protogalaxy above a critical stellar mass forms a nuclear star cluster with specified initial mass, gas fraction, and a log-uniform initial half-mass radius prior, and that stellar-mass black holes in these clusters grow through runaway collisions, repeated black-hole mergers, gas accretion, and stellar consumption to form massive black hole seeds. These seeds then grow and merge hierarchically during major galaxy mergers, producing predictions for the local MBH-Mstar relation, black hole spins, occupation fractions, number densities of IMBHs and SMBHs, NSC properties, and rates of massive black hole mergers, stellar-mass black hole captures, extreme-mass-ratio inspirals, and tidal disruption events. The headline claim is that the seeds of supermassive black holes formed in nuclear star clusters via stellar black hole mergers at early epochs, and that this scenario can reproduce local scaling relations and yields a LISA-detectable massive black hole merger rate.
Significance. If the central scenario is correct, the paper offers a concrete pathway from ordinary stellar-mass black holes to SMBH seeds without invoking exotic massive seed formation, and it produces a rich set of falsifiable predictions: local IMBH number densities, occupation fractions, TDE rates, EMRI rates, and LISA-detectable merger rates. The manuscript is transparent about its simplifications, provides a public code, and makes an explicit connection to a high-resolution cosmological simulation via merger trees. However, the model's headline outputs are not fully independent predictions: the essential initial conditions for high-redshift NSCs are unconstrained, and the authors state in the Conclusions that many free parameters were chosen to give a crude fit to local scaling laws. These issues must be addressed before the claim that the NSC paradigm explains SMBH seeding can be regarded as robust.
major comments (4)
- [§5 and §4.2] The log-uniform prior on the initial half-mass radius, rh,0 in [0.1, 10] pc, is load-bearing but observationally unconstrained at high redshift. The seed-formation thresholds in §2.3.1 (0.2 tau_rh,0 < 3 Myr for runaway collisions and v_esc > 300 km/s for repeated BH mergers) imply that only the compact tail, roughly rh,0 < 0.2-0.5 pc for N ~ 1e6-1e7, contributes to seed formation. Under the fiducial prior, this tail contains roughly 15-35% of the clusters, but if the true high-redshift NSC radius distribution is centered at a few pc with only a small compact tail, the seed-forming fraction drops by about an order of magnitude, and the predicted z=0.25 IMBH density of 0.13 Mpc^-3 falls below the 0.02 Mpc^-3 lower end of the observationally allowed range from Greene et al. (2020). The authors themselves flag this in §4.1: 'Compact NSCs on sub-pc scales could theoretically form... This assumption is essential for assembling SMBH seeds in the NSC paradigm.' I request a quantitative sensitivity test that reweights or replaces the radius prior with a distribution informed by local NSC radii and by the few available high-z cluster observations, and that reports how the IMBH density, occupation fraction, and merger rates change.
- [§2.2 and §3.1.3] The Conclusions state explicitly that 'Our code has many free parameters, that we have chosen to give a crude fit to local scaling laws involving NSCs.' Consequently, headline outputs such as the MBH-Mstar slope in Eq. (1), the occupation fractions in Fig. 6, and the IMBH number density in Table 2 are partly calibrated rather than independent predictions. The hyperparameter tests in §4.2 show that the slope of the MBH-Mstar relation varies from about 0.35 (fg,max = 0.1) to about 1.13 (fstar,max = 0.1), so the fiducial slope of 1.10 is not robust to plausible parameter choices. Please specify which hyperparameters were tuned, against which observables and with what tolerance, and then separate calibrated statements from genuinely predicted quantities. A small grid or emulator-based exploration of the parameter region that still matches the calibration targets would greatly strengthen the paper.
- [§2.4 and §3.3.1] The assumption that every protogalaxy forms one NSC in a single burst when its stellar mass reaches Mcr = 1e7 Msun directly produces the near-unity NSC and IMBH occupation fractions at stellar masses below about 1e8 Msun. The authors acknowledge in §3.1.3 that 'Our predictions for high occupation fractions are impacted by our choice to form an NSC in every galaxy with stellar mass larger than 1e7 Msun.' Because the occupation fraction is a key comparison with Greene et al. (2020) and Nguyen et al. (2019), this assumption needs to be relaxed or bracketed, for example by introducing an NSC formation efficiency with a physically motivated dependence on galaxy mass or environment, and by modeling the ex-situ contribution from inspiraling globular clusters, which is currently neglected.
- The model ignores time delays between galaxy mergers, NSC mergers, and massive BH binary mergers. The estimate in §4.1, tau_har + tau_gw ~ 44 Myr, addresses only the hardening and gravitational-wave phase and explicitly drops the dynamical-friction delay tau_df. Since the predicted LISA merger-rate redshift distribution in Fig. 10 and the detectability statement in §3.3.1 depend on when massive BH binaries actually coalesce, the instantaneous-merger approximation could bias the z-distribution of events. A simple log-normal delay prescription, or at least a comparison of the 44 Myr estimate with the median inter-merger time in the NewHorizon trees, would show whether this simplification is benign for the paper's LISA predictions.
minor comments (4)
- [§1] The phrase 'strong HHEII emission' contains a typo and should read 'strong He II emission.'
- [§3.2.1] The sentence 'Most of these NSCs have half-mass radii larger than 50M⊙' should read 'larger than 50 pc,' consistent with the detectability threshold defined in §3.2.
- [§4.1] The word 'super-Eddigton' should be 'super-Eddington.'
- [§2.4] The description of the major/minor merger classification is confusing: 'if the asymmetry in the galaxy masses is smaller' appears to describe a minor merger, but the sentence then states the NSC becomes an ultra-compact dwarf. Please rephrase to make the threshold Qth and the two outcomes unambiguous.
Circularity Check
Partially circular: several headline 'predictions' (MBH-Mstar relation, NSC/BH occupation fractions) are explicitly partly fitted to, or built into, the assumed NSC scaling laws and universal NSC formation, so the claimed consistency with observations is partly by construction.
-
fitted input called prediction
[Section 5 (Conclusions), paragraph beginning 'Our code has many free parameters']
"Our code has many free parameters, that we have chosen to give a crude fit to local scaling laws involving NSCs. We caution the reader that this is in the same spirit as other discussions of semi-analytical models of galaxy formation, e.g., Nadler et al. (2023)."
The paper presents the MBH-Mstar relation (Eq. 1), the MBH-MNSC correlation, the occupation fractions, and the local number densities as model predictions in Sec. 3, and then compares them with observed scaling relations. However, the model's hyperparameters were explicitly chosen to 'give a crude fit to local scaling laws involving NSCs.' Since those local NSC scaling laws overlap with the observational benchmarks used in Sec. 3.1-3.2, the agreement is partly the target of the fit rather than an independent confirmation. The paper also declines a comprehensive hyperparameter study in Sec. 4.2, so the robustness of the 'prediction' over the parameter space is not quantified.
-
self definitional
[Section 3.1.3 (Occupation fraction), final paragraph]
"Our predictions for high occupation fractions are impacted by our choice to form an NSC in every galaxy with stellar mass larger than 10^7 M_sun."
The NSC occupation fraction fNSC is presented as a model output and compared with Virgo observations, but by construction every galaxy with Mstar > Mcr = 10^7 M_sun hosts an NSC. Consequently, the high massive-BH occupation fractions in dwarfs (e.g., ~80% with MBH > 10^3 M_sun) largely follow from this input assumption, because galaxies without NSCs cannot host nuclear IMBHs in the model. The quoted sentence concedes that the high occupation fractions are 'impacted by our choice' to assume universal NSC formation; this makes the occupation-fraction 'prediction' partly self-definitional rather than an emergent test of the model.
full rationale
The core Nuce evolution plus NewHorizon merger-tree calculation is not itself circular: the seed-formation thresholds (runaway collisions requiring 0.2 tau_rh,0 < 3 Myr, repeated BH mergers requiring vesc > 300 km/s), the Bondi/Eddington accretion prescription, and the GW-recoil retention criterion are stated as explicit physical inputs, and the resulting growth curves in Fig. 2-3 are computed rather than assumed. The self-citations to Kritos et al. (2024a,b) for Nuce and spin evolution are not load-bearing in a circular sense here, because the present paper spells out the relevant equations and assumptions, and the code is public. The compact-NSC prior (log-uniform rh,0 over 0.1-10 pc) is admittedly essential for seed formation, but it is an openly flagged physical assumption rather than a derived prediction; the paper even labels it 'essential' in Sec. 4.1. The genuine circularity is narrower and located in the presentation of results: the paper states that its free parameters were chosen to give a crude fit to local NSC scaling laws, and that the high NSC occupation fraction follows from a deliberate assumption that every galaxy above 10^7 M_sun forms an NSC. Therefore, the agreement of the model's MBH-Mstar slope, occupation fractions, and NSC abundance with observed scaling relations is partly a consequence of those choices, not an independent confirmation of the seeding scenario. This is partial circularity in the headline 'prediction' claims, not a full reduction of the derivation to its inputs; the seed-formation mechanism itself retains independent dynamical content.
Assumptions & free parameters
free parameters (8)
- Mcr (critical stellar mass for NSC formation) =
1e7 solar masses
- Qth (major merger mass ratio threshold) =
0.1
- tau_ge (gas expulsion timescale) =
100 Myr
- fg,max (maximum initial gas fraction in NSC) =
1.0
- fstar,max (maximum initial NSC mass fraction of galaxy) =
1.0
- Initial half-mass radius prior =
log-uniform 0.1 to 10 pc
- Runaway collision seed mass fraction =
1e-3 of initial stellar mass
- Direct collapse seed spin =
0
assumptions (8)
- domain assumption Henon's principle and balanced evolution govern NSC core evolution.
- ad hoc to paper Each protogalaxy forms one NSC in a single burst when its stellar mass reaches Mcr.
- domain assumption Initial NSC gas is ionized hydrogen with cs = 10 km/s and is removed exponentially with timescale tau_ge.
- domain assumption Runaway stellar collisions form a seed BH of mass 1e-3 Mstar,0 if 0.2 tau_rh,0 < 3 Myr.
- domain assumption Gas accretion is Bondi-like, capped at Eddington, and coherently aligned with BH spin.
- domain assumption The final parsec problem is assumed solved, so BH binaries merge promptly after galaxy mergers.
- ad hoc to paper Massive BHs and NSCs do not feed back onto the galaxy merger tree.
- domain assumption Merger trees retain only galaxies above 1e7 solar masses and stop at z = 0.25.
Cite this review
Pith. "Pith review of Supermassive black hole growth in hierarchically merging nuclear star clusters." pith.science (2026). https://pith.science/paper/BQPPJDEI
@misc{pith2026241215334,
author = {Pith},
title = {Pith review of: Supermassive black hole growth in hierarchically merging nuclear star clusters},
year = {2026},
howpublished = {\url{https://pith.science/paper/BQPPJDEI}},
note = {Machine review of arXiv:2412.15334}
}
read the original abstract
Supermassive black holes are prevalent at the centers of massive galaxies, and their masses scale with galaxy properties, increasing evidence suggesting that these trends continue to low stellar masses. Seeds are needed for supermassive black holes, especially at the highest redshifts explored by the James Webb Space Telescope. We study the hierarchical merging of galaxies via cosmological merger trees and argue that the seeds of supermassive black holes formed in nuclear star clusters via stellar black hole mergers at early epochs. Observable tracers include intermediate-mass black holes, nuclear star clusters, and early gas accretion in host dwarf galaxies, along with a potentially detectable stochastic gravitational wave background, ejection of intermediate and supermassive black holes, and consequences of a significant population of tidal disruption events and extreme-mass ratio inspirals.
Figures
Figures from the paper (8 more)
Forward citations
Cited by 4 Pith papers
-
Spatially Offset Active Galactic Nuclei in the Very Large Array Sky Survey: Tracers of Galaxy Mergers and Wandering Massive Black Holes
A new catalog of 328 offset AGN candidates shows that the fraction of galaxies hosting an off-center accreting black hole is 30-70 times lower than expected, hinting that massive black hole seeding is inefficient.
-
The diffraction-lensing interpretation of GW231123 with astrophysical priors
GW231123 is better fit by a lower-mass merger diffracted by an isolated ~1000 M_sun lens, but astrophysical priors and a frequentist rate estimate make this lensing interpretation unlikely.
-
From nuclear star clusters to Little Red Dots: black hole growth, mergers, and tidal disruptions
A semi-analytic model shows dense early star clusters can grow 1000-solar-mass black hole seeds to ~10^7 solar masses within ~100 million years, predicting tens of tidal disruptions per year from Little Red Dots at z=4-6.
-
Progenitor of the recoiling super-massive black hole RBH-1 identified using HST/JWST imaging
The 954 km/s runaway speed of RBH-1 implies its progenitor was a precessing, nearly equal-mass binary with a rapidly spinning primary black hole.
Reference graph
Works this paper leans on
-
[1]
Abbott, R., et al. 2023, Phys. Rev. X, 13, 011048, doi: 10.1103/PhysRevX.13.011048 Adamo, A., Bradley, L. D., Vanzella, E., et al. 2024, Nature, 632, 513, doi: 10.1038/s41586-024-07703-7 Alexander, T., & Natarajan, P. 2014, Science, 345, 1330, doi: 10.1126/science.1251053 Amaro-Seoane, P., et al
-
[4]
https://arxiv.org/abs/2407.09949 Chabrier, G. 2005, Astrophys. Space Sci. Libr., 327, 41, doi: 10.1007/978-1-4020-3407-7 5 Dekel, A., Stone, N. C., Dutta Chowdhury, D., et al. 2024, . https://arxiv.org/abs/2409.18605 Devecchi, B., Volonteri, M., Colpi, M., & Haardt, F. 2010, Mon. Not. Roy. Astron. Soc., 409, 1057, doi: 10.1111/j.1365-2966.2010.17363.x Dev...
work page Pith review arXiv 2005
-
[5]
https://arxiv.org/abs/2409.06527 Li, K., Volonteri, M., Dubois, Y., Beckmann, R., & Trebitsch, M. 2024, . https://arxiv.org/abs/2410.07856 Liu, B., Gurian, J., Inayoshi, K., et al. 2024, Mon. Not. R. Astron. Soc., 534, 290, doi: 10.1093/mnras/stae2066 Lupi, A., Colpi, M., Devecchi, B., Galanti, G., & Volonteri, M. 2014, Mon. Not. Roy. Astron. Soc., 442, 3...
arXiv 2024
-
[7]
2013, Dynamics and Evolution of Galactic Nuclei Miller, B
https://arxiv.org/abs/2409.08326 Merritt, D. 2013, Dynamics and Evolution of Galactic Nuclei Miller, B. P., Gallo, E., Greene, J. E., et al. 2015, Astrophys. J., 799, 98, doi: 10.1088/0004-637X/799/1/98 Mouri, H., & Taniguchi, Y. 2002, Astrophys. J. Lett., 566, L17, doi: 10.1086/339472 Mowla, L., Iyer, K. G., Desprez, G., et al. 2022, Astrophys. J. Lett.,...
arXiv 2013
-
[9]
https://arxiv.org/abs/2409.18096 Pfister, H., Dai, J., Volonteri, M., et al. 2020, Mon. Not. Roy. Astron. Soc., 500, 3944, doi: 10.1093/mnras/staa3471 Polkas, M., Bonoli, S., Bortolas, E., et al. 2024, Astron. Astrophys., 689, A204, doi: 10.1051/0004-6361/202449470 Portegies Zwart, S. F., & McMillan, S. L. W. 2002, Astrophys. J., 576, 899, doi: 10.1086/34...
arXiv 2020
-
[10]
https://arxiv.org/abs/2410.06141 Quinlan, G. D. 1996, New Astron., 1, 35, doi: 10.1016/S1384-1076(96)00003-6 Quinlan, G. D., & Shapiro, S. L. 1987, Astrophys. J., 321, 199, doi: 10.1086/165624 Rantala, A., Naab, T., & Lah´ en, N
arXiv 1996
-
[11]
A., Freitag, M., & Atakan Gurkan, M
https://arxiv.org/abs/2403.10602 Rasio, F. A., Freitag, M., & Atakan Gurkan, M. 2003, in Carnegie Observatories Centennial Symposium
arXiv 2003
-
[12]
Formation of Massive Black Holes in Dense Star Clusters
Coevolution of Black Holes and Galaxies. https://arxiv.org/abs/astro-ph/0304038 Reinoso, B., Schleicher, D. R. G., Fellhauer, M., Klessen, R. S., & Boekholt, T. C. N. 2018, Astron. Astrophys., 614, A14, doi: 10.1051/0004-6361/201732224 Reynolds, C. S. 2021, Ann. Rev. Astron. Astrophys., 59, 117, doi: 10.1146/annurev-astro-112420-035022 Rizzuto, F. P., Naa...
work page Pith review arXiv 2018
Show all 14 references
-
[13]
S., & Dopita, M
https://arxiv.org/abs/2405.05333 Sutherland, R. S., & Dopita, M. A. 1993, Astrophys. J. Supp. S., 88, 253, doi: 10.1086/191823 Teyssier, R. 2002, Astron. Astrophys., 385, 337, doi: 10.1051/0004-6361:20011817 Trebitsch, M., Blaizot, J., Rosdahl, J., Devriendt, J., & Slyz, A. 20...
1993 arXiv
-
[14]
https://arxiv.org/abs/2411.03448
-
[145]
2021, Mon
https://arxiv.org/abs/1104.4797 —. 2021, Mon. Not. Roy. Astron. Soc., 501, 1413, doi: 10.1093/mnras/staa3724 Neumayer, N., Seth, A., & Boeker, T. 2020, Astron. Astrophys. Rev., 28, 4, doi: 10.1007/s00159-020-00125-0 Neumayer, N., Seth, A., & B¨ oker, T. 2020, Astron. Astrophys...
2021 arXiv
-
[2017]
2019, Mon
https://arxiv.org/abs/1702.00786 Antonini, F., Gieles, M., & Gualandris, A. 2019, Mon. Not. Roy. Astron. Soc., 486, 5008, doi: 10.1093/mnras/stz1149 Aubert, D., Pichon, C., & Colombi, S. 2004, Mon. Not. Roy. Astron. Soc., 352, 376, doi: 10.1111/j.1365-2966.2004.07883.x Babak, ...
2019 arXiv
-
[2023]
2024, Nature, 627, 59, doi: 10.1038/s41586-024-07052-5 Mandel, I., & Broekgaarden, F
https://arxiv.org/abs/2308.01230 SMBHs grow in hierarchically merging NSCs 19 —. 2024, Nature, 627, 59, doi: 10.1038/s41586-024-07052-5 Mandel, I., & Broekgaarden, F. S. 2022, Living Rev. Rel., 25, 1, doi: 10.1007/s41114-021-00034-3 Mapelli, M. 2016, Mon. Not. Roy. Astron. Soc...
2024 arXiv
-
[2024]
2008, Astrophys
https://arxiv.org/abs/2410.02875 Berti, E., & Volonteri, M. 2008, Astrophys. J., 684, 822, doi: 10.1086/590379 Bhowmick, A. K., Blecha, L., Torrey, P., et al. 2024a, Mon. Not. Roy. Astron. Soc., 529, 3768, doi: 10.1093/mnras/stae780 Bhowmick, A. K., et al. 2024b. https://arxiv...
2008 arXiv
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.