A&A 486, 347-358 (2008)
DOI: 10.1051/0004-6361:200809459
S. Giacintucci1,2 - T. Venturi1 - G. Macario1,3 - D. Dallacasa1,3 - G. Brunetti1 - M. Markevitch2 - R. Cassano1,3 - S. Bardelli4 - R. Athreya5
1 -
INAF - Istituto di Radioastronomia, Via Gobetti 101, 40129 Bologna, Italy
2 -
Harvard-Smithsonian Centre for Astrophysics,
60 Garden Street, Cambridge, MA 02138, USA
3 -
Dipartimento di Astronomia, Università di Bologna,
via Ranzani 1, 40126 Bologna, Italy
4 -
INAF - Osservatorio Astronomico di Bologna, via Ranzani 1,
40126 Bologna, Italy
5 -
Tata Institute of Fundamental Research,
National Centre for Radio Astrophysics,
Ganeshkhind, Pune 411007, India
Received 25 January 2008 / Accepted 22 March 2008
Abstract
Aims. We present new high sensitivity observations of the radio relic in A 521 carried out using the Giant Metrewave Radio Telescope at 327 MHz and with the Very Large Array at 4.9 and 8.5 GHz.
Methods. We imaged the relic at these frequencies and carried out a detailed spectral analysis, based on the integrated radio spectrum between 235 MHz and 4.9 GHz, and on the spectral index image in the frequency range 327-610 MHz. In our present analysis we use our new GMRT observations in addition to proprietary and archival data. We search for a possible shock front co-located with the relic on a short archival Chandra X-ray observation of the cluster.
Results. The integrated spectrum of the relic is consistent with a single power law; the spectral index image shows a clear trend of steepening going from the outer portion of the relic toward the cluster centre. We discuss the origin of the source in the light of theoretical models for the formation of cluster radio relics. Our results on the spectral properties of the relic are consistent with acceleration of relativistic electrons by a shock in the intracluster medium. This scenario is supported by our detection of an X-ray surface brightness edge coincident with the outer border of the radio relic. This edge is probably a shock front.
Key words: radio continuum: galaxies - galaxies: clusters: general - galaxies: clusters: individual: A 521
The radio emission from clusters of galaxies comes into two flavours. In addition to the individual radio sources associated with cluster galaxies, a fraction of massive and X-ray luminous clusters with clear signature of ongoing mergers, hosts large scale diffuse radio sources, known as radio halos, if located at the cluster centre, and relics, if located in peripheral regions. The synchrotron emission in these diffuse sources arises directly within the intracluster medium (ICM), and probes the existence of non-thermal components spread over the cluster scale (e.g., the review paper by Feretti 2003).
Particle acceleration by means of turbulence injected into the cluster volume during mergers represents a promising possibility to understand the origin of radio halos (see the review papers by Brunetti 2003, 2004; Sarazin 2004; Petrosian & Bykov 2008). Statistical analyses (e.g., Kuo et al. 2004; Cassano & Brunetti 2005; Cassano et al. 2006, 2008) and deep radio observations of samples of galaxy clusters (Brunetti et al. 2007; Venturi et al. 2007, 2008) have provided further support to this scenario. On the other hand, our understanding of the origin of relics is still limited. Only a handful of radio relics have been observed in some detail (e.g., A 2256, Clarke & Enßlin 2006; A 3667, Röttgering et al. 1997; A 2744, Orrú et al. 2007). Integrated radio spectra over a wide range of frequencies are available for a few relics only (e.g., 1257+275 in the Coma cluster; Andernach et al. 1984; Thierbach et al. 2003, and references therein). All models proposed so far for the relic formation invoke the presence of a shock within the X-ray gas (Enßlin et al. 1998; Roettiger et al. 1999; Enßlin & Gopal-Krishna 2001; Hoeft & Brüggen 2007).
We focus our attention on the radio relic in A 521 (Ferrari et al.
2006; Giacintucci et al. 2006, hereinafter GVB06), an X-ray luminous and
massive galaxy cluster (
erg s-1; virial mass
) at
redshift z=0.247. Multiple merging episodes are known to occur in
this disturbed cluster whose properties are indicative of an object
in a complex dynamical state, which is still accreting a number of smaller
mass concentrations (e.g., Maurogordato et al. 2000; Ferrari
et al. 2003, 2006; see also Fig. 1 in GVB06 for a sketch of
the multiple optical and X-ray substructures in the cluster).
Table 1: Summary of the radio observations.
A radio study of A 521 based on 610 MHz Giant Metrewave Radio Telescope
(GMRT) observations (GVB06) shows that the relic is a diffuse elongated
structure located in the southeastern periphery of A 521, on the edge
of a dynamically active region where galaxy groups are infalling into the main
cluster (Maurogordato et al. 2000; Ferrari et al.
2003). The relic is at the boundary of the X-ray emission from the intracluster gas, at a projected distance of 930 kpc from the cluster
X-ray centre, and is apparently connected to the most powerful radio
galaxy in A 521 (J0454-1016a) by a faint bridge of radio emission. Even though
projection effects should be taken into account, this situation is similar to what is observed in the Coma cluster, where a bridge of radio emission connects the tails
of the radio galaxy NGC 4789 to the prototype relic source 1253+275 (Giovannini et al. 1991).
We present a study of the origin of the relic based on new GMRT 327 MHz high sensitivity observations of A 521 and new high frequency and high resolution images of the relic region obtained from Very Large Array (VLA) observations at 4.9 and 8.5 GHz. The amount of radio information available allows us to study in detail the spectral properties of the source (integrated and point-to-point) over a relatively wide range of frequencies, and compare them to the model predictions of its origin. The paper is organised as follows: Sect. 2 describes the radio observations and data reduction; the new GMRT 327 MHz images of the A 521 field and relic region are presented in Sect. 3; in Sect. 4, we report on the analysis of the VLA 4.9 and 8.5 GHz images; Sect. 5 deals with the study of the spectral index image; the source integrated radio spectrum is analysed in Sect. 6; the proposed scenarios for the relic origin are discussed in Sect. 7; finally a summary of our results is given in Sect. 8.
We adopt the CDM cosmology with H0=70 km s-1 Mpc-1,
and
.
At the redshift of A 521, this cosmology
leads to a linear scale of 1'' = 3.87 kpc. The spectral index
is defined according to
.
![]() |
Figure 1:
GMRT 327 MHz contours of the
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We carried out high sensitivity observations of A 521 using the
GMRT at 327 MHz. To image the sky region close to the relic at
higher frequency and resolution, and resolve the inner structure
of the radio galaxy J0454-1016a, we performed VLA observations
at 4.9 GHz in the hybrid BnA and CnB configurations, and at 8.5 GHz
in the BnA configuration. In Table 1, we summarise the
details of all the radio observations presented in this paper. Columns
in the table provide the following information: telescope/array; J2000 coordinates of the pointing centre; observing day; frequency; total
bandwidth; total time on source; half power bandwidth (HPBW) and rms
level (1)
in the full resolution image.
A 521 was observed using the GMRT at 327 MHz in November 2006 for a total integration time of 5.5 h (Table 1). These observations were performed using both the upper and lower side band (USB and LSB, respectively) for a total observing bandwidth of 32 MHz. The data were collected in spectral-line mode with 128 channels/band, and a spectral resolution of 125 kHz/channel.
The USB and LSB datasets were calibrated and analysed individually using the NRAO Astronomical Image Processing System (AIPS) package. The bandpass calibration was performed using the flux density calibrator. A RFI-free channel was chosen to normalise the bandpass for each antenna. The calibration solutions were applied to the data by running the AIPS task FLGIT, which subtracts a continuum from the channels in the u-v plane, determined on the basis of the bandpass shape and using a specified set of channels. The data whose residuals exceed a chosen threshold are then flagged. Despite this flagging procedure, both the USB and LSB datasets were still affected by strong residual radio frequency interferences (RFI). Hence, an accurate editing of the visibility data was carried out to identify and remove data affected by RFI.
To find a compromise between the size of the dataset and the need to
minimize bandwidth smearing effects within the primary beam, the
central channels were averaged to 6 channels of 2 MHz each after
bandpass calibration. Given the large field of view of the GMRT,
in each step of the data reduction we implemented the wide-field
imaging technique to minimise the errors due to the non-planar nature
of the sky. We used 25 facets covering a total field of view of
square degrees. After a number of phase
self-calibration cycles, the final USB and LSB datasets were
averaged further from 6 channels to 1 single channel
.
Due to residual phase errors in the LSB dataset, the USB-LSB data
combination led to images with a quality worse than those obtained
from the USB alone. For this reason, only the USB dataset was used for
the analysis presented in this present paper. A very high sensitivity
(1)
was achieved in our final full resolution image
(Table 1): from
90
Jy b-1 in the region
of the relic to
100
Jy b-1 in the outer parts
of the A 521 field, where the quality of the image is limited by the
presence of strong radio sources. The residual amplitude errors are
approximately
5%.
The 4.9 GHz observations of the relic region were carried
out using the VLA in the hybrid BnA and CnB configurations.
The 8.5 GHz observations were performed in the BnA configuration
(see Table 1 for details). A bandwidth of 50 MHz was
used for each of the two IF channels at each frequency. All
observations included full polarisation information.
The data were calibrated and reduced using the pilot semi-automatic
pipeline for the VLA data processing, a facility developed
at NRAO and implemented in the AIPS package. The pipeline provided
high quality calibrated datasets at each frequency. These datasets
were further phase self-calibrated in AIPS, to correct for residual
phase variations, and used to produce the final images. The rms noise
level achieved in the final images are reported in Table 1.
The average residual amplitude errors in the data are approximately 5
both at 4.9 and 8.5 GHz.
In Fig. 1, we present the GMRT full resolution
image at 327 MHz that covers the region within the cluster virial
radius (
Mpc; GVB06), which is delimited by the
solid circle. The figure shows the same
field presented at 610 MHz in Fig. 2 of GVB06. The
cross marks the centre of the cluster X-ray emission as detected by
ROSAT HRI (Arnaud et al. 2000). The dashed circle, which
has a radius corresponding to
1 Mpc, indicates the area
covered by the analysis of the optical substructures in Ferrari et al.
(2003). The radio contours are plotted starting from
0.5 mJy b-1, which corresponds to the 5
detection
significance level in the region with the highest noise.
Three extended radio sources are visible in the image: the radio relic, in the southeastern outskirts of the cluster, and the two radio galaxies J0453-0957 and J0454-1006, located North of A 521, and analysed at 610 MHz in GVB06.
In addition to discrete point sources, positive residuals of radio
emission are detected within the dashed circle in Fig. 1,
which suggest that diffuse emission is present at the cluster centre.
The investigation of this point is beyond the purpose of the present
paper, and will be addressed in a forthcoming paper (Brunetti et al., to be submitted).
![]() |
Figure 2:
Left panel: full resolution GMRT contours at 327 MHz
of the relic, overlaid on the red POSS-2 optical image. The resolution
is
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
As observed at higher frequencies, the 327 MHz radio emission within
the inner 1 Mpc radius is clearly dominated by the relic.
Figure 2 zooms into the 327 MHz image of the relic.
In the left panel, we show the full resolution contours overlaid on
the optical POSS-2 frame (grey scale). Labels A, B and C indicate
the position of the radio galaxies embedded in the diffuse relic
emission, and optically identified in GVB06. In the right panel,
we show an image at the resolution of
with grey scale and contours overlaid to highlight
the distribution of the radio surface brightness across the source.
The relic exhibits a highly elongated and arc-shaped structure with an
angular size of
,
which corresponds to a linear size
of
1 Mpc. The overall morphology and total extent in
Fig. 2 are in good agreement with the images at
610 MHz (GVB06; also reported in Fig. 5) and at 1.4 GHz (Ferrari
et al. 2006) of similar resolution. The relic emission along the minor
axis appears, on average, to be slightly wider in the 327 MHz image
(
,
i.e.,
230 kpc) than at higher frequencies
(
200 kpc at 610 MHz and
160 kpc at 1.4 GHz).
![]() |
Figure 3:
Panel a) GMRT 610 MHz full resolution contours of the
relic on the POSS-2 image. The HPBW is
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Source A in Fig. 2 is the cluster radio galaxy J0454-1016a. A faint bridge of emission is clearly detected between the northern part of the relic and this radio galaxy, which is located at a projected distance of approximately 350 kpc from the relic. This emission was also observed at 610 MHz, and suggested the presence of a physical link between the relic and J0454-1016a (GVB06). We studied the possible connection between these sources using VLA 4.9 and 8.5 GHz observations (Table 1), which resolve the inner structure of J0454-1016a and enable us to search for connections with the nearby relic in the form, for example, of bent jets and/or extended emission in that direction.
J0454-1016a is the most powerful radio galaxy in A 521 (GVB06;
Table 2), and is identified with the galaxy
(v=74 282 km s-1, I=17.00) in the optical catalogue by
Ferrari et al. (2003). Figure 3
shows the region of the relic and J0454-1016a (labelled as A)
in increasing frequency (from 610 MHz to 8.5 GHz) and resolution
order (
to
)
going from panels a) to d). Panel a) shows
the full resolution image at 610 MHz. The VLA-CnB 4.9 GHz image
of J0454-1016a is presented in panel b); the BnA full resolution
images at 4.9 and 8.5 GHz are shown in panels c) and d),
respectively. J0454-1016a appears extended in all images with a largest
linear size of
30 kpc. Its radio structure is consistent with
a head-tail morphology. The compact component, detected at 4.9 and
8.4 GHz (panels c and d), is coincident with the nucleus of
the host galaxy.
Table 2: Properties of the radio galaxy J0454-1016a.
The extended emission is entirely located southwest of the compact component, and there is no evidence of any emission in the direction of the radio bridge detected at 327 MHz (Fig. 2) and 610 MHz (panel a). This result might rule out a physical connection between the diffuse emission from the relic and this radio galaxy, whose tail extends (at least in projection) in almost the opposite direction.
![]() |
Figure 4:
Integrated radio spectrum of J0454-1016a between
327 MHz and 8.5 GHz. The solid line is the best fit of the
CI model to the data. The value of
![]() ![]() |
Open with DEXTER |
The flux densities available for J0454-1016a are presented in Table 2, where we also report the 1.4 GHz radio power (from VLA archival data of program AF 0390, data which were reanalysed by Giacintucci 2007), and the spectral index in the 327 MHz-8.5 GHz interval. The flux density measurements on both images at 4.9 GHz (panels b and c in Fig. 3) are consistent within the errors. Figure 4 shows the integrated radio spectrum of J0454-1016a between 327 MHz and 8.4 GHz, derived using the values in Table 2.
We fitted the spectrum with a a continuous injection model
(CI; Kardashev 1962), using the Synage++ package (Murgia 2001).
The best fit is shown as a solid line in Fig. 4,
and provides an injection spectral index
.
Even though there is an indication of a spectral
steepening above 4.9 GHz, the spectrum is consistent with a single power
law of slope
.
The spectral shape in Fig. 4 is similar to that observed in other active and low luminosity radio galaxies (e.g., Parma et al. 2002).
The most interesting result of the VLA 4.9 GHz observations is the
detection of the radio relic, shown in Fig. 5.
This is the second detection of a radio relic at a frequency as high
as 4.9 GHz, after the relic source 1253+275 in the Coma cluster
(Andernach et al. 1984; Thierbach et al. 2003). The
image was obtained from the VLA-CnB data (Table 1),
tapered to a resolution of
.
In the left panel, the relic is shown as contours, while in the right panel
the 4.9 GHz emission is reported in grey scale with the GMRT 610 MHz
contours overlaid (GVB06). The faintest features of the relic emission
are significant at the 3
level (1
Jy b-1),
and the peaks at the level of 12
.
We point out that, given the
short observing time (Table 1) and lack of short baselines (minimum
baseline
1 k
), the u-v coverage of the 4.9 GHz
observations is inadequate to properly image the entire diffuse emission
associated with the relic. The details of the structure of the source in
Fig. 5 may not be therefore fully reliable, and
deeper 4.9 GHz observations using a more appropriate array and u-v coverage are required to determine more accurately the relic
morphology at this frequency.
Spectral index imaging of cluster radio relics is a powerful tool to attempt to understand the origin, evolution, and connection of the relics with the merging activity of the host cluster. In particular, the spectral index images provide important information on the energy spectrum of the radio emitting electrons and magnetic field distribution in these sources (e.g., Clarke & Enßlin 2006).
We made an image of the spectral index distribution
in the A 521 relic, for the frequency range 327-610 MHz,
by comparing the GMRT image at 327 MHz, shown here (Fig. 2;
right panel), with the GMRT 610 MHz image obtained from the
observations presented in GVB06. The images were produced
using the same cell size, u-v range, and restoring beam.
For details, we refer to Table 3, where we provide
the u-v range, beam, and noise level (1)
of the two images. The images were aligned, pixels with
brightness below the 3
level were blanked, and the
image combination was carried out to create the spectral
index image using the Synage++ package. The resulting spectral
index image is shown in Fig. 6 (colour), with
the 610 MHz contours overlaid.
![]() |
Figure 5:
VLA 4.9 GHz image of the radio relic as contours
( left panel), and grey scale on the GMRT 610 MHz
image ( right panel). The resolution of the 4.9 GHz image
is
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Table 3: Details of the 327-610 MHz spectral index image.
The distribution of the spectral index in Fig. 6 shows
different features along the two axes of the relic. Along
the major axis, a number of irregularities provide a
rather patchy appearance of the spectral index image.
Such patchiness might arise from irregularities in
the u-v coverage that occur at different spacings
for the two frequencies. Along the minor axis, a steepening of
the spectral index from the eastern edge of the relic toward
the western border is visible. This trend is real and not driven
by a misalignment between the images at 327 and 610 MHz. We checked
carefully that the images were aligned correctly using point sources
in the relic region. The spectral index distribution appears
uniform within each of the point sources and consistent
with the value of
derived from their total flux
density at the two frequencies. This is clear
from Fig. 6, where, for example, the
radio galaxy J0454-1016a (labelled as A) has an average
spectral index
,
which agrees with the
source integrated spectrum shown in Fig. 4
(see also Table 2).
To check the significance of the spectral steepening visible in
Fig. 6, we determined the average spectral index in
3 independent strips of dimensions
in the northern region of the relic and 3 independent
strips of dimensions
in the southern part. These two regions are labelled respectively N and S
in the left panel of Fig. 7. The central portion of
the relic was excluded from the analysis because of the presence of the
radio galaxies B and C (Fig. 7), which may affect the
spectral trend. The strips were set parallel to the edge of the relic,
i.e., at a position angle of
and
in
the N and S regions, respectively. For each strip, we integrated the
flux density in the 327 MHz and 610 MHz images individually,
and calculated the corresponding spectral index values.
The average 327-610 MHz spectral index in each strip is
shown in the right panel of Fig. 7.
The spectral index trend shows a gradual steepening going
westwards both in the N and S regions. In particular,
ranges from
to
in the southern part (empty triangles), and
from
to
across
the northern region (filled triangles).
A steepening from the outer to the inner edge of the
relic with available spectral index images has been
observed in a few other cases, for example
in A 3667 (Röttgering et al. 1997) and A 2744
(Orrú et al. 2007). A more detailed study of the spectral index
gradient, similar to the analysis presented in this paper,
was carried out for the relic in A 2256 by
Clarke & Enßlin (2006), who found a significant
steepening of
between 1369 MHz and 1703 MHz from
the external edge of the relic toward the cluster core, i.e.,
a gradient similar to that found for the A 521 relic.
The wealth of multifrequency radio observations available for the relic in A 521 enables its integrated spectrum to be determined over almost two orders of magnitude in frequency.
In Table 4, we report all available flux
densities of the relic, along with the resolution of the images
used for the measurement.
To obtain a consistent measurement of the total flux
densities, we integrated over the same region in all
images, and subtracted the flux density of the embedded point-sources,
as measured in the corresponding full resolution images.
The source subtraction could not be applied to the 74 MHz data,
since the low resolution (80
)
of the Very Low-frequency Sky Survey (VLSS
)
image did not allow such an operation. Furthermore, the extended
structure visible in the VLSS image is not fully consistent with the relic
morphology at higher frequencies; this is probably due to
the different angular resolution and much lower sensitivity of the
VLSS image (which has an average rms noise of approximately
1
Jy b-1). For these reasons, the flux density
measurement at 74 MHz is very uncertain.
![]() |
Figure 6:
Spectral image of the relic in A 521 in the 327-610 MHz
range, with the GMRT 610 MHz contours at levels
0.12, 0.24, 0.48, 0.96, 2, 4 mJy b-1 overlaid. In both images
(spectral index and total intensity) the resolution is
![]() ![]() |
Open with DEXTER |
The integrated spectrum of the relic is shown in
Fig. 8. The spectrum does not
show any steepening up to 4.9 GHz and is well fitted by
a single power law with slope
(solid
line) between 235 MHz and 4.9 GHz. Given its large uncertainty,
the 74 MHz data point was not included in the fit. We note that
the 4.9 GHz flux density should be considered a lower limit, given
the array, u-v coverage, and resolution which led to its detection
(see Sect. 4.2).
New deep observations carried out with the VLA at 74 and 327 MHz (October 2007 and February 2008) will allow us to better constrain the low frequency end of the radio spectrum.
A preliminary discussion of the formation of the radio relic in A 521 was carried out in GVB06, where a number of possible theoretical frameworks were taken into account, all related to the assessed ongoing merging activity in this cluster. Two possible scenarios were considered, both invoking a tight connection with the presence of a merger-driven shock front at the location of the relic. Such a shock may accelerate electrons to ultra-relativistic energies (Enßlin et al. 1998; Roettiger et al. 1999; Hoeft & Brüggen 2007), or it may revive fossil radio plasma through adiabatic compression of the magnetic field (Enßlin & Gopal-Krishna 2001).
A third alternative scenario was also proposed, based on the hypothesis of a physical link between the diffuse emission of the relic and the nearby radio galaxy J0454-1016a, as suggested by the faint radio bridge of emission between the two sources observed at 610 MHz (e.g., Fig. 5, right panel). The high frequency images of J0454-1016a (Fig. 3) and the spectral analysis of the radio relic presented in this paper (Sects. 5 and 6) allow us to carry out a more detailed comparison between the expectations of the models and the observed properties of the relic.
![]() |
Figure 7:
Left - Grids used to derive the spectral
index trend overlaid on the 327-610 MHz spectral index image
(grey scale; same as Fig. 6). Right -
Trend of the spectral index
![]() |
Open with DEXTER |
GVB06 proposed that the relic might be the result of
the ram pressure stripping of the radio lobes of
J0454-1016a by either i) group merging activity
in the southern cluster region; or ii) the infall
of the radio galaxy itself into the cluster. On the basis
of pressure arguments, and given the projected distance
of J0454-1016a from the relic, they estimated that
the infall velocity of the merging group, or of the
galaxy itself, should be 3000 km s-1
(leading to a shock with Mach number
)
to allow the electrons in the radio lobes to still emit
in the radio band.
The high resolution and high frequency observations presented in Sect. 4.1 do not appear to support this scenario, since there is no obvious morphological link between the relic and J0454-1016a. The radio galaxy has a head-tail morphology, whose tails extend (in projection) in a direction that is almost opposite to that towards the relic region and the faint bridge of emission connecting the galaxy and the relic (Fig. 3). Even though the origin of the faint bridge remains unclear, it appears to be unrelated to the current AGN activity of J0454-1016a.
Table 4: Flux density values of the relic.
![]() |
Figure 8: Radio spectrum of the radio relic between 74 MHz and 4.9 GHz. The solid circle is the relic flux density measured on the new GMRT image at 327 MHz. The empty triangles represent the 610 MHz data point from GVB06, and the 1.4 GHz value from VLA archival data. The solid triangle is the source flux density measured on the new VLA image at 4.9 GHz. The empty circle is the flux density from a preliminary GMRT image at 235 MHz (see Appendix). The empty square is the 74 MHz flux density from the VLSS image. The solid line represents the linear fit to the data by using the points between 235 MHz and 4.9 GHz. |
Open with DEXTER |
The shock scenarios for the origin of radio relics provide a number of expectations, which can be tested by means of comparison with the observed properties of the source.
The results of the spectral analysis presented in this
paper provide insightful information to discriminate
between the models described above. In the case
of the relic in A 521, a significant steepening of
the spectral index was found going from the eastern
edge of the source toward the western border (Fig. 7).
This behaviour is in line with the expectations of all
the scenarios. It implies that the shock is expected
to be moving outwards with respect to the cluster
centre, and its current location should be approximately
coincident with the eastern edge of the relic.
In Sect. 6, we showed that
the integrated spectrum of the relic is well reproduced
by a single power law with steep spectral index
(
), and with no evidence of high frequency
steepening up to 4.9 GHz (Fig. 8).
Such a spectral shape is expected only in the framework of
the shock acceleration and re-acceleration
scenario.
![]() |
Figure 9:
Left panel - Smoothed Chandra image of A 521
in the 0.5-4 keV energy band. Right panel - Same image as left
panel with the GMRT 610 MHz contours of the A 521 field overlaid.
The resolution of the radio image is 13.1
![]() ![]() |
Open with DEXTER |
The adiabatic compression scenario should
produce curved spectra with a high frequency cut-off in
the case of moderate or weak shocks, such as those expected
to be developed during cluster mergers (e.g.,
Gabici & Blasi 2003; Ryu et al. 2003; Pfrommer et al. 2006,
and references therein). In principle, a cut-off just above
5 GHz cannot be ruled out. The compression enhances the magnetic
field and the particle energy density, moving the break
frequency in the synchrotron spectrum to higher
frequencies. Assuming thermal and relativistic particles
mixed into the fossil plasma, the ratio of the post-compression
and initial break frequencies (
and
,
respectively) is (e.g., Markevitch et al. 2005):
This value is not sufficiently high to explain the lack of
a cut-off in the observed spectrum of the relic.
This factor can be strongly increased if the shock
compression acts on a ghost of purely relativistic plasma (i.e.,
not mixed with the thermal ICM). In this case (Enßlin & Gopal-Krishna 2001):
An additional point is that the injection spectrum of the fossil
radio plasma in this scenario should be roughly equal to the
observed one (
), which is much steeper
than the typical spectra of radio galaxies.
To summarise, our spectral analysis suggests that the
origin of the relic in A 521 is consistent with the
shock acceleration scenario. According to Eq. (1),
we can estimate the Mach number of the shock responsible
for the electron acceleration. The spectral fit provided a
total spectral index
(Sect. 6),
which corresponds to
,
and thus
.
Such a Mach number is
in reasonable agreement with the values expected
for the cluster merger shocks, and indeed
observed in merging clusters (e.g., Markevitch & Vikhlinin
2007). As argued in GVB06, the relic is located in a
peripheral region of A 521, which is expected
to be dynamically active (Maurogordato et al.
2000; Ferrari et al. 2003, 2006; see also
Fig. 2 in GVB06). Thus the presence of a shock
front at the relic location is likely.
![]() |
Figure 10:
Left panel: X-ray 0.5-4.0 keV brighntess profile across the edge in
the relic region, extracted in the sector indicated in the left panel of Fig. 9.
Error bars are ![]() ![]() |
Open with DEXTER |
Using the available Chandra X-ray archival data of A 521, we checked if there is indeed a shock front in the cluster gas at the position of the relic. Shock fronts in clusters are observed very rarely. Only two unambigous examples, exhibiting both a sharp gas density edge and a clear gas temperature jump, have been discovered by Chandra so far, those in the merging clusters 1E 0657-56 (Markevitch et al. 2002) and A 520 (Markevitch et al. 2005). Their rarity is due to the fact that the viewing geometry and the moment of our observation must be favourable: a merger shock quickly moves to the cluster outskirts where it cannot be detected.
For A 521, two observations are available in the Chandra
public archive, one performed using ACIS-I and another using
ACIS-S, with 40 ks exposure each (Ferrari et al. 2006).
Unfortunately, the relic lies right at the S3 chip boundary in the
ACIS-S observation, so we could use only the ACIS-I observation
(OBSID 901) for our purpose. We cleaned the data and modelled the
detector background and instrumental spatial response as described
in Vikhlinin et al. (2005).
In Fig. 9, we show the resulting Chandra image in
the 0.5-4.0 keV energy range. The background was subtracted and the
image divided by the exposure map, and then smoothed with a
Gaussian. In the right panel, the
GMRT 610 MHz contours are overlaid on the same X-ray image. The image
reveals a clear brightness edge coincident with the outer
edge of the radio relic.
The radial brightness profile in Fig. 10 (left panel)
shows this edge more clearly. The profile was extracted from the unsmoothed
image in the sector shown in the left panel of Fig. 9, which
is centred on the centre of curvature of the relic and spans the angle covered
by the source. Discrete X-ray sources were excluded for the profile
derivation. Figure 10 (left panel) shows a
brightness edge at
,
which coincides
with the outer edge of the radio relic; we note that the drop
in the profile and increased error bar at
corresponds to an interchip gap. For comparison, in the right panel we plot
the X-ray brightness profiles extracted from
three other sectors of the cluster. None of these profiles show such
an edge, which is consistent with the visual inspection of the image
shown in Fig. 9.
The X-ray edge in Fig. 10 (left panel) has
the characteristic shape that corresponds to a projection of
a spherical density discontinuity. To quantify the
discontinuity, we fitted the brightness profile close to the
edge using a gas-density model of two power laws and an
abrupt jump, where the power laws and position and amplitude of the
jump were free parameters (as in Markevitch et al. 2000). A continuous
density profile (i.e., no jump, but a possible break), shown
by the red line in Fig. 11, is inconsistent
with the data at about .
If one assumes a 6 keV
gas temperature just inside the edge (Ferrari et al. 2006,
determined from this and the ACIS-S Chandra
observations) and converts brightness in the Chandra 0.5-4 keV band into plasma emission measure self-consistently, the
best-fit density jump (the model shown by green line)
corresponds to
.
However, as mentioned above,
such a shock would imply a pre-shock temperature of
0.4 keV, which is unlikely close to the centre of such a
hot cluster. A shock with M=2.3, predicted from the radio
spectrum (Sect. 7.2), corresponds to an edge shown by the
blue line in Fig. 11; this fit is only
away from the best fit, which we consider good
agreement.
Unfortunately, the accuracy of the existing Chandra observation is insufficient to measure the gas temperature on the faint side of the edge. Thus, we cannot rule out other interpretations of this feature, e.g., a cold front (Markevitch & Vikhlinin 2007). However, a cold front with such a density contrast would correspond to an outer gas temperature in excess of 15 keV - quite unlikely in a cluster of this X-ray luminosity. We conclude that the X-ray edge is most likely a shock front - found at a location and with an amplitude exactly as required to produce the radio relic by means of the shock acceleration mechanism.
![]() |
Figure 11: X-ray brighntess profile (same as left panel of Fig. 10) compared to different models. The green and red lines are the fit with and without a shock discontinuity, respectively. The blue line is the profile expected for a shock with M=2.3. |
Open with DEXTER |
In the following, we summarise the main results of our analysis:
Acknowledgements
We thank the staff of the GMRT for their help during the observations. GMRT is run by the National Centre for Radio Astrophysics of the Tata Institute of Fundamental Research. Thanks are due to E. Fomalont for his help with the use of the NRAO pipeline. We thank the anonymous referee for useful comments. S.G. and T.V. acknowledge partial support from the Italian Ministry of Foreign Affairs. This work has been partially supported by contracts ASI-INAF I/088/06/0, PRIN-MUR 2005, and PRIN-INAF 2005.
![]() |
Figure A.1:
Preliminary GMRT image (contours and grey scale) of radio
relic at 235 MHz. The resolution is 15.6
![]() ![]() ![]() ![]() |
Open with DEXTER |
The full
resolution image shown in Fig. A.1 was
produced using only the data from the first day
(i.e., 9 h). The image has a very high
sensitivity level (1
Jy b-1), and
the relic is clearly detected and imaged in its whole
extent. Given the high quality of this image we are
confident that the flux density given in Table 4 is
definitely reliable, with an error of the order of 5%.
The analysis of the entire dataset (i.e., from the
combination of the two observing days) will be presented
in a forthcoming paper.