Search

Many bright radio relics in the outskirts of galaxy clusters have low inferred Mach numbers, defying expectations from shock acceleration theory and heliospheric observations that the injection efficiency of relativistic particles plummets at low Mach numbers. With a suite of cosmological simulations, we follow the diffusive shock acceleration as well as radiative and Coulomb cooling of cosmic ray electrons during the assembly of a cluster. We find a substantial population of fossil electrons. When reaccelerated at a shock (through diffusive shock acceleration), they are competitive with direct injection at strong shocks and overwhelmingly dominate by many orders of magnitude at weak shocks, Mach < 3, which are the vast majority at the cluster periphery. Their relative importance depends on cooling physics and is robust to the shock acceleration model used. While the abundance of fossils can vary by a factor of ~10, the typical reaccelerated fossil population has radio brightness in excellent agreement with observations. Fossil electrons with 1 < gamma < 100 (10 < gamma < 10^4) provide the main seeds for reacceleration at strong (weak) shocks; we show that these are well-resolved by our simulation. We construct a simple self-similar analytic model which assumes steady recent injection and cooling. It agrees well with our simulations, allowing rapid estimates and physical insight into the shape of the distribution function. We predict that LOFAR should find many more bright steep-spectrum radio relics, which are inconsistent with direct injection. A failure to take fossil cosmic ray electrons into account will lead to erroneous conclusions about the nature of particle acceleration at weak shocks; they arise from well-understood physical processes and cannot be ignored., Published in MNRAS, 25 pages, 20 figures. Cluster region of interest changed to (0.8-1.0) R_200, new figure added, figures and references updated. Conclusions unchanged

About 1/3 of X-ray-luminous clusters show smooth, Mpc-scale radio emission, known as giant radio haloes. One promising model for radio haloes is Fermi-II acceleration of seed relativistic electrons by compressible turbulence. The origin of these seed electrons has never been fully explored. Here, we integrate the Fokker-Planck equation of the cosmic ray (CR) electron and proton distributions when post-processing cosmological simulations of cluster formation, and confront them with radio surface brightness and spectral data of Coma. For standard assumptions, structure formation shocks lead to a seed electron population which produces too centrally concentrated radio emission. Matching observations requires modifying properties of the CR population (rapid streaming; enhanced CR electron acceleration at shocks) or turbulence (increasing turbulent-to-thermal energy density with radius), but at the expense of fine-tuning. In a parameter study, we find that radio properties are exponentially sensitive to the amplitude of turbulence, which is inconsistent with small scatter in scaling relations. This sensitivity is removed if we relate the acceleration time to the turbulent dissipation time. In this case, turbulence above a threshold value provides a fixed amount of amplification; observations could thus potentially constrain the unknown CR seed population. To obtain sufficient acceleration, the turbulent magneto-hydrodynamics cascade has to terminate by transit time damping on CRs, i.e., thermal particles must be scattered by plasma instabilities. Understanding the small scatter in radio halo scaling relations may provide a rich source of insight on plasma processes in clusters., 18 pages, 7 figures, 1 table. This article supersedes arXiv:1503.07870, and has a significantly enlarged analysis and presentation. To be published in MNRAS, ref. MN-15-2504-MJ.R3

Entering a new era of high-energy gamma-ray experiments, there is an exciting quest for the first detection of gamma-ray emission from clusters of galaxies. To complement these observational efforts, we use high-resolution simulations of a broad sample of galaxy clusters, and follow self-consistent cosmic ray (CR) physics using an improved spectral description. We study CR proton spectra as well as the different contributions of the pion decay and inverse Compton emission to the total flux and present spectral index maps. We find a universal spectrum of the CR component in clusters with surprisingly little scatter across our cluster sample. The spatial CR distribution also shows approximate universality; it depends however on the cluster mass. This enables us to derive a semi-analytic model for both, the distribution of CRs as well as the pion-decay gamma-ray emission that results from hadronic CR interactions with ambient gas protons. In addition, we provide an analytic framework for the inverse Compton emission that is produced by shock-accelerated CR electrons and valid in the full gamma-ray energy range. Combining the complete sample of the brightest X-ray clusters observed by ROSAT with our gamma-ray scaling relations, we identify the brightest clusters for the gamma-ray space telescope Fermi and current imaging air Cherenkov telescopes (MAGIC, HESS, VERITAS). We reproduce the result in Pfrommer (2008), but provide somewhat more conservative predictions for the fluxes in the energy regimes of Fermi and imaging air Cherenkov telescopes when accounting for the bias of `artificial galaxies' in cosmological simulations. We find that it will be challenging to detect cluster gamma-ray emission with Fermi after the second year but this mission has the potential of constraining interesting values of the shock acceleration efficiency after several years of surveying.

Galaxy clusters are the largest and most massive gravitationally bound structures known in the Universe. Cosmic-Ray (CR) hadrons accelerated at structure formation shocks and injected by galaxies, are confined in galaxy clusters where they accumulate for cosmological times. The presence of diffuse synchrotron radio emission in several clusters proves the existence of high-energy electrons, and magnetic fields. However, a direct proof of CR proton acceleration is missing. The presence of CR protons can be probe through the diffuse gamma-ray emission induced by their hadronic interaction with the Intra-Cluster Medium (ICM). The Perseus cluster, a nearby cool-core cluster, has been identified to be among the best candidates to detect such emission. We present here the results of a very deep observation of the Perseus cluster with the MAGIC telescopes, accumulating about 250 hours of data from 2009 to 2014. No evidence of large-scale very-high-energy gamma-ray emission from CR-ICM interactions has been detected. The derived flux upper limits in the TeV regime allow us to put stringent constraints on the physics of cluster CRs, in particular on the CR-to-thermal pressure, the CR acceleration efficiency at formation shocks and the magnetic field of the central cluster region., Proceedings of The 34th International Cosmic Ray Conference (ICRC 2015), 8 pages, 8 figures

We report on the detection of very-high energy (VHE, E>100 GeV) gamma-ray emission from NGC 1275, the central radio galaxy of the Perseus cluster of galaxies. The source has been detected by the MAGIC telescopes with a statistical significance of 6.6 sigma above 100 GeV in 46 hr of stereo observations carried out between August 2010 and February 2011. The measured differential energy spectrum between 70 GeV and 500 GeV can be described by a power law with a steep spectral index of \Gamma=-4.1+/-0.7stat+/-0.3syst, and the average flux above 100 GeV is F_{gamma}=(1.3+/-0.2stat+/-0.3syst) x 10^-11 cm^-2 s^-1. These results, combined with the power-law spectrum measured in the first two years of observations by the Fermi-LAT above 100 MeV, with a spectral index of Gamma ~ -2.1, strongly suggest the presence of a break or cut-off around tens of GeV in the NGC 1275 spectrum. The light curve of the source above 100 GeV does not show hints of variability on a month time scale. Finally, we report on the nondetection in the present data of the radio galaxy IC 310, previously discovered by the Fermi-LAT and MAGIC. The derived flux upper limit F^{U.L.}_{gamma} (>300 GeV)=1.2 x 10^-12 cm^-2 s^-1 is a factor ~ 3 lower than the mean flux measured by MAGIC between October 2009 and February 2010, thus confirming the year time-scale variability of the source at VHE., Comment: Corresponding authors are S. Lombardi, P. Colin, D. Hildebrand, and F. Zandanel. Published in A&A; 4 pages, 4 figures

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 4-9 microG, depending on the rate of decline of the magnetic field strength toward larger radii., Corresponding authors are F. Zandanel, C. Pfrommer, P. Colin, A. Pinzke and S. Lombardi. Accepted for publication in A&A on 06/03/2012

We study the possibility for detecting gamma-ray emission from galaxy clusters. We consider 1) leptophilic models of dark matter (DM) annihilation that include a Sommerfeld enhancement (SFE), 2) different representative benchmark models of supersymmetric DM, and 3) cosmic ray (CR) induced pion decay. Among all clusters/groups of a flux-limited X-ray sample, we predict Virgo, Fornax and M49 to be the brightest DM sources and find a particularly low CR-induced background for Fornax. For a minimum substructure mass given by the DM free-streaming scale, cluster halos maximize the substructure boost for which we find a factor above 1000. Since regions around the virial radius dominate the annihilation flux of substructures, the resulting surface brightness profiles are almost flat. This makes it very challenging to detect this flux with imaging atmospheric Cherenkov telescopes. Assuming cold dark matter with a substructure mass distribution down to an Earth mass and using extended Fermi upper limits, we rule out the leptophilic models in their present form in 28 clusters, and limit the boost from SFE in M49 and Fornax to be < 5. This corresponds to a limit on SFE in the Milky Way of < 3, which is too small to account for the increasing positron fraction with energy as seen by PAMELA and challenges the DM interpretation. Alternatively, if SFE is realized in Nature, this would imply a limiting substructure mass of M_lim > 10^4 M_sol - a problem for structure formation. Using individual cluster observations, it will be challenging for Fermi to constrain our selection of DM benchmark models without SFE. The Fermi upper limits are, however, closing in on our predictions for the CR flux using an analytic model based on cosmological hydrodynamical cluster simulations. We limit the CR-to-thermal pressure in nearby bright galaxy clusters of the Fermi sample to < 10% and in Norma and Coma to < 3%., 43 pages, 23 figures, 10 tables. Accepted for publication in Phys. Rev. D: streamlined paper, added a paragraph about detectability to introduction, few references added, and few typos corrected

We report on the detection of very high energy emission (> 300 GeV) from IC310 with the MAGIC telescopes. IC 310 is a head-tail radio galaxy located in the Perseus galaxy cluster, and has been recently detected by the Fermi-LAT instrument. The spectrum observed by the MAGIC telescopes is hard (spectral index −2.00± 0.14). Strong hints of variability are seen in the data. We discuss the possible origin of the very high energy gamma-ray emission.

Recently, it has been shown that electrons and positrons from dark matter (DM) annihilations provide an excellent fit to the Fermi, PAMELA, and HESS data. Using this DM model, which requires an enhancement of the annihilation cross section over its standard value to match the observations, we show that it immediately implies an observable level of gamma-ray emission for the Fermi telescope from nearby galaxy clusters such as Virgo and Fornax. We show that this DM model implies a peculiar feature from final state radiation that is a distinctive signature of DM. Using the EGRET upper limit on the gamma-ray emission from Virgo, we constrain the minimum mass of substructures within DM halos to be > 5x10^-3 M_sun -- about four orders of magnitudes larger than the expectation for cold dark matter. This limits the cutoff scale in the linear matter power spectrum to k < 35/kpc which can be explained by e.g., warm dark matter. Very near future Fermi observations will strongly constrain the minimum mass to be > 10^3 M_sun: if the true substructure cutoff is much smaller than this, the DM interpretation of the Fermi/PAMELA/HESS data must be wrong. To address the problem of astrophysical foregrounds, we performed high-resolution, cosmological simulations of galaxy clusters that include realistic cosmic ray (CR) physics. We compute the dominating gamma-ray emission signal resulting from hadronic CR interactions and find that it follows a universal spectrum and spatial distribution. If we neglect the anomalous enhancement factor and assume standard values for the cross section and minimum subhalo mass, the same model of DM predicts comparable levels of the gamma-ray emission from DM annihilations and CR interactions. This suggests that spectral subtraction techniques could be applied to detect the annihilation signal., 5 pages, 2 figures (published version; minor corrections to figures and result, equation added)