Cosmic-ray protons accumulate for cosmological times in clusters of galaxies as their typical radiative and diffusive escape times are longer than the Hubble time. Their hadronic interactions with protons of the intra-cluster medium generate secondar
y electrons, gamma-rays and neutrinos. We here estimate the contribution from clusters to the diffuse gamma-ray and neutrino backgrounds. We model the cluster population by means of their mass function, using a phenomenological luminosity-mass relation applied to all clusters, as well as a detailed semi-analytical model. Additionally, we consider observationally-motivated values for the cluster magnetic field. This is a crucial parameter since the observed radio counts due to synchrotron emission by secondary electrons need to be respected. For a choice of parameters respecting all current constraints, and assuming a spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the extragalactic gamma-ray background observed by Fermi and to the IceCube flux. They account for less than 1% for spectral indices <=-2. The neutrino flux observed by IceCube can be reproduced without violating constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with energy. In the case of proton-photon interactions, we find that very likely protons do not reach sufficiently high energies to produce neutrinos in clusters. We argue that our results are optimistic due to our assumptions, and that clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds. Finally, we find that the cluster contribution to the angular fluctuations in the gamma-ray background is subdominant, less than 10%. [abridged]
The underlying physics of giant and mini radio halos in galaxy clusters is still an open question. We find that mini halos (such as in Perseus and Ophiuchus) can be explained by radio-emitting electrons that are generated in hadronic cosmic ray (CR)
interactions with protons of the intracluster medium. By contrast, the hadronic model either fails to explain the extended emission of giant radio halos (as in Coma at low frequencies) or would require a flat CR profile, which can be realized through outward streaming and diffusion of CRs (in Coma and A2163 at 1.4 GHz). We suggest that a second, leptonic component could be responsible for the missing flux in the outer parts of giant halos within a new hybrid scenario and we describe its possible observational consequences. To study the hadronic emission component of the radio halo population statistically, we use a cosmological mock galaxy cluster catalog built from the MultiDark simulation. Because of the properties of CR streaming and the different scalings of the X-ray luminosity (L_X) and the Sunyaev-Zeldovich flux (Y) with gas density, our model can simultaneously reproduce the observed bimodality of radio-loud and radio-quiet clusters at the same L_X as well as the unimodal distribution of radio-halo luminosity versus Y; thereby suggesting a physical solution to this apparent contradiction. We predict radio halo emission down to the mass scale of galaxy groups, which highlights the unique prospects for low-frequency radio surveys (such as the LOFAR Tier 1 survey) to increase the number of detected radio halos by at least an order of magnitude.
Cosmological hydrodynamical simulations of galaxy clusters are still challenged to produce a model for the intracluster medium that matches all aspects of current X-ray and Sunyaev-Zeldovich observations. To facilitate such comparisons with future si
mulations and to enable realistic cluster population studies for modeling e.g., non-thermal emission processes, we construct a phenomenological model for the intracluster medium that is based on a representative sample of observed X-ray clusters. We create a mock galaxy cluster catalog based on the large collisionless N-body simulation MultiDark, by assigning our gas density model to each dark matter cluster halo. Our clusters are classified as cool-core and non cool-core according to a dynamical disturbance parameter. We demonstrate that our gas model matches the various observed Sunyaev-Zeldovich and X-ray scaling relations as well as the X-ray luminosity function, thus enabling to build a reliable mock catalog for present surveys and forecasts for future experiments. In a companion paper, we apply our catalogs to calculate non-thermal radio and gamma-ray emission of galaxy clusters. We make our cosmologically complete multi-frequency mock catalogs for the (non-)thermal cluster emission at different redshifts publicly and freely available online through the MultiDark database (www.multidark.org).
Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to
be luminous gamma-ray emitters since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven galactic winds. To detect the gamma-ray emission from CR interactions with the ambient cluster gas, we conducted the deepest to date observational campaign targeting a galaxy cluster at very high-energy gamma-rays and observed the Perseus cluster with the MAGIC Cherenkov telescopes for a total of ~85 hr of effective observing time. This campaign resulted in the detection of the central radio galaxy NGC 1275 at energies E > 100 GeV with a very steep energy spectrum. Here, we restrict our analysis to energies E > 630 GeV and detect no significant gamma-ray excess. This constrains the average CR-to-thermal pressure ratio to be <= 1-2%, depending on assumptions and the model for CR emission. Comparing these gamma-ray upper limits to predictions from cosmological cluster simulations that include CRs constrains the maximum CR acceleration efficiency at structure formation shocks to be < 50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Finally, we derive lower limits on the magnetic field distribution assuming that the Perseus radio mini-halo is generated by secondary electrons/positrons that are created in hadronic CR interactions: assuming a spectrum of E^-2.2 around TeV energies as implied by cluster simulations, we limit the central magnetic field to be > 4-9 microG, depending on the rate of decline of the magnetic field strength toward larger radii.
The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{
-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies.
