evaporating on the Floor

Roberta Calabrese,1 Damiano F.G. Fiorillo,1 Gennaro Miele,1, 2 Stefano Morisi,1 and Antonio Palazzo3, 4, ∗ 1Universit`adegli Studi di Napoli “Federico II”, Complesso Univ. Monte S. Angelo, I-80126 Napoli, Italy, and INFN - Sezione di Napoli, Complesso Univ. Monte S. Angelo, I-80126 Napoli, Italy 2Scuola Superiore Meridionale, Universit`adi Napoli Federico II, Largo San Marcellino 10, 80138 Napoli, Italy 3Dipartimento Interateneo di Fisica “Michelangelo Merlin”, Universit`adegli Studi di Bari, Via Amendola 173, 70126 Bari, Italy 4Istituto Nazionale di Fisica Nucleare, Sezione di Bari, Via Orabona 4, 70126 Bari, Italy Primordial black holes (PBHs) hypothetically generated in the first instants of life of the Universe are potential dark matter (DM) candidates. Focusing on PBHs masses in the range [5 × 1014 − 5 × 1015]g, we point out that the emitted by PBHs evaporation can interact through the coherent elastic neutrino nucleus scattering (CEνNS) producing an observable signal in multi-ton DM direct detection experiments. We show that with the high exposures envisaged for the next- generation facilities, it will be possible to set bounds on the fraction of DM composed by PBHs improving the existing neutrino limits obtained with Super-Kamiokande. We also quantify to what extent a signal originating from a small fraction of DM in the form of PBHs would modify the so-called “neutrino floor”, the well-known barrier towards detection of weakly interacting massive particles (WIMPs) as the dominant DM component.

Introduction. The identity of dark matter (DM) is been successfully observed by COHERENT [45], where a one of the most puzzling mysteries in contemporary as- few kilograms detector was exposed to an intense neu- troparticle and cosmology. In spite of enormous trino flux of artificial origin. The very same process efforts, no uncontroversial non-gravitational signal of DM involving neutrinos of natural origin, such as the so- has emerged so far. In this context, the recent detec- lar, DSNB and atmospheric ones, constitutes an irre- tion of gravitational waves from binary black hole merg- ducible background [46–52] (forming the so-called “neu- ers by LIGO/Virgo [1,2] has strongly revamped the at- trino floor” [52]) towards detection and identification of tention [3–7] towards the hypothesis that DM may be WIMPs as DM candidates in next-generation direct de- composed of primordial black holes (PBHs). As first tection experiments. Here we show that neutrinos from recognized in the seventies, these objects can be gen- PBHs with masses in the range [5 1014, 5 1015]g, which erated in the early Universe from the collapse of large emit neutrinos with peak energy× between× 10 MeV and overdensities [8–14], and may constitute a fraction of the 100 MeV, may emerge as a signal on top of such a famil- observed amount of DM [10, 13] (see [15–19] for recent iar background. As a consequence, it is possible to set reviews). PBHs emit Hawking radiation [20], and for prospective bounds on the PBHs fraction fPBH of DM 14 large enough masses (MPBH & 5 10 g), have a lifetime in this mass range. As an interesting byproduct of our longer than the age of the Universe.× The evaporation pro- study, we show how the neutrino floor gets modified by cess can give rise to observable signals. In fact, bounds the presence of a hypothetical signal from PBHs. (present or prospective) on PBHs have been obtained Neutrinos emitted by PBHs. A black hole is be- from X-rays [21, 22], γ-rays [21, 23–27], 511 keV γ-ray lieved to quantum evaporate [20], emitting radiation akin line from galactic center [28–30], galactic e± [31], cos- to a hot body. For a neutral non-rotating (Schwarzschild) mic microwave background (CMB) [32–34], radio signals black hole, the Hawking temperature is given by [20, 53– from inverse Compton scattering on CMB photons [35] 55] and synchrotron radiation [36], and heating of the in- 3  16  terstellar medium [37, 38]. The possibility to constrain ~c 10 g kBTPBH = 1.06 MeV , (1) arXiv:2106.02492v1 [hep-ph] 4 Jun 2021 PBHs using the emitted neutrinos was discussed long ago 8πGN MPBH ' MPBH in [39–41]. More recently, limits on PBHs have been ob- tained in [29] exploiting the null searches of the diffuse where four fundamental constants appear, kB (Boltz- supernova neutrino background (DSNB) performed by mann), GN (gravitational), ~ (Planck), and c (speed of Super-Kamiokande [42]. Also, prospective bounds from light). The differential flux per unit time of emitted par- the experiment JUNO have been discussed in [43]. ticles depends on their spin. For spin 1/2 particles with In this Letter, we entertain a novel possibility, never mass negligible with respect to TPBH like neutrinos, it is addressed before in the literature, proposing to detect given by the emitted neutrinos from PBHs by coherent elastic neu- 2 trino nucleus scattering (CEνNS). It is only recently that d Nν 1 Γν (Eν ,MPBH) the CEνNS process, predicted long time ago [44], has = , (2) dEν dt 2π exp [Eν /kBTPBH] + 1

where Eν is the neutrino energy and Γν is a graybody factor [53–55] encoding the imprint of the space-time ge- ∗ Corresponding author; email: [email protected] ometry intervening between the event horizon and the 2

asymptotic observer. In our analysis, we employ the pub- ]

1 pp licly available BlackHawk code [56] for calculating the − 11 Background 10 15 2 M = 3 10 g, f = 5.2 10− energy spectra of the emitted neutrinos. BlackHawk pro- PBH · PBH · 13 15 3 N M = 2 10 g, f = 4.0 10− MeV 15 PBH PBH

1 · · vides the primary spectra for all fundamental Standard 8 O 14 5

− 10 MPBH = 7 10 g, fPBH = 7.0 10− Model particles using the Hawking evaporation spectrum s · · 2 17 8

− F B in Eq. (2). In addition, the code generates the spec- 5 tra of secondary neutrinos deriving from hadronization 10 of strongly-interacting constituents and decay of unsta- 2 hep ble particles. We sum up both kind of spectra in our 10 DSNB analysis. In [57], it has been shown that Dirac neutrinos, 1 having twice as many degrees of freedom as the Majo- 10− Atm rana ones, would affect the PBHs evaporation making it faster (see also [58, 59]). Although the additional Dirac Neutrino Flux [cm 4 10− 1 0 1 2 3 degrees of freedom are sterile for the electroweak inter- 10− 10 10 10 10 actions and not detectable in the CEνNS process, their Eν [MeV] existence can alterate indirectly the emission rate of the FIG. 1. Neutrino fluxes from PBHs. The black dashed active degrees of freedom. For the PBHs masses con- contours represent the background fluxes originating from so- sidered in our work, this effect can be quantified around lar, DSNB and atmospheric neutrinos. The colored solid lines 10% [57]. For definiteness, we assume that neutrinos ∼ correspond to neutrinos from PBHs evaporation for three rep- have Majorana nature. We ignore neutrino oscillations resentative values of their mass and fraction of total DM con- being irrelevant for the flavor-blind CEνNS process. We tent. These benchmark values lie on the 90% C.L. exclusion take into account both the contributions coming from the curve obtainable from a liquid experiment with 200 t galactic and extragalactic PBHs. The galactic differen- yr exposure (corresponding to the black stars in Fig.3). tial neutrino flux is given by

MW Z 2 Z equality to t , the minimum between the PBH lifetime dφν dΩ d N fPBH ρMW [r(l, ψ)] max = dl , (3) and the age of the Universe. The overall neutrino flux dEν 4π dEν dt MPBH from the sum of galactic and extragalactic contributions where ρMW(r) is the DM density of the Milky Way is plotted in Fig.1 for three benchmark values of MPBH (MW), r denotes the galactocentric distance and fPBH. As it will be discussed in the next section, these values are excludable at 90% C.L. from a xenon q 2 2 experiment with 200 t yr exposure, assuming a measure- r(l, ψ) = r 2lr cos ψ + l , (4) − ment compatible with the ordinary background. In the same plot, we represent the background which is formed with r being the distance of the Sun from the galactic from solar [63], DSNB [64] and low-energy atmospheric center, l the line-of-sight distance to the PBH, ψ the angle neutrinos [65, 66]. As expected from Eq. (1), smaller between these two directions, and Ω the solid angle under PBHs masses correspond to harder spectra of the emit- consideration. For definiteness, we employ a Navarro- ted neutrinos with peak located at 4.2 T [67]. We Frenk-White (NFW) profile [60] PBH see that PBHs neutrinos can be visible∼ above the abrupt  −1  2 fall-off of the solar hep neutrinos, where the dominant r 1 + r /rs ρMW(r) = ρ , (5) contribution to the background is provided by DSNB and r 1 + r/rs atmospheric neutrinos. Coherent scattering of neutrinos. Coherent elas- −3 where we take ρ = 0.4 GeV cm for the DM density in tic scattering of a neutrino ν (or antineutrinoν ¯) on a the solar neighborhood, r = 8.5 kpc, and rs = 20 kpc nucleus can occur if qR 1, where q = ~q is the three- for the scale radius. We stress that the value ρ = momentumN transfer and Ris the nuclear radius| | [44]. The −3 0.4 GeV cm lies at the lower end of its allowed range differential cross section can be expressed as [44] according to the latest estimates (see for example [61]). Therefore, the bounds we are going to derive will be dσ G2 m  m E  νN (E ,E ) = F N Q2 1 N R F 2(q) , conservative in this respect. For the extragalactic com- ν R w 2 dER 4π − 2Eν ponent, the differential flux integrated over the full sky (7) is [16] where GF is the Fermi constant, mN is the nucleus mass,  2  Qw = N Z(1 4 sin θW ) is the weak vector nuclear dφEG Z tmax f ρ d2N − − ν PBH DM ν charge, Z and N are the number of protons and neu- = dt [1 + z(t)] , 2 ˜ E˜ =[1+z(t)]E dEν tmin MPBH dEν dt ν ν trons, sin θW = 0.231 [68] is the Weinberg angle, ER (6) is the nucleus recoil energy, and Eν is the neutrino en- ρDM being the average DM density of the Universe at the ergy. The recoil energy can assume the maximum value max 2 present epoch, as determined by Planck [62]. The time ER = 2Eν /(mN + 2Eν ). For the nuclear form fac- integral runs from tmin set to the era of matter-radiation tor F (q), encoding the loss of coherence for qR > 1, we 3

2 10 defined as Background 15 2 2  2 T −1  M = 3 10 g, f = 5.2 10− PBH · PBH · χ (θ) = min χ (θ, α) + (1 α) Σα (1 α) , (9) ] 15 3 α − − 1 1 M = 2 10 g, f = 4.0 10− 10 PBH · PBH · − 14 5 MPBH = 7 10 g, fPBH = 7.0 10− T T · · where θ = [MPBH, fPBH], and α = [α1, α1, α3] rep- resent respectively the vector of the model parameters 100 and that of the nuisance parameters associated to the normalization of the three backgrounds components (so- lar hep, DSNB, atmospheric) with respect to their best 1 theoretical estimates. The uncertainties on the back-

Event Rate10 [keV − grounds are encoded in the covariance matrix Σα = thr diag(σ2 , σ2 , σ2 ), which we take diagonal because the ER = 5 keV 200 t yr α1 α2 α3 10 2 three fluxes have completely independent origin. For the − 101 uncertainties, we have assumed 30%, 50% and 20%, re- ER [keV] spectively for solar hep [63], DSNB [64] and atmospheric neutrinos [65, 66]. The first term in Eq. (9) is defined as FIG. 2. Differential neutrino events rate. The black dashed contour represents the total background rate (solar + 2 L0 DSNB + atmospheric). The colored solid lines correspond χ (θ, α) = 2 ln , (10) L to neutrinos emitted from PBHs evaporation for three repre- − 1 sentative values of their mass and fraction of total DM con- with likelihoods given by tent. These benchmark values are excludable at the 90% C.L. by a liquid xenon experiment with 200 t yr exposure (corre- Y i i i sponding to the black stars in Fig.3). The grey vertical line L0 = P (x = N Bck; λ = NPBH(θ) + NBck(α)) , (11) indicates the threshold recoil energy used in the analysis. i and Y i i employ the Helm parametrization [69] using the recipes L1 = P (x = N Bck; λ = N Bck) , (12) provided in [70]. The differential rate of events is given i by where P (x, λ) is the Poisson distribution, i is the bin in- i dex, N Bck is the nominal number of background events dRνN i = NT (ER) expected in the i-th bin, whereas NPBH is the number dE dt i R (8) of PBHs events, and N i = P α N is the float- Z dσ dφ Bck j j Bck,j dE νN ν Θ(Emax E ) , ing background. The second contribution in Eq. (9) is a ν dE dE R R × R ν − penalty term gauging the statistical weight of the devi- ation of the background fluxes from their central values. where NT is the number of target nuclei in the detector, We have neglected other possible backgrounds, in par- (ER) is the detector efficiency (assumed to be equal to ticular those arising from electron recoils of solar pp and one), dφν /dEν is the differential neutrino flux, and Θ is 7Be neutrinos, since as shown in [71], these can be ef- the Heaviside step function. In Fig.2, we show the dif- fectively suppressed through a statistical discrimination ferential rate of events as a function of the recoil energy of the photon and ionization signals. In our statistical for the background (dashed line), and for three bench- analysis, we employ ten bins in the recoil energy window mark values of the PBHs parameters (the same used in thr [5 50] keV. The choice of the threshold ER = 5 keV Fig.1). The plots refers to an exposure of 200 t yr. One is dictated− by the position of the sharp cutoff of the hep can observe that the shape of the spectrum induced by neutrino event rate, which in this way has a marginal im- PBHs neutrinos appreciably changes with the value of pact in our results. In view of the most recent theoretical the mass of the PBH. In particular, for smaller values evaluations of the DSNB [72–75], which point towards of M , which correspond to more energetic neutrino PBH somewhat larger fluxes of the νe component with respect fluxes, the event spectrum is similar to the atmospheric to previous findings, as a check, we have increased by a background. In contrast, for larger values of MPBH with factor of two the nominal value [64] of the DSNB flux for flux peaked at lower energies, the PBHs event spectrum both ν and ν and for all three flavors. We find that the ex- is quite different with respect to the background. For this clusion limits are basically unchanged because in the re- reason, in the statistical treatment presented in the next gion of interest the DSNB has a subleading role, with the section we employ a binned likelihood analysis, so as to dominant contribution arising from the atmospheric neu- exploit the information contained in the spectral shape. trinos. Figure3 displays the 90% C.L. exclusion limits PBHs at next-generation detectors. In order to obtainable from a xenon detector with the three bench- derive prospective upper limits on the fraction of DM mark exposures reported in the legend. The 20 t yr ex- composed of PBHs, we implement the χ2 test statistic posure should be attainable in LZ and XENONnT, while 4

0 10 the solar 8B neutrino background, which is dominant for WIMP masses around mχ 6 GeV. The next-generation 1 10− facilities with very high exposures∼ will unavoidably en- counter the neutrino floor also for larger WIMP masses, 2 10− where the background from DSNB and atmospheric neu- trinos is relevant. Going below such a floor will re- 3

PBH 10− quire exploiting timing structure of the signal [80], com- f bining different targets [81], or using directional infor- 4 mation [82–84]. Interestingly, the neutrino background 10− 20 t yr (LZ, XENONnT) 200 t yr (DARWIN) can be affected by exotic neutrino interactions [85–88], 5 with consequent modification of the standard neutrino 10− 2000 t yr Super K floor [89–94], which can be influenced also by neutrinos 6 originating from decay of massive particles [95, 96]. 10− 15 16 10 10 As illustrated in the previous sections, the neutrinos MPBH [g] emitted by PBHs would lie on top of such a familiar FIG. 3. Upper bounds on PBHs. The exclusion contours, background. Therefore, the existence of even a minute all drawn at the 90% C.L. for 1 d.o.f., refer to a xenon experi- fraction of PBHs in the DM content would modify the ment with three different exposures. Two of them correspond neutrino floor, making it higher. Here we quantify this to the planned experiments (LZ/XENONnT and DARWIN), effect by calculating the floor following the prescription while the third one refers to an ideal setup. The black stars in- of [52]. For definiteness and consistency with the pre- dicate the three benchmark points used for Figures1,2 and4. vious sections, we consider the case of a xenon target The shaded grey region represents the upper bounds obtained nucleus and adopt the same NFW profile as in Eq. (5). in [29] from the DSNB searches of Super-Kamiokande [42]. Following [70], we employ a galactic Maxwell-Boltzmann velocity distribution with most probable speed v0 = 220 km/s, truncated at the escape velocity v = 544 km/s the higher value of 200 t yr refers to the more ambitious esc and boosted into the laboratory frame with v = 232 project DARWIN. The third contour corresponds to 2000 lab km/s. For a fixed value of the WIMP mass we calcu- t yr, which probably may be considered as an ultimate late the cross-section σ that can be excluded at the goal for liquid xenon detectors. Such high exposures seem χn 90% confidence level (corresponding to 2.3 DM events) to be less extreme for Argon based detectors [76]. In the selecting the exposure which leads to 1 CEνNS count. same plot, for comparison, we report as a shaded grey re- Then, varying over the energy threshold Ethr, we keep gion, the upper bounds obtained in [29] from the DSNB R the value which minimizes the cross-section. By repeat- searches of Super-Kamiokande [42].1 We see that an im- ing such a procedure for a dense grid of WIMP masses in provement of one order of magnitude is attainable with the range [1, 1000] GeV, we construct the contours shown the most high exposures. Further amelioration will be in Fig.4. The black upper border of the yellow region possible with a better theoretical knowledge of the atmo- corresponds to the well-known ordinary floor. The up- spheric neutrino fluxes. per edge of the colored regions corresponds to the case Impact of PBHs on the neutrino floor. Neu- of PBHs with parameters reported in the legend, where trinos originating from the Sun, diffuse supernovae and for the sake of clarity, we use the same values and color Earth’s atmosphere constitute an irreducible background convention adopted in Figs.1 and2. in DM direct searches [46–51]. This background gives rise to the so-called “neutrino floor”, an ultimate limit to the Here, an important remark is in order. For suffi- 21 −12 discovery potential in the plane spanned by the WIMP ciently high PBH masses MPBH 10 g 10 M & ∼ mass mχ and the spin-independent WIMP-nucleon cross (M being the solar mass), the existence of WIMPs as section σχn. Of course, such a limitation holds for any a dominant or subdominant component of DM, is in- kind of interaction, and can be generalized to other types compatible even with a small fraction of DM made of of WIMP-nucleus effective field theory operators [77, 78], PBHs [97–104]. This occurs because WIMPs are ac- but is customary to adopt the case of spin-independent creted by the PBH potential well, forming spiked ultra- interaction as a benchmark. As a matter of fact, the compact mini-halos, whose existence can be excluded by running experiment XENON1T [79] is already surfacing the non-observation of the expected products of WIMPs self-annihilation. However, for the low PBHs masses con- sidered in our work, falling in the so-called evaporation 16 −17 range (MPBH . 10 g 10 M ), this phenomenon is 1 In [29], which considers both spinning and non-spinning PBHs, negligible. In addition,∼ one should note that DM may 15 the upper limits are reported for masses above 10 g because be asymmetric [105–107], and this may occur also for rotating black holes evaporate faster than the non-rotating ones WIMPs [108], in which case self-annihilation is reduced and can contribute to the present DM content for MPBH & 7 × 1015g. In our work, focused on non-rotating PBHs, we show the (albeit not necessarily absent [109–112]). Therefore, we 15 limits down to MPBH = 5 × 10 g, as it is customary in the deem it interesting to consider how the neutrino floor literature. would get modified by the presence of a minute DM frac- 5

] 40 2 10− for the next-generation facilities, it will be possible to Standard Neutrino Floor set bounds on the fraction of DM composed of PBHs, 15 2 M = 3 10 g, f = 5.2 10− 42 PBH · PBH · improving the existing neutrino limits obtained with 10− 15 3 M = 2 10 g, f = 4.0 10− PBH · PBH · 14 5 Super-Kamiokande. In addition, we have quantified how M = 7 10 g, f = 7.0 10− PBH · PBH · 44 much a signal originating from PBHs would heighten the 10− so-called “neutrino floor”, the ultimate barrier towards detection of WIMPs as the dominant DM component. 46 10− For definiteness we have focused our study on liquid xenon detectors such as DARWIN [113], LZ [114], and 48 XENONnT [115], but other targets such as liquid argon 10− employed in DarkSide-20k [76] and ARGO [76], or archeological lead in RES-NOVA [116], should offer a 50 WIMP-nucleon cross-section10 [cm − 0 1 2 3 similar opportunity provided that very high exposures 10 10 10 10 are reached (see [117] for an extensive review of the WIMP mass [GeV] experimental program of the direct detection facilities). Finally, we underline that, in the context of PBHs FIG. 4. Impact of PBHs on the Neutrino floor. The searches, the direct DM experiments would rather oper- black contour delimiting the yellow region represents the or- dinary neutrino floor, while the upper border of the colored ate as indirect DM observatories. From this perspective, bands correspond to the modifications induced by neutrinos our study lends further support to the emerging role of from PBHs with masses and DM fractions in the legend. such underground facilities as multi-purpose low-energy These benchmark values lie on the 90% C.L. exclusion curve neutrino telescopes complementary to their high-energy obtainable from a liquid xenon experiment with 200 t yr ex- “ordinary” counterparts, IceCube and KM3NeT. posure (corresponding to the black stars in Fig.3). ACKNOWLEDGMENTS tion of PBHs having mass in the range considered in this work. We thank B. Dasgupta and D. Montanino for help- Conclusions. We have explored a new avenue in ful discussions. We acknowledge partial support by constraining PBHs as Dark Matter with the neutrinos the research grant number 2017W4HA7S “NAT-NET: emitted in the Hawking radiation. Specifically, we have Neutrino and Astroparticle Theory Network” under the pointed out that neutrinos from PBHs can interact program PRIN 2017 funded by the Italian Ministero via the coherent elastic neutrino nucleus scattering dell’Istruzione, dell’Universit`ae della Ricerca (MIUR) (CEνNS) in multi-ton DM direct detection experiments. and by the research project TAsP funded by the Instituto We have shown that with the high exposures envisaged Nazionale di Fisica Nucleare (INFN).

[1] LIGO Scientific, Virgo collaboration, GW190521: [6] S. Blinnikov, A. Dolgov, N.K. Porayko and K. Postnov, A Binary Black Hole Merger with a Total Mass of Solving puzzles of GW150914 by primordial black 150M , Phys. Rev. Lett. 125 (2020) 101102 holes, JCAP 11 (2016) 036[ 1611.00541]. [2009.01075]. [7] T. Nakama, J. Silk and M. Kamionkowski, Stochastic [2] LIGO Scientific, Virgo collaboration, Binary Black gravitational waves associated with the formation of Hole Population Properties Inferred from the First and primordial black holes, Phys. Rev. D 95 (2017) 043511 Second Observing Runs of Advanced LIGO and [1612.06264]. Advanced Virgo, Astrophys. J. Lett. 882 (2019) L24 [8] Y.B..N. Zel’dovich, I. D., The Hypothesis of Cores [1811.12940]. Retarded during Expansion and the Hot Cosmological [3] S. Bird, I. Cholis, J.B. Mu˜noz,Y. Ali-Ha¨ımoud, Model, Soviet Astron. AJ (Engl. Transl. ), 10 (1967) M. Kamionkowski, E.D. Kovetz et al., Did LIGO 602. detect dark matter?, Phys. Rev. Lett. 116 (2016) [9] S. Hawking, Gravitationally collapsed objects of very 201301[ 1603.00464]. low mass, Mon. Not. Roy. Astron. Soc. 152 (1971) 75. [4] S. Clesse and J. Garc´ıa-Bellido, The clustering of [10] B.J. Carr and S.W. Hawking, Black holes in the early massive Primordial Black Holes as Dark Matter: Universe, Mon. Not. Roy. Astron. Soc. 168 (1974) 399. measuring their mass distribution with Advanced [11] S.W. Hawking, Black hole explosions, Nature 248 LIGO, Phys. Dark Univ. 15 (2017) 142[ 1603.05234]. (1974) 30. [5] M. Sasaki, T. Suyama, T. Tanaka and S. Yokoyama, [12] B.J. Carr, The Primordial black hole mass spectrum, Primordial Black Hole Scenario for the Astrophys. J. 201 (1975) 1. Gravitational-Wave Event GW150914, Phys. Rev. Lett. [13] G.F. Chapline, Cosmological effects of primordial black 117 (2016) 061101[ 1603.08338]. holes, Nature 253 (1975) 251. 6

[14] M. Khlopov, B.A. Malomed and I.B. Zeldovich, [1610.10051]. Gravitational instability of scalar fields and formation [33] S. Clark, B. Dutta, Y. Gao, L.E. Strigari and of primordial black holes, Mon. Not. Roy. Astron. Soc. S. Watson, Planck Constraint on Relic Primordial 215 (1985) 575. Black Holes, Phys. Rev. D 95 (2017) 083006 [15] M. Sasaki, T. Suyama, T. Tanaka and S. Yokoyama, [1612.07738]. Primordial black holes—perspectives in gravitational [34] H. Poulter, Y. Ali-Ha¨ımoud,J. Hamann, M. White wave astronomy, Class. Quant. Grav. 35 (2018) and A.G. Williams, CMB constraints on ultra-light 063001[ 1801.05235]. primordial black holes with extended mass [16] B. Carr, K. Kohri, Y. Sendouda and J. Yokoyama, distributions, 1907.06485. Constraints on Primordial Black Holes, 2002.12778. [35] B. Dutta, A. Kar and L.E. Strigari, Constraints on [17] A.M. Green and B.J. Kavanagh, Primordial Black MeV dark matter and primordial black holes: Inverse Holes as a dark matter candidate, J. Phys. G 48 Compton signals at the SKA, JCAP 03 (2021) 011 (2021) 4[ 2007.10722]. [2010.05977]. [18] B. Carr and F. Kuhnel, Primordial Black Holes as [36] M.H. Chan and C.M. Lee, Constraining Primordial Dark Matter: Recent Developments, Ann. Rev. Nucl. Black Hole Fraction at the Galactic Centre using radio Part. Sci. 70 (2020) 355[ 2006.02838]. observational data, Mon. Not. Roy. Astron. Soc. 497 [19] P. Villanueva-Domingo, O. Mena and (2020) 1212[ 2007.05677]. S. Palomares-Ruiz, A brief review on primordial black [37] R. Laha, P. Lu and V. Takhistov, Gas Heating from holes as dark matter, 2103.12087. Spinning and Non-Spinning Evaporating Primordial [20] S.W. Hawking, Particle Creation by Black Holes, Black Holes, 2009.11837. Commun. Math. Phys. 43 (1975) 199. [38] H. Kim, A constraint on light primordial black holes [21] G. Ballesteros, J. Coronado-Bl´azquezand D. Gaggero, from the interstellar medium temperature, 2007.07739. X-ray and gamma-ray limits on the primordial black [39] F. Halzen, B. Keszthelyi and E. Zas, Neutrinos from hole abundance from Hawking radiation, Phys. Lett. B primordial black holes, Phys. Rev. D 52 (1995) 3239 808 (2020) 135624[ 1906.10113]. [hep-ph/9502268]. [22] J. Iguaz, P.D. Serpico and T. Siegert, Isotropic X-ray [40] E.V. Bugaev and K.V. Konishchev, Constraints on bound on Primordial Black Hole Dark Matter, diffuse neutrino background from primordial black 2104.03145. holes, Phys. Rev. D 65 (2002) 123005 [23] B.J. Carr, K. Kohri, Y. Sendouda and J. Yokoyama, [astro-ph/0005295]. Constraints on primordial black holes from the [41] E.V. Bugaev and K.V. Konishchev, Cosmological Galactic gamma-ray background, Phys. Rev. D 94 constraints from evaporations of primordial black holes, (2016) 044029[ 1604.05349]. Phys. Rev. D 66 (2002) 084004[ astro-ph/0206082]. [24] A. Arbey, J. Auffinger and J. Silk, Constraining [42] Super-Kamiokande collaboration, Supernova Relic primordial black hole masses with the isotropic gamma Neutrino Search at Super-Kamiokande, Phys. Rev. D ray background, Phys. Rev. D 101 (2020) 023010 85 (2012) 052007[ 1111.5031]. [1906.04750]. [43] S. Wang, D.-M. Xia, X. Zhang, S. Zhou and Z. Chang, [25] R. Laha, J.B. Mu˜nozand T.R. Slatyer, INTEGRAL Constraining primordial black holes as dark matter at constraints on primordial black holes and particle dark JUNO, Phys. Rev. D 103 (2021) 043010[ 2010.16053]. matter, Phys. Rev. D 101 (2020) 123514[ 2004.00627]. [44] D.Z. Freedman, Coherent Neutrino Nucleus Scattering [26] A. Coogan, L. Morrison and S. Profumo, Direct as a Probe of the Weak Neutral Current, Phys. Rev. D Detection of Hawking Radiation from Asteroid-Mass 9 (1974) 1389. Primordial Black Holes, Phys. Rev. Lett. 126 (2021) [45] COHERENT collaboration, Observation of Coherent 171101[ 2010.04797]. Elastic Neutrino-Nucleus Scattering, Science 357 [27] A. Ray, R. Laha, J.B. Mu˜nozand R. Caputo, Closing (2017) 1123[ 1708.01294]. the gap: Near future MeV telescopes can discover [46] B. Cabrera, L.M. Krauss and F. Wilczek, Bolometric asteroid-mass primordial black hole dark matter, Detection of Neutrinos, Phys. Rev. Lett. 55 (1985) 25. 2102.06714. [47] A.K. Drukier, K. Freese and D.N. Spergel, Detecting [28] R. Laha, Primordial Black Holes as a Dark Matter Candidates, Phys. Rev. D 33 (1986) Candidate Are Severely Constrained by the Galactic 3495. Center 511 keV γ -Ray Line, Phys. Rev. Lett. 123 [48] J. Monroe and P. Fisher, Neutrino Backgrounds to (2019) 251101[ 1906.09994]. Dark Matter Searches, Phys. Rev. D 76 (2007) 033007 [29] B. Dasgupta, R. Laha and A. Ray, Neutrino and [0706.3019]. positron constraints on spinning primordial black hole [49] J.D. Vergados and H. Ejiri, Can Solar Neutrinos be a dark matter, Phys. Rev. Lett. 125 (2020) 101101 Serious Background in Direct Dark Matter Searches?, [1912.01014]. Nucl. Phys. B 804 (2008) 144[ 0805.2583]. [30] C. Keith and D. Hooper, The 511 keV Excess and [50] L.E. Strigari, Neutrino Coherent Scattering Rates at Primordial Black Holes in our Solar System, Direct Dark Matter Detectors, New J. Phys. 11 (2009) 2103.08611. 105011[ 0903.3630]. [31] M. Boudaud and M. Cirelli, Voyager 1 e± Further [51] A. Gutlein et al., Solar and atmospheric neutrinos: Constrain Primordial Black Holes as Dark Matter, Background sources for the direct dark matter search, Phys. Rev. Lett. 122 (2019) 041104[ 1807.03075]. Astropart. Phys. 34 (2010) 90[ 1003.5530]. [32] V. Poulin, J. Lesgourgues and P.D. Serpico, [52] J. Billard, L. Strigari and E. Figueroa-Feliciano, Cosmological constraints on exotic injection of Implication of neutrino backgrounds on the reach of electromagnetic energy, JCAP 03 (2017) 043 next generation dark matter direct detection 7

experiments, Phys. Rev. D 89 (2014) 023524 [72] A. Priya and C. Lunardini, Diffuse neutrinos from [1307.5458]. luminous and dark supernovae: prospects for upcoming [53] D.N. Page, Particle Emission Rates from a Black Hole: detectors at the O(10) kt scale, JCAP 11 (2017) 031 Massless Particles from an Uncharged, Nonrotating [1705.02122]. Hole, Phys. Rev. D 13 (1976) 198. [73] S. Horiuchi, K. Sumiyoshi, K. Nakamura, T. Fischer, [54] D.N. Page, Particle Emission Rates from a Black Hole. A. Summa, T. Takiwaki et al., Diffuse supernova 3. Charged Leptons from a Nonrotating Hole, Phys. neutrino background from extensive core-collapse Rev. D 16 (1977) 2402. simulations of 8-100M progenitors, Mon. Not. Roy. [55] J.H. MacGibbon and B.R. Webber, Quark and gluon Astron. Soc. 475 (2018) 1363[ 1709.06567]. jet emission from primordial black holes: The [74] K. Møller, A.M. Suliga, I. Tamborra and P.B. Denton, instantaneous spectra, Phys. Rev. D 41 (1990) 3052. Measuring the supernova unknowns at the [56] A. Arbey and J. Auffinger, BlackHawk: A public code next-generation neutrino telescopes through the diffuse for calculating the Hawking evaporation spectra of any neutrino background, JCAP 05 (2018) 066 black hole distribution, Eur. Phys. J. C 79 (2019) 693 [1804.03157]. [1905.04268]. [75] D. Kresse, T. Ertl and H.-T. Janka, Stellar Collapse [57] C. Lunardini and Y.F. Perez-Gonzalez, Dirac and Diversity and the Diffuse Supernova Neutrino Majorana neutrino signatures of primordial black Background, Astrophys. J. 909 (2021) 169 holes, JCAP 08 (2020) 014[ 1910.07864]. [2010.04728]. [58] D. Hooper, G. Krnjaic and S.D. McDermott, Dark [76] DarkSide-20k collaboration, DarkSide-20k: A 20 Radiation and Superheavy Dark Matter from Black tonne two-phase LAr TPC for direct dark matter Hole Domination, JHEP 08 (2019) 001[ 1905.01301]. detection at LNGS, Eur. Phys. J. Plus 133 (2018) 131 [59] I. Masina, Dark matter and from [1707.08145]. evaporating primordial black holes, Eur. Phys. J. Plus [77] J.B. Dent, B. Dutta, J.L. Newstead and L.E. Strigari, 135 (2020) 552[ 2004.04740]. Effective field theory treatment of the neutrino [60] J.F. Navarro, C.S. Frenk and S.D.M. White, A background in direct dark matter detection experiments, Universal density profile from hierarchical clustering, Phys. Rev. D 93 (2016) 075018[ 1602.05300]. Astrophys. J. 490 (1997) 493[ astro-ph/9611107]. [78] G.B. Gelmini, V. Takhistov and S.J. Witte, Casting a [61] M. Benito, A. Cuoco and F. Iocco, Handling the Wide Signal Net with Future Direct Dark Matter Uncertainties in the Galactic Dark Matter Distribution Detection Experiments, JCAP 07 (2018) 009 for Particle Dark Matter Searches, JCAP 03 (2019) [1804.01638]. 033[ 1901.02460]. [79] XENON collaboration, Search for Coherent Elastic [62] Planck collaboration, Planck 2018 results. VI. Scattering of Solar 8B Neutrinos in the XENON1T Cosmological parameters, Astron. Astrophys. 641 Dark Matter Experiment, Phys. Rev. Lett. 126 (2021) (2020) A6[ 1807.06209]. 091301[ 2012.02846]. [63] N. Vinyoles, A.M. Serenelli, F.L. Villante, S. Basu, [80] J.H. Davis, Dark Matter vs. Neutrinos: The effect of J. Bergstr¨om,M.C. Gonzalez-Garcia et al., A new astrophysical uncertainties and timing information on Generation of Standard Solar Models, Astrophys. J. the neutrino floor, JCAP 03 (2015) 012[ 1412.1475]. 835 (2017) 202[ 1611.09867]. [81] F. Ruppin, J. Billard, E. Figueroa-Feliciano and [64] J.F. Beacom, The Diffuse Supernova Neutrino L. Strigari, Complementarity of dark matter detectors Background, Ann. Rev. Nucl. Part. Sci. 60 (2010) 439 in light of the neutrino background, Phys. Rev. D 90 [1004.3311]. (2014) 083510[ 1408.3581]. [65] G. Battistoni, A. Ferrari, T. Montaruli and P.R. Sala, [82] P. Grothaus, M. Fairbairn and J. Monroe, Directional The atmospheric neutrino flux below 100-MeV: The Dark Matter Detection Beyond the Neutrino Bound, FLUKA results, Astropart. Phys. 23 (2005) 526. Phys. Rev. D 90 (2014) 055018[ 1406.5047]. [66] M. Honda, T. Kajita, K. Kasahara and S. Midorikawa, [83] C.A.J. O’Hare, A.M. Green, J. Billard, Improvement of low energy atmospheric neutrino flux E. Figueroa-Feliciano and L.E. Strigari, Readout calculation using the JAM nuclear interaction model, strategies for directional dark matter detection beyond Phys. Rev. D 83 (2011) 123001[ 1102.2688]. the neutrino background, Phys. Rev. D 92 (2015) [67] K. Kohri and J. Yokoyama, Primordial black holes and 063518[ 1505.08061]. primordial nucleosynthesis. 1. Effects of hadron [84] C.A.J. O’Hare, Can we overcome the neutrino floor at injection from low mass holes, Phys. Rev. D 61 (2000) high masses?, Phys. Rev. D 102 (2020) 063024 023501[ astro-ph/9908160]. [2002.07499]. [68] Particle Data Group collaboration, Review of [85] R. Harnik, J. Kopp and P.A.N. Machado, Exploring nu Particle Physics, Chin. Phys. C 40 (2016) 100001. Signals in Dark Matter Detectors, JCAP 07 (2012) [69] R.H. Helm, Inelastic and Elastic Scattering of 187-Mev 026[ 1202.6073]. Electrons from Selected Even-Even Nuclei, Phys. Rev. [86] D.G. Cerde˜no,M. Fairbairn, T. Jubb, 104 (1956) 1466. P.A.N. Machado, A.C. Vincent and C. Bœhm, Physics [70] J.D. Lewin and P.F. Smith, Review of mathematics, from solar neutrinos in dark matter direct detection numerical factors, and corrections for dark matter experiments, JHEP 05 (2016) 118[ 1604.01025]. experiments based on elastic nuclear recoil, Astropart. [87] D.K. Papoulias, R. Sahu, T.S. Kosmas, V.K.B. Kota Phys. 6 (1996) 87. and B. Nayak, Novel neutrino-floor and dark matter [71] J.L. Newstead, R.F. Lang and L.E. Strigari, searches with deformed shell model calculations, Adv. Atmospheric neutrinos in a next-generation xenon High Energy Phys. 2018 (2018) 6031362[ 1804.11319]. dark matter experiment, 2002.08566. 8

[88] A.M. Suliga and I. Tamborra, Astrophysical constraints [102] B. Carr, F. Kuhnel and L. Visinelli, Black Holes and on nonstandard coherent neutrino-nucleus scattering, WIMPs: All or Nothing or Something Else, Phys. Rev. D 103 (2021) 083002[ 2010.14545]. 2011.01930. [89] J.B. Dent, B. Dutta, J.L. Newstead and L.E. Strigari, [103] K. Kadota and J. Silk, Boosting small-scale structure Dark matter, light mediators, and the neutrino floor, via primordial black holes and implications for Phys. Rev. D 95 (2017) 051701[ 1607.01468]. sub-GeV dark matter annihilation, Phys. Rev. D 103 [90] E. Bertuzzo, F.F. Deppisch, S. Kulkarni, (2021) 043530[ 2012.03698]. Y.F. Perez Gonzalez and R. Zukanovich Funchal, Dark [104] H. Tashiro and K. Kadota, Constraining Mixed Matter and Exotic Neutrino Interactions in Direct Dark-Matter Scenarios of WIMPs and Primordial Detection Searches, JHEP 04 (2017) 073[ 1701.07443]. Black Holes from CMB and 21-cm observations, [91] M.C. Gonzalez-Garcia, M. Maltoni, 2104.09738. Y.F. Perez-Gonzalez and R. Zukanovich Funchal, [105] D.E. Kaplan, M.A. Luty and K.M. Zurek, Asymmetric Neutrino Discovery Limit of Dark Matter Direct Dark Matter, Phys. Rev. D 79 (2009) 115016 Detection Experiments in the Presence of [0901.4117]. Non-Standard Interactions, JHEP 07 (2018) 019 [106] K. Petraki and R.R. Volkas, Review of asymmetric [1803.03650]. dark matter, Int. J. Mod. Phys. A 28 (2013) 1330028 [92] C. Bœhm, D.G. Cerde˜no,P.A.N. Machado, [1305.4939]. A. Olivares-Del Campo, E. Perdomo and E. Reid, How [107] K.M. Zurek, Asymmetric Dark Matter: Theories, high is the neutrino floor?, JCAP 01 (2019) 043 Signatures, and Constraints, Phys. Rept. 537 (2014) [1809.06385]. 91[ 1308.0338]. [93] W. Chao, J.-G. Jiang, X. Wang and X.-Y. Zhang, [108] M.L. Graesser, I.M. Shoemaker and L. Vecchi, Direct Detections of Dark Matter in the Presence of Asymmetric WIMP dark matter, JHEP 10 (2011) 110 Non-standard Neutrino Interactions, JCAP 08 (2019) [1103.2771]. 010[ 1904.11214]. [109] T. Cohen and K.M. Zurek, Leptophilic Dark Matter [94] S. Sadhukhan and M.P. Singh, Neutrino floor in from the Lepton Asymmetry, Phys. Rev. Lett. 104

leptophilic U(1) models: Modification in U(1)Lµ−Lτ , (2010) 101301[ 0909.2035]. Phys. Rev. D 103 (2021) 015015[ 2006.05981]. [110] Y. Cai, M.A. Luty and D.E. Kaplan, Leptonic Indirect [95] Y. Cui, M. Pospelov and J. Pradler, Signatures of Dark Detection Signals from Strongly Interacting Radiation in Neutrino and Dark Matter Detectors, Asymmetric Dark Matter, 0909.5499. Phys. Rev. D 97 (2018) 103004[ 1711.04531]. [111] I. Baldes and K. Petraki, Asymmetric thermal-relic [96] M. Nikolic, S. Kulkarni and J. Pradler, The dark matter: Sommerfeld-enhanced freeze-out, neutrino-floor in the presence of dark radation, annihilation signals and unitarity bounds, JCAP 09 2008.13557. (2017) 028[ 1703.00478]. [97] B.C. Lacki and J.F. Beacom, Primordial Black Holes [112] I. Baldes, M. Cirelli, P. Panci, K. Petraki, F. Sala and as Dark Matter: Almost All or Almost Nothing, M. Taoso, Asymmetric dark matter: residual Astrophys. J. Lett. 720 (2010) L67[ 1003.3466]. annihilations and self-interactions, SciPost Phys. 4 [98] S.M. Boucenna, F. Kuhnel, T. Ohlsson and (2018) 041[ 1712.07489]. L. Visinelli, Novel Constraints on Mixed Dark-Matter [113] DARWIN collaboration, DARWIN: towards the Scenarios of Primordial Black Holes and WIMPs, ultimate dark matter detector, JCAP 11 (2016) 017 JCAP 07 (2018) 003[ 1712.06383]. [1606.07001]. [99] J. Adamek, C.T. Byrnes, M. Gosenca and [114] LUX-ZEPLIN collaboration, Projected WIMP S. Hotchkiss, WIMPs and stellar-mass primordial sensitivity of the LUX-ZEPLIN dark matter black holes are incompatible, Phys. Rev. D 100 (2019) experiment, Phys. Rev. D 101 (2020) 052002 023506[ 1901.08528]. [1802.06039]. [100] G. Bertone, A.M. Coogan, D. Gaggero, B.J. Kavanagh [115] XENON collaboration, Projected WIMP sensitivity of and C. Weniger, Primordial Black Holes as Silver the XENONnT dark matter experiment, JCAP 11 Bullets for New Physics at the Weak Scale, Phys. Rev. (2020) 031[ 2007.08796]. D 100 (2019) 123013[ 1905.01238]. [116] L. Pattavina, N. Ferreiro Iachellini and I. Tamborra, [101] R.-G. Cai, Y.-C. Ding, X.-Y. Yang and Y.-F. Zhou, Neutrino observatory based on archaeological lead, Constraints on a mixed model of dark matter particles Phys. Rev. D 102 (2020) 063001[ 2004.06936]. and primordial black holes from the galactic 511 keV [117] J. Billard et al., Direct Detection of Dark Matter – line, JCAP 03 (2021) 057[ 2007.11804]. APPEC Committee Report, 2104.07634.