Issue 
A&A
Volume 618, October 2018



Article Number  A139  
Number of page(s)  7  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201833483  
Published online  30 October 2018 
Constraints on primordial black hole dark matter from Galactic center Xray observations
^{1}
National Institute of Chemical Physics and Biophysics, Rävala 10, 10143 Tallinn, Estonia
email: andi.hektor@cern.ch, martti.raidal@cern.ch
^{2}
Tartu Observatory, University of Tartu, Tõravere, 61602 Tartumaa, Estonia
email: gert.hutsi@to.ee
Received:
23
May
2018
Accepted:
23
July
2018
Context. Surprisingly high masses of the black holes inferred from the Laser Interferometer GravitationalWave Observatory (LIGO) and Virgo gravitational wave measurements have lead to speculations that the observed mergers might be due to 𝒪(10) M_{⊙} primordial black holes (PBHs). Furthermore, it has been suggested that the whole amount of dark matter (DM) might be in that exotic form.
Aims. We investigate constraints on the PBH DM using NuSTAR Galactic center (GC) Xray data.
Methods. We used a robust Monte Carlo approach in conjunction with a radiatively inefficient PBH accretion model with commonly accepted model parameters. Compared to previous studies we allowed for multiple forms of DM density profiles. Most importantly, our study includes treatment of the gas turbulence, which significantly modifies the relative velocity between PBHs and gas.
Results. We show that inclusion of the effects of gas turbulence and the uncertainties related to the DM density profile reduces significantly the gas accretion onto PBHs compared to the claimed values in previous papers. It is highly improbable to obtain accreting PBHs brighter than the NuSTAR point source limit using observationally determined gas velocities.
Conclusions. One can safely conclude that GC Xray observations cannot rule out 𝒪(10) M_{⊙} PBH DM.
Key words: cosmology: theory / dark matter / black hole physics
© ESO 2018
1. Introduction
The first direct detection of gravitational waves by LIGO^{1} collaboration and several subsequent detections by LIGO and Virgo^{2}, have revolutionized observational astronomy by opening a new window into the cosmos. Thus far five binary black hole (BH) and a single binary neutron star mergers have been observed (Abbott et al. 2016a,b, Abbott et al. 2017a,b,c). Although the measurements have been fully consistent with general relativity, and thus have severely constrained modified gravity models, the large inferred BH masses have been somewhat surprising^{3}.
This has initiated vigorous speculations that the observed mergers might be between 𝒪(10) M_{⊙} primordial BHs (PBHs; Hawking 1971; Carr & Hawking 1974; Carr 1975; Meszaros 1975; Chapline 1975) and, furthermore, that the whole dark matter (DM) might be in the form of PBHs in that mass range (Kashlinsky 2016; Bird et al. 2016; Clesse & GarcíaBellido 2016; Sasaki et al. 2016; for a review see, e.g. Sasaki et al. 2018). However, several recent works have shown that this interpretation is strongly constrained. In the case of DM in the form of 𝒪(10) M_{⊙} PBHs, the increased twobody interaction rate has potentially significant impact on dynamics of the stellar systems (Brandt 2016; Koushiappas & Loeb 2017). Also the gravitational lensing gets modified (Mediavilla et al. 2017; Diego et al. 2017; Oguri et al. 2018; Zumalacarregui & Seljak 2018). The present constraints on PBH abundances from various lensing experiments and dynamical observations are rederived in Carr et al. (2017) for various nonmonochromatic PBH mass functions. Those constraints must be supplemented with bounds arising from nonobservation of the gravitational wave background (Raidal et al. 2017) and from the measured LIGO rate (AliHaïmoud et al. 2017; AliHaïmoud 2018).
Thus the accumulated experimental data collectively constrain the fraction of PBH DM, f_{PBH}, to be below unity, barring scenarios in which the radiation induced by PBH is strongly modified (Raidal et al. 2018) or the results of lensing experiments are misinterpreted (GarcíaBellido & Clesse 2018; Clesse & GarcíaBellido 2017; GarcíaBellido et al. 2017).
The aim of this work is to revisit the constraints on PBH abundance arising from the photon flux created by accreting PBHs in our Galaxy. Our approach, based on the Galactic measurements, is complementary to those previously adopted in the literature to constrain the PBH abundance from the cosmic microwave background observations (Ricotti et al. 2008; Horowitz 2016; AliHaimoud & Kamionkowski 2017; Poulin et al. 2017) and from recent global 21 cm measurements (Hektor et al. 2018). The authors of Gaggero et al. (2017) have claimed very strong constraints on 𝒪(10) M_{⊙} PBH DM using Xray and radio observations of the Galactic center (GC). In this paper we reconsider the GC Xray constraints on 𝒪(10) M_{⊙} PBH DM taking into account physical effects that were overlooked in the latter work. Compared to Gaggero et al. (2017) we include in our analysis turbulent gas motions inside molecular clouds and also allow for uncertainties in the DM density profile of our Galactic halo. We show that inclusion of those ingredients, in particular the measured velocities of the gas motion, strongly suppress the PBH accretion and remove the previously claimed strong bounds arising from the NuSTAR GC Xray data.
Bounds on PBHs based on accretion arguments using Xray data have also been found in Inoue & Kusenko (2017). There the authors use the extragalactic luminosity function of Xray binaries as determined by Mineo et al. (2012) and demand that the accreting PBH population should not overshine this astrophysical component.
The paper is structured as follows. In Sect. 2 we briefly sketch our physical assumptions and calculation method, our main results are presented in Sect. 3, and the discussion and summary are in Sect. 4.
2. Model details and parameter assumptions
Black holes are only visible in electromagnetic (EM) radiation if they accrete a sufficient amount of baryonic material. The only EMvisible stellar mass BHs we know are in binary systems, where the companion star can provide sufficient mass transfer. However, owing to the relative diluteness of 𝒪(10) M_{⊙} PBH DM and thus rather negligible twobody scattering with stars in Milky Way (MW) size galaxies, it is usually not an option for PBHs to get incorporated into binary systems with ordinary star as a companion. Thus PBHs have to accrete matter directly from the interstellar medium (ISM).
In this case it is relevant to apply the Bondi accretion model (Bondi 1952) as a useful starting point. However, the Bondi mass accretion rate Ṁ_{B} has to be reduced by a factor λ to be consistent with the nonobservation of significant population of isolated neutron stars^{4}. The present upper bound on λ is ∼10^{−3} − 10^{−2} (Perna et al. 2003). To obtain model for the emitted EM radiation, the description of the mass accretion rate has to be augmented with a description of the radiative efficiency. As common in cases with low mass accretion rate and small opacity, we assume an advection dominated accretion flow (ADAF) model (Narayan & Yi 1994). The radiative efficiency η (bolometric luminosity L_{bol} = ηṀc^{2}) is often approximated as(1)
where ṁ is the mass accretion rate in units of Eddington rate, i.e. ṁ ≡ Ṁ/Ṁ_{Edd}, and ṁ_{crit} ≃ 0.01 (e.g. Narayan & McClintock 2008). Taking Ṁ = λṀ_{B} and assuming the above radiative efficiency we write the accreting PBH luminosity as^{5} (2)
where M_{BH} is a BH mass, f is a fraction of energy going to a selected energy band, n_{H} is the hydrogen number density, and v_{eff} is the total relative velocity of a BH with respect to the gas. This can be modeled to contain various components: (i) motion of the BHs, (ii) thermal motion of the gas, and (iii) turbulent motion of the gas. We also denote a dangerous divergence due to L ∝ v^{−6}, which needs careful treatment in any realistic physical system.
The highest chance of observing a bright PBH is toward a dense molecular region of the GC, or the socalled central molecular zone (CMZ), where conditions for both the high number density of PBHs and a dense surrounding medium are satisfied. Because of the measured coldness of the dense medium, we can safely neglect thermal motions of the gas to estimate v_{eff} in Eq. (2). However, turbulent gas motions in CMZ are far from negligible. For example, the 1D velocity dispersions of 2.6–53 km s^{−1} with the median value of 9.8 km s^{−1} have been inferred in Henshaw et al. (2016). Under these assumptions v_{eff} in Eq. (2) can be approximated as . The BH velocity v_{BH} is assumed to follow the Maxwell–Boltzmann (MB) distribution, where the MB scale parameter (1D velocity dispersion) is obtained by solving spherically symmetrized Jeans equation (e.g. Binney & Tremaine 1987) with two components: the baryonic bulge and the DM halo.
We assume an isotropic velocity distribution, i.e. σ_{r} = σ_{θ} = σ_{ϕ}. The PBH velocity dispersion is obtained by solving(3)
where the tracer density is that of the PBHs, n_{PBH} (∝ρ_{DM}), but the mass inside radius r should contain all the assumed components, M_{tot}(< r)=M_{DM}(< r)+M_{bulge}(< r).
To be more precise, a subdominant contribution from the baryonic disk should also be included. For simplicity, we neglect it in our study, which makes our values for v_{BH} somewhat underestimated and thus more conservative. The density profile of the DM halo is allowed to have two analytic forms: NavarroFrenkWhite (NFW; Navarro et al. 1997) and Einasto (Einasto 1965). For the GC gas distribution we use an analytic fitting form from Ferrière et al. (2007). There the analytic spatial model for the mean gas densities is provided, but ISM is supposed to have a broad hierarchy of densities. In the case where v_{turb} = 0, it turns out that the Bondi accretion in conjunction with the low velocity tail of the MB distribution has a particular scaling property: the number of bright accreting PBHs above a fixed luminosity is independent of the hierarchy^{6}. However, if v_{turb} > 0 this scaling property gets broken and for more reliable treatment there should be a model for the smallscale gas density distribution. We performed our model calculations for two separate cases where the smallscale gas density distribution is assumed, first to follow a uniform distribution within a resolution element with the mean density given by the analytic model of Ferrière et al. (2007), and second has a powerlaw probability distribution function with a mean given by the model of Ferrière et al. (2007). We take β = 2.8, which is the value typical for the giant molecular clouds, e.g. Berkhuijsen et al. (1999); Agol & Kamionkowski (2002).
3. Main results
It is well known that because of the low level of contaminating backgrounds and absorption, a search for accreting compact objects is especially efficient in the Xray band of the EM spectrum. As such we compare our PBH DM model predictions with the NuSTAR ^{7} GC survey, which has lead to a detection of ∼70 Xray point sources. The NuSTAR GC survey has point source limits 4 × 10^{33} and 8 × 10^{33} erg s^{−1} in the 3 − 10 keV and 10 − 40 keV bands, respectively (Hong et al. 2016). Alternatively, it would be possible to use Chandra observations of the GC (Muno et al. 2009), which has approximately an order of magnitude higher sensitivity in softer Xray band 0.5 − 8 keV. However, in softer bands, the number of Xray sources is steeply rising (mostly due to cataclysmic variables), and thus it is significantly easier to hide a subdominant accreting PBH population. Thus in the following we present our results based on the NuSTAR observations.
Figure 1 shows an approximate mask we used for the NuSTAR GC survey derived from the results presented in Hong et al. 2016). In Fig. 2 we show the ADAF model spectra (M_{BH} = 10 M_{⊙}) for different values of the specific accretion rate ṁ ≡ Ṁ/Ṁ_{Edd} taken from Yuan & Narayan 2014 (2014; see Fig. 1 therein). The fractions of energy going to the soft and hard NuSTAR bands (see the gray shaded regions) are also shown for each of the model curves. We note that the ADAF model depends on various input parameters such as (i) viscosity parameter, (ii) magnetization parameter, (iii) electron heating parameter, and (iv) wind parameter, which in this particular case have values α = 0.1, β = 9, δ = 0.5, and s = 0.4, respectively. In case a jet with a nonthermal population of electrons appears or if the hot flow itself contains a nonthermal component, the prominent inverse Compton bumps of Fig. 2 get significantly smoothed out. For example, such a nonthermal electron component is required for the ADAF model to provide a satisfactory fit to the available SgrA* data (Yuan et al. 2003).
Fig. 1. Approximate spatial mask of the NuSTAR GC survey shown in Galactic coordinates (l, b). 

Open with DEXTER 
Fig. 2. Spectra from ADAF model for various values of the specific accretion rate ṁ ≡ Ṁ/Ṁ_{Edd} (Yuan & Narayan 2014). The parameters f denote the fractions of energy going to the soft and hard NuSTAR bands shown with the grayshaded regions. 

Open with DEXTER 
Figure 3 shows the PBH DM velocity dispersions derived from spherically symmetric Jeans equation. The three solid lines show the Galactic model consisting only of (i) NFW DM halo, (ii) NFW DM halo + baryonic bulge, and (iii) Einasto DM halo + baryonic bulge. The profile parameters relevant for the Galactic DM halo are taken from Cirelli et al. (2011). For the baryonic bulge an analytic density distribution from McMillan (2011) is assumed, while for simplicity spherically symmetric approximation is taken by replacing the scale radius r_{0} with the following geometric mean ; the Galactic bulge has an approximate aspect ratio of q ≃ 1/2. The value r_{cut} is also replaced similarly. For comparison, the triangles and squares in Fig. 3 denote the DM velocity dispersions taken from Ferrer & Hunter (2013), in which more detailed modeling is performed. Compared to Ferrer & Hunter (2013) we assumed spherical symmetry and completely ignored contribution from the baryonic disk, which can be seen to lead to ∼25% underestimation of the DM velocity dispersion, i.e. a factor of ∼3 overestimation of the accreting PBH luminosity. Also, replacement of NFW with the Einasto halo leads to a factor of ∼1.5 increase in the velocity dispersion in the very central region. However, at extremities of the CMZ the effect is far less pronounced, but it is important to keep in mind that, according to Eq. (2), velocity rises as L ∝ v^{−6}; the approximate size of the CMZ is shown as a light gray shaded region.
Even though the velocity dispersions for the NFW and Einasto profiles shown in Fig. 3 at r ∼ 0.1 kpc differ only by ∼10%, it turns out that the more central regions, where the velocity dispersions and the DM density profiles differ more strongly, are responsible for very large differences we obtain between the NFW and Einasto cases (see Figs. 4 and 5).
Fig. 4. Distributions for the number of bright PBHs above the NuSTAR point source detection limit for different levels of turbulence in the CMZ (0–5 km s^{−1}) and for different DM profiles (NFW, Einasto). Each distribution is derived from 1000 independent Monte Carlo realizations. The orange histograms represent the model with only a bulge and DM halo. The green histograms correspond to the case in which the MB 1D velocity dispersion has been increased by 25% to correct for the missing Galactic disk contribution (see Fig. 3). The uniform gas distribution within resolution element is assumed. 

Open with DEXTER 
Fig. 5. Analog of Fig. 4 for the model with a smallscale, powerlaw gas density distribution. 

Open with DEXTER 
Now we have all the ingredients available to proceed with our main calculations. We ran several Monte Carlo simulations, each containing 1000 independent realizations, the results of which are shown in Figs. 4 and 5. We assume that all the DM is in the form of PBHs. To speed up our computations we implemented a gridbased Monte Carlo by dividing the lineofsight cone, which covers the full NuSTAR GC survey footprint into small cubic cells of size ΔV = 2 × 2 × 2 pc^{3}. We looped over the cells and calculate the number of PBHs in each of these by drawing numbers from the Poisson distribution with mean , where ρ_{DM} is the assumed DM density and M_{PBH} the PBH mass. For each object drawn this way we generated their 3D velocity components by drawing numbers from the Gaussian with zero mean and dispersion taken from Eq. (3). By using a finer grid we confirm that the above cell size is sufficient for obtaining converging results.
In case the smallscale gas density distribution is taken to be a power law, , we draw the relevant number density as follows: where . The value u is a random variable following the standard uniform distribution, i.e. u ∼ U(0, 1), ⟨n_{H}⟩ is the mean gas density taken from Ferrière et al. (2007); the powerlaw index is assumed to be β = 2.8.
The results of our Monte Carlo calculations for the model with uniform gas density within the resolution element are shown in Fig. 4 in which the left and righthand panels correspond to the NFW and Einasto profiles, respectively. The red histograms represent the model with only the bulge and DM halo. The blue histograms correspond to the case in which the MB 1D velocity dispersion was increased by 25% to correct for the missing disk contribution (see Fig. 3). To facilitate a comparison with the results of Gaggero et al. (2017) the analogous green and yellow histograms assume instead a fixed spectral factor of f = 0.3. From top to bottom the effective velocity for turbulent gas motions is allowed to vary in the range 0 − 8 km s^{−1} with 2 km s^{−1} step size. In these calculations, PBH mass of 30 M_{⊙} is assumed. The numbers of bright objects should be compared with the NuSTAR point source count of ∼70. Only in the case of NFW halo plus negligible turbulent gas motions, is it possible to obtain (depending on a particular model) at most ∼10 − 30 PBHs above the NuSTAR point source limit. More realistically, stochastic treatment for turbulent gas motions should also be included, i.e. the resulting effective distribution should be calculated as a weighted average over the distributions shown in Fig. 4. In case turbulent gas motions are modeled with MB distribution with a realistic dispersion of 10 km s^{−1}, the total weight for the distributions of Fig. 4 (velocity range 0 − 8 km s^{−1}) is only ∼14%. As a result, in this case PBHs brighter than the NuSTAR limit are only rarely obtained, and it is safe to conclude, contrary to the claims made in Gaggero et al. (2017), that the GC Xray observations cannot rule out 𝒪(10) M_{⊙} PBH DM.
The analogous results for the model with a powerlaw gas density distribution are shown in Fig. 5. In order to ease comparison we kept the axes scales identical to those used in Fig. 4. As expected, the results in case of small v_{turb} are practically the same as obtained for the model with a constant smallscale gas density. For v_{turb} ≳ 5 km s^{−1} the results start to differ more visibly, namely, we obtain somewhat larger number of bright objects. However, this increase is very moderate, and thus our previous conclusions remain unchanged.
The above modeling has completely ignored possible twobody relaxation effects due to the granularity of the DM distribution. According to Koushiappas & Loeb (2017), in the case of dwarf galaxies these effects are not negligible. However, it turns out that the relevant relaxation time for the MW size halo inside the central region comparable to the size of the CMZ is significantly larger than the Hubble time. Thus, the twobody scattering effects can be safely neglected. The relaxation only becomes noticeable in the very central (few tens of parsecs) region. In this region PBHs, which are significantly more massive than typical stars, would start to migrate to the center. The corresponding density profile would steepen by trying to approach the isothermal ρ ∝ r^{−2}. However this does not increase the chances of obtaining bright PBHs, since the central subsystem would heat up and the v^{−6} scaling of Eq. (2) would dominate above the rise from the steepening of the density profile.
4. Discussion and summary
In this paper we investigated how well the Galactic Xray measurements are able to constrain 𝒪(10) M_{⊙} PBH DM. The probability of seeing bright accreting PBHs is highest for the GC, where the conditions for large PBH number density and high gas density, are simultaneously met. As such, we used the GC data from the NuSTAR Xray telescope in our analysis.
A similar study has previously been performed by Gaggero et al. (2017). Compared to this work, rather than having a single fixed form for the DM density profile, we investigate how much the results change by allowing a reasonable level of profile variability. The other, arguably the most important difference, is that our model includes treatment for the ISM gas turbulence. In addition, our treatment for the spectral factor f of Eq. (2) is more complete: we calculate f selfconsistently using the ADAF spectral templates taken from Yuan & Narayan (2014), whereas Gaggero et al. (2017) assume a fixed value f = 0.3.
In Gaggero et al. (2017) 𝒪(10) M_{⊙} PBH DM is claimed to be ruled out by up to 40σ using the GC Xray data, i.e. this would correspond to around ∼1500 Xray visible PBHs. In comparison, only in our most optimistic case with negligible gas turbulence and the NFW density profile, we obtain on average ∼10 − 15 Xray visible PBHs above the NuSTAR point source limit. For almost all the other cases in which a reasonable level of gas turbulence is allowed, hardly any bright PBHs are found. Thus the extraordinary strong claims made in Gaggero et al. (2017) are artifacts of the unphysical v^{−6} singularity in the PBH luminosity estimation arising from their assumption that the accreted gas is standing still. We demonstrate that the inclusion of measured gas velocities removes this singularity, and as a result refutes their claims regarding strong bounds on the PBH abundance.
We also note that the use of the experimental upper bound for the accretion efficiency parameter of Eq. (2), i.e. λ ∼ 10^{−2}, as done in Gaggero et al. (2017) in order to rule out PBH DM, does not sound correct. For this purpose actually a lower bound is needed. The latter cannot be reliably estimated since the accretion physics in a highly turbulent and magnetized ISM has very large uncertainties.
In this paper we opted to use a specific ADAF accretion model taken from Yuan & Narayan (2014) to model the energy distribution of the emitted radiation. To compare our results directly with the results presented in Gaggero et al. (2017) we also used a simple model in which 30% of the emitted radiation is assumed to fall into NuSTAR Xray energy band. We stress that beyond ADAFs, where most of the energy gained by the accreting gas is simply carried beyond the BH event horizon, there are many other possible models for the accretion from the relatively dilute ISM. Such models include adiabatic inflowoutflow solution (ADIOS) and convectiondominated accretion flow (CDAF) type solutions (Blandford & Begelman 1999; Quataert & Gruzinov 2000), which are inefficient because most of the accreting mass is simply driven out, and thus the mass actually gained by the BH is much smaller than the largescale mass accretion scale as estimated by Ṁ_{B}. Also depending on a particular model, the resulting energy distributions for the emitted radiation can vary substantially. It is not the topic of this paper to study all these possibilities. The main message we would like to convey is that under these very large model uncertainties, together with the absence of the relevant phenomenology for the isolated stellarmass BHs accreting from the ISM, it is impossible to convincingly rule out a possible existence of the 𝒪(10) M_{⊙} PBH DM by accretion arguments alone. In general we agree that under the specific assumptions made in Gaggero et al. (2017) it is possible to obtain a measurable population of Xray sources above the NuSTAR sensitivity limit; however the above result is very different when allowing for different DM profiles and taking the significant level of turbulence of the gas in the CMZ into account. To conclude, it is fair to say that with the current state of knowledge about the isolated stellarmass BHs accreting from the ISM, the present GC Xray observations cannot rule out 𝒪(10) M_{⊙} PBH DM.
Considering the complexity of star formation and binary stellar evolution, and the fact that this is the first time we are able to probe BHBH systems directly, one should not immediately draw too strong conclusions; see, e.g. Belczynski et al. (2016); Fishbach & Holz (2017); Dvorkin et al. (2018); Lamberts et al. (2018).
We note that in the context of the accretion flow described by the ADAF model, it is natural to try to use a population of isolated neutron stars to gain access to the possible values for the λ parameter. In particular, neutron stars are almost as compact as BHs. These stars have, for their limited range of masses, radii that are only a factor of ∼3 larger than those of the wouldbe BHs with similar masses (i.e. comparable to the size of the last stable orbit of the nonrotating BH). Compared to BHs, however they have a surface that lights up, and thus provides a direct way to probe the inflow of advected material.
We have assumed that we are always in a regime where the mass accretion rate is below 1% of the Eddington’s and thus radiative efficiency is given by the first line of Eq. (1).
Assuming independence of the velocity and gas density distributions (f_{MB}(v) and f(n_{H}), respectively) and that the PBH luminosity follows a deterministic relation given by Eq. (2), we can write the joint probability density function as , where κ is a constant factor and δ is the Dirac delta function. By integrating over n_{H} we obtain the joint luminosity and velocity distribution ℱ(L, v) from which the luminosity probability function can be obtained as as . To a good approximation f_{MB}(v) is given as , since even for the very dense molecular cloud cores with n_{H} ∼ 10^{5} cm^{−3} the threshold velocity below which the luminosity exceeds the NuSTAR limit is v ∼ 35 km s^{−1}, i.e. comfortably smaller than the 1D velocity dispersion shown in Fig. 3. As such, it is easy to see that f(L) indeed does not depend on gas density distribution, . The number of bright objects in volume ΔV follows as N = n_{PBH}ΔVf(> L_{NuSTAR}), i.e. indeed independent of the gas number density distribution function.
Fig. 3. Velocity dispersion of the PBH DM. The solid curves show our estimate from the Jeans equation. Three cases are shown: the Galaxy consisting only of (i) NFW DM halo, (ii) NFW DM halo + baryonic bulge, and (iii) Einasto DM halo + baryonic bulge. The triangles and squares denote the DM velocity dispersions taken from Ferrer & Hunter (2013). Neglecting the baryonic disk leads to a ∼25% underestimation of the DM velocity dispersion, i.e. a factor of ∼3 overestimation of the accreting PBH luminosity. Replacing NFW with the Einasto density profile leads to a factor of ∼1.5 increase in the velocity dispersion in the very central parts. The light gray shaded region shows an approximate extent of the CMZ. 

Open with DEXTER 
Acknowledgments
This work was supported by the grants IUT236, IUT262, PUT808, and by EU through the ERDF CoE program grant TK133 and by the Estonian Research Council via the Mobilitas Plus grant MOBTT5. AH thanks the Horizon 2020 program as the project has received funding from the program under the Marie SklodowskaCurie grant agreement No 661103.
References
 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016a, Phys. Rev. Lett., 116, 241103 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016b, Phys. Rev. Lett., 116, 061102 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017a, ApJ, 851, L35 [NASA ADS] [CrossRef] [Google Scholar]
 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017b, Phys. Rev. Lett., 119, 141101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017c, Phys. Rev. Lett., 119, 161101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Agol, E., & Kamionkowski, M. 2002, MNRAS, 334, 553 [NASA ADS] [CrossRef] [Google Scholar]
 AliHaïmoud, Y. 2018, Phys. Rev. Lett., 121, 081304 [NASA ADS] [CrossRef] [Google Scholar]
 AliHaimoud, Y., & Kamionkowski, M. 2017, Phys. Rev. D, 95, 043534 [NASA ADS] [CrossRef] [Google Scholar]
 AliHaïmoud, Y., Kovetz, E. D., & Kamionkowski, M. 2017, Phys. Rev. D, 96, 123523 [NASA ADS] [CrossRef] [Google Scholar]
 Belczynski, K., Holz, D. E., Bulik, T., & O’Shaughnessy, R. 2016, Nature, 534, 512 [NASA ADS] [CrossRef] [Google Scholar]
 Berkhuijsen, E. M. 1999, in Plasma Turbulence and Energetic Particles in Astrophysics, ed. M. Ostrowski, 61 [Google Scholar]
 Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton University Press). [Google Scholar]
 Bird, S., Cholis, I., Muñoz, J. B., et al. 2016, Phys. Rev. Lett., 116, 201301 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Begelman, M. C. 1999, MNRAS, 303, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Bondi, H. 1952, MNRAS, 112, 195 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Brandt, T. D. 2016, ApJ, 824, L31 [NASA ADS] [CrossRef] [Google Scholar]
 Carr, B. J. 1975, ApJ, 201, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Carr, B. J., & Hawking, S. W. 1974, MNRAS, 168, 399 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Carr, B., Raidal, M., Tenkanen, T., Vaskonen, V., & Veermäe, H. 2017, Phys. Rev. D, 96, 023514 [NASA ADS] [CrossRef] [Google Scholar]
 Chapline, G. F. 1975, Nature, 253, 251 [NASA ADS] [CrossRef] [Google Scholar]
 Cirelli, M., Corcella, G., Hektor, A., et al. 2011, Phys. Dark Univ., 3, 051 [Google Scholar]
 Clesse, S., & GarcíaBellido, J. 2016, Phys. Dark Univ., 10, 002 [Google Scholar]
 Clesse, S., & GarcíaBellido, J. 2017, ArXiv eprints [arXiv:1711.10458] [Google Scholar]
 Diego, J. M., Kaiser, N., Broadhurst, T., et al. 2017, ApJ, 857, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Dvorkin, I., Uzan, J.P., Vangioni, E., & Silk, J. 2018, MNRAS, 479, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Einasto, J. 1965, Trudy Astrofizicheskogo Instituta AlmaAta, 5, 87 [NASA ADS] [Google Scholar]
 Ferrer, F., & Hunter, D. R. 2013, Phys. Dark Univ., 9, 005 [Google Scholar]
 Ferrière, K., Gillard, W., & Jean, P. 2007, A&A, 467, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fishbach, M., & Holz, D. E. 2017, ApJ, 851, L25 [NASA ADS] [CrossRef] [Google Scholar]
 Gaggero, D., Bertone, G., Calore, F., et al. 2017, Phys. Rev. Lett., 118, 241101 [NASA ADS] [CrossRef] [Google Scholar]
 GarcíaBellido, J., & Clesse, S. 2018, Phys. Dark Univ., 19, 144 [CrossRef] [Google Scholar]
 GarcíaBellido, J., Clesse, S., Fleury, P. 2017, ArXiv eprints [arXiv:1712.06574] [Google Scholar]
 Hawking, S. 1971, MNRAS, 152, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Hektor, A., Hütsi, G., Marzola, L., et al. 2018, Phys. Rev. D, 98, 023503 [NASA ADS] [CrossRef] [Google Scholar]
 Henshaw, J. D., Longmore, S. N., Kruijssen, J. M. D., et al. 2016, MNRAS, 457, 2675 [NASA ADS] [CrossRef] [Google Scholar]
 Hong, J., Mori, K., Hailey, C. J., et al. 2016, ApJ, 825, 132 [NASA ADS] [CrossRef] [Google Scholar]
 Horowitz, B. 2016, ArXiv eprints [arXiv:1612.07264]. [Google Scholar]
 Inoue, Y., & Kusenko, A. 2017, Phys. Dark Univ., 10, 034 [Google Scholar]
 Kashlinsky, A. 2016, ApJ, 823, L25 [NASA ADS] [CrossRef] [Google Scholar]
 Koushiappas, S. M., & Loeb, A. 2017, Phys. Rev. Lett., 119, 041102 [NASA ADS] [CrossRef] [Google Scholar]
 Lamberts, A., GarrisonKimmel, S., Hopkins, P., et al. 2018, MNRAS, 480, 2704 [NASA ADS] [CrossRef] [Google Scholar]
 McMillan, P. J. 2011, MNRAS, 414, 2446 [NASA ADS] [CrossRef] [Google Scholar]
 Mediavilla, E., JiménezVicente, J., Muñoz, J. A., VivesArias, H., & CalderónInfante, J. 2017, ApJ, 836, L18 [NASA ADS] [CrossRef] [Google Scholar]
 Meszaros, P. 1975, A&A, 38, 5 [NASA ADS] [Google Scholar]
 Mineo, S., Gilfanov, M., & Sunyaev, R. 2012, MNRAS, 419, 2095 [NASA ADS] [CrossRef] [Google Scholar]
 Muno, M. P., Bauer, F. E., Baganoff, F. K., et al. 2009, ApJS, 181, 110 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & McClintock, J. E. 2008, New Astron. Rev., 51, 733 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I.S. 1994, ApJ, 428, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Oguri, M., Diego, J. M., Kaiser, N., Kelly, P. L., & Broadhurst, T. 2018, Phys. Rev. D, 97, 023518 [NASA ADS] [CrossRef] [Google Scholar]
 Perna, R., Narayan, R., Rybicki, G., Stella, L., & Treves, A. 2003, ApJ, 594, 936 [NASA ADS] [CrossRef] [Google Scholar]
 Poulin, V., Serpico, P. D., Calore, F., Clesse, S., & Kohri, K. 2017, Phys. Rev. D, 96, 083524 [NASA ADS] [CrossRef] [Google Scholar]
 Quataert, E., & Gruzinov, A. 2000, ApJ, 539, 809 [NASA ADS] [CrossRef] [Google Scholar]
 Raidal, M., Vaskonen, V., & Veermäe, H. 2017, JCAP, 1709, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Raidal, M., Solodukhin, S., Vaskonen, V., & Veermäe, H. 2018, Phys. Rev. D, 97, 123520 [NASA ADS] [CrossRef] [Google Scholar]
 Ricotti, M., Ostriker, J. P., & Mack, K. J. 2008, ApJ, 680, 829 [NASA ADS] [CrossRef] [Google Scholar]
 Sasaki, M., Suyama, T., Tanaka, T., & Yokoyama, S. 2016, Phys. Rev. Lett., 117, 061101 [NASA ADS] [CrossRef] [Google Scholar]
 Sasaki, M., Suyama, T., Tanaka, T., & Yokoyama, S. 2018, Class. Quant. Grav., 35, 063001 [NASA ADS] [CrossRef] [Google Scholar]
 Yuan, F., & Narayan, R. 2014, Ann. Rev. Astron. Astrophys., 52, 529 [NASA ADS] [CrossRef] [Google Scholar]
 Yuan, F., Quataert, E., & Narayan, R. 2003, ApJ, 598, 301 [NASA ADS] [CrossRef] [Google Scholar]
 Zumalacarregui, M., & Seljak, U. 2018, Phys. Rev. Lett., 121, 141101 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1. Approximate spatial mask of the NuSTAR GC survey shown in Galactic coordinates (l, b). 

Open with DEXTER  
In the text 
Fig. 2. Spectra from ADAF model for various values of the specific accretion rate ṁ ≡ Ṁ/Ṁ_{Edd} (Yuan & Narayan 2014). The parameters f denote the fractions of energy going to the soft and hard NuSTAR bands shown with the grayshaded regions. 

Open with DEXTER  
In the text 
Fig. 4. Distributions for the number of bright PBHs above the NuSTAR point source detection limit for different levels of turbulence in the CMZ (0–5 km s^{−1}) and for different DM profiles (NFW, Einasto). Each distribution is derived from 1000 independent Monte Carlo realizations. The orange histograms represent the model with only a bulge and DM halo. The green histograms correspond to the case in which the MB 1D velocity dispersion has been increased by 25% to correct for the missing Galactic disk contribution (see Fig. 3). The uniform gas distribution within resolution element is assumed. 

Open with DEXTER  
In the text 
Fig. 5. Analog of Fig. 4 for the model with a smallscale, powerlaw gas density distribution. 

Open with DEXTER  
In the text 
Fig. 3. Velocity dispersion of the PBH DM. The solid curves show our estimate from the Jeans equation. Three cases are shown: the Galaxy consisting only of (i) NFW DM halo, (ii) NFW DM halo + baryonic bulge, and (iii) Einasto DM halo + baryonic bulge. The triangles and squares denote the DM velocity dispersions taken from Ferrer & Hunter (2013). Neglecting the baryonic disk leads to a ∼25% underestimation of the DM velocity dispersion, i.e. a factor of ∼3 overestimation of the accreting PBH luminosity. Replacing NFW with the Einasto density profile leads to a factor of ∼1.5 increase in the velocity dispersion in the very central parts. The light gray shaded region shows an approximate extent of the CMZ. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.