2
just the recorded black hole mergers during the third
observing run (O3) in the framework of the PBH sce-
nario, while more realistic halo models (e.g., those with
ellipsoidal-collapse) can generate consistent PBH merg-
ers with gravitational wave observations, see our previous
studies in Refs. [32–35] for more details.
On the other side, it is believed that the large-scale vis-
ible Universe appears to have a web-like structure called
the cosmic web. In some parts of the cosmic web, clus-
ters are separated by large, almost empty regions called
cosmic voids where the density of such regions is lower
than the average density of the Universe. Being a ma-
jor part of the cosmic web, cosmic voids are exception-
ally underdense regions containing matter, which evac-
uates them toward other regions. Researches in recent
years show that there are galaxies and dark matter halos
in cosmic voids [36–39]. According to some theoretical
models, dark matter particles are expected to emit de-
tectable gamma-ray signals as a result of their decay and
annihilation [40–44]. Additionally, a diffuse background
of gamma rays can be detected across the sky by the
present gamma-ray observatories [45]. This background
consists of unknown signals that remain after subtracting
the contribution from all possible astrophysical sources,
such as supermassive black holes and pulsars. Besides,
such signals have a non-uniform distribution at different
spatial angles, which is adjustable to what is expected
from dark matter emission [46]. Accordingly, the way
of emitting signals related to the dark matter from over-
dense and underdense structures of the Universe has been
simulated [47]. The results indicate that although overall
dark matter signals from cosmic voids are weaker, they
are less contaminated by astrophysical sources, making
them easier to detect. In light of everything discussed
so far, the existence of PBHs in cosmic voids is not far
from expected. Therefore, it seems interesting to calcu-
late their merger rate in underdense environments with
minimal contamination by astrophysical sources.
In this work, we propose to calculate the merger rate
of PBHs in the medium of cosmic voids. In this respect,
the outline of the work is as follows. In Sec. II, we briefly
discuss cosmic voids and the need to include them in
cosmological studies. Then, in Sec. III, we present a con-
venient dark matter halo model in cosmic voids and de-
scribe some related quantities like halo density profile,
halo concentration parameter, and halo mass function.
Also, in Sec. IV, we calculate the merger rate of PBHs
in cosmic voids, as well as those in other structures. Fi-
nally, we discuss the results and summarize the findings
in Sec. V.
II. COSMIC VOIDS
The visible Universe on large scales seems to have a
web-like structure called the cosmic web. Such a struc-
ture is the result of the time evolution of primordial den-
sity fluctuations. There are smaller structures inside the
cosmic web, including knots, filaments, sheets, and voids,
within which matter is distributed differently. Since the
initial density field is a Gaussian random field defined
by the power spectrum of density fluctuations, and since
density fluctuations evolve due to gravity, the gravita-
tional field increases the density contrast in the Universe.
As a result, parts of the Universe with stronger gravita-
tional fields become denser over time, while parts with
less density become even more empty [48].
In the light these arguments, most of the matter can
be distributed in knots, filaments, and sheets. However,
clusters are separated by large, almost empty regions
called cosmic voids. The density of such regions is lower
than the average density of the Univers. Actually, cos-
mic voids are underdense regions with a low density of
matter. According to recent cosmological studies, cosmic
voids can provide clues about cosmic mass distribution
and serve as a convenient medium for constraining cos-
mological parameters [49–51]. In addition, by taking ad-
vantage of the dynamics governing cosmic voids, many
studies have been conducted on baryon acoustic oscilla-
tions [52,53], dark matter-dark energy interaction [54],
cosmic microwave background [55], and many other hot
topics.
In Ref. [56], a comparison between the properties of
various galaxies within cosmic voids has been performed,
in which the number density of galaxies is less than 10%
of the mean density of the Universe. Also, to find out the
properties of cosmic voids, several studies have been car-
ried out via N-body simulations [57–59]. By considering
the size, shape, and structure of cosmic voids, it can be
found that such regions might play a prominent role in
cosmological evolution because they cover a large volume
of the Universe.
Based on the cosmological perturbation theory in the
formation and evolution of large-scale structures, a vast
majority of matter budget is distributed in dark mat-
ter halos [60]. Numerical simulations and explanatory
studies indicate that the gravitational enhancement of
density fluctuations leads to the formation of large-scale
structures [61]. Moreover, a few studies have mapped the
population of cosmic voids within the local Universe [62–
64]. Large regions of cosmic voids are related to the fore-
most noticeable viewpoint of the megaparsec-scale Uni-
verse. Cosmic voids are attributed to enormous regions
with sizes about (20 ∼50) Mpc h−1that have a lower
distribution of matter compared to other structures and
possess a significant share of the Universe [65].
Although some studies have been performed on black
holes in cosmic voids, e.g., [66–68], no attention has been
paid to the merger rate PBHs and their evolution in such
regions. As mentioned earlier, it has been suggested that
there are galaxies and dark matter halos in cosmic voids.
Hence, it is likely that PBHs can be clustered in cos-
mic voids. Furthermore, research continues on develop-
ing new methods for finding galaxies, molecular gas, and
star formation in cosmic voids [69,70]. Cosmic voids
are interesting from a theoretical perspective because of