Deep observation of the NGC 1275 region with MAGIC: search of diffuse -ray emission from cosmic rays in the Perseus cluster
Key Words.:gamma rays: galaxies: clusters – acceleration of particles – galaxies: clusters: individual: Perseus – galaxies: individual: NGC 1275, NGC 1265
Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse -ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about from the centre, have been detected as point-like VHE -ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is described well by a power law with a photon index between 90 GeV and 1200 GeV. We do not detect any diffuse -ray emission from the cluster and so set stringent constraints on its CR population. To bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes . For , the CR-to-thermal pressure within the cluster virial radius is constrained to be %, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to %. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, , that depend on the CR distribution. For , G, which is below the 25 G inferred from Faraday rotation measurements, whereas for , the hadronic interpretation of the diffuse radio emission contrasts with our -ray flux upper limits independently of the magnetic field strength.
Clusters of galaxies represent the latest stage of structure formation,
and are presently assembling through mergers of smaller groups of
galaxies and gas accretion. They are powerful cosmological tools
for testing the evolution of the Universe (see Voit, 2005 for a review).
Cosmic-ray (CR) protons can accumulate in clusters of galaxies for cosmological times,
accelerated by structure formation shocks,
and outflows from active galactic nuclei (AGNs) and galaxies (see, e.g.,
Völk et al., 1996; Berezinsky et al., 1997;
see Brunetti & Jones, 2014 for a review). These CR protons can
interact hadronically with the protons of the intra-cluster
medium (ICM), a hot thermal plasma ( 1–10 keV) filling the cluster volume, and
generate pions. While the charged pions decay
to secondary electrons and neutrinos, the neutral pions decay directly
to high-energy rays.
Despite many observational efforts in the past decade, -ray
emission from clusters of galaxies remains elusive.
Non-thermal emission is observed at radio frequencies in many clusters of galaxies in the form of diffuse synchrotron radiation (see Feretti et al., 2012 for a review). This probes for the presence of relativistic CR electrons and magnetic fields in the cluster environment. However, a conclusive proof of CR-proton acceleration has yet to be found. The observed CR electrons can also produce hard X-rays by inverse-Compton (IC) scattering of cosmic microwave background (CMB) photons. Several claims of IC detection have been made in the past (see Rephaeli et al., 2008 for a review), but more recent observations do not confirm them (Ajello et al., 2009, 2010; Wik et al., 2011, 2012, 2014; Gastaldello et al., 2015), and the possible diffuse IC emission from clusters remains elusive, too.
The observed diffuse radio emission in clusters can be divided in two main
categories: peripheral radio relics and central radio halos
(e.g., Feretti et al., 2012; Brunetti & Jones, 2014). The latter are
usually divided in two other categories: giant-halos hosted in merging
non-cool-core clusters (e.g., the Coma cluster;
Deiss et al., 1997; Brown & Rudnick, 2011), and mini-halos hosted in
relaxed cool-core clusters (e.g., the Perseus cluster;
Pedlar et al., 1990; Sijbring, 1993; Gitti et al., 2002).
While radio relics can be roughly related to merger shocks due to their spatial coincidence and morphology,
the explanation for the origin of radio halos is more challenging.
The generation mechanism of radio halos has been historically debated
The key questions are the following. What is the origin of the radio-emitting electrons? What is the role of CR protons and how do they affect the cluster environment? Upcoming X-ray observations have the potential to detect IC emission in clusters and, hopefully, to break the degeneracy between the electron and the magnetic-field distributions (Kitayama et al., 2014; Bartels et al., 2015), providing an alternative estimate for the magnetic field in clusters with respect to Faraday rotation (FR) measurements (Kim et al., 1991; Clarke et al., 2001; Carilli & Taylor, 2002; Vogt & Enßlin, 2005; Bonafede et al., 2010; Kuchar & Enßlin, 2011; Bonafede et al., 2013). However, the presence and role of CR protons in clusters can only be probed directly through the rays and neutrinos induced by hadronic interactions. The high-energy astronomy window is then crucial for understanding non-thermal phenomena in clusters of galaxies. (This is also true, albeit even more challenging, for neutrinos; see, e.g, Murase et al., 2008; Zandanel et al., 2015.)
The Perseus cluster of galaxies (a.k.a. Abell 426) is a relaxed cool-core cluster located at a distance of about Mpc (redshift ). It hosts the brightest thermal ICM in X-rays (Reiprich & Böhringer, 2002) and a very luminous radio mini-halo (Pedlar et al., 1990; Sijbring, 1993; Gitti et al., 2002). The high ICM density at the centre of the cluster (Churazov et al., 2003) implies a high density of target protons for hadronic interactions with CR protons. Therefore, Perseus is the best cluster for searching for CR-induced -ray emission (we refer the reader to Pinzke & Pfrommer, 2010; MAGIC Collaboration, 2010b, 2012a; Pinzke et al., 2011 for a detailed argumentation). The Perseus cluster also hosts three bright radio galaxies (Ryle & Windram, 1968): NGC 1275, the central dominant galaxy of the cluster, NGC 1265, archetype of a head-tail radio galaxy, in which the jets are bent by their interaction with the ICM, and IC 310, a peculiar object that shows properties of different classifications and which could be an intermediate state between a BL Lac and a radio galaxy. The AGNs of both NGC 1275 and IC 310 show a bright and variable -ray emission in the energy ranges of both the Fermi-Large-Area-Telescope (LAT) (Neronov et al., 2010; Fermi-LAT Collaboration, 2009) and the atmospheric-Cherenkov telescopes (MAGIC Collaboration, 2010a, 2012b, 2014a, 2014b). The NGC 1275 AGN emission obscures the expected diffuse cluster emission over most of the -ray band, in particular in the GeV region, where the spectral energy distribution is peaking (MAGIC Collaboration, 2014a).
Since 2008, the Perseus cluster is intensively observed by the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes. In this paper we present the results obtained with 250 hours of effective observation time taken in stereoscopic mode from 2009 to 2014 and derive tight constraints on the CR population in the cluster. We discuss neither the AGN physics nor the possible indirect detection of dark matter from clusters (Palacio et al., 2015). After describing the observation and data analysis in Section 2, we report our observational results and searches for a CR-induced signal in Sects. 3 and 4, respectively. In Section 5, we discuss the interpretation of our constraints on the CR physics in galaxy clusters, and, finally, in Section 6 we present our conclusions. Throughout the paper, we assume a standard CDM cosmology with km s Mpc.
2 MAGIC observations and data analysis
MAGIC is a system of two 17 m diameter imaging atmospheric Cherenkov telescopes located on the Canary island of La Palma, which observes the -ray sky from 50 GeV to more than 50 TeV. The observation of the Perseus cluster started in 2008 with 24 hr of observation with a single telescope, which did not allow any detection (MAGIC Collaboration, 2010b). Since 2009, MAGIC operates in stereoscopic mode, which provides much better sensitivity (about a factor 2 above 300 GeV). A major upgrade of the telescopes then occurred during the northern-hemisphere summers of 2011 and 2012 (MAGIC Collaboration, 2016a). The improvement in the MAGIC performance after this upgrade is reported by the MAGIC Collaboration (2016b). Here, we combine all the stereoscopic data taken from October 2009 to November 2014. The observations carried out before the upgrade, from October 2009 to February 2011, led to the detection of IC 310 (MAGIC Collaboration, 2010a) and NGC 1275 (MAGIC Collaboration, 2012b). These data were taken solely during dark time at a low zenith angle (from 12 to 36), ensuring the lowest possible energy threshold (100 GeV at the analysis level). After data-quality selection, this first sample consists of 85 hr of effective time. The observations carried out after the upgrade from November 2012 to November 2014 were performed under heterogeneous conditions, including moonlight and large-zenith angles (from 12 to 60), in order to accumulate the largest amount of data possible and improve the performance at TeV energies. After data-quality selection, the second sample consists of 168 h of effective time.
All the observations were taken in so-called “wobble” mode, pointing to alternative sky positions lying 0.4 from NGC 1275, which is at the centre of the cluster. Most of the pointings were also separated by 0.4 from IC 310 to allow the survey of both AGNs simultaneously. The standard MAGIC Analysis and Reconstruction Software (MARS, Zanin et al., 2013) was used to analyse the data. The results of the 2009–2011 campaign were already reported by the MAGIC Collaboration (2012a). Here, we re-analyse these data with same calibration, image cleaning, and gamma-hadron separation as previously and combine them with the more recent data only at the last stage of the analysis. For the 2012–2014 data sample, we use the so-called “sum image cleaning” with cleaning thresholds that are 33% higher than the standard cleaning (reported in MAGIC Collaboration, 2016b) to properly handle the data taken under moonlight conditions. This increases the energy threshold from 100 GeV to about 150 GeV but does not affect the performance at TeV energies. Additionally, only events with more than 150 photo electrons in the camera image of both telescopes are kept. Below this cut our simulation of the instrument response with fixed pixel discriminator thresholds does not correctly describe the actual MAGIC trigger, which uses adaptive thresholds, during Moon observation.
The gamma-hadron separation is performed by the standard MAGIC method using the Random Forest algorithm (MAGIC Collaboration, 2008). The remaining background from the CR-induced air showers is estimated from background control regions (OFF regions) lying at 0.4 from the camera centre. To prevent contamination by the strong IC 310 signal, the OFF regions closer than 0.4 to IC 310 have been excluded. Our analysis assumes different source extensions that only differ on the signal-region radius, , and the number of OFF regions. For the point-like source analysis, an average of 5 OFF regions in the field of view are used. In the case of extended source analysis, only the most distant OFF region, lying 0.8 from NGC 1275, is used.
3 MAGIC results
The 253 hr of stereo observation from different
periods were combined to provide the deepest view of the Perseus cluster at
very high energy. The upper left- and right-hand panels of Figure 1
show the relative-flux (i.e., signal-to-background ratio) sky maps for an energy
threshold of 250 GeV and 1 TeV, respectively.
A clear signal is detected from the two previously discovered AGNs.
The bright and hard source IC 310 is visible in both maps with a high
significance, and it could mask smaller signals. To search for weak emissions,
we modelled the point-like emission from IC 310 as an additional background component. The lower panels of
Figure 1 show the resulting significance sky maps. Above
250 GeV, the emission from NGC 1275 is detected with a significance of
8.5 and a signal-to-background ratio greater than 20%. Above 1 TeV, however, no
source other than IC 310 is detected.
Figure 2 compares the excess (-ray) event
distribution above 250 GeV as a function of the squared distance from NGC 1275
() with the MAGIC point spread function (PSF) obtained from
contemporaneous Crab nebula data at similar zenith angles. In this plot, the PSF is described by a double exponential
(corresponding to a double Gaussian in a 2D map), which matches the excess shape of the Crab nebula
The average energy spectrum of NGC 1275 obtained with the full data set (August
2009 – November 2014) is shown in Figure 3, together with
the previous results from the first two years of observation (MAGIC Collaboration, 2014a).
The new spectrum starts at higher energy because the data include moonlight and
large-zenith-angle observations. With increased photon statistics, we get
better precision and reach higher energies up to 880 GeV.
This last data point is only marginally significant (2)
and in agreement with the upper limits discussed later.
The spectrum between 90 GeV and 1200 GeV can be described well by a simple
with a photon index and a normalisation constant at 200 GeV of . The systematic errors on the flux normalisation and spectral slope take the signal-to-noise ratio into account as explained in MAGIC Collaboration (2016b). Additionally, the uncertainty on the energy scale is 15%.
The average NGC 1275 spectrum from August 2009 to November 2014 is in good agreement with previous measurements (MAGIC Collaboration, 2014a). The transition between a flat spectrum, , measured by Fermi-LAT in 0.1–10 GeV range, and the soft spectrum observed above 100 GeV is confirmed to be smooth, and better fitted by a log parabola or a broken power law compared to a cut-off hypothesis (see discussion in MAGIC Collaboration, 2014a). The AGN-physics interpretation is not substantially modified with respect to our previous results.
Figure 3 also reports the differential flux upper limits of four energy bins above 1 TeV for a point-like source with a soft power-law spectrum () as measured in the TeV range, and for a hard spectral index (), as approximately expected from the CR-induced signal (see next section). Upper limits are calculated using the method of Rolke et al. (2005) for a confidence level (c.l.) of 95% and with a total systematic uncertainty of 30%. These results assume a point-like emission, while the CR-induced signal should be spatially extended. Upper limits on such diffuse emission depend on the assumption of the surface brightness shape. In the next section, we discuss several CR-induced emission models and report the associated flux upper limits.
The non-detection of NGC 1265 in Figure 1 allows us to derive flux upper limits for this radio galaxy. Also, rays could be expected from the central AGN or from the bowshock (Sijbring & de Bruyn, 1998) as speculated earlier for IC 310 before the flux variability was established, confirming the AGN nature of the emission (Neronov et al., 2010). The NGC 1265 position is slightly off-centred, and the exposure at this position is about 20% lower than for NGC 1275. Assuming a point-like source with the same spectral shape as NGC 1275 (a power law with photo index ), the upper limit of the integral flux above 250 GeV is estimated to be cms.
4 Search for cosmic-ray induced emission
As discussed in Section 3, the -ray emission from the central
galaxy NGC 1275 is consistent with a point-like source, and no diffuse component
seems to be present. The measured flux around 100 GeV is much larger than what is
expected from the CR-induced emission, and it could be outshining the latter.
However, the NGC 1275 AGN spectrum is very steep, and no signal is detected
above about 1 TeV. The CR-induced emission should be harder with no
intrinsic cut-off in the MAGIC energy range.
In this section, we discuss the expected energy spectrum and spatial distribution of the CR-induced -ray signal as it should be observed by MAGIC, when taking the absorption during the travel to Earth and the instrument response into account. Then, we present the optimisation of our analysis to detect such emissions and the derived flux upper limits.
4.1 Spectrum expected on Earth
The hydrodynamical simulations of Pinzke & Pfrommer (2010) suggest a universal CR-proton momentum spectrum with at the energies of interest here. The very high-energy -ray spectrum induced by pion decays from CR-ICM interactions should have approximately the same spectral index as the CR spectrum. However, these rays can interact by pair production with the extragalactic background light (EBL) during their travel to Earth, reducing the observed flux. To convert the intrinsic spectrum from the Perseus cluster at into the spectrum observed on Earth, we use the EBL model of Domínguez et al. (2011). Below 300 GeV, the EBL absorption is negligible (¡4%). Between 300 GeV and 10 TeV, the effect of the absorption can be approximated by an increase (softening) in the power-law index of about and a reduction of the differential-flux normalisation at 1 TeV of 17%. Above 10 TeV, the absorption increases dramatically. A source with an intrinsic power-law spectrum index of would appear above 300 GeV as a power law with an index and a cut-off above 10 TeV. About 20% of the flux above 1 TeV and 60% above 10 TeV is absorbed during the travel to Earth.
There are three different scenarios for the fate of the ultra-relativistic electron-positron pairs that are produced by -rayEBL interactions. Each possibility predicts a different signal for the energy range of interest here.
The pairs can interact with CMB photons that are Compton up-scattered in energy. As a result, an initial TeV -ray is reprocessed to GeV energies. If the initial ray has a much higher energy, this process can be repeated to produce the so-called inverse Compton cascade (ICC) effect (see, e.g., Neronov & Vovk, 2010). To first approximation, the ICC emission into the energy range of interest here, 1–10 TeV, is induced by primary rays from the energy range 30–100 TeV. Assuming the most favourable case, when all the absorbed energy is re-emitted, the ICC emission from a power-law spectrum with would represent about 30% of the intrinsic emission of the 1–10 TeV range.
The compact ICC signal could be diluted by the presence of extragalactic magnetic field (EGMF), which deflects the electron and positron tracks and, consequently, induces angular dispersion, as well as a time delay of the signal (see, e.g., Neronov & Vovk, 2010; Tavecchio et al., 2010, 2011; Dermer et al., 2011; Taylor et al., 2011; Dolag et al., 2011; Takahashi et al., 2012). Being constant, the CR-induced emission is not affected by the time delay. According to the approximation of Neronov & Vovk (2010), the ICC emission induced by 50 TeV rays from 78 Mpc should have a spatial extension of the order of for an EGMF G. For a weaker EGMF, the spatial distribution of the ICC would be almost the same as the intrinsic CR-induced emission, and one could expect up to 30% higher signal in the 1–10 TeV range (assuming the intrinsic power-law spectrum extends to 100 TeV). For an EGMF G, however, the ICC would be diluted in the diffuse extragalactic -ray background. Thus, depending on the EGMF level, the -ray loss due to the EBL absorption in the 1–10 TeV range could be compensated for by the ICC emission, up to a full compensation for a very low EGMF.
A competing mechanism exists that could modify the evolution of the pairs on a faster time scale than the ICC, namely powerful plasma instabilities driven by the anisotropy of the ultra-relativistic pair beams (Broderick et al., 2012; Schlickeiser et al., 2012a, b; Chang et al., 2012; Pfrommer et al., 2012; Miniati & Elyiv, 2013; Schlickeiser et al., 2013; Sironi & Giannios, 2014; Supsar & Schlickeiser, 2014; Chang et al., 2014). This picture is interesting because it can match the observed extragalactic -ray background spectrum above 3 GeV and flux distributions of TeV blazars simultaneously, using a unified model of AGN evolution (Broderick et al., 2014a, b). In contrast to the set-up studied by Broderick et al. (2012), here, we have to compare the oblique instability growth rate, (where is density of beam pairs), to the IC cooling rate, , at a fixed distance to the source. The mean-free-path to pair production, , of primary rays with energy 15 TeV is smaller than Mpc, the distance of the Perseus cluster. Hence, the density of beam pairs is lower at in comparison to . As the expected cluster luminosity is smaller than the minimum luminosity of (at and from the source) needed for the oblique instability to grow faster than IC cooling rate, at (a fortiori at ), the absorbed 30–100 TeV -ray flux is very likely reprocessed via ICCs to our energy range of interest.
The EBL effect was neglected in all previous papers on this topic, for which, therefore, the constraints derived from -ray flux upper limits above 1 TeV can be considered as optimistic. They correspond roughly to the case of a strong ICC emission (EGMF G), which compensates for the EBL absorption. Here, in contrast, we include EBL absorption in our reference case and assume the most conservative case without ICC emission. Also, CR propagation effects, such as streaming and diffusion, could cause a softening of the CR spectrum. In the following we do not consider this possibility because of the rather uncertain modelling at this stage (e.g., Wiener et al., 2013).
4.2 Flux upper limits for different spatial distributions
To account for our limited knowledge of the spatial shape of the CR-induced emission,
we adopted three different models that describe the CR density as a function of the distance from the cluster centre:
i) the isobaric model assuming a constant CR-to-thermal pressure (Pfrommer & Enßlin, 2004a); ii) the
semi-analytical model of Pinzke & Pfrommer (2010) derived from hydrodynamical
simulations of clusters; and iii) the extended
hadronic model of Zandanel et al. (2014), where the possibility of CR
propagation out of the cluster core is considered,
resulting in a significantly flatter CR profile.
Since the predicted CR-induced signal extension is significantly larger than the PSF, the optimal used to detect the emission is different than for a point-like source. Comparing the predicted signal in the ON region () with the corresponding background level estimated from our data, and taking the expected signal leakage inside the OFF region at 0.8 from the cluster centre into account, we optimised the for each model. In practice, is optimal on a relatively wide range. The same cut can be used for both the isobaric and extended models, which have very similar optimum values. Since the MAGIC PSF is relatively stable above 630 GeV (MAGIC Collaboration, 2016b), we use the same value for all energies for any given model. The resulting optimum cuts are presented in Table 1 and shown in Figure 5.
|[deg]||ON (OFF)||real (recon.)|
|point-like||0.0075||77% (0.0%)||100% (96%)|
|isobaric||0.0325||41% (0.7%)||39% (35%)|
|semi-analytical||0.0175||46% (0.2%)||55% (50%)|
|extended||0.0325||26% (1.4%)||22% (20%)|
Concerning the CR-proton spectrum, we adopt for both the semi-analytical and extended models, the universal power-law momentum spectrum , with for the energies of interest here, found in hydrodynamical simulations by Pinzke & Pfrommer (2010). For the less predictive isobaric model, the spectral index is free to vary, and we assume a range of values, . In all cases, we also consider the EBL absorption, which results in a slightly softer observed -ray spectrum, , in the energy range 300 GeV–10 TeV.
Table 2 presents the 95%-c.l. upper limits of the integral flux
between several energy thresholds () and 10 TeV for different spatial models
assuming a power-law spectral index .
The upper limits are converted to the corresponding flux contained within the reference radius and
the cluster virial radius
The effective area of MAGIC is relatively flat above 630 GeV, and the integral flux upper limits do not depend strongly on the assumed spectral shape. Therefore, the upper limits reported in Table 2 are valid, within 2%, for an observed spectral index range of , which corresponds to an intrinsic index range of about .
5 Interpretation and discussion
The flux upper limits reported in Section 4 allow to constrain the CR content in the Perseus cluster. We discuss the implications for the CR-to-thermal pressure for the three adopted models: isobaric, semi-analytical, and extended. Additionally, for the semi-analytical model of Pinzke & Pfrommer (2010), our constraints can be translated into constraints on the CR acceleration efficiency at structure formation shocks. Finally, assuming that the Perseus radio mini-halo is induced by secondary electrons (i.e., assuming that the hadronic model for the radio mini-halo is valid), we discuss the expected minimum -ray flux and derive constraints on the cluster magnetic field.
5.1 CR acceleration efficiency
A major uncertainty in modelling CR physics in clusters of galaxies is the CR-acceleration efficiency, i.e., the percentage of the energy dissipated in structure formation shocks, which goes into particle acceleration. We define the CR-acceleration efficiency as in Pinzke & Pfrommer (2010)as the ratio of the energy density in freshly injected CRs to the total dissipated energy density in the downstream region of a shock. We emphasise that it is particularly difficult to put meaningful constraints on this fundamental quantity because of the lack of a general theory of particle acceleration. In fact, our current knowledge of the dependence of the acceleration efficiency with the shock Mach number is mainly phenomenological and numerical. Nevertheless, we assess the impact of our upper limits in the context of the state-of-the-art model by Pinzke & Pfrommer (2010) and caution that our results are also subject to the limitations stated in that work.
Our semi-analytical model is based on predictions by Pinzke & Pfrommer (2010), which are derived from hydrodynamical simulations of galaxy clusters assuming a maximum CR-proton acceleration efficiency %. We parametrize with a flux multiplier , which is equal to for %, and we assume a linear scaling between and (Fermi-LAT Collaboration, 2014). The hadronic-induced emission is proportional to , which provides the overall normalisation of the CR distribution. In the left-hand panel of Figure 6, we show the integral -ray fluxes predicted by Pinzke & Pfrommer (2010) for the Perseus cluster with within from the centre and the highest fluxes allowed by our observational upper limits. The value of is constrained to be when neglecting the EBL absorption and when including the EBL absorption from Domínguez et al. (2011). This corresponds to 28% and 37%, respectively.
Compared to our previous paper (MAGIC Collaboration, 2012a), we accumulated three times more data and derived upper limits that are significantly lower (except for the upper limit at GeV that suffers contamination from the NGC 1275 AGN). Our previous constraint, , derived from the upper limit at TeV, was underestimated because both the EBL absorption and the signal leakage in the OFF regions were neglected. These effects relax the constraint by about 25%. As discussed in Section 4.1, this loss can be compensated by ICC emission. Moreover, the uncertainty on the EBL and ICC effects stays within the 30% systematics considered in deriving the upper limits. Here, considering the most conservative case, we can constrain %, which is even slightly below our previous results. We note that Fermi-LAT observations suggest values of % for the case of the merging Coma cluster (Fermi-LAT Collaboration, 2014; Zandanel & Ando, 2014), although all these results assume a negligible active CR transport, hence represent model-dependent upper limits.
We are able to constrain the CR acceleration efficiency only in the context of the hydrodynamical simulations of Pinzke & Pfrommer (2010) because a full modelling of the formation history of a galaxy cluster and of the CR acceleration at the corresponding structure formation shocks is needed for this. Therefore, our constraint on is strictly valid only in this context, which assumes in particular i) no CR transport relative to the plasma rest frame, and ii) a simplified model for CR acceleration in which rises steeply for shocks with weak Mach numbers and already saturates at (Enßlin et al., 2007). We assess the first assumption by adopting the extended model, shown in the right-hand panel of Figure 6, which was constructed independently of to match the radio emission from the Perseus mini-halo (Zandanel et al., 2014; see Section 5.3). Concerning the second assumption, non-linear diffusive shock acceleration models (Kang & Ryu, 2011, 2013) predict a slower rise in with increasing , which then saturates at higher Mach numbers with respect to Enßlin et al. (2007). However, Cherenkov telescopes probe the high-momentum part of the CR-proton spectrum, which is generated by intermediate-strength shocks where the acceleration efficiency is close to saturation for all models.
Our result on the maximum CR acceleration efficiency is a useful proxy for calibrating the total expected emission in the context of the adopted semi-analytical model, but as stressed above, it should be taken with a grain of salt as general constraint in the context of CR acceleration at shocks. Additionally, the latest supernova remnant observations tend to suggest values below 27% (Helder et al., 2013), with a commonly accepted value around 10% to be able to explain Galactic CRs (e.g., Morlino & Caprioli, 2012). The acceleration conditions at structure formation shock may be different from those at supernova remnants, complicating any direct comparison.
5.2 CR-to-thermal pressure
Depending on the spectral and spatial distribution of the CR protons, the MAGIC upper limits from Table 2 can be translated into constraints on the cluster CR population. In Table 3, we report the constraints on the CR-to-thermal pressure obtained with the isobaric model for different CR spectral indexes, , using the most constraining integral flux upper limit, i.e., the one above TeV. The CR-to-thermal pressure must be below about 1% for , below about 2% for , and below 15% for . Including the EBL absorption effectively causes a worsening of the derived constraints of about %. Therefore, the improvement of the constraints is modest, but the results are more robust with respect to MAGIC Collaboration (2012a).
The constraints on the volume-averaged CR-to-thermal pressure profiles for the semi-analytical and extended models are below 1.4% and 1.9% within (), respectively, as shown in Figure 7 and Table 3. When considering the full galaxy cluster volume up to , is below 2% for the semi-analytical model, but significantly less constrained since below 19% for the extended model. This last weak constraint is expected because the CR distribution is significantly flatter in the extended model, leading to high values in the cluster outskirts characterised by lower thermal pressure.
While in both simulations and analytical models, the CR distribution in galaxy clusters is assumed to roughly scale with the thermal gas, the real CR distribution is unknown and could be significantly flatter if CRs propagate out of the cluster core (Enßlin et al., 2011; Wiener et al., 2013; Zandanel et al., 2014). Indeed, phenomenological evidence from observations of giant radio halos indicate that the CR distribution appears to be flatter than the ICM distribution independently of the generation mechanism of the observed radio emission (Brunetti et al., 2012; Zandanel et al., 2014; Pinzke et al., 2015). This should be kept in mind when using galaxy clusters for cosmological purposes because, depending on the exact amount of CR protons, this can induce a bias in the estimation of cluster masses. While this effect is limited to a few percent for the standard assumptions (isobaric and semi-analytical models), a flat CR distribution could generate up to a 20% bias (e.g., in the specific case considered here) in hydrostatic mass estimates, which are potentially relevant in the current era of precision cosmology.
Notes. Upper limits on the volume-averaged obtained within derived from the integral-flux upper limit in the energy range TeV. is the constraint obtained neglecting the EBL absorption, given for comparison with previous results. includes the EBL-absorption correction.
5.3 Radio mini-halo and magnetic fields
As mentioned in Section 1, the Perseus cluster hosts the brightest radio mini-halo known to date (Pedlar et al., 1990; Sijbring, 1993; Gitti et al., 2002). Assuming that the observed radio emission has a hadronic origin, i.e., is generated by secondary electrons produced in hadronic interactions between the CR and ICM protons, the pion-decay -ray emission is directly linked to the radio signal. Since the intensity of the synchrotron radio emission depends on the amount of secondary electrons, proportional to the hadronically-induced rays, and the local magnetic field, our -ray upper limits can be turned into lower limits on the cluster magnetic field (MAGIC Collaboration, 2012a).
The -ray and synchrotron luminosities can be expressed as (adapted from Pfrommer, 2008)
where is the synchrotron spectral index (), and are the CR and ICM densities, and are constants that depend on the physics of the hadronic interactions, and , and denote the energy density of the magnetic field (=), the CMB, and the star-and-dust light in the cluster, respectively. The magnetic field in galaxy clusters can be parametrised as
where is the central magnetic field, the ICM electron density, and a parameter describing the radial decline of the magnetic field. Such a parametrisation is favoured both by FR measurements in clusters and by hydrodynamical simulations (e.g., Dubois & Teyssier, 2008; Bonafede et al., 2010; Kuchar & Enßlin, 2011).
Our radio synchrotron modelling includes energy losses due to IC scattering of ambient photons. We consider the CMB, as well as the light from stars and dust (SD) in the cluster, , according to the model of Pinzke et al. (2011), which has the advantage of including the average contribution from the intra-cluster light obtained from a stacking of cluster observations. We stress that the SD energy density dominates both the CMB and magnetic field energy densities in the very centre of the cluster, typically , for low values of the magnetic field, i.e., G. Therefore, for the case of Perseus, where the observed radio emission arises from within about (see next section), including this term in the synchrotron losses can significantly affect the modelling.
In Figure 6, we show the EBL-corrected
-ray emission within from the centre for both the
semi-analytical and the extended models, corresponding to the parameters,
for which secondary electrons reproduce the observed radio surface brightness radial
profile of the Perseus mini-halo at 1.4 GHz
(Pedlar et al., 1990).
At parity of , the extended model always shows a higher -ray emission than the semi-analytical model because of the flatter radial profile. In fact, while our -ray upper limits imply G for the semi-analytical model, for the extended model we obtain G (for the choices of as above). The FR measurements of cluster central magnetic field strengths range from G for merging clusters (Bonafede et al., 2010, 2013) to about G for cool-core objects (Kuchar & Enßlin, 2011). For the Perseus cluster, FR measurements are only available on very small scales, i.e., a few tens of pc, and this suggests magnetic field strengths of about G (Taylor et al., 2006). There are, however, large uncertainties related to this measurement, and we refer the reader to more extensive discussions on this topic in our previous publications on the Perseus cluster (MAGIC Collaboration, 2010b, 2012a). In conclusion, when assuming the radio mini-halo in Perseus has a hadronic origin, the lower limits obtained on the magnetic field strength in the Perseus cluster using our -ray flux upper limits are consistent with FR measurements.
5.4 Minimum -ray fluxes
For clusters that host diffuse radio emission, such as the radio mini-halo in Perseus, we can estimate a theoretical minimum -ray flux in the hadronic scenario, which assumes that the observed radio emission has a secondary origin. The idea is that if the magnetic field is strong enough in all the radio-emitting region, i.e., , a stationary distribution of CR electrons loses all its energy to synchrotron radiation (Pfrommer et al., 2008; Pfrommer, 2008; MAGIC Collaboration, 2010b, 2012a). In this case, the ratio of -raytosynchrotron luminosity, , becomes independent of the spatial distribution of CRs, of the ICM, and of the magnetic field, if the observed synchrotron spectral index is , as can be seen from Eqs. (2) and (3). Therefore, a minimum theoretical -ray flux, , can be derived as
A lower magnetic field value would require a higher secondary CR electron density, hence a higher CR proton density, to reproduce the observed radio synchrotron luminosity and would therefore result in a higher -ray flux. This is why we can consider the above -ray flux as a theoretical minimum in the hadronic scenario.
The measured spectral index of the Perseus radio mini-halo ranges from to (Sijbring, 1993; Gitti et al., 2002). Assuming, for example, G, the CR protons responsible for the GHz-synchrotron-emitting secondary electrons have an energy of 20 GeV, which is about 400 times lower than the energy of the CR protons of 8 TeV responsible for the TeV -ray emission. Here, again, we consider spectral indexes of which are consistent with the concavity in the CR spectrum that connects the low-energy CR population around GeV energies (characterised by ) with the harder CR population at TeV energies () (Pinzke & Pfrommer, 2010; see also the discussion in MAGIC Collaboration, 2012a). While Eq. (5) is only approximately valid in the adopted range of spectral indexes, we estimate that the deviations in from the exact equation with are within 10%.
We take as reference the total radio luminosity of the mini-halo measured at MHz by Sijbring (1993) of erg s Hz cm, and at GHz by Pedlar et al. (1990) of erg s Hz cm. The measured maximum emission radius of the Perseus radio mini-halo ranges from about 100 kpc at 1.4 GHz (Pedlar et al., 1990) to about 200 kpc at 327 MHz (Sijbring, 1993; Gitti et al., 2002), which correspond to (0.05) and (0.1), respectively. The shape of the radio surface brightnesses at MHz and at GHz is more compact than the MAGIC PSF. The minimum -ray flux of the hadronic scenario is therefore close to a point-like source for MAGIC. About 67% of the signal would be within the used for the point-like assumption instead of 77% (see Table 1). An appropriate correction factor of 1.15 must be applied to the point-like upper limits presented in Table 2 to correct for the expected emission not being perfectly point-like. In doing so, we implicitly assumed that no -ray emission is coming from beyond the observed extension in radio frequencies (see MAGIC Collaboration, 2012a for a discussion).
Figure 8 shows the EBL-corrected minimum -ray emission derived as described above for , and adopting both and as synchrotron luminosity. We compare them with the properly scaled () point-like flux upper limits. The most striking result is that the minimum -ray flux for conflicts with our upper limits. This implies that, if , the observed diffuse radio emission in Perseus cannot be uniquely hadronic in origin, independently of the magnetic field strength in the cluster. For the softer spectral indexes considered here, the current -ray upper limits cannot exclude the hadronic origin of the diffuse radio emission in Perseus. The other spectral index cases should be in the reach of future ground-based -ray observations with the Cherenkov Telescope Array (CTA; e.g., Doro et al., 2013) with one order of magnitude better sensitivity.
Clusters of galaxies are expected not only to contain -ray-bright AGNs, but also to host diffuse -ray emission due to neutral pion decays, induced by CR-ICM hadronic interactions. Indeed, CR protons should be accelerated in clusters by structure formation shocks and injected by outflows from galaxies and AGNs, and then should hadronically interact with the ICM protons, generating pions. The most promising galaxy cluster to search for such diffuse -ray emission is the Perseus cluster, which has been intensively observed with the MAGIC telescopes since 2008. These observations resulted in the detection of two -ray-bright AGNs in the central galaxy NGC 1275 and in the peculiar galaxy IC 310, both already reported in previous MAGIC publications. Here, we report the search of diffuse -ray emission, using 253 hr of MAGIC observation in stereoscopic mode, accumulated from 2009 to 2014.
The region of the hard point-like source IC 310, located from the cluster centre, can be easily excluded from our search for extended emission as this latter is expected to be centred on the cluster core. The emission from the NGC 1275 AGN, however, overlies the searched signal region. We derived the most precise NGC 1275 spectrum in the range 90 GeV–1.2 TeV ever done. It is described well by a simple power law with a very steep photon index and a differential-flux normalisation at 200 GeV of , in agreement with previous measurements. No signal is detected above approximately 1 TeV. Since the CR-induced emission is expected to have a harder spectrum, we preferred the high energies for our search. No other point-like emission is detected in the cluster, in particular from the radio galaxy NGC 1265, for which we derived a 95%-c.l. flux upper limit above 250 GeV of cms.
To bracket the uncertainty on the CR spatial and spectral distribution in Perseus, we considered three different models. First, the isobaric model, in which the CR-to-thermal pressure is constant and the CR-spectrum index ranges between 2.1 and 2.5. Second, the semi-analytical model of Pinzke & Pfrommer (2010) was derived from hydrodynamical simulations of clusters, for which the CR spectrum follows a universal spectrum with at the energies of interest here. Finally, the extended hadronic model of Zandanel et al. (2014), in which CRs propagate out of the cluster core and generate a significantly flatter radial distribution with respect to the previous two models. In this last model, the CR spectrum is the same as in the semi-analytical model. In this work we have not considered any softening of the CR-proton spectrum induced by possible CR propagation effects (e.g., Wiener et al., 2013).
We optimised our analysis for the different considered CR models. No diffuse -ray emission or large-scale structures were detected in Perseus. We derived 95%-c.l. integral flux upper limits, in different energy ranges and compared to the signal expected from the models over the same range. For the first time, we included the effect of the EBL absorption, which reduces the -ray flux above 1 TeV coming from Perseus by 20%. We discuss the fate of the produced electron pairs, including ICC and possibility of plasma instabilities driven by the anisotropy of the pair beams. We concluded that the absorbed rays with TeV are very likely reprocessed via ICC to our energy range of interest. In the most optimistic scenario (EGMF G), the ICC rays could fully compensate for the effect of the EBL absorption in 1–10 TeV range. The strongest constraints on the CR models come from the 1.6–10 TeV integral-flux upper limits of about in a central region of 0.15 radius.
The comparison with the semi-analytical model sets a constraint on the maximum CR-proton acceleration efficiency, , as defined in Pinzke & Pfrommer (2010). The derived constraint, %, is not much below our previous result obtained with 85 hr of data because the EBL absorption was not taken into account in that early work. Our new study is, therefore, more conservative and more robust. We stress, however, that this constraint is only valid in the context of the Pinzke & Pfrommer (2010) model.
More model-independent constraints were set on the CR-to-thermal pressure ratio in the cluster. In the context of the isobaric model, must be % for , % for , and % for . When considering the semi-analytical model, within is constrained to be %. In the extended model, % within , but only % within because the volume-averaged pressure ratio builds up to significant values in the cluster outskirts where the ICM pressure drops. The actual CR distribution in clusters is unknown, and if it deviates significantly from the ICM distribution, as for example, in the extended model, it could induce a bias on the estimates of the cluster hydrostatic mass, where its contribution is usually neglected, at a level that is potentially important in the current era of precision cosmology.
The Perseus cool-core cluster hosts the brightest known radio mini-halo. Assuming that this diffuse radio emission is generated by synchrotron radiation of secondary electrons from CR hadronic interactions with the ICM, we can turn our -ray flux upper limits into lower limits on the central magnetic field strength in the cluster. For the first time, we included in our modelling energy losses due to IC scattering of ambient photons from stars and dust in the cluster, in addition to the commonly considered CMB. We found G and G for the semi-analytical and extended models, respectively. These constraints are consistent with FR measurements in clusters. Additionally, assuming that CR electrons lose all their energy by synchrotron emission in the radio emitting region (), the derived -ray flux becomes independent of the CR, magnetic field, and ICM distributions. This represents a theoretical flux lower limit in the hadronic scenario because lower magnetic fields would imply a higher -ray emission. With this approach, we found that for the hadronic interpretation of the Perseus radio mini-halo is in conflict with our upper limits. For more realistic , the minimum -ray flux is several times below our upper limits, hence out of reach with MAGIC.
The large amount of data presented in this work, about 250 hr of observations, implies that it would be difficult to significantly improve upon our constraints with the current generation of Cherenkov telescopes. Therefore, this five-year-long campaign represents one of the reference results, together with the Fermi-LAT observations, for the cluster physics in the -ray energy regime until the planned CTA observatory becomes operational in a few years from now.
Acknowledgements.The MAGIC collaboration would like to thank the Instituto de Astrofísica de Canarias for the excellent working conditions at the Observatorio del Roque de los Muchachos in La Palma. The financial support of the German BMBF and MPG, the Italian INFN and INAF, the Swiss National Fund SNF, the ERDF under the Spanish MINECO (FPA2012-39502), and the Japanese JSPS and MEXT is gratefully acknowledged. This work was also supported by the Centro de Excelencia Severo Ochoa SEV-2012-0234, CPAN CSD2007-00042, and MultiDark CSD2009-00064 projects of the Spanish Consolider-Ingenio 2010 programme, by grant 268740 of the Academy of Finland, by the Croatian Science Foundation (HrZZ) Project 09/176 and the University of Rijeka Project 18.104.22.168.02, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, and by the Polish MNiSzW grant 745/N-HESS-MAGIC/2010/0. F. Zandanel acknowledges the support of the Netherlands Organisation for Scientific Research (NWO) through a Veni grant.
- footnotetext: * Corresponding authors: P. Colin
(firstname.lastname@example.org) and F. Zandanel (email@example.com).
- For space-based cluster observations in the GeV-band, see Reimer et al. (2003); Fermi-LAT Collaboration (2010a, b); Jeltema & Profumo (2011); Han et al. (2012); Ando & Nagai (2012); Huber et al. (2013); Zandanel & Ando (2014); Fermi-LAT Collaboration (2014); Prokhorov & Churazov (2014); Vazza & Brüggen (2014); Griffin et al. (2014); Selig et al. (2015); Vazza et al. (2015); Fermi-LAT Collaboration (2015b, a). For ground-based observations in the energy band above 100 GeV, see Perkins et al. (2006); Perkins (2008); HESS Collaboration (2009a, b); Domainko et al. (2009); Galante et al. (2009); Kiuchi et al. (2009); VERITAS Collaboration (2009); MAGIC Collaboration (2010b, 2012a); VERITAS Collaboration (2012); HESS Collaboration (2012).
- See, e.g., Schlickeiser et al. (1987); Brunetti et al. (2001); Petrosian (2001); Gitti et al. (2002); Ohno et al. (2002); Fujita et al. (2003); Brunetti et al. (2004); Brunetti & Blasi (2005); Cassano & Brunetti (2005); Brunetti & Lazarian (2007, 2011); Brunetti et al. (2012); Donnert et al. (2013); ZuHone et al. (2013); Pinzke et al. (2015); Miniati (2015); Bravi et al. (2016) .
- See, e.g., Dennison (1980); Vestrand (1982); Enßlin et al. (1997); Blasi & Colafrancesco (1999); Dolag & Enßlin (2000); Miniati et al. (2001); Miniati (2003); Pfrommer & Enßlin (2003); Gabici & Blasi (2003); Pfrommer & Enßlin (2004a, b); Blasi et al. (2007); Pfrommer et al. (2008); Pfrommer (2008); Kushnir et al. (2009); Donnert et al. (2010a, b); Keshet & Loeb (2010); Pinzke & Pfrommer (2010); Pinzke et al. (2011); Enßlin et al. (2011); Fujita & Ohira (2012); Zandanel et al. (2014); ZuHone et al. (2015).
- The Crab nebula has a much harder spectrum than NGC 1275, so the PSF for NGC 1275 could be slightly larger. The PSF shown in Figure 2 has been normalised to fit the NGC 1275 signal. It is just illustrative and no quantitative result on the intrinsic extension can be derived.
- Power-law fit obtained with the forward-unfolding method over 7 reconstructed-energy bins ().
- A cut-off is expected at PeV energies owing to the escape of high-energy CR protons that are no longer confined in the cluster volume (see, e.g., Völk et al., 1996; Berezinsky et al., 1997; Pinzke & Pfrommer, 2010).
- The extended model adopted in this work corresponds to the model with in Zandanel et al. (2014), where is a parameter that indicates the dominant CR transport mechanism.
- The flux within 0.15 is used as reference for all models.
- The cluster virial radius is defined here with respect to an average density that is 200 times the critical density of the Universe.
- Number of events in the signal () and background () regions.
- Number of events in the signal () and background () regions.
- Significance of the measured excess in standard deviations.
- 95%-c.l. flux upper limits in units of 10 cm s within a radius of 0.15 and 1.4 from the cluster centre.
- 95%-c.l. flux upper limits in units of 10 cm s within a radius of 0.15 and 1.4 from the cluster centre.
- For negative measured excess, we provide conservative upper limits assuming zero excess.
- We take as reference here the radio surface brightness radial profile from Pedlar et al. (1990) at 1.4 GHz rather than the one at MHz from Gitti et al. (2002), as the latter may be affected by residual point-source contamination as pointed out in Sijbring (1993) where the MHz data was taken.
- Ajello, M., Rebusco, P., Cappelluti, N., et al. 2009, ApJ, 690, 367
- Ajello, M., Rebusco, P., Cappelluti, N., et al. 2010, ApJ, 725, 1688
- Ando, S. & Nagai, D. 2012, JCAP, 7, 17
- Bartels, R., Zandanel, F., & Ando, S. 2015, A&A, 582, A20
- Berezinsky, V. S., Blasi, P., & Ptuskin, V. S. 1997, ApJ, 487, 529
- Blasi, P. & Colafrancesco, S. 1999, Astroparticle Physics, 12, 169
- Blasi, P., Gabici, S., & Brunetti, G. 2007, International Journal of Modern Physics A, 22, 681
- Bonafede, A., Feretti, L., Murgia, M., et al. 2010, A&A, 513, A30+
- Bonafede, A., Vazza, F., Brüggen, M., et al. 2013, MNRAS, 433, 3208
- Bravi, L., Gitti, M., & Brunetti, G. 2016, MNRAS, 455, L41
- Broderick, A. E., Chang, P., & Pfrommer, C. 2012, ApJ, 752, 22
- Broderick, A. E., Pfrommer, C., Puchwein, E., & Chang, P. 2014a, ApJ, 790, 137
- Broderick, A. E., Pfrommer, C., Puchwein, E., Chang, P., & Smith, K. M. 2014b, ApJ, 796, 12
- Brown, S. & Rudnick, L. 2011, MNRAS, 412, 2
- Brunetti, G. & Blasi, P. 2005, MNRAS, 363, 1173
- Brunetti, G., Blasi, P., Cassano, R., & Gabici, S. 2004, MNRAS, 350, 1174
- Brunetti, G., Blasi, P., Reimer, O., et al. 2012, MNRAS, 426, 956
- Brunetti, G. & Jones, T. W. 2014, International Journal of Modern Physics D, 23, 30007
- Brunetti, G. & Lazarian, A. 2007, MNRAS, 378, 245
- Brunetti, G. & Lazarian, A. 2011, MNRAS, 410, 127
- Brunetti, G., Setti, G., Feretti, L., & Giovannini, G. 2001, MNRAS, 320, 365
- Carilli, C. L. & Taylor, G. B. 2002, ARA&A, 40, 319
- Cassano, R. & Brunetti, G. 2005, MNRAS, 357, 1313
- Chang, P., Broderick, A. E., & Pfrommer, C. 2012, ApJ, 752, 23
- Chang, P., Broderick, A. E., Pfrommer, C., et al. 2014, ApJ, 797, 110
- Churazov, E., Forman, W., Jones, C., & Böhringer, H. 2003, ApJ, 590, 225
- Clarke, T. E., Kronberg, P. P., & Böhringer, H. 2001, ApJl, 547, L111
- Deiss, B. M., Reich, W., Lesch, H., & Wielebinski, R. 1997, A&A, 321, 55
- Dennison, B. 1980, ApJl, 239, L93
- Dermer, C. D., Cavadini, M., Razzaque, S., et al. 2011, ApJl, 733, L21
- Dolag, K. & Enßlin, T. A. 2000, A&A, 362, 151
- Dolag, K., Kachelriess, M., Ostapchenko, S., & Tomàs, R. 2011, ApJl, 727, L4
- Domainko, W., Nedbal, D., Hinton, J. A., & Martineau-Huynh, O. 2009, International Journal of Modern Physics D, 18, 1627
- Domínguez, A., Primack, J. R., Rosario, D. J., et al. 2011, MNRAS, 410, 2556
- Donnert, J., Dolag, K., Brunetti, G., & Cassano, R. 2013, MNRAS, 429, 3564
- Donnert, J., Dolag, K., Brunetti, G., Cassano, R., & Bonafede, A. 2010a, MNRAS, 401, 47
- Donnert, J., Dolag, K., Cassano, R., & Brunetti, G. 2010b, MNRAS, 407, 1565
- Doro, M., Conrad, J., Emmanoulopoulos, D., et al. 2013, Astroparticle Physics, 43, 189
- Dubois, Y. & Teyssier, R. 2008, A&A, 482, L13
- Enßlin, T., Pfrommer, C., Miniati, F., & Subramanian, K. 2011, A&A, 527, A99+
- Enßlin, T. A., Biermann, P. L., Kronberg, P. P., & Wu, X.-P. 1997, ApJ, 477, 560
- Enßlin, T. A., Pfrommer, C., Springel, V., & Jubelgas, M. 2007, A&A, 473, 41
- Feretti, L., Giovannini, G., Govoni, F., & Murgia, M. 2012, A&Ar, 20, 54
- Fermi-LAT Collaboration. 2009, ApJ, 699, 31
- Fermi-LAT Collaboration. 2010a, JCAP, 5, 25
- Fermi-LAT Collaboration. 2010b, ApJl, 717, L71
- Fermi-LAT Collaboration. 2014, ApJ, 787, 18
- Fermi-LAT Collaboration. 2015a, ApJ, 812, 159
- Fermi-LAT Collaboration. 2015b, ArXiv e-prints [\eprint[arXiv]1507.08995]
- Fujita, Y. & Ohira, Y. 2012, ApJ, 746, 53
- Fujita, Y., Takizawa, M., & Sarazin, C. L. 2003, ApJ, 584, 190
- Gabici, S. & Blasi, P. 2003, Astroparticle Physics, 19, 679
- Galante, N., Acciari, V. A., Aliu, E., et al. 2009, ArXiv e-prints [\eprint[arXiv]0907.5000]
- Gastaldello, F., Wik, D. R., Molendi, S., et al. 2015, ApJ, 800, 139
- Gitti, M., Brunetti, G., & Setti, G. 2002, A&A, 386, 456
- Griffin, R. D., Dai, X., & Kochanek, C. S. 2014, ApJL, 795, L21
- Han, J., Frenk, C. S., Eke, V. R., et al. 2012, MNRAS, 427, 1651
- Helder, E. A., Vink, J., Bamba, A., et al. 2013, MNRAS, 435, 910
- HESS Collaboration. 2009a, A&A, 502, 437
- HESS Collaboration. 2009b, A&A, 495, 27
- HESS Collaboration. 2012, A&A, 545, A103
- Huber, B., Tchernin, C., Eckert, D., et al. 2013, A&A, 560, A64
- Jeltema, T. E. & Profumo, S. 2011, ApJ, 728, 53
- Kang, H. & Ryu, D. 2011, ApJ, 734, 18
- Kang, H. & Ryu, D. 2013, ApJ, 764, 95
- Keshet, U. & Loeb, A. 2010, ApJ, 722, 737
- Kim, K.-T., Kronberg, P. P., & Tribble, P. C. 1991, ApJ, 379, 80
- Kitayama, T., Bautz, M., Markevitch, M., et al. 2014, ArXiv e-prints [\eprint[arXiv]1412.1176]
- Kiuchi, R., Mori, M., Bicknell, G. V., et al. 2009, ApJ, 704, 240
- Kuchar, P. & Enßlin, T. A. 2011, A&A, 529, A13+
- Kushnir, D., Katz, B., & Waxman, E. 2009, JCAP, 9, 24
- MAGIC Collaboration. 2008, Nuclear Instruments and Methods in Physics Research A, 588, 424
- MAGIC Collaboration. 2010a, ApJl, 723, L207
- MAGIC Collaboration. 2010b, ApJ, 710, 634
- MAGIC Collaboration. 2012a, A&A, 541, A99
- MAGIC Collaboration. 2012b, A&A, 539, L2
- MAGIC Collaboration. 2014a, A&A, 564, A5
- MAGIC Collaboration. 2014b, A&A, 563, A91
- MAGIC Collaboration. 2016a, Astroparticle Physics, 72, 61
- MAGIC Collaboration. 2016b, Astroparticle Physics, 72, 76
- Miniati, F. 2003, MNRAS, 342, 1009
- Miniati, F. 2015, ApJ, 800, 60
- Miniati, F. & Elyiv, A. 2013, ApJ, 770, 54
- Miniati, F., Ryu, D., Kang, H., & Jones, T. W. 2001, ApJ, 559, 59
- Morlino, G. & Caprioli, D. 2012, A&A, 538, A81
- Murase, K., Inoue, S., & Nagataki, S. 2008, ApJl, 689, L105
- Neronov, A., Semikoz, D., & Vovk, I. 2010, A&A, 519, L6
- Neronov, A. & Vovk, I. 2010, Science, 328, 73
- Ohno, H., Takizawa, M., & Shibata, S. 2002, ApJ, 577, 658
- Palacio et al. 2015, ArXiv e-prints [\eprint[arXiv]1509.03974]
- Pedlar, A., Ghataure, H. S., Davies, R. D., et al. 1990, MNRAS, 246, 477
- Perkins, J. S. 2008, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 569–572
- Perkins, J. S., Badran, H. M., Blaylock, G., et al. 2006, ApJ, 644, 148
- Petrosian, V. 2001, ApJ, 557, 560
- Pfrommer, C. 2008, MNRAS, 385, 1242
- Pfrommer, C., Chang, P., & Broderick, A. E. 2012, ApJ, 752, 24
- Pfrommer, C. & Enßlin, T. A. 2003, A&A, 407, L73
- Pfrommer, C. & Enßlin, T. A. 2004a, A&A, 413, 17
- Pfrommer, C. & Enßlin, T. A. 2004b, MNRAS, 352, 76
- Pfrommer, C., Enßlin, T. A., & Springel, V. 2008, MNRAS, 385, 1211
- Pinzke, A., Oh, S. P., & Pfrommer, C. 2015, ArXiv e-prints [\eprint[arXiv]1503.07870]
- Pinzke, A. & Pfrommer, C. 2010, MNRAS, 409, 449
- Pinzke, A., Pfrommer, C., & Bergström, L. 2011, Phys. Rev. D,, 84, 123509
- Prokhorov, D. A. & Churazov, E. M. 2014, A&A, 567, A93
- Reimer, O., Pohl, M., Sreekumar, P., & Mattox, J. R. 2003, ApJ, 588, 155
- Reiprich, T. H. & Böhringer, H. 2002, ApJ, 567, 716
- Rephaeli, Y., Nevalainen, J., Ohashi, T., & Bykov, A. M. 2008, Space Science Reviews, 134, 71
- Rolke, W. A., López, A. M., & Conrad, J. 2005, Nuclear Instruments and Methods in Physics Research A, 551, 493
- Ryle, M. & Windram, M. D. 1968, MNRAS
- Schlickeiser, R., Elyiv, A., Ibscher, D., & Miniati, F. 2012a, ApJ, 758, 101
- Schlickeiser, R., Ibscher, D., & Supsar, M. 2012b, ApJ, 758, 102
- Schlickeiser, R., Krakau, S., & Supsar, M. 2013, ApJ, 777, 49
- Schlickeiser, R., Sievers, A., & Thiemann, H. 1987, A&A, 182, 21
- Selig, M., Vacca, V., Oppermann, N., & Enßlin, T. A. 2015, A&A, 581, A126
- Sijbring, D. & de Bruyn, A. G. 1998, A&A, 331, 901
- Sijbring, L. G. 1993, PhD thesis, Groningen University
- Sironi, L. & Giannios, D. 2014, ApJ, 787, 49
- Supsar, M. & Schlickeiser, R. 2014, ApJ, 783, 96
- Takahashi, K., Mori, M., Ichiki, K., & Inoue, S. 2012, ApJl, 744, L7
- Tavecchio, F., Ghisellini, G., Bonnoli, G., & Foschini, L. 2011, MNRAS, 414, 3566
- Tavecchio, F., Ghisellini, G., Foschini, L., et al. 2010, MNRAS, 406, L70
- Taylor, A. M., Vovk, I., & Neronov, A. 2011, A&A, 529, A144
- Taylor, G. B., Gugliucci, N. E., Fabian, A. C., et al. 2006, MNRAS, 368, 1500
- Vazza, F. & Brüggen, M. 2014, MNRAS, 437, 2291
- Vazza, F., Eckert, D., Brüggen, M., & Huber, B. 2015, MNRAS, 451, 2198
- VERITAS Collaboration. 2009, ApJl, 706, L275
- VERITAS Collaboration. 2012, ApJ, 757, 123
- Vestrand, W. T. 1982, AJ, 87, 1266
- Vogt, C. & Enßlin, T. A. 2005, A&A, 434, 67
- Voit, G. M. 2005, Reviews of Modern Physics, 77, 207
- Völk, H. J., Aharonian, F. A., & Breitschwerdt, D. 1996, Space Science Reviews, 75, 279
- Wiener, J., Oh, S. P., & Guo, F. 2013, MNRAS, 434, 2209
- Wik, D. R., Hornstrup, A., Molendi, S., et al. 2014, ApJ, 792, 48
- Wik, D. R., Sarazin, C. L., Finoguenov, A., et al. 2011, ApJ, 727, 119
- Wik, D. R., Sarazin, C. L., Zhang, Y.-Y., et al. 2012, ApJ, 748, 67
- Zandanel, F. & Ando, S. 2014, MNRAS, 440, 663
- Zandanel, F., Pfrommer, C., & Prada, F. 2014, MNRAS, 438, 124
- Zandanel, F., Tamborra, I., Gabici, S., & Ando, S. 2015, A&A, 578, A32
- Zanin et al. 2013, in International Cosmic Ray Conference, id: 0773
- ZuHone, J. A., Brunetti, G., Giacintucci, S., & Markevitch, M. 2015, ApJ, 801, 146
- ZuHone, J. A., Markevitch, M., Brunetti, G., & Giacintucci, S. 2013, ApJ, 762, 78
- ZuHone, J. A., Markevitch, M., & Johnson, R. E. 2010, ApJ, 717, 908
- ZuHone, J. A., Markevitch, M., & Lee, D. 2011, ApJ, 743, 16