ORIGINAL RESEARCH article

Front. Astron. Space Sci., 28 May 2021

Sec. Cosmology

Volume 8 - 2021 | https://doi.org/10.3389/fspas.2021.681084

A Brief Review on Primordial Black Holes as Dark Matter

  • Instituto de Física Corpuscular (IFIC), CSIC-Universitat de València, Paterna, Spain

Abstract

Primordial black holes (PBHs) represent a natural candidate for one of the components of the dark matter (DM) in the Universe. In this review, we shall discuss the basics of their formation, abundance and signatures. Some of their characteristic signals are examined, such as the emission of particles due to Hawking evaporation and the accretion of the surrounding matter, effects which could leave an impact in the evolution of the Universe and the formation of structures. The most relevant probes capable of constraining their masses and population are discussed.

1 Introduction

The hypothesis of the formation of black holes (BHs) in the early Universe was first suggested in 1967 by Zeldovich and Novikov (Zel’dovich and Novikov, 1967), and independently by Hawking in 1971 (Hawking, 1971). Soon after, the possibility that primordial black holes (PBHs) could account for at least part of the DM became obvious (; Meszaros, 1975). At that time, the DM question started to be outlined as one of the fundamental problems in cosmology [see, e.g., the reviews in Sanders (2010); ]. DM (partly) composed by PBHs constitutes an exciting possibility, presenting an enormous number of observable signatures which can constrain its parameter space, as shall be detailed in Sec. 6. The variety of phenomenological effects produced by PBHs allows placing stringent bounds on the abundance of PBHs, usually indicated by the energy fraction of DM as PBHs, , with and the ratios of energy densities and critical density. Moreover, since PBHs are usually expected to be formed before nucleosynthesis, BBN constraints on the baryon abundance do not apply to them, as they do not intervene in the nucleosynthesis of elements, and thus, can be regarded as non-baryonic DM (). Several recent reviews are devoted to discuss PBHs in great detail (see, e.g., Sasaki et al., 2018; ; ; Green and Kavanagh, 2021).

Shortly after the first detection of gravitational waves from a merger of BHs by LIGO (), the question whether these could be of primordial nature was raised (). Analysis of posterior data from the gravitational wave detectors LIGO and Virgo showed that the detected mergers are compatible with the hypothesis of their components being of primordial nature, although there is no strong preference over stellar BHs (Sasaki et al., 2016; ; Clesse and García-Bellido, 2017; ; Garcia-Bellido et al., 2020; De Luca et al., 2020b; De Luca et al., 2021; Wong et al., 2021). PBHs with a lognormal mass function have been claimed to better fit data than BHs from astrophysical origin (Dolgov et al., 2020), although this is in contrast to the results ofHall et al. (2020), and a mixed population seems compatible or even favored (Hütsi et al., 2021).

Unlike stellar BHs, formed from the collapse of a massive star, which can present masses only above [the Tolman–Oppenheimer–Volkoff limit (Oppenheimer and Volkoff, 1939; Tolman, 1939)], PBHs could be produced with any mass. Thus, a positive measurement of a BH with a mass lower than would be a confirmation of the existence of primordial, non-stellar BHs (). PBHs could also conform intermediate mass BHs, with masses between and , too massive to be originated from a single star. It is the case of the merger event of BHs with masses and , producing a remnant BH of , in a so far mostly unobserved range of masses (). Furthermore, PBHs could constitute the seeds for super massive black holes (SMBHs), present in the nuclei of most galaxies, with masses ranging from to , and already existing at redshifts (). Such massive objects can hardly be produced from accreting stellar remnant BHs (Volonteri, 2010). However, the existence of massive enough PBHs may act as seeds for the SMBHs, from which they could have grown by accretion.

PBH formation is already present in standard cosmologies, although extremely unlikely. However, their production usually requires some exotic inflationary scenarios or physics beyond the Standard Model (BSM) in order to obtain a large enough abundance. The typically considered formation mechanism of PBHs arises from the direct collapse of primordial fluctuations, whose power is enhanced at small scales as a consequence of some inflationary potential, as we shall further comment on below. There are, however, other scenarios which naturally predict a population of PBHs as a result of phase transitions in the early Universe and by the collapse of topological defects (Hawking et al., 1982; Hawking, 1989; Polnarev and Zembowicz, 1991; Garriga et al., 2016; Deng et al., 2017; Deng and Vilenkin, 2017). Hence, the existence of PBHs would provide valuable hints about the still unknown physics of the very early Universe (see, e.g., Polnarev and Khlopov, 1985; Khlopov, 2010; ), and may allow to probe high-energy scales and supersymmetric theories (Ketov and Khlopov, 2019).

In this review, the most relevant aspects of PBHs as DM are briefly discussed, such as the mechanism of formation in Section 2, the initial abundance and mass distribution in Section 3, the process of accretion in Section 4 and other PBH features, which could leave imprints on different observables, in Section 5. Current observational constraints on their population are summarized in Section 6, and we conclude in Section 7.

2 Formation and Conditions of Collapse

The mass of a BH which collapsed in the early Universe depends on its formation time. A BH can be characterized by an extremely dense amount of matter in a very compact region, i.e., within the well-known Schwarzschild radius, . Thus, the mean density inside that region can be estimated as . On the other hand, the mean density of the universe in the radiation era scales with time as . In order to have PBH formation, densities at least of the order of the mean inside the BH horizon, , are needed. Therefore, the mass of the resulting PBHs should be of the order of the horizon mass at that time, i.e., the mass within a region of the size of the Hubble horizon, (), which is defined as

PBHs with masses of would have been formed at around the QCD phase transition, at . Since the PBH mass is roughly given by the mass within the horizon, it means that fluctuations entering the horizon can collapse into PBHs. A detailed calculation shows that , where the proportionality factor γ depends on the details of gravitational collapse, and gets values lower than 1. Early estimates showed that it can be approximated as , with the sound speed at the radiation epoch () More refined results show that the PBH mass is given by a scaling relation with the overdensityδ, (Niemeyer and Jedamzik, 1998; Niemeyer and Jedamzik, 1999), where κ and α are order unity constants, which depend on the background cosmology and on the shape of the perturbation (Niemeyer and Jedamzik, 1998; Niemeyer and Jedamzik, 1999; Musco et al., 2005), and is the collapse threshold [see, e.g., Escrivà et al., 2021; Musco et al., 2021 for recent accurate computations]. Figure 1 depicts a sketch of the process of PBH formation.

FIGURE 1

Since PBHs are formed when fluctuations cross the horizon by the time of formation, , their mass can be related to the wavelength of perturbations. When the mode of wavenumber k crosses the horizon, the condition holds. Since the mass of PBHs is proportional to the horizon mass at the moment of formation, , at the radiation dominated era (Sasaki et al., 2018),

Hence, probing a given scale k could constrain a PBH population of its corresponding mass. Furthermore, an enhancement in the power spectrum around that scale would result in a large number of PBHs of such masses. In this review, we shall only consider PBHs formed during the radiation era. Those produced before inflation ends would have been diluted due to their negligible density during the inflationary accelerated expansion. PBHs formed during the matter-dominated era, or in an early matter domination era previous to the radiation era, have also been considered in the literature, and may have different imprints, since the conditions of collapse are less restrictive, and could start from smaller inhomogeneities (see, e.g., Green and Kavanagh, 2021).

3 Abundance and Mass Function OF PBHS

It is possible to estimate the initial abundance of PBHs at the moment of formation, taking into account all overdensities above the threshold for collapse, . Assuming a Gaussian probability distribution, , for the overdensities with variance at a mass scale M, the initial abundance, defined as , is given by (Sasaki et al., 2018; Green and Kavanagh, 2021)where in the last equality, has been assumed. For the standard cosmological scenario with an initial scale invariant power spectrum, at CMB scales, the amplitude of the fluctuations is around , leading to , which is completely negligible (Green, 2014). Therefore, in order to have a relevant population of PBHs, larger values of the initial power spectrum are needed. On the other hand, the assumption of a Gaussian distribution may not be consistent with enhanced fluctuations and the presence of PBHs, so deviations are unavoidable (Franciolini et al., 2018; De Luca et al., 2019b) [except for specific inflation models presenting an inflection point ]. Non-gaussianities could have a great impact on the initial fraction and lead to a larger population, as well as leaving detectable signatures in gravitational waves (). Finally, note that the above formula follows the simple Press-Schechter approach, whereas to account for the non-universal nature of the threshold, the use of peak theory provides more accurate results (Green et al., 2004; Germani and Musco, 2019). Its validity, however, is limited to relatively narrow power spectra, while broader spectra require the use of non-linear statistics (Germani and Sheth, 2020).

Nonetheless, although the initial fraction β is a very small quantity, since matter and radiation densities scale differently with redshift [ and respectively], the PBH contribution can become relevant at current times. The fraction can be related to the current density parameters of PBH and radiation, and , as (),

For initial fractions as low as of solar mass BHs, the fraction of energy in PBHs could be, thus, of order unity.

Depending on the specific mechanism of formation, a population of PBHs with different masses could be generated. The specific shape of the enhancement of fluctuations determines the mass distribution function. Sharp peaks in the power spectrum imply approximately monochromatic distributions. For instance, chaotic new inflation may give rise to relatively narrow peaks (Yokoyama, 1998; Saito et al., 2008). However, inflation models with an inflection point in a plateau of the potential (García-Bellido, 2017; García-Bellido and Ruiz Morales, 2017), or hybrid inflation (), predict, instead, extended mass functions, which can span over a large range of PBHs masses. In this review, we focus on monochromatic distributions for simplicity, although it is possible to translate constraints to extended mass functions, which can be very stringent despite the fact of having more parameters to fit (; ). Nonetheless, even if there are bounds that exclude in the monochromatic case, there are choices for the mass function which allow PBHs to constitute most or all of the DM abundance (Lehmann et al., 2018; ; Sureda et al., 2020).

4 Accretion Onto PBHS

One of the consequences of the existence of PBHs with greater impact on different observables is the process of accretion. Infalling matter onto PBHs would release radiation, injecting energy into the surrounding medium, and strongly impacting its thermal state, leaving significant observable signatures. The physics of accretion is highly complex, but one can attempt a simplified approach considering the spherical non-relativistic limit, following the seminal work by . In this framework, the BH is treated as a point mass surrounded by matter, embedded in a medium which tends to constant density far enough from the BH. The relevant scale is the so-called Bondi radius, defined bywhere is the speed of sound far enough from the PBH. At distances around , the accretion process starts to become relevant, and the velocity of the accreted matter reaches , when the density is still close to the boundary value, . Taking into account the velocity of the BH relative to the medium, , one can write the accretion rate as ()where the parameter λ is the dimensionless accretion rate and depends on the equation of state, but it is order for the cases of interest (). In the early Universe, can be estimated as the baryon-DM relative velocity, computed in linear theory. Its root-mean-square value is approximately constant before recombination, dropping linearly with at later times ().

On the other hand, real BHs spin and form an accreting disk. Thus, the spherical symmetric case may not be applicable. Even though PBH spins are expected to be small, an accretion disk would form if the angular momentum is large enough to keep matter orbiting at Keplerian orbits at distances much larger than the innermost stable orbits, which are roughy given by the Schwarzschild radius (). Applying this criterion, the formation of accretion disks around PBHs has been suggested to occur if the condition is fulfilled (Poulin et al., 2017b). This is satisified for and large enough abundances at the epoch of CMB decoupling or at later times. Some of the results outlined above are still valid if the dimensionless accretion rate λ is suppresed by roughly two orders of magnitude, after accounting for viscosity effects and matter outflows through jets (Poulin et al., 2017b).

The matter falling onto BHs is greatly accelerated, which gives rise to radiative emission of high-energy photons by processes such as bremsstrahlung. The luminosity of accreting BHs is proportional to its accretion rate, and can be written as (Xie and Yuan, 2012)where denotes the radiative efficiency, which is in general a complicated function of the accretion rate , and depends upon the geometry, viscosity and other hydrodynamical considerations.

If the accretion disk is optically thin, most of the energy released through viscous dissipation is radiated away, and the luminosities obtained can be close to the Eddington luminosity, , explaining the extreme brightness of many far Active Galactic Nuclei (AGN) (Shakura and Sunyaev, 1973). However, PBHs, like nearby astrophysical BHs, may radiate in a much less efficient way, through the Advective-Dominated-Accretion-Flow (ADAF) (Ichimaru, 1977; Rees et al., 1982; Narayan and Yi, 1994; Narayan and Yi, 1995) (see, e.g., Yuan and Narayan, 2014 for a review). In this scenario, the dynamics is ruled by advective currents, forming a hot thick disk. Most of the emitted energy is deposited in the same accretion disk, heating it up. Thus, only a small portion of energy is released to the surrounding medium, the radiative process being inefficient. In the ADAF scenario, the efficiency function can be fitted by a broken power-law formula, with the slopes and amplitudes dependent on the mass range and the specific modeling of viscosity effects (Xie and Yuan, 2012).

Finally, the energy emitted in the accretion processes is deposited through different channels into the medium. The energy deposition rate for each channel reads (; Poulin et al., 2017b)where the subscript c denotes the channel in which energy is deposited, namely: ionization, heating, or atomic excitations (where the transitions are the most relevant ones). The energy deposition factors quantify the fraction of energy which goes to the different channels, and has been computed numerically (Slatyer, 2016). This energy injection would affect different types of observables, as we briefly outline below.

5 Other Generic Features

A noteworthy phenomenon of BHs is that of evaporation. Due to quantum effects in curved spacetimes, BHs may emit particles at their event horizon, as was noticed in Hawking (1974). The emitted radiation would have a nearly thermal black body spectrum, with a temperature given by ()

which is known as Hawking temperature. Due to this emission, BHs would slowly lose mass until completely evaporate. The lifetime of a PBH of initial mass M is (Lopresto, 2003)

Thus, the lower the PBH mass, the earlier it evaporates. Those with masses of 1015 g or below would have already evaporated by now, having lifetimes shorter than the age of the Universe (Page, 1976), so they cannot contribute to the current DM abundance. These evaporation products or the effects they produce in different observables can be search for in a variety of experiments, probing different mass ranges. Detailed computations of the emitted spectra can be performed by codes such as BlackHawk .

Another important feature is that of clustering. If fluctuations are originally Gaussian distributed and around a relatively narrow peak, PBHs are not expected to be originated in clusters, being initially randomly distributed on small scales (; Desjacques and Riotto, 2018). However, either primordial non-gaussianities or a broad peak in the power spectrum could lead to a significant initial clustering (; Suyama and Yokoyama, 2019) [although broad spectra have also been argued not to produce appreciable clustering (Moradinezhad Dizgah et al., 2019)]. Anyway, PBHs could become bounded as the Universe evolves. A proper determination of their clustering properties at later times is of great importance, for instance, in order to estimate their merger rates (Desjacques and Riotto, 2018; De Luca et al., 2020a). Indeed, the formation of clusters could alleviate some constraints on the PBH abundance ().

On another hand, since PBHs would be formed from the collapse of high density peaks relatively spherically symmetric, their torques and angular momentum are expected to be small (; Mirbabayi et al., 2020). It is usually quantified with the dimensionless spin parameter, , where S is the spin. Estimations of for PBHs show that it is a small quantity, equal or lower than 0.01 (De Luca et al., 2019a). In contrast, astrophysical BHs are expected to have substantially larger spins, since angular momentum must be conserved during the collapse of their stars of origin, which are often rotating. Hence, the spin can serve as a good proxy to distinguish the nature of a population of BHs. The measurement of low spin parameters could represent a hint for the detection of PBHs. The latest Bayesian analyses of LIGO/Virgo mergers suggest that low values of the spin parameter are strongly preferred by data, regardless of the assumed priors (Garcia-Bellido et al., 2020). Note, however, that the PBH mass and spin depend on the accretion mechanism and their time evolution is correlated (De Luca et al., 2020c).

Furthermore, due to the discrete nature of PBHs, a Poisson shot noise contribution to the matter power spectrum, constant in wavenumber, , would be expected (; Gong and Kitajima, 2017). PBHs fluctuations give rise to isocurvature modes (; ; Inman and Ali-Haïmoud, 2019), and thus, affect only at scales smaller than the horizon at the epoch of matter-radiation equality (Peacock, 1999). Therefore, this leads to an enhancement in the matter power spectrum, increasing the population of low-mass halos, which can be constrained by large scale structure and Lyα forest analyses. This effect is different from the one induced by other non-CDM candidates, such as warm DM or fuzzy DM, which suppress small scale fluctuations, washing out small structures. This contribution becomes relevant for low-mass halos not large enough to cool and collapse to form stars, which are commonly known as minihalos. It has been argued that this enhancement may produce a non-negligible 21 cm signal from the neutral hydrogen in minihalos (Gong and Kitajima, 2017, Gong and Kitajima, 2018). However, consistently accounting for the heating of the IGM due to PBH accretion increases the Jeans mass and suppresses the minihalo 21 cm signal, making it almost negligible (Mena et al., 2019).

6 Observational Constraints on PBHS as DM

PBHs can impact cosmology and astrophysics in a wide range of ways, leaving different observational effects which allow to constrain their properties. In this section, we review the most important bounds on the current fraction of PBHs as DM, , for a wide range of masses , for monochromatic mass functions. A collection of limits from the different probes is depicted in Figure 2. For a more comprehensive list of constraints, see Green and Kavanagh (2021); .

FIGURE 2

) and on the CMB spectrum (); non-observation of microlensing events (blue) from the MACHO (), EROS (Tisserand et al., 2007), Kepler (Griest et al., 2014), Icarus (Oguri et al., 2018), OGLE (Niikura et al., 2019b) and Subaru-HSC (Croon et al., 2020) collaborations; PBH accretion signatures on the CMB (orange), assuming spherical accretion of PBHs within halos (Serpico et al., 2020); dynamical constraints, such as disruption of stellar systems by the presence of PBHs (green), on wide binaries (Monroy-Rodríguez and Allen, 2014) and on ultra-faint dwarf galaxies (); power spectrum from the Lyα forest (cyan) (Murgia et al., 2019); merger rates from gravitational waves (purple), either from individual mergers (Kavanagh et al., 2018; ) or from searches of stochastic gravitational wave background (). Gravitational waves limits are denoted by dashed lines, since they could be invalidated (). Dotted brown line corresponds to forecasts from the 21 cm power spectrum with SKA sensitivities (Mena et al., 2019) and from 21 cm forest prospects (Villanueva-Domingo and Ichiki, 2021). Figure created with the publicly available Python code PBHbounds (Kavanagh, 2019).

Evaporation

Since BHs emit energy due to Hawking radiation, those with a lifetime shorter than the age of the Universe must have disintegrated nowadays, a fact which excludes PBHs with to form part of the current DM (Page, 1976). Moreover, PBHs with masses small enough, although still present, should emit a strong γ ray and cosmic ray background which could be observed. Absence of its detection strongly constrains the range of masses . In particular, the maximum fraction allowed is , with (). Comparable limits have been found from INTEGRAL and COMPTEL observations of the Galactic Center (DeRocco and Graham, 2019; Laha, 2019; Coogan et al., 2020; Laha et al., 2020). Furthermore, bounds from the isotropic X-ray and soft gamma-ray background have also been recently updated (; Iguaz et al., 2021). Additionally, data on the diffuse supernova neutrino background at Super-Kamiokande are also able to set constraints (Dasgupta et al., 2020).

Microlensing

If a compact object crosses the line of sight of a star, it may produce a so-called microlensing effect, which implies a transient and achromatic amplification of its flux. The range of masses of the objects which can produce it span from to (Paczynski, 1986; Green and Kavanagh, 2021). The non-detection of these events leads to bounds on the maximum abundance of PBHs about by the MACHO () and EROS (Tisserand et al., 2007) surveys in the Large and Small Magellanic Clouds, the Subaru Hyper Suprime Cam (HSC) in M31 (Andromeda) (Niikura et al., 2019a) and the Optical Gravitational Lensing Experiment (OGLE) in the Galactic bulge (Niikura et al., 2019b). Nonetheless, the existence of Earth-mass PBHs with a fraction could explain the observed excess of six microlensing events found in the OGLE data (Mróz et al., 2017), which is consistent with other constraints in this range of PBHs masses (Niikura et al., 2019b). Although this may constitute a hint of their existence, it cannot be regarded as a detection of PBHs, since these microlensing observations could also be explained by free-floating planets. There are some caveats regarding the results of the MACHO collaboration (Hawkins, 2015), since the limits reported are model dependent and could be biased by the assumption of an over-massive halo. Moreover, the results of the MACHO and EROS projects have been found to be statistically incompatible. Therefore, these bounds are not completely reliable, and PBHs could not be definitely ruled out within this range of masses. Variations of the lensing effect in type Ia Supernovae (Zumalacárregui and Seljak, 2018) and gamma-ray bursts () have also been proposed to constrain the PBH parameter space, although these limits have been later invalidated (Garcia-Bellido et al., 2017; Katz et al., 2018).

Gravitational Waves

The observation of BH mergers by LIGO and Virgo collaborations can be employed to constrain the allowed number of PBHs. Demanding that the predicted merger rates of PBH binaries cannot exceed the ones measured by gravitational waves, tight upper bounds of have been found for PBHs masses between 1 and 300 (; Kavanagh et al., 2018). The non-observation of a stochastic gravitational wave background of mergers expected from a population of PBHs has also been used to constrain their abundace (). However, to derive these limits, BHs have been treated as Schwarzschild BHs, while it would be more appropiate to use cosmological BH solutions embedded in a FLRW metric, such as the Thakurta metric (Thakurta, 1981). This implies a time-dependent mass, and that PBH binaries created before galaxy formation would have merged at much earlier times, allowing to obtain merger rates consistent with the LIGO observations, and completely avoiding these constraints ().

Dynamical Constraints

Due to two-body interactions, kinetic energies of systems of different masses usually become balanced and match. If a stellar system also has an additional MACHO population, its stars would gain kinetic energy and, due to the virial theorem, the system would expand. Therefore, the presence of PBHs would dynamically heat star clusters, making them larger and with higher velocity dispersions, leading to an eventual dissolution into its host galaxy. Populations with high mass to luminosity ratios are more sensitive to this effect, as happens with ultra faint dwarf galaxies (UFDW), which would be disrupted by the presence of PBHs. Making use of these effects, tight bounds have been obtained at for , weakening at lower masses down to for (). In a similar way, wide binary stellar systems may be perturbed by compact objects, potentially being disrupted after multiple encounters. The separation distribution of wide binaries restricts the PBH fraction from for down to at (Monroy-Rodríguez and Allen, 2014).

CMB

Radiation emitted either by accretion or from Hawking evaporation may affect the CMB spectrum in two ways: producing spectral distortions and modifying temperature anisotropies. The energetic radiation can enhance the ionization rate, delay recombination and shift the peaks of the CMB anisotropy spectrum, as well as induce more diffusion damping. The polarization spectrum could also be modified, since the increase of the free electron fraction would increase the Thomson optical depth and enhance the reionization bump at large angular scales. Although early CMB analyses (Ricotti et al., 2008) found very stringent bounds on the allowed abundance of accreting PBHs, later works revisited these computations and found much milder constraints (Horowitz, 2016; ). On another hand, while the former constraints rely on the assumption of spherical accretion, accreting disks have been argued to be more realistic for PBHs, resulting in tighter limits (Poulin et al., 2017b). Taking into account that PBHs could be immersed in DM halos with higher densities than the background, their accretion rates would be increased, also leading to more stringent constraints (Serpico et al., 2020). CMB limits from accretion are currently the most stringent ones for masses . The main caveat is their dependence on some details of the accretion mechanisms, such as the effective velocity and the accretion rate, which may not be very well understood yet.

On the other hand, the energy injection from PBHs evaporation would produce anisotropies and spectral distortions in the CMB spectrum, which would also limit the maximum abundance, leading to similar costraints to those obtained from the extra-galactic γ-ray background commented above (; Poulin et al., 2017a; ). Besides energy injection from accretion of BH evaporation, spectral distortions can also be produced by other means, such as the diffusion of photons due to Silk damping at small scales. This allows the translation of constraints on spectral distortions from FIRAS to stringent upper bounds on the PBH abundance, for masses (Nakama et al., 2018).

Lyα Forest

The discrete nature of PBHs would lead to a shot-noize contribution to the matter power spectrum, enchancing small scale fluctuations. Observations of the Lyα forest, which traces matter distribution at small scales, have been employed to extract limits on the maximum allowed fraction of PBHs (). The shot-noize contribution to the power spectrum depends on the joint product , for which the upper bound has been obtained (Murgia et al., 2019). The drawback of this method is on the priors of the reionization modeling and, as any Lyα forest analysis, is model dependent.

21 cm Cosmology

The 21 cm line signal from the hyperfine structure of the hydrogen is highly sensitive to the thermal state of the IGM, and thus, energy injection from PBH accretion or evaporation may leave strong observable signatures. The first claimed measurement of a global absorption dip by the EDGES collaboration () may lead to competitive bounds on the PBH abundance, either from accretion processes (Hektor et al., 2018) or from energy injection by evaporation (; Halder and Banerjee, 2021; Halder and Pandey, 2021). It must be noted, however, that the EDGES signal has not been confirmed yet by other experiments, and it has been argued that it could be explained by alternative mechanisms (Hills et al., 2018; ; Sims and Pober, 2020). On the other hand, although the cosmological 21 cm power spectrum has not been detected yet, forecasts with future experiments such as HERA and SKA have shown that 21 cm power spectrum data from these two radiotelescopes could potentially improve the bounds up to for masses above (Mena et al., 2019). The 21 cm forest observed as absorption troughs in the spectra of radio sources at could also provide similar limits on the PBH abundance, due to the Poisson shot noise and to the accretion heating effect (Villanueva-Domingo and Ichiki, 2021).

7 Conclusion

The extremely rich physics involved in the formation, evolution and distribution of PBHs implies a large number of observable effects which allow probing them. A myriad of constraints are present for a large range of PBHs masses. In recent years some of these limits, as those from microlensing, femtolensing, CMB accretion or BH mergers, have been revisited. Some of these bounds have been significantly weakened or have even disappeared, opening windows in the parameter space where PBHs could still form a substantial part of the DM, if not all. On the other hand, future experiments with better sensitivities may be able to reach yet unexplored regions in the parameter space and tighten up current limits. New probes, such as the 21 cm line pursued in radio interferometers like SKA, will present a promising and powerful way to proof or refute the existence of PBHs formed in the early universe and their potential contribution to the DM in the Universe.

Statements

Author contributions

PV-D is the main author of the text and the figures. OM and SP-R have substantially collaborated in the editing process.

Funding

PV and OM are supported by the Spanish grants FPA2017-85985-P and PROMETEO/2019/083. SP is supported by the Spanish FEDER/MCIU-AEI grant FPA2017-84543-P, and partially, by the Portuguese FCT (UID/FIS/00777/2019 and CERN/FIS-PAR/0004/2019). All authors acknowledge support from the European ITN project HIDDeN (H2020-MSCA-ITN-2019//860881-HIDDeN).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AbbottB. P.AbbottR.AbbottT.AbbottD.AbernathyM. R.AcerneseK.et al (2016). Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett.116, 061102. 10.1103/PhysRevLett.116.061102

  • 2

    AbbottB. P.AbbottR.AbbottT. D.AbbottS.AbrahamF.AcerneseK.et al (2019). Search for Subsolar Mass Ultracompact Binaries in Advanced LIGO’s Second Observing Run. Phys. Rev. Lett.123, 161102. 10.1103/PhysRevLett.123.161102

  • 3

    AbbottR.AbbottS.AbrahamF.AcerneseK.AckleyC.AdamsR. X.et al (2020). GW190521: a Binary Black Hole Merger with a Total Mass of 150 M⊙. Phys. Rev. Lett.125, 101102. 10.1103/PhysRevLett.125.101102

  • 4

    AcharyaS. K.KhatriR. (2020). CMB Spectral Distortions Constraints on Primordial Black Holes, Cosmic Strings and Long Lived Unstable Particles Revisited. JCAP02, 010. 10.1088/1475-7516/2020/02/010

  • 5

    AfshordiN.McDonaldP.SpergelD. N. (2003). Primordial Black Holes as Dark Matter: the Power Spectrum and Evaporation of Early Structures. ApJ594, L71L74. 10.1086/378763

  • 6

    AgolE.KamionkowskiM. (2002). X-rays from Isolated Black Holes in the Milky Way. Mon. Not. Roy. Astron. Soc.334, 553. 10.1046/j.1365-8711.2002.05523.x

  • 7

    AlcockC.AllsmanR. A.AlvesD. R.AxelrodT. S.BeckerA. C.BennettD. P.et al (2001). MACHO Project Limits on Black Hole Dark Matter in the 1-30 M Range. Astrophys. J. Lett.550, L169. 10.1086/319636

  • 8

    Ali-HaïmoudY. (2018). Correlation Function of High-Threshold Regions and Application to the Initial Small-Scale Clustering of Primordial Black Holes. Phys. Rev. Lett.121, 081304. 10.1103/PhysRevLett.121.081304

  • 9

    Ali-HaïmoudY.KamionkowskiM. (2017). Cosmic Microwave Background Limits on Accreting Primordial Black Holes. Phys. Rev.D95, 043534. 10.1103/PhysRevD.95.043534

  • 10

    Ali-HaïmoudY.KovetzE. D.KamionkowskiM. (2017). Merger Rate of Primordial Black-Hole Binaries. Phys. Rev. D96, 123523. 10.1103/PhysRevD.96.123523

  • 11

    ArbeyA.AuffingerJ. (2019). BlackHawk: a Public Code for Calculating the Hawking Evaporation Spectra of Any Black Hole Distribution. Eur. Phys. J. C79, 693. 10.1140/epjc/s10052-019-7161-1

  • 12

    AshoorioonA.RostamiA.FirouzjaeeJ. T. (2020). “Charting the Landscape in Our Neighborhood from the PBHs Mass Distribution and GWs,”.

  • 13

    AtalV.GermaniC. (2019). The Role of Non-gaussianities in Primordial Black Hole Formation. Phys. Dark Univ.24, 100275. 10.1016/j.dark.2019.100275

  • 14

    BallesterosG.Coronado-BlázquezJ.GaggeroD. (2020). X-ray and Gamma-ray Limits on the Primordial Black Hole Abundance from Hawking Radiation. Phys. Lett. B808, 135624. 10.1016/j.physletb.2020.135624

  • 15

    BallesterosG.SerpicoP. D.TaosoM. (2018). On the Merger Rate of Primordial Black Holes: Effects of Nearest Neighbours Distribution and Clustering. J. Cosmology Astropart. Phys.2018, 043. 10.1088/1475-7516/2018/10/043

  • 16

    BarnackaA.GlicensteinJ. F.ModerskiR. (2012). New Constraints on Primordial Black Holes Abundance from Femtolensing of Gamma-ray Bursts. Phys. Rev. D86, 043001. 10.1103/PhysRevD.86.043001

  • 17

    BellomoN.BernalJ. L.RaccanelliA.VerdeL. (2018). Primordial Black Holes as Dark Matter: Converting Constraints from Monochromatic to Extended Mass Distributions. J. Cosmology Astropart. Phys.2018, 004. 10.1088/1475-7516/2018/01/004

  • 18

    BelotskyK. M.DokuchaevV. I.EroshenkoY. N.EsipovaE. A.KhlopovM. Y.KhromykhL. A.et al (2019). Clusters of Primordial Black Holes. Eur. Phys. J. C79, 246. 10.1140/epjc/s10052-019-6741-4

  • 19

    BelotskyK. M.DmitrievA. D.EsipovaE. A.GaniV. A.GrobovA. V.et al (2014). Signatures of Primordial Black Hole Dark Matter. Mod. Phys. Lett. A.29, 1440005. 10.1142/S0217732314400057

  • 20

    BertoneG.HooperD. (2018). History of Dark Matter. Rev. Mod. Phys.90, 045002. 10.1103/RevModPhys.90.045002

  • 21

    BirdS.CholisI.MuñozJ. B.Ali-HaïmoudY.KamionkowskiM.et al (2016). Did LIGO Detect Dark Matter?Phys. Rev. Lett.116, 201301. 10.1103/PhysRevLett.116.201301

  • 22

    BoehmC.KobakhidzeA.O’hareC. A. J.PickerZ. S. C.SakellariadouM. (2021). Eliminating the LIGO Bounds on Primordial Black Hole Dark Matter. JCAP03, 078. 10.1088/1475-7516/2021/03/078

  • 23

    BondiH. (1952). On Spherically Symmetrical Accretion. MNRAS112, 195. 10.1093/mnras/112.2.195

  • 24

    BowmanJ. D.RogersA. E. E.MonsalveR. A.MozdzenT. J.MaheshN. (2018). An Absorption Profile Centred at 78 Megahertz in the Sky-Averaged Spectrum. Nature555, 6770. 10.1038/nature25792

  • 25

    BradleyR. F.TauscherK.RapettiD.BurnsJ. O. (2019). A Ground Plane Artifact that Induces an Absorption Profile in Averaged Spectra from Global 21 Cm Measurements, with Possible Application to EDGES. ApJ874, 153. 10.3847/1538-4357/ab0d8b

  • 26

    BrandtT. D. (2016). Constraints on MACHO Dark Matter from Compact Stellar Systems in Ultra-faint dwarf Galaxies. ApJ824, L31. 10.3847/2041-8205/824/2/L31

  • 27

    CaiR.-g.PiS.SasakiM. (2019). Gravitational Waves Induced by Non-gaussian Scalar Perturbations. Phys. Rev. Lett.122, 201101. 10.1103/PhysRevLett.122.201101

  • 28

    CarrB. J. (1975). The Primordial Black Hole Mass Spectrum. ApJ201, 119. 10.1086/153853

  • 29

    CarrB. J.KohriK.SendoudaY.YokoyamaJ. (2010). New Cosmological Constraints on Primordial Black Holes. Phys. Rev. D81, 104019. 10.1103/PhysRevD.81.104019

  • 30

    CarrB.KohriK.SendoudaY.YokoyamaJ. (2020). “Constraints on Primordial Black Holes,”.

  • 31

    CarrB.KühnelF. (2020). Primordial Black Holes as Dark Matter: Recent Developments. Annu. Rev. Nucl. Part. Sci.70, 355394. 10.1146/annurev-nucl-050520-125911

  • 32

    CarrB.KühnelF.SandstadM. (2016). Primordial Black Holes as Dark Matter. Phys. Rev. D94, 083504. 10.1103/PhysRevD.94.083504

  • 33

    CarrB.RaidalM.TenkanenT.VaskonenV.VeermäeH. (2017). Primordial Black Hole Constraints for Extended Mass Functions. Phys. Rev. D96, 023514. 10.1103/PhysRevD.96.023514

  • 34

    CarrB.SilkJ. (2018). Primordial Black Holes as Generators of Cosmic Structures. MNRAS478, 37563775. 10.1093/mnras/sty1204

  • 35

    ChaplineG. F. (1975). Cosmological Effects of Primordial Black Holes. Nat. J. C253, 251252. 10.1038/253251a0

  • 36

    ChenZ.-C.HuangQ.-G. (2020). Distinguishing Primordial Black Holes from Astrophysical Black Holes by Einstein Telescope and Cosmic Explorer. J. Cosmology Astropart. Phys.2020, 039. 10.1088/1475-7516/2020/08/039

  • 37

    ChibaT.YokoyamaS. (2017). Spin Distribution of Primordial Black Holes. PTEP2017, 083E01. 10.1093/ptep/ptx087

  • 38

    ChisholmJ. R. (2006). Clustering of Primordial Black Holes: Basic Results. Phys. Rev. D73, 083504. 10.1103/PhysRevD.73.083504

  • 39

    ClarkS.DuttaB.GaoY.MaY.-Z.StrigariL. E. (2018). 21 Cm Limits on Decaying Dark Matter and Primordial Black Holes. Phys. Rev. D98, 043006. 10.1103/PhysRevD.98.043006

  • 40

    ClarkS.DuttaB.GaoY.StrigariL. E.WatsonS. (2017). Planck Constraint on Relic Primordial Black Holes. Phys. Rev. D95, 083006. 10.1103/PhysRevD.95.083006

  • 41

    ClesseS.García-BellidoJ. (2015). Massive Primordial Black Holes from Hybrid Inflation as Dark Matter and the Seeds of Galaxies. Phys. Rev. D92, 023524. 10.1103/PhysRevD.92.023524

  • 42

    ClesseS.García-BellidoJ. (2018). Seven Hints for Primordial Black Hole Dark Matter. Phys. Dark Universe22, 137146. 10.1016/j.dark.2018.08.004

  • 43

    ClesseS.García-BellidoJ. (2017). The Clustering of Massive Primordial Black Holes as Dark Matter: Measuring Their Mass Distribution with Advanced LIGO. Phys. Dark Univ.15, 142147. 10.1016/j.dark.2016.10.002

  • 44

    CooganA.MorrisonL.ProfumoS. (2020). Direct Detection of Hawking Radiation from Asteroid-Mass Primordial Black Holes. Phys. Rev. Lett.126, 171101. 10.1103/PhysRevLett.126.171101

  • 45

    CroonD.McKeenD.RajN.WangZ. (2020). Subaru-HSC through a Different Lens: Microlensing by Extended Dark Matter Structures. Phys. Rev. D102, 083021. 10.1103/PhysRevD.102.083021

  • 46

    DasguptaB.LahaR.RayA. (2020). Neutrino and Positron Constraints on Spinning Primordial Black Hole Dark Matter. Phys. Rev. Lett.125, 101101. 10.1103/PhysRevLett.125.101101

  • 47

    De LucaV.DesjacquesV.FrancioliniG.MalhotraA.RiottoA. (2019a). The Initial Spin Probability Distribution of Primordial Black Holes. J. Cosmology Astropart. Phys.2019, 018. 10.1088/1475-7516/2019/05/018

  • 48

    De LucaV.DesjacquesV.FrancioliniG.RiottoA. (2020a). The Clustering Evolution of Primordial Black Holes. JCAP11, 028. 10.1088/1475-7516/2020/11/028

  • 49

    De LucaV.FrancioliniG.KehagiasA.PelosoM.RiottoA.ÜnalC. (2019b). The Ineludible Non-gaussianity of the Primordial Black Hole Abundance. J. Cosmology Astropart. Phys.2019, 048. 10.1088/1475-7516/2019/07/048

  • 50

    De LucaV.FrancioliniG.PaniP.RiottoA. (2021). Bayesian Evidence for Both Astrophysical and Primordial Black Holes: Mapping the GWTC-2 Catalog to Third-Generation Detectors, J. Cosmology Astroparticle Phys., 2021, 003. 10.1088/1475-7516/2021/05/003

  • 51

    De LucaV.FrancioliniG.PaniP.RiottoA. (2020b). Primordial Black Holes Confront LIGO/Virgo Data: Current Situation. JCAP06, 044. 10.1088/1475-7516/2020/06/044

  • 52

    De LucaV.FrancioliniG.PaniP.RiottoA. (2020c). The Evolution of Primordial Black Holes and Their Final Observable Spins. JCAP04, 052. 10.1088/1475-7516/2020/04/052

  • 53

    DengH.GarrigaJ.VilenkinA. (2017). Primordial Black Hole and Wormhole Formation by Domain walls. J. Cosmology Astropart. Phys.2017, 050. 10.1088/1475-7516/2017/04/050

  • 54

    DengH.VilenkinA. (2017). Primordial Black Hole Formation by Vacuum Bubbles. J. Cosmology Astropart. Phys.2017, 044. 10.1088/1475-7516/2017/12/044

  • 55

    DeRoccoW.GrahamP. W. (2019). Constraining Primordial Black Hole Abundance with the Galactic 511 keV Line. Phys. Rev. Lett.123, 251102. 10.1103/PhysRevLett.123.251102

  • 56

    DesjacquesV.RiottoA. (2018). Spatial Clustering of Primordial Black Holes. Phys. Rev. D98, 123533. 10.1103/PhysRevD.98.123533

  • 57

    DolgovA. D.KuranovA. G.MitichkinN. A.PoreyS.PostnovK. A.SazhinaO. S.et al (2020). On Mass Distribution of Coalescing Black Holes. JCAP12, 017. 10.1088/1475-7516/2020/12/017

  • 58

    EscrivàA.GermaniC.ShethR. K. (2021). Analytical Thresholds for Black Hole Formation in General Cosmological Backgrounds. JCAP01, 030. 10.1088/1475-7516/2021/01/030

  • 59

    FrancioliniG.KehagiasA.MatarreseS.RiottoA. (2018). Primordial Black Holes from Inflation and Non-gaussianity. J. Cosmology Astropart. Phys.2018. 01610.1088/1475-7516/2018/03/016

  • 60

    García-BellidoJ. (2017). Massive Primordial Black Holes as Dark Matter and Their Detection with Gravitational Waves. J. Phys. Conf. Ser.840, 012032. 10.1088/1742-6596/840/1/012032

  • 61

    Garcia-BellidoJ.ClesseS.FleuryP. (2017). LIGO Lo(g)Normal MACHO: Primordial Black Holes Survive SN Lensing Constraints. Phys. Dark Univ.20, 95100. 10.1016/j.dark.2018.04.005

  • 62

    Garcia-BellidoJ.Nuño SilesJ. F.Ruiz MoralesE. (2020). Bayesian Analysis of the Spin Distribution of LIGO/Virgo Black Holes. Phys. Dark Universe31, 100791. 10.1016/j.dark.2021.100791

  • 63

    García-BellidoJ.Ruiz MoralesE. (2017). Primordial Black Holes from Single Field Models of Inflation. Phys. Dark Universe18, 4754. 10.1016/j.dark.2017.09.007

  • 64

    GarrigaJ.VilenkinA.ZhangJ. (2016). Black Holes and the Multiverse. J. Cosmology Astropart. Phys.2016, 064. 10.1088/1475-7516/2016/02/064

  • 65

    GermaniC.MuscoI. (2019). Abundance of Primordial Black Holes Depends on the Shape of the Inflationary Power Spectrum. Phys. Rev. Lett.122, 141302. 10.1103/PhysRevLett.122.141302

  • 66

    GermaniC.ShethR. K. (2020). Nonlinear Statistics of Primordial Black Holes from Gaussian Curvature Perturbations. Phys. Rev. D101, 063520. 10.1103/PhysRevD.101.063520

  • 67

    GongJ.-O.KitajimaN. (2018). Distribution of Primordial Black Holes and 21cm Signature. JCAP1811, 041. 10.1088/1475-7516/2018/11/041

  • 68

    GongJ.-O.KitajimaN. (2017). Small-scale Structure and 21cm Fluctuations by Primordial Black Holes. JCAP1708, 017. 10.1088/1475-7516/2017/08/017

  • 69

    GreenA. M.KavanaghB. J. (2021). Primordial Black Holes as a Dark Matter Candidate. J. Phys. G48, 4. 10.1088/1361-6471/abc534

  • 70

    GreenA. M.LiddleA. R.MalikK. A.SasakiM. (2004). A New Calculation of the Mass Fraction of Primordial Black Holes. Phys. Rev. D70, 041502. 10.1103/PhysRevD.70.041502

  • 71

    GreenA. M. (2014). Primordial Black Holes: Sirens of the Early Universe. Quan. Aspects Black Holes178, 129149. 10.1007/978-3-319-10852-0_5

  • 72

    GriestK.CieplakA. M.LehnerM. J. (2014). Experimental Limits on Primordial Black Hole Dark Matter from the First 2 Yr of Kepler Data. Astrophys. J.786, 158. 10.1088/0004-637X/786/2/158

  • 73

    HalderA.BanerjeeS. (2021). Bounds on Abundance of Primordial Black Hole and Dark Matter from EDGES 21-cm Signal. Phys. Rev. D103, 063044. 10.1103/PhysRevD.103.063044

  • 74

    HalderA.PandeyM. (2021). Investigating the Effect of PBH, Dark Matter – Baryon and Dark Matter – Dark Energy Interaction on EDGES in 21cm Signal

  • 75

    HallA.GowA. D.ByrnesC. T. (2020). Bayesian Analysis of LIGO-Virgo Mergers: Primordial vs. Astrophysical Black Hole Populations. Phys. Rev. D102, 123524. 10.1103/PhysRevD.102.123524

  • 76

    HawkingS. (1989). Black Holes from Cosmic Strings. Phys. Lett. B231, 237239. 10.1016/0370-2693(89)90206-2

  • 77

    HawkingS. (1971). Gravitationally Collapsed Objects of Very Low Mass. MNRAS152, 75. 10.1093/mnras/152.1.75

  • 78

    HawkingS. W. (1974). Black Hole Explosions?Nat. J. C248, 3031. 10.1038/248030a0

  • 79

    HawkingS. W.MossI. G.StewartJ. M. (1982). Bubble Collisions in the Very Early Universe. Phys. Rev. D26, 26812693. 10.1103/PhysRevD.26.2681

  • 80

    HawkinsM. R. S. (2015). A New Look at Microlensing Limits on Dark Matter in the Galactic Halo. A&A575, A107. 10.1051/0004-6361/201425400

  • 81

    HektorA.HütsiG.MarzolaL.RaidalM.VaskonenV.VeermäeH. (2018). Constraining Primordial Black Holes with the EDGES 21-cm Absorption Signal. Phys. Rev. D98, 023503. 10.1103/PhysRevD.98.023503

  • 82

    HillsR.KulkarniG.MeerburgP. D.PuchweinE. (2018). Concerns about Modelling of the EDGES Data. Nat. J. C564, E32E34. 10.1038/s41586-018-0796-5

  • 83

    HorowitzB. (2016). Revisiting Primordial Black Holes Constraints from Ionization History

  • 84

    HütsiG.RaidalM.VaskonenV.VeermäeH. (2021). Two Populations of LIGO-Virgo Black Holes. JCAP03, 068. 10.1088/1475-7516/2021/03/068

  • 85

    IchimaruS. (1977). Bimodal Behavior of Accretion Disks - Theory and Application to Cygnus X-1 Transitions. Astrophys. J.214, 840855. 10.1086/155314

  • 86

    IguazJ.SerpicoP. D.SiegertT. (2021). Isotropic X-ray Bound on Primordial Black Hole Dark Matter

  • 87

    InmanD.Ali-HaïmoudY. (2019). Early Structure Formation in Primordial Black Hole Cosmologies. Phys. Rev. D100, 083528. 10.1103/PhysRevD.100.083528

  • 88

    KatzA.KoppJ.SibiryakovS.XueW. (2018). Femtolensing by Dark Matter Revisited. J. Cosmology Astropart. Phys.2018, 005. 10.1088/1475-7516/2018/12/005

  • 89

    KavanaghB. J.[Dataset](2019). Bradkav/Pbhbounds: Release version. 10.5281/zenodo.3538999

  • 90

    KavanaghB. J.GaggeroD.BertoneG. (2018). Merger Rate of a Subdominant Population of Primordial Black Holes. Phys. Rev. D98, 023536. 10.1103/PhysRevD.98.023536

  • 91

    KetovS. V.KhlopovM. Y. (2019). Cosmological Probes of Supersymmetric Field Theory Models at Superhigh Energy Scales. Symmetry11, 511. 10.3390/sym11040511

  • 92

    KhlopovM. Y. (2010). Primordial Black Holes. Res. Astron. Astrophys.10, 495528. 10.1088/1674-4527/10/6/001

  • 93

    LahaR.MuñozJ. B.SlatyerT. R. (2020). INTEGRAL Constraints on Primordial Black Holes and Particle Dark Matter. Phys. Rev. D101, 123514. 10.1103/PhysRevD.101.123514

  • 94

    LahaR. (2019). Primordial Black Holes as a Dark Matter Candidate Are Severely Constrained by the Galactic center 511 keV γ -ray Line. Phys. Rev. Lett.123, 251101. 10.1103/PhysRevLett.123.251101

  • 95

    LehmannB. V.ProfumoS.YantJ. (2018). The Maximal-Density Mass Function for Primordial Black Hole Dark Matter. JCAP04, 007. 10.1088/1475-7516/2018/04/007

  • 96

    LoprestoM. C. (2003). Some Simple Black Hole Thermodynamics. Phys. Teach.41, 299301. 10.1119/1.1571268

  • 97

    MenaO.Palomares-RuizS.Villanueva-DomingoP.WitteS. J. (2019). Constraining the Primordial Black Hole Abundance with 21-cm Cosmology. Phys. Rev.D100, 043540. 10.1103/PhysRevD.100.043540

  • 98

    MeszarosP. (1975). Primeval Black Holes and Galaxy Formation. A&A38, 513.

  • 99

    MirbabayiM.GruzinovA.NoreñaJ. (2020). Spin of Primordial Black Holes. JCAP03, 017. 10.1088/1475-7516/2020/03/017

  • 100

    Monroy-RodríguezM. A.AllenC. (2014). The End of the MACHO Era Revisited: New Limits on MACHO Masses from Halo Wide Binaries. ApJ790, 159. 10.1088/0004-637X/790/2/159

  • 101

    Moradinezhad DizgahA.FrancioliniG.RiottoA. (2019). Primordial Black Holes from Broad Spectra: Abundance and Clustering. J. Cosmology Astropart. Phys.2019, 001. 10.1088/1475-7516/2019/11/001

  • 102

    MrózP.UdalskiA.SkowronJ.PoleskiR.KozlowskiS.SzymanskiM. K.et al (2017). No Large Population of Unbound or Wide-Orbit Jupiter-mass Planets. Nat. J. C548, 183186. 10.1038/nature23276

  • 103

    MurgiaR.ScelfoG.VielM.RaccanelliA. (2019). Lyman-α forest Constraints on Primordial Black Holes as Dark Matter. Phys. Rev. Lett.123, 071102. 10.1103/PhysRevLett.123.071102

  • 104

    MuscoI.De LucaV.FrancioliniG.RiottoA. (2021). Threshold for Primordial Black Holes. II. A Simple Analytic Prescription. Phys. Rev. D103, 063538. 10.1103/PhysRevD.103.063538

  • 105

    MuscoI.MillerJ. C.RezzollaL. (2005). Computations of Primordial Black Hole Formation. Class. Quant. Grav.22, 14051424. 10.1088/0264-9381/22/7/013

  • 106

    NakamaT.CarrB.SilkJ. (2018). Limits on Primordial Black Holes from μ Distortions in Cosmic Microwave Background. Phys. Rev. D97, 043525. 10.1103/PhysRevD.97.043525

  • 107

    NarayanR.YiI. (1995). Advection Dominated Accretion: Underfed Black Holes and Neutron Stars. Astrophys. J.452, 710. 10.1086/176343

  • 108

    NarayanR.YiI. (1994). Advection-dominated Accretion: a Self-Similar Solution. ApJ428, L13. 10.1086/187381

  • 109

    NiemeyerJ. C.JedamzikK. (1999). Dynamics of Primordial Black Hole Formation. Phys. Rev. D59, 124013. 10.1103/PhysRevD.59.124013

  • 110

    NiemeyerJ. C.JedamzikK. (1998). Near-critical Gravitational Collapse and the Initial Mass Function of Primordial Black Holes. Phys. Rev. Lett.80, 54815484. 10.1103/PhysRevLett.80.5481

  • 111

    NiikuraH.TakadaM.YasudaN.LuptonR. H.SumiT.MoreS.et al (2019a). Microlensing Constraints on Primordial Black Holes with Subaru/HSC Andromeda Observations. Nat. Astron.3, 524534. 10.1038/s41550-019-0723-1

  • 112

    NiikuraH.TakadaM.YokoyamaS.SumiT.MasakiS. (2019b). Constraints on Earth-Mass Primordial Black Holes from OGLE 5-year Microlensing Events. Phys. Rev. D99, 083503. 10.1103/PhysRevD.99.083503

  • 113

    OguriM.DiegoJ. M.KaiserN.KellyP. L.BroadhurstT. (2018). Understanding Caustic Crossings in Giant Arcs: Characteristic Scales, Event Rates, and Constraints on Compact Dark Matter. Phys. Rev. D97, 023518. 10.1103/PhysRevD.97.023518

  • 114

    OppenheimerJ. R.VolkoffG. M. (1939). On Massive Neutron Cores. Phys. Rev.55, 374381. 10.1103/PhysRev.55.374

  • 115

    PaczynskiB. (1986). Gravitational Microlensing by the Galactic Halo. ApJ304, 1. 10.1086/164140

  • 116

    PageD. N. (1976). Particle Emission Rates from a Black Hole. II. Massless Particles from a Rotating Hole. Phys. Rev. D14, 32603273. 10.1103/PhysRevD.14.3260

  • 117

    PeacockJ. A. (1999). Cosmological Physics. Cambridge University Press. 10.1017/CBO9780511804533

  • 118

    PolnarevA. G.KhlopovM. Y. (1985). Cosmology, Primordial Black Holes, and Supermassive Particles. Sov. Phys. Usp.28, 213232. 10.1070/PU1985v028n03ABEH003858

  • 119

    PolnarevA.ZembowiczR. (1991). Formation of Primordial Black Holes by Cosmic Strings. Phys. Rev. D43, 11061109. 10.1103/PhysRevD.43.1106

  • 120

    PoulinV.LesgourguesJ.SerpicoP. D. (2017a). Cosmological Constraints on Exotic Injection of Electromagnetic Energy. JCAP03, 043. 10.1088/1475-7516/2017/03/043

  • 121

    PoulinV.SerpicoP. D.CaloreF.ClesseS.KohriK. (2017b). CMB Bounds on Disk-Accreting Massive Primordial Black Holes. Phys. Rev.D96, 083524. 10.1103/PhysRevD.96.083524

  • 122

    ReesM. J.PhinneyE. S.BegelmanM. C.BlandfordR. D. (1982). Ion Supported Tori and the Origin of Radio Jets. Nature295, 1721. 10.1038/295017a0

  • 123

    RicottiM.OstrikerJ. P.MackK. J. (2008). Effect of Primordial Black Holes on the Cosmic Microwave Background and Cosmological Parameter Estimates. ApJ680, 829845. 10.1086/587831

  • 124

    SaitoR.YokoyamaJ.NagataR. (2008). Single-field Inflation, Anomalous Enhancement of Superhorizon Fluctuations and Non-gaussianity in Primordial Black Hole Formation. J. Cosmology Astropart. Phys.2008, 024. 10.1088/1475-7516/2008/06/024

  • 125

    SandersR. H. (2010). The Dark Matter Problem: A Historical Perspective. Cambridge University Press.

  • 126

    SasakiM.SuyamaT.TanakaT.YokoyamaS. (2016). Primordial Black Hole Scenario for the Gravitational-Wave Wvent GW150914. Phys. Rev. Lett.117, 061101. 10.1103/PhysRevLett.117.061101

  • 127

    SasakiM.SuyamaT.TanakaT.YokoyamaS. (2018). Primordial Black Holes—Perspectives in Gravitational Wave Astronomy. Class. Quant. Grav.35, 063001. 10.1088/1361-6382/aaa7b4

  • 128

    SerpicoP. D.PoulinV.InmanD.KohriK. (2020). Cosmic Microwave Background Bounds on Primordial Black Holes Including Dark Matter Halo Accretion. Phys. Rev. Res.2, 023204. 10.1103/PhysRevResearch.2.023204

  • 129

    ShakuraN. I.SunyaevR. A. (1973). Black Holes in Binary Systems. Observational Appearance. Astron. Astrophys.24, 337355.

  • 130

    SimsP. H.PoberJ. C. (2020). Testing for Calibration Systematics in the EDGES Low-Band Data Using Bayesian Model Selection. MNRAS492, 2238. 10.1093/mnras/stz3388

  • 131

    SlatyerT. R. (2016). Indirect Dark Matter Signatures in the Cosmic Dark Ages II. Ionization, Heating and Photon Production from Arbitrary Energy Injections. Phys. Rev.D93, 023521. 10.1103/PhysRevD.93.023521

  • 132

    SuredaJ.MagañaJ.ArayaI. J.PadillaN. D. (2020). Press-Schechter Primordial Black Hole Mass Functions and Their Observational Constraints

  • 133

    SuyamaT.YokoyamaS. (2019). Clustering of Primordial Black Holes with Non-gaussian Initial Fluctuations. Prog. Theor. Exp. Phys.2019, 103E02. 10.1093/ptep/ptz105

  • 134

    ThakurtaS. N. G. (1981). Kerr Metric in an Expanding Universe. Indian J. Phys.55B, 304310.

  • 135

    TisserandP.Le GuillouL.AfonsoC.AlbertJ. N.AndersenJ.AnsariE.et al (2007). Limits on the MACHO Content of the Galactic Halo from the EROS-2 Survey of the Magellanic Clouds. Astron. Astrophys.469, 387404. 10.1051/0004-6361:20066017

  • 136

    TolmanR. C. (1939). Static Solutions of Einstein’s Field Equations for Spheres of Fluid. Phys. Rev.55, 364373. 10.1103/PhysRev.55.364

  • 137

    Villanueva-DomingoP.IchikiK. (2021). 21 Cm forest Constraints on Primordial Black Holes(Preprint) Arxiv

  • 138

    VolonteriM. (2010). Formation of Supermassive Black Holes. Astron. Astrophysics Rev.18, 279315. 10.1007/s00159-010-0029-x

  • 139

    WongK. W. K.FrancioliniG.De LucaV.BaibhavV.BertiE.PaniP.et al (2021). Constraining the Primordial Black Hole Scenario with Bayesian Inference and Machine Learning: the GWTC-2 Gravitational Wave Catalog. Phys. Rev. D103, 023026. 10.1103/PhysRevD.103.023026

  • 140

    XieF.-G.YuanF. (2012). The Radiative Efficiency of Hot Accretion Flows. Mon. Not. Roy. Astron. Soc.427, 1580. 10.1111/j.1365-2966.2012.22030.x

  • 141

    YokoyamaJ. (1998). Chaotic New Inflation and Formation of Primordial Black Holes. Phys. Rev. D58, 083510. 10.1103/PhysRevD.58.083510

  • 142

    YuanF.NarayanR. (2014). Hot Accretion Flows Around Black Holes. Ann. Rev. Astron. Astrophys.52, 529588. 10.1146/annurev-astro-082812-141003

  • 143

    Zel’dovichY. B.NovikovI. D. (1967). The Hypothesis of Cores Retarded during Expansion and the Hot Cosmological Model. Soviet Ast10, 602.

  • 144

    ZumalacárreguiM.SeljakU. (2018). Limits on Stellar-Mass Compact Objects as Dark Matter from Gravitational Lensing of Type Ia Supernovae. Phys. Rev. Lett.121, 141101. 10.1103/PhysRevLett.121.141101

Summary

Keywords

primordial black holes, dark matter, cosmology, accretion, 21 cm cosmology, gravitational waves, cosmic microwave background, microlensing

Citation

Villanueva-Domingo P, Mena O and Palomares-Ruiz S (2021) A Brief Review on Primordial Black Holes as Dark Matter. Front. Astron. Space Sci. 8:681084. doi: 10.3389/fspas.2021.681084

Received

15 March 2021

Accepted

11 May 2021

Published

28 May 2021

Volume

8 - 2021

Edited by

Fabrizio Nesti, University of L'Aquila, Italy

Reviewed by

Maxim Yurievich Khlopov, Southern Federal University, Russia

Vyacheslav Ivanovich Dokuchaev, Institute for Nuclear Research (RAS), Russia

Updates

Copyright

*Correspondence: Pablo Villanueva-Domingo,

This article was submitted to Cosmology, a section of the journal Frontiers in Astronomy and Space Sciences

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics