Issue |
A&A
Volume 660, April 2022
|
|
---|---|---|
Article Number | A78 | |
Number of page(s) | 40 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202143020 | |
Published online | 13 April 2022 |
The Planck clusters in the LOFAR sky
I. LoTSS-DR2: New detections and sample overview⋆
1
Leiden Observatory, Leiden University, PO Box 9513 2300 RA Leiden, The Netherlands
e-mail: botteon@strw.leidenuniv.nl
2
ASTRON, the Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA Dwingeloo, The Netherlands
3
INAF – IRA, Via P. Gobetti 101, 40129 Bologna, Italy
4
Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany
5
SRON Netherlands Institute for Space Research, Niels Bohrweg 4, 2333 CA Leiden, The Netherlands
6
Dipartimento di Fisica e Astronomia, Università di Bologna, Via P. Gobetti 93/2, 40129 Bologna, Italy
7
INAF – IASF Milano, Via A. Corti 12, 20133 Milano, Italy
8
Dipartimento di Fisica, Università degli Studi di Milano, Via Celoria 16, 20133 Milano, Italy
9
Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of Tokyo, Kashiwa, Chiba 277-8583, Japan
10
Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield AL10 9AB, UK
11
INAF – Astronomical Observatory of Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy
12
Thüringer Landessternwarte, Sternwarte 5, 07778 Tautenburg, Germany
Received:
30
December
2021
Accepted:
3
February
2022
Context. Relativistic electrons and magnetic fields permeate the intra-cluster medium (ICM) and manifest themselves as diffuse sources of synchrotron emission observable at radio wavelengths, namely radio halos and radio relics. Although there is broad consensus that the formation of these sources is connected to turbulence and shocks in the ICM, the details of the required particle acceleration, the strength and morphology of the magnetic field in the cluster volume, and the influence of other sources of high-energy particles are poorly known.
Aims. Sufficiently large samples of radio halos and relics, which would allow us to examine the variation among the source population and pinpoint their commonalities and differences, are still missing. At present, due to the physical properties of the sources and the capabilities of existing facilities, large numbers of these sources are easiest to detect at low radio frequencies, where they shine brightly.
Methods. We examined the low-frequency radio emission from all 309 clusters in the second catalog of Planck Sunyaev Zel’dovich detected sources that lie within the 5634 deg2 covered by the Second Data Release of the LOFAR Two-meter Sky Survey (LoTSS-DR2). We produced LOFAR images at different resolutions, with and without discrete sources subtracted, and created overlays with optical and X-ray images before classifying the diffuse sources in the ICM, guided by a decision tree.
Results. Overall, we found 83 clusters that host a radio halo and 26 that host one or more radio relics (including candidates). About half of them are new discoveries. The detection rate of clusters that host a radio halo and one or more relics in our sample is 30 ± 11% and 10 ± 6%, respectively. Extrapolating these numbers, we anticipate that once LoTSS covers the entire northern sky it will provide the detection of 251 ± 92 clusters that host a halo and 83 ± 50 clusters that host at least one relic from Planck clusters alone. All images and results produced in this work are publicly available via the project website.
Key words: galaxies: clusters: general / galaxies: clusters: intracluster medium / catalogs / radiation mechanisms: non-thermal / radiation mechanisms: thermal
Data are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/660/A78
© ESO 2022
1. Introduction
Radio emission associated with galaxy clusters and their member galaxies is mainly related either to radio galaxies that are powered by a central active galactic nucleus (AGN) or to nonthermal components residing in the intra-cluster medium (ICM). While AGN are often bright sources that contribute the majority of the radio flux from a cluster at gigahertz frequencies, the diffuse sources generated by relativistic electrons (Lorentz factors of γL > 1000) that propagate in the ICM magnetic field (∼μG level) have remained somewhat elusive, despite many extensive searches, due to their lower surface brightness and rapidly declining flux density with increasing frequency (e.g., Feretti et al. 2012; van Weeren et al. 2019, for reviews). The observed levels of diffuse radio emission from clusters suggest that a few percent of the energy of a cluster merger is dissipated by shocks and turbulence in the ICM and transferred to nonthermal components (the largest amount goes into ICM heating; see Markevitch & Vikhlinin 2007). However, the details of the particle acceleration and magnetic field amplification mechanisms on cluster scales are still poorly understood (e.g., Brunetti & Jones 2014, or a review). Thus, by studying the emission associated with the ICM we can probe the fundamental physics of particle acceleration in highly rarefied plasmas that are beyond the reach of those that can be studied in laboratories, and more generally we can provide insights into large-scale structure formation and evolution.
Diffuse cluster sources are typically classified as radio halos, mini-halos, relics, and revived fossil plasma sources (or phoenixes) according to their location in the cluster, morphology, size, and radio spectral properties. Observations with many facilities, such as the Very Large Array (VLA; Thompson et al. 1980), the Westerbork Synthesis Radio Telescope (WSRT; Hogbom & Brouw 1974), and the Giant Metrewave Radio Telescope (GMRT; Swarup et al. 1991), have played a crucial role in the discovery of new cluster radio sources and in constraining their main properties (e.g., Giovannini et al. 1999, 2006; Giovannini & Feretti 2000; Kempner & Sarazin 2001; Venturi et al. 2007, 2008; Rudnick & Lemmerman 2009; van Weeren et al. 2009, 2011). These instruments, in combination with X-ray observations, have provided conclusive evidence that diffuse (up to megaparsec-scale) radio sources in the ICM are connected to the dynamical motions of the ICM. Proton-proton collisions in the ICM represent an alternative process for producing (secondary) electrons in clusters (e.g., Dennison 1980; Blasi & Colafrancesco 1999); however, their contribution is likely not dominant enough to explain extended emission on megaparsec scales (e.g., Jeltema & Profumo 2011; Zandanel & Ando 2014; Brunetti et al. 2017; Adam et al. 2021). The search for correlations between diffuse radio sources and host cluster properties (e.g., Liang et al. 2000; Cassano et al. 2007, 2008; Brunetti et al. 2009; de Gasperin et al. 2014; Yuan et al. 2015) as well as their connection with the cluster dynamical state (e.g., Buote 2001; Cassano et al. 2010a; Wen & Han 2013; Cuciti et al. 2015; Giacintucci et al. 2017) is fundamental to unveiling the origin of these objects. However, until recently, many of these studies were hampered by the sensitivity of the observations, which has limited the number of detections of diffuse emission to about a hundred and the statistical analysis to very massive systems. This has challenged the overall interpretation of the population of these sources through theoretical models (e.g., Cassano & Brunetti 2005; Cassano et al. 2006; Nuza et al. 2012, 2017; Brüggen & Vazza 2020).
Thanks to the increased sensitivity to these diffuse sources that has been made possible due to upgrades to facilities such as the Karl G. Jansky Very Large Array (Perley et al. 2011) and the upgraded Giant Metrewave Radio Telescope (uGMRT; Gupta et al. 2017), as well as the advent of new-generation interferometers, such as the LOw Frequency ARray (LOFAR; van Haarlem et al. 2013), the Murchison Widefield Array (MWA; Tingay et al. 2013), the Australian Square Kilometre Array Pathfinder (ASKAP; Hotan et al. 2021), and MeerKAT (Jonas 2009), it is now possible to search for diffuse radio sources in clusters with a number of complementary and sensitive instruments. In particular, the leap forward in the capabilities of low-frequency interferometers allows us to study diffuse cluster sources in a regime where they are brighter due to their steep synchrotron spectra (α > 1, with Sν ∝ ν−α, where Sν is the flux density at frequency ν and α is the spectral index). In this respect, LOFAR has recently enabled the first detailed observations of galaxy clusters at frequencies of < 200 MHz thanks to the unprecedented high sensitivity and high resolution in its operational frequency range. This potential has already been demonstrated as LOFAR has proved to be very fruitful in investigating different aspects of nonthermal phenomena in the ICM, allowing us: to discover new instances of diffuse sources in clusters (e.g., Shimwell et al. 2016; Savini et al. 2018a, 2019; Wilber et al. 2019), including ultra-steep spectrum emission (e.g., Brüggen et al. 2018; Wilber et al. 2018; Mandal et al. 2020; Biava et al. 2021a) and very large-scale emission outside the central cluster region (e.g., Govoni et al. 2019; Botteon et al. 2019a, 2020a; Bonafede et al. 2021; Hoeft et al. 2021; Hoang et al. 2021a), as well as new faint halos and relics (e.g., Botteon et al. 2019a, 2021a; Locatelli et al. 2020; Hoang et al. 2021b) and high-z systems (e.g., Cassano et al. 2019; Di Gennaro et al. 2021a); to pinpoint the complex interplay between tailed cluster AGN and ICM motions (e.g., de Gasperin et al. 2017; Clarke et al. 2019; Hardcastle et al. 2019; Botteon et al. 2020b, 2021b; Ignesti et al. 2022a); to study the central cluster AGN duty cycle, structure, and interaction with the hot ICM (e.g., Brienza et al. 2020; Bîrzan et al. 2020; Biava et al. 2021b; Timmerman et al. 2022); and to detect extended, extraplanar emission from star-forming galaxies infalling into clusters (e.g., Ignesti et al. 2020, 2022b; Roberts et al. 2021a,b, 2022).
Sensitive wide-area searches for diffuse cluster sources require significant observational time and are arguably most efficient at low frequencies due to the higher survey speed resulting from the larger field-of-view (FoV) of the instruments. In addition, much of the undiscovered population is thought to have very steep spectra (α > 1.5; e.g., Cassano & Brunetti 2005; Cassano et al. 2006). This implies that, of the planned wide-area surveys, those at low frequencies are anticipated to make the largest number of discoveries (e.g., Cassano et al. 2010b, 2012; Nuza et al. 2012). In this respect, LOFAR is carrying out wide and deep surveys, and it is currently observing the entire northern sky at 120−168 MHz and 42−66 MHz in the context of the LOFAR Two-meter Sky Survey (LoTSS; Shimwell et al. 2017) and the LOFAR LBA Sky Survey (LoLSS; de Gasperin et al. 2021), respectively. These surveys have enormous discovery potential in many fields of astrophysics, including galaxy cluster science, offering the opportunity to study large samples of objects in synergy with surveys performed at other wavelengths. As it is currently believed that the cluster mass is a key parameter for the formation of the most extended radio sources in the ICM (namely, halos and relics), catalogs of clusters detected via the Sunyaev-Zel’dovich (SZ) effect (Sunyaev & Zel’dovich 1972) are particularly interesting as they provide unbiased samples that are almost mass-selected and that are ideal to be cross-matched with the LOFAR surveys.
This is Paper I of a series dedicated to the study of diffuse radio emission in the ICM of galaxy clusters selected from the second Planck catalog of SZ sources (PSZ2; Planck Collaboration XXVII 2016) that have also been covered by the Second LoTSS Data Release (LoTSS-DR2; Shimwell et al. 2022). It represents an extension of our previous work (van Weeren et al. 2021), which was based on the galaxy clusters covered by the First LoTSS Data Release (LoTSS-DR1; Shimwell et al. 2019). Here, we present the new sample (Sect. 2), describe the methods and data used (Sect. 3), classify the cluster radio sources (Sect. 4), provide the quantities used for the analysis that will be performed in subsequent papers, and present the new detections and the results of our study (Sects. 5 and 6). In Cassano et al. (in prep.) and Cuciti et al. (in prep.), we discuss the occurrence and the scaling relations of radio halos in the sample, while in Jones et al. (in prep.) we focus on radio relics. Other papers dedicated to the study of the X-ray properties of the sample (Zhang et al., in prep.) and to the methods developed to derive upper limits to the diffuse cluster radio emission (Bruno et al., in prep.) are also forthcoming.
Hereafter, we adopt a Λ cold dark matter cosmology, with ΩΛ = 0.7, Ωm = 0.3, and H0 = 70 km s−1 Mpc−1.
2. Cluster sample
The PSZ2 catalog (Planck Collaboration XXVII 2016) contains 1653 SZ sources detected over the entire sky. Here, we focus on the 309 entries listed in Table A.1 that lie in the LoTSS-DR2 footprint. This comprises two regions covering 5634 deg2 that are centered at approximately 12h45m00s +44°30′00″ and 01h00m00s +28°00′00″. In the PSZ2 catalog, 63 entries out of 309 are without redshift and mass estimates, meaning that they were not confirmed detections at the time of their publication. However, in follow-up optical studies by Buddendiek et al. (2015), Burenin (2017), Burenin et al. (2018), Barrena et al. (2018), Streblyanska et al. (2018, 2019), Aguado-Barahona et al. (2019), Boada et al. (2019), and Zohren et al. (2019), redshifts have been obtained for 35 of these 63 Planck detections. From these redshifts we computed M500 by interpolating the M500 versus z curves provided in the PSZ2 individual algorithm catalogs for each detection (see Appendix D in Planck Collaboration XXVII 2016). There are 28 remaining PSZ2 detections without redshift confirmation in the LoTSS-DR2 area. For simplicity, in the paper we refer to all 309 entries in Table A.1 as “galaxy clusters”, even if 28 of them should formally be referred to as “SZ detections”. In the end, our sample consists of galaxy clusters that are known to span at least the redshift and mass ranges of 0.016 < z < 0.9 (median of 0.280) and 1.1 × 1014 M⊙ < M500 < 11.7 × 1014 M⊙ (median of 4.9 × 1014 M⊙). As shown in Fig. 1, the distribution of redshift and mass in the sample of clusters included in our study provides a qualitatively good representation of the full PSZ2 population. Nonetheless, we note that our sample was selected only based on right ascension and declination cuts. For a better assessment of the similarity between the two samples, we performed a two-sample Kolmogorov-Smirnov test on M500 and z, and found p-values of the null hypothesis (that the two samples are drawn from the same distribution) of 0.25974 and 0.00148, respectively. The absolute differences between the median values of the two samples are 0.17 × 1014 M⊙ (for the mass) and 0.056 (for the redshift). These numbers indicate that the mass distributions of the two samples are in agreement. Concerning the redshift distributions, the low p-value and the slightly higher median value of the LoTSS-DR2 sample are related to the fact that our sample includes clusters that were confirmed with optical follow-ups after the publication of the Planck catalog. These clusters did not have a redshift in the original PSZ2 catalog and they are mostly clusters at high z (30 out of the 35 clusters confirmed by optical follow-ups have redshift higher than the median value of the full PSZ2 sample). The distribution of the clusters within the LoTSS-DR2 area is shown in Fig. 2.
![]() |
Fig. 1. Redshift-mass distribution of PSZ2 sources. Clusters that are located in the LoTSS-DR2 area are indicated in blue. The histograms show the number of clusters at various redshifts and masses; the dashed and dotted lines mark the median values of the LoTSS-DR2 sample and the full PSZ2 sample, respectively. Similarly to the full PSZ2 sample, our sample spans a wide range of redshifts and masses. |
![]() |
Fig. 2. Position of the PSZ2 clusters in the RA-13 (top) and RA-1 (bottom) regions covered by LoTSS-DR2. The color code indicates the redshift of the cluster. The radius of the circle is proportional to M500. Clusters without redshift and mass are reported as black crosses. The background image represents the noise variations in LoTSS-DR2 (darker colors denote higher noise values) and is reproduced from Shimwell et al. (2022). |
3. Methods and data analysis
3.1. Data reduction
LoTSS is an ongoing radio survey that employs LOFAR High Band Antennas (HBA) to observe the entire northern sky in the frequency range 120−168 MHz. LoTSS observations are generally 8 h long, the nominal central frequency of the survey is 144 MHz, and the typical root-mean-square (rms) noise σ is ∼0.1 mJy beam−1. More details on LoTSS, such as its design and scientific goals, can be found in Shimwell et al. (2017, 2019). Here, we use the data from the LoTSS-DR2 (Shimwell et al. 2022), which covers an area that is a factor of ∼13 larger than LoTSS-DR1 and has additional improvements to image fidelity particularly for faint diffuse structures. LoTSS-DR2 pointings are processed with fully automated pipelines developed by the LOFAR Surveys Key Science Project team that aim to correct for direction-independent and direction-dependent effects that are present in the data. These pipelines are PREFACTOR1 (van Weeren et al. 2016; Williams et al. 2016; de Gasperin et al. 2019) and ddf-pipeline2 (Tasse et al. 2021). The latter employs KILLMS (Tasse 2014a,b; Smirnov & Tasse 2015) and DDFACET (Tasse et al. 2018) to perform direction-dependent self-calibration of the entire LOFAR FoV, and has been significantly improved compared to the version used to process LoTSS-DR1 (Shimwell et al. 2017). We refer the reader to Tasse et al. (2021) and Shimwell et al. (2022) for more details.
In order to further improve the image quality toward the targets in our sample while also allowing for more flexible imaging, we adopted the “extraction + recalibration” scheme described by van Weeren et al. (2021), which was also used for the analysis of the galaxy clusters in the LoTSS-DR1 region. This method consists of the subtraction of the sources outside a small square region of the sky (typically, ∼0.3−0.7 deg2) containing the target from the uv data and using the direction-dependent calibration solutions and sky model derived from ddf-pipeline. The extracted data sets are then phase-shifted to the center of the region, averaged, and corrected for the LOFAR station beam in this direction. Finally, the calibration of the data is refined by performing a series of typically 4 phase and 6 phase and amplitude calibration loops. LoTSS pointings have a full width at half maximum of 3.96° at 144 MHz and are separated by ∼2.6°, so usually a specific target is covered by multiple pointings, which are combined and analyzed together. We typically extract the visibility data from pointings that are < 2.2° from the center of the extracting region.
Among the 309 PSZ2 sources in the LoTSS-DR2 area, we were not able to apply this method to 5 targets. This included the Coma cluster (PSZ2 G057.80+88.00) whose radio emission is too large for us to approximate the ionospheric and beam errors with single solutions as is done in the extraction + recalibration scheme, requiring a special treatment (see Bonafede et al. 2021). Embedded within the Coma cluster radio halo, there are a further two clusters (PSZ2 G056.62+88.42 and PSZ2 G061.75+88.11) where we are unable to differentiate their emission from that of Coma. Finally, PSZ2 G060.10+15.59 and PSZ2 G075.08+19.83 are located in regions where the direction-dependent calibration with ddf-pipeline failed likely due to very poor ionospheric conditions. These 5 targets were excluded from the analysis. A collection of the LOFAR images of our PSZ2 sample is shown in Fig. B.1.
3.2. Alignment of the flux scale
Because of inaccuracies in the LOFAR beam model, transferring amplitude solutions from the calibrator field data to the target field data may introduce offsets in the flux density scale of the target field (e.g., Hardcastle et al. 2016). For this reason, as described by Hardcastle et al. (2021) and Shimwell et al. (2022), when constructing final LoTSS-DR2 catalogs and mosaics the images are scaled to align the flux density scale with the Roger et al. (1973) scale. This procedure involves cross-matching catalogs derived from each LoTSS-DR2 observation with the NRAO VLA Sky Survey (NVSS; Condon et al. 1998) catalog and assuming a global scaling relationship between NVSS and the 6C catalog (Hales et al. 1988, 1990), which is thought to be consistent with Roger et al. (1973) to 5%. In LoTSS-DR2 the derived scaling factors are applied to the images during the mosaicing and not to the visibilities directly (Shimwell et al. 2022). Hence, for our processing, which uses the archived uv data (Sect. 3.1), we instead adopt a procedure where we align catalogs created from the images obtained from the extracted data sets with the final LoTSS-DR2 catalog in which the scaling factors have been applied. For this we perform a simple cross-match between the two catalogs (5 arcsec) and use the criteria given in Shimwell et al. (2022) to select only compact sources and remove those that are not (nearest neighbor within 30 arcsec) and those at low signal-to-noise (less than 7). As outliers can still exist in this cross-matched catalog, we used three different fitting methods (Sen 1968; Huber 1981, and regular linear regression). All three are available within the scikit-learn package3 (Pedregosa et al. 2012), and all have different outlier rejection criteria, ranging from a more robust median calculation to a less robust simple linear regression. About 90% of the time the derived values from the different methods give results that are consistent within 10%. The remaining cases are those where outliers are more prominent and are rejected differently by the adopted fitting methods. Thus, for each fit we calculated the mean absolute error and selected the method with the lowest value, which we then used to scale our images for that particular object and align it with the LoTSS-DR2 scale.
3.3. Radio images
For each cluster, we produced images at different resolutions to search for diffuse radio emission in the ICM and perform the subsequent analysis. The imaging was done with WSCLEAN v2.8 (Offringa et al. 2014) adopting the Briggs (1995) weighting scheme with robust = −0.5, and applying Gaussian uv tapers in arcsec equivalent approximately to 25, 50, and 100 kpc at the cluster redshift. For the PSZ2 entries without redshift, the Gaussian uv taper was set to 15, 30, and 60 arcsec to span a wide range of resolutions. The 60 arcsec tapered images were produced only with discrete sources subtracted (see below). For each cluster we also produced a higher-resolution image by using robust = −1.25, which leads to a resolution typically of 5.0 arcsec × 3.5 arcsec. The multi-scale multifrequency deconvolution option (Offringa & Smirnov 2017) was enabled in WSCLEAN adopting fixed scales (-multiscale-scales 0,4,8,16,32,64) and subdividing the bandwidth into 6 channels for all imaging runs.
We used the images obtained with robust = −0.5 and no uv taper as reference images to assess the quality of the data sets, which were visually inspected and graded according to: 1 for high quality images; 2 for images that are partially affected by calibration artifacts or higher rms levels; and 3 for low quality images where the scientific analysis is not possible due to strong calibration artifacts or very high rms noise levels (i.e., > 0.3 mJy beam−1). This image quality is reported in Table A.1. Radio contours from the reference images were overlaid on optical Panoramic Survey Telescope and Rapid Response System (Pan-STARRS; Chambers et al. 2016) mosaics using the g, r, i filters to verify the presence of optical counterparts.
To better study the diffuse emission, we removed the contribution of discrete sources by imaging the data sets with a uv cut corresponding to a physical scale of 250 kpc at the cluster redshift (for PSZ2 entries without redshift we arbitrarily adopted a uv cut of 2722λ, corresponding to 2.82 arcmin or 250 kpc at z = 0.2) and subtracting their clean components from the visibility data. The new visibility data were then imaged with the same values of uv taper adopted for the images obtained before discrete source subtraction. The low-resolution radio contours with discrete sources subtracted were overlaid onto the Chandra and/or XMM-Newton X-ray images (when available) smoothed to 30 kpc by a Gaussian function (for more details on the X-ray images, see Sect. 3.4).
An example of the set of images produced for each cluster is shown in Fig. 3. All images are available for download in PNG and FITS format on the project website4.
![]() |
Fig. 3. Example of the set of images that we produced for each object in our sample. The displayed cluster is PSZ2 G149.75+34.68, and the reported panels show (from left to right and from top to bottom) the reference image (robust = −0.5), the high-resolution image (robust = −1.25), taper 25, 50, and 100 kpc images with and without discrete sources, the clean model components used for the source subtraction, contours from the reference radio image starting from 3σ overlaid on an optical image (Pan-STARRS g, r, i), and taper 50 and 100 kpc discrete-source-subtracted contours starting from 2σ overlaid on an X-ray image (Chandra or XMM-Newton). Contours are always spaced by a factor of 2. The beam is shown in the bottom-left corner, and the mass (M500, in ×1014 M⊙ units), redshift (z), and image noise (rms, in mJy beam−1 units) are reported in the top-left corner. The circle denotes r500 and is centered at the coordinates reported in the PSZ2 catalog. The images are available at full resolution as well as in FITS format on the project website (https://lofar-surveys.org/planck_dr2.html). |
3.4. X-ray images and morphological parameters
As we will describe in Sect. 4, for the purpose of classifying the detected cluster diffuse radio sources, it is crucial to compare the position of the extended radio emission with the other cluster components and especially the ICM, which can be traced by its X-ray emission and by the SZ effect. Since all clusters in our sample have been detected by Planck, two-dimensional maps of the SZ signal are available for all of them (Planck Collaboration XXII 2016) but their use is hampered by their 10 arcmin resolution, which does not allow us to spatially resolve most of the targets of the sample. We thus decided to map the ICM distribution with the X-ray images obtained by the current generation X-ray telescopes (Chandra and XMM-Newton), whose spatial resolution is higher or comparable to that of our radio images (Sect. 3.3). We searched the Chandra and XMM-Newton archives for observations of the targets in our sample and we retrieved the data for 115 and 100 clusters, respectively (72 targets have been observed both by Chandra and XMM-Newton). The procedures used to prepare the X-ray images for each instrument are described in Sects. 3.4.1 and 3.4.2.
The ICM distribution is very sensitive to the dynamical history of the clusters, and therefore quantitative measurements of the morphology of the X-ray emission of galaxy clusters have proved to be an effective way to characterize the dynamical state of large samples of galaxy clusters (e.g., Buote 2001; Santos et al. 2008; Cassano et al. 2010a; Rasia et al. 2013; Parekh et al. 2015; Rossetti et al. 2017; Lovisari et al. 2017, and references therein). The use of a set of morphological parameters of the X-ray emission of clusters is only an approximation to the daunting task of assessing the dynamical state of the cluster. However, the combination of two morphological parameters is effective to provide a relatively robust classification. In particular the optimal choice is to combine a parameter sensitive to the presence of substructure, the centroid shift, with a parameter more sensitive to the core properties, the concentration parameter (Lovisari et al. 2017). In fact this has been the usual choice made in previous studies of the classification of radio sources (for recent examples, see Cuciti et al. 2015, 2021). The physical scale over which the morphological parameters are measured is also an important factor: here, following previous studies (Cassano et al. 2010a; Cuciti et al. 2021), we analyze an aperture of Rap = 500 kpc centered on the peak of the X-ray emission.
The concentration parameter has been introduced by Santos et al. (2008) as the ratio of the flux within two circular apertures to effectively identify cool cores even at high redshift. Here we adopt the choice of apertures made by Cassano et al. (2010a)
where F(r < 100 kpc) is the flux within 100 kpc and F(r < Rap) is the flux within the aperture of 500 kpc. The error on this parameter is obtained by taking into account the Poisson noise in both the source and background images.
The centroid shift (Mohr et al. 1993; Poole et al. 2006) is defined as the variance of the separation between the X-ray peak and the centroid of the emission obtained within a number N of apertures of increasing radius out to Rap,
where Δi is the distance between the X-ray peak and the centroid of the ith aperture. It traces the variation in the position of the centroid introduced by the presence of substructures in the X-ray emission. The number N of apertures is fixed at 20 in the XMM-Newton analysis and it is given by the number of annuli of fixed 5 arcsec width within 500 kpc for the Chandra analysis. The error on this parameter is obtained by using a Monte Carlo approach: for the Chandra analysis we simulated 100 realizations of the X-ray images obtained by resampling the counts per pixel according to their Poisson error, performed the measurement on the simulated image and estimated the standard deviation of the distribution of w thus obtained; for the XMM-Newton analysis we simulated 10 000 realizations of the centroids of the 20 apertures, sampled within their statistical errors.
We measured the concentration parameter and the centroid shift for 105 PSZ2 objects with Chandra and for 98 PSZ2 objects with XMM-Newton as a result of the following four selections: (i) a low redshift cut to accommodate the aperture of 500 kpc within the FoV of each respective detector (z > 0.065 for ACIS-I, z > 0.072 for ACIS-S, and z > 0.035 for XMM-Newton); (ii) PSZ2 G165.95+41.01 does not have a measurement because of a possible incorrect redshift estimate: the X-ray and optical images suggest that this object is at a higher z compared to the value of z = 0.062 reported in the PSZ2 catalog (see the note in the catalog that discusses a superposition with a z = 0.21 object); (iii) PSZ2 G067.52+34.75 does not have a Chandra morphological measurement because the observation is performed in a sub-array mode that does not allow a 500 kpc aperture to be accommodated; and (iv) PSZ2 G126.27+51.61 does not have a measurement because the available Chandra observation is too shallow for the faint emission of this high redshift (z = 0.815) cluster.
There are 63 objects that have both Chandra and XMM-Newton measurements and the total number of PSZ2 clusters in our sample for which we have X-ray morphological parameters is 140. For objects for which different clumps of X-ray emission could be clearly distinguished we measured morphological parameters for each component, labeling them according to their position in the sky. This explains why Table A.2, where we list the c and w, has 150 entries. For the 63 objects (65 measurements including the objects with multiple clumps) with both Chandra and XMM-Newton measurements we provide combined morphological parameters according to the following equations:
where 𝒫 is either c or w and the error σ𝒫 on this combined parameter is given by the sum of the statistical σ𝒫, stat and systematic σ𝒫, sys error,
where we also take into account the covariance with the term 2σ𝒫, xmmσ𝒫, chandra. These values are used to discuss the occurrence of diffuse radio emission with the cluster dynamical state in Cassano et al. (in prep.). A thorough comparison between parameters derived from XMM-Newton and Chandra as well as the detailed analysis on the X-ray data will be presented in Zhang et al. (in prep.).
We describe the reduction and analysis steps used for the Chandra and XMM-Newton data in the following subsections.
3.4.1. Chandra data reduction and analysis
We analyzed Chandra data with the Chandra Interactive Analysis of Observations (CIAO) software v4.13 using CALDB v4.9.4 (Fruscione et al. 2006), reprocessing data from the level 1 event files and following the standard data reduction threads5. We reprocessed event files using the chandra_repro tool and soft proton flares were excluded with the deflare task with the lc_clean routine analyzing the light curves extracted from the S2 chip when in ACIS-I configuration and from the S3 chip when in ACIS-S configuration. We used the fluximage tool to produce images in the 0.5−7.0 keV bands and the appropriate exposure and point spread function maps. For the purpose of background subtraction we used the blanksky and blankskyimage tools to provide a corresponding background image to be subtracted. We detected the point sources using the wavdetect tool, and by means of dmfilth we replaced their emission with the mean count rate in a surrounding annulus. Images were smoothed to a resolution of 30 kpc at the cluster redshift to minimize the effect of having different physical sizes for the same pixel scale given the broad redshift range of our sample (see the discussion in Yuan & Han 2020). The peak of the X-ray image used as the center of the cluster has been selected as the brightest pixels in the smoothed image after masking the point sources. When multiple observations were available for the same object, we used the observation with the longest available exposure time. For all cases, our single ObsID images have at least 500 counts, which is a safe limit to have a robust measurement of the X-ray morphological parameters (e.g., Nurgaliev et al. 2013).
3.4.2. XMM-Newton data reduction and analysis
We used XMM-Newton Science Analysis System (SAS) v18.0.0 for XMM-Newton European Photon Imaging Camera (EPIC) data reduction (Gabriel et al. 2004). MOS and pn event files are obtained from the observation data files with the tasks emproc and epproc. The out-of-time (OoT) event files of pn are produced by epproc as well. We extracted count images in the 0.5−2.0 keV band. The OoT count maps are directly subtracted from each pn count image following the user guide6. The corresponding exposure maps were generated using task eexpmap with parameter withvignetting = yes. Each exposure map was then multiplied by the on-axis effective area calculated by arfgen.
We used stacked filter wheel closed (FWC) event files as non-X-ray background (NXB). For each ObsID, the FWC event files were re-projected using the task evproject to match the observation. For the two MOS detectors, the NXB count maps were scaled using the blank regions out-of-FoV. For the pn detector, because of the contamination in out-of-FoV regions, we estimated the scaling factor based on the long-term NXB variation due to solar activities, which will be detailed described in Zhang et al. (in prep.).
Point source detection and removal procedures are the same as the Chandra data analysis. For each object, we stacked the point-source-removed NXB-subtracted count images and exposure maps, respectively. The stacked net count image was then divided by the stacked exposure map to obtain the final point source free flux map for morphological analysis. The X-ray peak of each object is determined from the 30 kpc Gaussian smoothed flux map.
4. Classification of radio sources
For each object listed in Table A.1 we searched for diffuse radio sources in the ICM that are not clearly associated with any AGN by visually inspecting the set of LOFAR images at different resolutions (with and without source subtraction) together with the optical and X-ray overlay images. To make the classification of the radio emission as objective as practical at present and easily reproducible, we created a decision tree (Fig. 4) that we followed during the inspection of the images and classified each cluster. Below we define the six classes of objects that form the end points of our decision tree.
![]() |
Fig. 4. Decision tree used to classify the diffuse radio sources in the ICM. The classes of objects that form the end points of the decision tree are described in Sect. 4. |
“Radio halos” (RH) are extended sources that occupy the region where the bulk of the X-ray emission from the ICM is detected. Historically, they were divided according to their sizes: giant halos extended on cluster-scale and mini-halos covering the cluster central region. Since in this work we are dealing with a large sample of clusters with masses spanning over one order of magnitude of difference, we prefer to not separate mini-halos from giant-radio halos based on the size of the radio emission; instead, we refer generically to radio halos.
“Radio relics” (RR) are elongated sources whose position is offset from the bulk of the X-ray emission from the ICM and consistent with lying in cluster outskirts. We also checked for a sharp edge in the radio emission and a largest-linear size (LLS) ≳300 kpc.
The classification “candidate radio halo” or “candidate radio relic” (cRH or cRR) is used when Chandra or XMM-Newton X-ray observations are not available and as such the presence of a halo or relic cannot be firmly claimed because we are uncertain if the emission is in the central or outer region of the ICM. However, it is possible to make an assessment based on the position of the radio emission with respect to the apparent overdensity of galaxies in the optical image. We thus classify the emission as being a candidate radio halo or radio relic if it is consistent with the other properties of these types of sources but colocated with or offset from the overdensity of galaxies rather than the X-ray emission from the ICM.
The “uncertain” (U) classification is for objects whose emission was significantly affected by calibration and/or subtraction artifacts or which did not show a morphology, size, and/or position consistent with the categories of radio halos and relics.
“No diffuse emission” (NDE) indicates that these objects do not show the presence of diffuse emission that could be attributed to the ICM (although they may show the presence of lobes or tails from AGN in the field).
Finally, “not applicable” (N/A) is used if the image quality is 3 (see Sect. 3.3) and as such the object cannot be adequately classified because of the poor data quality.
The classification for each target in the sample performed following the decision tree of Fig. 4, where the endpoints reflect the definitions above, is reported in Table A.1. If a target showed more than one diffuse source not related with any AGN, we repeated the decision tree for each source separately. All the answers to the questions of the decision tree are available on the project website. A gallery highlighting the wide variety of different types of emission is shown in Fig. 5.
![]() |
Fig. 5. Collection of clusters that show several types of radio emission. PSZ2 G053.53+59.52 has a central radio halo and a number of sources of uncertain origin. PSZ2 G071.21+28.86 has a double radio relic. PSZ2 G088.53+41.18 is a system without diffuse cluster emission as the only emission detected is associated with an optical galaxy. PSZ2 G097.72+38.12 has a radio halo. PSZ2 G113.91–37.01 has a radio halo and two relics. PSZ2 G114.99+70.36 has emission of uncertain origin. PSZ2 G143.44+53.66 has emission of uncertain origin. PSZ2 G148.36+75.23 has a radio halo. PSZ2 G166.62+42.13 has a radio halo and multiple relics. We remark that the classification was done by inspecting all the sets of images available, while the images displayed in this gallery are only the reference ones (i.e., those with robust = −0.5 and no uv taper). For a more complete picture of the large variety of radio structures observed in our sample, we refer the reader to Fig. B.1 and to the project website. |
We note that during our classification we did not attempt to identify radio phoenixes. These complex sources trace AGN radio plasma that has an ultra-steep (α ≳ 1.8) spectrum and that is thought to have been reenergized through processes in the ICM, unrelated to the radio galaxy itself. New instances of this class of objects are rapidly emerging thanks to sensitive observations at low frequencies (e.g., de Gasperin et al. 2017; Kale et al. 2018; Mandal et al. 2019, 2020; Duchesne et al. 2020, 2021, 2022; Botteon et al. 2021b; Hodgson et al. 2021, for recent works). Our choice of not yet attempting to classify radio phoenixes is twofold. First, their definition is tightly connected to the spectral index, requiring the analysis of multifrequency observations. Second, their amorphous morphology and small sizes (at most a few hundred kiloparsecs) make a robust classification even more challenging. That said, in Table A.1 we do highlight several instances where a suspected radio phoenix was noted to form a very prominent part of the cluster emission (for example, when it is the dominant emission in the cluster volume).
5. Flux density measurements
In the following subsections, we describe how we computed the total flux density Sν of the radio halos and relics in our sample. These values were used to calculate the k-corrected radio power of the sources at 150 MHz as
where DL is the luminosity distance and . Because our integrated flux density measurements are obtained at 144 MHz, S150 is only marginally affected by the adopted (unknown) radio spectral index. In this respect, we assumed α = 1.3 for radio halos and α = 1.0 for radio relics, which are typical values found in the literature for these kinds of sources (see, e.g., Feretti et al. 2012; van Weeren et al. 2019). Instead, the rest-frame radio power is more sensitive to the spectral index, especially for sources at higher redshift, due to the k-correction term appearing in Eq. (5). The interested reader should take P150 with more caution particularly for the radio halos in the sample at high z, which could be characterized by values of α > 1.3 (Di Gennaro et al. 2021b).
The quantities derived for the radio halos and relics in our sample are reported in Tables A.3 and A.4, respectively. In the tables, uncertainties on S150 (hence, on P150) take into account the statistical (σstat), systematic (σsys), and subtraction (σsub) errors, which were summed in quadrature. The statistical error
is related to the rms of the image in the integration area (and, in the case of halos, to the fitting process; see Sect. 5.1). The systematic error,
is given by the uncertainty of the flux scale calibration, δcal, which is set to 10% for LoTSS-DR2 (see Sect. 3.2 and Shimwell et al. 2022). The subtraction error,
takes into account the presence of possible residuals in the diffuse emission region from the discrete sources that were subtracted in the uv-plane. To determine the fraction of residual contaminating emission, ξres, we visually inspected the images used for the subtraction for a subsample of clusters characterized by different values of flux density in discrete sources, and found that the following percentages,
provide a good approximation to quantify the level of contamination in our measurements.
5.1. Radio halos
We employed the Halo-Flux Density CAlculator7 (HALO-FDCA; Boxelaar et al. 2021) to measure the integrated flux density from the observed radio halos. This code fits the two-dimensional surface brightness profile with a Markov chain Monte Carlo (MCMC) method that estimates the best-fit parameters and associated uncertainties. As proposed by Murgia et al. (2009), for the fitting we assume exponential profiles in the form
where the fitted parameters are I0, which is the central brightness, and G(r), which is a function that determines the model morphology (i.e., circular, elliptical or skewed; see Boxelaar et al. 2021 for more details). As done for the LoTSS-DR1 cluster sample (van Weeren et al. 2021), we primarily used a simple circular exponential model to fit the discrete-source-subtracted images that were obtained with a Gaussian uv taper corresponding to 50 kpc at the cluster redshift. If the signal to noise in these images was low or the emission was clearly non circular, we instead used the discrete-source-subtracted images with uv taper of 100 kpc or of the elliptical exponential model. In total, the circular model has four free parameters: I0, the coordinates of the center (x0 and y0), and a single e-folding radius (r1). The elliptical model has two additional free parameters: a second e-folding radius (r2) and a rotation angle (ϕ).
Prior to fitting the radio halos we carefully examined the images for contaminating extended sources that had been poorly subtracted, such as tailed radio galaxies and radio relics. We also identified regions affected by residual calibration errors. These problematic regions were manually masked during the fitting. In order to reduce the processing time and the size of the area to manually inspect for masking, we provided as input to HALO-FDCA images with a FoV reduced to approximately 1.5r500 × 1.5r500.
A qualitative assessment of the fit quality can be done by inspecting the residual images and corner plots produced by HALO-FDCA, which we have made publicly available for each cluster on the project website. An example of these plots for the prominent radio halo that is hosted in the cluster PSZ2 G149.75+34.68 (see Fig. 3) is reported in Fig. 6. We note that HALO-FDCA returns the signal-to-noise ratio (S/N) and the of the fit, which can also be used to assess the fit quality. However, we urge caution in interpreting these values. The S/N in HALO-FDCA is defined as the integrated flux density value divided by the uncertainty introduced by the fit (which is based on the image noise and the number of data points related to the area of the halo), and hence it is not a parameter that determines the significance of the detection. The
instead is calculated from the difference between the fitted exponential model and the data, and high
values may suggest that a halo has a lot of substructure. For example, this is the case of the halo shown in Fig. 6, where
. We note a positive trend between S/N and
for the radio halos and candidate radio halos in our sample (Fig. 7), indicating that the exponential model, while being a reasonable description of the data, may be affected by the presence of asymmetries and substructures that become evident only when the data are of sufficient statistical quality.
![]() |
Fig. 6. Results obtained by fitting the radio halo shown in Fig. 3 with HALO-FDCA (Boxelaar et al. 2021). Top figures: image used for the fit with overlaid: the contours (white circles) of the best-fit circular model drawn at [1, 2, 4, 8, …] × σ (left panel), the image of the best-fit model (central panel), and the residual image of the fit, with the 2σ model contour denoted by the black circle (right panel). Contaminating sources are masked out and are highlighted by the green and red regions (left and right panels, respectively). Bottom figure: MCMC corner plot presenting the distribution of the posteriors of each fitted parameter (Foreman-Mackey et al. 2013). |
![]() |
Fig. 7.
|
To demonstrate the performance of HALO-FDCA, in Appendix C we compare the integrated flux density obtained by integrating the flux density within a circle (or ellipse, depending on the model used in HALO-FDCA) that roughly encompasses the 2σ contour of the radio halo (S2σ) with that obtained with HALO-FDCA (Sfit). We generally see good agreement between the two quantities, but, as detailed in the appendix, we found that the fitting could not be done reliably due to the low significance of the emission for 10 out of the 83 fitted radio halos and candidate halos. These sources are reported as RH*/cRH* in Table A.1 and their integrated flux densities cannot be determined accurately with current data. For all halos except PSZ2 G107.10+65.32 and PSZ2 G139.18+56.37 we use as reference the HALO-FDCA derived integrated flux density. For these two clusters there is significant substructure and we instead derived the integrated flux density by summing the pixels within the 2σ contour level.
As suggested by Murgia et al. (2009), when calculating the HALO-FDCA derived integrated flux densities we integrated the best-fit models up to a radius of three times the e-folding radius. This choice leads to a flux density that is ∼80% of the one that would be obtained by integrating the model up to infinity and is motivated by the fact that halos do not extend indefinitely.
The quantities derived for the radio halos and candidate radio halo in our sample are reported in Table A.3.
5.2. Radio relics
The integrated flux density of radio relics was computed in polygonal regions encompassing the 2σ contour of the diffuse emission. In the process, particular care was devoted to excluding possible artifacts due to calibration errors and/or other contaminating extended sources in the cluster. As done for halos, we primarily used the images obtained with a Gaussian uv taper corresponding to 50 kpc at the cluster redshift with discrete sources removed. In a handful of relics (PSZ2 G089.52+62.34 N2, PSZ2 G091.79–27.00, PSZ2 G113.91–37.01 S, PSZ2 G166.62+42.13 E, and PSZ2 G205.90+73.76 N/S), we instead used the discrete-source-subtracted images with uv taper of 100 kpc, where the integrated flux density of the relics is > 10% than that obtained in the 50 kpc images. From inspection of the images of PSZ2 G190.61+66.46, we found that the relic was partially included in the model for source subtraction and hence that a small fraction of the relic emission was removed. Since there are no compact sources in the region of the relic, only for this target we chose to use the 50 kpc image prior to the source subtraction to measure its properties.
We estimated the LLS of each relic by computing the distance between the pixels above 2σ with the largest separation in the diffuse emission. Owing to the reliability of the flux density within a beam, the error in the measured LLS corresponds to the size of the restoring beam of the image.
The downstream extension of radio relics, or relic width, often varies significantly across the extent of the relic. Additionally, the irregular shapes of some radio relics make it challenging to decide at which position to define and measure the width. We therefore decided to take a statistical approach to measuring the width of a relic. The straight line joining the LLS pixel pair typically lies approximately perpendicular to the width direction. This was verified by eye. Therefore, by measuring the largest separation between relic pixels along the line perpendicular to the LLS at each pixel along the LLS, we obtained an estimate of the width at all points along the relic. We then took the median of these width measurements as the characteristic relic width and adopt as error the standard deviation of the measurements, which reflects the spread of values obtained.
We used the coordinates of X-ray centroids as reference points for calculating the projected distance of the radio relics from the clusters. The X-ray centroid of each cluster was calculated within a region of r500 centered at the coordinates reported in the PSZ2 catalog. Only for PSZ2 G107.10+65.32 S, we used a region of 500 kpc centered at the X-ray peak due to the PSZ2 coordinates being outside this subcluster.
Due to possible projection effects, the shock front associated with the relic may not be located at the sharp surface brightness edge of the relic, but at the brightest region of the relic. We therefore chose to define the position of the relic as the center of the brightest 10% of pixels within the relic, weighted by their flux density. The distance between the relic and the centroid of the cluster X-ray emission (DRR − c) as well as the separation between double relics (DRR − RR) were computed from these coordinates. We include an additional error to account for the projection between the relic and X-ray centroid axis. Since the merger axis of clusters hosting double radio relics is approximately on the plane of the sky, the additional error corresponds to a projection of 10°. For all other relics we take an offset of 30°.
The quantities derived for the radio relics and candidate radio relics in our sample are reported in Table A.4.
6. Results and discussion
6.1. Number and distribution of sources
With 309 objects, the LoTSS-DR2/PSZ2 sample represents the largest statistical sample of galaxy clusters observed with highly sensitive low-frequency observations that has been ever used to search for and study diffuse synchrotron emission in the ICM. The number of clusters divided per category as defined in Sect. 4 is summarized in Fig. 8. We found 73 clusters hosting a radio halo (i.e., 53 RH and 20 cRH), 26 clusters hosting one or more relics (i.e., 20 RR and 6 cRR), and 47 clusters with diffuse emission of uncertain origin. Additional 10 clusters are found to host a radio halo from visual inspection but since the surface brightness profile fittings were of poor quality they were reported with an asterisk in Fig. 8 and Table A.1 (these are the 6 RH* and 4 cRH* discussed in Appendix C). No diffuse radio emission from the ICM was found for 140 targets while for 36 objects it was not possible to investigate the presence of diffuse emission either due to the impossibility of applying the extraction + recalibration method described in Sect. 3.1 (5 out of 36) or due to the bad image quality (31 out of 36). Notes on individual clusters (including references on previous studies employing radio observations) are reported in Appendix D.
![]() |
Fig. 8. Summary of the number of clusters classified in our sample divided per category (NDE = no diffuse emission; U = uncertain; RH = radio halo; RR = radio relic; c = candidate; * = flux density measurement not reliable; N/A = not applicable). We note that a cluster can be classified under multiple categories, such as both RH and RR. The percentage on the y axis is with respect to the total number of PSZ2 in LoTSS-DR2 (309 targets). |
The number of clusters with newly discovered radio halos and relics in our sample is 50 (i.e., 50% of the total number of halo and relic detections) and these clusters are highlighted in Table A.1. These new discoveries consist of 35 clusters hosting a radio halo (i.e., 21 RH and 14 cRH, namely 48% of the number of halo detections) and 15 clusters hosting one or more radio relics (i.e., 9 RR and 6 cRR, namely 58% of the number of relic detections). In addition, all the 10 halos classified with an asterisk (i.e., 6 RH* and 4 cRH*) except PSZ2 G118.34+68.79 (see van Weeren et al. 2021) are reported in this work for the first time. In Fig. 9 we show the mass and redshift distribution for the clusters hosting halos, (one or more) relics, and uncertain sources in our sample.
![]() |
Fig. 9. Histograms of the mass (top panels) and redshift (bottom panels) distribution for the main categories of objects classified in this work. |
We detected radio halos over the entire M500 range of the PSZ2 sample we studied, which spans nearly one order of magnitude. The bulk of detections occur in the mass range 5 × 1014 M⊙ ≲ M500 ≲ 7 × 1014 M⊙. In Cassano et al. (in prep.) we discuss the occurrence of halos with the cluster mass, properly taking into account the selection effects and completeness of the PSZ2 catalog. Of particular interest is the large number of detections of halos in clusters with a mass M500 < 5 × 1014 M⊙ (i.e., 17 RH, 4 cRH, 3 cRH*, and 1 cRH*). This mass regime is poorly studied and it has been disclosed only very recently thanks to sensitive observations with new generation instruments (e.g., Hlavacek-Larrondo et al. 2018; Hoang et al. 2019; van Weeren et al. 2021; Botteon et al. 2019a, 2021a,b; Duchesne et al. 2021, for recent works). We note that the candidate radio halo in PSZ2 G192.77+33.14 [M500 = (1.66 ± 0.20)×1014 M⊙], if confirmed, would be the least massive system presently known to host a radio halo. Radio halos in our sample cover also a wide redshift range, from z = 0.05 (PSZ2 G192.77+33.14) up to z = 0.888 (PSZ2 G160.83+81.66; see Di Gennaro et al. 2021a,b), with the bulk located at 0.25 ≲ z ≲ 0.4. The discovery of bright radio halos in systems (partially overlapping with our sample) with z > 0.6 has been recently made possible thanks to LoTSS observations (Cassano et al. 2019; Di Gennaro et al. 2021a). The evolution of the radio halo properties with redshift will be discussed in Cassano et al. (in prep.). Finally, we note that in Fig. 9 we have not reported the candidate radio halo in PSZ2 G144.23–18.19, which is our only detection in a cluster without redshift and mass, while the double radio halo cluster in PSZ2 G107.10+65.32 (see Botteon et al. 2018) was counted only once.
The radio relic sample has a broad distribution both in terms of M500 and z (despite being more limited in size with respect to that of halos). Unlike halos, we find only one radio relic in clusters at z > 0.6, namely in PSZ2 G069.39+68.05 (z = 0.762). This relic is claimed here for the first time. The least massive cluster in our sample is also that with the lowest redshift: PSZ2 G089.52+62.34 [M500 = (1.83 ± 0.19)×1014 M⊙ at z = 0.07], which was already reported in van Weeren et al. (2021). The total number of relic sources is 35, and the number of clusters hosting more than one relic in the ICM is 8. We only found 6 cases of classic, symmetric (i.e., located in diametrically opposite sides of the cluster), double radio relic systems (PSZ2 G071.21+28.86, PSZ2 G099.48+55.60, PSZ2 G113.91–37.01, PSZ2 G165.46+66.15, PSZ2 G181.06+48.47, PSZ2 G205.90+73.76). We remark that clusters hosting more than one relic are reported only once in Fig. 9. The statistical analysis of radio relics is presented in Jones et al. (in prep.).
Clusters with uncertain diffuse emission are mainly observed for z < 0.4 and M500 < 6 × 1014 M⊙. From visual inspection, we expect that a large fraction of these systems may host a radio halo (which, however, cannot be firmly claimed at the moment), while a smaller fraction may contain different kinds of sources in the ICM. Follow-up observations are required to confirm the nature of the emission observed. These clusters are not considered to derive scaling relations in Cuciti et al. (in prep.) and Jones et al. (in prep.), while for the study of the occurrence of radio halos in Cassano et al. (in prep.) we consider the two extreme cases where all or none of the uncertain sources are radio halos.
6.2. Classification
In this work and in the accompanying papers of the series we focus on the detection and characterization of radio halos and relics in the PSZ2 sample. These types of emissions are the most widely studied classes of diffuse synchrotron sources in the ICM, and are arguably the best defined and understood.
Diffuse sources in our sample were classified by visually inspecting LOFAR 120−168 MHz images at different resolutions, with and without discrete sources removed. We also made use of overlays of our radio data with optical and X-ray data when available. The most challenging aspect of the inspections was to disentangle the sources of interest (halos and relics) from other contaminating sources such as AGN and phoenixes and also from calibration artifacts. This separation relied upon the high resolution of LOFAR and availability and high quality of the auxiliary data sets. However, in some cases it still remained very challenging to conclusively classify the emission. Indeed, low-frequency observations are very sensitive to low-energy relativistic electrons in the ICM, which even if emitted from a single compact region can occupy a significant fraction of the cluster volume due to their long lifetimes. As a consequence, the filling factor of emission associated with AGN increases at low frequency, and diffuse emission of somewhat uncertain origin and low significance can be identified almost in all clusters. This highlights the role of seed electrons in the formation of cluster-scale diffuse sources, which are an important ingredient for (re)acceleration models.
Although we found the application of the decision tree (Fig. 4) useful to make the classification less biased and more reproducible, we acknowledge that the visual inspection of the 309 objects in our sample is a painstaking procedure. It is clear that the strategy used in this work needs further development to properly assess significantly larger cluster samples or to eventually conduct a blind search for cluster emission. In this respect, machine-based techniques represent an appealing solution to classify the emission in large object samples (e.g., Aniyan & Thorat 2017; Alhassan et al. 2018; Domínguez Sánchez et al. 2018; Lukic et al. 2018, 2019; Sadeghi et al. 2021; Vavilova et al. 2021). As part of this work, we have made public all our images and the detailed results of our decision-tree-based classification, which we hope can provide a good training set for algorithms that attempt either the full classification or to aid the automation at specific intersections in a decision-tree-type approach. Another approach would be establishing a citizen science project dedicated to galaxy clusters (i.e., a LOFAR Galaxy Clusters Zoo) with the goal of distributing the classification work among citizen scientists, similarly to what is currently done by the LOFAR Radio Galaxy Zoo8 project, which has the aim of joining different components of the same radio galaxies observed by LoTSS and identifying their optical counterparts.
As a final note, the unprecedented high resolution and high sensitivity at low frequencies provided by LOFAR have the potential of probing previously unseen features of halos and relics, but also to unveil new kinds of emission in the ICM. Collecting the large variety of source morphologies in clusters in an automated way (e.g., Gheller et al. 2018; Mostert et al. 2021) will then help us to study the nature of the sources and the role of the surrounding environment to shape the radio emission.
6.3. Prospects
The fraction of identifications in each of our classifications with respect to the full LoTSS-DR2 sample is shown Fig. 8. In the following, reported fractions were computed excluding the 36 targets of the sample for which we could not verify the presence of diffuse emission in the radio images (i.e., 5 not extracted and 31 N/A). The reported uncertainties are computed based on Poisson statistics.
The fractions of PSZ2 clusters in which a radio halo or one or more radio relics is claimed in this work are 19 ± 9% and 7 ± 5%, respectively, and increase to 27 ± 10% and 10 ± 6% if candidates are also included. If the 10 halos classified with an asterisk for which we could not provide a reliable flux density measurement (Appendix C) are also taken into account, the fraction of clusters with a radio halo detected rises further, to 30 ± 11%. The fraction of clusters hosting sources with uncertain origin is 17 ± 8%, while that of targets that do not show the presence of diffuse synchrotron emission in the ICM is 51 ± 14%.
In our earlier work (van Weeren et al. 2021) we presented the first statistical analysis of galaxy clusters observed with LoTSS and studied the 26 PSZ2 clusters residing in the ∼420 deg2 covered by LoTSS-DR1 (this region is contained within LoTSS-DR2). We found that 73 ± 15% of PSZ2 clusters hosted some type of diffuse ICM related radio emission. More specifically, 62 ± 15% and 27 ± 10% PSZ2 clusters hosted a radio halo or one or more relics (including candidates), respectively. When comparing the fraction of diffuse emission in LoTSS-DR1 with that of LoTSS-DR2 one should take into account three important factors. Firstly, the LoTSS-DR1 area is on average more sensitive than LoTSS-DR2, as it covers a sky area above declination +45°, while LoTSS-DR2 goes down to declination +16°, where the sensitivity of LOFAR is reduced (Shimwell et al. 2022). This observational limitation makes the detection of faint ICM emission more challenging and likely biases low the fraction of extended sources in LoTSS-DR2 compared to LoTSS-DR1. Secondly, the area covered by LoTSS-DR1 comprises a relatively small region of the sky, about 2% of the northern sky, and may therefore be affected by cosmic variance. LoTSS-DR2 instead spans about 27% of the northern sky, enabling the possibility of performing reliable statistical studies. Thirdly, 28 PSZ2 entries in the LoTSS-DR2 sample (i.e., 9% of the total) lack z and M500, and are not confirmed galaxy clusters at the moment, while all the PSZ2 detections in LoTSS-DR1 are confirmed galaxy clusters. If the sources without z and M500 in LoTSS-DR2 were Planck false detections, they should be rejected from the sample, resulting in an increase in the fraction of clusters hosting diffuse radio emission in the ICM (as all these 28 entries except PSZ2 G144.23-18.19 have been classified as NDE). We checked the quality flag Q_NEURAL reported in the PSZ2 catalog for these 28 entries and found that 15 have Q_NEURAL < 0.4, which the threshold below which Aghanim et al. (2015) classify low reliability detections. This suggests that about half of the clusters for which we do not have z and M500 are likely spurious Planck detections.
The area covered by LoTSS-DR2 is ideal to study the extragalactic sky as it avoids regions at low Galactic latitude and low-declination fields (i.e., between declination 0° and +15°), where the sensitivity of LoTSS is generally a factor of 2−3 lower than the nominal survey noise9. Compared to LoTSS-DR1, it has a sky coverage ∼13 times larger (5634 deg2 vs. 424 deg2) and it samples a broader declination and right ascension range. This makes the new data release more representative of the quality of LoTSS for extragalactic studies. We thus use our findings to refine the predictions on the number of PSZ2 sources that will be found to contain relics and halos at the completion of LoTSS.
We consider that there are 835 detections in the PSZ2 catalog above a declination of 0° and assume a uniform sensitivity for LoTSS for simplicity. We note that the presence of the Galactic plane is taken into account in the PSZ2 selection due to the lack of Planck detections in the zone of avoidance. The results found in LoTSS-DR2 indicate that we will find 251 ± 92 clusters that host a halo and 83 ± 50 clusters that host one or more relics (including candidates), and, as in our study, approximately half of them should be new discoveries. When extrapolating the number of halos, we considered the fraction of 30 ± 11%, which takes into account also the ten halos reported with an asterisk, as in these systems the visual inspection and classification tree led to the class of RH or cRH. If we assume that also all uncertain sources trace a halo, we obtain a conservative upper limit on the number of clusters hosting a radio halo that LoTSS should find in PSZ2 clusters at its completion, namely < 401 ± 117. By considering the number of relic sources found in LoTSS-DR2, we instead predict that at the completion of the survey LoTSS will have observed 109 ± 58 radio relics in PSZ2 clusters.
In the context of the turbulent re-acceleration models for radio halos, Cassano et al. (2010b, 2012) predicted the observation of 350−450 radio halos at z ≤ 0.6 in LoTSS, while the number of radio relics expected to be observed in LoTSS according to Nuza et al. (2012) is ∼2500. As these model expectations were not tailored for PSZ2 clusters, the comparison between the number of halos and relics at completion of LoTSS and the extrapolations obtained from our analysis of the LoTSS-DR2/PSZ2 sample is not straightforward. In this respect, the thoroughly statistical analysis of the results obtained for the PSZ2 clusters in LoTSS-DR2 and implications on our understanding of halos and relics will be presented in the forthcoming papers of the series, where expectations will be further refined by considering: the sensitivity of present observations, the PSZ2 selection functions, and the new P150−M500 correlations. Moreover, as radio halos and relics can be found in LoTSS even in non-PSZ2 clusters, a census of diffuse radio emission in non-PSZ2 clusters in LoTSS-DR2 is also currently ongoing (Hoang et al. 2022).
An important parameter to test the models of particle acceleration in the ICM is the spectrum of diffuse emission. For example, turbulent re-acceleration models predict the existence of a population of radio halos with steep spectra while spectral gradients in relics depend on the mechanisms of particle (re)acceleration at cluster shocks. Interestingly, we note that about half of the radio halos reported in this work are new discoveries (see Sect. 6.1) and the turbulent re-acceleration scenario predicts that about half of the radio halos in LoTSS should have ultra-steep spectra (Cassano et al. 2012). If follow up observations of the radio halos in LoTSS-DR2 will confirm that about 50% of them have very steep spectra, this would be an important corroboration of such a class of models. Still, since LoTSS is currently exploring an uncharted territory in terms of resolution and sensitivity compared to other completed surveys, it is not yet possible to perform spectral analysis for the full cluster sample. In this respect, we have planned a series of targeted observations for a selected number of objects with MeerKAT, VLA, and uGMRT. For a systematic study, ongoing and future sensitive radio surveys covering the northern sky, such as LoLSS (de Gasperin et al. 2021), LoDeSS10, APERture Tile In Focus (Apertif; Hess et al., in prep.), and VLA Sky Survey (VLASS; Lacy et al. 2020), will supplement LoTSS providing a multifrequency view of clusters that will be critical to investigate the statistics of the spectral properties of halos and relics as well as to properly disentangle diffuse sources from the emission of AGN, allowing us to determine the role of seed relativistic electrons in the ICM.
The presence of X-ray data is critical for the classification of diffuse sources in the ICM, and therefore the ongoing eROSITA All-Sky Survey (eRASS; Merloni et al. 2012; Predehl et al. 2021) will be fundamental for providing the X-ray detection of many clusters that currently do not have pointed XMM-Newton or Chandra observations (i.e., about half of the targets in our sample). This instrument will also allow the detection of low-mass and/or high-z clusters (Pillepich et al. 2012) that are currently not confirmed or missed by Planck due to its low sensitivity. Indeed, combining X-ray and SZ surveys is particularly efficient to find new galaxy clusters and derive their mass (see, e.g., ComPRASS; Tarrío et al. 2019). The potential of the joint analysis of LOFAR and eROSITA observations of galaxy clusters has been demonstrated by recent papers (e.g., Ghirardini et al. 2021; Pasini et al. 2022; Brienza et al. 2021). Until the eRASS data are publicly released, we are planning to follow-up a subsample of clusters with Chandra and/or XMM-Newton.
Finally, we note that there are opportunities to further improve the image quality with respect to that presented in this paper. This could be achieved by reapplying the extraction + self-calibration step and fine-tuning the parameters of the self-calibration (van Weeren et al. 2021). The subtraction of discrete sources can also be more careful than that obtained during the analysis of this sample. Moreover, the addition of LOFAR international baselines will help to improve the calibration of the targets affected by the sidelobes of a central bright AGN, enabling the search for diffuse emission in regions that cannot be thoroughly examined with the current calibration. In this respect, we note that a number of targets presented in this sample will be the subject of focused studies in the future (see Appendix D).
7. Conclusions
We have presented the largest statistical sample of galaxy clusters used to date to search for and study diffuse synchrotron sources in the ICM. We examined the 120−168 MHz radio emission from 309 galaxy clusters selected from PSZ2 that span a redshift and mass range of 0.016 < z < 0.9 and 1.1 × 1014 M⊙ < M500 < 11.7 × 1014 M⊙, respectively, and have been covered by LoTSS-DR2. We produced radio images with different resolutions and with or without the discrete source subtracted as well as overlays with Pan-STARRS optical images for all the targets in our sample. When available, we also used targeted Chandra and XMM-Newton observations to compare the radio emission with that of the X-rays and derive morphological parameters of the ICM. All these images are publicly available on the project website11.
We divided the diffuse synchrotron emission into the classes of halos, relics, and uncertain sources. The physical properties of halos and relics have been collected into tables (also available on the website) and are used in Cassano et al. (in prep.), Cuciti et al. (in prep.), and Jones et al. (in prep.) to discuss their statistical properties, such as occurrence and scaling relations. Overall, we found 83 clusters that host a radio halo and 26 clusters that host one or more radio relic (including candidates), of which about half are new discoveries. These numbers correspond to a detection fraction in our sample of 30 ± 11% and 10 ± 6%, respectively. Based on these results, we expect to find 251 ± 92 cluster that host a halo and 83 ± 50 clusters that host at least one relic in the PSZ2 catalog at the completion of LoTSS. Other searches are being made to examine different cluster samples and try to gauge how many more halos and relics can be found in LoTSS in non-PSZ2 clusters (Hoang et al. 2022).
In the future, LoTSS will benefit from the synergy of complementary radio surveys (e.g., LoLSS, LoDeSS, and Apertif), which will be fundamental for studying the spectral properties of the observed sources. The all-sky survey with eROSITA will enable a systematic comparison of the radio and X-ray properties of the PSZ2 clusters as well as the discoveries of new galaxy clusters (especially with low mass and at high z).
The sensitivity of LoTSS observations scales with the elevation of the target as A × cos(90° −elevation)−2, where A = 63 μJy beam−1 and the dependence on elevation is fixed according to the projected size of the LOFAR stations (Shimwell et al. 2022).
Acknowledgments
We thank the anonymous referee for constructive comments that helped to improve the presentation of our results. A.B. and R.J.v.W. acknowledge support from the VIDI research programme with project number 639.042.729, which is financed by the Netherlands Organisation for Scientific Research (NWO). R.C., F.G., M.R. and G.B. acknowledge support from INAF mainstream project ‘Galaxy Clusters Science with LOFAR’ 1.05.01.86.05. V.C. and G.D.G. acknowledge support from the Alexander von Humboldt Foundation. X.Z. acknowledges support from Chinese Scholarship Council (CSC) and thanks Vittorio Ghirardini for useful discussions. A.S. is supported by the Women In Science Excel (WISE) programme of the NWO, and acknowledges the World Premier Research Center Initiative (WPI) and the Kavli IPMU for the continued hospitality. R.J.v.W. and C.G. acknowledge support from the ERC-Stg ClusterWeb 804208. M.B. and F.d.G. acknowledge support from the Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy – EXC 2121 “Quantum Universe” – 390833306. D.N.H. and A.B. acknowledge support from ERC-Stg DRANOEL 714245. M.J.H. acknowledges from the UK Science and Technology Facilities Council [ST/V000624/1]. A.I. acknowledges the Italian PRIN-Miur 2017 (PI A. Cimatti). A.Bon. acknowledges support from MIUR FARE grant “SMS”. A.D. acknowledges support by the BMBF Verbundforschung under the grant 05A20STA. LOFAR (van Haarlem et al. 2013) is the LOw Frequency ARray designed and constructed by ASTRON. It has observing, data processing, and data storage facilities in several countries, which are owned by various parties (each with their own funding sources), and are collectively operated by the ILT foundation under a joint scientific policy. The ILT resources have benefitted from the following recent major funding sources: CNRS-INSU, Observatoire de Paris and Université d’Orléans, France; BMBF, MIWF-NRW, MPG, Germany; Science Foundation Ireland (SFI), Department of Business, Enterprise and Innovation (DBEI), Ireland; NWO, The Netherlands; The Science and Technology Facilities Council, UK; Ministry of Science and Higher Education, Poland; Istituto Nazionale di Astrofisica (INAF), Italy. This research made use of the Dutch national e-infrastructure with support of the SURF Cooperative (e-infra 180169) and the LOFAR e-infra group, and of the LOFAR-IT computing infrastructure supported and operated by INAF, and by the Physics Dept. of Turin University (under the agreement with Consorzio Interuniversitario per la Fisica Spaziale) at the C3S Supercomputing Centre, Italy. The Jülich LOFAR Long Term Archive and the German LOFAR network are both coordinated and operated by the Jülich Supercomputing Centre (JSC), and computing resources on the supercomputer JUWELS at JSC were provided by the Gauss Centre for Supercomputing e.V. (grant CHTB00) through the John von Neumann Institute for Computing (NIC). This research made use of the University of Hertfordshire high-performance computing facility and the LOFAR-UK computing facility located at the University of Hertfordshire and supported by STFC [ST/P000096/1]. The scientific results reported in this article are based in part on data obtained from the Chandra Data Archive and on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. SRON Netherlands Institute for Space Research is supported financially by the NWO. This research made use of APLpy, an open-source plotting package for Python (Robitaille & Bressert 2012), Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018), Matplotlib (Hunter 2007), NumPy (Harris et al. 2020), and SciPy (Virtanen et al. 2020).
References
- Adam, R., Goksu, H., Brown, S., Rudnick, L., & Ferrari, C. 2021, A&A, 648, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Aghanim, N., Hurier, G., Diego, J. M., et al. 2015, A&A, 580, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Aguado-Barahona, A., Barrena, R., Streblyanska, A., et al. 2019, A&A, 631, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Alhassan, W., Taylor, A., & Vaccari, M. 2018, MNRAS, 480, 2085 [NASA ADS] [CrossRef] [Google Scholar]
- Aniyan, A., & Thorat, K. 2017, ApJS, 230, 20 [NASA ADS] [CrossRef] [Google Scholar]
- Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
- Bacchi, M., Feretti, L., Giovannini, G., & Govoni, F. 2003, A&A, 400, 465 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barrena, R., Streblyanska, A., Ferragamo, A., et al. 2018, A&A, 616, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Biava, N., de Gasperin, F., Bonafede, A., et al. 2021a, MNRAS, 508, 3995 [NASA ADS] [CrossRef] [Google Scholar]
- Biava, N., Brienza, M., Bonafede, A., et al. 2021b, A&A, 650, A170 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bîrzan, L., Rafferty, D., Cassano, R., et al. 2019, MNRAS, 487, 4775 [CrossRef] [Google Scholar]
- Bîrzan, L., Rafferty, D., Brüggen, M., et al. 2020, MNRAS, 496, 2613 [CrossRef] [Google Scholar]
- Blasi, P., & Colafrancesco, S. 1999, Astropart. Phys., 12, 169 [Google Scholar]
- Boada, S., Hughes, J. P., Menanteau, F., et al. 2019, ApJ, 871, 188 [NASA ADS] [CrossRef] [Google Scholar]
- Bonafede, A., Giovannini, G., Feretti, L., Govoni, F., & Murgia, M. 2009, A&A, 494, 429 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bonafede, A., Feretti, L., Murgia, M., et al. 2010, A&A, 513, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bonafede, A., Brüggen, M., van Weeren, R. J., et al. 2012, MNRAS, 426, 40 [CrossRef] [Google Scholar]
- Bonafede, A., Brunetti, G., Vazza, F., et al. 2021, ApJ, 907, 32 [Google Scholar]
- Bonafede, A., Brunetti, G., Rudnick, L., et al. 2022, ApJ, submitted [arXiv:2203.01958] [Google Scholar]
- Botteon, A., Gastaldello, F., Brunetti, G., & Dallacasa, D. 2016, MNRAS, 460, L84 [Google Scholar]
- Botteon, A., Shimwell, T. W., Bonafede, A., et al. 2018, MNRAS, 478, 885 [Google Scholar]
- Botteon, A., Cassano, R., Eckert, D., et al. 2019a, A&A, 630, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Botteon, A., Shimwell, T. W., Bonafede, A., et al. 2019b, A&A, 622, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Botteon, A., van Weeren, R. J., Brunetti, G., et al. 2020a, MNRAS, 499, L11 [Google Scholar]
- Botteon, A., Brunetti, G., van Weeren, R. J., et al. 2020b, ApJ, 897, 93 [Google Scholar]
- Botteon, A., Cassano, R., van Weeren, R. J., et al. 2021a, ApJ, 914, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Botteon, A., Giacintucci, S., Gastaldello, F., et al. 2021b, A&A, 649, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Boxelaar, J. M., van Weeren, R. J., & Botteon, A. 2021, Astron. Comput., 35, 100464 [NASA ADS] [CrossRef] [Google Scholar]
- Brienza, M., Morganti, R., Harwood, J. J., et al. 2020, A&A, 638, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brienza, M., Shimwell, T. W., de Gasperin, F., et al. 2021, Nat. Astron., 5, 1261 [NASA ADS] [CrossRef] [Google Scholar]
- Briggs, D. 1995, BAAS, 27, 1444 [NASA ADS] [Google Scholar]
- Brown, S., & Rudnick, L. 2011, MNRAS, 412, 2 [Google Scholar]
- Brüggen, M., & Vazza, F. 2020, MNRAS, 493, 2306 [Google Scholar]
- Brüggen, M., Rafferty, D., Bonafede, A., et al. 2018, MNRAS, 477, 3461 [CrossRef] [Google Scholar]
- Brunetti, G., & Jones, T. W. 2014, Int. J. Mod. Phys. D, 23, 30007 [Google Scholar]
- Brunetti, G., Cassano, R., Dolag, K., & Setti, G. 2009, A&A, 507, 661 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brunetti, G., Zimmer, S., & Zandanel, F. 2017, MNRAS, 472, 1506 [Google Scholar]
- Buddendiek, A., Schrabback, T., Greer, C., et al. 2015, MNRAS, 450, 4248 [NASA ADS] [CrossRef] [Google Scholar]
- Buote, D. A. 2001, ApJ, 553, L15 [Google Scholar]
- Burenin, R. 2017, Astron. Lett., 43, 507 [NASA ADS] [CrossRef] [Google Scholar]
- Burenin, R., Bikmaev, I., Khamitov, I., et al. 2018, Astron. Lett., 44, 297 [NASA ADS] [CrossRef] [Google Scholar]
- Cassano, R., & Brunetti, G. 2005, MNRAS, 357, 1313 [Google Scholar]
- Cassano, R., Brunetti, G., & Setti, G. 2006, MNRAS, 369, 1577 [Google Scholar]
- Cassano, R., Brunetti, G., Setti, G., Govoni, F., & Dolag, K. 2007, MNRAS, 378, 1565 [Google Scholar]
- Cassano, R., Brunetti, G., Venturi, T., et al. 2008, A&A, 480, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cassano, R., Ettori, S., Giacintucci, S., et al. 2010a, ApJ, 721, L82 [Google Scholar]
- Cassano, R., Brunetti, G., Röttgering, H. J., & Brüggen, M. 2010b, A&A, 509, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cassano, R., Brunetti, G., Norris, R. P., et al. 2012, A&A, 548, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cassano, R., Botteon, A., Di Gennaro, G., et al. 2019, ApJ, 881, L18 [Google Scholar]
- Chambers, K. C., Magnier, E., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560] [Google Scholar]
- Clarke, T. E., Randall, S. W., Sarazin, C. L., Blanton, E. L., & Giacintucci, S. 2013, ApJ, 772, 84 [NASA ADS] [CrossRef] [Google Scholar]
- Clarke, A., Scaife, A., Shimwell, T. W., et al. 2019, A&A, 627, A176 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cohen, A., & Clarke, T. E. 2011, AJ, 141, 149 [Google Scholar]
- Condon, J., Cotton, W., Greisen, E., et al. 1998, AJ, 115, 1693 [NASA ADS] [CrossRef] [Google Scholar]
- Cuciti, V., Cassano, R., Brunetti, G., et al. 2015, A&A, 580, A97 [CrossRef] [EDP Sciences] [Google Scholar]
- Cuciti, V., Cassano, R., Brunetti, G., et al. 2021, A&A, 647, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Gasperin, F., van Weeren, R. J., Brüggen, M., et al. 2014, MNRAS, 444, 3130 [NASA ADS] [CrossRef] [Google Scholar]
- de Gasperin, F., Ogrean, G., van Weeren, R. J., et al. 2015, MNRAS, 448, 2197 [NASA ADS] [CrossRef] [Google Scholar]
- de Gasperin, F., Intema, H. T., Shimwell, T. W., et al. 2017, Sci. Adv., 3, e1701634 [Google Scholar]
- de Gasperin, F., Dijkema, T. J., Drabent, A., et al. 2019, A&A, 622, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Gasperin, F., Williams, W. L., Best, P. N., et al. 2021, A&A, 648, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dennison, B. 1980, ApJ, 239, L93 [Google Scholar]
- Di Gennaro, G., van Weeren, R. J., Brunetti, G., et al. 2021a, Nat. Astron., 5, 268 [Google Scholar]
- Di Gennaro, G., van Weeren, R. J., Cassano, R., et al. 2021b, A&A, 654, A166 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Domínguez Sánchez, H., Huertas-Company, M., Bernardi, M., Tuccillo, D., & Fischer, J. 2018, MNRAS, 476, 3661 [CrossRef] [Google Scholar]
- Drabent, A., Hoeft, M., Pizzo, R. F., et al. 2015, A&A, 575, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Duchesne, S. W., Johnston-Hollitt, M., Zhu, Z.-H., Wayth, R., & Line, J. 2020, PASA, 37, e037 [NASA ADS] [CrossRef] [Google Scholar]
- Duchesne, S. W., Johnston-Hollitt, M., & Bartalucci, I. 2021, PASA, 38, e053 [CrossRef] [Google Scholar]
- Duchesne, S., Johnston-Hollitt, M., Riseley, C. J., Bartalucci, I., & Keel, S. R. 2022, MNRAS, 511, 3525 [NASA ADS] [CrossRef] [Google Scholar]
- Farnsworth, D., Rudnick, L., Brown, S., & Brunetti, G. 2013, ApJ, 779, 189 [Google Scholar]
- Feretti, L., & Neumann, D. 2006, A&A, 450, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Feretti, L., Orrù, E., Brunetti, G., et al. 2004, A&A, 423, 111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Feretti, L., Giovannini, G., Govoni, F., & Murgia, M. 2012, A&ARv, 20, 54 [Google Scholar]
- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
- Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, SPIE Conf. Ser., 6270, 62701V [Google Scholar]
- Gabriel, C., Denby, M., Fyfe, D., et al. 2004, in Astronomical Data Analysis Software and Systems XIII, eds. F. Ochsenbein, M. G. Allen, & D. Egret, ASP Conf. Ser., 314, 759 [NASA ADS] [Google Scholar]
- George, L. T., Kale, R., & Wadadekar, Y. 2021, MNRAS, 507, 4487 [NASA ADS] [CrossRef] [Google Scholar]
- Gheller, C., Vazza, F., & Bonafede, A. 2018, MNRAS, 480, 3749 [NASA ADS] [CrossRef] [Google Scholar]
- Ghirardini, V., Bulbul, E., Hoang, D. N., et al. 2021, A&A, 647, A4 [EDP Sciences] [Google Scholar]
- Giacintucci, S., Venturi, T., Murgia, M., et al. 2007, A&A, 476, 99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Giacintucci, S., Venturi, T., Cassano, R., Dallacasa, D., & Brunetti, G. 2009, ApJ, 704, L54 [NASA ADS] [CrossRef] [Google Scholar]
- Giacintucci, S., Markevitch, M., Venturi, T., et al. 2014, ApJ, 781, 9 [Google Scholar]
- Giacintucci, S., Markevitch, M., Cassano, R., et al. 2017, ApJ, 841, 71 [Google Scholar]
- Giovannini, G., & Feretti, L. 2000, New Astron., 5, 335 [Google Scholar]
- Giovannini, G., Feretti, L., & Stanghellini, C. 1991, A&A, 252, 528 [NASA ADS] [Google Scholar]
- Giovannini, G., Feretti, L., Venturi, T., Kim, K.-T., & Kronberg, P. P. 1993, ApJ, 406, 399 [NASA ADS] [CrossRef] [Google Scholar]
- Giovannini, G., Tordi, M., & Feretti, L. 1999, New Astron., 4, 141 [CrossRef] [Google Scholar]
- Giovannini, G., Feretti, L., Govoni, F., Murgia, M., & Pizzo, R. F. 2006, Astron. Nachr., 327, 563 [NASA ADS] [CrossRef] [Google Scholar]
- Giovannini, G., Bonafede, A., Feretti, L., et al. 2009, A&A, 507, 1257 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gitti, M. 2013, MNRAS, 436, L84 [NASA ADS] [CrossRef] [Google Scholar]
- Gitti, M., Brunetti, G., Feretti, L., & Setti, G. 2004, A&A, 417, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Govoni, F., Feretti, L., Giovannini, G., et al. 2001, A&A, 376, 803 [CrossRef] [EDP Sciences] [Google Scholar]
- Govoni, F., Murgia, M., Markevitch, M., et al. 2009, A&A, 499, 371 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Govoni, F., Murgia, M., Giovannini, G., Vacca, V., & Bonafede, A. 2011, A&A, 529, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Govoni, F., Ferrari, C., Feretti, L., et al. 2012, A&A, 545, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Govoni, F., Orrù, E., Bonafede, A., et al. 2019, Science, 364, 981 [Google Scholar]
- Gupta, Y., Ajithkumar, B., Kale, H., et al. 2017, Curr. Sci., 113, 707 [NASA ADS] [CrossRef] [Google Scholar]
- Hales, S., Baldwin, J., & Warner, P. 1988, MNRAS, 234, 919 [NASA ADS] [CrossRef] [Google Scholar]
- Hales, S., Masson, C., Warner, P., & Baldwin, J. 1990, MNRAS, 246, 256 [NASA ADS] [Google Scholar]
- Hallman, E. J., Alden, B., Rapetti, D., Datta, A., & Burns, J. O. 2018, ApJ, 859, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Hardcastle, M. J., Gürkan, G., van Weeren, R. J., et al. 2016, MNRAS, 462, 1910 [Google Scholar]
- Hardcastle, M. J., Croston, J., Shimwell, T. W., et al. 2019, MNRAS, 488, 3416 [Google Scholar]
- Hardcastle, M. J., Shimwell, T. W., Tasse, C., et al. 2021, A&A, 648, A10 [EDP Sciences] [Google Scholar]
- Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [NASA ADS] [CrossRef] [Google Scholar]
- Hlavacek-Larrondo, J., Gendron-Marsolais, M.-L., Fecteau-Beaucage, D., et al. 2018, MNRAS, 475, 2743 [NASA ADS] [CrossRef] [Google Scholar]
- Hoang, D. N., Shimwell, T. W., van Weeren, R. J., et al. 2018, MNRAS, 478, 2218 [Google Scholar]
- Hoang, D. N., Shimwell, T. W., van Weeren, R. J., et al. 2019, A&A, 622, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hoang, D. N., Zhang, X., Stuardi, C., et al. 2021a, A&A, 656, A154 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hoang, D. N., Shimwell, T. W., Osinga, E., et al. 2021b, MNRAS, 501, 576 [Google Scholar]
- Hoang, D. N., Brüggen, M., Botteon, A., et al. 2022, A&A, submitted [Google Scholar]
- Hodgson, T., Bartalucci, I., Johnston-Hollitt, M., et al. 2021, ApJ, 909, 198 [CrossRef] [Google Scholar]
- Hoeft, M., Dumba, C., Drabent, A., et al. 2021, A&A, 654, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hogbom, J., & Brouw, W. 1974, A&A, 33, 289 [NASA ADS] [Google Scholar]
- Hotan, A., Bunton, J., Chippendale, A., et al. 2021, PASA, 38, e009 [NASA ADS] [CrossRef] [Google Scholar]
- Huber, P. J. 1981, Robust Statistics (New York: John Wiley and Sons) [CrossRef] [Google Scholar]
- Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
- HyeongHan, K., Jee, M. J., Rudnick, L., et al. 2020, ApJ, 900, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Ignesti, A., Gitti, M., Brunetti, G., Feretti, L., & Giovannini, G. 2017, A&A, 604, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ignesti, A., Shimwell, T. W., Brunetti, G., et al. 2020, A&A, 643, A172 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ignesti, A., Brunetti, G., Shimwell, T. W., et al. 2022a, A&A, 659, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ignesti, A., Vulcani, B., Poggianti, B. M., et al. 2022b, ApJ, 924, 64 [NASA ADS] [CrossRef] [Google Scholar]
- Jeltema, T. E., & Profumo, S. 2011, ApJ, 728, 53 [Google Scholar]
- Jonas, J. 2009, IEEE Proc., 97, 1522 [NASA ADS] [CrossRef] [Google Scholar]
- Kale, R., & Gitti, M. 2017, MNRAS, 466, L19 [NASA ADS] [CrossRef] [Google Scholar]
- Kale, R., Venturi, T., Giacintucci, S., et al. 2013, A&A, 557, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kale, R., Parekh, V., & Dwarakanath, K. 2018, MNRAS, 480, 5352 [NASA ADS] [CrossRef] [Google Scholar]
- Kempner, J. C., & Sarazin, C. L. 2001, ApJ, 548, 639 [Google Scholar]
- Kim, K.-T., Kronberg, P. P., Giovannini, G., & Venturi, T. 1989, Nature, 341, 720 [NASA ADS] [CrossRef] [Google Scholar]
- Kokotanekov, G., Wise, M. W., de Vries, M., & Intema, H. T. 2018, A&A, 618, A152 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lacy, M., Baum, S. A., Chandler, C., et al. 2020, PASP, 132, 35001 [Google Scholar]
- Liang, H., Hunstead, R. W., Birkinshaw, M., & Andreani, P. 2000, ApJ, 544, 686 [Google Scholar]
- Locatelli, N., Rajpurohit, K., Vazza, F., et al. 2020, MNRAS, 496, L48 [NASA ADS] [CrossRef] [Google Scholar]
- Lovisari, L., Forman, W. R., Jones, C., et al. 2017, ApJ, 846, 51 [Google Scholar]
- Lukic, V., Brüggen, M., Banfield, J., et al. 2018, MNRAS, 476, 246 [NASA ADS] [CrossRef] [Google Scholar]
- Lukic, V., Brüggen, M., Mingo, B., et al. 2019, MNRAS, 487, 1729 [Google Scholar]
- Macario, G., Venturi, T., Brunetti, G., et al. 2010, A&A, 517, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Macario, G., Venturi, T., Intema, H. T., et al. 2013, A&A, 551, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mandal, S., Intema, H. T., Shimwell, T. W., et al. 2019, A&A, 622, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mandal, S., Intema, H. T., van Weeren, R. J., et al. 2020, A&A, 634, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Markevitch, M., & Vikhlinin, A. 2007, Phys. Rep., 443, 1 [Google Scholar]
- Merloni, A., Predehl, P., Becker, W., et al. 2012, ArXiv e-prints [arXiv:1209.3114] [Google Scholar]
- Moffet, A., & Birkinshaw, M. 1989, AJ, 98, 1148 [NASA ADS] [CrossRef] [Google Scholar]
- Mohr, J. J., Fabricant, D. G., & Geller, M. J. 1993, ApJ, 413, 492 [Google Scholar]
- Mostert, R. I., Duncan, K. J., Röttgering, H. J., et al. 2021, A&A, 645, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Murgia, M., Govoni, F., Markevitch, M., et al. 2009, A&A, 499, 679 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nurgaliev, D., McDonald, M., Benson, B., et al. 2013, ApJ, 779, 112 [NASA ADS] [CrossRef] [Google Scholar]
- Nuza, S., Hoeft, M., van Weeren, R. J., Gottlöber, S., & Yepes, G. 2012, MNRAS, 420, 2006 [NASA ADS] [CrossRef] [Google Scholar]
- Nuza, S., Gelszinnis, J., Hoeft, M., & Yepes, G. 2017, MNRAS, 470, 240 [NASA ADS] [CrossRef] [Google Scholar]
- Offringa, A. R., & Smirnov, O. M. 2017, MNRAS, 471, 301 [Google Scholar]
- Offringa, A. R., McKinley, B., Hurley-Walker, N., et al. 2014, MNRAS, 444, 606 [Google Scholar]
- Orrù, E., Murgia, M., Feretti, L., et al. 2007, A&A, 467, 943 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Osinga, E., van Weeren, R. J., Boxelaar, J. M., et al. 2021, A&A, 648, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Owen, F. N., & Ledlow, M. J. 1997, ApJS, 108, 41 [Google Scholar]
- Owen, F. N., Morrison, G. E., & Voges, W. 1999, in Diffuse Thermal and Relativistic Plasma in Galaxy Clusters, eds. H. Boehringer, L. Feretti, & P. Schuecker, 9 [Google Scholar]
- Parekh, V., van der Heyden, K., Ferrari, C., Angus, G., & Holwerda, B. 2015, A&A, 575, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pasini, T., Brüggen, M., Hoang, D. N., et al. 2022, A&A, in press, https://doi.org/10.1051/0004-6361/202141211 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Paul, S., Salunkhe, S., Sonkamble, S., et al. 2020, A&A, 633, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Paul, S., Gupta, P., Salunkhe, S., et al. 2021, MNRAS, 506, 5389 [NASA ADS] [CrossRef] [Google Scholar]
- Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2012, ArXiv e-prints [arXiv:1201.0490] [Google Scholar]
- Perley, R. A., Chandler, C., Butler, B., & Wrobel, J. 2011, ApJ, 739, L1 [NASA ADS] [CrossRef] [Google Scholar]
- Pillepich, A., Porciani, C., & Reiprich, T. H. 2012, MNRAS, 422, 44 [Google Scholar]
- Planck Collaboration XXII. 2016, A&A, 594, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XXVII. 2016, A&A, 594, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Poole, G., Fardal, M., Babul, A., et al. 2006, MNRAS, 373, 881 [Google Scholar]
- Predehl, P., Andritschke, R., Arefiev, V., et al. 2021, A&A, 647, A1 [EDP Sciences] [Google Scholar]
- Rasia, E., Meneghetti, M., & Ettori, S. 2013, Astron. Rev., 8, 40 [Google Scholar]
- Roberts, I. D., van Weeren, R. J., McGee, S. L., et al. 2021a, A&A, 650, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Roberts, I. D., van Weeren, R. J., McGee, S. L., et al. 2021b, A&A, 652, A153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Roberts, I. D., van Weeren, R. J., Timmerman, R., et al. 2022, A&A, 658, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Robitaille, T. P., & Bressert, E. 2012, Astrophysics Source Code Library [record ascl:1208.017] [Google Scholar]
- Roger, R., Costain, C., & Bridle, A. 1973, AJ, 78, 1030 [NASA ADS] [CrossRef] [Google Scholar]
- Rossetti, M., Gastaldello, F., Eckert, D., et al. 2017, MNRAS, 468, 1917 [Google Scholar]
- Rudnick, L., & Lemmerman, J. 2009, ApJ, 697, 1341 [NASA ADS] [CrossRef] [Google Scholar]
- Sadeghi, M., Javaherian, M., & Miraghaei, H. 2021, AJ, 161, 94 [NASA ADS] [CrossRef] [Google Scholar]
- Santos, J., Rosati, P., Tozzi, P., et al. 2008, A&A, 483, 35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Savini, F., Bonafede, A., Brüggen, M., et al. 2018a, MNRAS, 478, 2234 [Google Scholar]
- Savini, F., Bonafede, A., Brüggen, M., et al. 2018b, MNRAS, 474, 5023 [NASA ADS] [CrossRef] [Google Scholar]
- Savini, F., Bonafede, A., Brüggen, M., et al. 2019, A&A, 622, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schellenberger, G., David, L. P., O’Sullivan, E., Vrtilek, J. M., & Haines, C. P. 2019, ApJ, 882, 59 [NASA ADS] [CrossRef] [Google Scholar]
- Sehgal, N., Hughes, J. P., Wittman, D., et al. 2008, ApJ, 673, 163 [CrossRef] [Google Scholar]
- Sen, P. K. 1968, J. Am. Stat. Assoc., 63, 1379 [CrossRef] [Google Scholar]
- Shimwell, T. W., Luckin, J., Brüggen, M., et al. 2016, MNRAS, 459, 277 [Google Scholar]
- Shimwell, T. W., Röttgering, H. J., Best, P. N., et al. 2017, A&A, 598, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shimwell, T. W., Tasse, C., Hardcastle, M. J., et al. 2019, A&A, 622, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shimwell, T. W., Tasse, C., Hardcastle, M. J., et al. 2022, A&A, 659, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Smirnov, O. M., & Tasse, C. 2015, MNRAS, 449, 2668 [Google Scholar]
- Sommer, M. W., Basu, K., Intema, H. T., et al. 2017, MNRAS, 466, 996 [NASA ADS] [CrossRef] [Google Scholar]
- Streblyanska, A., Barrena, R., Rubiño-Martín, J., et al. 2018, A&A, 617, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Streblyanska, A., Aguado-Barahona, A., Ferragamo, A., et al. 2019, A&A, 628, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sunyaev, R., & Zel’dovich, Y. 1972, Comments Astrophys. Space Phys., 4, 173 [Google Scholar]
- Swarup, G., Ananthakrishnan, S., Kapahi, V., et al. 1991, Curr. Sci., 60, 95 [Google Scholar]
- Tarrío, P., Melin, J.-B., & Arnaud, M. 2019, A&A, 626, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tasse, C. 2014a, A&A, 566, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tasse, C. 2014b, ArXiv e-prints [arXiv:1410.8706] [Google Scholar]
- Tasse, C., Hugo, B., Mirmont, M., et al. 2018, A&A, 611, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tasse, C., Shimwell, T. W., Hardcastle, M. J., et al. 2021, A&A, 648, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Terni de Gregory, B., Feretti, L., Giovannini, G., et al. 2017, A&A, 608, A58 [EDP Sciences] [Google Scholar]
- Thierbach, M., Klein, U., & Wielebinski, R. 2003, A&A, 397, 53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Thompson, A., Clark, B., Wade, C., & Napier, P. 1980, ApJS, 44, 151 [NASA ADS] [CrossRef] [Google Scholar]
- Timmerman, R., van Weeren, R. J., Callingham, J., et al. 2022, A&A, 658, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tingay, S. J., Goeke, R., Bowman, J. D., et al. 2013, PASA, 30, e007 [Google Scholar]
- Vacca, V., Murgia, M., Govoni, F., et al. 2010, A&A, 514, A71 [CrossRef] [EDP Sciences] [Google Scholar]
- van Haarlem, M., Wise, M. W., Gunst, A., et al. 2013, A&A, 556, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Weeren, R. J., Röttgering, H. J., Brüggen, M., & Cohen, A. 2009, A&A, 508, 75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Weeren, R. J., Brüggen, M., Röttgering, H. J., et al. 2011, A&A, 533, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Weeren, R. J., Bonafede, A., Ebeling, H., et al. 2012, MNRAS, 425, L36 [NASA ADS] [CrossRef] [Google Scholar]
- van Weeren, R. J., Williams, W. L., Hardcastle, M. J., et al. 2016, ApJS, 223, 2 [Google Scholar]
- van Weeren, R. J., de Gasperin, F., Akamatsu, H., et al. 2019, Space Sci. Rev., 215, 16 [Google Scholar]
- van Weeren, R. J., Shimwell, T. W., Botteon, A., et al. 2021, A&A, 651, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vavilova, I., Dobrycheva, D., Vasylenko, M. Y., et al. 2021, A&A, 648, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Venturi, T., Giacintucci, S., Brunetti, G., et al. 2007, A&A, 463, 937 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Venturi, T., Giacintucci, S., Dallacasa, D., et al. 2008, A&A, 484, 327 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Venturi, T., Giacintucci, S., Dallacasa, D., et al. 2011, MNRAS, 414, L65 [NASA ADS] [CrossRef] [Google Scholar]
- Venturi, T., Giacintucci, S., Dallacasa, D., et al. 2013, A&A, 551, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nat. Meth., 17, 261 [Google Scholar]
- Wen, Z., & Han, J. 2013, MNRAS, 436, 275 [NASA ADS] [CrossRef] [Google Scholar]
- Wilber, A. G., Brüggen, M., Bonafede, A., et al. 2018, MNRAS, 473, 3536 [Google Scholar]
- Wilber, A. G., Brüggen, M., Bonafede, A., et al. 2019, A&A, 622, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Williams, W. L., van Weeren, R. J., Röttgering, H. J., et al. 2016, MNRAS, 460, 2385 [Google Scholar]
- Yuan, Z., & Han, J. 2020, MNRAS, 497, 5485 [NASA ADS] [CrossRef] [Google Scholar]
- Yuan, Z., Han, J., & Wen, Z. 2015, ApJ, 813, 77 [NASA ADS] [CrossRef] [Google Scholar]
- Zandanel, F., & Ando, S. 2014, MNRAS, 440, 663 [Google Scholar]
- Zohren, H., Schrabback, T., van der Burg, R. F., et al. 2019, MNRAS, 488, 2523 [CrossRef] [Google Scholar]
Appendix A: Tables
In this appendix we collect the tables listing: the main properties of the targets in the full PSZ2/LoTSS-DR2 sample (Table A.1), the X-ray morphological parameters for the clusters observed with Chandra/XMM-Newton (Table A.2), and the quantities derived for the radio halos (Table A.3) and radio relics (Table A.4) observed. All tables are available in FITS format on the project website https://lofar-surveys.org/planck_dr2.html.
Full sample of PSZ2 clusters in LoTSS-DR2.
Subsample of galaxy clusters with Chandra/XMM-Newton data available for which we computed the X-ray morphological parameters.
Sample of radio halos and candidate radio halos found in this work.
Sample of radio relics and candidate radio relics found in this work.
Appendix B: Image gallery
Figure B.1 shows a collection of LOFAR images of the targets in our LoTSS-DR2/PSZ2 sample.
![]() |
Fig. B.1. LOFAR image gallery of the targets in our sample, excluding the five objects for which it was not possible to apply the extraction + recalibration scheme (see end of Sect. 3.1). The collage is available at full resolution on the project website, https://lofar-surveys.org/planck_dr2.html. |
![]() |
Fig. B.1. continued. |
![]() |
Fig. B.1. continued. |
![]() |
Fig. B.1. continued. |
Appendix C: Comparison between flux density measurements for radio halos
Here we compare the flux density for the radio halos obtained by fitting their surface brightness profile with HALO-FDCA (Sfit) and by manual integration of a circular/elliptical region encompassing the 2σ level contour (S2σ). In both cases, regions affected by residual calibration errors and/or contaminating extended sources were masked out and the flux density beneath them was extrapolated assuming the best-fit model (in the case of HALO-FDCA) or the average surface brightness of the halo in the non-masked region (in the case of the manual measurement). The left panel of Fig. C.1 shows that there is a general agreement between the two measurements, although a few outliers are present.
![]() |
Fig. C.1. Flux density measurements for radio halos. Left: Comparison between the flux density derived with HALO-FDCA (Sfit) and that obtained manually within the 2σ contour (S2σ). The dashed line indicates the linear correlation as a reference. Right: Diagnostic plot used to determine when HALO-FDCA does not provide a reliable flux density measurement due to the low significance of the radio halo. The dashed lines indicate the thresholds of Sfit/S2σ > 1.5 and I0/3σ < 2 that were used to identify the sources that were classified as radio halos or candidate radio halos with the decision tree of Fig. 4, but whose flux density was not considered reliable. These sources are reported in red and in the manuscript are referred to as RH* and cRH*. |
We inspected the radio images and results of the fitting for all clusters and found that the outliers may indicate cases in which HALO-FDCA provides a unreliable fit because of the low significance of the radio halo. To verify this, we produce the plot shown in the right panel of Fig. C.1, which shows the ratio between the best-fit value of the central surface brightness I0 and the 3σ noise of the map as a function of the Sfit/S2σ ratio. We observe a trend between the plotted quantities, noting indeed that the halos where the two flux density measurements are in most disagreement (large Sfit/S2σ values) are typically those with lowest I0/3σ. Our interpretation is that when the radio halo is observed with low significance, HALO-FDCA does not converge on the diffuse emission but it provides a nonphysical best fit with a low I0 and a large r1, basically fitting the image noise. Importantly, these nonphysical best-fit models have still good S/N and values and corner plots as the fit is still converged; thus, they cannot be simply rejected by using these tools. We therefore used the right panel of Fig. C.1 as diagnostic plot, adopting the arbitrary thresholds of Sfit/S2σ > 1.5 and I0/3σ < 2 to identify the bad cases. These radio halos and candidate radio halos are marked with an asterisk in Table A.1 (i.e., RH*/cRH*) and represent 10 out of the 83 fitted halos. They are collected in Table C.1. Deeper LOFAR observations are required to reliably determine the integrated flux density of the halo emission in these clusters.
Sample of radio halos and candidate radio halos that were classified with an asterisk (the RH* and cRH* in Table A.1) because of their low significance. The integrated flux density of these sources cannot be determined accurately with current data.
Appendix D: Notes on individual clusters
In the following, we make specific remarks for some clusters of our sample.
PSZ2 G023.17+86.71 The radio halo is elongated in the E-W direction and it is claimed in this paper for the first time.
PSZ2 G031.93+78.71 (A1775) It hosts a prominent head-tail radio galaxy that has been studied in the past (Owen & Ledlow 1997; Giovannini & Feretti 2000; Giacintucci et al. 2007; Terni de Gregory et al. 2017). The presence in the cluster center of a radio halo bounded by a cold front was recently claimed by Botteon et al. (2021b) using LoTSS data, who referred to this emission as a “slingshot” radio halo. Toward the SE, a diffuse source elongated along the N-S direction and with uncertain nature is also observed.
PSZ2 G033.81+77.18 (A1795) Past observations with the VLA and GMRT found the presence of diffuse emission in the cluster center, which may be related to a mini-halo or AGN activity (Giacintucci et al. 2014; Kokotanekov et al. 2018). Due to the bad data quality, we cannot perform the analysis of this system. Bîrzan et al. (2020) used LOFAR data to search for an association between the central radio source and the X-ray cavity in the NW, but found none. A paper presenting the analysis of this system using the LOFAR international baselines is currently in preparation (Timmerman et al., in preparation).
PSZ2 G040.58+77.12 (A1800) The radio halo is claimed here for the first time. A roundish region of emission without clear optical counterpart is observed in the S.
PSZ2 G045.87+57.70 The presence of diffuse emission in this cluster was originally claimed by Di Gennaro et al. (2021a) using LoTSS data. We confirm the central emission and following our decision tree we classified it as radio halo.
PSZ2 G046.88+56.48 (A2069) This is a bimodal cluster in which the first detection of megaparsec-scale radio emission was claimed by Farnsworth et al. (2013). Drabent et al. (2015) found extended emission in both components. This is confirmed by our images, in which the radio emission seems to embed the whole system (including the region in-between the two subclusters). The detailed analysis of the LOFAR observations of this target will be presented by Drabent et al. (in preparation).
PSZ2 G048.10+57.16 (A2061) It hosts a radio halo and a relic in the SW direction, which have already been claimed (Rudnick & Lemmerman 2009; van Weeren et al. 2011; Farnsworth et al. 2013). In our images, the radio halo remarkably follows the X-ray emission, which is stretched along the SW-NE direction, and is connected with the radio relic with a trail of emission. The detailed analysis of the LOFAR observations of this target will be presented by Drabent et al. (in preparation).
PSZ2 G049.32+44.37 (A2175) The radio halo is claimed here for the first time.
PSZ2 G053.53+59.52 (A2034) Diffuse radio emission in the center of this cluster was noticed by different authors (Kempner & Sarazin 2001; Rudnick & Lemmerman 2009; Giovannini et al. 2009; van Weeren et al. 2011). Shimwell et al. (2016) provided the clearest view of the diffuse radio sources in this system to date using LOFAR observations. They claimed the presence of a radio halo with filamentary structures and a number of sources in the cluster outskirts with unclear origin. Our new LoTSS images confirm the structures observed by Shimwell et al. (2016). This cluster is displayed in Fig. 5.
PSZ2 G055.59+31.85 (A2261) Hints of diffuse emission in this system were first noticed with VLA 1.4 GHz observations by Venturi et al. (2008), but not confirmed with GMRT 235 and 610 MHz data by Kale et al. (2013). Sommer et al. (2017) reanalyzed these radio observations and claimed to detect the presence of megaparsec-scale emission in the cluster, which was classified as a radio halo. The radio halo was confirmed by LOFAR observations (Savini et al. 2019). Our images agree with these findings.
PSZ2 G056.77+36.32 (A2244) The radio halo is claimed here for the first time.
PSZ2 G057.61+34.93 (A2249) The presence of a radio relic in the E was recently reported with LOFAR and uGMRT data by Locatelli et al. (2020). The relic is confirmed in our image while we classified as uncertain the emission in the cluster center due to the presence of artifacts from the bright central AGN.
PSZ2 G057.80+88.00 (Coma cluster) It hosts the prototypical radio halo and relic, which have been studied in the past with numerous instruments (e.g., Kim et al. 1989; Giovannini et al. 1991, 1993; Thierbach et al. 2003; Feretti & Neumann 2006; Bonafede et al. 2010; Brown & Rudnick 2011). Due to its large angular size, the analysis of LOFAR data for this cluster requires a special treatment. For this reason, this cluster was excluded from our analysis. Results on this system using LOFAR observations have been recently reported by Bonafede et al. (2021, 2022).
PSZ2 G058.29+18.55 It is a double galaxy cluster, comprising RXC J1825.3+3026 (dubbed subcluster E in Table A.2) and CIZA J1824.1+3029 (dubbed subcluster W in Table A.2). It also known as Lyra complex, in which the presence of a radio halo in RXC J1825.3+3026 was claimed with a targeted LOFAR observation by Botteon et al. (2019a). Due to the bad quality of LoTSS data, we could not perform the analysis of this system.
PSZ2 G063.38+53.44 The candidate radio halo is claimed here for the first time.
PSZ2 G066.34+26.14 The candidate radio halo is claimed here for the first time.
PSZ2 G066.41+27.03 The radio halo is claimed here for the first time.
PSZ2 G067.17+67.46 (A1914) The intricate emission from this system was noticed by early observations with the VLA by Bacchi et al. (2003). Past LOFAR observations found the presence of very steep spectrum emission in the cluster due to very bright revived fossil plasma, a radio halo, and a head-tail radio galaxy (Mandal et al. 2019). Due to the complex emission of this system, we classified it as uncertain. We refer the reader to Mandal et al. (2019) for more details.
PSZ2 G068.36+81.81 We interpreted the diffuse emission from this cluster as the radio lobes of the central AGN. Therefore, we classified this cluster as NDE.
PSZ2 G069.39+68.05 A faint candidate radio halo is noticed in the lowest-resolution images. An elongated region of emission in the S is classified as candidate radio relic. These two sources are claimed here for the first time.
PSZ2 G071.21+28.86 (MACS J1752.0+4440) It is a textbook example of a double radio relic that has been reported in the past using GMRT and WSRT observations (van Weeren et al. 2012; Bonafede et al. 2012). A faint radio halo was claimed by van Weeren et al. (2012) and it is confirmed by our images. This cluster is displayed in Fig. 5.
PSZ2 G071.39+59.54 Central diffuse emission in form of a radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as RH*.
PSZ2 G072.62+41.46 (A2219) It hosts a radio halo detected in the NVSS (Giovannini et al. 1999) and confirmed in deeper VLA observation at 1.4 GHz (Bacchi et al. 2003) and 325 MHz (Orrù et al. 2007). Our LOFAR images of this cluster are affected by bad quality due to the presence of the strong AGN located in the cluster center. Therefore, we could not perform the analysis of this system.
PSZ2 G073.97-27.82 (A2390) Originally, this cluster was classified as a mini-halo by Bacchi et al. (2003) using VLA data at 1.4 GHz. Later, Sommer et al. (2017) analyzed deeper and wide-band VLA data at 1−2 GHz and reclassified the emission as radio halo. Past LOFAR observations revealed that the central radio galaxy may account for most or even all the radio flux that was attributed to the radio halo (Savini et al. 2019). In this paper, we interpret all the radio emission observed in this system as lobes of the central AGN. Another LOFAR work highlighted how these radio lobes are filling the X-ray cavities in the ICM (Bîrzan et al. 2020).
PSZ2 G076.55+60.29 The X-ray emission of the cluster is offset with respect to the position provided by the PSZ2 catalog. Faint emission in the form of a radio halo is noticed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as RH*.
PSZ2 G077.90-26.63 (A2409) Central diffuse emission in form of a radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as RH*.
PSZ2 G080.16+57.65 (A2018) A radio relic in the E direction was discovered together with a candidate radio halo by van Weeren et al. (2021) in the analysis of the LoTSS-DR1 cluster sample. More recently, Paul et al. (2021) confirmed the presence of the two sources using uGMRT data at 300−500 MHz. In this paper, we confirm only the presence of the radio relic.
PSZ2 G080.41-33.24 (A2443) This cluster is known to host complex radio emission with an ultra-steep spectrum in the form of revived fossil plasma (Cohen & Clarke 2011; Clarke et al. 2013), which is well recovered in our LOFAR images.
PSZ2 G080.70+48.31 (A2136) The candidate radio halo is claimed here for the first time.
PSZ2 G081.02+50.57 The radio halo is claimed here for the first time.
PSZ2 G081.72+70.15 (A1838) This system hosts two elongated sources with uncertain origin that are located outside the bulk of the X-ray emission. We did not classified these sources as radio relics as they do not show a sharp surface brightness edge in the LOFAR images.
PSZ2 G083.29-31.03 The radio halo is claimed here for the first time.
PSZ2 G084.10+58.72 The presence of diffuse emission at the center of this cluster was originally noted by previous LOFAR studies (Di Gennaro et al. 2021b,a; van Weeren et al. 2021), belonging to the LoTSS-DR1 area. Here, we confirm the emission and claim the presence of a radio halo.
PSZ2 G084.13-35.41 (A2472) The radio halo is claimed here for the first time. An uncertain source with a wedge-shape is also observed to the S of the radio halo.
PSZ2 G085.23+39.42 The candidate radio halo is claimed here for the first time. The diffuse emission is embedded into two bright structures likely associated with cluster radio galaxies.
PSZ2 G086.58+73.11 The galaxy overdensity seems located toward the SE with respect to the cluster coordinates reported in the PSZ2 catalog. To the S of the optical overdensity, we observed an elongated emission that we classify as candidate radio relic. This relic is claimed here for the first time. The presence of diffuse emission in the central region of the cluster is uncertain due to the artifacts introduced by the source subtraction.
PSZ2 G086.93+53.18 The presence of a radio halo in this cluster has been recently studied with LOFAR by Di Gennaro et al. (2021b,a). This systems belongs to the LoTSS-DR1 area and was also reported by van Weeren et al. (2021). Here, we confirm the emission.
PSZ2 G088.53+41.18 (A2208) This cluster hosts a peculiar emission that is likely associated with radio galaxies. This cluster is displayed in Fig. 5.
PSZ2 G089.39+69.36 The presence of a radio halo in this cluster was recently claimed and studied with LOFAR and uGMRT data (Di Gennaro et al. 2021b,a). Here, we confirm the emission.
PSZ2 G089.52+62.34 (A1904) Two radio relics located on the same side of the cluster (NE) have been firstly claimed by van Weeren et al. (2021) during the analysis of the LoTSS-DR1 cluster sample. We confirm the two detections. Recently, the innermost relic was detected also with the uGMRT at 300−500 MHz (Paul et al. 2021).
PSZ2 G091.79-27.00 The radio relic is claimed here for the first time.
PSZ2 G093.94-38.82 (A2572) This is a system with multiple cluster components in the X-rays, with the main component that we dubbed PSZ2 G093.94-38.82 W and two smaller sub-clumps that we dubbed PSZ2 G093.94-38.82 EN and PSZ2 G093.94-38.82 ES in Table A.2. We discover an extended diffuse source in the S of PSZ2 G093.94-38.82 W that has an uncertain origin.
PSZ2 G094.61-41.24 (A2589) Diffuse emission at the center of this cluster is noticed only in the image tapered at 100 kpc resolution. Due to the possible blend of partially subtracted sources at low resolution, we classified the emission as uncertain.
PSZ2 G095.22+67.41 Patches of emission with low significance and uncertain origin are observed in the cluster center and cluster outskirts. The peripheral emission on the E was classified as a candidate radio relic during the analysis of the LoTSS-DR1 sample (van Weeren et al. 2021).
PSZ2 G095.29+44.13 Central diffuse emission in form of a candidate radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as cRH*.
PSZ2 G096.43-20.89 The candidate radio relic is claimed here for the first time.
PSZ2 G096.83+52.49 (A1995) Extended emission at the center of this cluster was firstly reported by Giovannini et al. (2009) with VLA 1.4 GHz observations. Our images confirm the presence of a radio halo in this system.
PSZ2 G097.72+38.12 (A2218) Diffuse radio emission in the form of a radio halo in this system was originally reported by Moffet & Birkinshaw (1989) and later confirmed by other authors (Giovannini & Feretti 2000; Kempner & Sarazin 2001). In our image, the halo is slightly elongated in the E-W direction and it has an extension toward N. This cluster is displayed in Fig. 5.
PSZ2 G098.30-41.15 The candidate radio halo is claimed here for the first time.
PSZ2 G099.24+42.54 Radio emission with uncertain origin (due to the possible contribution of partially subtracted discrete sources) is noticed in the cluster center. Another emission in the SW with unclear origin is observed.
PSZ2 G099.48+37.72 (A2216) Extended emission is noticed only in the lowest-resolution radio emission. Due to the possible contribution of partially subtracted discrete sources, we classified it as uncertain.
PSZ2 G099.48+55.60 (A1925) This is a double radio relic system and it is claimed here for the first time.
PSZ2 G099.55+34.23 Central diffuse emission in form of a candidate radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as cRH*. Another source elongated in the N-S direction with uncertain origin is located W to the cluster center.
PSZ2 G099.86+58.45 The radio halo in this cluster was claimed by previous LOFAR observations (Cassano et al. 2019; Di Gennaro et al. 2021b,a; van Weeren et al. 2021) and is confirmed in our images.
PSZ2 G100.14+41.67 (A2146) This cluster hosts a double radio relic system that has been observed with the VLA and LOFAR (Hlavacek-Larrondo et al. 2018; Hoang et al. 2019). Due to the bad data quality, we could not perform the analysis of this system.
PSZ2 G100.45-38.42 (A2626) Peculiar emission in form of four arcs was observed in this system (Gitti et al. 2004; Gitti 2013; Kale & Gitti 2017; Ignesti et al. 2017), which led to the nickname of Kite cluster. Recent LOFAR observations suggest that the emission originate from the central radio galaxy whose fossil radio plasma has been compressed and revived as a consequence of motions of the ICM (Ignesti et al. 2020).
PSZ2 G100.96-24.07 The candidate radio halo is claimed here for the first time.
PSZ2 G105.55+77.21 (A1691) Extended emission is noticed only in the lowest-resolution radio image. Due to the possible contribution of partially subtracted discrete sources, we classified it as uncertain.
PSZ2 G106.41+50.82 (A1918) Possible diffuse radio emission is observed at the center of the cluster. However, due to an artifact introduced by the discrete source subtraction, we classified the diffuse emission as uncertain.
PSZ2 G106.61+66.71 Central diffuse emission in form of a radio halo was discovered by van Weeren et al. (2021) in the analysis of the LoTSS-DR1 cluster sample and was later confirmed by Paul et al. (2021). Our images agree with previous results.
PSZ2 G107.10+65.32 (A1758) This is a double galaxy cluster (dubbed PSZ2 G107.10+65.32 N and PSZ2 G107.10+65.32 S in Table A.2) that has been intensively studied in the past with different instruments (Kempner & Sarazin 2001; Giovannini et al. 2009; Venturi et al. 2013; Botteon et al. 2018, 2020a; Schellenberger et al. 2019; van Weeren et al. 2021). It hosts a double radio halo system (Botteon et al. 2018), and the two cluster components are connected with a bridge of radio emission (Botteon et al. 2020a). The northern cluster hosts revived fossil plasma emission and the southern system hosts a radio relic. Our images agree with previous results.
PSZ2 G108.27+48.66 The candidate radio halo is claimed here for the first time.
PSZ2 G109.22-44.01 The candidate radio relic is claimed here for the first time.
PSZ2 G109.97+52.84 The radio halo is claimed here for the first time.
PSZ2 G111.75+70.37 (A1697) The relic and halo in this system were firstly reported by Paul et al. (2020). This cluster belongs to the LoTSS-DR1 sample presented by van Weeren et al. (2021). Our images confirm the presence of an elongated radio halo and a relic in the NE with a long trail of emission.
PSZ2 G112.48+56.99 (A1767) The radio halo is claimed here for the first time.
PSZ2 G113.29-29.69 Diffuse emission is observed at the center of the cluster. Due to the artifacts introduced by the subtraction of discrete sources, we classified it as uncertain.
PSZ2 G113.91-37.01 The radio halo and relics are claimed here for the first time. This is a double radio relic system. This cluster is displayed in Fig. 5.
PSZ2 G114.31+64.89 (A1703) Hints of diffuse radio emission in this cluster were noticed with the VLA by Owen et al. (1999) and with LOFAR by Savini et al. (2018b). Recently, van Weeren et al. (2021) confirmed the presence of a radio halo during the analysis of the LoTSS-DR1 sample. Our images agree with previous findings.
PSZ2 G114.99+70.36 (A1682) This systems hosts complex diffuse emission associated with the ICM and radio galaxies that has been studied in the past with GMRT and LOFAR (Venturi et al. 2008, 2011, 2013; Macario et al. 2013; Clarke et al. 2019; van Weeren et al. 2021). Due to the difficulty of disentangling the diffuse emission, we classified it as uncertain. This cluster is displayed in Fig. 5.
PSZ2 G116.32-36.33 There are two separate X-ray components for this cluster, RX J0027.8+2616 and WHL J002727.0+260707. They are dubbed PSZ2 G116.32-36.33 N and PSZ2 G116.32-36.33 S in Table A.2 respectively. RX J0027.8+2616 is the target of the Chandra observation and it is within the FoV of the XMM-Newton observation targeted at WHL J002727.0+260707.
PSZ2 G116.50-44.47 Central diffuse emission in form of a radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as RH*. We also report for the first time a radio relic located in the N.
PSZ2 G118.34+68.79 The presence of a candidate radio halo was claimed during the analysis of the LoTSS-DR1 sample (van Weeren et al. 2021). Our images confirm the presence of central diffuse emission in form of a candidate radio halo. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as cRH*.
PSZ2 G121.03+57.02 The radio relic in the S of the cluster is located where there is an apparent jump of the X-ray surface brightness and it is claimed here for the first time. The emission in the center has uncertain origin due to the possible residual of discrete source subtraction.
PSZ2 G123.00-35.52 Central diffuse emission in form of a radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as RH*.
PSZ2 G124.20-36.48 (A115) A prominent and bright radio relic in the N periphery of the system was originally claimed with VLA 1.4 GHz observations by Govoni et al. (2001). Follow-up studies focused on the connection between relic and shock (Botteon et al. 2016), also exploiting GMRT observations at 610 MHz (Hallman et al. 2018). Due to the bad quality of the LoTSS data, we could not study the radio properties of this system. The PSZ2 source comprises two subclusters, which are dubbed PSZ2 G124.20-36.48 N and PSZ2 G124.20-36.48 S in Table A.2.
PSZ2G126.28+65.62 A radio halo in this system was claimed by Di Gennaro et al. (2021a) with LOFAR data at 144 MHz. Follow-up uGMRT observations at 550−900 MHz did not show the presence of diffuse radio emission in the cluster possibly due to its steep spectrum (Di Gennaro et al. 2021b). Owing to the difficulty of separating the contribution of discrete sources from the diffuse emission, we classified this cluster as uncertain.
PSZ2 G133.59+50.68 Central diffuse emission in form of a radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as RH*.
PSZ2 G133.60+69.04 (A1550) The radio halo in this cluster was originally claimed by Govoni et al. (2012) and recently confirmed during the analysis of the LoTSS-DR1 sample (van Weeren et al. 2021). We confirm the presence of a radio halo elongated in the E-W direction and the presence of possible revived fossil plasma. This system will be analyzed in detail in a forthcoming publication (Pasini et al., in preparation).
PSZ2 G135.17+65.43 This cluster was recently studied by van Weeren et al. (2021) as part of the LoTSS-DR1 cluster sample. They found the presence of a halo and relic (source G in van Weeren et al. 2021). We confirm the previous classification.
PSZ2 G138.32-39.82 The radio halo in this cluster was originally claimed by Kale et al. (2013) in the context of the Extended GMRT Radio Halo Survey. Our images confirm the detection.
PSZ2 G139.00+50.92 (A1351) The existence of a halo source in this cluster was firstly reported by Owen et al. (1999) and later studied by different authors (Giacintucci et al. 2009; Giovannini et al. 2009). Our images confirm the presence of radio halo emission that is contaminated by two bright cluster AGN.
PSZ2 G141.05-32.61 The candidate radio halo is claimed here for the first time.
PSZ2 G143.26+65.24 (A1430) The presence of a radio halo in this cluster was claimed by van Weeren et al. (2021) during the analysis of the LoTSS-DR1 sample, and later confirmed by Paul et al. (2021). A dedicated work on this system, composed of two cluster components, using LOFAR data has been recently published by Hoeft et al. (2021).
PSZ2 G143.44+53.66 A number of diffuse radio sources are noted in the cluster center. These sources are claimed here for the first time, and their nature is still uncertain. This cluster is displayed in Fig. 5.
PSZ2 G144.23-18.19 The candidate radio halo is claimed here for the first time. This is the only target in our sample in which the indication of diffuse emission in the cluster center is noticed in a cluster without redshift/mass.
PSZ2 G144.99-24.64 The candidate radio relic is claimed here for the first time. Another diffuse emission with uncertain origin is observed at the center of our images.
PSZ2 G145.65+59.30 A candidate radio halo was claimed in this system during the analysis of the LoTSS-DR1 sample (van Weeren et al. 2021). We classified this emission as uncertain as its location is offset with respect to the bulk of the X-ray emission.
PSZ2 G147.88+53.24 The presence of central diffuse emission in the form of a radio halo was recently claimed by Di Gennaro et al. (2021a) and studied with deeper observations in the context of the LoTSS Deep Fields (Osinga et al. 2021) and with uGMRT data (Di Gennaro et al. 2021b). Our images confirm the emission (which is classified as candidate radio halo due to the lack of X-ray data).
PSZ2 G148.36+75.23 The radio halo is claimed here for the first time. The diffuse emission surrounds the central cluster AGN. This cluster is displayed in Fig. 5.
PSZ2 G149.22+54.18 (A1132) The radio halo in this system was claimed with LOFAR observations by Wilber et al. (2018) and later studied with deeper observations in the context of the LoTSS Deep Fields (Osinga et al. 2021). Our images agree with previous results.
PSZ2 G149.75+34.68 (A665) The presence of diffuse emission coincident with the cluster center was detected in early observations (Moffet & Birkinshaw 1989; Giovannini & Feretti 2000; Kempner & Sarazin 2001), and it was further studied by other authors (Feretti et al. 2004; Vacca et al. 2010; George et al. 2021). Our images show that the emission entirely fills the central volume of the cluster and has a sharp surface brightness edges toward N. This cluster is displayed in Fig. 3.
PSZ2 G150.56+46.67 The candidate radio halo is claimed here for the first time.
PSZ2 G150.56+58.32 The presence of a radio halo elongated in the NW-SE direction was recently claimed by van Weeren et al. (2021) during the analysis of the LoTSS-DR1 sample. Our images agree with previous results.
PSZ2 G151.19+48.27 (A959) Diffuse radio emission in this system was reported in Owen et al. (1999). More recently, this system was studied in detail with LOFAR and GMRT observations by Bîrzan et al. (2019) who confirmed the radio halo and identified a radio relic in the SE. Our images agree with these results.
PSZ2 G154.13+40.19 The candidate radio halo is claimed here for the first time. The diffuse emission surrounds the central cluster AGN.
PSZ2 G156.26+59.64 The presence of a candidate radio halo in this cluster was reported by van Weeren et al. (2021) during the analysis of the LoTSS-DR1 sample. Here we confirm the previous classification.
PSZ2 G160.83+81.66 The presence of a radio halo in this system was recently claimed by Di Gennaro et al. (2021a) and followed-up with the uGMRT (Di Gennaro et al. 2021b). Our images confirm the emission. This is the highest-z radio halo cluster observed to date.
PSZ2 G163.61+34.30 Central diffuse emission in form of a candidate radio halo is observed here for the first time. Since the fitting of the radio halo surface brightness profile does not provide reliable results, we classified this emission as cRH*.
PSZ2 G164.65+46.37 The radio halo is claimed here for the first time. The emission is elongated along the NW-SE direction.
PSZ2 G165.06+54.13 (A990) The radio halo in this cluster was recently claimed by Hoang et al. (2021b) using LOFAR observations. Our images agree with previous results.
PSZ2 G165.46+66.15 (A1240) This cluster hosts a double radio relic system. The two relics are located to the N and S, and were firstly hinted by Kempner & Sarazin (2001), and later firmly claimed by Bonafede et al. (2009). This cluster has been the subject of a dedicated study by Hoang et al. (2018), who employed LOFAR, GMRT, and VLA observations. Our images agree with previous findings.
PSZ2 G166.09+43.38 (A773) The existence of diffuse emission in the central region of the cluster was suggested by NVSS and WEsterbork Northern Sky Survey (WENSS) observations (Giovannini et al. 1999; Kempner & Sarazin 2001). The cluster was then targeted with the VLA, and the presence of a radio halo was firmly claimed (Govoni et al. 2001; Cuciti et al. 2021). In our images, the halo appears more extended than previously observed.
PSZ2 G166.62+42.13 (A746) The prominent radio relic in the W of this system was firstly claimed by van Weeren et al. (2011). In addition, our images highlight the presence of other two fainter and smaller relics in the E and NE directions. An elongated structure with high surface brightness surrounded by extended diffuse emission is observed toward the cluster center. We classify the extended emission a radio halo, but did not provide a flux density due to the difficulty of disentangling the emission from the bright elongated structure, which possibly traces a radio relic projected on the center of the cluster. A detailed study of this system will be presented in a forthcoming paper. This cluster is displayed in Fig. 5.
PSZ2 G172.63+35.15 (A655) The radio halo is claimed here for the first time and will be the focus of a forthcoming work (Groeneveld et al., in preparation).
PSZ2 G176.27+37.54 The candidate radio halo is claimed here for the first time. The diffuse emission surrounds the central cluster AGN.
PSZ2 G179.09+60.12 (A1068) The radio halo is claimed here for the first time, although hints of diffuse emission in the cluster center were noticed in the past by Govoni et al. (2009) using VLA observations at 1.4 GHz. The detailed analysis of this system will be presented in a forthcoming paper (Biava et al., in preparation).
PSZ2 G181.06+48.47 This is a double radio relic system and it is claimed here for the first time.
PSZ2 G183.30+34.98 The candidate radio halo is claimed here for the first time. A diffuse emission elongated in the N-S direction with uncertain origin is also observed to the E of the cluster center.
PSZ2 G184.68+28.91 (A611) The candidate radio halo is claimed here for the first time. The diffuse emission surrounds the central cluster AGN.
PSZ2 G186.37+37.26 (A697) Hints of diffuse emission located in the cluster center and were firstly claimed by Kempner & Sarazin (2001) using WENSS observations. The radio halo was later confirmed and studied in detail with the VLA, GMRT and WSRT (Venturi et al. 2008; van Weeren et al. 2011; Macario et al. 2010, 2013). Our images agree with previous results.
PSZ2 G186.61+62.94 The candidate radio halo is claimed here for the first time.
PSZ2 G186.99+38.65 The radio halo and radio relic are claimed here for the first time. The radio halo is slightly elongated in the N-S direction while the relic lays in the NE outskirts of the cluster. Curiously, this relic has a mildly convex morphology; to the best of our knowledge, such an unusual curvature has been reported only for the relics in SPT-CL J2023-5535 (HyeongHan et al. 2020), the Ant Cluster (Botteon et al. 2021a), and Abell 3266 (Riseley et al., in preparation).
PSZ2 G189.31+59.24 (A1033) This cluster hosts revived fossil plasma (de Gasperin et al. 2015) and the prototype of Gently Re-Energized Tail (GReET; de Gasperin et al. 2017). In our images, diffuse radio emission in form of a radio halo is observed to the N of the GReET. This radio halo is claimed here for the first time. A dedicated paper on this target is currently in preparation (Edler et al., in preparation).
PSZ2 G190.61+66.46 The radio halo and radio relic are claimed here for the first time. The radio halo is elongated in the E-W direction while the relic is located to the W of the halo emission.
PSZ2 G192.18+56.12 (A961) The radio halo is claimed here for the first time.
PSZ2 G192.77+33.14 (A671) The candidate radio halo is claimed here for the first time. If confirmed, this would be the radio halo in the lowest-mass cluster observed to date.
PSZ2 G192.90+29.63 The candidate radio halo is claimed here for the first time.
PSZ2 G195.60+44.06 (A781) This is a complex system with multiple galaxy cluster components (e.g., Sehgal et al. 2008) with two components at redshift z ∼ 0.3 (dubbed PSZ2 G195.60+44.06 E2 and PSZ2 G195.60+44.06 W1 in Table A.2) and the other two (dubbed PSZ2 G195.60+44.06 E1 and PSZ2 G195.60+44.06 W1 in Table A.2) at redshift z ∼ 0.43. The presence of a radio halo in the main component of A781 (PSZ2 G195.60+44.06 E2) was disputed in the literature, being observed with the VLA at 1.4 GHz (Govoni et al. 2011) but not with the GMRT at lower frequency (Venturi et al. 2008, 2011, 2013). Deeper observations performed with LOFAR did not confirm the presence of a radio halo in this system (Botteon et al. 2019b). Our images agree with the latter results. The peripheral diffuse emission in the SE was originally believed to trace a radio relic (Venturi et al. 2008; Govoni et al. 2011), while recently it has been proposed that it results from the interaction between a weak shock and a radio galaxy (Botteon et al. 2019b).
PSZ2 G198.46+46.01 The candidate radio relic is claimed here for the first time. It is located in the S and it is elongated in the NE-SW direction. Another diffuse, elongated, emission with uncertain origin is observed toward the cluster center.
PSZ2 G205.90+73.76 The radio halo and double radio relic system are claimed here for the first time. The radio halo is elongated in the N-S direction. The two relics are located in symmetric directions (N-S) with respect to the radio halo emission.
All Tables
Subsample of galaxy clusters with Chandra/XMM-Newton data available for which we computed the X-ray morphological parameters.
Sample of radio halos and candidate radio halos that were classified with an asterisk (the RH* and cRH* in Table A.1) because of their low significance. The integrated flux density of these sources cannot be determined accurately with current data.
All Figures
![]() |
Fig. 1. Redshift-mass distribution of PSZ2 sources. Clusters that are located in the LoTSS-DR2 area are indicated in blue. The histograms show the number of clusters at various redshifts and masses; the dashed and dotted lines mark the median values of the LoTSS-DR2 sample and the full PSZ2 sample, respectively. Similarly to the full PSZ2 sample, our sample spans a wide range of redshifts and masses. |
In the text |
![]() |
Fig. 2. Position of the PSZ2 clusters in the RA-13 (top) and RA-1 (bottom) regions covered by LoTSS-DR2. The color code indicates the redshift of the cluster. The radius of the circle is proportional to M500. Clusters without redshift and mass are reported as black crosses. The background image represents the noise variations in LoTSS-DR2 (darker colors denote higher noise values) and is reproduced from Shimwell et al. (2022). |
In the text |
![]() |
Fig. 3. Example of the set of images that we produced for each object in our sample. The displayed cluster is PSZ2 G149.75+34.68, and the reported panels show (from left to right and from top to bottom) the reference image (robust = −0.5), the high-resolution image (robust = −1.25), taper 25, 50, and 100 kpc images with and without discrete sources, the clean model components used for the source subtraction, contours from the reference radio image starting from 3σ overlaid on an optical image (Pan-STARRS g, r, i), and taper 50 and 100 kpc discrete-source-subtracted contours starting from 2σ overlaid on an X-ray image (Chandra or XMM-Newton). Contours are always spaced by a factor of 2. The beam is shown in the bottom-left corner, and the mass (M500, in ×1014 M⊙ units), redshift (z), and image noise (rms, in mJy beam−1 units) are reported in the top-left corner. The circle denotes r500 and is centered at the coordinates reported in the PSZ2 catalog. The images are available at full resolution as well as in FITS format on the project website (https://lofar-surveys.org/planck_dr2.html). |
In the text |
![]() |
Fig. 4. Decision tree used to classify the diffuse radio sources in the ICM. The classes of objects that form the end points of the decision tree are described in Sect. 4. |
In the text |
![]() |
Fig. 5. Collection of clusters that show several types of radio emission. PSZ2 G053.53+59.52 has a central radio halo and a number of sources of uncertain origin. PSZ2 G071.21+28.86 has a double radio relic. PSZ2 G088.53+41.18 is a system without diffuse cluster emission as the only emission detected is associated with an optical galaxy. PSZ2 G097.72+38.12 has a radio halo. PSZ2 G113.91–37.01 has a radio halo and two relics. PSZ2 G114.99+70.36 has emission of uncertain origin. PSZ2 G143.44+53.66 has emission of uncertain origin. PSZ2 G148.36+75.23 has a radio halo. PSZ2 G166.62+42.13 has a radio halo and multiple relics. We remark that the classification was done by inspecting all the sets of images available, while the images displayed in this gallery are only the reference ones (i.e., those with robust = −0.5 and no uv taper). For a more complete picture of the large variety of radio structures observed in our sample, we refer the reader to Fig. B.1 and to the project website. |
In the text |
![]() |
Fig. 6. Results obtained by fitting the radio halo shown in Fig. 3 with HALO-FDCA (Boxelaar et al. 2021). Top figures: image used for the fit with overlaid: the contours (white circles) of the best-fit circular model drawn at [1, 2, 4, 8, …] × σ (left panel), the image of the best-fit model (central panel), and the residual image of the fit, with the 2σ model contour denoted by the black circle (right panel). Contaminating sources are masked out and are highlighted by the green and red regions (left and right panels, respectively). Bottom figure: MCMC corner plot presenting the distribution of the posteriors of each fitted parameter (Foreman-Mackey et al. 2013). |
In the text |
![]() |
Fig. 7.
|
In the text |
![]() |
Fig. 8. Summary of the number of clusters classified in our sample divided per category (NDE = no diffuse emission; U = uncertain; RH = radio halo; RR = radio relic; c = candidate; * = flux density measurement not reliable; N/A = not applicable). We note that a cluster can be classified under multiple categories, such as both RH and RR. The percentage on the y axis is with respect to the total number of PSZ2 in LoTSS-DR2 (309 targets). |
In the text |
![]() |
Fig. 9. Histograms of the mass (top panels) and redshift (bottom panels) distribution for the main categories of objects classified in this work. |
In the text |
![]() |
Fig. B.1. LOFAR image gallery of the targets in our sample, excluding the five objects for which it was not possible to apply the extraction + recalibration scheme (see end of Sect. 3.1). The collage is available at full resolution on the project website, https://lofar-surveys.org/planck_dr2.html. |
In the text |
![]() |
Fig. B.1. continued. |
In the text |
![]() |
Fig. B.1. continued. |
In the text |
![]() |
Fig. B.1. continued. |
In the text |
![]() |
Fig. C.1. Flux density measurements for radio halos. Left: Comparison between the flux density derived with HALO-FDCA (Sfit) and that obtained manually within the 2σ contour (S2σ). The dashed line indicates the linear correlation as a reference. Right: Diagnostic plot used to determine when HALO-FDCA does not provide a reliable flux density measurement due to the low significance of the radio halo. The dashed lines indicate the thresholds of Sfit/S2σ > 1.5 and I0/3σ < 2 that were used to identify the sources that were classified as radio halos or candidate radio halos with the decision tree of Fig. 4, but whose flux density was not considered reliable. These sources are reported in red and in the manuscript are referred to as RH* and cRH*. |
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.