A GMRT Narrowband vs. Wideband Analysis of the ACT-CLJ0034.4+0225 Field Selected from the ACTPol Cluster Sample

Low frequency radio observations of galaxy clusters are a useful probe of the non-thermal intracluster medium (ICM), through observations of diffuse radio emission such as radio halos and relics. Current formation theories cannot fully account for some of the observed properties of this emission. In this study, we focus on the development of interferometric techniques for extracting extended, faint diffuse emissions in the presence of bright, compact sources in wide-field and broadband continuum imaging data. We aim to apply these techniques to the study of radio halos, relics and radio mini-halos using a uniformly selected and complete sample of galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect by the Atacama Cosmology Telescope (ACT) project, and its polarimetric extension (ACTPol). We use the upgraded Giant Metrewave Radio Telescope (uGMRT) for targeted radio observations of a sample of 40 clusters. We present an overview of our sample, confirm the detection of a radio halo in ACT-CL J0034.4+0225, and compare the narrowband and wideband analysis results for this cluster. Due to the complexity of the ACT-CL J0034.4+0225 field, we use three pipelines to process the wideband data. We conclude that the experimental SPAM wideband pipeline produces the best results for this particular field. However, due to the severe artefacts in the field, further analysis is required to improve the image quality.


Introduction
Non-thermal diffuse emission in galaxy clusters was first discovered in the Coma cluster [1][2][3]. This discovery was a confirmation of the existence of relativistic electrons and magnetic fields in the intracluster medium (ICM) [4,5]. Traditionally, diffuse radio emissions have been categorized into three groups based on morphology, size, and cluster dynamics; giant radio halos (GRHs), radio mini-halos (RMHs), and radio relics (RRs). The size of these non-thermal diffuse structures, which extend over kpc to Mpc scales, raised questions on their formation mechanisms. One prominent question is how the cosmic ray electrons (CRes) are (re)accelerated and transported given their diffusion timescale limitations [6]. Many studies allow electrons to reach relativistic speeds via the shock drift acceleration mechanism (SDA; [47,48]). However, the role of magnetic fields and its amplification by low Mach number shocks is still poorly constrained [35].
A new generation of telescopes such as the LOw-Frequency ARray (LOFAR; [49]), MeerKAT [50], and the uGMRT [51] have opened a window into a new group of ultra-steep radio sources which were previously undetectable due to frequency and sensitivity constraints. The new studies have resulted in a different class of diffuse emission, such as; radio phoenices and gently re-energized tails (GReETs; [26,52]). The sensitivity of these telescopes have opened up a new observational window which probes lower mass (<4 × 10 14 M ) and higher redshift (z > 0.3) clusters [53][54][55]. Large cluster samples of this nature will help in the refinement of currently existing formation theories. However, the broad bandwidths and large fields of view of these telescopes also present multiple data reduction challenges.
In this paper, we introduce an SZ-selected sample from the ACTPol observations. We will also discuss the challenges we faced with the uGMRT data reduction and present a case study of ACT−CL J0034.4+0225 (hereafter J0034). In Section 2, we introduce the cluster sample. In Section 3, we discuss the data reduction challenges and the pipelines we explored. In Section 4, we study J0034 and compare the narrowband and wideband results. Finally, we summarise our findings and future outlook in Section 5. We adopt a ΛCDM flat cosmology with H 0 = 70 km s −1 Mpc −1 , Ω m = 0.3, and Ω Λ = 0.7. For our radio spectral index calculations, we assume S ν ∝ ν −α , where S ν is the flux density at frequency ν and α is the spectral index.

The ACTPol Sample
Until recently, statistical studies of diffuse radio emission in galaxy clusters have been constrained to observations of high mass and low redshift systems selected using X-ray telescopes [20,56]. The X-ray selected samples had relatively low diffuse emission detection rates (∼20%). Andrade-Santos et al. [57] studied the fraction of cool-core clusters in an X-ray selected sample from Chandra versus an SZ selected sample from Planck. Their study revealed that the X-ray sample had a higher fraction (∼44%) of cool-core clusters in comparison to the SZ sample (∼28%). The majority of cool-core clusters do not host GRHs and RRs; and hence would attribute to the low detection rates. These findings in conjunction with the poor resolution (∼1.8 ) of ROSAT, which was mainly used for the X-ray samples, may account for the low detection rates of diffuse emission in X-ray selected samples. Such constraints have since led to preferentially using SZ-selected samples for scaling relation studies.
The thermal SZ effect is a powerful probe for high-redshift clusters due to the fact that it is not affected by dimming. Hence, the SZ survey catalogues offer an almost redshift independent selection function. Essentially, cluster samples from SZ surveys are restricted by cluster masses based on the sensitivity of the observing instrument. The three main microwave telescopes that have been used to detect galaxy clusters are the ACT, the Planck satellite [58], and the South Pole Telescope (SPT; [59]).
For our project, we use the galaxy cluster catalogue from the ACT's Polarimetric extension (ACTPol; [60]) to select our sample. The ACTPol catalogue was constructed using the E−D56 region, shown in Figure 1, which covers an area of 987.5 deg 2 [61]. The final catalogue consists of 182 optically confirmed clusters. For a signal-to-noise (SNR) > 4 the sample spans a mass and redshift range of 1.6 < 10 14 M < 9.1, and 0.1 < z < 1.4, respectively. The 90% sample completeness cut-off is M 500c > 4.5 × 10 14 M for SNR > 5 and a redshift range of 0.2 < z < 1.0. M 500c is the mass measured within a radius that encloses a region with an average density that is 500 times the critical density at the cluster redshift, and assuming the SZ-signal scales with mass as described in Arnaud et al. [62].
For our sample, we considered all the clusters with an SNR > 5. We then applied a mass and redshift cut of M 500,SZ > 4 × 10 14 M and 0.1 < z < 0.8. This results in a sample of 40 clusters, which formed the basis of our sample. Given that the ACTPol catalogue is 90% complete for clusters within 0.2 < z < 1.0 and M 500c > 4.5 × 10 14 M , we derive our sample completeness to be 68%. This completeness is derived from the mass completeness and the available radio information (30/40). In Figure 2, we compare our sample of 40 clusters to Planck and SPT samples, which are the most recent statistical studies using SZ-selected cluster samples. Our sample covers a wider redshift range compared to both samples. We also cover lower mass clusters compared to the Planck sample. The SPT sample overlaps with ours at higher redshifts; however, it does not cover lower redshifts (z < 0.33). This study is a precursor to studies of larger samples exploring lower mass and higher redshift cluster sample studies such as the MeerKAT Extended Relics, Giant Halos, and Extragalactic Radio Sources (MERGHERS) survey [63] and the MeerKAT Absorption Line Survey (MALS; [64]).
We chose to observe these clusters at low-frequency (250-500 MHz) using the wideband uGMRT in addition to using existing archival narrowband GMRT data. As of November 2020, there were 17 clusters with pre-existing GMRT legacy software backend (GSB) observations. The GSB has a bandwidth of 32 MHz. We obtained uGMRT wideband (GWB) observations for 13 clusters, taken over three observing semesters. These clusters were observed in GMRT's observing cycles 32, 33, and 36 (proposal IDs: 32_012, 33_010, 36_050). In total, we have 30 clusters with narrowband and/or wideband observations. For this paper, we focus on the data reduction and analysis of J0034, which has both narrowband and wideband observations. We selected this field because it was most severely affected by RFI and has multiple bright sources near the cluster region. Hence, successfully calibrating and imaging this field would imply that we could apply the techniques on the rest of the sample.  [25] and SPT [55] sample used for statistical studies of diffuse emission.

Data Reduction
For this paper, we focus on the case study of the complex field J0034. We start by reducing the narrowband observations using the source peeling and atmospheric modelling (SPAM; [65,66]). The final narrowband image indicates that this field has bright sources near the cluster's region. Such sources complicate the calibration procedure. Such a field is a good case study for testing calibration and imaging strategies. Hence, when reducing the wideband data, we use three pipelines and compare the resulting images for J0034. In the following sections, we discuss the narrowband and wideband data reduction strategies.

Legacy GMRT GSB
We use SPAM's calibration pipeline version 19.11.07 to reduce the GSB data. SPAM is a fully automated Python pipeline that uses the PARSELTONGUE [67] interface to access and execute the Astronomical Image Processing System (AIPS; [68]) tasks. The pipeline is divided into two steps; the 'pre-calibrate targets' step and the 'process target' step. The data reduction steps are as follows.

1.
The 'pre-calibrate target' step performs the cross-calibration step of the standard calibration procedure. The pipeline uses the primary calibrator(s) to determine the channels affected by RFI, determine the flux scale, and to produce cross calibration tables. The flags and calibration tables are then applied to the target source. Finally, the calibrator and target visibilities are split into separate UVFITS files.

2.
The 'process target' step begins by taking the cross-calibrated (1GC) target data and applying 2GC (or self calibration). For 2GC, the target visibilities are imaged using the facet-based method established by Cornwell and Perley [69]. The point sources covering the primary beam are obtained from a sky model extracted from the VLA low-frequency sky survey (VLSS; [70]) and the NRAO VLA sky survey (NVSS; [71]). The observed field is then faceted based on the sky model. The deconvolution is done using the Cotton-Schwab CLEAN algorithm [71]. The CLEANED visibilities are calibrated using the sky model. Then the calibrated visibilities are imaged to produce a better sky model and improved calibration solutions. This self calibration cycle is applied three times. For the wide-field imaging, the pipeline performs a single-scale CLEAN deconvolution down to 3 times the central background noise (σ), using automated CLEAN boxes placed at positive peaks of at least 5σ.

3.
The second part of the 'process target' step is to correct for DDEs (3GC), which include ionospheric effects [72]. For each 3GC cycle, the brightest source in the field is peeled. The ionospheric effects of the brightest source are modelled and phase corrected solutions for the peeled model are obtained. SPAM uses a phase screen model to extend the ionospheric phase solutions of the single source to the full field. Then the solutions are applied to each facet during imaging. This cycle is repeated six times, and each time the brightest source is selected and corrected for DDEs. This repeats until all the bright sources are corrected for DDEs. Finally, the astrometric corrections are applied to the DDEs corrected image and the image is primary beam corrected thereafter.
We used default data reduction parameters throughout the different stages of the pipeline. The default setting of the number of pixels is ∼3780, with a pixel size of ∼1.9 and the robust parameter is 1. The resulting image, shown in the bottom right panel of Figure 3, has a rms of 0.10 mJy/beam near the cluster centre, and the beam size is 18.1 × 14.2 , p.a. 60.3 • .

uGMRT Data Reduction
Thirteen clusters in our sample, including J0034, have uGMRT GWB band 3 (250-500 MHz) observations. The wide bandwidth of the upgraded GMRT simultaneously provides increased sensitivity and opportunity for in-band spectral index studies. However, the calibration and imaging of GWB data is a complex procedure. The main challenges with reducing the wideband, wide-field data are as follows. The primary beam pattern is dependent on the observing frequency, hence, for wide bandwidth observations, this pattern varies across the band. The flux density of radio sources correlates with the frequency (S ν ∝ ν −α ). As a result of this correlation, the spectral models of the point sources need to be taken into account during imaging. The uGMRT observations were carried out at low-frequencies, which resulted in various directional dependent effects, such as ionospheric effects. The GWB observations of J0034 are particularly affected by radio frequency interference (RFI), which corrupts data if not fully modelled and removed. The level of RFI in our data makes the reduction process onerous and leads to images with significantly compromised quality. The resulting images had rms noise values three to seven times higher than the theoretically predicted noise floor of 10 µJy/beam. Another challenge for our sample was that it is in the equatorial region (−7.2 • < δ < 4 • ), so despite the wide bandwidth and hours of integration time, uv coverage is more sparse because aperture synthesis is not as effective for equatorial sources (see Figure 4). The poor sampling of the visibility space results in north-south artefacts around bright sources, which are a reflection of the poor sampling function. To overcome these challenges, we explored various pipelines that deal with the intensive RFI flagging and calibration direction-dependent effects. Once the data reduction pipelines are refined, we plan to reduce the remaining wideband data to determine if we can detect lower surface brightness diffuse emission in clusters with non-detections in the GSB data, and enhanced features of the existing diffuse emission detections in other clusters. The pipelines we explored are summarised below, including the image comparison of the J0034 results.

CASA Pipeline
The CASA pipeline we adopt is a preliminary version of the CAPTURE pipeline [73]. For this pipeline, flagging, calibration and imaging is done using in-built Common Astronomy Software Applications (CASA; [74]) tasks. Firstly, we flag RFI using the manual and automated flagdata tasks. We apply the auto-flag separately for the different fields. For auto-flagging, we use the tfcrop mode which identifies and removes outliers on the time-frequency plane. We set the flagging threshold parameters higher for the calibrator fields in comparison to the target field because these fields are usually much brighter and are detected at a higher SNR. The timecutoff and freqcutoff deviation parameters for the calibrators were set to 5.0, while for the target they were set to 6.0. These values place constraints on the deviation of data points from the fitted time and frequency polynomial. Data points with higher deviation values are regarded as RFI.
We then apply a cycle of cross-calibration. We begin by setting the flux scale using the standard models from [75]. Thereafter, we produce calibration tables. The delay calibration solution interval is 10 min, and we only have one solution interval for bandpass and gain calibration. We transfer the flux and phase calibration solutions to the target. After the first cycle of cross-calibration, we apply flagging once more. We do this in order to excise low-level RFI only discernible after calibration. This time we use tfcrop and rflag mode. The rflag mode calculates statistics per time chunk and set thresholds which indicate the outliers that need to be flagged. We apply rflag post calibration because it tends to result in higher flag percentages if applied on uncalibrated data. Finally, we produce calibration solutions from the second cycle of cross-calibration and apply these solutions to the target.
We then separate the target visibilities and begin with 2GC. We apply multiple cycles of phase-only calibration on the data, we stop once the quality of the image is no longer improving. For most data sets, the four cycles were sufficient. We then apply multiple cycles of phase and amplitude calibration until the noise quality of the images reaches a plateau, three cycles were sufficient for most datasets. For both phase-only and phase and amplitude calibration, we begin with a solution interval of 16 min and decrease the interval per cycle by dividing the initial solution interval by a factor of two times the number of the cycle. Our imaging during self-calibration is done using tclean's mtmfs deconvolver, with nterms = 2 and robust parameter = 0. This deconvolver accounts for the varying spectral indices of sources across the wide bandwidth. The number of iterations is set to 2500 and increases per cycle by a factor of 2 n , where n is the cycle number. The mask is set to auto-multithreshold. The CLEANing threshold is 0.01 mJy while the sidelobe threshold is set to 3σ.
The resulting image, shown in the top left panel of Figure 3, has a rms of 89.2 µJ/beam and the beam size is 8.0 × 4.9 , p.a. 60.3 • . The GSB image is included in the panels for comparison purposes. After rigorous RFI excision strategies, including AOFLAGGER [76], the north-south artefacts around the bright sources were still visible. This led us to conclude that these artefacts and phase variations in the bright sources are due to DDEs. The results of this pipeline were unsatisfactory for our science goals, since we needed to correctly remove bright point sources in the visibilities. Hence, we explored other pipelines which apply DDEs calibration (3GC).

DDFacet/killMS Pipeline
The KILLMS 2 (the pipeline was accessed in August 2019) and DDFACET [77] based pipeline attempts to correct for DDEs by solving the full Jones matrix. This pipeline starts with a measurement set (MS) that has been calibrated using the CASA pipeline. From our investigation of the CASA pipeline, we found that more cycles of 2GC calibration resulted in images with lower peak fluxes. Hence, for this pipeline, 2GC is only applied once for phase-only calibration. We use the image produced after 2GC to create a mask using a threshold of 10σ to ensure we create a sky model that does not contain residual artefacts. Thereafter, we image the visibilities using DDFACET. The imaging parameters used are as follows. We apply three major cycles and 90,000 minor cycles with a deconvolution peak factor of 0.001. These parameters are fixed even for the post-KILLMS imaging cycle.
We then examine the 2GC image and locate the brightest sources (>0.1 Jy) in DS9 [78], and create a bright-source region file. A minimum of three sources is tagged for each field, and a maximum of 6 sources are tagged for fields with numerous bright sources. We use the region file to divide the target field into facets equal to the number of bright sources in the region file, with the tagged sources being at the phase centre of each facet. We use KILLMS to obtain a set of direction-independent solutions for each facet separately, which is combined into a set of direct-dependent solutions for the field. We apply the COHJONES solver with a time solution interval of 5 min and frequency interval of 8 channels per solution.
The number of directions we solve for is equivalent to the number of the bright sources that are tagged per field. Finally, we use DDFACET to apply the KILLMS solutions and produce the DDEs corrected image. DDFACET also applies direction-dependent PSF deconvolution, which explicitly accounts for the time variation and bandwidth fluctuation effects. We perform the KILLMS and DDFACET loops iteratively, improving the mask and increasing the time solution intervals with each loop until we get an image with significantly reduced artefacts.
The resulting image, shown in the top right panel of Figure 3, has a rms of 0.14 mJy/beam and the beam size is 7.9 × 4.5 , p.a. 52.4 • . We note that the noise floor is higher than theCASA result. The noise for the J0034 field increases by ∼4%. However, the sources in this pipeline have a better defined structure compared to CASA. The point sources are circular and have less phase variation. This indicated that the pipeline improved the phase corrections for the fields. Although the phase corrections had improved in comparison to the CASA pipeline, the north-south artefacts were still not significantly reduced. Such artefacts would have been problematic for the point source subtraction, which is often required when extracting measurements for faint extended emission. We tried various calibration solution intervals (30 s-2 min) and facet numbers (3)(4)(5)(6)(7)(8); however, these did not improve our results. The range of calibration solution intervals produced the same results, while increasing the facet numbers resulted in some facets having higher noise levels due to fewer sources in each facet.

Experimental SPAM Wideband Pipeline
The final pipeline we explored was the experimental wideband SPAM 3 (the pipeline was accessed in October 2020). This pipeline begins by splitting the GWB data into ∼7 sub-bands of ∼30 MHz each. Thereafter, it follows the conventional GSB narrowband data reduction for the individual sub-bands (see Section 3). The calibrated sub-band visibilities are then converted into MS files and concurrently imaged using WSCLEAN. For the WSCLEAN wide bandwidth imaging step, we use a Briggs robust = 0, number of iterations = 150,000, threshold = 1 µJy, auto-mask = 9σ, amd multiscale scales of (10,20,30).
The resulting image, shown in the bottom left panel of Figure 3, has a rms of 99.0 µJ/beam and the beam size is 14.0 × 7.6 , p.a. 65.3 • . As seen in the bottom left panel of Figure 3, this pipeline produces the best phase calibration, which results in an improvement in the structure of the brightest sources. The flag percentages in Table 1 also show that the SPAM pipeline results in the least flag percentages. These low percentages result in higher sensitivity, and this is indicated by the higher peak flux recovered for the double-lobed source shown in Figure 3.  Although the noise levels are higher compared to the CASA pipeline for both images, the improvement of the phase calibrations meant we would be able to extract the bright point sources, which is a necessary step when extracting faint diffuse emission. The peak fluxes of the sources in the SPAM-reduced images were much higher in than those in the CASA-reduced and DDFACET/KILLMS-reduced images. This indicates that the phase and amplitude calibration is improved, despite the higher noise floor, so that less of the real signal is being fractured into artefacts or side lobes. However, further steps are needed to improve the algorithm so that the 3GC strategy is more robust. These steps include more comprehensive ionospheric corrections in the sub-band data reduction and primary beam correction for the full bandwidth image.

Flux Density and P 1.4GHz Calculations
We use the 2GC calibrated data to search for extended diffuse radio emission in the cluster region. To enhance the brightness of the diffuse emission, we produce a low resolution image. For the low resolution imaging process, we use the SPAM pipeline. We first image the compact sources by applying a uv-range cut at baselines ≥ 3 kλ. At this uv-range, all the point source emission is captured, while higher uv cuts result in residual point source emission and lower uv cuts also capture the extended diffuse structure. We then use the high resolution imaging step to create model visibilities that only contain point sources. For the high resolution image, we use a Briggs robust parameter of −0.8 and a uv-range 6 kλ. We subtract the point sources in the visibility plane. We then produce a full resolution point-source-subtracted image to ensure that the subtraction is successful. Finally, we image the visibilities in the uv-range ≤ 8 kλ whilst applying a uv-taper at 5 kλ and using a Briggs robust parameter of 0.8. These tapering and uv-cut values are ideal for capturing the full extension of the emission while producing an image with decent resolution and rms noise levels for the J0034 field.
We followed the same procedure for the GWB dataset. The point sources are subtracted in the sub-band datasets, and the point-source-subtracted MS files are imaged concurrently using WSCLEAN. Finally, we proceed to extract the flux density of the detected diffuse radio source. We use the statistical method derived in Knowles et al. [79] to account for the contamination from the point source subtraction procedure. We select 100 random off-source positions and create a source catalogue using the high-resolution image. We use the off-source and source positions to calculate flux densities on beam-sized regions in the low resolution image. We calculate the mean (µ) and standard deviation (σ) of the flux densities. The flux contamination resulting from the point source subtraction is the mean of the on-source positions (µ srcs ). This results in a flux measurement given by where S meas is the measured flux density and N is the number of beams within the region that the flux density is measured. We use a polygon region guided by 3σ contours to measure the flux density of the detected diffuse emission in the low resolution image. The systematic error due to point source subtraction is where σ srcs and σ o f f −src are the standard deviations of the on-source and off-source populations, respectively. We determined the uncertainty associated with the flux measurement as where δS is the calibration uncertainty (∼10% for GMRT [80]), σ rms is the rms noise of the low resolution images, measured using a background region tagged in DS9, σ syst is the systematic error from the point source subtraction, and N is the number of beams within the region of flux measurement. Since J0034 has archival multi-frequency radio observations, we compute the integrated spectral index value using S ν ∝ ν −α , where ν is the central frequency of the observation and S ν is the measured flux at that frequency. The error associated with the spectral index is calculated arithmetically, with the assumption that the flux-density measurement errors at the two frequencies are independent. We use the spectral index value to extrapolate the radio power at 1.4 GHz. We use the following equation to calculate the k-corrected radio power at 1.4 GHz where D L is the luminosity distance of a cluster at redshift z, S ν is the flux density at frequency ν, and α is the power law spectral index. We also measure the largest angular size (LAS) and the largest linear size (LLS) of the diffuse emission. The LLS is calculated by dividing the LAS by the angular size distance (D A ). In the following section, we present a case study of J0034 which was reported to host diffuse radio emission by Knowles et al. [81] (hereafter K20), and we compare the GSB and the GWB results for this cluster.

Case Study: ACT-CL J0034.4+0225
J0034 is at redshift 0.382 and has a mass of M 500c = 9.0 × 10 14 M [61]. At this cluster redshift, 1 corresponds to 5.223 kpc and the luminosity distance is 2057.4 Mpc. Carrasco et al. [82] used data from the FOcal Reducer and low dispersion Spectrograph 2 (FORS2; [83]), which is mounted at Very Large Telescope (VLT) to study the galaxy cluster. From their observations they derived that the cluster has a rest-frame velocity dispersion of 713 ± 179 kms −1 and a dynamical mass of (3.37 ± 2.19) × 10 14 M . MeerKAT L-band observations by K20 indicated that the cluster hosts a candidate radio halo. From their observations, they measured the flux density of the radio halo to be S 1.16GHz = 1.26 ± 0.2 mJy, and it has LLS of 348 kpc.
From our 325 MHz GSB data (PI: Kenda Knowles, ID: 32_016), we observe the diffuse emission for the first time at this frequency, and confirm its presence as tentatively detected at 1.16 GHz by K20 (see Figure 5). The emission we detect has a larger extent than that detected in K20 due to the observations being at a lower frequency. The flux density of the detected diffuse emission is S 323MHz = 16.77 ± 2.34 mJy, and the LAS is 2.32 × 1.46 , corresponding to a LLS of 726 kpc × 459 kpc. We note that the flux density measurements might contain residual emission from the BCG that is embedded in the radio halo. This contamination is accounted for as explained in Section 3.3. Using our results and the 1.16 GHz flux density measurement from K20, we calculated the spectral index of the radio halo and found it to be 2.26 ± 0.32. This indicates that the radio halo is an ultra steep spectrum. We used the spectral index to extrapolate the halo radio power at 1.4 GHz and found it to be (0.60 ± 0.15) × 10 24 WHz −1 . Given the extent of the diffuse emission, the central location, and its regular morphology, we confirm that it falls in the category of radio halos. We also note that the irregular morphology of the emission might be an indication of merger activity, given that this cluster is paired with ACT−CL J0034.9+0233 which is at the same redshift with a projected separation of 3.5 Mpc [61].
From the GWB data (PI: Kenda Knowles, ID: 32_016), we detect extended diffuse emission in the central region of the cluster. We measure the LAS of the diffuse emission to be 2.31 × 1.48 , corresponding to a LLS of 723 kpc × 463 kpc. This is similar is size to the GSB detection (LLS ∼ 726 × 459 kpc). The sizes and morphology of the diffuse emission are similar for both the GWB and the GSB images (See Figure 5). Hence, we retain the radio halo classification from the narrowband data analysis.
The flux density of the radio halo for the GWB observations is S 398MHz = 6.22 ± 1.64 mJy. We use this value and the 1.16 GHz flux density measurement from K20 to calculate the spectral index of the radio halo. We determine a spectral index of α 1160 398 = 1.75 ± 0.36, consistent within the uncertainties of that determined using the GSB measurement. Using our GWB flux density and measured spectral index, we determine a 1.4 GHz radio power of (0.59 ± 0.12) × 10 24 WHz −1 , in agreement with the GSB extrapolated value. We recover the radio halo at a higher SNR in the GWB data. The halo radio power is in agreement for both the GWB and GSB derived values. The morphology of the radio halo in the narrowband image is slightly elongated compared to the wideband image. The difference in morphology could be an indication that the emission at the edge of the halo is faint and hence picked up by the GWB observations, which goes down to frequencies of ∼250 MHz. For this cluster, the GWB data did not significantly add any new scientific information for the radio halo. However, we note that this dataset was severely affected by RFI and DDEs. The bright sources make this field an extreme data reduction challenge and is therefore a good candidate to test future SKA pipelines.

Summary and Future Outlook
The sensitivity of new generation telescopes has allowed for the study of diffuse emission in previously unexplored parameter spaces. Recent statistical studies of diffuse radio emission now consist of cluster samples that target lower mass and higher redshift clusters. Such studies are crucial for understanding the cosmological evolution of the diffuse radio sources, their connection to the dynamical state of the host cluster, and for the refinement of the currently existing formation theories.
In this paper, we presented an overview of an SZ-selected ACTPol sample of 40 clusters that spans a wide mass (4.5 < 10 14 M < 10.5) and redshift range (0.15 < z < 1.0). We followed up this sample using the GMRT. Seventeen clusters have existing archival GSB observation, and we obtained uGMRT GWB band 3 data for thirteen clusters. We reduced the J0034 GSB data using SPAM. We then presented the challenges we faced with reducing the GWB data for the complex J0034 field, and the three pipelines we explored. The experimental SPAM wideband pipeline produced the most improved results. Finally, we used J0034 as a case study to compare the GSB and GWB data analysis results. We found that halo radio power is in agreement for both the GWB and GSB derived values. Although the GWB data did not significantly add any new scientific features for the radio halo, the flux is detected at a higher σ level. We note that the uGMRT data reduction still needs to be refined to perform in-band spectral index studies.
We note that the current best pipeline for the GWB data is still in the experimental stage and could be optimised to produce better results. The complex J0034 field indicates the need for advanced 3GC technique for the next generation of low frequency and wide bandwidth telescopes. Finally, we aim to carry out statistical studies on the full sample once the narrowband and wideband analysis is completed.

Institutional Review Board Statement: Not applicable
Informed Consent Statement: Not applicable Data Availability Statement: The uGMRT data reported in this paper are available through the GMRT online archive. https://naps.ncra.tifr.res.in/goa. The data were last accessed November 2020. Besides the raw UVFITS data, SPAM pipeline processed GSB images are also available in the GMRT archive. The observation IDs for the gsb and gwb datasets are 8651 and 9434, respectively.