A Study of Abell 3667 with KAT-7

Early Science with the Karoo Array Telescope: a Mini-Halo Candidate in Galaxy Cluster Abell 3667

Abstract

Abell 3667 is among the most well-studied galaxy clusters in the Southern Hemisphere. It is known to host two giant radio relics and a head-tail radio galaxy as the brightest cluster galaxy. Recent work has suggested the additional presence of a bridge of diffuse synchrotron emission connecting the North-Western radio relic with the cluster centre. In this work, we present full-polarization observations of Abell 3667 conducted with the Karoo Array Telescope at 1.33 and 1.82 GHz. Our results show both radio relics as well as the brightest cluster galaxy. We use ancillary higher-resolution data to subtract the emission from this galaxy, revealing a localised excess, which we tentatively identify as a radio mini-halo. This mini-halo candidate has an integrated flux density of mJy beam at 1.37 GHz, corresponding to a radio power of P W Hz, consistent with established trends in mini-halo power scaling.

keywords:
radio continuum: general – galaxies: clusters: individual (Abell 3667) – magnetic fields – polarization
12

1 Introduction

Galaxy clusters are among the largest gravitationally-bound structures in the Universe. They form following the density profile of the dark matter distribution, from the merger of smaller galaxy groups into larger structures. As these structures grow, they undergo many dynamical interactions which make them excellent physical laboratories. Clusters grow through both accretion of gas from their environment as well as merger events with other galaxy groups and clusters. Major mergers are dramatic, releasing tremendous amounts of energy ( ergs; Ferrari et al. 2008) into the intra-cluster medium (ICM). One tell-tale indicator of clusters that have undergone recent merger events is the presence of large-scale, diffuse radio emission, examples of which are radio haloes and radio relics. Haloes are usually located in the cluster centre, whereas relics are typically observed toward the periphery at megaparsec distances from the cluster centre. Approximately 42 haloes and 36 relics3 are currently known to be hosted by galaxy clusters (Feretti et al. 2012).

Common characteristics of radio relics include low surface-brightness (typically 1 Jy arcsec at 1.4 GHz) and a steep spectral profile of the form , where is the spectral index taking typical values of . The emission from such structures is synchrotron, indicative the presence of large-scale magnetic fields and relativistic electrons in the intracluster medium (ICM, see Feretti 2005; Ferrari et al. 2008). Observations have shown that relics are highly-polarized (Feretti et al., 2012) making them good targets for studying the physics of the magnetised ICM, as the magnetic field strength can be reconstructed by, for example, rotation measure (RM) synthesis (Brentjens & de Bruyn 2005) or equipartition arguments (Beck & Krause 2005).

Observational data on relics is broadly consistent with an interpretation of the diffusive shock acceleration (DSA) model of electron acceleration (Feretti et al. 2012, and references therein). DSA models predict that the magnetic field within radio relics will be aligned with the shock front, and observations concur with this prediction (e.g. van Weeren et al. 2010; van Weeren et al. 2012). Shocks have been confirmed using X-ray data (through detected surface brightness discontinuity and/or temperature jumps) in only a handful of clusters to-date, including for example Abell 754 (Macario et al., 2011), Abell 3376 (Akamatsu et al., 2012b), and Abell 3667 (e.g. Finoguenov et al. 2010; Akamatsu et al. 2012a). Typically observed shock Mach numbers lie in the range ; from simulations we expect a Mach number approximately in the range (Feretti et al. 2012, and references therein).

Relics are also known to be hosted by some cool-core clusters which suggests both major and minor/off-axis merger events can be energetic enough to generate relics. Although they exist on vastly different scales, examples of this are the small-scale (a few kpc in extent) relic in Abell 85 (Slee et al. 2001, Schenck et al. 2014) and the giant ( Mpc size) radio relic in the merging cluster Abell 115 (Govoni et al. 2001, Gutierrez & Krawczynski 2005). Additionally, X-ray observations of Abell 115 suggest that while the (off-axis) merger event has significantly disturbed the morphology of the two subclusters, both possess intact cool cores (Gutierrez & Krawczynski, 2005).

Additionally, some cool-core clusters are known to host smaller-scale diffuse radio emission known as mini-haloes. Radio mini-haloes (MH) are diffuse radio sources on a smaller scale than the giant radio haloes (GRH) and relics seen in dynamically-active clusters (typical MH scale size kpc; GRH scale size Mpc). MH are typically co-located with the brightest cluster galaxy (BCG) and although few are known (15 confirmed; Giacintucci et al. 2014) a general scaling trend is observed between BCG radio power and MH radio power (Giacintucci et al., 2014). MH typically exhibit spectra that are slightly steeper than those of relics or giant haloes (). Two main models exist for the mechanism responsible for MH electron acceleration: magnetohydrodynamic (MHD) turbulence within the cool core and secondary electron (or ‘hadronic’) models.

Two principal sources of MHD turbulence have been suggested: (1) turbulence generated by a cooling flow (for example Gitti, Brunetti & Setti 2002) and (2) turbulence induced by gas sloshing within the cool core (ZuHone, Markevitch & Brunetti, 2011). Evidence for each source of turbulence is provided by, for example, (1) the correlation between cooling flow power and MH power (Gitti et al., 2004) and (2) the detection of spiral-shaped cold fronts in observations (for example Owers et al. 2009) and simulations (ZuHone et al., 2013). The observed trends in MH/BCG co-location and radio power suggest a potential source of the electrons responsible for the MH radio emission (Cassano, Gitti & Brunetti 2008; Giacintucci et al. 2014).

The hadronic group of models predicts that the population of relativistic electrons responsible for MH arises from interactions between cosmic ray protons (CRp) and protons of the ICM (Dennison, 1980). This group of models can reproduce some properties of MH, but purely-hadronic models generally fail to replicate a number of observed characteristics, including the radial steepening of the spectral index and the detected high-frequency cutoff in some MH spectra (Giacintucci et al., 2014).

1.1 Abell 3667

Located at a redshift of (Owers, Couch & Nulsen, 2009) Abell 3667 (hereinafter A3667) is a massive double-cluster of galaxies with (Sodre et al., 1992). It is believed to have undergone a recent major merger (Finoguenov et al., 2010) and is known to host two large-scale, diffuse radio relics (for example Röttgering et al. 1997). The brightest cluster galaxy (BCG) is head-tail radio galaxy MRC B2007-569 (per the nomenclature of Goss et al. 1982).

Optical observations have yielded evidence of substructure in A3667: the cluster is believed to consist of two main sub-clusters of 242 (164) cluster members in the primary (secondary) sub-cluster and one tertiary sub-group with 27 members (Owers et al. 2009). The primary sub-cluster density distribution lies coincident with the peak of X-ray emission; the secondary sub-cluster is to the North-West of the cluster centre and the tertiary is to the South-East. Optical observations also suggest that the merger event is occurring close to the plane of the sky (Johnston-Hollitt, Hunstead & Corbett 2008, Owers et al. 2009).

A3667 hosts two giant radio relics (Röttgering et al. 1997; Johnston-Hollitt 2003) North-West (NW) and South-East (SE) of the cluster centre. Observations show the NW relic has an integrated flux density in the range 2–4 Jy at 1.4 GHz (Röttgering et al. 1997: 2.40.2 Jy; Johnston-Hollitt 2003: 3.70.3 Jy) and a spectral index of approximately between 1.4 and 2.4 GHz (Röttgering et al., 1997). The SE relic has an integrated flux density of 0.300.02 (Johnston-Hollitt, 2003).

There is tentative evidence of a radio halo in A3667, with an integrated flux density of 33 mJy at 1.4 GHz (ATCA: Johnston-Hollitt 2003) and 446 mJy at 3.3 GHz (Parkes: Carretti et al. 2013). However, due to the difference between angular scales recovered by the two instruments, the recovered flux density at 1.4 GHz should be considered a lower limit only. In addition, recent single-dish work has suggested the presence of an unpolarized radio ‘bridge’ of synchrotron emission (from electrons re-accelerated by turbulence in the post-shock region) connecting the NW radio relic with the cluster centre, with surface brightness at 2.3 GHz (Carretti et al., 2013).

From observations at radio wavelengths, the magnetic field strength in the NW relic has been estimated. RM analysis of two point sources background to the NW relic yields estimates of G (Johnston-Hollitt, 2004); for comparison, equipartition arguments yield values of G (Johnston-Hollitt, 2004). The magnetic field strength in A3667 has also been estimated from X-ray observations: lower limits of (Nakazawa et al., 2009) and G (Finoguenov et al., 2010) in the NW relic; central magnetic field strength G (Vikhlinin, Markevitch & Murray, 2001). These estimates of the magnetic field strength are summarised in Table 1.

Jumps in X-ray hardness and surface brightness are detected at the outer edge of the NW relic, consistent with predictions based on the shock (re-)acceleration model, with a Mach number of and shock speed (Finoguenov et al., 2010). A front of cold gas is detected coincident with the BCG (Markevitch, Sarazin & Vikhlinin, 1999). Hudson et al. (2010) apply 16 cool-core diagnostics to the Chandra HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS; Reiprich & Böhringer 2002) to identify the best parameters for characterising cool-core clusters; following which A3667 is classified as a weak cool-core cluster.

Using the Murchison Widefield Array (MWA; Tingay et al. 2013) at 149 MHz, Hindson et al. (2014) recover an integrated flux density of () Jy for the NW (SE) relics, and an average spectral index of for both. From the spectral index of the NW relic, the data imply a shock Mach number of 2.4; within the range expected from DSA and consistent with the Mach number derived from X-ray observations by Finoguenov et al. (2010). The NW relic exhibits a spectral index gradient that steepens toward the cluster centre, indicative of an ageing electron population; the radial profile is consistent with that of an expanding shock travelling out from the cluster centre, as has been observed in a number of other relics (for example the relic in the ’Sausage cluster’; van Weeren et al. 2010). Additionally, Hindson et al. recover flux from both the radio bridge region and radio halo at lower significance than expected from previous work; the interpretation is that previous detections of a bridge and halo in A3667 are the result of poorly-sampled Galactic emission that contaminates the region.

Paper Estimation method Region B
()
Vikhlinin et. al. (2001) X-Ray (Chandra) Centre 7.0-16.0
Johnston-Hollitt (2004) RM estimates NW relic 3.0-5.1
Johnston-Hollitt (2004) Equipartition NW relic 1.5-2.5
Nakazawa et al. (2009) X-Ray (Suzaku) NW relic
Finoguenov et al. (2010) X-Ray (XMM) NW relic
Table 1: Summary of previous magnetic field strength derivations for A3667.

In this work, A3667 was observed with the Karoo Array Telescope, a precursor to MeerKAT and the Square Kilometre Array, to investigate the nature of the diffuse radio emission from the relics and ‘bridge’ of the cluster. This paper is structured as follows: we discuss the Karoo Array Telescope, the observations, and the reduction process in §2. Our results are presented in §3 and in §4, while the implications are discussed in §5. We draw our conclusions in §6. Throughout this work, we assume a concordance cosmology of H km s Mpc, , . All errors are quoted to 1-. We adopt the spectral index convention that and we take the redshift of A3667 to be 0.0553 (Owers et al. 2009). At this redshift, an angular distance of 1 corresponds to 1.031 kpc (Wright, 2006).

2 The Karoo Array Telescope

The Karoo Array Telescope (KAT-7) is a synthesis telescope consisting of seven 12-m diameter dishes with baselines in the range  m, located in the Republic of South Africa. It has dual linear-feed receivers capable of observing in the range 1.20–1.95 GHz, with an IF bandwidth of 400 MHz. However, due to analog filters in the IF and baseband system, only the central 256 MHz of the IF bandwidth is useable.

This work was conducted as part of KAT-7 commissioning; as such central observing frequencies were selected to be 128 MHz from the bottom and top end of the observable bandwidth, yielding useable slices that cover the extremes of the KAT-7 frequency range. Our reference frequencies for this work are therefore 1382 and 1826 MHz. With a longest baseline of 819 (1124)  at 1328 (1822) MHz, KAT-7 has modest angular resolution of  arcmin, beneficial to recovering diffuse low surface-brightness emission.

KAT-7 is capable of observing using a number of correlator modes (see http://public.ska.ac.za/kat-7 for further details). In this work, wideband mode was used, with the IF divided into 1024 channels; after removal of the bandpass edges, the central 601 channels were selected for a useable bandwidth of 235 MHz. Table 2 provides a summary of telescope characteristics. Although primarily a testbed for new materials and engineering for the upcoming MeerKAT and Square Kilometre Array (SKA) projects, early science is underway with KAT-7 (see for example Armstrong et al. 2013).

Number of Antennas 7
Dish Diameter (m) 12.0
uv-min (m) 25.6
uv-max (m) 185.0
System Temperature (K) 30
IF Bandwidth (MHz) 400.0
Useable Bandwidth (MHz) 256.0
Effective Bandwidth in this Work (MHz) 157.0 235.0
Effective Number of Channels in this Work 401 601
Polarization XX,YY,XY,YX
Reference Frequency (MHz) 1328 1822
Effective Reference Frequency (MHz) 1372 1826
Resolution (arcmin) 4.19 3.06
Largest detectable size (arcmin) 30.31 22.09
  • Reduced bandwidth due to RFI environment – see §2.2.1 for further details.

Table 2: Technical Details of KAT-7

2.1 Observation Details

Kat-7

A3667 was observed in full-polarization in both upper and lower bands. Observations at 1822 MHz were conducted in 2013 March, and 1328 MHz data was taken during 2013 May. Between these observing periods, KAT-7 was upgraded with new receivers to reduce intrinsic radio frequency interference (RFI). All results in this paper were produced with a single observing run at 1822 MHz, and two runs at 1328 MHz. Table 3 summarises the details of observations. At 1822 MHz, 9 close-packed pointings were used to cover the cluster region, while 7 pointings were used at 1328 MHz. The KAT-7 primary beam has full-width at half maximum (FWHM) of approximately at 1822 MHz. The uv-coverage of the central pointing on A3667 at each frequency is shown in Figure 1.

{subfigure}

[t]0.45 {subfigure}[b]0.45

Figure 1: uv-coverage of A3667. top: 1822 MHz, bottom: 1328 MHz. Both figures plot the uv-coverage of the central pointing in the mosaic.

During each observing run, three sources were used as calibrators: PKS B1934-638, 3C138/3C286 and PMN J2040-5735. PKS B1934-638 (hereafter 1934-638) was used as both a bandpass calibrator and flux calibrator, 3C138 (3C286) was the polarization calibrator at 1822 (1328) MHz, and PMN J2040-5735 was to be phase calibrator. However, the local RFI environment of PMN J2040-5735 proved too strong to derive good calibration solutions, so 1934-638 was substituted as phase calibrator, as it was interleaved sufficiently often to be suitable. In the lower band, 1934-638 was interleaved every 30 minutes. In the upper band, it was interleaved every 60 minutes. Phase solutions in both bands were smoothly-varying, with less than variation over the entire observation, and variation of less than 3 between scans.

Thermal Measured
Pointing ID RA Dec
(J2000) (J2000) (s) (Jy beam) (mJy beam)
ACO3667-H01 20h11m16.s60 -56\textdegree5957.00 4110 61.18 1.38
ACO3667-H02 20h11m16.s90 -56\textdegree4957.00 4280 59.96 1.35
ACO3667-H03 20h11m17.s19 -56\textdegree3957.00 4230 60.31 1.43
ACO3667-H04 20h12m29.s99 -56\textdegree5958.99 4210 60.45 0.94
ACO3667-H05 20h12m29.s99 -56\textdegree4958.99 4250 60.17 1.14
ACO3667-H06 20h12m29.s99 -56\textdegree3958.99 4310 59.75 1.60
ACO3667-H07 20h13m43.s39 -56\textdegree5957.00 4240 60.24 1.52
ACO3667-H08 20h13m43.s10 -56\textdegree4957.00 4130 61.04 1.18
ACO3667-H09 20h13m42.s80 -56\textdegree3957.00 4190 60.60 1.60
ACO3667-L01 20h12m06.s13 -56\textdegree4912.30 5380 67.23 2.53
ACO3667-L02 20h11m02.s14 -56\textdegree4912.30 5550 66.19 1.99
ACO3667-L03 20h11m34.s14 -56\textdegree3520.90 5480 66.61 1.85
ACO3667-L04 20h12m29.s56 -56\textdegree3520.90 5310 67.67 2.79
ACO3667-L05 20h13m10.s14 -56\textdegree4912.30 5360 67.35 2.38
ACO3667-L06 20h12m29.s56 -57\textdegree0303.70 5420 66.98 1.93
ACO3667-L07 20h11m02.s14 -57\textdegree0303.70 5220 68.25 2.54
Table 3: KAT-7 pointings on A3667 at 1822 MHz (ACO3667-H*) and 1328 MHz (ACO3667-L*), with their individual integration times and both expected thermal noise component and measured rms based on individual images made for diagnostic purposes.

It has been suggested that 3C138 exhibits variation in its polarized flux emission, resulting from a flare that began in 2003 (Perley & Butler, 2013) so 3C286 was chosen as polarization calibrator at 1328 MHz. The polarization calibrator was observed for 4 scans during each observing run, for a total integration time of order 3000s.

2.2 Reduction Process

RFI Excision

Initial pre-processing was performed on the data with an in-house automatic RFI flagging pipeline (priv. comm., Mauch T.). Subsequently, the data were visually inspected to remove strong sources of RFI, and then automatically flagged using flagdata4 in casa 4.15.

For the upper band ( MHz) observations, KAT-7 was equipped with its first-generation receivers; both low-level and sharp RFI peaks were present in the data. Additionally, the first-generation receivers experienced strong, broad-band internal RFI – an issue which has since been corrected with the upgrade to the second-generation receivers. Following visual inspection approximately 50 per cent of the data was flagged out across the entire bandwidth, resulting in an effective reference frequency of 1826 MHz.

Prior to the lower band ( MHz) observations being conducted, KAT-7 was upgraded with improved second-generation receivers to eliminate instrumental RFI. In the lower band highly-significant RFI from a known source dominated the signal below approximately 1280 MHz. Hence a cut was made at 1294 MHz, resulting in an effective bandwidth of 157 MHz and effective reference frequency (adopted henceforth) of 1372 MHz. Aside from this cut, approximately 16 per cent of the data was flagged across the bandwidth. These effective reference frequencies and effective bandwidths are also summarised in Table 2.

Calibration

The results presented in this work have been calibrated using casa 4.1. The flux scale was set using the primary calibrator (1934-638) on the Perley-Butler (2010) model which yields a Stokes I flux density of 14.795 Jy at 1450 MHz and 13.159 Jy at 1944 MHz. Parallel-hand delays, complex gains and bandpass solutions were derived for each antenna using 1934-638, following standard practice. For calibration, 1934-638 was assumed to have zero polarization; although 1934-638 is known to have no linear polarization, it does have low-level circularly-polarized emission ( per cent; Rayner, Norris & Sault 2000). This assumption of zero polarization will lead to a slight offset in Stokes V, which should be corrected if analysis of that parameter is to be made. However, in this work we neglect this offset.

Cross-hand delays and phase solutions were determined using the polarization calibrator (3C286 at 1328 MHz; 3C138 at 1822 MHz) assuming a non-zero polarization model. Stokes parameters were derived directly from the measured complex visibilities and then compared with the source model from Perley & Butler (2013) (see Table 4) in order to correct for the 180 ambiguity present in the XY phase solution6. With these additional calibration solutions, the complex gains were solved for the polarization calibrator using the full-Stokes source model from Perley & Butler (2013) and used to calibrate the polarization leakage. Typical leakages were of the order of 3 per cent in the upper band and 3–6 per cent in the lower band.

From these data we recover a polarization fraction of per cent and electric vector position angle (EVPA) of for 3C286; for 3C138 the polarization fraction is per cent and EVPA , in reasonable agreement with those determined by Perley & Butler (2013) for these sources. See Table 4 for details.

Source Freq. Fractional Pol. EVPA
(MHz) () ()
3C286 Model 1050 8.6 33
1450 9.5 33
1640 9.9 33
1950 10.1 33

Recovered 1372
3C138 Model 1050 5.6
1450 7.5
1640 8.4
1950 9.0

Recovered 1826
Table 4: Polarization properties of 3C286 and 3C138 from Perley & Butler (2013) and those derived from data taken with KAT-7 during these observations.

Imaging

Deconvolution was performed in all Stokes parameters using multi-scale clean in casa. The deconvolution process used natural weighting, with scales set to [0 (point sources), 1beam, 2beam] to improve sensitivity to diffuse emission. Maps were produced at reference frequencies of 1826 and 1372 MHz. The fields were mosaicked during deconvolution, cleaning initially down to a threshold of 1.5 (where is the theoretical confusion limit – see section §2.3) before inspection. Subsequently, further cleaning was performed down to a threshold of . Finally, the images were corrected for the primary beam response. In both cases, the cluster was imaged out to 10 per cent of the primary beam response. See Figures 2, 3 and 4 for the maps.

2.3 Confusion and Thermal Noise

At the angular resolution of KAT-7, we expect our results to be limited by confusion noise rather than thermal noise in Stokes I. KAT-7 has a system temperature of K and we expect thermal noise values of Jy at 1822 (1328) MHz for the full observation. Pointings at each frequency, integration times and both predicted thermal noise and measured rms per pointing are listed in Table 3. The given rms values are measured from images of each pointing made individually for diagnostic purposes.

We derive the confusion limit using a 1.4 GHz differential source counts model which unifies contributions from flat-spectrum radio quasars, BL Lac objects and steep-spectrum objects (de Zotti et al., 2010). Defining confusion as when a level of one source per 20 beam areas is reached, we interpolate this 1.4 GHz model to find the corresponding flux density cutoff for confusion

(1)

where is the number of sources per beam area, is the flux density cutoff, is the beam area in steradians and is the differential source counts model (in Jy sr). We subsequently use to integrate the differential source counts model to derive the confusion limit

(2)

where is the confusion limit and is the source flux density. KAT-7 has beam size 3.05 (4.19) at 1822 (1328) MHz, therefore we derive values of = 1.28 (2.15) mJy beam. These confusion limit values dominate over the predicted thermal noise. We also note that KAT-7 rapidly becomes confusion-limited, as the thermal noise drops below the confusion limit after s at 1822 MHz and s at 1328 MHz.

2.4 Archival ATCA Data

Archival Australia Telescope Compact Array (ATCA) data was used to provide a model of the head-tail radio galaxy MRC B2007-569 (Goss et al. 1982; hereafter B2007-569) for subtraction purposes. B2007-569 was observed with the ATCA in 1.5D configuration, with baselines in the range 107 m – 4.4 km, providing a nominal resolution of 10 arcsec and maximum detectable scale size of 7 arcmin at 1.4 GHz. Observations were taken on 2011 November 26 (project code CX224, P.I. Carretti E.) for a total integration time of 9 h.

The data were calibrated with the Multichannel Image Reconstruction, Image Analysis and Display (miriad; Sault, Teuben & Wright 1995) package following standard techniques. The flux density scale was set using the Reynolds (1994) model for 1934-638, which yields flux densities consistent with those set by the casa task setjy. Therefore we can be confident that both the KAT-7 and ATCA observations will have consistent flux scaling. Following calibration, the full 2 GHz Compact Array Broadband Backend (CABB; Wright et al. 1994) bandwidth was divided into subbands of 200 MHz width. The calibrated data were imported into casa for further flagging and initial imaging using the task clean, followed by 5 rounds of phase-only self-calibration and imaging. Images were made in all subbands above 1180 MHz to enable modelling of the spectral variation of B2007-569. The 1180 MHz band was excluded as it was severely flagged (due to RFI excision and edge channel removal) resulting in a useable frequency coverage of 1288 – 3124 MHz.

To obtain a clean-component model (CCM) for subtraction from our KAT-7 data, we performed multi-frequency synthesis using the msmfs algorithm (Rau & Cornwell, 2011) as implemented in the casa task clean. This algorithm uses the large fractional bandwidth available to deconvolve the image using multiple spatial scales and Taylor coefficients (as a function of frequency) to describe the intensity distribution. After deconvolution, a spectral index map and its associated error map are produced. We performed mfs-clean with the reference frequency set to 1400 MHz in order to provide a model for subtraction from the KAT-7 data most sensitive to large-scale structure. See §4.5 for further details on the subtraction, maps and analysis.

3 Results

We present total intensity (Stokes I) images at 1826 MHz (Figure 2) and 1372 MHz (Figure 3). Throughout all figures, total intensity contours mark [-1, 1, 2, 3, 4, 6, 8, 12, 16, 32, 64, 128] x unless otherwise stated.

Figure 2 has image noise mJy beam, a value comparable to the confusion limit at this frequency (see §2.3). Large-scale diffuse emission from both radio relics is clearly detected, as are many point sources in the field of view, and the head-tail radio galaxy B2007-569. With the resolution of KAT-7, we do not resolve the structure of B2007-569. The NW relic extends for an angular size of 33 arcmin, corresponding to a physical distance of 2.04 Mpc at this redshift; the SE radio relic covers an extent of 23 arcmin or 1.42 Mpc. Figure 3 has image noise level mJy beam, comparable to the confusion limit derived in §2.3. With the marginally-lower resolution, the diffuse sources exhibit smoother morphologies than in Figure 2.

We present linearly-polarized images at 1826 and 1372 MHz in Figure 4. Contours are Stokes I at the corresponding frequency, and vectors display the polarization angle of emission from the relics for regions where the flux density exceeds a 10 threshold, where mJy beam at 1826 (1372) MHz.

The polarisation vectors have been both rotated by 90\degree  to represent the direction of the magnetic field and corrected for Galactic Faraday rotation using the map of Oppermann et al. (2012). We estimate the Galactic Faraday depth contribution in the direction of the NW (SE) relic to be () rad m; we correct by subtracting a rotation of 8.17\degree (5.14\degree) at 1826 MHz and 14.47\degree (9.11\degree) at 1372 MHz.

Figure 2: Total intensity map of A3667 at 1826 MHz. Image rms is 1.30 mJy beam. Contours mark [-1, 1, 2, 3, 4, 6, 8, 12, 16, 32, 64, 128] x . Beam size (unfilled circle, bottom-left) is arcsec.
Figure 3: Total intensity map of A3667 at 1372 MHz. Image rms is 2.18 mJy beam. Contours mark [-1, 1, 2, 3, 4, 6, 8, 12, 16, 32, 64, 128] x . Beam size (unfilled circle, bottom-left) is arcsec.
{subfigure}

[t]0.745 {subfigure}[b]0.7

Figure 4: Linearly-polarized maps of A3667 at 1826 (top) and 1372 (bottom) MHz. Contours are Stokes I at the corresponding frequency, as per Figures 2 and 3. Vectors are the polarization angle rotated by 90 to represent the direction of the magnetic field, derived where the polarized flux density exceeds 10 and (1.25) mJy beam at 1826 (1372) MHz.
Integrated Flux Density
KAT-7 SUMSS
SUMSS Catalog Name RA DEC 1826 MHz 1372 MHz 843 MHz C13 ID Notes
(J2000) (J2000) (mJy) (mJy) (mJy)
SUMSS J200930–570058 20h09m30s 82.44.8 93.26.3 127.24.0 - -
SUMSS J200940–565354 20h09m40s 38.33.1 37.33.1 50.21.8 13 -
SUMSS J201127-564358 20h11m27s 471.223.7 589.129.9 1030.8 1 B2007-569
SUMSS J201155–571655 20h11m56s 80.95.3 76.45.9 66.32.2 - -
SUMSS J201309–560842 20h13m10s - 129.08.3 142.84.5 - Index is
SUMSS J201518–563440 20h15m16s 49.73.4 149.610.5 81.82.7 - -
SUMSS J201524–563805 20h15m21s 38.34.2 58.92.0 - -
SUMSS J200925–563327 20h09m25s 54.8 74.5 108.0 5 Within NW relic.
SUMSS J201330–570552 20h13m32s 50.3 73.9 122.13.8 - Within SE relic.
  • This column lists source ID numbers from Table 1 of Carretti et al. (2013) where catalogued.

  • This source appears embedded in a region of extended emission at 1826 MHz (see Figure 2) so this spectral index should be considered an upper limit.

  • The SUMSS catalog lists two components for B2007-569; we quote the sum of the flux densities from both and derive the integrated spectrum from this value.

  • Primary beam response 50 per cent.

  • Due to the lower resolution at 1372 MHz, these sources cannot be resolved separately.

  • Fits could not converge, so we present peak flux densities read directly from the image.

Table 5: Source population of the A3667 field as observed with KAT-7. Corresponding integrated flux densities at 843 MHz have been retrieved from the SUMSS catalogue (Mauch et al., 2003) for each source. RA and DEC values quoted from the fit to the 1826 MHz map. Flux densities were only fitted to sources where the primary beam response is per cent and source flux density rises to S.
upper panel: Field sources. lower panel: Sources present within extended emission from relics.

4 Analysis

4.1 Source Identification

In these data, we have shown strong detections of both large-scale diffuse emission and compact sources at both frequencies. The two relics exhibit similar morphologies to those presented in Röttgering et al. (1997) and Johnston-Hollitt (2003). Additionally, we detect many of the brighter sources seen in other work, as revealed by (for example) comparison with the Sydney University Molonglo Sky Survey (SUMSS; Bock, Large & Sadler 1999) catalogue. With the resolution available to KAT-7, we have angular scale sensitivity that is intermediate to that of Röttgering et al. (1997) and Carretti et al. (2013).

The source populations of Figures 2 and 3 were catalogued and had positions and flux densities extracted using the casa task imfit7. A catalogue of sources common to both figures is presented in Table 5, containing integrated flux densities from both KAT-7 maps as well as the integrated flux density of the corresponding source in the SUMSS catalogue (Mauch et al., 2003). The spectral index was fitted using the values at 1826 and 843 MHz; values derived are consistent with non-thermal extragalactic source populations.

Integrated flux densities for each relic were computed using the fitflux flux density fitting routine (Green, 2007) deriving values of () Jy for the NW (SE) relic. At 1372 MHz, the NW (SE) relic has an integrated flux density of () Jy – see Table 6 for a summary. All integrated flux measurements were conducted after subtracting contributions from point sources identified by Carretti et al. (2013) – see §4.5 for details on the subtraction process.

With the resolution of KAT-7, the compact sources and extended relic in the SE form one structure. The compact sources are PMN J2014-5701 (a head-tail radio source) and SUMSS J201330-570552 (a suspected FRII; Johnston-Hollitt 2003). With the resolution of ATCA however, previous work (Johnston-Hollitt 2003, 2004) has resolved the separate features in the SE relic. Hence for the KAT-7 integrated flux density, the measurements were restricted to the elongated outer part of the structure.

The integrated flux density of the NW relic is consistent with the work of Röttgering et al. (1997) where an integrated flux of 2.4 Jy at 1.4 GHz was derived, although there is discrepancy with the integrated flux density of 3.7 Jy derived by Johnston-Hollitt (2003). The negative bowl around the NW relic in Figure 3 suggests some missing flux on the largest angular scales, which may partially account for this discrepancy. For the SE relic, the integrated flux densities of () Jy at 1826 (1372) MHz agree reasonably well with the integrated flux density of 0.300.02 Jy from Johnston-Hollitt (2003) at 1.4 GHz.

4.2 Polarized Structure

As can be seen in Figure 4, we detect linearly-polarized emission from both relics at each frequency. There is little ambient emission detected elsewhere at 1826 MHz, whereas at 1372 MHz there is a greater level of polarized emission outside the relics. We attribute this higher level of ambient emission to off-axis leakage effects. The linear polarization vectors in Figure 4 have been rotated by 90 to represent the direction of the magnetic field.

From Figure 4, at 1372 MHz the magnetic field vectors lie approximately along the long axis of the SE relic, a result consistent with DSA theory and observations (e.g. van Weeren et al. 2010; van Weeren et al. 2012). For the NW relic however the picture is harder to interpret as a result of the contorted morphology. Nevertheless, some departure from the expected alignment is observed, which may indicate the merger axis is slightly tilted, similar to the case of PSZ1 G096.89+24.17 (de Gasperin et al., 2014). In such a case, the ICM would produce an observable rotation in the polarization angle measured in the more distant relic, but would not be seen in the closer relic. This deviation is not inconsistent with optical observations that suggest a merger event close to the plane of the sky (Johnston-Hollitt et. al 2008; Owers et al. 2009)

At 1826 MHz however, the magnetic field vectors lie approximately perpendicular to the long axis of the relics, suggesting they are experiencing Faraday rotation. Given the relatively small frequency difference between the datasets, an rotation of approximately 60 in magnetic field direction might imply a high rotation measure (RM). However, although there is low fractional bandwidth difference between the two maps, in wavelength-squared-space the difference is more apparent. The magnetic field vector rotation of implies a RM of rad m, accounting for the ambiguity in RM measurements (see for example Heald 2009). Such a value is consistent with previous RM estimates for the NW relic (Johnston-Hollitt, 2004).

When deriving polarization fractions for the relics, fitflux was used to determine the integrated polarized flux. At 1826 MHz we recover an integrated flux density of () mJy for the NW (SE) relic; corresponding to a polarization fraction of () per cent. At 1372 MHz we recover an integrated flux density of () mJy for the NW (SE) relic; a polarization fraction of () per cent. For comparison, Carretti et al. (2013) derived fractional polarization values of 8.5 (14) per cent for the NW (SE) relic at 2.3 GHz. However, they note that – as a result of the resolution of the single-dish Parkes radio telescope – the NW relic proved difficult to constrain, as the polarized emission blends smoothly with the weaker background emission. Johnston-Hollitt (2003) measured a polarization fraction of 10–40 per cent for the NW relic at 1.4 GHz with ATCA.

Additionally, we detect no polarized emission from the bridge region above the noise level. This is consistent with the observations of Carretti et al. (2013) and the suggestion that turbulence generated in the wake of the passing merger shock would disrupt the ordered structures and large-scale alignment necessary for polarized emission on these scales.

(MHz) 1372 1826
NW relic S(I) (Jy)
S(P) (Jy)
P/I (per cent)
SE relic S(I) (Jy)
S(P) (Jy)
P/I (per cent)
(arcsec) 251 180
(I) (mJy beam) 2.18 1.30
Table 6: Integrated flux density of the total intensity and linearly-polarized intensity of diffuse radio emission from the relics of Abell 3667 at 1826 and 1372 MHz. Also listed are the corresponding sensitivity and resolution at each frequency. Note that we measure the Stokes I integrated flux density after performing point source subtraction as described in §4.5.

4.3 Spectral Index Profile

From Figures 2 and 3 is is possible to map the spectral index variation between 1826 and 1372 MHz. To improve the fractional bandwidth over which the spectrum is mapped, however, we employ ancillary data from SUMSS to map the spectral index between 1826 and 843 MHz. We initially convolve the SUMSS map to match the resolution of KAT-7 at 1826 MHz ( arcsec) before masking emission from the relics at 10 mJy beam () to mitigate contributions from noise.

The resultant spectral index map (along with the associated uncertainty) is shown in Figure 5. For the NW relic, the spectral index is typically in the range to , marginally flatter than the profile derived by Röttgering et al. (1997) where , which may indicate curvature in the spectrum. The spectrum steepens to toward the cluster centre, flattening to approximately toward the leading edge of the relic; such variation may indicate differences in the ages of the electrons responsible for the radio emission.

{subfigure}

[tl]0.47   {subfigure}[tr]0.46

Figure 5: Spectral index profile (left) and associated uncertainty (right) of the radio relics of A3667 between 843 and 1826 MHz (36 and 16cm) at a resolution of arcsec (resolution of KAT-7 at 1826 MHz). Contours are full-resolution Stokes I at 843 MHz, from SUMSS, starting at 3 mJy beam and scaling by a factor of 2. (A colour version of this Figure is available in the online journal).

While the spectrum inverts for the Eastern spur of the NW relic, this is likely due to differences in sensitivity to faint diffuse emission resulting from resolution differences between SUMSS ( arcsec with sensitivity 1 mJy beam) and KAT-7 ( arcsec with sensitivity 1.30 mJy beam). The resulting difference in sensitivity to large-scale diffuse emission is visible in the breaking up of the SUMSS contours in this region (see Figure 5).

The SE relic exhibits a flatter spectrum, typically in the range to . This is consistent with DSA being the acceleration mechanism (the flattest spectrum possible according to DSA is ) and in agreement with earlier evidence of DSA from the alignment of the magnetic field vectors with the long axis of the relic (see Figure 4 and §4.2). Both spectral index ranges are broadly consistent with the spectral index variation derived by Johnston-Hollitt (2003) although with our moderate-resolution data we do not resolve the spectral index structure of the SE relic. Given that MOST data at 843 MHz are common to both this work and Johnston-Hollitt (2003) such consistency is not unexpected.

Hindson et al. (2014) find integrated spectral index values of for both the NW and SE relics. From our KAT-7 map, the spectral index of the NW is typically , consistent with the work of Hindson et al. Additionally, Hindson et al. find the 843 MHz flux of the NW relic inconsistent with their 4-point spectral index map; their interpretation is that MOST may be missing some flux on scales of order 30 arcmin, comparable to the size of the NW relic. Such missing flux may be responsible for the spectral index shape recovered previously. With KAT-7, we are sensitive to structures up to 22 arcmin in size; we cannot discount the possibility that some flux is missing on the largest angular scales. For the SE relic, we find a typical spectral index of , inconsistent with the integrated spectrum derived by Hindson et al.

4.4 Equipartition Magnetic Field

The magnetic field strength of A3667 can be estimated using the equipartition field strength equation from Beck & Krause (2005). The equipartition magnetic field strength is given by

(3)

where is the spectral index, is the synchrotron intensity at frequency ; is the proton rest energy and is the ratio of the number densities of protons and electrons (Beck & Krause 2005). , and are constants given by the following:

(4a)
(4b)
(4c)
(4d)

where is the spectral index of the electron energy spectrum, related to the spectral index by ; is the inclination angle with respect to the sky plane and is assumed to be constant for to be valid as presented (Beck & Krause, 2005); m is the mass of the electron and c is the speed of light. A value of denotes a face-on view, which is a reasonable approximation for A3667 as the merger event is thought to be occurring approximately in the plane of the sky (for example Finoguenov et al. 2010). Typical ratios of proton-to-electron number density are (Pfrommer & Enßlin, 2004); we take a mid-range value of for our derivation. is the line-of-sight path length, which must be assumed for the calculation; we assume a cylindrical geometry, measuring a path length of 720 arcsec for the NW relic and 350 arcsec for the SE relic. Substituting in Equation 4b, we derive:

Figure 6: Estimated equipartition magnetic field strength in the radio relics of A3667. Contours are KAT-7 Stokes I at 1826 MHz as per Figure 2. Estimates for the SE relic are largely contaminated by the background FRII. (A colour version of this Figure is available in the online journal).

The resultant map of equipartition magnetic field strength () is presented in Figure 6. Large regions have no estimate due to the fact that Equation 3 diverges for . For the SE relic, the regions where estimates exist are largely contaminated by flux from either of the secondary sources8 which have previously been resolved separately (e.g. Johnston-Hollitt 2003, 2004). Therefore we restrict our discussion of the equipartition magnetic field strength to the NW relic.

For the NW relic, the flat spectrum of the Eastern spur of this relic prohibits estimation of . However, the estimates for the remainder of the relic yield values of in the range ; consistent with equipartition estimates from previous work (e.g. Johnston-Hollitt 2004; see Table 1).

4.5 Radio Bridge

With the modest resolution of KAT-7, we are able to probe the diffuse low surface-brightness emission from A3667 on scales up to 22 (30) arcminutes at 1826 (1372) MHz. Given the extent suggested for the proposed radio bridge (20 arcmin; Carretti et al. 2013) KAT-7 should possess sufficient sensitivity to recover this feature. However, as seen in Figure 2, only patchy low-level flux is recovered between the central galaxy and the NW relic; this emission is coincident with a number of point sources that appear in the SUMSS catalogue (Mauch et al., 2003).

Extrapolating from the surface brightness of the bridge as presented in Carretti et al. (2013) of 0.21 – 0.25 Jy arcsec at 2.3 GHz, using the spectral index of derived therein, we calculate 1826 MHz (1372 MHz) surface brightness in the range 0.29 – 0.33 (0.39 – 0.45) Jy arcsec. From these surface brightness values, we should expect to recover flux density of the order of 12 – 14 (35 – 40) mJy beam with KAT-7, corresponding to the order of 10 (17) – a level not recovered in Figures 2 and 3. At low frequency (226 MHz) Hindson et al. (2014) recover a surface brightness from the bridge region that is below the significance expected from extrapolation of the surface brightness as presented by Carretti et al. (2013).

In order to further investigate the nature of the radio ‘bridge’, we extrapolate the flux densities of the 13 sources above 3 mJy identified by Carretti et. al. (2013; Table 1) to the observing frequencies covered in this work and subtract these sources from the data. All sources are modelled as a simple power-law fit, with the exception of head-tail radio galaxy B2007-569, which requires more careful modelling.

Freq. Telescope Flux Density Flux Scale Scale Factor
(MHz) (mJy)
120 MWA 7300 BT
149 MWA 6500 BT
180 MWA 5100 BT
226 MWA 4200 BT
408 MOST 1920 W69 0.97
843 MOST 1031 (c)B77 n/a
1372 KAT-7 589 B77 n/a
1383 ATCA 497 B77 n/a
1588 ATCA 416 B77 n/a
1793 ATCA 365 B77 n/a
1826 KAT-7 471 B77 n/a
2202 ATCA 304 B77 n/a
2300 ATCA 280 B77 n/a
2407 ATCA 284 B77 n/a
2817 ATCA 246 B77 n/a
3022 ATCA 229 B77 n/a
3300 ATCA 219 B77 n/a
4850 Parkes 173 (c)B77 n/a
Flux scale notation:
BT: Bootstrap method using known sources in the field – see
      Hindson et al. (2014) for details.
W69: Wyllie (1969)
B77: Baars et al. (1977)
(c) implies the flux scale is consistent with the given scale.
  • Hindson et al. (2014)

  • Molonglo Reference Catalog; Large et al. (1981)

  • SUMSS; Mauch et al. (2003)

  • This work

  • Carretti et al. (2013)

  • Parkes-MIT-NRAO (PMN) Survey; Wright et al. (1994)

  • These error values taken as 5 percent of the integrated flux density.

Table 7: Integrated flux density measurements for B2007-569 from the literature and derived in this work, with their respective flux scales and the scale factor required to convert to the scale of Baars et al. (1977). Integrated flux density is quoted before any scaling factor is applied.

Table 7 presents the integrated flux density of B2007-569 between 408 and 4850 MHz, both derived in this work and sourced from the literature. The flux densities have been measured using a number of different instruments (listed in Table 7) sensitive to emission on several different scales. The ATCA flux densities derived from our analysis of the archival CABB data were obtained by running fitflux on images made in each 200 MHz subband. We include the measured MWA flux densities from Hindson et al. (2014) in Table 7 although we omit them from the fitting as uncertainties remain due to inconsistencies between the analytic beam model (Williams et al., 2012) and beam response as revealed during commissioning (Hurley-Walker et al., 2014).

With the exception of the Molonglo Reference Catalog (MRC; Large et al. 1981) flux density, all values in Table 7 were derived in work where the flux scale was either set according to the Baars et al. (1977) scale or shown to be consistent with it. The MRC employed the Wyllie (1969) flux scale. Baars et al. (1977) note that the Wyllie (1969) scale overestimates the flux density by 3 per cent compared to the Baars et al. (1977) scale; hence this can be accounted for by a scaling factor 0.97.

The SUMSS flux scale was set through both measurements at 843 MHz and interpolation between the MRC at 408 MHz and the Parkes catalog at 2700 MHz (Bock, Large & Sadler, 1999) following the method of Hunstead (1991). This flux scale is primarily tied to the historically-determined flux density of PKS B1934-638 (among others) using the Reynolds (1994) scale, which broadly agrees with the Baars et al. (1977) scale (Reynolds, 1994). Hence no scaling factor is necessary for the SUMSS flux densities used in this work. The PMN catalog flux density values have been shown to be consistent with the Baars et al. (1977) scale (Wright et al., 1994).

We present the scaled flux densities in Figure 7. From inspection, the integrated flux density of B2007-569 appears to follow a power-law between 408 and 1826 MHz; we derive a spectral index of . While the integrated flux densities derived from our analysis of the ATCA data agree well with those of Carretti et al. (2013) the flux densities derived from the MOST and KAT-7 data are in excess of the flux densities derived using the ATCA.

The flux density of B2007-569 as measured with ATCA in this work is well-described by a power-law fit with spectral index . We also fit a polynomial of the form from which we derive a spectral index of the form where is measured in MHz. Hindson et al. (2014) find a spectral index fit of between 120 and 4850 MHz, consistent with the spectral index values derived for B2007-569 from this work.

{subfigure}

[t]0.9 {subfigure}[b]0.9

Figure 7: Integrated spectrum of B2007-569 between 408 MHz and 4850 MHz. Both panels present identical data with different fits shown. Top panel presents the power-law fits to the MOST/KAT-7 data between 408–1826 MHz (dot-dashed, ) and ATCA data of Carretti et al. (2013; solid line, ). Lower panel presents fits to the ATCA data as analysed in this work. Dashed slope (dotted line) is a power-law (second-order polynomial) fit to the ATCA data (hollow triangles) analysed in this work, with the power-law fit yielding a spectral index . Also presented is the fit in the range 408–1826 MHz for reference.
{subfigure}

[l]0.45 {subfigure}[r]0.45 {subfigure}[c]0.45

Figure 8: Stokes I
Figure 9: Spectral Index
Figure 10: Spectral Index Error
Figure 11: ATCA maps of B2007-569 produced using MFS-clean on the CABB band in the frequency range 1288–3124 MHz. top left: Stokes I image at a reference frequency of 1400 MHz. top right: Spectral index. lower: Spectral index error. Contours are the Stokes I ATCA data starting at 1 mJy/beam and scaling by a factor 2. Beam size is arcsec. mJy beam. (A colour version of this Figure is available in the online journal).

B2007-569 is not point-like at the resolution of the instruments considered here. We use the archival ATCA data to obtain a higher-resolution clean-component model (CCM) for subtraction. The wide bandwidth available with CABB yields suitable frequency coverage for multi-frequency synthesis (MFS; Rau & Cornwell 2011) clean in casa. To obtain the CCM, we performed mfs-clean using two Taylor terms, with the reference frequency set to the 1400 MHz. Following mfs-clean, a spectral index map and its associated error map are produced; these are presented in Figure 11 alongside the Stokes I map at 1400 MHz.

From Figure 11, B2007-569 has a bright central component with the extended tails fading in surface-brightness. The spectral index map indicates that the spectrum of the central component is relatively flat (typically ), with a marked steepening toward the fainter tails. Both of these results are consistent with previously published images of B2007-569 (for example Röttgering et al. 1997). Imaging B2007-569 across the full CABB bandwidth has indicated that the tail components rapidly fade below the detectable flux limit, yielding images morphologically consistent with those of Carretti et al. (2013).

In the upper band, only the point sources were subtracted, in order to investigate their contribution to the flux from the bridge region. In the lower band, with the improved sensitivity to larger angular scales, B2007-569 was subtracted as well as the point sources. The frequency-dependent ATCA CCM was scaled to account for the position of B2007-569 relative to the phase centre of each KAT-7 pointing and subtracted from the uv data.

{subfigure}

[tl]0.45 {subfigure}[tr]0.45 {subfigure}[bl]0.45 {subfigure}[br]0.45

Figure 12: Full-resolution image of A3667 at 1826 MHz, prior to source subtraction. mJy beam.
Figure 13: Full-resolution image of A3667 at 1826 MHz, with point sources (as identified in Carretti et al. 2013) subtracted. mJy beam.
Figure 14: Full-resolution image of A3667 at 1372 MHz, with head-tail radio galaxy B2007-569 and point sources (as identified in Carretti et al. (2013) subtracted. mJy beam.
Figure 15: uv-tapered image A3667 at 1372 MHz, with resolution , to match the resolution of the the higher-frequency observations in Carretti et al. (2013). mJy beam.
Figure 16: Stokes I maps of A3667 observed with KAT-7. Grey contours are positive Stokes I as observed with KAT-7 at [1,2,3,4,6,8,12,16,32], where is stated in the image caption. Black contours are SUMSS data at 843 MHz, starting at 3 mJy beam and scaling by a factor 2. Note that the X-shaped emission extending from the location of B2007-569 has been identified as imaging artefact. Maps at matching frequencies (i.e. Figures 16 and 16; Figures 16 and 16) are set to matching colour scales to facilitate comparison.

The results of subtraction are presented in Figure 16. From the upper band maps – Figures 16 and 16 – the point sources in the bridge region appear to dominate the contribution to the detected flux, as minimal residual emission remains above the noise level. In the lower band maps, subtraction of the ATCA model of B2007-569 results in strong residual emission co-located with the head-tail radio galaxy, shown in Figure 16. We also apply a uv-taper to our 1372 MHz map to match the resolution of the Parkes telescope at 3.3 GHz, yielding a map that exhibits strong morphological resemblance to maps presented in Carretti et al. (2013; Figure 1). This is presented in Figure 16.

There is an extension of emission toward the cluster centre seen in the uv-tapered map (see Figure 16) which could indicate the presence of a radio bridge. However, examination of the wide-field ATCA image at 1400 MHz reveals 3 point sources in the region of this extension, all exceeding significance. These were either not detected in the 2.3–3.0 GHz frequency range imaged in Carretti et al. (2013) or below the 3 mJy threshold of the catalog. It is possible that these sources contribute the flux recovered by KAT-7, as the emission in this region is only detected at the level.

5 Discussion

The residual emission seen in Figure 16 lies spatially co-located with B2007-569. We do not have enough beams across the source to fully determine its size, but it appears extended in nature with an apparent major axis of approximately 500 arcsec; the deconvolved size is 270 arcsec. This would correspond to a physical size of 278 kpc, consistent with the scale size of known mini-haloes (typically  kpc in size, e.g. Kale et al. 2013). One explanation is that the residual emission might be attributed to flux from the tails of B2007-569 that exists on scales unrecovered by the ATCA observations. However, the archival ATCA data were taken in the 1.5D configuration, with baselines in the range 105 m to 4.4 km. At 1400 MHz therefore, the maximum scale size recovered by these data will be of the order of 7 arcmin, which overlaps well with the scale size range recovered by KAT-7. Hence the ATCA observations should be missing little flux from B2007-569 unless there exists some power on very large angular scales. An alternative explanation is that this residual emission is a new candidate radio mini-halo.

Hindson et al. (2014) have shown that A3667 lies coincident with a region of very extended Galactic emission. From this work we cannot rule out the possibility that the recovered emission we tentatively identify as a mini-halo is the result of poorly-sampled Galactic emission. However, from the Galactic synchrotron analysis of La Porta et al. (2008) and following Carretti et al. (2013) we estimate the contamination at the 1372 MHz beam scale of KAT-7 (4.19 arcmin) to be mJy. This contribution is small compared to the flux density recovered from our MHc (which peaks at mJy beam).

The suggestion of a radio halo in A3667 is not a new idea. Johnston-Hollitt (2003 and 2004) report a candidate radio halo with an integrated flux density of mJy at 1.4GHz. Additionally, Carretti et al. (2013) suggest a radio halo of integrated flux mJy at 3.3 GHz. However, both these candidate haloes are coincident with the peak of X-ray emission, whereas this mini-halo candidate (MHc) lies approximately coincident with B2007-569, the brightest cluster galaxy (BCG). The co-location of the MHc and BCG is consistent with previous mini-halo (MH) detections.

Giacintucci et al. (2014) present the most recent catalog of known MH, listing a total of 15 confirmed MH, with 3 additional MHc and a further 3 whose classification as MH is uncertain. Taking radio power from Giacintucci et al. (2014; Table 5, confirmed/candidates only) and bolometric X-ray luminosity from the Archive of Chandra Cluster Entropy Profile Tables (ACCEPT; Cavagnolo et al. 2008) we investigate whether the proposed MHc in A3667 is consistent with trends in the larger MH population.

{subfigure}

[l]0.497 {subfigure}[r]0.497

Figure 17: Scaling trends for clusters with confirmed or candidate MH. left panel shows the trend in MH and BCG radio power scaling; right panel shows the scaling trend between cluster X-ray bolometric luminosity and MH radio power. Unfilled circles are clusters with confirmed MH, unfilled triangles are clusters with candidate MH (Giacintucci et al., 2014); filled symbols correspond to values for A3667. All cluster L values are from ACCEPT (Cavagnolo et al., 2008); cluster MH/BCG radio power are quoted from Giacintucci et al. (2014) with the exception of A3667, where the MHc/BCG radio power were determined in this paper and have here been scaled to match the cosmology of Giacintucci et al. (2014). Errors in measured radio power are included in the plot.

This new MHc possesses an integrated flux density of 67.2 mJy at 1372 MHz, which corresponds to a radio power of P W Hz with our cosmology. Our 1400 MHz ATCA map of B2007-569 has an integrated flux density of mJy, corresponding to a radio power of P W Hz. Despite large scatter, from Figure 17, A3667 appears consistent with the trend that clusters which host more powerful MH tend to host more powerful BCGs (e.g. Govoni et al. 2009; Giacintucci et al. 2014). Additionally, A3667 is approximately consistent with the trend that the more powerful MH tend to be hosted by more X-ray bright clusters. Note that in Figure 17 we have scaled the radio power of the MHc in A3667 to match the cosmology of Giacintucci et al. (2014).

Hindson et al. (2014) subtract B2007-569 from their 226 MHz data by modelling it as a point source. After subtraction, Hindson et al. recover emission at the level Jy beam . In estimating this error we quote the image rms (50 mJy beam) and the approximate calibration uncertainty ( per cent) from Hindson et al. (2014). The flux density recovered from this new MHc peaks at mJy beam, equivalent to a surface brightness of approximately 1.1 Jy arcsec. A spectral index of would yield a 226 MHz flux density consistent with that recovered by Hindson et al., although this is slightly flatter than expected from observations.

A small sample of cool-core clusters host radio relics, but these are generally believed to arise as a result of weak shocks or off-axis mergers (Feretti et al., 2012). On the other hand, while it is a confirmed cool-core cluster (albeit a weak cool-core; Hudson et al. 2010) the merger event in A3667 is a major merger, with the sub-clusters colliding approximately head-on. It has been suggested however, that A3667 is observed in a post-merger state (Finoguenov et al., 2010). Two questions arise: (1) could a cool-core in A3667 survive a major merger event and (2) if a cool-core cannot survive a major merger, has sufficient time elapsed for re-establishment of a cool-core?

Based on a selection of clusters from the Giant Metrewave Radio Telescope (GMRT) Radio Halo Survey (GRHS; see Venturi et al. 2007; Venturi et al. 2008 for details) Rossetti et al. (2011) derive the time-scale () on which non-cool core clusters relax and re-establish cool cores. They derive an upper limit of 1.7 Gyr, and time-scales in the range 1.22.7 Gyr with the inclusion of clusters in the HIFLUGCS subsample of the ACCEPT catalog that were not observed during the GRHS. This is consistent with expectations based on the estimations of radiative lifetime of radio haloes, Gyr (Brunetti et al., 2009). This relaxation time-scale is significantly shorter than the expected cooling lifetime of non-cool-core clusters based on thermodynamic arguments, which is of the order of 10 Gyr (Rossetti & Molendi, 2010). Based on X-ray observations, the structure in the SE cold front in A3667 bears close resemblance to the simulated (Heinz et al. 2003) structure at 3–4 Gyr post core-passage (Briel, Finoguenov & Henry, 2004). While this would be consistent with such a timescale, a newly-formed cool core would likely lead to a detectable cooling flow; however observations have not suggested the presence of a cooling flow in A3667.

Additionally, in clusters which possess cooling flows the cool-core is typically found coincident with the peak of the X-ray emission, whereas the MHc and BCG in A3667 are co-located with the secondary sub-cluster, and significantly offset from the X-ray emission peak. Previous X-ray observations indicate the presence of a cool region of gas co-located with the BCG (Markevitch et al. 1999) which could indicate a relic cool-core, though further observations are required to confirm this. Some simulations have suggested that major merger events occurring in the local Universe have difficulty destroying established cool cores (Burns et al., 2008); however the question of the extent to which merger events disrupt cool cores is still open.

Given the co-location of the BCG (and MHc) in A3667 with the secondary sub-cluster (travelling in a NW direction; Owers, Couch & Nulsen 2009) offset from the peak of the X-ray emission, and lack of observed cooling flow, we suggest that this would likely not be a newly-formed cool-core, but a relic cool-core which has survived the cluster merger. In this scenario, the likely source of electron acceleration would be gas sloshing caused by perturbation of the cool-core by the merger event (e.g. ZuHone et al. 2011).

6 Conclusions

In this work, we have mapped the dual radio relics of the galaxy cluster Abell 3667 in both total intensity and linearly-polarized intensity at 1826 and 1372 MHz with KAT-7. We recover integrated flux densities and polarization fractions that are broadly consistent with previous work at similar frequencies.

Combining our new KAT-7 data with ancillary data from SUMSS we have mapped the spectral profile of both radio relics between 1826 and 843 MHz. We find typical spectral index values of in the North-West (South-East) relic. The North-West relic exhibits a rough transverse gradient, steepening toward the cluster centre and flattening toward the outer edge, indicative of differences in the age of the electrons responsible for the emission.

The typical spectral index of the South-East relic is consistent with the predicted spectral index according to the diffusive shock acceleration mechanism. Our polarization analysis is additionally consistent with this view, as the magnetic field vectors align along the long axis of the South-East relic. The polarization analysis also suggests that the merger axis may be slightly inclined, as the magnetic field vectors for the North-West relic do not quite lie along the long relic axis, indicative of Faraday rotation along the line-of-sight to this relic that does not occur for the South-East relic.

Using the spectral index map in conjunction with our KAT-7 map at 1826 MHz, we have estimated the equipartition magnetic field strength () for the NW relic using the formalism of Beck & Krause (2005). We find typical values of G, consistent with previous estimates from equipartition arguments.

We have subtracted point sources identified in previous work from the North-West relic region to isolate flux density contribution from the diffuse emission. Additionally, we have subtracted the head-tail radio galaxy MRC B2007-569 by using archival ATCA observations to produce a model using multi-frequency synthesis across the wide bandwidth available with the CABB. Localised residual emission remains, spatially co-located with MRC B2007-569, which we have tentatively identified as a MHc. This MHc has an integrated flux density of mJy at 1.37 GHz, corresponding to P W Hz. Subsequent investigation has shown this to be consistent with previously-established trends.

We suggest that the turbulence required for electron reacceleration may be generated by gas sloshing set off by the major merger perturbing (but not destroying) the weak cool-core in A3667 (see Hudson et al. 2010 for the classification of a weak cool-core cluster). This scenario is consistent with some simulations which have shown that established cool-cores in the local Universe are difficult to destroy with major merger events (Burns et al., 2008).

6.1 Acknowledgements

This work was carried out as part of commissioning efforts of KAT-7. All KAT-7 data and images are provided courtesy of SKA SA - http://www.ska.ac.za. We would like to thank the Square Kilometre Array South Africa (SKA SA) collaboration, and recognise the efforts of the KAT-7 observers and engineers for making this work possible. We extend our thanks to our referee M. Johnston-Hollitt for detailed comments and suggestions that assisted in the interpretation of our results and improved our scientific output. CJR gratefully acknowledges both support from the UK Science & Technology Facilities Council (STFC) and the support from SKA SA that enabled this work to be conducted. AMS gratefully acknowledges support from the European Research Council under grant ERC-2012-StG-307215 LODESTONE.

CJR wishes to thank T. Mauch and L. Richter for helpful discussions during the data reduction process; and M. Johnston-Hollitt for informative discussions during results analysis. This work made use of the cosmology calculator developed by Wright (2006) and the Archive of Chandra Cluster Entropy Profile Tables (ACCEPT; Cavagnolo et al. 2008). The colour schemes for Figures 5, 6 and 11 were made using the cubehelix colour scheme generator developed by D. Green (Green, 2011) to generate a custom colour scheme specifically for this work.

Footnotes

  1. pagerange: LABEL:firstpageLABEL:lastpage
  2. pubyear: 2014
  3. Counting only relics classified as ‘elongated’ per the nomenclature of Feretti et al. (2012); these are analogous to the ‘radio gischt’ of Kempner et al. (2004). We exclude ‘roundish’ relics as these are more akin to the radio ‘phoenix’ classification of source.
  4. http://casaguides.nrao.edu/index.php?title=Flagdata
  5. http://casa.nrao.edu
  6. We note that this ambiguity can be avoided by use of an injected XY phase calibration signal, such as that utilised by the ATCA instrument.
  7. http://casaguides.nrao.edu/index.php?title=Imfit
  8. These sources are SUMSS J201330-570552; a suspected FRII source (Johnston-Hollitt, 2003) and the head-tail radio source PMN J2014-5701.

References

  1. Akamatsu H., de Plaa J., Kaastra J., Ishisaki Y., Ohashi T., Kawaharada M., Nakazawa K., 2012a, PASJ, 64, 49
  2. Akamatsu H., Takizawa M., Nakazawa K., Fukazawa Y., Ishisaki Y., Ohashi T., 2012b, PASJ, 64, 67
  3. Armstrong R. P. et al., 2013, MNRAS, 433, 1951
  4. Baars J. W. M., Genzel R., Pauliny-Toth I. I. K., Witzel A., 1977, A&A, 61, 99
  5. Beck R., Krause M., 2005, Astronomische Nachrichten, 326, 414
  6. Bock D. C.-J., Large M. I., Sadler E. M., 1999, AJ, 117, 1578
  7. Brentjens M. A., de Bruyn A. G., 2005, A&A, 441, 1217
  8. Briel U. G., Finoguenov A., Henry J. P., 2004, A&A, 426, 1
  9. Brunetti G., Cassano R., Dolag K., Setti G., 2009, A&A, 507, 661
  10. Burns J. O., Hallman E. J., Gantner B., Motl P. M., Norman M. L., 2008, ApJ, 675, 1125
  11. Carretti E. et al., 2013, MNRAS, 430, 1414
  12. Cassano R., Gitti M., Brunetti G., 2008, A&A, 486, L31
  13. Cavagnolo K. W., Donahue M., Voit G. M., Sun M., 2008, ApJ, 683, L107
  14. de Gasperin F., van Weeren R. J., Bruggen M., Vazza F., Bonafede A., Intema H. T., 2014, ArXiv e-prints
  15. de Zotti G., Massardi M., Negrello M., Wall J., 2010, A&A Rev., 18, 1
  16. Dennison B., 1980, ApJ, 239, L93
  17. Feretti L., 2005, Advances in Space Research, 36, 729
  18. Feretti L., Giovannini G., Govoni F., Murgia M., 2012, A&A Rev., 20, 54
  19. Ferrari C., Govoni F., Schindler S., Bykov A. M., Rephaeli Y., 2008, Space Sci.Rev., 134, 93
  20. Finoguenov A., Sarazin C. L., Nakazawa K., Wik D. R., Clarke T. E., 2010, ApJ, 715, 1143
  21. Giacintucci S., Markevitch M., Venturi T., Clarke T. E., Cassano R., Mazzotta P., 2014, ApJ, 781, 9
  22. Gitti M., Brunetti G., Feretti L., Setti G., 2004, A&A, 417, 1
  23. Gitti M., Brunetti G., Setti G., 2002, A&A, 386, 456
  24. Goss W. M., Ekers R. D., Skellern D. J., Smith R. M., 1982, MNRAS, 198, 259
  25. Govoni F., Enßlin T. A., Feretti L., Giovannini G., 2001, A&A, 369, 441
  26. Govoni F., Murgia M., Markevitch M., Feretti L., Giovannini G., Taylor G. B., Carretti E., 2009, A&A, 499, 371
  27. Green D. A., 2007, Bulletin of the Astronomical Society of India, 35, 77
  28. Green D. A., 2011, Bulletin of the Astronomical Society of India, 39, 289
  29. Gutierrez K., Krawczynski H., 2005, ApJ, 619, 161
  30. Heald G., 2009, in IAU Symposium, Vol. 259, IAU Symposium, Strassmeier K. G., Kosovichev A. G., Beckman J. E., eds., pp. 591–602
  31. Heinz S., Churazov E., Forman W., Jones C., Briel U. G., 2003, MNRAS, 346, 13
  32. Hindson L. et al., 2014, MNRAS, 445, 330
  33. Hudson D. S., Mittal R., Reiprich T. H., Nulsen P. E. J., Andernach H., Sarazin C. L., 2010, A&A, 513, A37
  34. Hunstead R. W., 1991, Australian Journal of Physics, 44, 743
  35. Hurley-Walker N. et al., 2014, ArXiv e-prints
  36. Johnston-Hollitt M., 2003, PhD thesis, University of Adelaide
  37. Johnston-Hollitt M., 2004, in The Riddle of Cooling Flows in Galaxies and Clusters of galaxies, Reiprich T., Kempner J., Soker N., eds., p. 51
  38. Johnston-Hollitt M., Hunstead R. W., Corbett E., 2008, A&A, 479, 1
  39. Kale R., Venturi T., Giacintucci S., Dallacasa D., Cassano R., Brunetti G., Macario G., Athreya R., 2013, A&A, 557, A99
  40. Kempner J. C., Blanton E. L., Clarke T. E., Enßlin T. A., Johnston-Hollitt M., Rudnick L., 2004, in The Riddle of Cooling Flows in Galaxies and Clusters of galaxies, Reiprich T., Kempner J., Soker N., eds., p. 335
  41. La Porta L., Burigana C., Reich W., Reich P., 2008, A&A, 479, 641
  42. Large M. I., Mills B. Y., Little A. G., Crawford D. F., Sutton J. M., 1981, MNRAS, 194, 693
  43. Macario G., Markevitch M., Giacintucci S., Brunetti G., Venturi T., Murray S. S., 2011, ApJ, 728, 82
  44. Markevitch M., Sarazin C. L., Vikhlinin A., 1999, ApJ, 521, 526
  45. Mauch T., Murphy T., Buttery H. J., Curran J., Hunstead R. W., Piestrzynski B., Robertson J. G., Sadler E. M., 2003, MNRAS, 342, 1117
  46. Nakazawa K. et al., 2009, PASJ, 61, 339
  47. Oppermann N. et al., 2012, A&A, 542, A93
  48. Owers M. S., Couch W. J., Nulsen P. E. J., 2009, ApJ, 693, 901
  49. Owers M. S., Nulsen P. E. J., Couch W. J., Markevitch M., 2009, ApJ, 704, 1349
  50. Perley R. A., Butler B. J., 2013, ApJS, 206, 16
  51. Pfrommer C., Enßlin T. A., 2004, MNRAS, 352, 76
  52. Rau U., Cornwell T. J., 2011, A&A, 532, A71
  53. Rayner D. P., Norris R. P., Sault R. J., 2000, MNRAS, 319, 484
  54. Reiprich T. H., Böhringer H., 2002, ApJ, 567, 716
  55. Reynolds J. E., 1994, A Revised Flux Scale for the AT Compact Array, Tech. Rep. 39.3/040, Australia Telescope National Facility – . 1997, ATCA Calibrator Source Catalogue, ftp://ftp.atnf.csiro.au/pub/atnfdo cs/guides/at.cat
  56. Rossetti M., Eckert D., Cavalleri B. M., Molendi S., Gastaldello F., Ghizzardi S., 2011, A&A, 532, A123
  57. Rossetti M., Molendi S., 2010, A&A, 510, A83
  58. Röttgering H. J. A., Wieringa M. H., Hunstead R. W., Ekers R. D., 1997, MNRAS, 290, 577
  59. Sault R. J., Teuben P. J., Wright M. C. H., 1995, in Astronomical Society of the Pacific Conference Series, Vol. 77, Astronomical Data Analysis Software and Systems IV, Shaw R. A., Payne H. E., Hayes J. J. E., eds., p. 433
  60. Schenck D. E., Datta A., Burns J. O., Skillman S., 2014, AJ, 148, 23
  61. Slee O. B., Roy A. L., Murgia M., Andernach H., Ehle M., 2001, AJ, 122, 1172
  62. Sodre, Jr. L., Capelato H. V., Steiner J. E., Proust D., Mazure A., 1992, MNRAS, 259, 233
  63. Tingay S. J. et al., 2013, PASA, 30, 7
  64. van Weeren R. J., Röttgering H. J. A., Brüggen M., Hoeft M., 2010, Science, 330, 347
  65. van Weeren R. J., Röttgering H. J. A., Intema H. T., Rudnick L., Brüggen M., Hoeft M., Oonk J. B. R., 2012, A&A, 546, A124
  66. Venturi T., Giacintucci S., Brunetti G., Cassano R., Bardelli S., Dallacasa D., Setti G., 2007, A&A, 463, 937
  67. Venturi T., Giacintucci S., Dallacasa D., Cassano R., Brunetti G., Bardelli S., Setti G., 2008, A&A, 484, 327
  68. Vikhlinin A., Markevitch M., Murray S. S., 2001, ApJ, 549, L47
  69. Williams C. L. et al., 2012, ApJ, 755, 47
  70. Wright A. E., Griffith M. R., Burke B. F., Ekers R. D., 1994, ApJS, 91, 111
  71. Wright E. L., 2006, PASP, 118, 1711
  72. Wyllie D. V., 1969, MNRAS, 142, 229
  73. ZuHone J., Markevitch M., Brunetti G., 2011, MmSAI, 82, 632
  74. ZuHone J. A., Markevitch M., Brunetti G., Giacintucci S., 2013, ApJ, 762, 78
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
130040
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description