Recognition: 2 theorem links
· Lean TheoremSmoluchowski Coagulation Equation and the Evolution of Primordial Black Hole Clusters
Pith reviewed 2026-05-13 21:17 UTC · model grok-4.3
The pith
Simulations of the Smoluchowski coagulation equation show PBH clusters reach runaway merger on timescales set by initial number and density.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We solve the Smoluchowski coagulation equation for primordial black hole clusters by Monte Carlo methods after deriving a coagulation kernel that includes gravitational dynamics both with and without mass segregation. The simulations yield the runaway timescales of the clusters and the evolution of the PBH mass distribution across a wide range of cosmic redshifts, with both quantities depending on the number of PBHs inside the cluster and the associated density.
What carries the argument
The Smoluchowski coagulation equation with a derived gravitational kernel, solved by Monte Carlo sampling under the full-conditioning scheme.
If this is right
- Clusters containing more PBHs or higher density reach runaway merger at earlier redshifts.
- The PBH mass population shifts systematically toward higher masses as mergers proceed.
- The same initial conditions produce distinct final mass spectra at different redshifts.
- Runaway growth can occur within the cosmic time available before z approximately 10 for a range of cluster parameters.
Where Pith is reading between the lines
- The predicted mass spectra could be compared directly with future JWST or Roman telescope counts of high-redshift black holes.
- Early mergers inside clusters would generate a stochastic gravitational-wave background whose amplitude depends on the same cluster parameters.
- Adding external tidal fields or gas drag would shift the runaway timescales in a testable way.
Load-bearing premise
The derived coagulation kernel accurately models the physical merger rates inside PBH clusters without significant external perturbations.
What would settle it
A measured high-redshift black-hole mass function that lies outside the range of final distributions produced by the simulated runs for any plausible cluster number and density.
read the original abstract
In arXiv:2507.07171, we demonstrate that the high-redshift supermassive black holes in the so-called "little red dots" discovered by James Webb Space Telescope (JWST) can be explained by the primordial black hole (PBH) clustering on small scales. In this paper, we present a comprehensive simulation of the successive PBH mergers within a cluster by solving the Smoluchowski coagulation equation. We derive the coagulation kernel considering both cases with and without the effects of mass segregation. Then we employ the Monte Carlo method to solve the equation, implementing the full-conditioning scheme using the discrete inverse transformation method. Our simulations determine the runaway timescales of clusters and the mass population evolution of PBHs across a wide range of cosmic redshifts, depending on the number of PBHs within the cluster and the associated density.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to model successive mergers of primordial black holes (PBHs) in clusters by deriving coagulation kernels (with and without mass segregation) from gravitational dynamics and solving the Smoluchowski equation via a Monte Carlo method with full-conditioning and discrete inverse transformation. It reports runaway timescales and mass-population evolution across redshifts as functions of initial PBH number and cluster density, motivated by explaining JWST 'little red dots' via the clustering scenario of arXiv:2507.07171.
Significance. If the central results hold, the work supplies a quantitative Monte Carlo framework for PBH cluster evolution that could link small-scale clustering to early supermassive black hole growth, with potential falsifiable predictions for mass spectra at different redshifts.
major comments (2)
- [Kernel derivation (abstract and §2)] Coagulation kernel derivation: the kernel incorporates only internal gravitational dynamics and mass segregation. At the high redshifts and densities relevant to the little-red-dots scenario, Hubble drag, cosmic expansion, and external tidal fields are omitted; these can alter encounter velocities and effective merger rates by an O(1) factor, directly affecting the reported runaway timescales.
- [Monte Carlo solution and results] Monte Carlo implementation: the method is described as standard, yet the manuscript provides no validation against known analytic limits of the Smoluchowski equation, no error analysis or convergence tests, and no comparison to limiting cases (e.g., constant kernel or monodisperse initial conditions). This leaves quantitative support for the specific timescales and mass spectra unclear.
minor comments (1)
- [Initial conditions and redshift dependence] Clarify how the initial cluster density is mapped to physical scales at each redshift and whether the reported timescales include or exclude the Hubble time.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions we plan to make.
read point-by-point responses
-
Referee: [Kernel derivation (abstract and §2)] Coagulation kernel derivation: the kernel incorporates only internal gravitational dynamics and mass segregation. At the high redshifts and densities relevant to the little-red-dots scenario, Hubble drag, cosmic expansion, and external tidal fields are omitted; these can alter encounter velocities and effective merger rates by an O(1) factor, directly affecting the reported runaway timescales.
Authors: We acknowledge that the coagulation kernel is derived from internal gravitational dynamics within the cluster. For the high-density regimes considered, internal velocity dispersions are expected to dominate over the Hubble flow on cluster scales, but we agree that external effects merit discussion. In the revised manuscript we will add a dedicated paragraph estimating the magnitude of Hubble drag and tidal contributions, showing that they introduce at most an O(1) correction to the runaway timescales for the parameter ranges explored, and we will explicitly state this approximation as a limitation of the present model. revision: partial
-
Referee: [Monte Carlo solution and results] Monte Carlo implementation: the method is described as standard, yet the manuscript provides no validation against known analytic limits of the Smoluchowski equation, no error analysis or convergence tests, and no comparison to limiting cases (e.g., constant kernel or monodisperse initial conditions). This leaves quantitative support for the specific timescales and mass spectra unclear.
Authors: We agree that explicit validation strengthens the quantitative claims. In the revised version we will add a new subsection that compares the Monte Carlo results to known analytic solutions of the Smoluchowski equation for the constant kernel and for monodisperse initial conditions. We will also present convergence tests with respect to the number of Monte Carlo realizations and the number of particles, together with error estimates on the derived runaway timescales and mass spectra. revision: yes
Circularity Check
No significant circularity in the derivation chain
full rationale
The paper derives the coagulation kernel from gravitational dynamics (with and without mass segregation) and solves the Smoluchowski equation using Monte Carlo methods with the discrete inverse transformation scheme. These steps are self-contained once the initial number and density of PBHs are specified as inputs; the resulting runaway timescales and mass spectra follow directly from the numerical integration without reduction to fitted outputs or self-referential definitions. The citation to arXiv:2507.07171 supplies only contextual motivation for the existence of PBH clusters and does not bear the load of the kernel derivation or the simulation results. No ansatz is smuggled via self-citation, no uniqueness theorem is invoked from prior work, and no known empirical pattern is merely renamed. The derivation chain therefore stands independently of its own outputs.
Axiom & Free-Parameter Ledger
free parameters (1)
- initial PBH number and cluster density
axioms (1)
- domain assumption The coagulation kernel accurately captures gravitational merger rates in PBH clusters
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We derive the coagulation kernel considering both cases with and without the effects of mass segregation... Kij = ⟨σmerg(mi,mj)vrel⟩Fmsij (Eq. 2.12, 2.20)
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The Smoluchowski coagulation equation... ∂n(m,t)/∂t = ½∫K(m′,m−m′)n(m′)n(m−m′)dm′ − ... (Eq. 1.1)
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.
Reference graph
Works this paper leans on
-
[1]
B. Zhang, W.-X. Feng and H. An,Little Red Dots from Small-Scale Primordial Black Hole Clustering,2507.07171
-
[2]
Supermassive black holes from primordial black hole seeds
N. Duechting,Supermassive black holes from primordial black hole seeds,Phys. Rev. D70(2004) 064015, [astro-ph/0406260]
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[3]
Primordial seeds of supermassive black holes
M. Kawasaki, A. Kusenko and T. T. Yanagida,Primordial seeds of supermassive black holes,Phys. Lett. B711(2012) 1–5, [1202.3848]. – 29 –
work page Pith review arXiv 2012
- [4]
- [5]
-
[6]
B. J. Carr and S. W. Hawking,Black holes in the early Universe,Mon. Not. Roy. Astron. Soc. 168(1974) 399–415
work page 1974
-
[7]
B. J. Carr,The Primordial black hole mass spectrum,Astrophys. J.201(1975) 1–19
work page 1975
- [8]
-
[9]
Primordial Black Holes as Generators of Cosmic Structures
B. Carr and J. Silk,Primordial Black Holes as Generators of Cosmic Structures,Mon. Not. Roy. Astron. Soc.478(2018) 3756–3775, [1801.00672]
work page Pith review arXiv 2018
-
[10]
D. Inman and Y. Ali-Ha¨ ımoud,Early structure formation in primordial black hole cosmologies, Phys. Rev. D100(2019) 083528, [1907.08129]
- [11]
- [12]
-
[13]
H.-L. Huang, J.-Q. Jiang and Y.-S. Piao,High-redshift JWST massive galaxies and the initial clustering of supermassive primordial black holes,Phys. Rev. D110(2024) 103540, [2407.15781]
-
[14]
H.-L. Huang, Y.-Y. Li, Y.-T. Wang and Y.-S. Piao,Supermassive primordial black holes as seeds for the galaxies GHZ9 and UHZ1 recently observed by the JWST,Phys. Rev. D113(2026) 023519, [2410.05891]
-
[15]
P. Dayal and R. Maiolino,The properties of primordially-seeded black holes and their hosts in the first billion years: implications for JWST,Astron. Astrophys.706(2026) A72, [2506.08116]
-
[16]
Y. B. Zel’dovich and I. D. Novikov,The Hypothesis of Cores Retarded during Expansion and the Hot Cosmological Model,Sov. Astron.10(1967) 602
work page 1967
-
[17]
Hawking,Gravitationally collapsed objects of very low mass,Mon
S. Hawking,Gravitationally collapsed objects of very low mass,Mon. Not. Roy. Astron. Soc.152 (1971) 75
work page 1971
-
[18]
B. Carr, F. Kuhnel and M. Sandstad,Primordial Black Holes as Dark Matter,Phys. Rev. D94 (2016) 083504, [1607.06077]
work page Pith review arXiv 2016
-
[19]
B. Carr, K. Kohri, Y. Sendouda and J. Yokoyama,Constraints on primordial black holes,Rept. Prog. Phys.84(2021) 116902, [2002.12778]
work page internal anchor Pith review arXiv 2021
-
[20]
M. Gorton and A. M. Green,How open is the asteroid-mass primordial black hole window?, SciPost Phys.17(2024) 032, [2403.03839]
- [21]
-
[22]
B. Carr, A. J. Iovino, G. Perna, V. Vaskonen and H. Veerm¨ ae,Primordial black holes: constraints, potential evidence and prospects,2601.06024. [23]LIGO Scientific, Virgocollaboration, B. P. Abbott et al.,Observation of Gravitational Waves from a Binary Black Hole Merger,Phys. Rev. Lett.116(2016) 061102, [1602.03837]. [24]KAGRA, VIRGO, LIGO Scientificcoll...
-
[23]
S. Bird, I. Cholis, J. B. Mu˜ noz, Y. Ali-Ha¨ ımoud, M. Kamionkowski, E. D. Kovetz et al.,Did LIGO detect dark matter?,Phys. Rev. Lett.116(2016) 201301, [1603.00464]
work page Pith review arXiv 2016
-
[24]
Primordial Black Hole Scenario for the Gravitational-Wave Event GW150914
M. Sasaki, T. Suyama, T. Tanaka and S. Yokoyama,Primordial Black Hole Scenario for the Gravitational-Wave Event GW150914,Phys. Rev. Lett.117(2016) 061101, [1603.08338]
work page Pith review arXiv 2016
- [25]
- [26]
- [27]
-
[28]
Takhistov,Transmuted Gravity Wave Signals from Primordial Black Holes,Phys
V. Takhistov,Transmuted Gravity Wave Signals from Primordial Black Holes,Phys. Lett. B782 (2018) 77–82, [1707.05849]
-
[29]
Y. Ali-Ha¨ ımoud, E. D. Kovetz and M. Kamionkowski,Merger rate of primordial black-hole binaries,Phys. Rev. D96(2017) 123523, [1709.06576]
-
[30]
N. Bellomo, J. L. Bernal, A. Raccanelli and L. Verde,Primordial Black Holes as Dark Matter: Converting Constraints from Monochromatic to Extended Mass Distributions,JCAP01(2018) 004, [1709.07467]
- [31]
- [32]
- [33]
- [34]
-
[35]
M. Andr´ es-Carcasona, A. J. Iovino, V. Vaskonen, H. Veerm¨ ae, M. Mart´ ınez, O. Pujol` as et al., Constraints on primordial black holes from LIGO-Virgo-KAGRA O3 events,Phys. Rev. D110 (2024) 023040, [2405.05732]
- [36]
- [37]
-
[38]
V. De Luca, L. Del Grosso, G. Franciolini, K. Kritos, E. Berti, D. D’Orazio et al.,A cosmologist’s take on Little Red Dots,2512.19666
-
[39]
M. Aljaf and I. Cholis,The Merger Rate of Primordial Black Holes,2512.12227
- [40]
-
[41]
V. De Luca, V. Desjacques, G. Franciolini and A. Riotto,The clustering evolution of primordial black holes,JCAP11(2020) 028, [2009.04731]
-
[42]
G. Franciolini, K. Kritos, E. Berti and J. Silk,Primordial black hole mergers from three-body interactions,Phys. Rev. D106(2022) 083529, [2205.15340]
-
[43]
Stasenko,Redshift evolution of primordial black hole merger rate,Phys
V. Stasenko,Redshift evolution of primordial black hole merger rate,Phys. Rev. D109(2024) 123546, [2403.11325]. – 31 –
-
[44]
Stasenko,Merger history of clustered primordial black holes,Phys
V. Stasenko,Merger history of clustered primordial black holes,Phys. Dark Univ.50(2025) 102138, [2506.11290]
- [45]
- [46]
- [47]
-
[48]
S. Young and C. T. Byrnes,Primordial black holes in non-Gaussian regimes,JCAP08(2013) 052, [1307.4995]
-
[49]
S. Young and C. T. Byrnes,Long-short wavelength mode coupling tightens primordial black hole constraints,Phys. Rev. D91(2015) 083521, [1411.4620]
-
[50]
S. Young and C. T. Byrnes,Signatures of non-gaussianity in the isocurvature modes of primordial black hole dark matter,JCAP04(2015) 034, [1503.01505]
-
[51]
Y. Tada and S. Yokoyama,Primordial black holes as biased tracers,Phys. Rev. D91(2015) 123534, [1502.01124]
-
[52]
G. Franciolini, A. Kehagias, S. Matarrese and A. Riotto,Primordial Black Holes from Inflation and non-Gaussianity,JCAP03(2018) 016, [1801.09415]
-
[53]
V. Desjacques and A. Riotto,Spatial clustering of primordial black holes,Phys. Rev. D98(2018) 123533, [1806.10414]
-
[54]
Y. Ali-Ha¨ ımoud,Correlation Function of High-Threshold Regions and Application to the Initial Small-Scale Clustering of Primordial Black Holes,Phys. Rev. Lett.121(2018) 081304, [1805.05912]
-
[55]
S. Young and C. T. Byrnes,Initial clustering and the primordial black hole merger rate,JCAP03 (2020) 004, [1910.06077]
-
[56]
T. Suyama and S. Yokoyama,Clustering of primordial black holes with non-Gaussian initial fluctuations,PTEP2019(2019) 103E02, [1906.04958]
- [57]
-
[58]
V. De Luca, G. Franciolini and A. Riotto,Heavy Primordial Black Holes from Strongly Clustered Light Black Holes,Phys. Rev. Lett.130(2023) 171401, [2210.14171]
-
[59]
V. De Luca, G. Franciolini, A. Riotto and H. Veerm¨ ae,Ruling Out Initially Clustered Primordial Black Holes as Dark Matter,Phys. Rev. Lett.129(2022) 191302, [2208.01683]
-
[60]
M. Y. Khlopov, S. G. Rubin and A. S. Sakharov,Primordial structure of massive black hole clusters,Astropart. Phys.23(2005) 265, [astro-ph/0401532]
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[61]
Quasars formation around clusters of primordial black holes
V. Dokuchaev, Y. Eroshenko and S. Rubin,Quasars formation around clusters of primordial black holes,Grav. Cosmol.11(2005) 99–104, [astro-ph/0412418]
work page internal anchor Pith review Pith/arXiv arXiv 2005
- [62]
-
[63]
H.-L. Huang and Y.-S. Piao,Toward supermassive primordial black holes from inflationary bubbles,Phys. Rev. D110(2024) 023501, [2312.11982]
- [64]
- [65]
- [66]
-
[67]
C. Animali and V. Vennin,Clustering of primordial black holes from quantum diffusion during inflation,JCAP08(2024) 026, [2402.08642]
-
[68]
V. De Luca, G. Franciolini and A. Riotto,Bubble correlation in first-order phase transitions,Phys. Rev. D104(2021) 123539, [2110.04229]
-
[69]
J. Matthee et al.,Little Red Dots: An Abundant Population of Faint Active Galactic Nuclei at z∼ 5 Revealed by the EIGER and FRESCO JWST Surveys,Astrophys. J.963(2024) 129, [2306.05448]
-
[70]
F. Pacucci, B. Nguyen, S. Carniani, R. Maiolino and X. Fan,JWST CEERS and JADES Active Galaxies at z = 4–7 Violate the Local M •–M⋆ Relation at>3σ: Implications for Low-mass Black Holes and Seeding Models,Astrophys. J. Lett.957(2023) L3, [2308.12331]
-
[71]
F. Gentile et al.,Not-so-little Red Dots: Two Massive and Dusty Starbursts at z∼5–7 Pushing the Limits of Star Formation Discovered by JWST in the COSMOS-Web Survey,Astrophys. J. Lett.973(2024) L2, [2408.10305]
- [72]
- [73]
-
[74]
A. Bogdan et al.,Evidence for heavy-seed origin of early supermassive black holes from a z≈10 X-ray quasar,Nature Astron.8(2024) 126–133, [2305.15458]
- [75]
- [76]
-
[77]
Limits on primordial black holes from $\mu$ distortions in cosmic microwave background
T. Nakama, B. Carr and J. Silk,Limits on primordial black holes fromµdistortions in cosmic microwave background,Phys. Rev. D97(2018) 043525, [1710.06945]
work page Pith review arXiv 2018
-
[78]
W.-H. Ruan, Z.-K. Guo, R.-G. Cai and Y.-Z. Zhang,Taiji program: Gravitational-wave sources, Int. J. Mod. Phys. A35(2020) 2050075, [1807.09495]. [82]TianQincollaboration, J. Mei et al.,The TianQin project: current progress on science and technology,PTEP2021(2021) 05A107, [2008.10332]
-
[79]
Detecting a gravitational-wave background with next-generation space interferometers
H. Kudoh, A. Taruya, T. Hiramatsu and Y. Himemoto,Detecting a gravitational-wave background with next-generation space interferometers,Phys. Rev. D73(2006) 064006, [gr-qc/0511145]
work page internal anchor Pith review Pith/arXiv arXiv 2006
-
[80]
Kawamura et al., PTEP2021, 05A105 (2021), arXiv:2006.13545 [gr-qc]
S. Kawamura et al.,Current status of space gravitational wave antenna DECIGO and B-DECIGO, PTEP2021(2021) 05A105, [2006.13545]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.