IPMU20-0076 RIKEN-iTHEMS-Report-20 Constraining Primordial Black Holes with Dwarf Galaxy Heating Philip Lu,1, ∗ Volodymyr Takhistov,1, 2, † Graciela B. Gelmini,1, ‡ Kohei Hayashi,3, 4, § Yoshiyuki Inoue,5, 6, 2, ¶ and Alexander Kusenko1, 2, ∗∗ 1Department of Physics and Astronomy, University of California, Los Angeles Los Angeles, California, 90095-1547, USA 2Kavli Institute for the Physics and Mathematics of the Universe (WPI), UTIAS The University of Tokyo, Kashiwa, Chiba 277-8583, Japan 3Institute for Cosmic Ray Research, The University of Tokyo, Kashiwa 277-8582, Japan 4Astronomical Institute, Tohoku University, Sendai 980-8582, Japan 5Department of Earth and Space Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan 6Interdisciplinary Theoretical & Mathematical Science Program (iTHEMS), RIKEN, 2-1 Hirosawa, Saitama 351-0198, Japan (Dated: April 16, 2021) Black holes formed in the early universe, prior to the formation of stars, can exist as dark matter and also contribute to the black hole merger events observed in gravitational waves. We set a new limit on the abundance of primordial black holes (PBHs) by considering interactions of PBHs with the interstellar medium, which result in the heating of gas. We examine generic heating mechanisms, including emission from the accretion disk, dynamical friction, and disk outflows. Using the data from the Leo T dwarf galaxy, we set a new cosmology-independent limit on the abundance of PBHs 7 in the mass range O(1)M − 10 M , relevant for the recently detected gravitational wave signals from intermediate-mass BHs. INTRODUCTION has the desired properties, such as a low baryon velocity dispersion, making it a sensitive probe of PBH heating. Primordial black holes (PBHs) can form in the early While constraints on BHs interacting with surrounding Universe through a variety of mechanisms and can ac- stars have been extensively discussed [48–50], gas heat- count for all or part of the dark matter (DM) (e.g. [1– ing has not been considered in detail. Other constraints 26]). PBHs surviving until present can span many orders focused on the X-ray emission, but not the heating of 15 10 of magnitude in mass, from 10 g to well over 10 M , the surrounding gas [51, 52]. ISM heating has been used and they can account for the entirety of the DM in the to constrain particle DM candidates [53–55], which have −16 −10 mass window ∼ 10 − 10 M , where there are no different heating mechanisms with a different velocity de- observational constraints [27–29]. PBHs with sublunar pendence compared to PBHs. masses can play a role in the synthesis of heavy elements, Additionally, emission from PBH accretion can result production of positrons, as well as other astrophysical in other observables and affect cosmological history, in- phenomena [30–33]. PBHs with larger masses can ac- cluding modifications to reionization and recombination. count for some of the gravitational wave events detected As discussed in [56], ionization from X-rays due to PBH by LIGO (e.g. [34–36]) as well as seeds of supermassive accretion will increase the amount of molecular hydrogen 3 black holes [37–39]. The mass window of 10 − 10 M and the resultant cooling will enhance early star forma- is particularly intriguing in connection with signals ob- tion. PBH emission will also result in spectral distor- served by LIGO [39–45]. Very recently, LIGO announced tions of the cosmic microwave background (CMB). Fur- first detection of an intermediate-mass 150M BH [46]. thermore, heating and ionization of intergalactic medium While a variety of constraints exist for this PBH mass arising from PBH emission will modify 21 cm power spec- range (see Ref. [47] for review), they often rely on multi- trum and will be probed by upcoming experiments [57]. arXiv:2007.02213v2 [astro-ph.CO] 15 Apr 2021 ple assumptions and are subject to significant uncertain- While we do not discuss these effects in detail within our ties. work, they provide complementary probes of PBHs. In this work, we set new cosmology-independent con- straints on PBH abundance based on the lack of gas heat- ing from PBH interactions with the interstellar medium (ISM). We consider several generic heating mechanisms, including dynamical friction, accretion disk emission as BLACK HOLES IN INTERSTELLAR MEDIUM well as mass outflows/winds from accretion disk. We then apply our analysis to dwarf DM-rich galaxies, fo- cusing on Leo T. Leo T is a transitional object between The accretion of gas onto freely floating BHs has been a dwarf irregular galaxy and a dwarf spheroidal galaxy analyzed in [58] and applied to PBHs in [52]. Bondi- that has been well studied and modeled theoretically. It Hoyle-Lyttleton accretion results in the mass accretion 2 rate1 [23, 60, 61] GAS HEATING 2 2 ˙ 2 4πG M nµmp For a gas system in thermal equilibrium, the total M = 4πrBvρ˜ = , (1) v˜3 amount of heating by PBHs of mass M is 2 where M is the PBH mass, rB = GM/v˜ is the Bondi H (M) = N (M)H(M) (4) radius, µ is the mean molecular weight, n is the ISM gas tot PBH Z nmax Z vmax 2 dfn dfv number density, mp is the proton mass andv ˜ ≡ (v + = dndv H(M, n, v) , 2 1/2 dn dv cs) . Here, v is the PBH speed relative to the gas nmin vmin and cs is the temperature-dependent sound speed in gas, where df /dn is the gas density distribution, df /dv is which we take to be cs ∼ 10 km/s [52]. n v The accretion rate can be related to the bolometric the PBH relative speed distribution and H(M, n, v) is emission luminosity as L = (M˙ )M˙ , where (M˙ ) is the the amount of heat deposited into the system from a sin- radiative efficiency which scales with accretion rate. The gle PBH. Here, H represents the cumulative contribution Eddington accretion rate, assuming a characteristic ra- from all heating processes. For photon emission and out- flows we perform an additional integration to treat the diative efficiency of 0 = 0.1, is defined in terms of 2 absorption efficiency. For gas of approximately constant the Eddington luminosity M˙ Edd = LEdd/0c . A con- venient parameter for characterizing the accretion flow is density, one can replace dfn/dn by a delta function. First, we consider gas heating due to photon emission m˙ = M/˙ M˙ Edd. With a sufficient angular momentum, the infalling gas from accretion. Emission in the X-ray band generally can form an accretion disk around the BH. The angular constitutes the dominant contribution and it becomes momentum necessary for a disk formation can be sup- more efficient at high mass accretion rates. plied by perturbations in the density or the velocity of Photon emission from accretion depends on the accre- the accreting gas. For a Schwarzschild BH, the inner tion flow. To characterize the accretion flow, we fol- radius of the disk is taken to be the innermost stable cir- low the scheme outlined in [62] and assume that the accretion flow results in a (geometrically) thin disk for cular orbit (ISCO) of a test particle rISCO = 3rs, where 2 m˙ > M˙ = 0.07α. The thin α-disk is the so-called stan- rs = 2GM/c is the Schwarzschild radius. Following the arguments of [52, 58], we have confirmed that an accre- dard disk [63], where α ∼ 0.1 is a phenomenological tion disk always forms for our parameters of interest. parameter describing viscosity. A thin disk is optically thick and efficiently emits blackbody radiation. Thin disk If PBHs constitute a fraction fPBH of the DM, the total number of PBHs of mass M within a volume V is emission allows for a fully analytic description, and we employ the scaling characterization of [64]. ρDMV For accretion rates with efficiency below the thin disk N (M) = f , (2) PBH PBH M regime, accretion is described by the advection domi- nated accretion flow (ADAF) [62, 65]. Here, the heat where ρDM is the DM density, assumed to be approx- generated by viscosity during accretion is not efficiently imately constant. We assume a monochromatic PBH radiated out, and much of the energy is advected via mat- mass function for definiteness and for presenting our re- ter heat capture into the BH event horizon along with sults in the form of a differential exclusion plot. The ve- the gas inflow. In contrast with the thin disk, the ADAF locity of PBHs contributing to the DM can be described “disk” is geometrically thick and optically thin. by a Maxwell-Boltzmann distribution An ADAF disk results in a complicated multi- component emission spectrum. We consider three com- r 2 2 2 v v ponents of the ADAF spectrum, arising from electron fv(v) = 3 exp − 2 , (3) π σv 2σv cooling: synchrotron radiation, inverse Compton scat- tering and bremsstrahlung. To describe the ADAF spec- where σv is the velocity dispersion in a given system. trum, we employ approximate analytic expressions ob- A distribution in gas number density fn(n) can also be tained in [66] in combination with the updated values for introduced, as in [52, 58]. the phenomenological input parameters consistent with recent numerical simulations and observations [62]. We take the ratio of direct viscous heating to electrons and ions δ = 0.3, and the ratio of gas pressure to total pres- 1 Recent 3D hydrodynamical simulations show that the accretion sure β = 10/11. rate at high Mach number would be limited to ∼ 10–20% of the We do not consider a slim disk or other solutions for canonical Bondi-Hoyle-Lyttleton accretion rate of Eq.