Free Access
Issue
A&A
Volume 595, November 2016
Article Number A116
Number of page(s) 11
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201628788
Published online 09 November 2016

© ESO, 2016

1. Introduction

Radio halos (RHs) are Mpc-scale diffuse synchrotron sources observed in the central regions of galaxy clusters (see e.g. Feretti et al. 2012, for a recent observational review). Radio halos are found in about one-third of massive clusters (see Cuciti et al. 2015, for an updated study) and are located in merging systems (e.g. Cassano et al. 2013; see Bonafede et al. 2014, for a case of radio halo in an apparently relaxed system).

It is generally believed that turbulence induced by mergers in galaxy clusters can reaccelerate the relativistic electrons responsible for the origin of radio halos (see e.g. Brunetti et al. 2001; Petrosian 2001; Brunetti & Jones 2014, for reviews). This scenario naturally explains the connection between radio halos and mergers; however, it also poses fundamental questions on the micro-physics of the mechanisms that are responsible for the acceleration and transport of these relativistic electrons (see Miniati 2015; Brunetti & Lazarian 2016, for recent discussions).

The study of the connection between radio halos and the thermal properties of the hosting clusters sheds light on these mechanisms and in general on the interplay between thermal and non-thermal components in these systems.

Several correlations between thermal and non-thermal properties of galaxy clusters have been found (P1.4LX, P1.4–Mass, P1.4Y500; where P1.4 and Y500 are, respectively, the radio power of halos at 1.4 GHz and the cluster integrated SZ signal within the radius at which the mean mass density is 500 times the critical density at the cluster redshift R500; see e.g. Basu 2012; Cassano et al. 2013). In this respect radio follow up of X-ray or mass-selected samples of galaxy clusters provide a unique way to probe the formation of radio halos and their connection with cluster mergers. Models predict that the bulk of radio halos should be generated at z = 0.2−0.4 (e.g. Cassano et al. 2006), yet current statistical studies are available only for a limited range of masses and redshifts (e.g. Cuciti et al. 2015).

In this framework, this work is part of the MACS-Planck Radio Halo Cluster Project (Macario et al. 2014) conceived to explore the origin and occurrence of RHs and their connection with the dynamical state of the host systems by extending previous studies to a higher redshift range. Feretti et al. (2012) show that the redshift distribution of clusters hosting RHs is homogeneous up to z = 0.35, but statistically incomplete at higher redshifts. Our sample includes 32 intermediate redshift clusters (0.3 < z < 0.45), which are being analysed through deep ~325 MHz GMRT or ~2.1 GHz ATCA observations, depending on the declination of the targets. In this paper we present the discovery of a RH in PLCK G285.0-23.7, which is a massive cluster (M500 = 8.39 × 1014MPlanck Collaboration IX 2011) located at redshift z = 0.39.

In the following, Sect. 2 describes the ATCA observations of PLCK G285.0-23.7 together with the data reduction and the image reconstruction strategy. The physical properties of the detected diffuse radio source and our analysis of its nature are presented in Sect. 3. Our study is discussed and concluded in Sect.  4. In the ΛCDM cosmology adopted throughout this paper (with H0 = 70 km s-1 Mpc-1, ΩM = 0.3, ΩΛ = 0.7), 1 arcsec corresponds to 5.29 kpc at the redshift of PLCK G285.0-23.7.

2. Radio observations and data reduction

Radio observations of PLCK G285.0-23.7 were undertaken on the Australia Telescope Compact Array (ATCA) in three separate array configurations (6D, 750A, 1.5A) using the Compact Array Broadband Backend (CABB) correlator with a central frequency of 2.1 GHz and spanning 1.13.1 GHz. Observations were carried out in continuum mode with the correlator set to produce 2000 × 1 MHz channels. Details of the observations can be found in Table 1. The primary flux scale was set relative to the unresolved source PKS B1934-638 for which the detailed spectral behaviour is well understood over the ATCA band. The phase calibration was initially performed relative to PKS B1036-697 for first of two observations in the 6D configuration which commenced on 8 June 2012. However, it was quickly noted that PKS B1036-697 was partially confused in this configuration and so this calibrator was replaced with PKS B0606-795 for the remaining observation run.

Table 1

Details of our ATCA observations towards PLCK G285.0-23.7 .

Radio frequency interference (RFI) and bad channels were manually excised from the target and primary and secondary calibrators using a combination of clipping algorithms and visual inspection via the PGFLAG task within MIRIAD (Sault et al. 1995). Owing to the nature of CABB data, it is necessary to perform calibration and cleaning on narrower frequency intervals. The precise intervals used depend on properties of the observations (phase stability) and the complexity of sources in the field, which affects the success of cleaning. We conducted a number of trials and determined that ~500 MHz × 4 sub-bands produce the optimal results for these data. Thus, the target, primary, and secondary calibrator data were divided into the required sub-bands centred at 1.381 (from 1.130 GHz to 1.630 GHz), 1.867 (from 1631 GHz to 2.130 GHz), 2.380 (from 2.131 GHz to 2.630 GHz), and 2.769 MHz (from 2.631 GHz to 3.100 GHz). Each sub-band was then calibrated in MIRIAD as per standard procedures outlined in the 2014 MIRIAD manual to produce four continuum images. During this procedure it became apparent that data from the 6D array suffered from a bandpass ripple on all baselines associated with antennas 5 and 6 which persisted for 75 min. This effect is a transient error of unknown origin seen rarely in ATCA data and was not evident in data from the 750A or 1.5A configuration. Attempts to correct for this proved unsuccessful and baselines to antennas 5 and 6 were removed for this configuration over this time. An example of the resultant uv-coverage for the combination of all three configurations for the 1.867 GHz sub-band image is shown in Fig. 1.

thumbnail Fig. 1

uv-coverage of PLCK G285.0-23.7 observations in the 500 MHz band centred at 1.867 MHz. The wide-band of the full 2 GHz observations completes the coverage.

thumbnail Fig. 2

Final full-resolution, wide-band image of PLCK G285.0-23.7. The outer, big circle denotes the boundary of the primary beam with a radius of ~0.22 deg. The central 1 Mpc diameter region is indicated by the inner, smaller red circle. The rms noise of this image is 11.3 μJy/beam and the resolution is 5′′.

Table 2

Properties of the full-resolution (first five rows) and tapered (last five rows) radio maps in different frequency bands centred at ν0.

Following calibration, sub-band images were created out to the second null in the point spread function using the task INVERT and a Steer CLEAN (Steer et al. 1984) was applied to all sources within the primary beam. Two bright sources (at coordinates RA = 07:19:25.285, Dec = 73:32:21.16 and RA = 07:23:30.799, Dec = 73:08:40.16) with fluxes of just over 100 mJy sitting just outside of the primary beam were difficult to clean. They were modelled and subtracted from all of the visibility data and the sub-bands were re-imaged. Details of the resultant sub-band images are presented in Table 2. The first 200 MHz of the band is strongly affected by RFI meaning the image at 1.3 GHz has the lowest sensitivity.

The sub-band images were convolved with a Gaussian to a common resolution, slightly lower than the lowest resolution sub-band image. The sub-band images were then added together to create a final wide-band image (see Table 2). As the noise levels in the sub-band images are mostly identical for the three highest bands and only differed marginally in the lowest frequency band, it was not necessary to weight the images relative to the sensitivities in the mosaicing process. The final deep ATCA image achieves a root mean squared (rms) noise of ~11 μJy/beam measured with AIPS TVSTAT at the field centre in regions inside the primary beam without any trace of point sources or diffuse emission. We get a dynamic range of 10 000:1, making it one of the deepest, highest dynamic range ATCA images yet published, equal to the recent image of the Bullet Cluster at 11 μJy/beam (Shimwell et al. 2014, 2015). The image is shown in Fig. 2.

A visual inspection of our full-resolution ATCA image of PLCK G285.0-23.7 shows no obvious presence of diffuse radio emission in the central area where the cluster lies (see Figs. 2 and 3 for a zoomed version on the cluster). In order to investigate the possible presence of low-surface-brightness diffuse sources in the cluster, visibilities were re-imaged with a robust = 0.5 weighting and a ~9 kλ Gaussian tapering was applied allowing a FWHM of 30 arcsec (see Table 2). Again, data were imaged in each of the four sub-bands with the taper applied and then mosaicked into a single deep map (see Table 2). A diffuse, low-surface-brightness source is detected coincident with the cluster core.

thumbnail Fig. 3

Contours of the ATCA tapered image centred at 1.867 GHz are overlaid on the ATCA wide-band tapered (left) and full-resolution maps (right). Dashed contours correspond to the islands of significant emission detected by PyBDSM on the Block 3 tapered map (cyan contours in Fig. A.1). The continuous curves trace instead “classical” 3, 3 ×2\hbox{$\times \sqrt{2}$}, 6σ contours (i.e. with the 51.1 μJy/beam rms value of the map estimated by hand on the Block 3 tapered map at ~30 arcsec resolution; see Table 2).

3. Results

3.1. Detection and characterisation of a radio halo in PLCK G285.0-23.7

Contours of the low-resolution radio map of PLCK G285.0-23.7 are overlaid on the X-ray XMM image of the same sky region in Fig. 4. The morphology of the diffuse radio source matches very closely the thermal ICM X-ray emission in terms of extension and of elongation (with a main axis along the east–west direction). This, together with the disturbed dynamical state of the cluster suggested by its X-ray morphology and density profile (Planck Collaboration IX 2011), allows us to confidently classify the diffuse source in PLCK G285.0-23.7 as a classical RH.

thumbnail Fig. 4

Radio contours of the ATCA tapered image centred at 1.867 GHz overlaid on the smoothed raw XMM image in the [0.32.0] keV energy band of PLCK G285.0-23.7 (Planck Collaboration IX 2011). Contour levels are the same as in Fig. 3.

Table 3

Physical properties of the diffuse radio source at 1.867 GHz from tapered images.

We now examine the radio flux density of the diffuse source on the tapered image relative to Block 3 (ID 2 in Table 3), i.e. the image centred at 1.867 GHz and characterised by an average rms noise level of 51.1 μJy/beam (Table 2), as it shows the best sensitivity to the diffuse radio emission. Owing to the typical spectral behaviour of radio halos (with synchrotron spectral index α ≳ 1)1, the rms sensitivity of our maps is comparable from Block 1 to Block 3; however, a lower luminosity is expected for the diffuse radio source at the higher frequencies sampled by Blocks 1 and 2. In terms of sensitivity to steep synchrotron sources, Block 4 would have been the optimal frequency range in which to measure the radio halo flux density, if it had not been badly affected by RFI (~60% of the data were flagged).

thumbnail Fig. 5

Zoom in the HST-ACS image of PLCK G285.0-23.7 central field (available in the HST archive). For reference, the 3σ contour of the ATCA tapered image centred at 1.867 GHz is shown in black. The blue contours indicate sources of significant radio emission in the full-resolution deep ATCA image of the cluster (i.e. 5 × rms of the final wide-band ATCA map shown in Fig. 2).

The total flux density was first measured on the Block 3 tapered image by integration over the source image within 3σ contours using the AIPS task TVSTAT. By overlaying the 5σ contours corresponding to the full-resolution wide-band radio map on the optical image of the cluster central region, we identify two point-like significant objects lying within the region occupied by the diffuse source on the low-resolution radio map. The two compact radio sources show quite clear optical counterparts (see Fig. 5). We measured their flux densities using the AIPS verb JMFIT and we subtracted them from the total flux density measured inside the 3σ contours of the tapered image.

To verify that large-scale diffuse emission is not caused by the blending of discrete sources and to get a complementary flux measurement, we then produced an image of the diffuse cluster emission by subtracting compact sources in the uv-plane. We first identified the clean components of the discrete sources by using only the longest baselines (uv-range 3.640.0 kλ). These components were then subtracted from the original data set in the uv-plane by using the MIRIAD task UVMODEL. At this point, we calculated the flux density of the residual diffuse emission by integrating the surface brightness down to the same region considered before, i.e. within the 3σ level of the non-point source subtracted tapered map (ID 1 in Table 3). The results of these flux density measurements are given in the top part of Table 3. The error in flux density is calculated following the formula ΔF=(σrms)2Nbeam+0.01F2\hbox{$\Delta F = \sqrt{(\sigma_{\rm rms})^2N_{\rm beam}+0.01F^{2}}$}, where F is the measured flux density and Nbeam is the number of beams contained in the measured area.

Table 3 shows the results of the different flux density measurements; all the values are consistent within the error bars. The flux measurement that we adopt in our following analysis is the one measured in the tapered map with point source subtracted from the uv-plane (ID1 in Table 3), for which we obtain a flux density of 2.02 ± 0.25 mJy at 1.867 GHz.

In order to get an estimate of the size of the diffuse radio source, we measured the smallest ellipse that fully contains the 3σ contours of the Block 3 tapered map without point source subtraction (ID2 in Table 3). We obtain an angular extent of 2.34′ × 1.34′, which, for our cosmology, corresponds to a major and minor axis of ~742 kpc and 425 kpc, respectively (see Table 3).

thumbnail Fig. 6

Synchrotron power of RHs at 1.4 GHz (P1.4 GHz) vs. cluster mass (M500 from Planck Collaboration XXVI 2016). Symbols are colourcoded depending on the redshift of RH hosting clusters. The red line is the fit derived in the present paper, including only the RHs with radio fluxes measured at 1.4 GHz (indicated by triangles; see also Table B.1). Analogously to Cassano et al. (2013), we exclude for the fitting ultra-steep spectrum RHs and halos detected in frequencies different from 1.4 GHz, both shown in dots (see Table B.1). The red star corresponds to PLCK G285.0-23.7 detected in this work. When the spectral index is unknown (including PLCK G285.0-23.7 ) we assume the value of α = 1.3.

Finally, we measured with TVSTAT the total flux density on the ID1 and ID2 maps by following the contours of the significant islands of emission identified by the source finder software PyBDSM around the diffuse source (dashed contours in Fig. 3; see Appendix A for details on the use of PyBDSM). The results obtained following this final method are given in the bottom part of Table 3.

3.2. Radio luminosity vs. cluster mass correlation

We calculated the power of the newly discovered diffuse radio source and compared it with radio powers of previously discovered halos. We selected objects firmly classified as giant radio halos (i.e. with diffuse emission extending beyond the cluster core), whose flux density has been measured through radio interferometric observations and with a point source identification and subtraction strategy very similar to ours, and with information about M500 in the PSZ2 cluster catalog (Planck Collaboration XXVI 2016). The list of sources is presented in Table B.1.

In order to convert all values to the cosmology adopted in this paper and to take into account some inconsistencies found for published radio powers, we compiled the total flux densities of RHs reported in the literature (see Table B.1). We then calculated radio powers using the formula Lν=4πDL2Sν(1+z)α1\hbox{$L_{\nu}=4\pi D_{L}^{2}S_{\nu}(1+z)^{\alpha-1}$}, where DL is the luminosity distance of the cluster. Whenever a measured value of the spectral index is not available, we adopted α = 1.3.

Traditionally RH powers are reported at 1.4 GHz; however, owing to the lower quality of the Block 4 image centred at 1.381 GHz, we considered the most reliable sub-band map (centred at 1.867 GHz) and then extrapolated to obtain the 1.4 GHz radio power. As the spectral index of the radio source is too uncertain over the ATCA band owing to the very low surface brightness of the object, we undertake this extrapolation assuming the value of α = 1.3. We obtain a value for the radio power at 1.4 GHz of 1.72 ± 0.22 × 1024 W/Hz.

In Fig. 6 we plot the radio power vs. cluster mass for the RH in PLCK G285.0-23.7 (shown as a red star) as compared to all halos included in Table B.1, which are indicated by triangles (colourcoded based on the redshift of their host clusters), except when they are classified as ultra steep spectrum (USS) radio halos and/or when their flux was not measured at 1.4 GHz (see Table B.1), in which case they appear as dots (keeping the same colour code). Analogously to Cassano et al. (2013), we fit a power-law relation using linear regression in the log–log space and adopt a BCES-bisector regression algorithm (Akritas & Bershady 1996). For the fitting, shown as a red line in Fig. 6 and characterised by a slope of 4.24 ± 0.14, we consider only the RHs represented as triangles in Fig. 6. The best fit derived by Cassano et al. (2013), traditionally reported in the literature, has a slope of 3.77 ± 0.57. For completeness we also derived the value for the slope using all radio powers of halos present in Table B.1 for the best fit, which is 4.05 ± 0.09. The three results are consistent within the error bars. The RH in PLCK G285.0-23.7 is slightly under-luminous with respect to the best-fit relation, but is within the scatter of the observed data points. We also point out that, based on the present radio data, we do not find a clear trend with redshift of the P1.4M500 correlation.

Following Cassano et al. (2007), we show in Fig. 7 the position of the newly detected radio halo on the plot of radio halo sizes as a function of the virial radius of their host clusters. There is evidence that powerful radio halos are also bigger and hosted by more massive systems (e.g. Giovannini & Feretti 2000; Kempner & Sarazin 2001; Cassano et al. 2007; Murgia et al. 2009). Although it is a massive cluster (M500 ~ 8.39 × 1014M), PLCK G285.0-23.7 hosts a quite small radio halo of radius ab~281\hbox{$\sqrt{a*b} \sim 281$} kpc, where a and b are the semi-major and semi-minor axis derived in Sect. 3.1. However, if we consider the correlation between the size of the radio halos and the virial radius of their hosting clusters (Fig. 7), we see that PLCK G285.0-23.7 (red point) follows the trend of other radio loud clusters. This occurs because PLCK G285.0-23.7 is at a relatively high redshift (z ~ 0.39), which means that its virial radius, and thus its radio halo, is smaller (owing to the cosmological growth and virialisation of galaxy clusters) than those of nearby clusters with the same mass.

4. Discussion and conclusions

We present high-sensitivity and wide-band (1.13.1 GHz) ATCA observations of the galaxy cluster PLCK G285.0-23.7, which reveals the presence of a diffuse radio source with an extension (major axis) of 700 kpc and a total 1.4 GHz power of 1.7 × 1024 W/Hz, as extrapolated from measurements in the 500 MHz-wide sub-band centred at 1.867 GHz.

Both the diffuse nature of the source, which is not related to blending of compact radio objects, and the comparison of our high-sensitivity tapered radio map to the X-ray image of the cluster allow us to classify the detected emission as a classical RH. First, the morphology and size of the diffuse radio source indicate a very similar distribution of the non-thermal ICM and the thermal component, traced by its X-ray bremsstrahlung radiation. Similarly to most known RHs, our newly detected halo is thus hosted by a massive cluster (8.4 × 1014M), and there are strong indications that it is a merging system.

Our wide-band full-resolution radio map of PLCK G285.0-23.7 is one of the deepest ATCA images yet published, reaching an rms level of ~11 μJy/beam with a synthesised beam size of few arcsec, i.e. the kind of sensitivity expected for the upcoming continuum ASKAP-EMU survey (Norris et al. 2011). Despite the depth and quality of our radio observations, the detection of the diffuse source required a tapered weighting (see Table 2), which has allowed the detection of one of the lowest luminosity RH in the sub-sample of known radio-loud clusters at z > 0.35 (Fig. 6). This analysis is thus a very good example of the impact that the 30 000 sq. deg (Dec < 30°) EMU survey will have on statistical studies of the non-thermal component of the ICM, if an appropriate data reduction strategy, which includes data tapering and compact source subtraction, is implemented.

thumbnail Fig. 7

Correlation between radio halo sizes and the virial radius of their host clusters (see Cassano et al. 2007). The newly detected radio halo in PLCK G285.0-23.7 is indicated in red.

The measurements of the halo radio power can be significantly biased by several factors. First of all, the quality of our map centred at ~1.4 GHz does not allow a proper detection and flux density measurement of the RH, which can instead be reliably measured in the sub-band image centred at ~1.8 GHz. Since an in-band spectral index cannot be derived from our data, a value α = 1.3 is assumed here to get the k-corrected 1.4 GHz luminosity plotted in Fig. 6. The other significant uncertainty in comparing our measurement to literature values is intrinsically related to the relative depth of the different published maps and to the method adopted for flux density measurements. Here we measure the flux density on images with point source subtraction both on the uv-data and on the image plane, showing that our different measurements are consistent within the error bars (see Table 3). Is it worth mentioning that the telescope only measures power on the spatial scales of the sampled uv-coverage, which could limit the total power measured from the radio halo. This would be especially the case if the uv-coverage happened to lack at the physical scales of the radio halo.

With these elements in mind, the newly detected RH seems to be slightly under-luminous compared to objects hosted by clusters in a similar mass range, as shown in the P1.4 GHz vs. M500 plot of Fig. 6; however, it is not a clear outlier like the RH recently discovered by Bonafede et al. (2015).

It is expected that radio halos and their hosting clusters develop complex patterns with time in the radio-mass (or X-ray luminosity) diagram as a result of the evolution of cluster dynamics and particle acceleration and spectra (e.g. Donnert et al. 2013). For most of their lifetimes these systems are expected to be “off-state” (i.e. under-luminous) or in the region spanned by the correlation leading to an apparent bimodality (e.g. Brunetti et al. 2009; Cassano et al. 2013). In this respect under-luminous systems are expected to be generally associated with both young or old mergers (Donnert et al. 2013); naively PLCK G285.0-23.7 could be in one of these stages. On the other hand the spectrum of radio halos also plays a role, with ultra-steep-spectrum radio halos, which are expected in less powerful mergers and at higher redshift (as a result of stronger IC losses), being statistically under-luminous (Cassano et al. 2006, 2010; Brunetti et al. 2008, see Fig. 6).

The new generation of more sensitive radio telescopes will allow us to discriminate between these two possible scenarios in the case of PLCK G285.0-23.7 and similarly weak radio loud clusters through deeper multi-frequency radio observations.


1

In this paper we use the convention S(ν) να, with S(ν) being the radio flux density.

Acknowledgments

We thank the referee for providing constructive comments and help in improving the contents of this paper. We warmly thank Gianluca Castignani for very useful discussions and Gabriel Pratt, Monique Arnaud, and Iacopo Bartalucci for supplying the X-ray images. G.M.A. is supported by the Erasmus Mundus Joint Doctorate Program by Grants Number 2013-1471 from the agency EACEA of the European Commission. G.M.A., C.F., and M.J.-H. acknowledge financial support from Programme National Cosmologie et Galaxies (PNCG) and Université de Nice-Sophia Antipolis – Programme Professeurs Invités 2015. M.J.-H. acknowledges the Marsden Fund administered by the Royal Society of New Zealand on behalf of the Ministry of Business, Innovation and Employment. L.P. is funded via Marsden Funding to M.J.-H. The Australia Telescope Compact Array is part of the Australia Telescope National Facility, which is funded by the Australian Government for operation as a National Facility managed by CSIRO. The results of this paper are partially based on data retrieved from the ESA Planck Legacy Archive, NASA SkyView, and the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA), and the Canadian Astronomy Data Centre (CADC/NRC/CSA).

References

  1. Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bacchi, M., Feretti, L., Giovannini, G., & Govoni, F. 2003, A&A, 400, 465 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Basu, K. 2012, MNRAS, 421, L112 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bonafede, A., Brüggen, M., van Weeren, R., et al. 2012, MNRAS, 426, 40 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bonafede, A., Intema, H. T., Brüggen, M., et al. 2014, MNRAS, 444, L44 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bonafede, A., Intema, H., Brüggen, M., et al. 2015, MNRAS, 454, 3391 [NASA ADS] [CrossRef] [Google Scholar]
  7. Brunetti, G., & Jones, T. W. 2014, Int. J. Mod. Phys. D, 23, 30007 [Google Scholar]
  8. Brunetti, G., & Lazarian, A. 2016, MNRAS, 458, 2584 [NASA ADS] [CrossRef] [Google Scholar]
  9. Brunetti, G., Setti, G., Feretti, L., & Giovannini, G. 2001, MNRAS, 320, 365 [NASA ADS] [CrossRef] [Google Scholar]
  10. Brunetti, G., Giacintucci, S., Cassano, R., et al. 2008, Nature, 455, 944 [NASA ADS] [CrossRef] [Google Scholar]
  11. Brunetti, G., Cassano, R., Dolag, K., & Setti, G. 2009, A&A, 507, 661 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Cassano, R., Brunetti, G., & Setti, G. 2006, MNRAS, 369, 1577 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cassano, R., Brunetti, G., Setti, G., Govoni, F., & Dolag, K. 2007, MNRAS, 378, 1565 [NASA ADS] [CrossRef] [Google Scholar]
  14. Cassano, R., Ettori, S., Giacintucci, S., et al. 2010, ApJ, 721, L82 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cassano, R., Ettori, S., Brunetti, G., et al. 2013, ApJ, 777, 141 [NASA ADS] [CrossRef] [Google Scholar]
  16. Clarke, T. E., & Ensslin, T. 2006, Astron. Nachr., 327, 553 [NASA ADS] [CrossRef] [Google Scholar]
  17. Cuciti, V., Cassano, R., Brunetti, G., et al. 2015, A&A, 580, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Dallacasa, D., Brunetti, G., Giacintucci, S., et al. 2009, ApJ, 699, 1288 [NASA ADS] [CrossRef] [Google Scholar]
  19. Donnert, J., Dolag, K., Brunetti, G., & Cassano, R. 2013, MNRAS, 429, 3564 [NASA ADS] [CrossRef] [Google Scholar]
  20. Farnsworth, D., Rudnick, L., Brown, S., & Brunetti, G. 2013, ApJ, 779, 189 [NASA ADS] [CrossRef] [Google Scholar]
  21. Feretti, L., Fusco-Femiano, R., Giovannini, G., & Govoni, F. 2001, A&A, 373, 106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Feretti, L., Giovannini, G., Govoni, F., & Murgia, M. 2012, A&ARv, 20, 54 [Google Scholar]
  23. Giacintucci, S., Venturi, T., Cassano, R., Dallacasa, D., & Brunetti, G. 2009, ApJ, 704, L54 [NASA ADS] [CrossRef] [Google Scholar]
  24. Giacintucci, S., Dallacasa, D., Venturi, T., et al. 2011, A&A, 534, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Giacintucci, S., Kale, R., Wik, D. R., Venturi, T., & Markevitch, M. 2013, ApJ, 766, 18 [NASA ADS] [CrossRef] [Google Scholar]
  26. Giovannini, G., & Feretti, L. 2000, New Astron., 5, 335 [NASA ADS] [CrossRef] [Google Scholar]
  27. Giovannini, G., Bonafede, A., Feretti, L., et al. 2009, A&A, 507, 1257 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Govoni, F., Feretti, L., Giovannini, G., et al. 2001, A&A, 376, 803 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Govoni, F., Murgia, M., Feretti, L., et al. 2005, A&A, 430, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Kempner, J. C., & Sarazin, C. L. 2001, ApJ, 548, 639 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kim, K.-T., Kronberg, P. P., Dewdney, P. E., & Landecker, T. L. 1990, ApJ, 355, 29 [NASA ADS] [CrossRef] [Google Scholar]
  32. Liang, H., Hunstead, R. W., Birkinshaw, M., & Andreani, P. 2000, ApJ, 544, 686 [NASA ADS] [CrossRef] [Google Scholar]
  33. Lindner, R. R., Baker, A. J., Hughes, J. P., et al. 2014, ApJ, 786, 49 [NASA ADS] [CrossRef] [Google Scholar]
  34. Macario, G., Markevitch, M., Giacintucci, S., et al. 2011, ApJ, 728, 82 [NASA ADS] [CrossRef] [Google Scholar]
  35. Macario, G., Intema, H. T., Ferrari, C., et al. 2014, A&A, 565, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Miniati, F. 2015, ApJ, 800, 60 [NASA ADS] [CrossRef] [Google Scholar]
  37. Mohan, N., & Rafferty, D. 2015, PyBDSM: Python Blob Detection and Source Measurement, Astrophysics Source Code Library [Google Scholar]
  38. Murgia, M., Govoni, F., Markevitch, M., et al. 2009, A&A, 499, 679 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Norris, R. P., Hopkins, A. M., Afonso, J., et al. 2011, PASA, 28, 215 [NASA ADS] [CrossRef] [Google Scholar]
  40. Norris, R. P., Afonso, J., Bacon, D., et al. 2013, PASA, 30, e020 [NASA ADS] [CrossRef] [Google Scholar]
  41. Petrosian, V. 2001, ApJ, 557, 560 [NASA ADS] [CrossRef] [Google Scholar]
  42. Planck Collaboration IX. 2011, A&A, 536, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Planck Collaboration XXVI. 2016, A&A, 594, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Data Analysis Software and Systems IV, eds. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, ASP Conf. Ser., 77, 433 [Google Scholar]
  45. Shakouri, S., Johnston-Hollitt, M., & Pratt, G. W. 2016, MNRAS, 459, 2525 [NASA ADS] [CrossRef] [Google Scholar]
  46. Shimwell, T. W., Brown, S., Feain, I. J., et al. 2014, MNRAS, 440, 2901 [NASA ADS] [CrossRef] [Google Scholar]
  47. Shimwell, T. W., Markevitch, M., Brown, S., et al. 2015, MNRAS, 449, 1486 [NASA ADS] [CrossRef] [Google Scholar]
  48. Steer, D. G., Dewdney, P. E., & Ito, M. R. 1984, A&A, 137, 159 [NASA ADS] [Google Scholar]
  49. Vacca, V., Govoni, F., Murgia, M., et al. 2011, A&A, 535, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. van Weeren, R. J., Röttgering, H. J. A., Brüggen, M., & Cohen, A. 2009, A&A, 505, 991 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. van Weeren, R. J., Brüggen, M., Röttgering, H. J. A., et al. 2011, A&A, 533, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. van Weeren, R. J., Fogarty, K., Jones, C., et al. 2013, ApJ, 769, 101 [NASA ADS] [CrossRef] [Google Scholar]
  53. Venturi, T., Giacintucci, S., Dallacasa, D., et al. 2013, A&A, 551, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Detection of the source using PyBDSM

The use of packages for the automatic detection of radio sources is becoming more and more necessary and common, especially in the era of the new generation of deep and wide radio surveys that will ultimately bring to the SKA a revolutionary view of the long-wavelength sky (see e.g. Norris et al. 2013). In this framework, and in order to get a complementary measure of the total flux density of the diffuse radio source, we ran the automatic source finder PyBDSM (Mohan & Rafferty 2015) on the tapered images with and without point source subtraction from the uv-data. It is worth noting that PyBDSM finds significant emission in both maps, corresponding to the diffuse radio source that we found by visual inspection and analysis on the maps (see e.g. Fig. A.1).

We used the threshold technique of the PyBDSM process_image task2, which locates islands of emission above some multiple of the noise in the image (thresh_isl parameter, set to 3σ here). This determines the region where fitting is done. In addition, we set to 5σ the source detection threshold in number of sigma above the mean (thresh_pix parameter). Finally, we activated the wavelet module of process_image (i.e. atrous_do=True), which improves the detection of diffuse sources by doing wavelet transforms at increasing scales of the residual image after subtraction of the initial fitted Gaussians. The modelled Gaussians are shown in the bottom right panel of Fig. A.1, while the top right panel shows the islands of significant emission (cyan) and the position of the fitted Gaussians. Violet ellipses indicate the sources identified on the input radio map (top left), while the red empty ellipse shows the extended source recovered through the wavelet analysis. As shown in the residual maps (bottom left), we nicely managed to fit most of the source components, in particular the central diffuse source.

The final output catalogs of PyBDSM give a list of all the Gaussian functions fitted to model the significant emission within the input radio maps and give a source list where different Gaussians are grouped together if they satisfy objective criteria to be considered as a single source. For each detected source, PyBDSM provides the values of the FWHM of the major and minor axisand the total integrated Stokes I flux density. To subtract the contribution of point sources in the ID2ss case reported in Table 3, we ran PyBDSM on the full-resolution Block 3 image of the cluster. The flux densities of compact sources obtained in this way and contained within the cyan region of Fig. A.1 were then subtracted from the total flux density of the diffuse source obtained in the ID 2 case.

The size and flux density of the newly detected radio halo are given as an output of PyBDSM, being 5.76 arcmin × 4.34 arcmin (914 kpc × 689 kpc) and 4.91 ± 0.04 mJy for the map with ID1. In the case of ID2 and ID2ss, we obtained a halo size of 4.48 arcmin × 2.70 arcmin, and flux densities of 4.08 ± 0.07 mJy and 3.72 ± 0.16 mJy, respectively.

We note here that, compared to classical measurements performed “by hand” (Sect. 3.1 and Table 3), PyBDSM gives systematically higher values of both source sizes and total flux densities. This is because with classical methods we integrate the surface brightness of the diffuse source within a region delimited by the 3σ contours, while this same region is used by PyBDSM as a support to fit one (or multiple) Gaussian function(s) giving the total flux density of the source(s) whose size is reconstructed based on moment analysis. This last method tends to include flux density coming from regions outside the original islands of significant emission, in particular when the wavelet module, which decomposes the residual image that results from the normal fitting of Gaussians into wavelet images of various scales, is activated3.

Both methods are based on different assumptions. In the case of classical “by-hand” measurements, for instance, it is assumed that the mean value of the surface brightness of all pixels within 3σ contours multiplied by the number of synthesised beams within the considered region is a proper measurements of the integrated flux density for the whole RH. Instead, in the case of PyBDSM, it is assumed that sources are correctly modelled by one or more Gaussian functions of increasing size. We do not conclude that one method is better than the other since they are based on different approaches, but we definitely recommend not directly comparing results obtained through a mix of the two different methods when producing plots such as our Fig. 6.

thumbnail Fig. A.1

Output of the PyBDSM source finder tool. The axis units are in pixels. From left to right and from top to bottom: original image on which PyBDSM is run (in this example, the ID2 map of Table 3); islands of significant emission (cyan contours) and fitted Gaussians that model identified sources overlaid on the ID2 map; residual after the subtraction of the fitted sources; model of the reconstructed sky. The big red ellipse corresponds to the diffuse emission in the centre of the cluster.

Appendix B: Additional table

Table B.1

Collection of clusters known to host a giant radio halo (see Sect. 3.2).

All Tables

Table 1

Details of our ATCA observations towards PLCK G285.0-23.7 .

Table 2

Properties of the full-resolution (first five rows) and tapered (last five rows) radio maps in different frequency bands centred at ν0.

Table 3

Physical properties of the diffuse radio source at 1.867 GHz from tapered images.

Table B.1

Collection of clusters known to host a giant radio halo (see Sect. 3.2).

All Figures

thumbnail Fig. 1

uv-coverage of PLCK G285.0-23.7 observations in the 500 MHz band centred at 1.867 MHz. The wide-band of the full 2 GHz observations completes the coverage.

In the text
thumbnail Fig. 2

Final full-resolution, wide-band image of PLCK G285.0-23.7. The outer, big circle denotes the boundary of the primary beam with a radius of ~0.22 deg. The central 1 Mpc diameter region is indicated by the inner, smaller red circle. The rms noise of this image is 11.3 μJy/beam and the resolution is 5′′.

In the text
thumbnail Fig. 3

Contours of the ATCA tapered image centred at 1.867 GHz are overlaid on the ATCA wide-band tapered (left) and full-resolution maps (right). Dashed contours correspond to the islands of significant emission detected by PyBDSM on the Block 3 tapered map (cyan contours in Fig. A.1). The continuous curves trace instead “classical” 3, 3 ×2\hbox{$\times \sqrt{2}$}, 6σ contours (i.e. with the 51.1 μJy/beam rms value of the map estimated by hand on the Block 3 tapered map at ~30 arcsec resolution; see Table 2).

In the text
thumbnail Fig. 4

Radio contours of the ATCA tapered image centred at 1.867 GHz overlaid on the smoothed raw XMM image in the [0.32.0] keV energy band of PLCK G285.0-23.7 (Planck Collaboration IX 2011). Contour levels are the same as in Fig. 3.

In the text
thumbnail Fig. 5

Zoom in the HST-ACS image of PLCK G285.0-23.7 central field (available in the HST archive). For reference, the 3σ contour of the ATCA tapered image centred at 1.867 GHz is shown in black. The blue contours indicate sources of significant radio emission in the full-resolution deep ATCA image of the cluster (i.e. 5 × rms of the final wide-band ATCA map shown in Fig. 2).

In the text
thumbnail Fig. 6

Synchrotron power of RHs at 1.4 GHz (P1.4 GHz) vs. cluster mass (M500 from Planck Collaboration XXVI 2016). Symbols are colourcoded depending on the redshift of RH hosting clusters. The red line is the fit derived in the present paper, including only the RHs with radio fluxes measured at 1.4 GHz (indicated by triangles; see also Table B.1). Analogously to Cassano et al. (2013), we exclude for the fitting ultra-steep spectrum RHs and halos detected in frequencies different from 1.4 GHz, both shown in dots (see Table B.1). The red star corresponds to PLCK G285.0-23.7 detected in this work. When the spectral index is unknown (including PLCK G285.0-23.7 ) we assume the value of α = 1.3.

In the text
thumbnail Fig. 7

Correlation between radio halo sizes and the virial radius of their host clusters (see Cassano et al. 2007). The newly detected radio halo in PLCK G285.0-23.7 is indicated in red.

In the text
thumbnail Fig. A.1

Output of the PyBDSM source finder tool. The axis units are in pixels. From left to right and from top to bottom: original image on which PyBDSM is run (in this example, the ID2 map of Table 3); islands of significant emission (cyan contours) and fitted Gaussians that model identified sources overlaid on the ID2 map; residual after the subtraction of the fitted sources; model of the reconstructed sky. The big red ellipse corresponds to the diffuse emission in the centre of the cluster.

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.