# Thermal Sunyaev-Zel’dovich anisotropy due to primordial black holes

###### Abstract

We investigate the thermal Sunyaev-Zel’dovich (SZ) effect caused by primordial black holes (PBHs) on the cosmic microwave background (CMB) temperature fluctuations. The gas accreting on a PBH heats up by the release of the gravitational energy. As a result, the heated gas in the vicinity of the PBH emits UV and x-ray photons. These photons can ionize and heat the intergalactic medium (IGM) around the PBH. Assuming the simple model of these emitting photons, we compute the profiles of the IGM ionization fraction and temperature around a PBH by using the numerical calculation of the radiative transfer. Using these profiles, we evaluate the Compton -parameter created by the IGM gas around a PBH. Finally, we estimate the CMB temperature angular power spectrum due to the PBH SZ effect in our model. We show that the SZ temperature anisotropy due to the PBHs has the flat angular power spectrum on small scale, and could dominate the primordial temperature spectrum on smaller scales than the Silk scale. This flat spectrum extends to the scale of the ionized region by the PBH emission. Therefore, in future SZ anisotropy measurements, the detection or nondetection of the flat spectrum gives useful information about the existence of PBHs.

## I Introduction

Many astronomical observations including the motion of stars in a galaxy, the large-scale structure distribution, and the cosmic microwave background (CMB) temperature anisotropy indicate the existence of dark matter (DM) in the Universe. However, its nature is still poorly understood. Currently, since the presence of DM can be probed only through its gravitational interaction, observational information about its nature is very limited. However, there are many works to study the nature of DM and many theoretical models have been suggested so far. Generally, the DM candidates are classified into two types; the nonbaryonic particle model and astrophysical compact object model. A lot of models for the nonbaryonic particle type are predicted in the physics beyond the standard model such as weakly interacting massive particles [1], axions [2, 3, 4, 5] and axionlike particles [6]. On the other hand, in the latter type, a primordial black hole (PBH) is the most potential candidate [7].

PBHs could have formed from high-density peaks in the very early Universe [8, 9]. In the inflation paradigm, the density fluctuations are generated from the quantum fluctuations. When an overdense region which exceeds a critical density threshold enters the horizon scale, the gravitational collapse of this region happens to form a PBH. The resultant mass of the formed PBH corresponds the horizon mass at the horizon-crossing epoch of the overdense region. As a result, the PBH mass range can span very widely. Therefore, the PBH abundance has been studied for a long time not only for DM but also as the relic of the primordial density fluctuations on small scales. Besides, the recent detection of gravitational wave (GW) events drew attention to PBHs as sources of GWs. The analysis of the observed GW data reveals that the detected GW events were produced by binary black hole mergers with masses larger than [10, 11, 12]. It would be difficult to produce a black hole with such large mass from stellar evolution in the standard solar metallicity environment [13]. On the other hand, the broad mass range of PBHs can cover the observed black hole masses. Therefore, the PBH abundance is also studied as a responsible source for the detected GW events [14, 15, 16].

The existence of PBHs affects various cosmological phenomena, depending on the PBH mass. For small mass PBHs, the abundance can be constrained by the effects of their evaporation. As first pointed out by Hawking [17], a black hole emits many kinds of particles with the thermal spectrum. As a result, PBHs with a mass smaller than g have evaporated by the present epoch. The abundance constraint on evaporated PBHs is obtained by investigating the effect of their evaporation on big bang nucleosynthesis [18], the CMB spectrum distortion [19], the recombination and reionization processes [20, 21], and the diffuse gamma-ray background [22, 23, 24].

For nonevaporated PBHs, the robust constraint is provided by gravitational lensing observations [25, 26, 27]. The black hole merger rate evaluated from the recent detection of GW events also constrains the PBH abundance [15, 16]. Recently several works focused on the gas accreting on massive PBHs. Because of the release of the gravitational energy during accretion, the gas becomes hot and emits x-ray and UV photons [28, 29, 30]. Resultantly the surrounding gas around a PBH is heated and ionized. Studying the cosmological effects of such heating and reionization provides the constraint on stellar mass PBHs with the recent CMB measurement [28, 29, 30, 31]. There are also some works suggesting that future 21-cm observations can probe these PBH heating and ionizing processes and give the strong constraint on PBHs in this mass range [32, 33].

In this paper, we study the thermal Sunyaev-Zel’dovich (SZ) effect due to PBHs. The thermal SZ effect is the distortion of the CMB energy spectrum through the inverse Compton scattering by high-energy electrons in hot plasma and this effect is parametrized by the Compton -parameter [34, 35]. As mentioned above, PBHs can heat and ionize the surrounding gas. The resultant ionized gas contributes not only to the global optical depth of the Thomson scattering as discussed in Refs. [28, 29, 30], but also to creating the nonzero Compton -parameter. The sky-averaged (global) Compton -paremter is already evaluated in Ref. [29]. They have obtained the global Compton -parameter induced by PBHs from evaluating the energy deposit of the PBH luminosity into a plasma. According to their estimation, the global Compton -parameter could be roughly , where describes the PBH abundance fraction to the total DM one. Since this estimated value is below the sensitivity of the Far-Infrared Absolute Spectrophotometer [36], they concluded that it is difficult to obtain a constraint on the PBH abundance from the global Compton -parameter with the current observation status.

Here, we focus on the spatial fluctuations of the Compton -parameter induced by PBHs and the resultant CMB temperature anisotropy due to them. The angular power spectrum of the CMB temperature anisotropy could provide additional information through its amplitude and shape depending on scales which we cannot access from the global value. Our main aim is to study the dependence of the SZ angular power spectrum on the PBH properties such as the abundance, mass and so on.

To evaluate the spatial fluctuations of the Compton -parameter, it is required to obtain the profiles of the ionization faction and the intergalactic medium (IGM) gas temperature around a PBH. Introducing the emission efficiency parameter of the PBH luminosity, we evaluate these profiles through the one-dimensional radiative transfer calculation. Then, assuming that the PBH distribution follows the dark matter density fluctuations, we calculate the CMB temperature anisotropy due to the SZ effect of PBHs. We investigate the dependence of the PBH fraction, the PBH mass and the emission efficiency parameter on the generated CMB anisotropy. From the comparison with the small-scale measurement of the CMB temperature anisotropy by South Pole Telescope (SPT) [37], we also discuss the suggestion of the SZ effect due to PBHs on the constraint on the PBH abundance with the emission efficiency parameter.

The rest of this paper is organized as follows. In Sec. II, we compute the gas temperature and ionization fraction around a PBH. Accordingly, we calculate the profile of the thermal SZ effect. In Sec. III, introducing the PBH fraction to dark matter, we evaluate the CMB temperature angular power spectrum due to the thermal SZ effect around PBHs. We also discuss the PBH abundance constraint from comparing our results with the SPT data. Finally, we summarize in Sec. IV. Throughout our paper, we take the flat CDM model with the Planck best fit parameters [38]: =.

## Ii Thermal Sunyaev-Zel’dovich Effect Due to a Pbh

In this section, we evaluate the gas temperature and ionization fraction of the IGM around a PBH, assuming the photon energy spectrum emitted from the hot gas in the vicinity of a PBH. Accordingly, we calculate the profile of the thermal SZ effect.

### ii.1 The luminosity from a PBH

Since a PBH creates a gravitational potential, the surrounding gas accretes on the PBH. During accretion, as it goes closer to a PBH, the gas becomes hot enough to emit x-ray and UV photons [28, 29, 30]. However, since the gas temperature highly depends on the astrophysical processes and environmental condition, there is a theoretical uncertainty in this luminosity. References [28, 29] evaluated the luminosity of a PBH considering spherical accretion on the PBH and showed that the luminosity is much lower than the Eddington luminosity. On the other hand, Ref. [30] pointed out that a PBH could have an accretion disk and that its luminosity is much higher than in the case of the spherical accretion. Moreover, when a massive PBH has an accretion disk, its luminosity might become nearly sub-Eddington luminosity similar to the active galactic nuclei as considered in Ref. [39].

Therefore, this emission efficiency can depend on the redshift, the PBH mass and other physical conditions around the PBH. Here, for simplicity, we introduce one free parameter for the emission efficiency, , which represents the total PBH luminosity, , in terms of the Eddington luminosity.

(1) |

where is the Eddington luminosity for mass

(2) |

### ii.2 IGM temperature and ionization profiles around a PBH

Now we consider the profiles of the ionization fraction and gas temperature around a PBH with . For simplicity, we consider only hydrogen as the IGM gas component.

The evolutions of the ionization fraction, and temperature in the IGM are given by the following equations,

(5) | |||

(6) |

where , , , , , and are the number density of hydrogen nucleus, the Boltzmann constant, the IGM gas adiabatic index, , the mean molecular weight, the proton mass, the gas mass density, and the case B recombination rate given in Ref. [41] respectively. In Eqs. (5) and (6), , and are the ionization, heating and cooling rates, respectively. Here, we adopt cooling mechanisms due to recombination cooling, collisional ionization cooling, collisional excitation cooling and Compton cooling [41]. For simplicity, we assume that the IGM gas is homogeneous and we set the IGM gas density to the cosmological mean (background) values at a given time .

We consider only photons emitted from a PBH as the ionization and heating source. We solve Eqs. (5) and (6) with the spherically symmetric assumption. Since we evaluate the ionization and temperature evolutions for the cosmological timescale, it is useful to introduce the comoving radial distance from a PBH for representing the spatial position, instead of the physical distance. The ionization and heating rates at the comoving radial distance can be written, respectively, as

(7) | ||||

(8) |

where , and are the Planck constant, the neutral hydrogen number density given in , and the absorbed cross section area of ionization photons, , respectively. In the above equations, represents the photon energy flux for a frequency at a comoving distance :

(9) |

Here is the scale factor normalized as at the present epoch, , and is the optical depth of HI gas from the central PBH to the comoving distance ,

(10) |

Solving Eqs. (5) and (6) consistently by a spherically symmetric radiative transfer code, we obtain the evolutions of and at each comoving distance and cosmic time .

Figure. 1 shows the radial profile of the ionization fraction as a function of the physical distance, , from a PBH. In the figure, we set . The difference of colors represents different redshifts. As the Universe evolves, the ionized region extends outward. This expansion is proportional to . This behavior is consistent with the redshift evolution of the Strömgren radius with the constant photon flux from a PBH. The Strömgren radius (in physical) can be related to the number rate of photons emitted from a PBH, , and has the redshift dependence as in

(11) |

where is obtained from

(12) |

Since we assume that the photon emission from a PBH does not have the redshift evolution, the redshift dependence in Eq. (11) comes from the (physical) number density of electrons and protons, .

Figure 2 provides the redshift evolution of the IGM gas temperature. We plot the radial profile of the temperature as a function of the physical radial distance from a PBH. The color difference shows the difference of the redshifts. Similarly to the ionization fraction in FIG. 1, the lower the redshift becomes, the more the volume of the heated region increases. In the fully ionized central region, the temperature profile is flat and the redshift dependence of the amplitude is very weak. Since the heating mechanism is photoionization, the heating becomes efficient in the neutral region. Therefore, the peak of the gas temperature locates at the spatial position where the ionization fraction becomes about . Thus, the heated region extends more than the ionized region. In a sufficiently distant region, the gas temperature becomes the background IGM temperature which is represented by the thin dashed lines in Fig. 2.

At the end of this subsection, we comment on the dependence of the ionization fraction and gas temperature on the PBH luminosity, . According to Eq. (11) with Eq. (12), the Strömgren radius depends on . As shown in Figs. 1 and 2, the profiles of and is sensitive to the evolution of the Strömgren radius. Therefore, as or increases, the ionization front expands towards outside roughly proportional to . However, since the temperature is determined by the balance between the heating and cooling, the amplitude of the temperature at the fully ionized region does not have the simple dependence on .

### ii.3 Compton -parameter induced by a PBH

Figures. 1 and 2 show that a PBH can make a hot gas plasma around itself. Passing through this plasma, CMB photons suffer from the SZ effect. As a result, the observed brightness temperature of the CMB is shifted by at a frequency from the background CMB temperature, [34]. At the comoving distance from the PBH on the sky, this observed shift of the brightness temperature can be written with the Compton -parameter,

(13) |

where is the frequency spectral function of the SZ effect,

(14) |

The Compton -parameter at the comoving distance on the sky, , can be obtained from the integral form along a line-of-sight with the impact parameter from the PBH,

(15) |

where is the Thomson cross section, is the electron mass and is the speed of light. In Eq. (15), is the comoving distance projected on the line-of-sight direction and is the comoving radial distance from the PBH satisfying, .

Using Eq. (15) with the profiles of the ionization fraction and temperature shown in Figs. 1 and 2, we can compute the -parameter. We plot the -parameter due to a PBH with as a function of the physical distance, , for different redshifts in Fig. 3. As expected, the -parameter profile depends on the ionization fraction profile. Therefore, when becomes larger than the Strömgren radius, the -parameter quickly falls down. Following the evolution of the Strömgren radius, as the Universe evolves, the flat region in the -parameter profile extends outward. In Eq. (15), most of the contribution to the integration come from the region inside the Strömgren radius (the fully ionized region).

The amplitude of -parameter decreases as the Universe evolves because the redshift dependence of the electron number density dominates that of the integral distance and gas temperature.

We also comment on the dependence of the -parameter profile on and . As mentioned above, the -parameter profile strongly depends on the Strömgren radius. Since the Strömgren radius is proportional to , the parameters and can affect the -parameter profile. As or becomes large, the amplitude of the -parameter increases and the tail of the profile moves toward large .

Based on Fig. 3, we can roughly estimate the total -parameter induced by PBHs on the line-of-sight direction between and . Figure 3 shows that one PBH with and provides . Since the comoving number density of such PBHs is about where is the PBH fraction to the total DM abundance, the resultant typical distance of PBHs is roughly . Since the radial distance between and is about , the PBH number on the line-of-sight is Therefore, the total -parameter integrated over the line-of-sight is . Consequently, this y-parameter can generate a SZ signature on the order of K. In the next section, we evaluate the anisotropic signal in more detail.

## Iii THERMAL SUNYAEV-ZEL’DOVICH ANISOTROPY DUE TO PBHs

As shown in the previous section, a hot plasma is generated around a PBH and can cause the SZ effect. If PBHs account for a significant fraction of the dark matter abundance, they can produce the observable CMB temperature anisotropy through the SZ effect. In this section, we evaluate the angular power spectrum of the CMB temperature fluctuations by the SZ effect. We also discuss the constraint on the PBH abundance from the comparison with the small-scale CMB measurement from the SPT.

### iii.1 The CMB temperature angular power spectrum by the SZ effect

To calculate the angular power spectrum of the CMB temperature due to the SZ effect caused by PBHs, we take a similar method to the one for galaxy clusters based on the halo formalism [42, 43]. Accordingly, the angular power spectrum can be described as the sum of the two contributions:

(16) |

where is the “one-PBH” term describing the Poisson contribution and is the “two-PBH” term arising from clustering of PBHs. Assuming that the PBH mass function is restricted to a single mass , we can write these two terms as

(17) | ||||

(18) |

where is the comoving volume, is the comoving distance to the redshift , is the matter power spectrum and is the comoving number density of PBHs with mass . Here we assume that PBHs contribute to the DM abundance with the fraction . Therefore, can be written as with the present critical density of the Universe, . We integrate these equations from the redshift to for whose values we will add small discussion later.

In Eqs. (17) and (18), taking the small angle approximation, we get as the two-dimension Fourier transform of the Compton -parameter for PBH mass at the redshift obtained in Sec. II.3:

(19) |

where is a vector describing the two-dimensional Fourier mode with , the comoving distance on the sky is given in and is the angular direction on the sky sphere.

Calculating Eqs. (17) and (18) with the profile of , we can obtain the CMB temperature anisotropy induced by the SZ effect due to PBHs. In Fig. 4, we plot the obtained angular power spectra of the CMB temperature anisotropy for the different PBH parameter sets . Our PBH parameter sets are summarized in Table 1. For comparison, we show the primordial CMB temperature anisotropy as the black dashed-dotted line. We also plot the SPT data with the error bars [37].

The shape of the spectrum is independent on the PBH parameters. We find out that the contribution from the two-PBH term dominates the one-PBH term. Since the typical scale of the -parameter profile in Fig. 3 is the Strömgren radius, , which is much smaller than the CMB observation scales, is constant on the CMB observation scales. Therefore, according to Eq. (18), the shape of the angular spectrum is determined by the matter power spectrum. In fact, independently on the PBH parameters, Fig. 4 shows that the spectrum has a flat shape on larger multipoles than . Because of this flat shape, the SZ temperature anisotropy due to PBHs can dominate the primordial temperature anisotropy on smaller scales than the Silk scale.

On the other hand, the amplitude depends on the PBH parameters. Because of the integration range in Eq. (15), is proportional to . Figure 3 tells us that the typical scale of non-negligible is also . Therefore, according to Eq. (19), the angular Fourier component of the -parameter, , is proportional to . The PBH number density is . As a result, has the dependence as and we found that the angular power spectrum on the flat region can be approximated in

(20) |

Figure 4 clearly shows this dependence of to us. Because of the degeneracy between and , the angular spectrum for the case C is coincident with that for the case D.

We also investigate the redshift contribution to the angular spectrum. The observed temperature anisotropy is a resultant effect which is integrated along the line-of-sight direction. As already mentioned, is proportional to which becomes large as the Universe evolves. In particular, we find out that the contribution from the redshift higher than is less than 1 %. On the other hand, in the low-redshift side, the Universe is gradually reionized and finally fully ionized until [38]. During the epoch of reionization, the regions ionized by PBHs are caught up in the reionization process by stars and galaxies. As a result, in lower redshifts around the epoch of reionization, the SZ effect due to PBHs becomes fainter as the reionization proceeds. Here, for simplicity, we take into account the PBH SZ effect between and in the integration of Eqs. (17) and (18).

In Ref. [29], they calculated the typical value of in the analytic model with the spherical accretion on a PBH. Their obtained is on the order of for the PBH mass at . Substituting this value into Eq. (20), we can see that the signal is much smaller than the current observation sensitivity such as the SPT even for . However, there exists the theoretical uncertainty in as mentioned in Sec. II. Depending on the environment around a PBH, could become larger than the above one. In that case, we might get the observable value in the present observation status. Therefore, in the next subsection, we discuss the possible constraint on the PBH abundance from the thermal SZ angular power spectrum with the free parameter .

(A) | 10 | ||
---|---|---|---|

(B) | 100 | ||

(C) | 1000 | ||

(D) | 10 | ||

(E) | 100 | ||

(F) | 1000 |

### iii.2 Application to the PBH constraint

Now we discuss the possible constraint on the PBH abundance from the CMB SZ power spectrum measurement. When we fix the PBH mass and the emission efficiency , the amplitude of the SZ spectrum is determined by the PBH abundance in our model.

Since the SZ effect due to PBHs can induce the CMB temperature anisotropy on small scales, the small-scale CMB measurement provides the constraint on the PBH abundance. The SPT data reach the minimum value, at . Therefore, the condition, at gives the limit on with fixing and , because of the flat shape of the spectrum. We found that the constraint from the SZ angular power spectrum is roughly independent on the PBH mass and given in

(21) |

for . Therefore, when we adopt , for , as in Refs. [28, 29], the current SPT CMB measurement on small scales Eq. (21) cannot provide the constraint on the PBH abundance. When the PBH has larger than as discussed in Refs. [30, 39], we can obtain the abundance constraint from the SPT data. However, according to Refs. [28, 29], we found that the CMB measurement of the global Thomson optical depth provides the PBH abundance constraint in terms of . Therefore, we conclude that the constraint from the SZ measurement is always weaker than the one from the global Thomson optical depth.

## Iv Conclusion

In this paper, we have investigated the CMB temperature angular power spectrum due to the SZ effect caused by PBHs. The gas accreting on a PBH heats up because of the release of the gravitational energy. As a result, the heated gas in the vicinity of the PBH emits the UV and x-ray photons. These photons can ionize and heat the IGM around the PBH. The ionized hot IGM causes the secondary CMB temperature anisotropy through the SZ effect.

First assuming the luminosity of photons emitted in the vicinity of a PBH with a free parameter, we have evaluated the profiles of the IGM ionization fraction and temperature around a PBH by solving the one-dimensional radiative transfer equations. Based on these profiles, we have obtained the SZ Compton -parameter profile around a PBH. Following the halo formalism, finally, we have calculated the CMB temperature angular power spectrum due to the PBH SZ effect with assuming the PBH abundance.

We have shown that the SZ spectrum due to PBHs could dominate the primordial temperature spectrum on smaller scales than the Silk scale. Since the correlation on such scales is made by distortions between two different PBHs, the shape of the angular power spectrum depends on the matter power spectrum. We have also found that the amplitude of the spectrum is sensitive to the PBH abundance and the emission efficiency of the gas accreting on PBHs. In the case of a scale-invariant primordial spectrum, the angular power spectrum has a flat shape on small scales . Using in Ref. [29], we cannot get the observable signal in current observations such as the SPT. However, has a theoretical uncertainty and Refs.[30, 39] suggest that PBHs could have larger than . In this case, PBHs could make the observable signature.

We also investigated the impact for the constraint of the PBH abundance from the SZ effect. Our obtained constraint on the PBH fraction to dark matter is for the PBH mass range . However, we should note that the CMB constraint from the optical depth provides a much tighter constraint with the same . Therefore we conclude that the SZ measurement does not give a new constraint on the PBH abundance.

Although the constraint from the SZ angular power spectrum is weaker than that from the optical depth, it is worth mentioning the impact of the future small-scale SZ measurements on the probe of the existence of PBHs. We have shown that the CMB temperature spectrum of the SZ effect due to PBHs has the flat spectrum on smaller scales than . This flat spectrum extends up to the scale of the ionized region by the PBH emission, e.g., roughly 1 kpc in the physical scale for with . The SZ effect in galaxy clusters can produce large CMB anisotropies on small scales. However, its spectrum has a peak around and the amplitude decays on higher . Therefore, the detection or nondetection of the flat SZ spectrum on higher than gives useful information about the existence of PBHs.

In this work, we have assumed the luminosity from the hot gas in the vicinity of a PBH with introducing a free parameter for the emission efficiency. References. [28, 29, 30] have studied the emission efficiency with the simple assumptions analytically. However, to obtain the emission efficiency consistently, we need to solve the dynamics of gas accreting on a PBH with considering the backreaction of the luminosity from hot gas. In the next study, we address this issue by numerical simulations and investigate the cosmological effects of the emission from gas accreting on a PBH in more detail.

###### Acknowledgements.

We are grateful to Kenji Hasegawa for his useful comments. This work is supported by JSPS KAKENHI Grants No. 15K17646 (H.T.) and No. 17H01110 (H.T.).## References

- [1] G. Steigman and M. S. Turner. Cosmological constraints on the properties of weakly interacting massive particles. Nuclear Physics B, 253:375–386, 1985.
- [2] R. D. Peccei and H. R. Quinn. Constraints imposed by CP conservation in the presence of pseudoparticles. Phys. Rev. D, 16:1791–1797, September 1977.
- [3] R. D. Peccei and H. R. Quinn. CP conservation in the presence of pseudoparticles. Physical Review Letters, 38:1440–1443, June 1977.
- [4] Steven Weinberg. A new light boson? Phys. Rev. Lett., 40:223–226, Jan 1978.
- [5] F. Wilczek. Problem of strong and invariance in the presence of instantons. Phys. Rev. Lett., 40:279–282, Jan 1978.
- [6] Y. Chikashige, R.N. Mohapatra, and R.D. Peccei. Are there real goldstone bosons associated with broken lepton number? Physics Letters B, 98(4):265 – 268, 1981.
- [7] B. Carr, F. Kühnel, and M. Sandstad. Primordial black holes as dark matter. Phys. Rev. D, 94(8):083504, October 2016.
- [8] Y. B. Zel’dovich and I. D. Novikov. The Hypothesis of Cores Retarded during Expansion and the Hot Cosmological Model. Soviet Ast., 10:602, February 1967.
- [9] S. Hawking. Gravitationally collapsed objects of very low mass. MNRAS, 152:75, 1971.
- [10] B. P. Abbott, R. Abbott, T. D. Abbott, and et al. Properties of the Binary Black Hole Merger GW150914. Physical Review Letters, 116(24):241102, June 2016.
- [11] B. P. Abbott, R. Abbott, T. D. Abbott, and et al. Gw170104: Observation of a 50-solar-mass binary black hole coalescence at redshift 0.2. Phys. Rev. Lett., 118:221101, Jun 2017.
- [12] B. P. Abbott, R. Abbott, T. D. Abbott, and et al. Gw170814: A three-detector observation of gravitational waves from a binary black hole coalescence. Phys. Rev. Lett., 119:141101, Oct 2017.
- [13] K. Belczynski, D. E. Holz, T. Bulik, and R. O’Shaughnessy. The first gravitational-wave source from the isolated evolution of two stars in the 40-100 solar mass range. Nature (London), 534:512–515, June 2016.
- [14] S. Bird, I. Cholis, J. B. Muñoz, Y. Ali-Haïmoud, M. Kamionkowski, E. D. Kovetz, A. Raccanelli, and A. G. Riess. Did LIGO Detect Dark Matter? Physical Review Letters, 116(20):201301, May 2016.
- [15] M. Sasaki, T. Suyama, T. Tanaka, and S. Yokoyama. Primordial Black Hole Scenario for the Gravitational-Wave Event GW150914. Physical Review Letters, 117(6):061101, August 2016.
- [16] S. Clesse and J. García-Bellido. The clustering of massive Primordial Black Holes as Dark Matter: Measuring their mass distribution with advanced LIGO. Physics of the Dark Universe, 15:142–147, March 2017.
- [17] S. W. Hawking. Particle creation by black holes. Communications in Mathematical Physics, 43:199–220, August 1975.
- [18] B. J. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama. New cosmological constraints on primordial black holes. Phys. Rev. D, 81(10):104019, May 2010.
- [19] H. Tashiro and N. Sugiyama. Constraints on primordial black holes by distortions of the cosmic microwave background. Phys. Rev. D, 78(2):023004, July 2008.
- [20] S. J. Clark, B. Dutta, Y. Gao, L. E. Strigari, and S. Watson. Planck constraint on relic primordial black holes. Phys. Rev. D, 95(8):083006, April 2017.
- [21] S. J. Clark, B. Dutta, Y. Gao, Y.-Z. Ma, and L. E. Strigari. 21 cm limits on decaying dark matter and primordial black holes. Phys. Rev. D, 98(4):043006, August 2018.
- [22] D. N. Page and S. W. Hawking. Gamma rays from primordial black holes. Astrophys. J. , 206:1–7, May 1976.
- [23] B. J. Carr. Some cosmological consequences of primordial black-hole evaporations. Astrophys. J. , 206:8–25, May 1976.
- [24] J. H. MacGibbon and B. J. Carr. Cosmic rays from primordial black holes. Astrophys. J. , 371:447–469, April 1991.
- [25] P. Tisserand, L. Le Guillou, C. Afonso, J. N. Albert, J. Andersen, R. Ansari, É. Aubourg, P. Bareyre, J. P. Beaulieu, X. Charlot, C. Coutures, R. Ferlet, P. Fouqué, J. F. Glicenstein, B. Goldman, A. Gould, D. Graff, M. Gros, J. Haissinski, C. Hamadache, J. de Kat, T. Lasserre, É. Lesquoy, C. Loup, C. Magneville, J. B. Marquette, É. Maurice, A. Maury, A. Milsztajn, M. Moniez, N. Palanque-Delabrouille, O. Perdereau, Y. R. Rahal, J. Rich, M. Spiro, A. Vidal-Madjar, L. Vigroux, S. Zylberajch, and EROS-2 Collaboration. Limits on the Macho content of the Galactic Halo from the EROS-2 Survey of the Magellanic Clouds. A&A, 469:387–404, July 2007.
- [26] H. Niikura, M. Takada, N. Yasuda, R. H. Lupton, T. Sumi, S. More, T. Kurita, S. Sugiyama, A. More, M. Oguri, and M. Chiba. Microlensing constraints on primordial black holes with the Subaru/HSC Andromeda observation. arXiv e-prints, January 2017.
- [27] M. Zumalacárregui and U. Seljak. Limits on Stellar-Mass Compact Objects as Dark Matter from Gravitational Lensing of Type Ia Supernovae. Physical Review Letters, 121(14):141101, October 2018.
- [28] M. Ricotti, J. Ostriker, and K. Mack. Effect of Primordial Black Holes on the Cosmic Microwave Background and Cosmological Parameter Estimates. page U5.004, January 2017.
- [29] Y. Ali-Haïmoud and M. Kamionkowski. Cosmic microwave background limits on accreting primordial black holes. Phys. Rev. D, 95(4):043534, February 2017.
- [30] Vivian Poulin, Pasquale D. Serpico, Francesca Calore, Sebastien Clesse, and Kazunori Kohri. CMB bounds on disk-accreting massive primordial black holes. Phys. Rev., D96(8):083524, 2017.
- [31] J. L. Bernal, N. Bellomo, A. Raccanelli, and L. Verde. Cosmological implications of primordial black holes. J. Cosmology Astropart. Phys., 10:052, October 2017.
- [32] H. Tashiro and N. Sugiyama. The effect of primordial black holes on 21-cm fluctuations. MNRAS, 435:3001–3008, November 2013.
- [33] J.-O. Gong and N. Kitajima. Small-scale structure and 21cm fluctuations by primordial black holes. J. Cosmology Astropart. Phys., 8:017, August 2017.
- [34] Y. B. Zeldovich and R. A. Sunyaev. The Interaction of Matter and Radiation in a Hot-Model Universe. Ap&SS, 4:301–316, July 1969.
- [35] R. A. Sunyaev and Y. B. Zeldovich. The Observations of Relic Radiation as a Test of the Nature of X-Ray Radiation from the Clusters of Galaxies. Comments on Astrophysics and Space Physics, 4:173, November 1972.
- [36] D. J. Fixsen, E. S. Cheng, J. M. Gales, J. C. Mather, R. A. ShaFer, and E. L. Wright. The cosmic microwave background spectrum from the full cobe firas data set. The Astrophysical Journal, 473(2):576, 1996.
- [37] R. Chown, Y. Omori, K. Aylor, B. A. Benson, L. E. Bleem, and et al. Maps of the Southern Millimeter-wave Sky from Combined 2500 deg SPT-SZ and Planck Temperature Data. ApJS, 239:10, November 2018.
- [38] Planck Collaboration and N. et al. Aghanim. Planck 2018 results. VI. Cosmological parameters. ArXiv e-prints, July 2018.
- [39] José Luis Bernal, Alvise Raccanelli, Licia Verde, and Joseph Silk. Signatures of primordial black holes as seeds of supermassive black holes. Journal of Cosmology and Astro-Particle Physics, 2018:017, May 2018.
- [40] D. Giannios, N. D. Kylafis, and D. Psaltis. Spectra and time variability of Galactic black-hole X-ray sources in the low/hard state. A&A, 425:163–169, October 2004.
- [41] M. Fukugita and M. Kawasaki. Reionization during Hierarchical Clustering in a Universe Dominated by Cold Dark Matter. MNRAS, 269:563, August 1994.
- [42] E. Komatsu and T. Kitayama. Sunyaev-Zeldovich Fluctuations from Spatial Correlations between Clusters of Galaxies. ApJ, 526:L1–L4, November 1999.
- [43] Uroš Seljak. Analytic model for galaxy and dark matter clustering. MNRAS, 318:203–213, October 2000.