
Heavy primordial black holes from strongly clustered light black holes
Valerio De Luca,
1, ∗
Gabriele Franciolini,
2, 3, †
and Antonio Riotto
4, ‡
1
Center for Particle Cosmology, Department of Physics and Astronomy,
University of Pennsylvania 209 S. 33rd St., Philadelphia, PA 19104, USA
2
Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, 00185, Roma, Italy
3
INFN, Sezione di Roma, Piazzale Aldo Moro 2, 00185, Roma, Italy
4
Département de Physique Théorique and Gravitational Wave Science Center (GWSC),
Université de Genève, CH-1211 Geneva, Switzerland CH-1211 Geneva, Switzerland
We show that heavy primordial black holes may originate from much lighter ones if the latter are
strongly clustered at the time of their formation. While this population is subject to the usual
constraints from late-time universe observations, its relation to the initial conditions is different from
the standard scenario and provides a new mechanism to generate massive primordial black holes even
in the absence of efficient accretion, opening new scenarios, e.g. for the generation of supermassive
black holes.
Introduction. Multiple detections of gravitational waves
(GWs) coming from black hole binary mergers [
1
–
4
]
have revived the interest in the physics of Primor-
dial Black Holes (PBHs) [
5
–
8
]. Indeed, some of the
LIGO/Virgo/KAGRA data may be of primordial origin
[
9
–
16
] and future GW experiments will shed light on the
possible existence of PBHs [17–22].
PBHs in the early universe are commonly born in the
radiation-dominated phase (see Ref. [
7
] for a review on
the various formation mechanisms). Given our ignorance
of the production mechanism giving rise to PBHs, if any,
we do not know if they are born randomly distributed or
with a strong correlation among them. In the standard
scenario where PBHs are generated by the collapse of large
overdensities created during inflation on small scales [
5
],
PBHs are Poisson distributed in space [
23
–
26
]. However,
large initial correlations are conceivable and not ruled out,
unless the PBH abundance in the universe (normalised to
the dark matter)
fPBH
is larger than
O
(0
.
1) in the stellar
mass range [27].
In this Letter we propose a new mechanism to generate
heavy PBHs from light ones making use of the large initial
PBH clustering, a mechanism that we dub “clusteringen-
esis”. The idea is quite simple: if PBHs are born close to
each other, the strong gravitational interactions among
them may result in the collapse of this clump into a more
massive PBH, even in the radiation dominated phase
of the early universe. This mechanism provides a novel
way to increase the mass of light PBHs in the primordial
epochs, which is alternative to the more standard process
relying on baryonic mass accretion, which is efficient only
for PBHs with mass larger than
O
(10)
M
[
28
] at lower
redshifts.
In the following we will describe the basics of this idea
and discuss some of its possible implications.
Collapse of large overdensities in the radiation-dominated
era. Independently from the formation mechanism and be-
ing discrete objects, the most generic initial two-point cor-
relator for the PBH density contrast
δPBH
=
δρPBH/ρPBH
acquires the form [23]
hδPBH(~r)δPBH(0)i=1
nPBH
δD(r) + ξPBH(r),(1)
in terms of their distance r, where
nPBH '30fPBH MPBH
M−1
kpc−3(2)
is the average PBH number density per comoving volume
for a monochromatic PBH population with mass MPBH .
We suppose that, at the time of formation, such two-
point correlator is dominated by the reduced correlation
function
ξPBH
(
r
)up to some comoving clustering scale
rcl
,
while on larger scales the Poisson shot noise, arising from
the discrete nature of PBHs, dominates. For simplicity
and for the sake of the argument, we assume an approx-
imately constant in space and large reduced two-point
correlation function up to rcl,
ξPBH(r)'(ξ01 for r∼
<rcl,
0 otherwise.(3)
We will be agnostic in the following regarding the origin
of such correlations. Let us just point out that large
clustering appears if the local properties of the PBH
overdensity field are space dependent. This effect might
either come from an actual field different from the over-
density field, or from a long wavelength modulation of the
overdensity field itself, resulting from a self-coupling of
long and short scales as happens, e.g., in local models of
non-Gaussianity [
29
]. Alternatively, large correlation may
arise if PBHs are formed by bubble collisions in first-order
phase transitions [
30
–
32
] or for PBHs generated thanks
to long-range scalar forces [
33
]. These models can be used
to construct explicit realisations of the scenario proposed
in this work. Once rescaled to the total dark matter den-
sity
δρPBH/ρDM ≡
Φ
'ξ0fPBH
, the corresponding density
contrast is suppressed by the PBH abundance, which we
will assume to be tiny in the following in order to avoid
arXiv:2210.14171v2 [astro-ph.CO] 15 Apr 2023