Dark fluxes from accreting black holes through several mechanisms

We discuss the possibility that accreting black hole systems can be sources for dark matter flux through several different mechanisms. We firstly discuss two types of systems‘: coronal thermal plasmas around supermassive black holes in active galactic nuclei (AGNs), and accretion disks of stellar-mass X-ray black hole binaries (BHBs). We explore how these black hole systems may produce keV light dark matter fluxes and find that the dark fluxes from those sources might be too weak to account for the current XENON1T excess. On the other hand, black holes can be good accelerators to accrete and boost heavy dark matter particles. If considering collisions or dark electromagnetism, those particles can then escape and reach the benchmark speed of 0.1c at the detector. We also extend the black hole mass region to primordial black holes (PBHs) and discuss the possibility of contributing to keV light dark flux via superradiance of PBHs.


Introduction
Dark matter direct search experiments have been very successfully developed to put constraints on the dark matter properties [1]. While the search for a few GeV dark matter is still going on without confirmed signals so far, lighter sub-GeV dark matter scenarios have received more attention in recent years [2][3][4][5][6][7][8][9][10][11], encouraged by the comic ray excesses [12][13][14]. The galactic center dark matter has been widely a e-mail: sichunssun@bit.edu.cn (corresponding author) considered to be the sources of the cosmic ray excesses as well as the possible explanation for the XENON1T excess [1,15,16]. The other sources for relativistic dark matter particles or axions are the Sun and stars, although such light dark matter models face many astrophysical constraints such as stellar cooling [17].
With the first detection of gravitational waves from the binary black hole mergers [18] and the Event Horizon Telescope (EHT) [19], we have new methods to probe black holes. A black hole can be a host of many processes involving the beyond the Standard Model (SM) physics [20,21]. Particles may be accreted and pairs may be produced around the event horizon, which can then be accelerated around the black hole without coupling to SM physics. In principle, the dark particles can also be thermally produced or thermalize around the black holes, and be ejected along with the standard keV X-ray radiations. The field of dark matter has very rich and diverse contents, here we just choose three difference types: light (m ∼ eV) dark photon/axion-like particles, intermediate (m ∼ keV) dark photon dark matter [73,74] and heavy (m ∼ MeV to GeV) heavier dark matter. We may further specify the details of these viable models, e.g. with/without self-interactions, and with/without the couplings to the SM models. Since all these models can be quite different, and there are so many different workable types in the literatures (there are even multiple models to explain Xenon1T), we try our best to organize these possibilities and be inclusive.
One can see related examples like how to produce dark matter relic in the early universe through black hole evap- Fig. 1 The schematic diagram of axion flux converted from the X-rays passing through the cluster of magnetic field oration in [22][23][24]. Black holes have also been discussed as the arbitrary high energy particle accelerators [25,26]. Moreover, the disks around accreting black holes can be the places for nucleosynthesis [27]. The visible accretion disk is heated through viscous dissipation of gravitational energy.
For light particles, such as dark photon or axion which may be thermally produced in the coronal thermal plasmas and accretion disks, a flux of them at keV temperature would be expected, similar to the solar source [28]. Astrophysical black holes are known to carry a large-scale magnetic field, which is the agent to tap the spin energy of black holes through the Blandford-Znajek process. The poloidal magnetic field of a black hole provides an environment for photon-to-axion conversion (γ → a) [29,30]. Besides, as depicted in Fig. 1, the cosmic magnetic field between the source and Earth can convert part of the X-ray photons to axions at the same energy scale, which is similar to neutrino oscillations. It has been found that up to 1/3 of photons can be converted, when the traveling distance inside the magnetic field is saturated [31][32][33].
If couplings between the dark sector and SM fields are quite small (For e.g. emergent dark sector in [34,35], there is only gravitational interaction between the dark sector and the standard model sector.), the thermal production of dark particle radiation from the visible thermal media, if any, is much smaller than the radiations of the Standard Model charged particles. For this case, we can look back at the traditional Penrose process [36] or superradiance, where energy can be extracted from the rotating black holes directly. Axions/dark photon vectors may also be accumulated around the black hole through gravitational superradiance effects and form clouds [62]. The process can also be independent of the black hole accretion of standard model matters.
The heavier MeV-GeV dark matter particles can also be accelerated to be relativistic, since the orbiting speed at the so-called innermost stable circular orbit (ISCO) is around the order of 0.1 c. These kinds of sources with keV kinetic energy are a lot less constrained compared to stellar sources, and provide a mechanism to accumulate large dark matter density. Moreover, black hole systems provide a good mechanism to heat up dark matter in this heavier dark matter case. The scattering process around the black hole can also heat up dark particles around it, in analogy to the visible process considering that black holes can be assumed to interact with the visible/dark sector equivalently. In this case, we need selfinteractions in the dark sector or couplings to the standard model to provide enough scattering. In this self interacting dark matter case [37], accretion disk model with the viscous hydrodynamical approximation can be applied, e.g. the double-disk dark matter [38,39]. The cross section of self interacting dark matter is bounded by σ χ /m χ 1 cm 2 /g [40]. It is however, comparable to σ T /m p 0.3 cm 2 /g for ionized hydrogens with the Thomson scattering crosssection. We will show how the dark luminosity can be as bright as the visible luminosity in a model made of dark photon and dark fermion later.
In the following sections, we discuss these processes in more details and relate them to the XENON1T excess explanations. In Sect. 2, we discuss the case with the light (m ∼ eV) dark matters, including axions and dark photons, which are converted from the X-rays passing through the cluster of magnetic field. In Sect. 3, we discuss other possibilities of the dark processes around black holes, including self-interaction dark matter from the accreting black holes, heavier dark matter boosted from accreting black holes, and superradiance from primordial black holes (PBHs). We summarize and discuss these issues more in Sect. 4.

Dark matter produced from keV environment near black holes
In this section, we discuss several environments associated to accreting black holes of different masses and characterized by keV-scale temperatures. We first consider active galactic nuclei, characterized by supermassive black holes heavier than 10 5 M . The other class is the stellar-mass black holes (3-10 2 M ) with companion stars, which are the so-called Xray black hole binaries (BHBs). There are also possibilities for primordial black holes (PBHs) such that the mass range can be extended to the region that is much smaller than three solar masses.

KeV environment around the black holes
The observations and models show that the radiations from plasmas around black holes are partly X-rays, around the keV scale. The light dark matter particles/axions can be thermally produced around visible matter and their energy is expected to be around keV, in analogy to the production in the Sun. The keV medium can produce axion/dark photon flux, as calculated in [28,59], through a few processes. These processes include dark Higgs strahlung, oscillations from visible photons to dark photons, Primakoff processes γ + Ze → Ze+a, and ABC processes for axion production. The ABC processes stand for Atomic axio-recombination and Atomic axio-deexcitation, axio-Bremsstrahlung in electron-Ion or electron-electron collisions, Compton scattering (see eg [59] for a sample of Feynman diagrams). The Eddington luminosity corresponds to the balance between gravity and radiation in the spherical plasmas i.e. [55]. It leads to G Mm where m is the test mass of a small part in the accreting matter at radius R. One then obtains the Eddington luminosity where κ E ≡ σ i /m i is the cross section due to photon scattering per unit mass. In the high energy approximation of accretion, the accreting matters are mostly composed on ionized hydrogen and the mass is dominated by the proton mass m p . The opacity is provided by Thomson scattering with the cross-section σ T . Observationally [44,45], astrophysical black holes radiate keV X-rays, which can be explained in theoretical models. If most of the emission is released in the keV band, one may estimate the X-ray flux from the Eddington luminosity in (1), which leads to the total energy flux The number flux can be obtained where E is the energy scale of each particle in the flux. For super massive black hole at kpc distances, the energy flux E E is comparable to the X-ray flux from the Sun detected at the Earth at 1 erg/(cm 2 · s) 10 8 keV/(cm 2 · s). However, most of those active AGNs are at Gpc distances. Since their luminosities are much higher, their fluxes can be comparable to those of X-ray BHBs as discussed below.

AGNs in supermassive black holes
One source of X-ray radiations is the active galactic nucleus (AGN), which is a compact region at the center of a galaxy that has extremely high luminosity. Broad-band emission from radio wave of O(10 −3 ) keV to gamma-ray of O(10 3 ) keV has been observed, which is powered by accretion of matters surrounding supermassive black holes ranging from 10 5 M to 10 9 M . For supermassive black holes with typical mass M ∼ 10 6 M , the luminosity is L ∼ 10 44 erg/s. At the galactic center of the Milky Way galaxy, there is a supermassive black hole with mass 4 × 10 6 M which powers the compact radio source Sgr A* and is about 8 kpc from Earth [174,175], which might be a possible source for the dark matter fluxes. In our Galactic center, the luminosity of Sgr A * [178,179] and other low luminosity AGNs is in the order of 10 −3 L E . For very active galaxies, the bolometric luminosities can be higher than 10 2 L E .

X-ray black hole binaries (BHBs)
Except for the matter accretion of AGNs, around stellar-mass black holes with stellar companions, the accretion of matter is also happening. They are so-called X-Ray black hole binaries (BHBs). The candidates include Cygnus X-1, XTE J1650-500 and GX 339-4 [56][57][58]. A density cusp of quiescent Xray binaries can also be found in the central parsec of the Galaxy [181]. There are two kinds of X-ray spectra of the accreting black holes. In the so-called high-soft state, the Xray spectrum is composed of a strong and narrow peak at a few keV and a soft power-law spectrum up to hundreds of keV. The low-hard state peaks below 1 keV while the Comptonized component is extended up to 100 keV. The luminosity of the high-soft state can be near or much higher than L E , while the low-hard state is usually two orders of magnitude below the Eddington luminosity L E . The widely considered accretion flow model is the socalled Shakura-Sunyaev model [46], which is used for the observed emission energies less than 10 keV. This thin disk model was also worked out by Lynden-Bell, Pringle and Rees [47], and a general relativistic treatment for the inner part of the disk can be found in Novikov and Thorne's [48,49]. For the purpose of order of magnitude estimation of this paper, we take the Shakura-Sunyaev model as an example to show that they can produce keV emission. There are many more variants of accretion models. As long as they can produce X-rays, our following discussion remains relevant. We use the notations r = R/R s and m B H = M/M . The local radiation energy flux from unit surface at the radius R of the disk is determined by the gravitational energy release  [50].
If considering the critical Eddington luminosity with L = f sṀ c 2 = f e L E , the typical energy flux is around Here the Eddington luminosity has been 3.8 × 10 33 erg/s is the solar luminosity. In the following calculation, we will simply take f s 0.1 and f e = 1 for the estimation [51]. The Shakura-Sunyaev disk is optically thick, which can be roughly approximated as blackbody. From the law of thermal radiation σ s T 4 The typical value of T 0 is calculate from T 0 ≡ (E 0 /σ s ) 1/4 . In Fig. 2, we make the contour plot of temperature T r in (3) as a function of locations in the disk r = R/R s and black hole masses m B H = M/M . One can see that the peak temperature of the disk is around T r 1.2 keV/m 1/4 at the radius r 4. And in Table 1, we list two benchmark points for black hole sources. The gravity of primordial black holes is too weak to have significant accretion, so we leave the discussion to section 3.3.
With the simple model of the thin accretion disk with the temperature distribution in (3), in the X-ray black hole binaries, the temperature of the disk peaks at energies around 1 keV. The temperature of disk around the supermassive black Table 1 Here we list two benchmark points for accreting black hole sources, with the Eddington luminosity and the peaked temperature in thin accreting disk in (3) Types of X-ray sources AGNs BHBs Benchmark masses (M ) hole peaks at around 10 −2 keV. For AGNs, more complicated models are required. From the observations of a couple of AGNs, e.g. 1H 0707-495 [52][53][54] and analysis, the X-ray emission comes from the corona surrounding the black hole of an AGN as well as the reflection and thermal emission of the rays from the accretion disk. It has multiple features including reflection from accretion disk from 0.5 to 1 keV, corona X-ray around 1-5 keV and Iron Kα emission line peaked at 6.8 keV. The number flux in the spectrum from sources [52][53][54]

Generating keV dark matter
The temperature of AGNs or X-ray binary disks can be as high as keV, which is close to the temperature in the core of the Sun. However, the density of the accretion disk for AGNs with mostly ions and electrons, ρ disk ∼ 10 −7 g/cm 3 [60], is much smaller than the ρ HB ∼ 10 4 -10 6 g/cm 3 (HB stands for the stars in the horizontal branch with the solar masses), and is model-dependent. The nuclear reactions take place in the core of the Sun and the reaction rates are highly temperature dependent. Black hole accretion disks can have high temperatures generated by the high accretion rates, such that dark particles/axions can be generated in this environment.
Similar to the solar environment where the axion is converted via γ → a [17], AGN corona (or solar corona) and X-ray BHB accretion disks can have a temperature of ∼ 10 7 Kelvin and produce keV dark matter. Thermal emission is a process where the collected excitations, plasmons are converted to the light-dark matter of the comparable energy. The frequency of the plasmons is dominated by the medium temperature at a high temperature. However, since the production rate is related to density so that the production rate from AGNs has been found to be quite low [61]. Only the solar corona or BHB accretion disks may reach an considerable axion flux, which is usually bound by the visible luminosity. Considering that the luminosity of dark particles/axions from the Sun is bound by 10% of visible luminosity [28,59], we assume that the dark luminosity L D is comparable to the visible luminosity.
For AGN corona, one needs to consider the alternative mechanism as shown in Fig. 1. It has been proposed that the cosmic magnetic field may convert part of the X-ray photons to high energy axions [31,32], where up to about 1/3 of photons can be converted, similar to neutrino oscillations. Thus, we simply start from the total energy flux (2) with Eddington luminosity in (1), and make the major assumption that the total dark luminosity is characterized by the Eddington luminosity L D L E .
Notice here that this assumption L D ∼ L E is rather ad hoc. If the dark particles are produced thermally, the couplings with the standard model particles are bound by the Fig. 3 The contour plot of the bound of total dark flux D in Eq. (4) with the Eddington luminosity and E D 3.5 keV, as a function of the distances R/kpc and masses M/M of black holes stellar cooling constraints which further constrains the production rate. If the dark flux is not produced thermally, then Eddington luminosity can be a benchmark point for the dark luminosity, assuming the dark sector is accreted by the black holes and balanced by some self-interaction-induced scatterings within the dark sector. If considering the axion fluxes converted from the X-ray photons under the cosmic magnetic field, this bound L D ∼ L E is also plausible.
In summary, the total number flux is given by In Fig. 3, we make the contour plot of dark flux D as a function of locations in the disk R/R s and black hole masses M/M , which will be helpful to estimate the benchmark parameters below.

Possible signals and benchmark parameters
Below we discuss the light dark matter candidates that can be produced around black holes and estimate their flux. Taking into account of known AGN or BHB sources, we find that the flux is too low to explain the XENON1T excess. However, it is interesting to estimate the benchmark parameter of the flux, which may be detected in the future experiments. Before that, we summarize the general result on the detection. One possibility is to consider the electron recoil, and the total number of signal events is given by N s = i × σ ie × T X C X . Here C X = Z e n X V X and i is the incoming flux with dimension [cm −2 s −1 ], σ e is the cross section with dimension [cm 2 ] and T X is the operation time. C X is the total number of effective electrons in the detector, which can be calculated by the product of the effective number Z e of elec-trons in each Xenon atom that undergo recoils, the number density of Xenon atoms n X and the total fiducial volume V X . On the XENON1T excess with N s ∼ O(100), e.g. in [69], the required product of the number flux and cross section has been found to be i × σ i O(10 −35 )/s. In the following, we will use this relation to estimate the required flux and then the distance to the sources with (4), considering the accreting black holes as the sources.

Axion/axion-like particles (ALP)
One possible explanation of the claimed XENON1T excess is the solar axion/ALP. It requires the coupling to electrons without/with photon coupling through the inverse Primakoff process [70], a+Ze → γ +Ze or inverse Compton scattering a + e → e + γ , with the couplings in the Lagrangian Since the solar axion luminosity is bounded by stellar cooling, the couplings have been constrained as g aγ <10 −10 GeV −1 or g ae <10 −12 [28]. The kinetic energy of axion/ALP considered for Xenon here is 2-4 keV with the axion mass m a 1eV. However, this explanation is in tension with astrophysical stellar cooling constraints [17], because the solar axion production is well studied.
What's more, the couplings bound from the stellar cooling constraints also lead to the cross sections bound. The minimum required number flux density of the incoming axions for the XENON1T excess is around a 10 11 /(cm 2 s). Here instead, it is possible to consider accreting black holes as the sources. If we consider the keV axion in the dark flux bound by (4), with the help of Table 1, one can obtain the minimum distance of the three types of black holes. It is tempting to explain the Xenon 1T axion with the black hole sources rather than the Sun, provided that the thermal plasma around the observed black hole is at the keV temperature.

Light dark photon
The dark photon production in the Sun is mostly from dark Higgstrahlung and visible photon oscillations with real or virtue photons. The dark photon and dark fermion sector is given by the Lagrangian density The typical solar flux is γ 10 10 /(cm 2 s) for mixing ∼ O(10 −14 ) [59]. It was recently pointed out in [73] that Table 2 Here we list two benchmark points with some known black hole binary sources combined in [71,72]. We assume the dark particle luminosity of the host black hole is comparable to the visible luminosity L D ∼ L E , and those sources are around kpc distance from the earth
Here if we replace the Sun by other sources such as two types of accreting black holes in our Table 1, with different couplings, the flux aimed at the XENON1T excess received on Earth is γ 10 11 /(cm 2 · s). Considering a few known keV X-ray sources combined in [71,72], and with the help of Fig. 3 or Eq. (4), one can summarize in Table 2 below with ∼kpc distance black holes with the assumed dark luminosity L D ∼ L E [71,72]. However, the maximal intensity of the flux D is around 10 3 /(cm 2 · s), which is a few orders lower for the Xenon1T excess, but may be considered for the future astrophysical observation.
Here we would also like to comment on the X-ray telescope constraints, similar to the discussions in [75]. In general, these benchmark points in Table 1 do not contradict the X-ray observation, since we basically assume the dark flux/axion flux is converted from the observed photons in Xray. So the dark flux is around the same order as the X-ray observations and does not exceed the observed limits. In the other words, the explanation that the Xenon1T signals come from the known X-ray sources of AGN and BHB conversions is basically excluded. In Fig. 5, we use the green dashed line to denote the bound from the X-ray background in [67]. In the next section, we will consider the dark fluxes from the accreting black holes directly.

Dark matter flux and black hole dynamics
In this section, we will discuss how heavier dark matter may accumulate, and accelerate around the black hole,as well as the benchmark parameters compared with the Xenon1T experiment.

Dark luminosity and accreting black hole
Now we consider a different mechanism to produce an outflow of dark matter particles. The heavier non-relativistic MeV-GeV dark matter cannot be thermally produced in the plasma. They can however get boosted around the event horizon or captured/accreted by the black hole. Especially the black holes have been considered to act like particle accelerators to boost the particles to arbitrary high energy and collide [25,26]. However, for the traditional collision-less cold dark matter, the dark matter particles can have a larger speed near the black hole, but that does not mean they can escape. If we assume the collision can happened, then some of the products may escape the near horizon region of black holes. That means that we need to consider self-interacting dark matter in order to account for a dark flux [37]. Another possibility is the charged DM model with the dark electromagnetism (mediated by the dark photon) as discussed in the lagrangian (6). One may be able to invoke shocks or magnetic reconnection processes to accelerate particles via the first-order Fermi acceleration mechanism. All these processes require charged dark matter particles.
At the galaxies scale, the self-interacting dark matter model is introduced to solve the problems in small scale structures [37], such as the core-cusp problem, the rotation curve's diversity problem, the missing satellites problem, etc. The double-disk dark matter is also motivated [38,39], such that the accretion disk model with viscous hydrodynamical approximation can be applied. For example in [30], the hypothetical dark photon was proposed to mediate the interactions of dark matter particles through (6). Neither of them interacts much with the standard model particles.
It is interesting to estimate the dark luminosity following the similar derivation of Eddington luminosity in (1). In [180], it has been found that that accretion of dissipative dark matter onto AGN contributes to the growth rate of the black hole. Here we assume the dark matter particles are scattered upon the dark photons via the cross section σ χγ due to the interaction in (6), we can have the dark luminosity where κ χ ≡ σ χγ /m χ . For ionized hydrogen, one has κ E = σ T /m p ∼ 0.3cm 2 /g. From numerical simulations of merging cluster galaxies [37,40], one bound of the self-interaction cross-section of dark matter particles is σ χχ m χ 1cm 2 /g. Here the cross sections σ χχ and σ χγ are both related to the dark photon fine-structure constant α , which is bounded by the relic abundance [30]. Hence, it is reasonable to assume κ χ ∼ κ E , and we can have the dark luminosity as bright as the Eddington luminosity L χ ∼ L E . If considering the case of the self-interacting dark matter without dark radiation, the gravitation could be balanced by the pressure of the dark particle gas κ χχ ≡ σ χχ /m χ , then again we can reach For the much lighter axion, the self interaction is very small [76]. If we take the cross section σ aγ σ aγ and usual Fig. 4 The schematic diagram of black hole and dark star(or axion clump) binary system. If dark matter particles have enough selfinteractions, they can clump together to form a dark star. If such dark star happens to be accreted to one black hole, it is possible to form the dark matter disc and dark matter flux QCD axion at μeV scale, then κ a ≡ σ aγ m a will be many orders of magnitude smaller than κ E , which means that the luminosity L a = L E κ E κ a can be much higher. However, this kind of system is usually unstable and may lead to bursts of photons or axions after a short time. This process may be an interesting source to search for axions and will be interested in further exploration. For example, one may consider a possible black hole dark star binary system which is depicted in Fig. 4. If dark matter particles have enough self-interactions, it can clump together to form a star. We can see the examples of e.g. boson stars, axion stars or moduli stars [21,41]. If such a dark star happens to be accreted by a black hole, a dark disk may form around the black hole. Similar to an X-ray binary, this object, if exists, could be a power source of dark matter flux, without necessarily being observed in the X-ray band, such that it could be quite close to Earth. Since these systems are mostly invisible, it is possible for the black hole -axion star binaries to produce the required relativistic flux of axions with self-interaction. More detailed modeling of black hole dynamic processes is required for a similar production rate calculation around the black holes. In the following, we list a few mechanisms to produce the dark matter flux.

Boosting dark matter by black holes
There are a few interesting circular orbits very close to the event horizons of the black holes. The photon sphere near a black hole is where photons travel in circular orbits. For a Schwarzschild black hole with the event horizon located at R s = 2G M c 2 , the photon sphere is located at 3R s /2. And the innermost stable circular orbit of the massive particles is located at 3R s . For a rotating black hole, it has corrections due to the rotational effects [42] which can be diagnosed using BH images obtained from EHT (see e.g. [19]). In general, for a rotating and charged black hole with angular momentum J and charge Q, the outer event horizon is R + = R s (1 + 1 − a 2 − q 2 )/2, where a = J/M and q = Q/M. Here when a 2 + q 2 = 1, R + = R s /2 is the extremal radius. Normally a 2 +q 2 < 1 is required by the cosmic censorship hypothesis to avoid naked singularities. The Hawking temperature sourced by the quantum radiation of a Schwarzschild black hole is around T H 6×10 −8 K M M B H [43]. For stellar-mass black holes around 3-10M , one has T H ∼ 10 −8 K . For super massive black holes ( 10 6 M ), one has T H 10 −13 K . Both are much lower than the observed temperature of the cosmic microwave background temperature 2.7 K. However, some black holes are surrounded by very hot coronal thermal plasmas and accretion disks, with the inner circle located around ISCO R I = 3R s . For the primordial black holes with much smaller masses, the Hawking temperature can be as high as keV∼ 10 7 K , and we leave the discussion of this case to the next section.
The heavier boosted dark matter from MeV to GeV mass with around 0.1c can also explain the excess as in [69,87] with a number flux of χ ∼ 10 −6 /(cm 2 · s) in an NFW profile [78] from the Galactic center to produce Xenon excess considering elastic scattering between the dark matter particle and an electron χ + e → χ + e . For the non-relativistic case, the transferred recoil energy Considering the mass of electron as m e 511keV and v e ∼ 10α e 10/137 as in [87], for the XENON1T excess between 2-4 keV, a required benchmark velocity of the heavy dark matter particles is v χ ∼ 0.1c, with the boost factor , where σ eχ is the required cross section and here we just use it as the benchmark flux for our dark sources [69,87].
Black holes can be good accelerators to accrete and boost heavy dark matter particles. For the boosted dark matter to escape and reach the benchmark speed of 0.1c at the Xenon detector, we need to assume extra ingredients. One may consider the semi-annihilation dark matter χ χ → χφ, where χ and χ are the heavy DM particles with the same mass, and φ is a light particle coupled to the SM sector [88]. The boost factor of the final states turns out to be where γ i is the boost factor of χ in the center of mass frame. When γ i = 1, we reach the limit γ χ = 1.25 and v χ = 0.6c. For the twocomponent DM model χ AχA → χ BχB studied in [69] with initial boost factor γ A , the boost factor for χ B is simply given by γ B = γ A m A /m B [2], which can go beyond 0.6c easily. Thus, for both cases, the velocities of the produced dark matter particles χ B can reach or go beyond the escape velocity outside ISCO. The dark accretion disk also provides a higher density environment to enhance semi-annihilation and Table 3 Here we list the other two benchmark points for the XENON1T explanations with black hole sources. The parameter choices are based on [73] and [69]. We assume the dark particle luminosity is L D 10 38 erg/s. The last row shows the benchmark distances that can produce the Xenon1T excess. For the dark photon dark matter(DPDM), there is no requirement. The BHs help to increase the flux γ and thus relax the benchmark choice of mixing parameter
Here we again assume that the dark luminosity is characterized by the Eddington luminosity. In Table 3, we list heavier dark matter candidates, such as the dark photon dark matter [74], where the dark photon is produced non-thermally e.g. from inflation decay and acting like cold dark mater, or boosted dark matter with a black hole source. The parameter choices are from [73] and [69], and we assume that the dark particle luminosity of the host black hole is L D 10 38 erg/s.
There are also a lot of discussions on the cosmological origin of the dark photons (see e.g. [77]), similar to the dark photon dark matter model. The mass of the dark photon dark matter is generally heavier than the dark photon case we discussed in Sect. 2. It was found in [73] that the local dark photon dark matter with mixing angles ∼ O(10 −16 ) can fit the recent Xenon1T excess well and satisfy the astrophysical constraints. Here one can find for the local 3 keV dark photon dark matter with the density ρ γ 0.3 GeV/cm 3 and 10 −3 c, the flux can be as high as γ ∼ 10 12 /(cm 2 · s). The intrinsic mass m γ ∼ 3keV is built into the model [73]. In fact, if we consider the black hole sources, the flux can be enhanced, and mixing parameters and m γ can then be relaxed accordingly. In Fig. 5, we plot the Xenon1T signals, background and the predicted signals from two scenarios discussed in Table 3.
With the benchmark parameters and Fig. 5, we would also like to comment on the dark matter constraints from direct detections. For the boosted dark matter case, the black hole accumulation and boosting mechanism does not depend on the DM-nucleon scattering cross section [69], different from solar capture mechanism, therefore the current cross section σ eχ ∼ 10 −29 cm 2 does not subject to the other DM direct detection experiments constraints in general.
For the dark photon dark matter in [73], the star cooling constraints are stronger than other direct detection limits e.g. XENON100 [79], XENON1T [80], SuperCDMS [81], Fig. 5 The data points and background curve B 0 are reproduced from the Xenon1T data and the background model in [1]. The red dashed line is produced from the boosted dark matter (BSD) model accelerated by the black holes. This line is schematically produced as in [69], with v χ 0.1c and other benchmark parameters in Table 3. Here we smear the recoil energy in equation (8) with the Rayleigh distribution, considering the detector resolution. The blue dashed line denotes the signals of the dark photon dark matter(DPDM) with m A = 2.8 keV with resolution of 1keV as in [73] and other benchmark parameters in Table 3. The green dashed line denotes the axion flux bounded by the X-ray background in [67]. The event number is shown in the plot after being multiplied by the factor of 10 7 CDEX [82], SENSEI [83,85,86] and CRESST-III [84]. But in our case, the flux of dark photon dark matter is enhanced by black holes with γ 10 12 /(cm 2 ·s) around m A ∼ 2.8keV. The mixing parameter can then be relaxed 10 −15 . Thus, our benchmark parameters choice in Table 3 is safe from direct detection constraints.

Dark superradiance of black holes
On the superradiance of black hole, one kind of interesting objects is the gravitationally bound "atoms" with light bosons (e.g. axions) or vectors [62]. Similar to the hydrogen atoms, those light particles occupy different energy levels of the rotating black holes. Although the kinetic energy in the superradiance clouds is non-relativistic, these accumulated light particles can be accelerated through different energy extraction processes around the black hole. Notice that in gravitational atoms, when the Compton wavelength of dark particle is comparable to the size of the black hole, the number of dark particles can grow exponentially from extracting energy of the rotating black holes. When the attractive selfinteractions become stronger than the gravitational binding energy, the cloud collapses and the explosion induces an outflow [62]. If considering the decay of axion through a → γ γ , it has been shown in [20,21] that the μeV axion clouds around 10 −5 M PBHs may induce fast radio bursts (FRBs), and the luminosity can be as high as 10 39 erg/ms, more than enough to power observed FRBs.
It is also interesting to consider the bursts of dark sectors with the similar mechanism. If considering the channel of φ → γ γ , the dark photon at 3.5 keV requires m φ 7 keV. Considering the reduced Compton wavelength λ φ =h/(m φ c), the dimensionless number describing the dark matter superradiance cloud is α s ≡ R s  3 , which is much longer than the age of the universe at 10 10 years. If the particles in Hawking radiation of PBHs are mainly in the dark sectors [64], dark radiations at the keV energy scale could also be one candidate source for dark fluxes.
The PBHs, the luminosity of Hawking radiation L B H ∼ 10 −25 L is quite small compared to the solar luminosity. To estimate the total flux, one needs to consider the integrated luminosity with the distribution of PBHs. The related X-ray constraints have already been discussed in the PBHs as the dark matter scenarios with the Hawking radiation in [65][66][67]. It is also interesting to consider whether we can find axion signals from PBHs in X-ray telescopes. See more recent discussions in [68]. The axion flux converted from the known X-ray observation is plotted as the dashed green line in Fig. 5, which is a few orders lower than the Xenon1T excess.

Discussion and summary
In this study, we propose that accreting black holes could be the natural sources to produce the keV-order light dark matter and boost the heavier dark matter. We discuss several mechanisms to produce, accumulate and accelerate the dark photon, axion, even heavier dark matter around the black holes, motivated by the fast dark matter explanations of XENON1T excess, e.g.  .
The X-ray black hole binares with stellar masses, or AGNs in the super massive black holes can be possible hosts (see also a recent event of an intermediate mass black hole [173]). However, the fluxes observed in AGNs and BHBs are many orders of magnitude smaller than the flux required for the XENON1T excess since their distances to Earth are > 1 kpc. Although the black hole corona and disk productions of axions and dark photons require some detailed modeling of accretion processes, the observed X-ray emissions provide benchmark points around the 1-8 keV range. There are many bright X-ray sources in the sky that originate from accreting or spindown-powered neutron stars. Supernova remnants, Gamma ray bursts (GRB) afterglows can also produce bright X-rays although the density in the emission region is extremely low and the radiation process is synchrotron radiation. These systems could also be interesting subjects of study in the future.
One viable mechanism for the current Xenon1T excess is black holes as the dark matter boosters. We discuss and calculate the boosted factor and dark luminosity here. The black hole accretion can also accumulate dark matter to achieve a higher dark matter local density. The black hole candidates such as primordial ones (e.g. the smaller "Planet 9" with 10 −5 M ) with closer locations may also be the sources to boost dark matter [176,177]. The other possible region is the micro PBHs around 10 −15 M which can emit X-rays through thermal radiation. Although the luminosity is quite small, it is plausible to assume the radiations partly belong to the dark sector.
The dark cosmic fluxes created around black holes can be interesting and less explored targets for direct detections, in light of the multi-messenger detection experiments for black holes such as gravitational waves and EHT. The constraints are also less severe compared to the stellar productions. Especially with the help of the multi-messenger detections, we hope the dark processes around black holes can be understood and modeled better in the near future.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .