Open Access
Issue
A&A
Volume 690, October 2024
Article Number A329
Number of page(s) 22
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/202451515
Published online 22 October 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

A variety of discrete and diffuse radio sources can be found in galaxy clusters, such as radio galaxies, radio halos, mini-halos, and radio relics (see Brunetti & Jones 2014; van Weeren et al. 2019 for reviews). These targets offer a chance to probe the evolution and dynamics of galaxy clusters and their members, along with the cosmic magnetism and the energy transfer processes on a range of scales.

Radio galaxies in clusters are typically of Fanaroff-Riley I-type (FRI; Fanaroff & Riley 1974) exhibiting two jets lacking hotspots at their tips, which are launched in opposite directions from the core at relativistic velocities and become subrelativistic at a distance of a few kpc. In galaxy clusters, FRIs can be reshaped due to the high-velocity (v ~ 1000 km s−1) motion of the host galaxy throughout the rarefied (nICM ~ 10−2−10−4 cm−3) thermal intracluster medium (ICM). Specifically, ram pressure (Pram ∝ v2nICM; Gunn & Gott 1972) can deflect the radio jets by angles ≳90° and generate narrow-angle tail (NAT) galaxies (e.g. Miley et al. 1972; Pfrommer & Jones 2011; Terni de Gregory et al. 2017; Pal & Kumari 2023). These galaxies exhibit a bright core (the head) and roughly parallel jets that rapidly diffuse as tails of length ~ 100–500 kpc. Ram pressure is also responsible for moderate deflections by angles ≪90° observed in some targets, possibly associated with merging subclusters; in this case, they form a wide-angle tail (WAT) galaxy (e.g. Owen & Rudnick 1976; Owen et al. 1978; Pinkney et al. 1993; Missaglia et al. 2019; O’Dea & Baum 2023). When the jets are not resolved (due to the intrinsic bending, projection effects, or resolution), NAT and WAT galaxies are generally referred to as head-tail (HT) galaxies. Owing to their origin, tailed radio galaxies are useful tracers of high-z and/or low-mass groups and clusters (e.g. Blanton et al. 2000, 2003; Giacintucci & Venturi 2009), which are scarcely detectable by current X-ray telescopes.

In typical tailed galaxies, the radio brightness decreases from the core along the tail due to particle ageing. Therefore, the spectral index steepens with the distance from the core (e.g. Cuciti et al. 2018; Botteon et al. 2021; Rudnick et al. 2021). However, departures from these trends have been observed in a number of targets (e.g. Jones & McAdam 1994; de Gasperin et al. 2017; Sebastian et al. 2017; van Weeren et al. 2017; Lal 2020; Müller et al. 2021; Rudnick et al. 2021; Edler et al. 2022; Riseley et al. 2022a; Lusetti et al. 2024). These observations suggest the presence of particle re-acceleration processes and/or amplification of magnetic fields. Furthermore, tailed galaxies, especially those extending throughout large fractions of the radii of clusters (e.g. Owen et al. 2014; Wilber et al. 2018; Srivastava & Singal 2020; Ignesti et al. 2022) are likely to release relativistic electrons and magnetic fields into the ICM (e.g. de Gasperin et al. 2017; van Weeren et al. 2017; Vazza & Botteon 2024). Therefore, tailed galaxies are interesting targets to probe the interplay between thermal and non-thermal components in galaxy clusters, and the complex re-acceleration mechanisms that are still poorly constrained (e.g. Brunetti & Jones 2014).

Abell 2142 (A2142) is a nearby and massive galaxy cluster in an intermediate dynamical state between relaxed and merging systems (Rossetti et al. 2013). It is characterised by complex dynamics (Liu et al. 2018) and is known to host peculiar discrete and diffuse radio sources (Venturi et al. 2017; Bruno et al. 2023; Riseley et al. 2024). In the present work, we focus on two morphologically interesting HT galaxies in A2142, labelled as T1 and T2 in Venturi et al. (2017). We aim to analyse their morphological and spectral properties, over a wide range of wavelengths (30 MHz-6.5 GHz), to search for possible interactions between the tails and the ICM and constrain their dynamics.

Throughout this paper, we adopt a standard ΛCDM cosmology with H0 = 70 km s−1 Mpc−1, ΩM = 0.3, and ΩΛ = 0.7. At the cluster redshift of z = 0.0894, we have 1″ = 1.669 kpc (or 1′ ~ 100 kpc). We adopted the convention on the spectral index α as defined from the flux density S(ν) ∝ να. The paper is organised as follows. In Sec. 2, we describe the galaxy cluster A2142. In Sec. 3, we present the radio and X-ray data, along with their processing. In Sec. 4, we report on the results of our analysis. In Sec. 5 we discuss scenarios to explain the properties of T1 and T2. In Sec. 6, we summarise our findings.

2 The galaxy cluster Abell 2142

A2142 (RAJ2000 = 15h58m20s, DecJ2000 = 27°14′00″) is a nearby (z = 0.0894) galaxy cluster of mass1 M500 = (8.8 ± 0.2) × 1014 M within a radius of R500 = 14.07 ± 0.70 arcmin (1408.5 ± 70.4 kpc at the cluster redshift; Planck Collaboration XXIV 2016; Tchernin et al. 2016). The galaxy cluster A2142 is located at the centre of the A2142 supercluster, which it is named after (Einasto et al. 2015; Gramann et al. 2015). In Fig. 1, we provide a multi-wavelength (optical, X-ray, radio) view of A2142 and label the various sources and features that we discuss below.

A2142 hosts about 900 member galaxies gathered in small groups and hierarchically organised in structures and substructures (Owers et al. 2011; Einasto et al. 2018; Liu et al. 2018). Furthermore, several groups are infalling towards the richest structure in the cluster centre (e.g. Owers et al. 2011; Eckert et al. 2017; Liu et al. 2018), which hosts the brightest cluster galaxy ‘BCG1’ (RAJ2000 = 239.5834, DecJ2000 = 27.2334, z = 0.09081, M* = 1.9 × 1011 M). The secondary brightest cluster galaxy, BCG2 (RAJ2000 = 239.5554, DecJ2000 = 27.2481, z = 0.0965, M* = 1.5 × 1011 M) is located at a projected distance of ~180 kpc from BCG1 and is thought to be the main member of a merging group. Minor mergers are likely responsible for the uncommon properties of the ICM, owing to their intermediate dynamical state in between that of relaxed cool cores and unrelaxed (major) mergers (e.g. Cavagnolo et al. 2009; Rossetti et al. 2013; Tchernin et al. 2016; Wang & Markevitch 2018; Cuciti et al. 2021).

Radio observations of A2142 with the Low Frequency Array (LOFAR) revealed diffuse emission from the ICM in the form of a hybrid radio halo (Fig. 1) consisting of three distinct components (Bruno et al. 2023). The brightest component (the ‘core’, ‘H1’) in the cluster centre has a roundish morphology and a diameter of ~200 kpc, the second component (the ‘ridge’, ‘H2’) is elongated towards south-east for ~400 kpc, and the third component (‘H3’) extends for ~2 Mpc embedding both H1 and H2 (Venturi et al. 2017; Bruno et al. 2023; Riseley et al. 2024). The origin of the hybrid halo is likely associated with turbulent re-acceleration triggered by mergers taking place on different spatial scales and/or timescales (see Bruno et al. 2023, for a detailed discussion).

Besides the hybrid halo, extended radio emission is associated with two prominent HT galaxies, T1 and T2, which are the focus of the present work. The host of T1 (RAJ2000 = 239.5596, Dec5J2000 = 27.2721, z = 0.09540, M* = 0.7 × 1011 M) is an elliptical galaxy likely being a member of the merging group associated with BCG2 (Einasto et al. 2018). Its optical spectrum2 shows the presence of weak emission lines (e.g. O[III], Hα) that indicate AGN activity. The core is active at radio wavelengths (e.g. Colla et al. 1972; Govoni et al. 2010; Terni de Gregory et al. 2017) and bright in the 0.5–10 keV band (L0.5–10 = 2.1 × 1042 erg s−1; Sun 2009), further supporting the association with an AGN. The host of T2 (RAJ2000 = 239.5870, DecJ2000 = 27.3337, z = 0.08953, M* = 2.1 × 1011 M) is an elliptical member galaxy of A2142. Its optical spectrum lacks emission lines that would confirm the presence of current AGN activity, while its radio and X-ray emission are discussed throughout this paper.

3 Observations and data reduction

In this section, we present the radio and X-ray data to study the HT radio galaxies in A2142. We first provide a brief overview of the radio data recently analysed in Bruno et al. (2023) and in Riseley et al. (2024) for the study of the cluster’s hybrid radio halo. The additional radio and X-ray data used in the present work are described in Sects. 3.13.5. The details of all radio and X-ray observations are reported in Tables 1 and 2, respectively.

Some of the radio data exploited in this work were recently presented in Bruno et al. (2023) to characterise the hybrid radio halo. These include data from LOFAR, the Giant Metrewave Radio Telescope (GMRT), the upgraded GMRT (uGMRT), and the Very Large Array (VLA). Furthermore, MeerKAT data have been recently used by Riseley et al. (2024) for a follow-up analysis. Here, we briefly summarise these observations, and refer to Bruno et al. (2023) and Riseley et al. (2024) for details on the telescope-specific setup and calibration strategies.

LOFAR observed A2142 for 16 and 32 hours with the Low Band Antenna (LBA) and High Band Antenna (HBA) arrays operating at 30–70 MHz and 120–168 MHz, respectively. The GMRT observed A2142 for 5 hours at 305–340 MHz (see also Venturi et al. 2017); furthermore, 3-hour observations were carried out with the uGMRT at 300–500 MHz (band 3). Mosaicked pointings (each having a field of view of ~30′) on A2142 were obtained with the Very Large Array (VLA) at 1–2 GHz (L-band) in C-array and D-array configurations, for a total of 2 hours. MeerKAT observed A2142 at 872–1712 MHz (L-band) for 5.5 hours as part of the ‘MeerKAT-meets-LOFAR mini-halo census’ project (Riseley et al. 2022b, 2023).

thumbnail Fig. 1

Composite RGB image of A2142: radio (LOFAR at 143 MHz, 9″ × 6″) in red, optical (DSS-2, r-filter) in green, and X-rays (XMM-Newton, 0.7–1.2 keV) in blue. Labels indicate the tailed radio galaxies (T1, T2), the primary and secondary brightest cluster galaxies (BCG1, BCG2), and the two components (H1, H2) of the radio halo. The dashed white contours represent the 3σ level of the 73″ × 66″ LOFAR image at 143 MHz reported in Bruno et al. (2023), showing the third component (H3) of the radio halo. The yellow arrows indicate the location of the NW, central, and SE cold fronts reported by Rossetti et al. (2013).

Table 1

Details of the radio data analysed in this work.

3.1 GMRT data

Archival GMRT observations of A2142, first presented by Venturi et al. (2017), are available in the 225–240 MHz and 590–625 MHz bands, for 6 and 5 hour on source, respectively. The total bandwidths of 16 and 32 MHz are split into 128 and 256 channels, respectively. The sources 3C286 and 3C48 were used as absolute flux density scale calibrators.

Following the same procedure described in Bruno et al. (2023) for the 305–340 MHz band data, we reprocessed the 225–240 and 590–625 MHz band data by means of the Source Peeling and Atmospheric Modeling (SPAM) automated pipeline (Intema et al. 2009), which corrects for ionospheric effects by deriving directional-dependent gains from bright sources across the field of view. We achieved noise levels of −180 μJy beam−1 at 13″ × 10″ and ~40 μJy beam−1 at 6″ × 4″ for the 234 and 608 MHz datasets, respectively, which are consistent with those reported in Venturi et al. (2017).

Table 2

Details of the Chandra X-ray data analysed in this work.

3.2 uGMRT band 5 data

New observations of A2142 were performed with the uGMRT in the 1050–1450 MHz (band 5) frequency range in March 2018. The total 400 MHz bandwidth is split into 8192 channels of width ~49 kHz each. The cluster was covered in its full extent with five different pointings (each having a field of view of ~27′), for a total of 7 hours. For the study presented in this paper, only the three pointings covering the northern part of the cluster (~40′ in total) are considered (namely ‘A2142_2′, ‘A2142_3′, ‘A2142_5′) for a total of 4 hours. The sources 3C286 and 1602+334 were used as amplitude and phase calibrators, respectively.

Data reduction of band 5 observations with SPAM has not been tested in depth. However, direction-dependent corrections are negligible for these gigahertz frequency and small field of view data; thus, we did not used SPAM. To process these data, we split the total bandwidth in 6 sub-bands of ~67 MHz each and carried out a standard data reduction with the Common Astronomy Software Applications (CASA v. 6.5; McMullin et al. 2007) by iteratively performing flagging of Radio Frequency Interference (RFI), and bandpass, amplitude, and phase calibrations for each sub-band and pointing. We then recombined the calibrated sub-bands of each pointing to perform rounds of phase and phase plus amplitude self-calibration. The self-calibrated pointings were imaged separately with WSClean v. 2.10 (Offringa et al. 2014; Offringa & Smirnov 2017) with multi-frequency and multi-scale synthesis options. After correcting each of the three images by the corresponding primary beam attenuation at the central frequency of 1250 MHz and convolving them to the same resolution, they were combined by means of the tool lm.makemosaic in CASA to produce a single mosaic image. At a resolution of 2.5″, the final noise level is in the range ~25–40 μJy beam−1.

3.3 VLA L-band and C-band radio data

As first presented by Terni de Gregory et al. (2017), the head-tail galaxy T1 was studied with the VLA in A configuration at 1–2 GHz (L-band) and 4.5–6.5 GHz (C-band) for 15 minutes on source in each band. In both observations, the sources 3C286 and 1602+3326 were used as absolute flux density and phase calibrators, respectively. The data were recorded in 16 spectral windows of 128 MHz each.

The field of view of these VLA observations include the core of T2, which we aim to study. We reprocessed the L-band and C-band data in CASA following standard calibration procedures (see Sec. 3.2), performing an additional cycle of phase self-calibration. Our data processing improved the quality of the images with respect to those reported in Terni de Gregory et al. (2017) in terms of noise (improvement by factors of ~1.6 and ~1.2 in L-band and C-band, respectively). We reached a noise level of ~27 μJy beam−1 in L-band at ~1″ resolution and ~11 μJu beam−1 in C-band at ~0.3″ resolution.

3.4 Radio imaging

The imaging of all the radio data used in this study was carried out with WSClean v. 2.10 (Offringa et al. 2014; Offringa & Smirnov 2017) to account for wide-field, multi-frequency, and multi-scale synthesis. For both VLA and uGMRT, mosaicked observations were imaged separately and then properly combined following the procedure described in Sect. 3.2 for uGMRT.

In the following, uncertainties on the reported radio flux densities S are computed as: ΔS=(σ2Nbeam )+(ξcal S)2,$\[\Delta S=\sqrt{\left(\sigma^2 \cdot N_{\text {beam }}\right)+\left(\xi_{\text {cal }} \cdot S\right)^2},\]$(1)

where Nbeam is the number of independent beams within the considered region, and ξcal is the calibration error. We assumed standard calibration errors of ξcal = 10% for LOFAR (Shimwell et al. 2022; de Gasperin et al. 2023), ξcal = 7%, 6%, 6%, 5%, 5% for GMRT at 234 MHz, 323 MHz, 407 MHz, 608 MHz, and 1250 MHz, respectively (Chandra et al. 2004), ξcal = 5% for VLA and MeerKAT in L-band (Perley & Butler 2013), and ξcal = 3% for VLA in the C-band (Perley & Butler 2013).

3.5 Chandra X-ray data

In Bruno et al. (2023), we analysed deep Chandra observations of A2142. These consist of four pointings of 187 ks in total that mainly cover the central regions of the clusters. In the present work, we considered three additional pointings covering the northern and north-eastern regions of A2142, in the direction of T2. These observations were carried out in 2014, in VFAINT mode, with both ACIS-I and ACIS-S CCDs, and were first presented in Eckert et al. (2017). A summary of all Chandra data considered in this work is reported in Table 2.

As for the other pointings, we reprocessed the additional data by means of CIAO v. 4.13, with CALDB v. 4.9.4. After extracting light curves in source-free regions, soft proton flares were filtered out with the lc_clean algorithm, leaving a clean time of 169.1 ks. Overall, the seven pointings provide a total clean time of 356 ks. In this work, we used these data to investigate the local conditions of the ICM towards T1 and T2 with a combination of resolution and signal-to-noise ratio (S/N) that depends on the depth of the pointings covering the same regions. We refer to Bruno et al. (2023) for the description of background treatment for imaging and spectral analysis.

thumbnail Fig. 2

Radio images and contours labelling the sub-regions and features discussed in the text. Left panels: 143 MHz images (see Figs. 3 and 4) of T1 (top) and T2 (bottom). Right panels: Pan-STARSS optical (composite i, r, g filters) images zoomed towards T1 (top) and T2 (bottom). Radio contours at 143 MHz (5σ, in red) and 1250 MHz (5σ for T1 and 3σ for T2, in cyan) are overlaid.

4 Results

4.1 Radio morphology and sub-regions

Within the tailed radio galaxies T1 and T2 we define sub-regions that are discussed throughout this section and labelled in Fig. 2, along with optical overlays from the Panoramic Survey Telescope & Rapid Response System (Pan-STARSS; Flewelling et al. 2020). Radio images of T1 and T2 are presented in Figs. 3 and 4 at different frequencies and resolutions (see details in Table 3).

The host of T1 is located at a projected distance of ~270 kpc from the brightest cluster galaxy. Starting from the head, T1 extends from east to west, for a total projected length of ~5.5′ (at 143 MHz), corresponding to ~550 kpc at the cluster redshift. The sub-region T1-A, which corresponds to the initial ~200 kpc, is the brightest part of the source and is well imaged at all frequencies. In this region, the width of tail is ~35 kpc, thus suggesting that the whole structure (made by the two merged jets into the tail) is highly collimated. For the subsequent ~150 kpc (T1-B), the tail doubles its width up to ~70 kpc and shows clear brightness fluctuations, which we will refer to as wiggles due to their oscillating pattern. We note the presence of a galaxy (RAJ2000 = 239.5267, DecJ2000 = 27.2742, z = 0.08684) at the transition from T1-A and T1-B (Fig. 2). At the end of sub-region T1-B, a bright compact component is visible at higher frequencies (its location is indicated by the purple circle in Fig. 2), which is likely a compact radio source seen in projection. For the last ~200 kpc (T1-C), the tail follows a straight path. Interestingly, additional emission (T1-D) is well detected at 143 MHz (and partly visible at 50 and 323 MHz in lower resolution images). A thin, ~50 kpc long filament connects T1-C with an arc-shaped structure extending for ~50 kpc and 200 kpc along east-west and north-south, respectively. If the arc were the termination of the tail, the total length of T1 would be ~650 kpc (in lower resolution images, more emission from T1-D is recovered, reaching ~700 kpc in total). We aim to shed light on this feature with a spectral analysis in Sect. 4.4.

The morphology of T2, which extends along the SE-NW axis for ~400 kpc, is more peculiar than classical tailed galaxies. The optical counterpart of T2 is located at a projected distance of ~600 kpc from the brightest cluster galaxy, and hosts a weak radio core (T2-D; see also Fig. 2, bottom right panel) that is resolved from the rest of the source at 608 and 1250 MHz only (this is likely due to a favourable combination of higher resolution and lower sensitivity to extended components than other images). The core is the base of the first sub-region (T2-A), which has a light bulb shape of width ~75 kpc and length ~100 kpc. As highlighted by our high resolution (2.5″) 1250 MHz data, the light bulb ends with a sharp edge (see also Fig. 2). The global morphology of T2-A may be presumably due to the backward bending of the jets by the ram pressure, but these are not resolved by any of our images. A second sub-region (T2-B) is defined by a choking, that is the abrupt shrinking of the width of the tail and a drop in the radio surface brightness. T2-B extends for ~150 kpc, has a fairly constant width of ~75 kpc, and exhibits a single bright spot. The last sub-region (T2-C) is defined by the spread of the tail into a diffuse and filamentary plume of length ~150 kpc, maximum width ~200 kpc, and non-uniform brightness. The plume is detected in our images only at low frequencies (ν < 407 MHz), thus suggesting a very steep spectral index for this region. Interestingly, the western part of the plume is bent towards SW (see also LBA images at lower resolution in Bruno et al. 2023).

Table 3

Summary of the radio images shown in Fig. 3.

thumbnail Fig. 3

Radio images of T1. Details on resolution and noise are reported in Table 3. The overlaid contour levels are [±5, 10, 20,...] × σ. From top left to bottom right: LOFAR LBA at 50 MHz, LOFAR HBA at 143 MHz, GMRT at 234 MHz, GMRT at 323 MHz, uGMRT at 407 MHz, GMRT at 608 MHz, uGMRT at 1250 MHz (cropped to a smaller area for inspection purposes), MeerKAT at 1284 MHz, VLA at 1380 MHz, and VLA at 1810 MHz.

thumbnail Fig. 4

Radio images of T2. Details on resolution and noise are reported in Table 3. The overlaid contour levels are [±5, 10, 20, ...] × σ. From top left to bottom right: LOFAR LBA at 50 MHz, LOFAR HBA at 143 MHz, GMRT at 234 MHz, GMRT at 323 MHz, uGMRT at 407 MHz, GMRT at 608 MHz, uGMRT at 1250 MHz, MeerKAT at 1284 MHz, VLA at 1380 MHz, and VLA at 1810 MHz.

thumbnail Fig. 5

Radio surface brightness profiles of T1 and T2 normalised to the peak value within T1-A and T2-A, respectively. Left panels: profiles of T1 (middle) and T2 (bottom) at 143 MHz (resolution 9″ × 6″). The sampling boxes are shown in the top panels and have a width of 10″. Right panels: comparison of profiles at 50 (brown), 143 (yellow), 323 (green), 608 (light blue), and 1284 (magenta) MHz. The resolution of the images and the width of the sampling boxes (not shown) are 14″. In middle and bottom panels, the vertical lines indicate the boundaries of sub-regions as in Fig. 2.

4.2 Radio surface brightness fluctuations

To investigate the brightness fluctuations within each sub-region, we computed the surface brightness profiles of T1 and T2. In the left panels of Fig. 5 we report the profiles (black data points) measured from the LOFAR HBA image (Figs. 3 and 4) in boxes of width 10″. Each data point of T1 and T2 is normalised by the peak value within T1-A and T2-A, respectively. Analogously, in the right panels of Fig. 5 we compare the normalised profiles at 50 (brown), 143 (yellow), 323 (green), 608 (light blue), and 1284 (magenta) MHz, as measured from images convolved at the same resolution of 14″. The surface brightness profiles are reflective of the sub-regions that we defined in Sect. 4.1 by visual inspection, and reveal interesting features discussed below.

The absolute peak value of T1 is coincident with the radio core (first sampling box). In T1-A, the brightness at 143 MHz rapidly decreases (down to a factor of ~10) with the distance from the core. A discontinuity is visible at ~100 kpc, where the tail first deviates from its straight path and is slightly compressed. In T1-B (the region of the wiggles), the brightness is enhanced instead of declining, but such growth becomes progressively shallower with the increasing frequency. In T1-C, the declining trend resumes, but we also report the presence of a moderate peak at ~500 kpc at low-ν. The last peak in T1-D is associated with the arc, which is detected only at 50, 143, and 323 MHz.

The surface brightness profile of T2 is largely unusual for tailed galaxies. The core of T2 is weak (see also Sec. 4.3) and is not coincident with the absolute peak of emission. Overall, the normalised profile exhibits three peaks of similar relative amplitude, corresponding to the light bulb (T2-A), the main body of the tail (T2-B), and the plume (T2-C). The first peak is associated with the bright spot at the edge of T2-A. The choke that separates T2-A and T2-B is identified as a sharp discontinuity in our profiles at a distance of ~100–150 kpc. Within T2-B, the surface brightness increases with the distance up to ~200 kpc. Finally, the emission of the plume produces the last peak, which becomes progressively shallower with the increasing frequency.

thumbnail Fig. 6

Radio spectra of T1 and T2. Left and right panels refer to T1 and T2, respectively. Data points are the flux densities reported in Table 4, measured from the regions (including the green circle used for the core of T2) shown in the top panels following the same colour-code. Flux densities were fitted with single (solid lines) or double (solid plus dotted lines) power laws (see Table 5).

4.3 Integrated radio spectra

To derive the integrated radio spectra of T1 and T2, we imaged all the datasets (except VLA data in A-array) with a common uv-range of 350λ − 16kλ. The chosen minimum baseline length provides more uniform uv-coverage of our data at short spacings. The obtained images were convolved at the same resolution of 14″. We report our flux density measurements in Table 4; the corresponding radio spectra are shown in Fig. 6, and the obtained spectral indices are summarised in Table 5.

We measured the total flux densities of T1 (excluding T1-D) in a box of size 6.0′ × 1.3′ (blue box in Fig. 6) and fitted the data points with a single power law. We found that a single power law of a slope of α = 0.87 ± 0.01 (blue line in Fig. 6) can well reproduce the integrated spectrum of T1 from 50 to 1810 MHz. In addition, we obtained the radio spectrum of each sub-region by measuring the flux densities within boxes of size 2.0′ × 1.0′, 1.5′ × 1.0′, 2.0′ × 1.0′, and 1.0′ × 2.0′ for T1-A (red), T1-B (orange), T1-C (cyan), and T1-D (green), respectively. Fig. 6 clearly shows that the total flux density of T1 is dominated by T1-A at all frequencies. For T1-A, data points can be described by a single power law of a slope of α = 0.68 ± 0.02. However, we find evidence of spectral breaks for T1-B and T1-C, as single power laws cannot fit our measurements. We thus considered double power laws, with a fixed break at 234 MHz, as this appears to be roughly the frequency where the spectrum steepens. The details on the fitted spectral indices between 50 and 234 MHz (solid lines) and between 234 and 1284 MHz (dotted lines) are reported in Table 5. The arc is detected (above 3σ) only by LOFAR and GMRT at 323 MHz; therefore, we did not attempt to fit two power laws, but the poor fit (χred 2=7.5)$\[\left(\chi_{\text {red }}^2=7.5\right)\]$ suggests the existence of a break for this sub-region as well.

Similarly to T1, we obtained the total flux density of T2 in a box of size 3.0′ × 4.5′ (blue box in Fig. 6), but a single power law does not reproduce all our data points (χred 2=5.9)$\[\left(\chi_{\text {red }}^2=5.9\right)\]$. The flux densities of T2-A (red), T2-B (orange), and T2-C (cyan) were measured in an ellipse of axis 0.8′ × 1.1′, a box of size 1.0′ × 1.5′, and an ellipse of axis 1.8′ × 2.4′, respectively, then fitted with two power laws having a fixed break at 407 MHz (see details in Table 5). Even though we notice that the low-ν spectrum of T2-B is not accurately fitted (χred 2=6.5)$\[\left(\chi_{\text {red }}^2=6.5\right)\]$, likely due to contamination from the plume in overlapping areas, our analysis shows that the spectrum of T2 progressively steepens from the inner to the outer sub-regions. Moreover, we produced additional radio images (including VLA data in A-array) with specific combination of uv-range and weighting schemes to maximise the resolution (<5″) of our images, and derive the spectrum of the core of T2 (T2-D, green line). We used the imfit task in CASA to derive the peak value of a Gaussian fit to the core at each frequency. These values provide a fitted spectral index of α = 0.69 ± 0.02 from 143 to 5500 MHz.

Table 4

Flux densities of T1, T2, and corresponding sub-regions.

Table 5

Integrated spectral index of T1 and T2.

4.4 Resolved spectral properties

The analysis of the integrated radio spectra of T1 and T2 in Sec. 4.3 has shown that the spectral index is not constant across our targets. Under the hypothesis of pure radiative losses, a gradual steepening of the spectral index is expected along the tail, whereas a sudden spectral flattening may suggest re-energising mechanisms. To determine the spectral trend along the tails, we produced spectral index maps by combining sets of radio images (see details in Table 6) and setting a minimum flux density threshold of 5σ at each frequency. In Fig. 7 we report the spectral index maps at 50–143 and 143–323 MHz at 30″ (~50 kpc), and at 143–323 MHz and 608–1284 MHz at 14″ (~25 kpc); the associated error maps are shown in Fig. A.1. By measuring the spectral index within boxes of width equal to the beam size, we derived the corresponding spectral index profiles that are shown in Fig. 8. As typically observed along the lobes of FRI galaxies due to ageing of the emitting particles, on average both T1 and T2 exhibit a progressive steepening of the spectral index along the tail. We discuss the spectral trends within each sub-region in the following paragraphs.

Within T1-A, the spectral index ranges from α − 0.5 (in the core) up to α ~ 1.5, with steeper values for higher frequency pairs. The profiles steepen with the distance with approximately constant slopes (even though the trend is flatter at lower frequencies). Along T1-B and T1-C, the spectra further steepen, reaching ultra-steep (a ≳ 1.5) values of α ~ 2 and α ~ 2.5 between 143–323 MHz and 608–1284 MHz, respectively. In the lower-frequency regime (<323 MHz), the constant steepening trend with the distance is retained, whereas deviations in the form of flatter and steeper features are clearly visible at 608–1284 MHz both in T1-B and T1-C. Despite the different resolutions (14″ and 30″), we notice that the trends of the 143–323 MHz maps are consistent. Beyond T1-C, emission is detected only below 323 MHz. In T1-D (the arc) the spectral index is α − 1.5 and α − 2.5 between 50–143 and 143–323 MHz, respectively. The constant steepening of the spectral index is still preserved (with moderate deviations at 143–323 MHz), thus allowing us to confidently conclude that the arc is the oldest part of T1. As a consequence, the total projected length of T1 is ~700 kpc.

In T2-A, the inner regions exhibit a spectral index that is flatter (α − 0.5) at 50–143 MHz than that at higher frequencies (α − 0.7). For both T2-A and T2-B, the spectral index steepens with the distance with roughly constant and smooth trends at all wavelengths. The spectrum of T2-B reaches values of α −0.8, α − 1.1, and α − 1.5 at 50–143 MHz, 143–323 MHz, and 608–1284 MHz, respectively. Studying the profile of T2-C is not trivial due to the expansion of the tail into the plume and its western bending at low frequencies. Along our sampling boxes, the measured spectral indices are ultra-steep (α ~ 1.5–2.5) for all maps. The smooth trend is retained below 323 MHz, while a rapid steepening is observed at 608–1284 MHz.

Table 6

Parameters used to produce images for the spectral index maps in Fig. 7.

thumbnail Fig. 7

Spectral index maps of T1 and T2 computed with a threshold of S ≥ 5σ. In all the panels, the contour levels are [±5, 10, 20, ...] × σ of the lowest-frequency image in each combination. Top: spectral index maps at 30 between 50–143 MHz (left) and 143–323 MHz (right). Bottom: spectral index maps at 14″ between 143–323 MHz (left) and 608–1284 MHz (right). The corresponding error maps are shown in Fig. A.1.

4.5 Testing the radiative ageing

Relativistic particles are ejected from the core of the radio galaxy and then radiatively age along the tail due to synchrotron and inverse Compton losses. By assuming a constant injection rate, we would expect to observe a progressive decline of the flux density and a steepening of the spectral index with the increasing distance. However (as described in Sects. 4.2, 4.4), we found complex distributions of the surface brightness and spectral index within each sub-region of T1 and T2. These suggest either local deviations from a pure ageing scenario or a varying particle injection rate. In this section, we aim to test ageing models in detail.

Ageing models depend on the form of the initial electron energy distribution: N(E)Eδinj $\[N(E) \propto E^{-\delta_{\text {inj }}}\]$, where δinj = 2αinj + 1 is the population injection index and αinj is the spectral index at age t = 0, based on the assumption of Fermi I acceleration mechanism. They also depend on the spatial distribution of the magnetic field (B(r)), which are not known a priori. Models assuming a single injection event are the Kardashev-Pacholczyk (KP; Kardashev 1962; Pacholczyk 1970), Jaffe-Perola (JP; Jaffe & Perola 1973), and Tribble-Jaffe-Perola (TJP; Tribble 1993). Both the KP and JP models assume a uniform magnetic field, but differ in terms of the treatment of the pitch angle θp. In the KP model, a constant and isotropic θp is assumed throughout the entire lifetime of the electrons; the JP model considers electron scattering, which leads to isotropic θp on short-time scales only, and thus assumes a time-averaged pitch angle. As a consequence of the assumption on θp, energetic electrons with small pitch angles can live indefinitely for the KP model if emitting synchrotron radiation only, whereas an exponential cut-off arises in the electron distribution at high energy for the JP model (e.g. Hardcastle 2013). The TJP model is based on the JP model, but introduces Gaussian spatial fluctuations of the magnetic field around a central value (B0). The high-energy cut-off is shallower in the TJP model than in the JP model, thus allowing particles to live longer in non-uniform magnetic fields.

In the following, we explain how we tested the KP, JP, and TJP ageing models. We considered a specific value for the magnetic field that minimises the radiative losses and maximises the lifetime of the source; namely: B0=BCMB/3$\[B_0=B_{\mathrm{CMB}} / \sqrt{3}\]$, where BCMB = 3.25(1 + z)2 μG is the equivalent magnetic field of the cosmic microwave background (CMB). This yields B0 ~ 2.2 μG for both T1 and T2 at the cluster redshift. The injection spectral index was derived by means of the findinject task of the Broadband Radio Astronomy ToolS (BRATS3; Harwood et al. 2013, 2015) software. It fits the radio spectrum of the target from multi-frequency radio images and outputs the particular value of αinj that minimises the distribution of χ2 for the fitted model; we obtained αinjT1=0.51±0.01$\[\alpha_{\mathrm{inj}}^{\mathrm{T} 1}=0.51 \pm 0.01\]$ for T1 and αinjT2=0.72±0.01$\[\alpha_{\mathrm{inj}}^{\mathrm{T} 2}=0.72 \pm 0.01\]$ for T2, in agreement with the measured spectral index of the radio cores (see Sects. 4.3, 4.4), thus suggesting that their spectra are representative of the injection distributions.

Radio colour-colour diagrams (RCCDs; Katz-Stone et al. 1993) are diagnostic plots meant to probe the local shape of radio spectra computed from two pairs of frequencies. These spectra are independent of the magnetic field and possible adiabatic expansion and compression (these can only shift the spectrum in frequency and affect the age), thus allowing us to test the theoretical ageing models. In RCCDs, the one-to-one line represents a power law spectrum with α = αinj, whereas data points lying below the one-to-one line (α > αinj) indicate particle ageing. By considering the spectral index maps and sampling boxes as in Sec. 4.4, we obtained the RCCDs shown in Fig. 9. We overlaid the theoretical ageing curves (dotted lines) as obtained from the kpdata, jpdata, and tribbledata tasks in BRATS under the assumptions discussed above.

For the low-resolution (30″) and low-ν (50–143–323 MHz) set of images of T1 (upper left panel in Fig. 9), both the JP and TJP models reproduce the emission of the source in each sub-region, whereas the KP model cannot describe data points beyond ~400 kpc (T1-C and T1-D). With 14″-resolution (upper right panel), the observed 143–323–608–1284 MHz spectral distribution can be scarcely described by the models in T1-A (inner ~100 kpc) and progressively increasing deviations are found in T1-B and T1-C. Although a steeper injection index (αinj ~ 0.6–0.8) would shift the ageing curves towards our measurements, models still fail to reproduce data points at large distance from the core. We will further discuss the results and implications of the RCCDs for T1 in Sects. 5.1 and 5.2.

The bottom panels in Fig. 9 report the complex distribution of the data points in the RCCD for T2. In the low-ν (50–143–323 MHz) set at 30″ (left), data points align roughly parallel to the one-to-one line (except for an outlier associated with the plume). In the high-ν (143–323–608–1284 MHz) set at 14″ (right), the KP, JP, and TJP models can approximately reproduce the data points associated with T2-A, but prominent deviations are found for T2-B and T2-C. In summary, none of the considered models can entirely describe the spectral behaviour of T2 (see further discussion in Sec. 5.5).

thumbnail Fig. 8

Spectral index profiles as measured from maps in Fig. 7 for T1 (left) and T2 (right). The sampling boxes have beam-size width and are shown in the upper panels. The vertical lines indicate the boundaries of sub-regions as in Fig. 2.

thumbnail Fig. 9

Radio colour-colour diagrams of T1 (top) and T2 (bottom). Data points are colour-coded based on their increasing distance from the core from blueish to reddish colours. The black dashed line (one-to-one line) indicates a power law with α = αinj. The dotted lines represent the theoretical KP (orange), JP (blue), and TJP (red) ageing curves obtained with αinj = 0.51 (T1), αinj = 0.72 (T2), B0 = 2.2 μG (T1, T2) as input. The frequency pairs and resolution are the same as in Fig. 8.

4.6 Local ICM conditions

Throughout this section, our aim is the search for direct evidence of interplay between radio emission and local ICM conditions. By means of XMM-Newton images in different energy bands, Rossetti et al. (2013) produced the projected temperature map that we show in the upper left panel of Fig. 10 as it provides a useful overview of the ICM temperature over the whole cluster. Nevertheless, this map does not allow us to probe the peripheral regions with sufficient resolution. In this respect, following the same procedure described in Bruno et al. (2023), we produced maps of projected thermodynamic quantities towards T1 and T2 by performing a spectral analysis with our Chandra data. CONTBIN4 (Sanders 2006) was used to bin the 0.5–2 keV exposure-corrected Chandra image in regions with a minimum S/N = 40. We extracted the spectra in each region of the event and blanksky files, subtracted the background from the ICM emission and jointly fitted the resulting spectra in XSPEC (Arnaud 1996) with an absorbed thermal plasma component (phabs × apec) by fixing the Galactic hydrogen column density in the direction of A2142 and the ICM metal abundance to values of NH = 3.8 × 1020 cm−2 and Z = 0.28 Z (Markevitch et al. 2000). The Cash statistics (Cstat; Cash 1979) was adopted for fitting.

Each spectrum provides values of temperature kT (in units of keV) and normalisation 𝒩 (in units of cm−5), which is proportional to the squared numerical density integrated over the volume. We derived the (projected) pressure as p=kT×(NA)12[keV cm5/2 arcmin1],$\[p=k T \times\left(\frac{\mathcal{N}}{A}\right)^{\frac{1}{2}} \quad\left[\mathrm{keV} \mathrm{~cm}^{-5 / 2} ~\operatorname{arcmin}^{-1}\right],\]$(2)

where A is the area of each region (in units of arcmin2). By propagating errors, uncertainties on p are computed as Δp=p(ΔkTkT)2+14(ΔNN)2.$\[\Delta p=p \sqrt{\left(\frac{\Delta k T}{k T}\right)^2+\frac{1}{4}\left(\frac{\Delta \mathcal{N}}{\mathcal{N}}\right)^2}.\]$(3)

In Fig. 10, we report the temperature and pressure maps in the direction of T1 (top right panels) and T2 (middle panels), respectively, and the corresponding profiles (bottom panels) as measured from regions that are labelled with green and magenta circles. While the chosen high S/N ensures an accurate spectral fitting, the extraction regions are relatively large and limit our analysis to minimum spatial scales of −30 kpc for T1 and ~100 kpc for T2.

For T1, we found drops in temperature and pressure from regions 4 to 5, which correspond to the transition from T1-A to T1-B. We further probed these trends by extracting and fitting the X-ray surface brightness with a broken power law by means of pyproffit5 (Eckert et al. 2020). As shown in Fig. 11, the profile is dominated by the presence of the prominent NW cold front, which we detect as a density jump of C = 1.80 ± 0.14 (Cstat/dof = 45/49); this is consistent with the value reported by Rossetti et al. (2013). Furthermore, we observe a peak at a distance of ~1′ (indicated by the red arrow in the profile), which is co-spatial with the transition from T1-A to T1-B. This feature is likely responsible for the observed kT and p drops and is suggestive of local ICM compression, but we were not able to constrain its nature through fitting procedures. While our analysis is inconclusive, it is unlikely that such compression is driven by a shock, as its passage along the tail would have left signatures in the spectral index distribution with a discontinuity between T1-A and T1-B. However, we did not observe any such features (Fig. 8).

For T2, both kT and p profiles are roughly continuous within errors. Even though we did notice a tentative (significance ≳1σ) drop in pressure between T2-A and T2-B, we cannot draw any solid conclusion on possible discontinuities from these plots. A complementary analysis (not shown) of the X-ray surface brightness profile in the direction of T2 with high spatial resolution (ranging from 2″ to 12″) is consistent with continuous trends of the thermodynamic properties. We notice that modelling of the background (e.g. Gastaldello et al. 2007; Bartalucci et al. 2014; Botteon et al. 2018; Bruno et al. 2021) may provide more accurate spectral results rather than subtraction, especially for the cluster outskirts (which is the case of T2). However, this is beyond the aim of the present work.

thumbnail Fig. 10

Projected thermodynamic maps. Top left: temperature map from XMM-Newton photometric analysis covering the whole cluster (adapted from Rossetti et al. 2013). Top right: temperature and pressure maps from Chandra spectral analysis towards T1. Bin ‘4’ roughly marks the transition T1-A/T1-B. Middle: temperature and pressure maps from Chandra spectral analysis towards T2. Bin ‘1’ includes T2-A and the choke, bin ‘2’ includes T2-B. Bottom: temperature and pressure profiles for T1 (green) and T2 (magenta), as measured from the bins that are labelled in the corresponding maps with circles. The white dots in the Chandra maps indicate point sources excluded from the analysis. The LOFAR HBA radio contours are overlaid in black.

thumbnail Fig. 11

X-ray surface brightness analysis in the direction of T1. Left: Chandra X-ray image in the 0.5–2 keV band with overlaid radio contours. The surface brightness was measured from the red sector in annuli of width 2″. The fitted position of the discontinuity is shown within the sector, and is coincident with the NW cold front. Right: data points fitted with a broken power law (blue line). The red arrow indicates the peak discussed in Sec. 4.6.

5 Discussion

In the previous sections, we analysed the morphological and spectral properties of T1 and T2. In the following, we discuss possible scenarios explaining the observed features.

5.1 On the discrepancy of RCCD analysis for T1

Through the RCCD analysis in Sec. 4.5 we showed that the observed spectral distribution of T1 at low frequency (≲300 MHz) can be reproduced by the JP and TJP models under simple assumptions with respect to the magnetic field (Fig. 9). Nevertheless, an inconsistent spectral trend is measured at higher frequencies. In RCCDs, data points are expected to follow the same spectral shape in a pure radiative ageing scenario, regardless of the considered frequency pairs. Therefore, the observed inconsistency requires a deeper understanding.

Offsets in the flux density scale and calibration artefacts can systematically shift the data points. The spectra of compact sources in the field do not reveal clear offsets in our datasets, and although images at 608 MHz show spurious emission around the tail (~4σ significance, likely generated by self-calibration errors), this is not driving the observed spectral distribution in the RCCD. A combination of projection effects and mixing of particles with different energies within each sampling region can broaden the intrinsic spectral distribution, mimicking our observed trend. While this hypothesis is plausible, the discrepancy among RCCDs is preserved when using images produced with different weighting schemes and resolutions, changing the considered frequency pairs, and sampling with beam-size circular regions to reduce possible transverse mixing across the boxes. Physical phenomena, such as compression/expansion and re-energising, can also alter the spectral distribution in RCCDs. However, compression and expansion are frequency-independent processes, while standard re-energising mechanisms via shock and turbulence predominantly affect the low-energy particles emitting at low-frequency. In this respect, unknown exotic processes favouring the re-acceleration of high-energy against low-energy particles would be necessary to solve the RCCD discrepancy.

In summary, we did not identify any obvious systematic effects or physical conditions explaining the observed spectral discrepancy. The possible role of subtle calibration effects should be investigated through independent reprocessing and methods. Mixed (either intrinsic or projected) energy distributions are plausible, but it is a challenge to further disentangle them.

thumbnail Fig. 12

Radiative age of T1. Top: radiative age map at 30″ obtained by fitting radio images at 50, 143, and 323 MHz with the TJP model (αinj = 0.51, B0 = 2.2 μG). Typical errors are ~20 Myr. The sub-regions are indicated by horizontal bars. Bottom: age profile measured in boxes of beam-size width. The boundaries of the sub-regions are indicated with dashed vertical lines. The blue solid line is the result of a linear fit to the data points.

5.2 Dynamics of T1

In Sect. 5.1, we discussed possible solutions to the RCCD discrepancy of T1. Throughout this section, we only consider the low frequency RCCD in Sect. 4.5, as we showed that the measured spectral trend is predicted by the JP and TJP ageing curves.

Following up on the results of the RCCD, we performed a pixel-by-pixel fitting of a TJP spectrum to our 50, 143, and 323 MHz images at 30″-resolution by means of BRATS, and obtained the radiative age map of T1 that is shown in the top panel of Fig. 12. The lower panel of Fig. 12 reports the corresponding age profile. Fitted ages are in the range t ~ 50–150 Myr in T1-A, t ~ 150–200 Myr in T1-B, t ~ 200–300 Myr in T1-C, and t ~ 350 Myr in T1-D, with typical errors of ~20 Myr.

We fit the age profile data points with a linear relation, which provides the tangential velocity vsky (in the plane of the sky) under the assumption that T1 is moving at a constant velocity, the radiative age coincides with the dynamical age (implying vsky ~ L/t), and there are not compression or expansion, re-energising, or bulk motions. We obtained a fitted velocity of vsky = 2100 ± 87 km s−1, which has to be considered as a lower limit because the value of B0 we used provides an upper limit to the radiative age and we ignored projection effects. The radial velocity (along the line of sight) of T1 is computed from the spectroscopic redshift of its host galaxy (zT1 = 0.0954) and that of A2142 (zA2142 = 0.0894) in units of the light speed c as (e.g. Davis & Scrimgeour 2014): vlos=c(zT1zA21421+zA2142)1650 km s1.$\[\mathrm{v}_{\mathrm{los}}=c\left(\frac{z_{\mathrm{T} 1}-z_{\mathrm{A} 2142}}{1+z_{\mathrm{A} 2142}}\right) \sim 1650 \mathrm{~km} \mathrm{~s}^{-1}.\]$(4)

The inferred velocity components provide constraints on the deprojected dynamics of T1, as we derive a 3D velocity of v=vlos2+vsky 2=2670 km s1$\[v=\sqrt{v_{\mathrm{los}}^2+v_{\text {sky }}^2}=2670 \mathrm{~km} \mathrm{~s}^{-1}\]$ and a viewing angle of i = arctan (vsky/vlos) = 52°. For comparison, the radial velocity dispersion of A2142 is σA2142 = 1193 km s−1 (Munari et al. 2014), meaning that v ~ 2.2σA2142, which is a reasonable value for head-tail galaxies.

thumbnail Fig. 13

Transitions along T1. Top: 143 MHz image (Fig. 3) stretched by a DEC to RA axis ratio of 2.5 to visually emphasise the transitions of each sub-region. Bottom: position (on the y-axis) of the emission peak along the tail (x-axis) with respect to the zero-level horizontal line shown in the inset (in black) and passing through the radio core.

5.3 Transitions along T1

In Sect. 5.2, we demonstrate that a simple TJP ageing scenario with a constant velocity is a good and physically reasonable representation of the global emission of T1. In line with this scenario, the absence of spectral flattening (Sect. 4.4) rules out ongoing large-scale re-acceleration processes; furthermore, our analysis in Sect. 4.6 does not conclusively support evidence of interaction with the ICM (Fig. 10). However, it is evident that local phenomena of unclear nature are taking place within its sub-regions. Indeed, we observed several major discontinuities in the surface brightness profile of T1 (Fig. 5), which we emphasise in Fig. 13. In this section, we briefly speculate on the possible origin of such discontinuities.

In T1-B, the tail shows a distinct set of wiggles, but they disappear as soon as the tail fades abruptly at the beginning of T1-C. A possible explanation to the wiggles is the precession of the jets before their downstream bending, but this scenario is unlikely. Indeed, we did not detect similar wiggles in T1-A. Thus, such a scenario would require the unlikely condition of sudden stabilising of the jets in a more recent epoch than the age of the plasma in T1-B.

Plausibly, the wiggles result from the development of Kelvin-Helmholtz (KH) instabilities, which are caused by the velocity difference across the interface between the jets and surrounding medium. These may be triggered either in T1-A, which is unresolved in our images, or at the transition in jet properties at the beginning of T1-B. With reference to the latter possibility, we recall the tantalising presence of a galaxy visible in the transition from T1-A to T1-B (Fig. 2). We stress that the growth rate of KH instabilities depends on the magnitude of the velocity gradient and on the effective viscosity and magnetic fields in the interface. Therefore, for the development of such instabilities, the jets have to move fast enough through the ICM, and the displacement of the jets downstream from the nucleus should result from the combination of the motion of the host galaxy through the ICM and the jet flow (O’Neill et al. 2019). As the jet decelerates and comes to rest in the ICM, the instabilities are not efficiently driven anymore, resulting in a turbulent tail which starts mixing with the ICM; this is a promising scenario for the plasma in T1-C and T1-D. Possibly, some local re-energising processes are occurring in T1-D, where the brightness increases and slight departures from the constant velocity scenario are found (Fig. 12). Higher resolution images of T1-A could further shed light on the onset of the instabilities.

5.4 Nuclear emission of T2

In this section, we discuss the properties of the host and core of T2 at different wavelengths. These are important to probe the overall origin of the source, as outlined in Sec. 5.5.

The host of T2 is an elliptical galaxy with no optical emission lines (Sec. 2). It emits in the soft X-ray band, but the analysis of its spectrum (see details in Appendix B) indicates that the observed emission comes from thermal gas. Therefore, the optical and X-ray spectra indicate that the accretion of the black hole is either currently off-state or radiatively inefficient.

Useful information is provided by the comparison of the core prominence at radio frequencies, which is CP=S1400T2D/S150T2=2×104$\[\mathrm{CP}=S_{1400}^{\mathrm{T} 2-\mathrm{D}} / S_{150}^{\mathrm{T} 2}=2 \times 10^{-4}\]$, when considering the ratio of 1.4 GHz flux density of the core and the 150 MHz total flux density, and the total 1.4 GHz radio power in logarithmic scale log10 P1400 = 23.9 W Hz−1. These values are particularly low (see e.g. Giovannini et al. 1988; de Ruiter et al. 1990; Jurlin et al. 2021, for works using these quantities or their combinations) and suggest that T2 is a remnant radio galaxy. Similar values were reported by Riseley et al. (2022a) for a remnant radio galaxy with a tailed morphology, whereas, for a comparison, we obtained CP = 4 × 10−2 and log10 P1400 = 24.4 W Hz−1 for T1 in line with active FRI galaxies.

5.5 Origin of T2

As highlighted in the previous sections, the global morphology of T2 is consistent with that of a HT galaxy, but it also exhibits unusual features. Moreover, its core prominence (Sec. 5.4) suggests that T2 is a remnant tailed galaxy. During its lifetime, the AGN may have experienced either a single or multiple outbursts, which we invoke to discuss possible scenarios explaining the radiatively old components of T2 (Sects. 4.3, 4.4) and some peculiar features.

In the context of a single AGN outburst, the core launched two radio jets in opposite directions, which were then bent by the ram pressure, and originated a tail as in classical HT galaxies. Afterwards, regions of the tail at different distances from the core passed through diverse ICM phases during the infalling towards the cluster centre. Such gas phases reshaped the subregions of T2 into the present T2-A and T2-B due to density and pressure gradients. This scenario is supported by the observed radial steepening of the spectral index (Fig. 8), in line with standard HT galaxies, but it relies on specific conditions of the ICM. Indeed, the relativistic plasma should be tremendously compressed at the location of the choke to explain the abrupt separation between T2-A and T2-B. Even though our projected thermodynamic maps (Fig. 10) do not provide solid conclusions owing to a combination of poor resolution and low ICM counts in the cluster outskirts, the existence of a thin layer where the thermal pressure is dramatically enhanced is disfavoured and appears unlikely.

A multiple AGN outburst scenario is more plausible and can be reconciled with the choke without invoking thermal pressure exerted by particular ICM layers. We can assume that a first AGN outburst was triggered when the host of T2 was far from its present position (beyond T2-B in projection). During the infall, the radio galaxy developed a tail, which we observe today as T2-B, and the core switched off at the location of the present choke. In proximity of its present location, the galaxy experienced a second AGN outburst, which would be responsible for the formation of T2-A. This scenario naturally explains the double-peaked surface brightness profile of T2 (Fig. 5) and the choke and is in line with the complex spectral distribution in the RCCDs that deviates from a single injection event6. As a first approximation, the integrated spectra of T2-A and T2-B (Fig. 6) suggest similar break frequencies that yield consistent radiative ages for a uniform magnetic field. This might be indicative of a short period between the two phases, possibly triggered by ram pressure itself (e.g. Poggianti et al. 2021; Peluso et al. 2022), but firm conclusions cannot be drawn (see e.g. Rudnick 2002 for a discussion on complex spectra).

In the context of the two scenarios described above, T2-A may represent the superposition of the two radio lobes in projection. The remnant scenario indicated by the core prominence analysis is supported by the non-detection of radio jets, which suggests that they are switched off on large scales (at least down to ~4 kpc). However, a completely different interpretation is also viable. Indeed, the conical morphology of the light bulb and the sharp transition at its edge (Fig. 2) are reminiscent of the structure of FRI-type jets reported in Laing & Bridle 2014 (see e.g. the case of M84 in their Fig. 3), which results from their deceleration on kiloparsec-scales. In other words, T2-A itself could potentially be a radio jet caught in its initial deceleration phase throughout the ambient medium. We notice that the sharp edge of T2-A suggests that the source approximately lies in the plane of the sky, implying that it has a one-sided jet. This is unlikely, therefore a possibility is that a combination of deceleration of the two jets and their downstream bending due to ram pressure during the infall is ongoing. In summary, the nature of T2-A remains unconfirmed, and higher-resolution radio data towards the nuclear regions could be helpful for further investigation.

Regardless of a single or double AGN outburst, the plume is interpreted as the oldest part of the tail. The aged relativistic plasma progressively expanded, diffused, mixed with the thermal ICM, and originated T2-C. Our high-resolution spectral index maps (Fig. 7) show flatter (but still ultra-steep) patches that deviate from the expected radial steepening of the tail. The morphology and spectral behaviour of T2-C may be indicative of some kind of interplay between thermal and non-thermal components and/or trace substructures of the magnetic field.

6 Summary and conclusions

In this work we reported on the study of two tailed radio galaxies, T1 and T2, in the galaxy cluster A2142. These targets show interesting morphological features that are suggestive of complex dynamics and interaction with the thermal ICM. By means of LOFAR, uGMRT, VLA, and MeerKAT radio data, we provided a detailed spectral analysis of T1 and T2. Auxiliary Chandra X-ray observations were used to investigate the local conditions of the ICM in the direction of the targets. In this section, we summarise our results and discuss future prospects.

T1 (Fig. 3) is a long HT galaxy extending for ~700 kpc and exhibiting clear surface brightness fluctuations and discontinuities that define four sub-regions (Figs. 2 and 13). A single power law of a slope of α = 0.87 ± 0.01 can fit the flux density measurements of T1 from 50 to 1810 MHz, but spectral breaks are found within its sub-regions (Fig. 6). The overall spectral index profile steepens with the increasing distance from the core (Fig. 8). In the low-frequency (50–143–323 MHz) regime that we consider here, standard ageing models (JP, TJP) can reproduce the observed spectral behaviour of T1 (Fig. 9) well. Under simple assumptions on the magnetic field, we produced a radiative age map (Fig. 12), which we used to constrain the tangential velocity of the target (Sec. 5.2). We computed a lower limit on the 3D velocity of the galaxy of v > 2670 km s−1, which is a factor of ~2 higher than the radial velocity dispersion of the cluster. Although we showed that a pure radiative ageing scenario and a constant velocity are reasonable approximations, each subregion shows signs of a complex phenomenology. Indeed, we detected spots, labelled as wiggles due to their oscillating pattern, where the surface brightness is locally enhanced (Fig. 5), sharp transitions along the tail, and fossil emission visible only at the lowest frequencies. All these features might be connected with the development and evolution of KH instabilities along the tail at different distances from the core. Moreover, we found that the spectral shape at high frequencies is surprisingly inconsistent with that at low frequencies (Fig. 4.4), hinting at subtle calibration effects, mixing of older and younger emitting electrons, or peculiar physical conditions (Sec. 5.1).

T2 (Fig. 4) is a HT galaxy extending for ~400 kpc and featuring morphological properties that are unusual for tailed sources. The whole structure of T2 can be decomposed into three sub-regions (Fig. 2) identified by distinct peaks in the surface brightness profile (Fig. 5). The radio core is not coincident with one of these peaks and our spectral analysis suggests that T2 is a remnant radio galaxy. While we did observe a spectral steepening along the tail that is in line with classical HT galaxies (Fig. 8), none of the standard (single injection) ageing models are able to reproduce the spectral measurements (Fig. 9), and non-trivial ad-hoc modelling would be required. The most intriguing feature of T2 is the choke, which is a sharp depletion of radio emission that separates the sub-regions T2-A and T2-B. We have discussed possible scenarios to explain such feature, which involve either a single or double AGN outburst event during the infall towards the cluster centre (Sec. 5.5). The simplest scenario that we propose is based on the assumption that two AGN outbursts are behind the tailed morphology of T2-B and the light bulb structure of T2-A, respectively. The X-ray analysis of the ICM disfavours exceptional thermal compression at the location of the choke (Sec. 4.6), but deeper data in this direction are necessary to definitely confirm our findings.

In conclusion, our work further highlights the increasing complexity in the phenomenology of tailed galaxies that has been shown by recent studies with the advent of deep and high-fidelity images from low (~100 MHz) to high (~1 GHz) frequencies. These works are providing valuable constraints on physical parameters, but theoretical modelling and numerical simulations are also necessary to shed light on the complex mechanisms that shape the properties of radio galaxies in clusters, such as re-energising processes, hydrodynamical instabilities, plasma mixing, magnetic field structure, and multiple AGN outbursts. As a follow-up work, we suggest looking towards the exploitation of LOFAR HBA data acquired with the international stations (available in the archive as part of the observations, with the core and remote stations) analysed here. These sensitive data could provide an insightful view of the sub-regions of T1 and T2 down to spatial scales of ~500 pc, possibly providing information on the onset of instabilities along T1 and the presence of jets and the choke in T2. Additionally, polarisation studies at gigahertz frequencies with MeerKAT are ongoing and will be used to constrain the structure of the magnetic field across T1 and T2.

Data availability

A copy of the reduced images is available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/690/A329

Acknowledgements

We thank the referee for their comments and suggestions. M.B. acknowledges support from the agreement ASI-INAF n. 2017-14-H.O and from the PRIN MIUR 2017PH3WAT ‘Blackout’. A.I. acknowledges the European Research Council (ERC) programme (grant agreement No. 833824, PI: B. Poggianti), and the INAF founding program ‘Ricerca Fondamentale 2022’ (project ‘Exploring the physics of ram pressure stripping in galaxy clusters with Chandra and LOFAR’, PI: A. Ignesti). C.J.R. acknowledges financial support from the ERC Starting Grant ‘DRANOEL’, number 714245. M.R., F.G., and G.B. acknowledge support from INAF mainstream project ‘Galaxy Clusters Science with LOFAR’. A.B. acknowledges financial support from the European Union – Next Generation EU. R.J.vW. acknowledges support from the ERC Starting Grant ClusterWeb 804208. D.V.L. acknowledges support of the Department of Atomic Energy, Government of India, under project No. 12-R&D-TFR-5.02-0700. The research leading to these results has received funding from the European Unions Horizon 2020 Programme under the AHEAD project (grant agreement No. 654215). LOFAR (van Haarlem et al. 2013) is the Low Frequency Array designed and constructed by ASTRON. It has observing, data processing, and data storage facilities in several countries, which are owned by various parties (each with their own funding sources), and that are collectively operated by the ILT foundation under a joint scientific policy. The ILT resources have benefited from the following recent major funding sources: CNRS-INSU, Observatoire de Paris and Université d’Orléans, France; BMBF, MIWF- NRW, MPG, Germany; Science Foundation Ireland (SFI), Department of Business, Enterprise and Innovation (DBEI), Ireland; NWO, The Netherlands; The Science and Technology Facilities Council, UK; Ministry of Science and Higher Education, Poland; The Istituto Nazionale di Astrofisica (INAF), Italy. This research made use of the Dutch national e-infrastructure with support of the SURF Cooperative (e-infra 180169) and the LOFAR e-infra group. The Jülich LOFAR Long Term Archive and the German LOFAR network are both coordinated and operated by the Jülich Supercomputing Centre (JSC), and computing resources on the supercomputer JUWELS at JSC were provided by the Gauss Centre for Supercomputing e.V. (grant CHTB00) through the John von Neumann Institute for Computing (NIC). This research made use of the University of Hertfordshire high-performance computing facility and the LOFAR-UK computing facility located at the University of Hertfordshire and supported by STFC [ST/P000096/1], and of the Italian LOFAR IT computing infrastructure supported and operated by INAF, and by the Physics Department of Turin University (under an agreement with Consorzio Interuniversitario per la Fisica Spaziale) at the C3S Supercomputing Centre, Italy. This research made use of the HOTCAT cluster (Bertocco et al. 2020; Taffoni et al. 2020) at Osservatorio Astronomico di Trieste. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. We thank the staff of the GMRT that made these observations possible. GMRT is run by the National Centre for Radio Astrophysics of the Tata Institute of Fundamental Research. The MeerKAT telescope is operated by the South African Radio Astronomy Observatory, which is a facility of the National Research Foundation, an agency of the Department of Science and Innovation. The scientific results reported in this article are based on observations made by the Chandra X-ray Observatory data obtained from the Chandra Data Archive. This research has made use of SAOImageDS9, developed by Smithsonian Astrophysical Observatory (Joye & Mandel 2003). This research has made use of the VizieR catalogue access tool, CDS, Strasbourg Astronomical Observatory, France (DOI: 10.26093/cds/vizier). This research made use of APLpy, an open-source plotting package for Python (Robitaille & Bressert 2012), Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013, Astropy Collaboration 2018), Matplotlib (Hunter 2007), Numpy (Harris et al. 2020).

Appendix A Spectral index error maps

In Fig. A.1, we report the error maps associated with the spectral index maps shown in Fig. 7. Errors are computed from the standard error propagation formula as Δα=|1ln(ν1ν2)|(ΔS1S1)2+(ΔS2S2)2,$\[\Delta \alpha=\left|\frac{1}{\ln \left(\frac{\nu_1}{\nu_2}\right)}\right| \sqrt{\left(\frac{\Delta S_1}{S_1}\right)^2+\left(\frac{\Delta S_2}{S_2}\right)^2},\]$(A.1)

where ΔS is obtained as in Eq. 1.

thumbnail Fig. A.1

Spectral index error maps of T1 and T2, corresponding to the maps shown in Fig 7.

Appendix B X-ray spectrum of T2 host

The host of T2 emits in the X-ray band (see Fig. 1). In this section we analyse its Chandra X-ray spectrum, which is shown in Fig. B.1.

We extracted the spectrum of the background from a circular annulus of width 25″ centred on the target. The spectrum of the source was extracted from a circular region of radius 5″ (8.5 kpc), chosen as the one that maximises the S/N, as: S/N=NcntfareaNbkgNcnt+fareaNbkg,$\[\mathrm{S} / \mathrm{N}=\frac{N_{\mathrm{cnt}}-f_{\text {area}} N_{\mathrm{bkg}}}{\sqrt{N_{\mathrm{cnt}}+f_{\text {area}} N_{\mathrm{bkg}}}},\]$(B.1)

where Ncnt is the total (target plus background) count number within the circle, Nbkg is the background count number within the annulus, and farea is the ratio of target to background extraction region areas. We extracted the corresponding background and source spectra from each pointing covering the region of T2, and we jointly fitted the background-subtracted spectra with an absorbed thermal component (phabs × apec). The Galactic hydrogen column density, redshift, and gas metallicity were kept fixed to the values considered for A2142 as in Sec. 4.6. The thermal component is sufficient to reproduce the spectrum of the target, as no hard X-ray (> 2 keV) emission, which is distinctive of AGN, is detected. An additional power law component (phabs × (apec+po)) to model the possible soft X-ray emission from AGN did not actually improve the fit (significance below 1σ) and was rejected.

The best-fit temperature, density (computed from the apec normalisation) are kT = 1.0 ± 0.1 keV and ne = (1.7 ± 0.4) × 10−3 cm−3. These values are consistent with those of the hot ionised medium in elliptical galaxies (e.g. Mathews & Brighenti 2003; Sun et al. 2007) and suggest that the observed X-ray emission comes from thermal gas rather than AGN activity.

thumbnail Fig. B.1

Chandra X-ray spectrum of the host of T2. The upper inset shows the extraction regions used for the source (red circle) and background (with annulus). Spectra were extracted for all available observations and jointly fitted, but a single spectrum is shown for clearer inspection. Data points were fitted with an absorbed thermal component (red curve). The subtracted background is shown in blue. The ratio of the data to the model is reported in the lower inset.

References

  1. Arnaud, K. A. 1996, ASP Conf. Ser., 101, 17 [Google Scholar]
  2. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  4. Bartalucci, I., Mazzotta, P., Bourdin, H., & Vikhlinin, A. 2014, A&A, 566, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bertocco, S., Goz, D., Tornatore, L., et al. 2020, ASP Conf. Ser., 527, 303 [NASA ADS] [Google Scholar]
  6. Blanton, E. L., Gregg, M. D., Helfand, D. J., Becker, R. H., & White, R. L. 2000, ApJ, 531, 118 [NASA ADS] [CrossRef] [Google Scholar]
  7. Blanton, E. L., Gregg, M. D., Helfand, D. J., Becker, R. H., & White, R. L. 2003, AJ, 125, 1635 [CrossRef] [Google Scholar]
  8. Botteon, A., Shimwell, T. W., Bonafede, A., et al. 2018, MNRAS, 478, 885 [Google Scholar]
  9. Botteon, A., Giacintucci, S., Gastaldello, F., et al. 2021, A&A, 649, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Brunetti, G., & Jones, T. W. 2014, Int. J. Mod. Phys. D, 23, 1430007 [Google Scholar]
  11. Bruno, L., Rajpurohit, K., Brunetti, G., et al. 2021, A&A, 650, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bruno, L., Botteon, A., Shimwell, T., et al. 2023, A&A, 678, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Cash, W. 1979, ApJ, 228, 939 [Google Scholar]
  14. Cavagnolo, K. W., Donahue, M., Voit, G. M., & Sun, M. 2009, ApJS, 182, 12 [Google Scholar]
  15. Chandra, P., Ray, A., & Bhatnagar, S. 2004, ApJ, 612, 974 [Google Scholar]
  16. Colla, G., Fanti, C., Fanti, R., et al. 1972, A&AS, 7, 1 [NASA ADS] [Google Scholar]
  17. Cuciti, V., Brunetti, G., van Weeren, R., et al. 2018, A&A, 609, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Cuciti, V., Cassano, R., Brunetti, G., et al. 2021, A&A, 647, A51 [EDP Sciences] [Google Scholar]
  19. Davis, T. M., & Scrimgeour, M. I. 2014, MNRAS, 442, 1117 [NASA ADS] [CrossRef] [Google Scholar]
  20. de Gasperin, F., Intema, H. T., Shimwell, T. W., et al. 2017, Sci. Adv., 3, e1701634 [Google Scholar]
  21. de Gasperin, F., Edler, H. W., Williams, W. L., et al. 2023, A&A, 673, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. de Ruiter, H. R., Parma, P., Fanti, C., & Fanti, R. 1990, A&A, 227, 351 [NASA ADS] [Google Scholar]
  23. Eckert, D., Gaspari, M., Vazza, F., et al. 2017, ApJ, 843, L29 [Google Scholar]
  24. Eckert, D., Finoguenov, A., Ghirardini, V., et al. 2020, Open J. Astrophys., 3, 12 [NASA ADS] [CrossRef] [Google Scholar]
  25. Edler, H. W., de Gasperin, F., Brunetti, G., et al. 2022, A&A, 666, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Einasto, M., Gramann, M., Saar, E., et al. 2015, A&A, 580, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Einasto, M., Deshev, B., Lietzen, H., et al. 2018, A&A, 610, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31P [Google Scholar]
  29. Flewelling, H. A., Magnier, E. A., Chambers, K. C., et al. 2020, ApJS, 251, 7 [Google Scholar]
  30. Gastaldello, F., Buote, D. A., Humphrey, P. J., et al. 2007, ApJ, 669, 158 [NASA ADS] [CrossRef] [Google Scholar]
  31. Giacintucci, S., & Venturi, T. 2009, A&A, 505, 55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Giovannini, G., Feretti, L., Gregorini, L., & Parma, P. 1988, A&A, 199, 73 [NASA ADS] [Google Scholar]
  33. Govoni, F., Dolag, K., Murgia, M., et al. 2010, A&A, 522, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Gramann, M., Einasto, M., Heinämäki, P., et al. 2015, A&A, 581, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Gunn, J. E., & Gott, J. Richard, I. 1972, ApJ, 176, 1 [NASA ADS] [CrossRef] [Google Scholar]
  36. Hardcastle, M. J. 2013, MNRAS, 433, 3364 [Google Scholar]
  37. Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [Google Scholar]
  38. Harwood, J. J., Hardcastle, M. J., Croston, J. H., & Goodger, J. L. 2013, MNRAS, 435, 3353 [NASA ADS] [CrossRef] [Google Scholar]
  39. Harwood, J. J., Hardcastle, M. J., & Croston, J. H. 2015, MNRAS, 454, 3403 [NASA ADS] [CrossRef] [Google Scholar]
  40. Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [Google Scholar]
  41. Ignesti, A., Brunetti, G., Shimwell, T., et al. 2022, A&A, 659, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Intema, H. T., van der Tol, S., Cotton, W. D., et al. 2009, A&A, 501, 1185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Jaffe, W. J., & Perola, G. C. 1973, A&A, 26, 423 [NASA ADS] [Google Scholar]
  44. Jones, P. A., & McAdam, W. B. 1994, PASA, 11, 74 [NASA ADS] [CrossRef] [Google Scholar]
  45. Joye, W. A., & Mandel, E. 2003, ASP Conf. Ser., 295, 489 [Google Scholar]
  46. Jurlin, N., Brienza, M., Morganti, R., et al. 2021, A&A, 653, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Kardashev, N. S. 1962, Sov. Astron., 6, 317 [NASA ADS] [Google Scholar]
  48. Katz-Stone, D. M., Rudnick, L., & Anderson, M. C. 1993, ApJ, 407, 549 [NASA ADS] [CrossRef] [Google Scholar]
  49. Komissarov, S. S., & Gubanov, A. G. 1994, A&A, 285, 27 [NASA ADS] [Google Scholar]
  50. Laing, R. A., & Bridle, A. H. 2014, MNRAS, 437, 3405 [NASA ADS] [CrossRef] [Google Scholar]
  51. Lal, D. V. 2020, AJ, 160, 161 [Google Scholar]
  52. Liu, A., Yu, H., Diaferio, A., et al. 2018, ApJ, 863, 102 [NASA ADS] [CrossRef] [Google Scholar]
  53. Lusetti, G., de Gasperin, F., Cuciti, V., et al. 2024, MNRAS, 528, 141 [NASA ADS] [CrossRef] [Google Scholar]
  54. Markevitch, M., Ponman, T. J., Nulsen, P. E. J., et al. 2000, ApJ, 541, 542 [Google Scholar]
  55. Mathews, W. G., & Brighenti, F. 2003, ARA&A, 41, 191 [NASA ADS] [CrossRef] [Google Scholar]
  56. McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, ASP Conf. Ser., 376, 127 [Google Scholar]
  57. Miley, G. K., Perola, G. C., van der Kruit, P. C., & van der Laan, H. 1972, Nature, 237, 269 [NASA ADS] [CrossRef] [Google Scholar]
  58. Missaglia, V., Massaro, F., Capetti, A., et al. 2019, A&A, 626, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Müller, A., Pfrommer, C., Ignesti, A., et al. 2021, MNRAS, 508, 5326 [CrossRef] [Google Scholar]
  60. Munari, E., Biviano, A., & Mamon, G. A. 2014, A&A, 566, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. O’Dea, C. P., & Baum, S. A. 2023, Galaxies, 11, 67 [CrossRef] [Google Scholar]
  62. Offringa, A. R., & Smirnov, O. 2017, MNRAS, 471, 301 [Google Scholar]
  63. Offringa, A. R., McKinley, B., Hurley-Walker, N., et al. 2014, MNRAS, 444, 606 [Google Scholar]
  64. O’Neill, B. J., Jones, T. W., Nolting, C., & Mendygral, P. J. 2019, ApJ, 884, 12 [CrossRef] [Google Scholar]
  65. Owen, F. N., & Rudnick, L. 1976, ApJ, 205, L1 [NASA ADS] [CrossRef] [Google Scholar]
  66. Owen, F. N., Burns, J. O., & Rudnick, L. 1978, ApJ, 226, L119 [NASA ADS] [CrossRef] [Google Scholar]
  67. Owen, F. N., Rudnick, L., Eilek, J., et al. 2014, ApJ, 794, 24 [Google Scholar]
  68. Owers, M. S., Nulsen, P. E. J., & Couch, W. J. 2011, ApJ, 741, 122 [Google Scholar]
  69. Pacholczyk, A. G. 1970, Radio Astrophysics. Nonthermal Processes in Galactic and Extragalactic Sources (USA: W.H. Freeman & Co Ltd) [Google Scholar]
  70. Pal, S., & Kumari, S. 2023, J. Astrophys. Astron., 44, 17 [Google Scholar]
  71. Peluso, G., Vulcani, B., Poggianti, B. M., et al. 2022, ApJ, 927, 130 [NASA ADS] [CrossRef] [Google Scholar]
  72. Perley, R. A., & Butler, B. J. 2013, ApJS, 204, 19 [Google Scholar]
  73. Pfrommer, C., & Jones, T. W. 2011, ApJ, 730, 22 [NASA ADS] [CrossRef] [Google Scholar]
  74. Pinkney, J., Rhee, G., Burns, J. O., et al. 1993, ApJ, 416, 36 [NASA ADS] [CrossRef] [Google Scholar]
  75. Planck Collaboration XXIV. 2016, A&A, 594, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Poggianti, B. M., Bellhouse, C., Deb, T., et al. 2021, IAU Symp., 359, 108 [NASA ADS] [Google Scholar]
  77. Riseley, C. J., Bonnassieux, E., Vernstrom, T., et al. 2022a, MNRAS, 515, 1871 [NASA ADS] [CrossRef] [Google Scholar]
  78. Riseley, C. J., Rajpurohit, K., Loi, F., et al. 2022b, MNRAS, 512, 4210 [NASA ADS] [CrossRef] [Google Scholar]
  79. Riseley, C. J., Biava, N., Lusetti, G., et al. 2023, MNRAS, 524, 6052 [NASA ADS] [CrossRef] [Google Scholar]
  80. Riseley, C. J., Bonafede, A., Bruno, L., et al. 2024, A&A, 686, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Robitaille, T., & Bressert, E. 2012, Astrophysics Source Code Library [record ascl:1208.017] [Google Scholar]
  82. Rossetti, M., Eckert, D., De Grandi, S., et al. 2013, A&A, 556, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Rudnick, L. 2002, New A Rev., 46, 95 [NASA ADS] [CrossRef] [Google Scholar]
  84. Rudnick, L., Cotton, W., Knowles, K., & Kolokythas, K. 2021, Galaxies, 9, 81 [NASA ADS] [CrossRef] [Google Scholar]
  85. Sanders, J. S. 2006, MNRAS, 371, 829 [Google Scholar]
  86. Sebastian, B., Lal, D. V., & Pramesh Rao, A. 2017, AJ, 154, 169 [Google Scholar]
  87. Shimwell, T. W., Hardcastle, M. J., Tasse, C., et al. 2022, A&A, 659, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Srivastava, S., & Singal, A. K. 2020, MNRAS, 493, 3811 [Google Scholar]
  89. Sun, M. 2009, ApJ, 704, 1586 [NASA ADS] [CrossRef] [Google Scholar]
  90. Sun, M., Jones, C., Forman, W., et al. 2007, ApJ, 657, 197 [NASA ADS] [CrossRef] [Google Scholar]
  91. Taffoni, G., Becciani, U., Garilli, B., et al. 2020, ASP Conf. Ser., 527, 307 [NASA ADS] [Google Scholar]
  92. Tchernin, C., Eckert, D., Ettori, S., et al. 2016, A&A, 595, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Terni de Gregory, B., Feretti, L., Giovannini, G., et al. 2017, A&A, 608, A58 [EDP Sciences] [Google Scholar]
  94. Tribble, P. C. 1993, MNRAS, 261, 57 [NASA ADS] [CrossRef] [Google Scholar]
  95. van Haarlem, M. P., Wise, M. W., Gunst, A. W., et al. 2013, A&A, 556, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. van Weeren, R. J., Ogrean, G. A., Jones, C., et al. 2017, ApJ, 835, 197 [Google Scholar]
  97. van Weeren, R. J., de Gasperin, F., Akamatsu, H., et al. 2019, Space Sci. Rev., 215, 16 [Google Scholar]
  98. Vazza, F., & Botteon, A. 2024, Galaxies, 12, 19 [NASA ADS] [CrossRef] [Google Scholar]
  99. Venturi, T., Rossetti, M., Brunetti, G., et al. 2017, A&A, 603, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Wang, Q. H. S., & Markevitch, M. 2018, ApJ, 868, 45 [CrossRef] [Google Scholar]
  101. Wilber, A., Brüggen, M., Bonafede, A., et al. 2018, MNRAS, 473, 3536 [Google Scholar]

1

M500 is the mass within a radius R500, which encloses 500ρc(z), where ρc(z) is the critical density of the Universe at a given redshift.

6

The complexity of the observed spectral distribution would require ad-hoc ageing modelling that we did not attempt in this work. A useful starting point might be the KGJP (Komissarov & Gubanov 1994) model, as it assumes a continuous injection of fresh electrons for a certain period, followed by a passive ageing.

All Tables

Table 1

Details of the radio data analysed in this work.

Table 2

Details of the Chandra X-ray data analysed in this work.

Table 3

Summary of the radio images shown in Fig. 3.

Table 4

Flux densities of T1, T2, and corresponding sub-regions.

Table 5

Integrated spectral index of T1 and T2.

Table 6

Parameters used to produce images for the spectral index maps in Fig. 7.

All Figures

thumbnail Fig. 1

Composite RGB image of A2142: radio (LOFAR at 143 MHz, 9″ × 6″) in red, optical (DSS-2, r-filter) in green, and X-rays (XMM-Newton, 0.7–1.2 keV) in blue. Labels indicate the tailed radio galaxies (T1, T2), the primary and secondary brightest cluster galaxies (BCG1, BCG2), and the two components (H1, H2) of the radio halo. The dashed white contours represent the 3σ level of the 73″ × 66″ LOFAR image at 143 MHz reported in Bruno et al. (2023), showing the third component (H3) of the radio halo. The yellow arrows indicate the location of the NW, central, and SE cold fronts reported by Rossetti et al. (2013).

In the text
thumbnail Fig. 2

Radio images and contours labelling the sub-regions and features discussed in the text. Left panels: 143 MHz images (see Figs. 3 and 4) of T1 (top) and T2 (bottom). Right panels: Pan-STARSS optical (composite i, r, g filters) images zoomed towards T1 (top) and T2 (bottom). Radio contours at 143 MHz (5σ, in red) and 1250 MHz (5σ for T1 and 3σ for T2, in cyan) are overlaid.

In the text
thumbnail Fig. 3

Radio images of T1. Details on resolution and noise are reported in Table 3. The overlaid contour levels are [±5, 10, 20,...] × σ. From top left to bottom right: LOFAR LBA at 50 MHz, LOFAR HBA at 143 MHz, GMRT at 234 MHz, GMRT at 323 MHz, uGMRT at 407 MHz, GMRT at 608 MHz, uGMRT at 1250 MHz (cropped to a smaller area for inspection purposes), MeerKAT at 1284 MHz, VLA at 1380 MHz, and VLA at 1810 MHz.

In the text
thumbnail Fig. 4

Radio images of T2. Details on resolution and noise are reported in Table 3. The overlaid contour levels are [±5, 10, 20, ...] × σ. From top left to bottom right: LOFAR LBA at 50 MHz, LOFAR HBA at 143 MHz, GMRT at 234 MHz, GMRT at 323 MHz, uGMRT at 407 MHz, GMRT at 608 MHz, uGMRT at 1250 MHz, MeerKAT at 1284 MHz, VLA at 1380 MHz, and VLA at 1810 MHz.

In the text
thumbnail Fig. 5

Radio surface brightness profiles of T1 and T2 normalised to the peak value within T1-A and T2-A, respectively. Left panels: profiles of T1 (middle) and T2 (bottom) at 143 MHz (resolution 9″ × 6″). The sampling boxes are shown in the top panels and have a width of 10″. Right panels: comparison of profiles at 50 (brown), 143 (yellow), 323 (green), 608 (light blue), and 1284 (magenta) MHz. The resolution of the images and the width of the sampling boxes (not shown) are 14″. In middle and bottom panels, the vertical lines indicate the boundaries of sub-regions as in Fig. 2.

In the text
thumbnail Fig. 6

Radio spectra of T1 and T2. Left and right panels refer to T1 and T2, respectively. Data points are the flux densities reported in Table 4, measured from the regions (including the green circle used for the core of T2) shown in the top panels following the same colour-code. Flux densities were fitted with single (solid lines) or double (solid plus dotted lines) power laws (see Table 5).

In the text
thumbnail Fig. 7

Spectral index maps of T1 and T2 computed with a threshold of S ≥ 5σ. In all the panels, the contour levels are [±5, 10, 20, ...] × σ of the lowest-frequency image in each combination. Top: spectral index maps at 30 between 50–143 MHz (left) and 143–323 MHz (right). Bottom: spectral index maps at 14″ between 143–323 MHz (left) and 608–1284 MHz (right). The corresponding error maps are shown in Fig. A.1.

In the text
thumbnail Fig. 8

Spectral index profiles as measured from maps in Fig. 7 for T1 (left) and T2 (right). The sampling boxes have beam-size width and are shown in the upper panels. The vertical lines indicate the boundaries of sub-regions as in Fig. 2.

In the text
thumbnail Fig. 9

Radio colour-colour diagrams of T1 (top) and T2 (bottom). Data points are colour-coded based on their increasing distance from the core from blueish to reddish colours. The black dashed line (one-to-one line) indicates a power law with α = αinj. The dotted lines represent the theoretical KP (orange), JP (blue), and TJP (red) ageing curves obtained with αinj = 0.51 (T1), αinj = 0.72 (T2), B0 = 2.2 μG (T1, T2) as input. The frequency pairs and resolution are the same as in Fig. 8.

In the text
thumbnail Fig. 10

Projected thermodynamic maps. Top left: temperature map from XMM-Newton photometric analysis covering the whole cluster (adapted from Rossetti et al. 2013). Top right: temperature and pressure maps from Chandra spectral analysis towards T1. Bin ‘4’ roughly marks the transition T1-A/T1-B. Middle: temperature and pressure maps from Chandra spectral analysis towards T2. Bin ‘1’ includes T2-A and the choke, bin ‘2’ includes T2-B. Bottom: temperature and pressure profiles for T1 (green) and T2 (magenta), as measured from the bins that are labelled in the corresponding maps with circles. The white dots in the Chandra maps indicate point sources excluded from the analysis. The LOFAR HBA radio contours are overlaid in black.

In the text
thumbnail Fig. 11

X-ray surface brightness analysis in the direction of T1. Left: Chandra X-ray image in the 0.5–2 keV band with overlaid radio contours. The surface brightness was measured from the red sector in annuli of width 2″. The fitted position of the discontinuity is shown within the sector, and is coincident with the NW cold front. Right: data points fitted with a broken power law (blue line). The red arrow indicates the peak discussed in Sec. 4.6.

In the text
thumbnail Fig. 12

Radiative age of T1. Top: radiative age map at 30″ obtained by fitting radio images at 50, 143, and 323 MHz with the TJP model (αinj = 0.51, B0 = 2.2 μG). Typical errors are ~20 Myr. The sub-regions are indicated by horizontal bars. Bottom: age profile measured in boxes of beam-size width. The boundaries of the sub-regions are indicated with dashed vertical lines. The blue solid line is the result of a linear fit to the data points.

In the text
thumbnail Fig. 13

Transitions along T1. Top: 143 MHz image (Fig. 3) stretched by a DEC to RA axis ratio of 2.5 to visually emphasise the transitions of each sub-region. Bottom: position (on the y-axis) of the emission peak along the tail (x-axis) with respect to the zero-level horizontal line shown in the inset (in black) and passing through the radio core.

In the text
thumbnail Fig. A.1

Spectral index error maps of T1 and T2, corresponding to the maps shown in Fig 7.

In the text
thumbnail Fig. B.1

Chandra X-ray spectrum of the host of T2. The upper inset shows the extraction regions used for the source (red circle) and background (with annulus). Spectra were extracted for all available observations and jointly fitted, but a single spectrum is shown for clearer inspection. Data points were fitted with an absorbed thermal component (red curve). The subtracted background is shown in blue. The ratio of the data to the model is reported in the lower inset.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.