Ruling out dark matter interpretation of the galactic GeV excess by gamma-ray data of galaxy clusters.

Recently, some very tight constraints of annihilating dark matter have been obtained from gamma-ray data of the Milky Way and Milky Way dwarf spheroidal satellite galaxies. In this article, we report that there are two excellent galaxy clusters (A2877 and Fornax) which can provide interesting constraints for annihilating dark matter. The lower limits of the dark matter mass for the thermal relic annihilation cross section are 25 GeV, 6 GeV, 130 GeV and 100 GeV respectively for the e + e -, μ + μ -, τ + τ - and [Formula: see text] channels. For some configuration of our working assumptions, our results improve the Fermi-LAT upper limits of annihilation cross sections by a factor of 1.3 - 1.8 for wide ranges of dark matter mass for e + e -, μ + μ - and [Formula: see text] channels, and a factor of 1.2-1.8 for τ + τ - channel with dark matter mass ≤100 GeV. These limits basically rule out most of the existing popular dark matter interpretation of the GeV excess in the Milky Way.

Introduction

In the past few years, gamma-ray observations of the Milky Way initiated a hot research area about dark matter annihilation. Many studies[1-5] point out that the GeV excess in the Milky Way can be best explained by dark matter annihilation with dark matter mass  GeV via τ + τ − or channels. However, later gamma-ray observations of the Milky Way dwarf spheroidal satellite (MW dSphs) galaxies by the Fermi-Large Area Telescope (Fermi-LAT) put strong constraints on the annihilation cross sections for different channels[6,7]. These limits give a certain tension for the dark matter interpretation of the GeV excess in the Milky Way. Besides, some recent studies suggest that the GeV excess in the Milky Way can be explained by pulsars’ emission[8-13] or depletion of gamma rays in molecular clouds[14]. Furthermore, recent results of Fermi-LAT gamma-ray detection strongly disfavour the dark matter interpretation of the GeV excess[15]. Therefore, the GeV excess in the Milky Way becomes a very controversial issue in dark matter astrophysics now.

Although most of the gamma-ray observations are now focusing on the Milky Way and the MW dSphs galaxies, some studies start to realize that galaxy clusters are also good candidates to constrain annihilating dark matter. For example, Ando & Nagai[16] find that the gamma-ray data from the Fornax cluster is quite good to constrain the parameters of dark matter annihilation because its cosmic ray emission is not very significant. The cross section constrained by using the Fornax cluster data is <συ>≤ (2–3) × 10 cm3 s−1  [16]. Besides, Huang et al.[17] use three years of Fermi-LAT data of several nearby galaxy clusters to constrain annihilating dark matter. The upper limits of annihilation cross section are also of the order 10−25 cm−3 s−1. Later studies of gamma-ray emission from galaxy clusters and galaxy groups give more stringent constraints on dark matter parameters. For example, Ackermann et al.[18] and Lisanti et al.[19] exclude the thermal relic annihilating dark matter via channel for m ≤ 40 GeV and m ≤ 30 GeV respectively. Also, using new Fermi-LAT data of nearby galaxy clusters can provide stringent constraints for gamma-ray line feature of annihilating dark matter[20]. Besides gamma-ray data, by using radio data of galaxy clusters, the annihilation cross sections for some channels can be constrained down to  cm3 s−1  [21-23], though the uncertainties due to magnetic field strength are somewhat significant. Overall speaking, these limits are quite close to or even tighter than some of the parameters for the dark matter interpretation of the GeV excess in the Milky Way.

In this article, by using the latest gamma-ray data of galaxy clusters from Fermi-LAT, we show that the A2877 and Fornax clusters can give very stringent constraints for annihilating dark matter. Our results improve the current limits of annihilation cross sections for wide ranges of dark matter mass. These limits are tighter than the parameters involved for the dark matter interpretation of the GeV excess in the Milky Way.

Method

Generally speaking, there are two different kinds of gamma-ray emissions for annihilating dark matter in galaxy clusters. One is ‘direct emission’ and the other is ‘indirect emission’. ‘Direct emission’ means the gamma-ray photons produced when dark matter particles annihilate each other. Since dark matter annihilation also produces electron and positron pairs, these electron and positron pairs can interact with the thermal electrons in galaxy clusters via Bremsstrahlung cooling. This process would emit gamma-ray photons if the energy of the electrons or positrons is high enough. This kind of emission is regarded as ‘indirect emission’.

Direct emission

The total number of photons for direct emission with energy greater than E 0 produced by annihilating dark matter within a galaxy cluster can be calculated bywhere ρ is the dark matter density profile, R 200 is the radius of a galaxy cluster, and dN /dE is the energy spectrum of gamma-ray produced per one annihilation, which can be obtained in[24]. Based on the sample of galaxy clusters in[25] and the observed limits of gamma-ray flux Φobs obtained in[26], there are two galaxy clusters (Fornax and A2877) which have the largest ratio I/Φobs, where . Larger the ratio of I/Φobs can give tighter constraints for annihilating dark matter. Therefore, we choose the Fornax and A2877 clusters to perform analysis. By assuming the dark matter density profiles of the A2877 and Fornax clusters follow the Navarro-Frenk-White (NFW) density profile ρ = ρ r 3[r(r + r )2]−1  [27] and using the corresponding parameters in[28], we can get the total photons produced with energy greater than E 0 for direct emission.

On the other hand, recent numerical simulations show that substructures’ contribution for dark matter annihilation is very significant. High dark matter density of the substructures can boost the annihilation signal to more than 10–100 times for galaxy clusters[29-32]. However, the actual value of the boost factor in a galaxy cluster is controversial. Some high-resolution simulations can get the boost factor (e.g. Phoenix simulation project)[29] while other studies point to a lower value of the order [30,31]. Moreover, the effect of tidal interaction in galaxy clusters is also controversial. Some studies point out that tidal interaction would give a lower boost[31] while some studies suggest that tidal stripping can increase the boost factor by 2–5 times[33]. Here, for a conservative estimation, we neglect the effect of tidal interaction and take the lower bound of the boost factor, which can be described by a simple parametrization[31]:where b = [−0.442, 0.0796, −0.0025, 4.77 × 10−6, 4.77 × 10−6, −9.69 × 10−8] and M 200 is the total mass of a galaxy cluster. The accuracy of this parametrization is better than 5% in the mass range [31]. For the A2877 and Fornax clusters, the total mass M 200 are and respectively[28]. Therefore, the parametrization can be applied and the corresponding boost factors are 35 and 30 respectively. As discussed in[32], the boost factor of the order are standard values for galaxy clusters. By taking these boost factors, we can calculate the direct emission flux for each galaxy cluster by , where D is the distance of the galaxy cluster. Here, we are using the point-source approximation, which can be justified because we will compare the calculated fluxes with the ‘point-source upper limits’ in[26].

Indirect emission

Dark matter annihilation also produces electron and positron pairs. These electrons and positrons would diffuse and cool down to a lower energy by synchrotron cooling, inverse Compton scattering, Bremsstrahlung cooling and Coulomb cooling. Since the diffusion coefficient for a typical galaxy cluster is about cm2 s−1  [34], the diffusion time scale for a 1 GeV electron is about s. However, the total cooling time scale for the above cooling processes is just s. Therefore, the effect of diffusion is nearly negligible in galaxy clusters. The equilibrium electron energy spectrum can be simply obtained by[21] where is the integrated energy spectrum of electrons produced per one annihilation and b(E) is the cooling function which is given by[34] where n is the thermal electron number density in cm−3, E GeV = E/1 GeV, B is the magnetic field strength in μG, b IC = 0.25 × 10−16 GeV s−1, b syn = 0.0254 × 10−16 GeV s−1, b cou = 6.13 × 10−16[1 + log(γ/n)/75] GeV s−1, b brem = 1.51 × 10−16[log(γ/n) + 0.36] GeV s−1 and γ is the Lorentz factor of a high-energy electron or positron. We follow the β-profile to model the thermal electron number density profile in a galaxy cluster:where n 0, r and β are fitted parameters, which can be obtained in[25]. For the magnetic field profile, we follow[21] to usewhere we have assumed B 0 = 4.7 μG and η = 0.5 for both galaxy clusters (the effect of the uncertainties in B will be discussed below)[21].

For Bremsstrahlung emission (indirect emission), the total gamma-ray photons emitted per second for indirect emission with energy greater than E 0 can be given by[35] where

By taking the same boost factors, we can calculate the indirect emission flux for each galaxy cluster by , with the point-source approximation.

Data availability

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Result

We can now obtain the total gamma-ray flux Φ = Φdir + Φind emitted from each of the galaxy clusters. Based on the observational results from Fermi-LAT[26], we obtain the fluxes for three lower limits of gamma-ray energy E 0 = 0.5 GeV, E 0 = 1 GeV and E 0 = 10 GeV. Also, we use the upper limits of the fluxes for these values of E 0 to constrain E 0 (see Table 1 for the point-source approximation upper limits for each E 0).

Table 1

The observed gamma-ray flux 95% C.L. upper limits for the A2877 and Fornax clusters (point-source approximation)[26]. All the units for the gamma-ray fluxes are in 10−11 cm−2 s−1.

	A2877	Fornax
Flux upper limits for E ≥ 0.5 GeV	8	11
Flux upper limits for E ≥ 1 GeV	3.6	4.9
Flux upper limits for E ≥ 10 GeV	0.19	0.26

In Fig. 1, we plot the upper limits for the annihilation cross section e + e − for 4 popular channels (e + e −, μ + μ −, τ + τ − and ). We can see that the gamma-ray flux upper limits can give very stringent constraints for the annihilation cross section. In standard cosmology, the thermal relic annihilation cross section is<συ> = 2.2 × 10−26 cm3 s−1  [36]. The lower limits of dark matter mass for the thermal relic annihilation cross section are 25 GeV, 6 GeV, 130 GeV and 100 GeV for the e + e −, μ + μ −, τ + τ − and channels respectively. Although these limits are nearly the same as the current one obtained by Fermi-LAT data of the MW dSphs[6,7,37], overall speaking, we get tighter constraints for the upper limits of annihilation cross sections, except for the τ + τ − channel with m > 100 GeV. Our results improve the upper limits of annihilation cross sections by a factor of 1.3–1.8 for wide ranges of m for e + e −, μ + μ − and channels, and a factor of 1.2–1.8 for τ + τ − channel (m ≤ 100 GeV) (see Fig. 2).

Figure 1

The upper limits of the annihilation cross section for 4 channels. The lines in red and blue represent the limits based on the data from A2877 and the Fornax cluster respectively (dotted line: E 0 = 0.5 GeV; dashed line: E 0 = 1 GeV; solid line: E 0 = 10 GeV). The black solid lines and dashed lines represent the upper limits of the annihilation cross section based on the Fermi-LAT data of the Milky Way dwarf spheroidal satellite galaxies[6,7]. The black dotted lines represent the thermal relic cross section <συ> = 2.2 × 10−26 cm3 s−1 [36].

Figure 2

The ratio of the upper limits of the annihilation cross section obtained in this study to the Fermi-LAT upper limits in[6,7] for 4 channels.

Furthermore, our results are more stringent than that obtained by other recent studies with the data of galaxy clusters. For example, recent studies using the data of Virgo cluster[18] and Coma cluster[32] obtain the lower limit of GeV for dark matter annihilating via channel. Our results improve this limit by a factor of 2.5. Also, our results improve the limits obtained in[17,21] (using the data of other clusters, such as Ophiuchus, A2199 and AWM7) by nearly an order of magnitude. Nevertheless, by using radio data of galaxy clusters, some recent studies can give more stringent constraints of annihilation cross sections. For example, Beck & Colafrancesco[22] obtain m ≥ 200 GeV for the thermal relic annihilation cross section via or τ + τ − by using the radio data of the Coma cluster. However, these studies constrain dark matter annihilation by synchrotron radiation which may suffer from large uncertainties of magnetic field strength. Although our study also includes the consideration of magnetic field strength, we consider the effect of Bremsstrahlung radiation due to magnetic field, but not synchrotron radiation. As discussed below, the uncertainties of magnetic field do not affect our result significantly. Therefore, using gamma-ray data is more reliable than using radio data for constraining annihilating dark matter, especially in using data of galaxy clusters. Based on the above discussion, we suggest that the Fornax and A2877 clusters are better candidates for constraining annihilating dark matter.

Moreover, our constraints here can safely rule out all of the existing popular dark matter interpretation of the GeV excess in our Milky Way (via τ + τ − or )[2,4,5] (see Fig. 3). Besides, some recent models of dark matter interpretation of the positron excess and antiproton excess in our Milky Way are also ruled out (e.g. τ + τ − and channels, see Fig. 3)[38,39].

Figure 3

The solid lines represent the upper limits of the annihilation cross sections for the e + e −, μ + μ −, τ + τ − and channels. The data points with 1 σ error bars are the results obtained in[2,4,5,38,39] for the dark matter interpretation of the GeV excess and positron excess. Here, different colors of solid lines and data points represent different annihilation channels. The dotted line is the thermal relic cross section.

We have also checked with the effects by varying the magnetic field strength, dark matter profile and boost factor for both galaxy clusters. First of all, the magnetic field strength mainly affects the indirect emission for low energy (E ≤ 10 GeV). The contribution of indirect emission varies with dark matter mass mildly. The largest values of the ratio Φind/Φdir for the e + e −, μ + μ −, τ + τ − and channels are 26%, 19%, 3% and 1% respectively (mainly for E ≤ 1 GeV). In the above calculations, we used a typical central magnetic field strength B 0 = 4.7 μG and a typical index η = 0.5. We have tried another extreme value B 0 = 25 μG[21] and change the index η to 0.4, we find that the overall upper limits of annihilation cross sections are just larger by E ≥ 10%. Therefore, the effect of the systematic uncertainties of magnetic field strength (or contribution from indirect emission) is not very significant. The reason is that the most stringent constraints mainly come from the data of gamma-ray flux E ≥ 10 GeV. Changing the magnetic field strength does not affect too much for the gamma-ray flux for E ≥ 10 GeV.

For the dark matter profile, we have used the NFW profile to model the density profile of galaxy clusters. We have tried another profile, Einasto profile[40], to model the Fornax and A2877 clusters as some recent studies suggest that the Einasto profile is a better description of dark matter halos on large scale[41]. We finally get tighter constraints after using the Einasto profile (see Fig. 4). Here, we only show the effect for the Fornax cluster because it gives the most stringent constraints. Generally speaking, using the NFW profile would give a more conservative result.

Figure 4

The upper limits of the annihilation cross section for 4 channels. The lines in blue and orange represent the limits using NFW and Einasto density profiles respectively (based on the data of the Fornax cluster only) (dashed line: E 0 = 1 GeV; solid line: E 0 = 10 GeV). The black solid lines and dashed lines represent the upper limits of the annihilation cross section based on the Fermi-LAT data of the Milky Way dwarf spheroidal satellite galaxies[6,7]. The black dotted lines represent the thermal relic cross section <συ> = 2.2 × 10−26 cm3 s−1  [36].

For the boost factor, as mentioned above, we used the most conservative values to model the boost factors of Fornax and A2877 clusters. These values depend on the virial mass of the clusters. By using an older value of virial mass for Fornax cluster ()[42], the boost factor of the Fornax cluster is 29, which gives only about 3% larger for the upper limits. We have also tried another extreme empirical relation to calculate the boost factor (based on the simulation results from the Phoenix simulation project[29]), the boost factor of the Fornax cluster is 525. The resultant constraints are much tighter (see Fig. 5).

Figure 5

The upper limits of the annihilation cross section for 4 channels. The lines in blue and orange represent the limits using B = 30 and B = 525 respectively (based on the data of the Fornax cluster only) (dashed line: E 0 = 1 GeV; solid line: E 0 = 10 GeV). The black solid lines and dashed lines represent the upper limits of the annihilation cross section based on the Fermi-LAT data of the Milky Way dwarf spheroidal satellite galaxies[6,7]. The black dotted lines represent the thermal relic cross section <συ> = 2.2 × 10−26 cm3 s−1  [36].

Overall speaking, the assumed parameters and the profiles used are conservative. Therefore, the upper limits obtained are conservative limits.

Discussion and Conclusion

Recently, it has been reported that the excess gamma-ray emission in the Milky Way can be best explained by dark matter annihilation. The proposed dark matter mass is about GeV. On the other hand, standard cosmology indicates that the thermal relic annihilation cross section should be <συ> = 2.2 × 10−26 cm3 s−1 so that it gives the correct amount of dark matter in our universe[36]. In this article, we obtain the upper limits of the annihilation cross sections for four popular channels by using the Fermi-LAT gamma-ray upper limits of two galaxy clusters (A2877 and Fornax). We improve the upper limits of annihilation cross sections by a factor of 1.3–1.8 for wide ranges of m for e + e −, μ + μ − and channels, and a factor of 1.2–1.8 for τ + τ − channel (for m ≤ 100 GeV). Tighter constraints are obtained because we have considered the ‘indirect emission’ (Bremsstrahlung emission) and the effect of substructures in galaxy clusters. In particular, the effect of substructures is significant in galaxy clusters as they are large and massive structures. Besides, we have used the most updated gamma-ray data for analysis, which provide more stringent limits for annihilating dark matter. Furthermore, among a large sample of nearby galaxy clusters, we found that the A2877 and Fornax clusters are the best candidates to constrain dark matter annihilation cross section because they have larger values of I/Φobs.

Besides, the lower limits of the dark matter mass for the thermal relic annihilation cross section are larger than the mass expected in dark matter interpretation of the gamma-ray excess. We show that most of the existing popular dark matter interpretation of the GeV excess in the Milky Way are ruled out. Our results basically support the latest Fermi-LAT analysis of the Milky Way GeV excess[15].

In fact, searching dark matter particles is one of the most important tasks in particle physics. However, the null detection in direct detection experiments may indicate that dark matter mass is probably larger than 1 TeV[43,44]. Therefore, our results are quite consistent with these experimental results that dark matter mass may be much larger than the previous expected values. We expect that future observations of TeV gamma rays[45,46] and anti-proton detection[47,48] might provide some new insight of dark matter annihilation.