Free Access
Issue
A&A
Volume 589, May 2016
Article Number A33
Number of page(s) 13
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201527846
Published online 11 April 2016

© ESO, 2016

1. Introduction

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 (kBT ≃ 1−10 keV) filling the cluster volume, and can 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 elusive1.

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 between re-acceleration2 and hadronic models3. In the re-acceleration model, a seed population of CR electrons can be re-accelerated by interacting turbulent waves, while in the hadronic scenario, the radio-emitting electrons are secondaries produced by CR protons interacting with the protons of the ICM. Currently, the re-acceleration scenario is favoured as the generation mechanism for giant radio halos, while for mini-halos, both the re-acceleration and hadronic models can explain the observed emission (see, e.g., Enßlin et al. 2011; and Brunetti & Jones 2014, for extensive discussions).

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, 2013; Kuchar & Enßlin 2011). 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 Dlum = 78 Mpc (redshift z = 0.018). 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,b). 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 has been intensively observed by the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes. In this paper we present the results obtained with ~250 h 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 Sect. 2, we report our observational results and searches for a CR-induced signal in Sects. 3 and 4, respectively. In Sect. 5, we discuss the interpretation of our constraints on the CR physics in galaxy clusters, and, finally, in Sect. 6 we present our conclusions. Throughout the paper, we assume a standard ΛCDM cosmology with H = 70 km s-1 Mpc-1.

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 h 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 h 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.

thumbnail Fig. 1

Perseus cluster sky maps for an energy threshold of 250 GeV (left-hand panels) and 1 TeV (right-hand panels) obtained from 253 h of MAGIC observation. Upper panels show the relative flux (colour code, expressed in signal-to-background ratio) and the excess significance (contour lines starting from 4σ with steps of 2σ). Lower panels show the significance maps where the signal from IC 310 is subtracted (see text). Symbols indicate the positions of the three brightest radio galaxies of the cluster.

Open with DEXTER

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, θcut, 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 h 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 Fig. 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 Fig. 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 (θ2) with the MAGIC point spread function (PSF) obtained from contemporaneous Crab nebula data at similar zenith angles. In this θ2 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 nebula4 well (MAGIC Collaboration 2016b). The shape of the detected signal is in perfect agreement with a point-like object such as the one expected for an AGN. At both energies, there is no sign of diffuse γ-ray structures inside the cluster.

The average energy spectrum of NGC 1275 obtained with the full data set (August 2009–November 2014) is shown in Fig. 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 power law5(1)with a photon index Γ = 3.6 ± 0.2stat ± 0.2syst and a normalisation constant at 200 GeV of f0 = (2.1 ± 0.2stat ± 0.3syst) × 10-11 cm-2 s-1 TeV-1. 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, Γ ≃ 2, measured by Fermi-LAT in 0.110 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 (Γ = 3.5) as measured in the 0.1−1 TeV range, and for a hard spectral index (Γ = 2.3), 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 Fig. 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 Γ = 3.6), the upper limit of the integral flux above 250 GeV is estimated to be 5.6 × 10-13 cm-2 s-1.

thumbnail Fig. 2

Excess signal above 250 GeV as a function of the squared distance to NGC 1275 (θ2). The solid red line is the MAGIC PSF fit, and the vertical dashed line shows the optimal θ2 cut used for a point-like source detection.

Open with DEXTER

thumbnail Fig. 3

Spectral energy distribution of NGC 1275 averaged over different periods with the newest results in red (~250 h). Arrows indicate the 95%-c.l. differential flux upper limits (5 bins/decade) for a power-law spectrum with a photon index of Γ = 3.5 (thick solid lines) and Γ = 2.3 (thin and dashed lines), respectively.

Open with DEXTER

4. Search for cosmic-ray induced emission

As discussed in Sect. 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 range6. Therefore, we use the energies above 1 TeV to search for the possible diffuse CR-induced component.

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 pα with α ≈ 2.2 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 z = 0.018 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 0.13 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 α = 2.2 would appear above 300 GeV as a power law with an index Γ = 2.33 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.

  • 1.

    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, ~110 TeV, is induced by primary γ rays from the energy range ~30100 TeV. Assuming the most favourable case, when all the absorbed energy is re-emitted, the ICC emission from a power-law spectrum with α = 2.2 would represent about 30% of the intrinsic emission of the ~110 TeV range.

  • 2.

    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 0.1° for an EGMF ≈ 10-14 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 ≫ 10-14 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 ~110 TeV range could be compensated for by the ICC emission, up to a full compensation for a very low EGMF.

  • 3.

    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, 2013; Chang et al. 2012; Pfrommer et al. 2012; Miniati & Elyiv 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, ΓM,knbeam (where nbeam is density of beam pairs), to the IC cooling rate, ΓIC, at a fixed distance to the source. The mean-free-path to pair production, dpp, of primary γ rays with energy >15 TeV is smaller than Dlum = 78 Mpc, the distance of the Perseus cluster. Hence, the density of beam pairs is lower at Dlum in comparison to dpp. As the expected cluster luminosity is smaller than the minimum luminosity of ELE,min ≈ 1042 ergs-1 (at z ≈ 0 and dpp from the source) needed for the oblique instability to grow faster than IC cooling rate, ΓM,k< ΓIC at dpp (a fortiori at Dlum), 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 < 10-14 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 XCR = PCR/Pth (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 profile7. In Fig. 4, we show the γ-ray surface brightness above 100 GeV for the three models. In all cases, fundamental input parameters are the Perseus ICM density and temperature as measured in X-rays (Churazov et al. 2003). The fraction of signal expected within a circular region of a given radius θ from the cluster centre is shown in Fig. 5, together with the MAGIC PSF above 630 GeV.

Since the predicted CR-induced signal extension is significantly larger than the PSF, the optimal θcut used to detect the emission is different than for a point-like source. Comparing the predicted signal in the ON region (θ<θcut) 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 θcut for each model. In practice, θcut 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 θcut 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 Fig. 5.

thumbnail Fig. 4

Surface brightness profiles, above 100 GeV, of the three tested spatial templates for the CR-induced emission in Perseus: isobaric with α = 2.2, semi-analytical, and extended. The normalisation of the isobaric model is set to respect our previous upper limits (MAGIC Collaboration 2012a), while the normalisations of the semi-analytical and extended models are as from Pinzke & Pfrommer (2010) and Zandanel et al. (2014), respectively. See main text for details.

Open with DEXTER

thumbnail Fig. 5

Cumulative fraction of signal within a given radius for different models: point-like in black (long-dashed line), isobaric in green, semi-analytical in red, and extended in blue. The distribution of the point-like model follows the MAGIC PSF above 630 GeV. The thin and thick coloured solid lines represent the real and reconstructed (i.e., smeared by the PSF) signal fractions, respectively. The dotted lines represent the fraction of these signals contained in the reference background region. The vertical dashed lines correspond to the used optimum θcut values.

Open with DEXTER

Table 1

Optimum and contained flux fractions for different models.

Concerning the CR-proton spectrum, we adopt for both the semi-analytical and extended models, the universal power-law momentum spectrum pα, with α ≈ 2.2 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, 2.1 ≤ α ≤ 2.5. In all cases, we also consider the EBL absorption, which results in a slightly softer observed γ-ray spectrum, Γ = α + 0.13, 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 (Eth) and 10 TeV for different spatial models assuming a power-law spectral index Γ = 2.33. The upper limits are converted to the corresponding flux contained within the reference radius 0.15° and the cluster virial radius8R200 ≃ 1.4° (Reiprich & Böhringer 2002). The conversion factors depend on the surface brightness distributions, also considering the signal contained in the OFF region. The latter is particularly important for the extended CR model, which has a flatter surface brightness profile. For the point-like assumption, we estimated the background level from an average of five OFF regions, while for the CR models we only used a single region at 0.8° from the cluster centre, as mentioned in Sect. 2.

Table 2

Integral flux upper limits between Eth and 10 TeV assuming observed power-law spectrum with an index Γ = 2.33.

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 2.1 ≲ Γ ≲ 2.6, which corresponds to an intrinsic index range of about 2.0 ≲ α ≲ 2.5.

5. Interpretation and discussion

The flux upper limits reported in Sect. 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 that 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 ζp,max = 50%. We parametrize ζp,max with a flux multiplier Aγ, which is equal to 1 for ζp,max = 50%, and we assume a linear scaling between Aγ and ζp,max (Fermi-LAT Collaboration 2014). The hadronic-induced emission is proportional to Aγ, which provides the overall normalisation of the CR distribution. In the left-hand panel of Fig. 6, we show the integral γ-ray fluxes predicted by Pinzke & Pfrommer (2010) for the Perseus cluster with Aγ = 1 within 0.15° from the centre and the highest fluxes allowed by our observational upper limits. The value of Aγ is constrained to be ≤ 0.56 when neglecting the EBL absorption and Aγ ≤ 0.75 when including the EBL absorption from Domínguez et al. (2011). This corresponds to ζp,max 28% and 37%, respectively.

thumbnail Fig. 6

Expected integral flux from the Perseus cluster core within 0.15° and the associated 95%-c.l. MAGIC upper limits for the semi-analytical model (Left) and the extended model (Right). For the semi-analytical model, we show the spectrum predicted by Pinzke & Pfrommer (2010) (Aγ = 1, see text) and the one constrained by our upper limits for both cases with (solid lines) and without (dashed lines) EBL absorption. Arrows in light grey show our previous flux upper limits (MAGIC Collaboration 2012a). For both the semi-analytical and extended models, the EBL-absorbed spectra obtained by assuming that the Perseus radio mini-halo has a hadronic origin and fixing a central magnetic field value B0 = 5,10, and 20μG are shown in colours. We used as reference the radio surface brightness radial profile of the Perseus mini-halo at 1.4 GHz from Pedlar et al. (1990). Fluxes are integrated up to 10 TeV.

Open with DEXTER

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 E> 630 GeV that suffers contamination from the NGC 1275 AGN). Our previous constraint, Aγ ≤ 0.8, derived from the upper limit at E> 1 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 Sect. 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 ζp,max ≤ 37%, which is even slightly below our previous results. We note that Fermi-LAT observations suggest values of ζp,max ≲ 15% 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 ζp,max 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 ζp,max rises steeply for shocks with weak Mach numbers M and already saturates at M ≳ 3 (Enßlin et al. 2007). We assess the first assumption by adopting the extended model, shown in the right-hand panel of Fig. 6, which was constructed independently of ζp,max to match the radio emission from the Perseus mini-halo (Zandanel et al. 2014; see Sect. 5.3). Concerning the second assumption, non-linear diffusive shock acceleration models (Kang & Ryu 2011, 2013) predict a slower rise in ζp,max with increasing M, 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.

Table 3

Constraints on the volume-averaged CR-to-thermal pressure ratio within R200 for different CR 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 XCR = PCR/Pth obtained with the isobaric model for different CR spectral indexes, 2.1 ≤ α ≤ 2.5, using the most constraining integral flux upper limit, i.e., the one above 1.6 TeV. The CR-to-thermal pressure must be below about 1% for α ≤ 2.2, below about 2% for α = 2.3, and below 15% for α = 2.5. Including the EBL absorption effectively causes a worsening of the derived constraints of about 25%. 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 0.15° (≈ 0.11 × R200), respectively, as shown in Fig. 7 and Table 3. When considering the full galaxy cluster volume up to R200, XCR is below 2% for the semi-analytical model, but significantly less constrained, 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 XCR 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.

thumbnail Fig. 7

Volume-averaged XCR within a given radius R as constrained by the upper limits presented in this work for the semi-analytical, extended, and isobaric models, the last one with α = 2.2 matching the effective spectral index of the other models in the energies of interest here. We recall that in the isobaric model, XCR(r) is constant by construction. The vertical grey line represents the radius R500, expressed with respect to an average density that is 500 times the critical density of the Universe, for an easier comparison with cosmological studies.

Open with DEXTER

5.3. Radio mini-halo and magnetic fields

As mentioned in Sect. 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 αν = α/ 2 is the synchrotron spectral index (Sνναν), nCR and nICM are the CR and ICM densities, Cγ and Cν are constants that depend on the physics of the hadronic interactions, and εB, εCMB and εSD denote the energy density of the magnetic field (=B2/ 8π), the CMB, and the star-and-dust light in the cluster, respectively. The magnetic field in galaxy clusters can be parametrised as (4)where B0 is the central magnetic field, ne the ICM electron density, and αB 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, εSD, 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 ≤ 0.03 × R200, for low values of the magnetic field, i.e., B0 ≲ 5μG. Therefore, for the case of Perseus, where the observed radio emission arises from within about 0.1 × R200 (see next section), including this term in the synchrotron losses can significantly affect the modelling.

In Fig. 6, we show the EBL-corrected γ-ray emission within 0.15° 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)9. We adopted three different values of the central magnetic field B0 = 5, 10, and 20μG with αB = 0.3 for the semi-analytical model and αB = 0.5 for the extended model as best-fit values. Values of αB ≥ 0.5 are theoretically preferred, for example from simulations of gas sloshing in cool-core clusters (ZuHone et al. 2010, 2011), and the low value of αB required for the semi-analytical model is due to the centrally peaked CR profile in the Pinzke & Pfrommer (2010) simulations, in which AGN feedback was not accounted for, causing enhanced radiative cooling.

At parity of B0, 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 B0 ≳ 5.5μG for the semi-analytical model, for the extended model we obtain B0 ≳ 8μG (for the choices of αB as above). The FR measurements of cluster central magnetic field strengths range from 5μG for merging clusters (Bonafede et al. 2010, 2013) to about 40μ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 25μ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., ϵBϵCMB + ϵSD, 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, Lγ/Lν, becomes independent of the spatial distribution of CRs, of the ICM, and of the magnetic field, if the observed synchrotron spectral index is αν = α/ 2 ≈ 1, as can be seen from Eqs. (2) and (3). Therefore, a minimum theoretical γ-ray flux, , can be derived as (5)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 αν = 1.1 to αν = 1.4 (Sijbring 1993; Gitti et al. 2002). Assuming, for example, B0 = 20μ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 2.1 ≤ α ≤ 2.5, which are consistent with the concavity in the CR spectrum that connects the low-energy CR population around GeV energies (characterised by α ≃ 2.5) with the harder CR population at TeV energies (α ≃ 2.2; 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 Lγ/Lν from the exact equation with α = 2 are within 10%.

We take as reference the total radio luminosity of the mini-halo measured at 327 MHz by Sijbring (1993) of L327 MHz = 17.57 × 10-23 erg s-1 Hz-1 cm-2, and at 1.4 GHz by Pedlar et al. (1990) of L1.4 GHz = 3.04 × 10-23 erg s-1 Hz-1 cm-2. 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.075° (~0.05 × R200) and 0.15° (~0.1 × R200), respectively. The shape of the radio surface brightnesses at 327 MHz and at 1.4 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 θcut 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).

thumbnail Fig. 8

Minimum γ-ray flux in the hadronic scenario obtained assuming that the Perseus radio mini-halo has a secondary origin and that ϵBϵCMB + ϵSD. We show the EBL-corrected fluxes for different spectral indexes α, adopting as reference the total synchrotron luminosity measured both at 327 MHz (Gitti et al. 2002) and at 1.4 GHz (Pedlar et al. 1990). We compare these with our 95%-c.l. point-like upper limits appropriately scaled for the radio surface brightness shape, i.e., multiplied by a factor of 1.15 with respect to what is reported in Table 2. The flux is integrated up to 10 TeV.

Open with DEXTER

Figure 8 shows the EBL-corrected minimum γ-ray emission derived as described above for 2.1 ≤ α ≤ 2.5, and adopting both L327 MHz and L1.4 GHz as synchrotron luminosity. We compare them with the properly scaled (× 1.15) point-like flux upper limits. The most striking result is that the minimum γ-ray flux for α = 2.1 conflicts with our upper limits. This implies that, if α ≤ 2.1, 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.

6. Conclusions

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 0.6° 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 f0(E/ 200 GeV)− Γ with a very steep photon index Γ = 3.6 ± 0.2stat ± 0.2syst and a differential-flux normalisation at 200 GeV of f0 = (2.1 ± 0.2stat ± 0.3syst) × 10-11 cm-2 s-1 TeV-1, 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 5.6 × 10-13 cm-2 s-1.

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 α = 2.2 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 E> 30 TeV are very likely reprocessed via ICC to our energy range of interest. In the most optimistic scenario (EGMF < 10-14 G), the ICC γ rays could fully compensate for the effect of the EBL absorption in 110 TeV range. The strongest constraints on the CR models come from the 1.6–10 TeV integral-flux upper limits of about 5 × 10-14 cm-2 s-1 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, ζp,max, as defined in Pinzke & Pfrommer (2010). The derived constraint, ζp,max ≤ 37%, 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, XCR must be ≲ 1% for α ≤ 2.2, ≲ 2% for α = 2.3, and ≲ 15% for α = 2.5. When considering the semi-analytical model, XCR within R200 is constrained to be ≲ 2%. In the extended model, XCR ⟩ ≲ 2% within 0.15°, but only ≲ 20% within R200 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 B0 ≳ 5.5μG and B0 ≳ 8μ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 (ϵBϵCMB + ϵSD), 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 α ≤ 2.1, the hadronic interpretation of the Perseus radio mini-halo is in conflict with our upper limits. For more realistic α> 2.1, 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 h 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.


4

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 Fig. 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.

5

Power-law fit obtained with the forward-unfolding method over 7 reconstructed-energy bins (χ2/nd.o.f. = 2.4/5).

6

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).

7

The extended model adopted in this work corresponds to the model with γtu = 3 in Zandanel et al. (2014), where γtu is a parameter that indicates the dominant CR transport mechanism.

8

The cluster virial radius is defined here with respect to an average density that is 200 times the critical density of the Universe.

9

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 327 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 327 MHz data was taken.

Acknowledgments

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 13.12.1.3.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. C. Pfrommer gratefully acknowledges support of the Klaus Tschira Foundation.

References

All Tables

Table 1

Optimum and contained flux fractions for different models.

Table 2

Integral flux upper limits between Eth and 10 TeV assuming observed power-law spectrum with an index Γ = 2.33.

Table 3

Constraints on the volume-averaged CR-to-thermal pressure ratio within R200 for different CR models.

All Figures

thumbnail Fig. 1

Perseus cluster sky maps for an energy threshold of 250 GeV (left-hand panels) and 1 TeV (right-hand panels) obtained from 253 h of MAGIC observation. Upper panels show the relative flux (colour code, expressed in signal-to-background ratio) and the excess significance (contour lines starting from 4σ with steps of 2σ). Lower panels show the significance maps where the signal from IC 310 is subtracted (see text). Symbols indicate the positions of the three brightest radio galaxies of the cluster.

Open with DEXTER
In the text
thumbnail Fig. 2

Excess signal above 250 GeV as a function of the squared distance to NGC 1275 (θ2). The solid red line is the MAGIC PSF fit, and the vertical dashed line shows the optimal θ2 cut used for a point-like source detection.

Open with DEXTER
In the text
thumbnail Fig. 3

Spectral energy distribution of NGC 1275 averaged over different periods with the newest results in red (~250 h). Arrows indicate the 95%-c.l. differential flux upper limits (5 bins/decade) for a power-law spectrum with a photon index of Γ = 3.5 (thick solid lines) and Γ = 2.3 (thin and dashed lines), respectively.

Open with DEXTER
In the text
thumbnail Fig. 4

Surface brightness profiles, above 100 GeV, of the three tested spatial templates for the CR-induced emission in Perseus: isobaric with α = 2.2, semi-analytical, and extended. The normalisation of the isobaric model is set to respect our previous upper limits (MAGIC Collaboration 2012a), while the normalisations of the semi-analytical and extended models are as from Pinzke & Pfrommer (2010) and Zandanel et al. (2014), respectively. See main text for details.

Open with DEXTER
In the text
thumbnail Fig. 5

Cumulative fraction of signal within a given radius for different models: point-like in black (long-dashed line), isobaric in green, semi-analytical in red, and extended in blue. The distribution of the point-like model follows the MAGIC PSF above 630 GeV. The thin and thick coloured solid lines represent the real and reconstructed (i.e., smeared by the PSF) signal fractions, respectively. The dotted lines represent the fraction of these signals contained in the reference background region. The vertical dashed lines correspond to the used optimum θcut values.

Open with DEXTER
In the text
thumbnail Fig. 6

Expected integral flux from the Perseus cluster core within 0.15° and the associated 95%-c.l. MAGIC upper limits for the semi-analytical model (Left) and the extended model (Right). For the semi-analytical model, we show the spectrum predicted by Pinzke & Pfrommer (2010) (Aγ = 1, see text) and the one constrained by our upper limits for both cases with (solid lines) and without (dashed lines) EBL absorption. Arrows in light grey show our previous flux upper limits (MAGIC Collaboration 2012a). For both the semi-analytical and extended models, the EBL-absorbed spectra obtained by assuming that the Perseus radio mini-halo has a hadronic origin and fixing a central magnetic field value B0 = 5,10, and 20μG are shown in colours. We used as reference the radio surface brightness radial profile of the Perseus mini-halo at 1.4 GHz from Pedlar et al. (1990). Fluxes are integrated up to 10 TeV.

Open with DEXTER
In the text
thumbnail Fig. 7

Volume-averaged XCR within a given radius R as constrained by the upper limits presented in this work for the semi-analytical, extended, and isobaric models, the last one with α = 2.2 matching the effective spectral index of the other models in the energies of interest here. We recall that in the isobaric model, XCR(r) is constant by construction. The vertical grey line represents the radius R500, expressed with respect to an average density that is 500 times the critical density of the Universe, for an easier comparison with cosmological studies.

Open with DEXTER
In the text
thumbnail Fig. 8

Minimum γ-ray flux in the hadronic scenario obtained assuming that the Perseus radio mini-halo has a secondary origin and that ϵBϵCMB + ϵSD. We show the EBL-corrected fluxes for different spectral indexes α, adopting as reference the total synchrotron luminosity measured both at 327 MHz (Gitti et al. 2002) and at 1.4 GHz (Pedlar et al. 1990). We compare these with our 95%-c.l. point-like upper limits appropriately scaled for the radio surface brightness shape, i.e., multiplied by a factor of 1.15 with respect to what is reported in Table 2. The flux is integrated up to 10 TeV.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.