Issue |
A&A
Volume 664, August 2022
|
|
---|---|---|
Article Number | A25 | |
Number of page(s) | 16 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202243174 | |
Published online | 03 August 2022 |
Seeing the forest and the trees: A radio investigation of the ULIRG Mrk 273⋆
1
Kapteyn Astronomical Institute, University of Groningen, Postbus 800, 9700 AV Groningen, The Netherlands
e-mail: kukreti@astro.rug.nl
2
ASTRON, The Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD, Dwingeloo, The Netherlands
3
INAF – Istituto di Radioastronomia, Via P. Gobetti 101, 40129 Bologna, Italy
4
Department of Physics and Astronomy, University of Sheffield, Sheffield S3 7RH, UK
5
Centre for Extragalactic Astronomy, Department of Physics, Durham University, Durham DH1 3LE, UK
6
Institute for Computational Cosmology, Department of Physics, University of Durham, South Road, Durham DH1 3LE, UK
7
Astronomisches Institut der Ruhr-Universität Bochum (AIRUB), Universitätsstrasse 150, 44780 Bochum, Germany
8
Department of Astronomy, University of Cape Town, Private Bag X3, Rondebosch 7701, South Africa
9
Hamburger Sternwarte, Universitat Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany
10
Thüringer Landessternwarte, Sternwarte 5, 07778 Tautenburg, Germany
11
Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía s/n, 18008 Granada, Spain
12
Department of Electrical Engineering, Chalmers University of Technology, Gothenburg, Sweden
13
Anton Pannekoek Institute, University of Amsterdam, Postbus 94249, 1090 GE Amsterdam, The Netherlands
14
Netherlands eScience Center, Science Park 140, 1098 XG Amsterdam, The Netherlands
15
Leiden Observatory, Leiden University, PO Box 9513 2300 RA Leiden, The Netherlands
16
University of Oslo Center for Information Technology, PO Box 1059 0316 Oslo, Norway
17
Astro Space Center of Lebedev Physical Institute, Profsoyuznaya Str. 84/32, 117997 Moscow, Russia
Received:
21
January
2022
Accepted:
4
May
2022
Galaxy mergers have been observed to trigger nuclear activity by feeding gas to the central supermassive black hole. One such class of objects are Ultra Luminous InfraRed Galaxies (ULIRGs), which are mostly late stage major mergers of gas-rich galaxies. Recently, large-scale (∼100 kpc) radio continuum emission has been detected in a select number of ULIRGs, all of which also harbour powerful Active Galactic Nuclei (AGN). This hints at the presence of large-scale radio emission being evidence for nuclear activity. Exploring the origin of this radio emission and its link to nuclear activity requires high sensitivity multi-frequency data. We present such an analysis of the ULIRG Mrk 273. Using the International LOFAR telescope (ILT), we detected spectacular large-scale arcs in this system. This detection includes, for the first time, a giant ∼190 kpc arc in the north. We propose these arcs are fuelled by a low power radio AGN triggered by the merger. We also identified a bright ∼45 kpc radio ridge, which is likely related to the ionised gas nebula in that region. We combined this with high sensitivity data from APERture Tile In Focus (Apertif) and archival data from the Very Large Array (VLA) to explore the spectral properties. The ILT simultaneously allowed us to probe the nucleus at a resolution of ∼0.3″, where we detected three components, and, for the first time, diffuse emission around these components. Combining this with archival high frequency VLA images of the nucleus allowed us to detect absorption in one component, and a steep spectrum radio AGN in another. We then extrapolate from this case study to the importance of investigating the presence of radio emission in more ULIRGs and what it can tell us about the link between mergers and the presence of radio activity.
Key words: galaxies: active / radio continuum: galaxies / galaxies: individual: Mrk 273 / galaxies: jets / galaxies: interactions / techniques: high angular resolution
Fits files of the radio maps are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/664/A25
© P. Kukreti et al. 2022
Open 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
Galaxy mergers are spectacular events that are considered to be the way massive galaxies form. They are also thought to be one of the mechanisms that trigger AGNs, including powerful radio AGNs, by helping the gas reach the nuclear regions and accrete onto the central supermassive black hole (SMBH; e.g., Hopkins et al. 2008; Yuan et al. 2010). Galaxy mergers go through different phases, which can show intense star formation or nuclear activity, and they can be identified using different diagnostics. A class of objects related to a merger phase are UltraLuminous InfraRed Galaxies (ULIRGs; e.g., Sanders & Mirabel 1996; Rigopoulou et al. 1999; De Lucia et al. 2006), which are understood to be formed by major mergers of gas-rich spiral galaxies. These are extremely bright IR sources with typical luminosities of LIR ≥ 1012 L⊙ (for a review see Lonsdale et al. 2006). ULIRGs are known to be powered by both starbursts and AGNs (e.g., Veilleux et al. 2009; Nagar et al. 2003; Nardini et al. 2010), and they can therefore represent different phases in the evolutionary sequence of a major merger. Although the properties of optical and radio AGNs in mergers and their feedback effects have been investigated by many studies (e.g., Silk & Rees 1998; Ramos Almeida et al. 2011; Cano-Díaz et al. 2012; Cresci et al. 2015; Carniani et al. 2016; Murphy 2013; Barcos-Muñoz et al. 2017), the presence and properties of large-scale radio emission in such systems are less known.
The AGN phase in mergers is understood to occur later than the starbursts. This is due to the fact that in order to be accreted to the central SMBH, the gas has to reach the central few parsecs and lose the majority of its angular momentum (Shlosman et al. 1990; Di Matteo et al. 2005; Hopkins et al. 2006), whereas the starburst can precede on galactic scales (Mihos & Hernquist 1996; Springel et al. 2005). However, optical imaging studies have found that AGNs are also triggered prior to the coalescence of the galaxy nuclei, that is at the pre-coalescence stage (Ramos Almeida et al. 2011; Bessiere et al. 2012; Pierce et al. 2022). Pierce et al. (in prep.) have found that this is true for the majority (61%) of nearby type 2 quasars. In general, phases of AGN activity are understood to peak either at the pre-coalescence stage or following the coalescence of the merging nuclei (Mihos & Hernquist 1996; Tadhunter et al. 2011).
Investigating the presence of radio AGNs is interesting because by using the radio spectrum, we can gain information about the timescales on which nuclear activity is triggered in mergers. This will help us understand the role mergers play in triggering radio galaxies. For instance, Emonts et al. (2006) studied the radio properties of B2 0648+27 and found a significant time delay between the merger-starburst event and the triggering of radio activity.
In their study of radio continuum emission from ULIRGs, Yun et al. (1999) (hereafter Y99) estimated that only ∼30% of the mergers host an AGN, on the basis of mid-IR excess. Out of these AGNs, ∼40% host a radio AGN. Among the merger systems with an AGN identified in mid-IR, Y99 observed large (hundreds of kiloparsecs) radio continuum tails in a small number of cases (including Mrk 273). They concluded that the luminosity, size, high degree of polarisation (20−50%), and the presence of a powerful AGN in these systems suggested that the power source for these tails were the AGNs and not the nuclear starburst. Similar radio features have been detected in more ULIRGs recently. Hayashi et al. (2021) discovered radio emission with a size of ∼100 kpc in IRAS 01004−2237. The two-sided morphology and lack of X-ray detection suggests that the radio emission was not powered by a cluster merger, but an AGN in that system. Nandi et al. (2021) also observed extended radio emission in IRAS 15001+1433, which hosts an AGN as well. They also conclude that ULIRGs harbour young (∼0.4−20 Myr) radio sources and are possible progenitors of radio galaxies. So these detections beg the question − are such large-scale tails signposts of radio AGN activity in ULIRGs?
To answer this question, detailed studies of the radio properties of ULIRGs are required. One such study investigated the ULIRG Mrk 231. Radio emission on scales of parsecs to tens of kiloparsecs was detected in this object (Baum et al. 1993; Ulvestad et al. 1999). Morganti et al. (2016) studied the large southern radio lobe of tens of kiloparsecs size and detected a plateau of emission near the nucleus and poorly collimated bridge structure a few kiloparsecs south of the nucleus (see their Fig. 6). They concluded that the radio lobe emission could either be due to fuelling of fresh electrons by the central AGN, or in situ acceleration of electrons, or both.
Although the current detection rate of large-scale radio emission from ULIRGs is low (≲10%), ongoing and upcoming high sensitivity surveys like the LOFAR LBA Sky Survey (LoLSS; de Gasperin et al. 2021), LOFAR Two-metre Sky Survey (LoTSS; Shimwell et al. 2017, 2019, 2022), APERture Tile In Focus (Apertif; van Cappellen et al. 2022) and Evolutionary Map of the Universe (EMU; Norris et al. 2011) provide a great opportunity to investigate its presence for a large sample.
In this paper, we present such a study of the merger/ULIRG Mrk 273 using LOFAR (down to 0.3″ using international stations) and Apertif. Since this is a well known system, archival data with the Very Large Array (VLA) were also available, which we used to explore the spectral properties. Mrk 273 is a good object for this study since it is the second nearest ULIRG (z = 0.0373, LIR = 1012.1 L⊙; Kim & Sanders 1998) and a late stage merger, with a starburst and multiple AGNs. Previous studies have found evidence for three nuclei in the system, two in the near-IR (Knapen et al. 1997) − northern (N) and south-western (SW) − and a third south-eastern (SE) in the radio (Condon et al. 1991) and [OIII] (Rodríguez Zaurín et al. 2014, hereafter RZ14). The presence of an AGN in these nuclei has been confirmed in the N component by Knapen et al. (1997) and Liu et al. (2019, hereafter L19), in the SW by Iwasawa et al. (2011) and L19, and a third AGN has been suggested in the SE component by RZ14.
The presence of a dense neutral and molecular gas disk (Condon et al. 1991; Cole et al. 1999; Carilli & Taylor 2000) and a collection of compact luminous supernovae remnants and radio supernovae (Bondi et al. 2005) have also been confirmed in the N component. Vardoulaki et al. (2015) had suggested the presence of a radio AGN in the nuclear region, on the basis of the steep spectral index (α ∼ 1.1) from 1400−8415 MHz, but they could not pin down its location to one of the three nuclear components. The nuclear region also has a high star formation rate (SFR) of ∼139 M⊙ yr−1 (Cicone et al. 2014).
The stellar population in the nuclear region and optical tail has an age of 0.7−2.0 Gyr, but there is a very young stellar population of age ≲50 Myr present as well (Rodríguez Zaurín et al. 2009). The gas kinematics in the nuclear region of Mrk 273 is complex, and emission line studies have revealed ionised gas outflows up to ∼5 kpc (Rupke & Veilleux 2013; RZ14; Spence et al. 2016, hereafter S16; Leung et al. 2021). On the kiloparsec scales, optical images of the system show a ∼30 kpc tidal tail, filaments, and diffuse emission to the south and north-east of the nucleus (RZ14, S16, also see Fig. 1).
![]() |
Fig. 1. Multi-wavelength emission from Mrk 273. The blue contours map the 144 MHz emission from the LOFAR 144 MHz image, red contours map the soft X-ray emission (0.4−2 keV) from L19, overlaid on the optical continuum image from S16 at λ = 6572 Å. The optical continuum tail seen towards the south of the nucleus is marked. The extended soft X-ray nebulae seen in south-west (SW) and north-east (NE) can be seen as well. The Chandra X-ray image has been smoothed with a Gaussian kernel of σ = 10 pixels, using the CIAO package. |
Although large-scale radio emission in Mrk 273 was reported by Y99, their low resolution (∼20″) did not allow them to resolve the different features of this emission. Using LOFAR, Apertif and VLA data, we have identified extended low surface brightness features in Mrk 273, and investigated their spectral properties. We have also studied the radio emission from the resolved nuclear components. The paper is structured in the following manner: In Sect. 2, we describe the data reduction and the procedure to make the spectral index maps; in Sect. 3 we present the morphology and spectral properties of the source, the spectral age estimates and the absorption analysis; and in Sect. 4, we discuss the origin of the radio emission and the spectral index trends of the system. Throughout the paper, we define the spectral index α as: S ∝ ν−α. The cosmology adopted in this work assumes a flat universe with H0 = 71 km s−1 Mpc−1, Ωm = 0.27 and Ωvac = 0.73. At the redshift of Mrk 273, 1″ corresponds to 0.73 kpc.
2. Data reduction
We have combined data from the International LOFAR Telescope (ILT) Low Band Antenna (LBA) at 42−68 MHz and High Band Antenna (HBA) at 120−168 MHz, Apertif at 1.36 GHz and VLA at 1.4 GHz and 4.7 GHz to trace the radio continuum and spectral properties of Mrk 273 at both high and low resolutions.
2.1. LOFAR 144 MHz
To trace the low frequency radio continuum emission, we have used the mosaic P205+55 from the LOFAR Two-metre Sky Survey (LoTSS; Shimwell et al. 2019). Mrk 273 lies at a distance of 1.2° from the phase centre. The dataset is part of LoTSS-DR2 (Shimwell et al. 2022) and the observations were carried out with the standard LoTSS survey setup, i.e. 8 h on-source observations bookended on either side with a flux density calibrator, in this case 3C 147. A total of 62 Dutch stations were used for this observation, which included 48 core stations with baselines up to 4 km and 14 remote stations with baselines up to 120 km. This dataset did not have all 13 international stations, which have baselines up to ∼2000 km. Therefore, for the sub-arcsecond image, we used data from Morabito et al. (2022). They used data from a re-observation of this field. This re-observation also used the standard LoTSS observational setup, and the same flux density calibrator but has 13 international stations.
Both datasets were recorded with an integration time of 1 s, a 48 MHz bandwidth centred at 144 MHz and a channel width of 3.05 kHz. The data were then passed through the standard LOFAR pre-processing pipeline (Heald et al. 2010) which performed the RFI flagging using the AOflagger (Offringa et al. 2010, 2012) and averaged down to a channel width of 12.2 kHz. Direction independent calibration was then performed using the PREFACTOR1 pipeline (de Gasperin et al. 2019; van Weeren et al. 2016; Williams et al. 2016).
2.1.1. Dutch array image
The direction dependent self-calibration of the low-resolution Dutch array image was done using the DDF-pipeline, described in Shimwell et al. (2019), Tasse et al. (2020). The LoTSS-DR2 data reduction is described in detail in Shimwell et al. (2022). We have used images of the field at both 6″ and 20″ resolution to trace as much of the faint large-scale emission as possible while also being able to map the details in the large-scale structure. The 6″ image has an rms noise of 80 μJy beam−1 and an integrated flux density of 466 ± 47 mJy measured using 3σ contours as reference and is shown in Fig. 2. The integrated flux density in the 20″ image is 451 ± 45 mJy. The error associated to the measured flux densities is dominated by the flux calibration uncertainty, that is typically assumed to be 10%. The residual image was noise like – i.e. the object was full deconvolved. To recover more of the extended structure in this object, we have also smoothed the 20″ image to a resolution of 30″, as shown in Fig. 2. We recover the same integrated flux density in the 6″, 20″ and 30″ image, within the errors. However, the structure of the northern arc is most clearly visible in the 30″ image.
![]() |
Fig. 2. LOFAR 144 MHz images showing emission at different scales in Mrk 273. Left panel: LoTSS-DR2 image at 30″ resolution shows the low surface brightness large-scale emission, with the northern arc marked. Black cross in the northern arc marks the nearby galaxy 2MFGC 11079 at that location. Middle panel: LoTSS-DR2 image at 6″ resolution, showing the southern arc and ridge in Mrk 273. Right panel: sub-arcsecond image of the nucleus of Mrk 273 at a resolution of 0.3″ × 0.2″. The northern (N), south-eastern (SE), south-western (SW) and diffuse north-east (NE) components are marked. The contour levels in all images are 3 |
2.1.2. Sub-arcsecond International array image
The calibration of the international stations data for this field was done using the LOFAR long baseline pipeline2. The pipeline and the data reduction for the LoTSS pointing P205+55 are described in detail using the same field in Morabito et al. (2022). The paper also described the procedure for the delay and self-calibration of the delay calibrators in the field (Jackson et al. 2016, 2021). After applying the self-calibration solutions from the delay calibrator to the target, we have used the procedure from van Weeren et al. (2021) to iteratively self-calibrate the phases and amplitudes and finally make the final image of the target. This procedure uses WSClean (Offringa et al. 2014; Offringa & Smirnov 2017). Using the Briggs weighting scheme and a robust = −1 along with a minimum uv length of 40kλ, we obtain a resolution of 0.2″ × 0.3″ and an rms noise of 90 μJy beam−1 for the target image, see Table 1. The image is shown in Fig. 2. To recover more extended emission around the nuclear components, we made another image with robust = 0 and a minimum uv length of 2kλ. This image has a resolution of 0.4″ × 0.5″, an rms noise of 78 μJy beam−1 and is shown in Fig. A.1e.
Image statistics summary.
2.2. LOFAR 54 MHz
We have used the image from the mosaic P205+55 of the LOFAR LBA Sky Survey (LoLSS; de Gasperin et al. 2021). The survey observations were carried out with the LOFAR Low Band Antenna (LBA) Dutch stations for ∼8 h and a 24 MHz bandwidth centred at 54 MHz. The time and frequency resolution were 1s and 3.052 kHz before Radio Frequency Interference (RFI) flagging, and were averaged to 2 s and 48.828 kHz after flagging. The data reduction is described in detail in the survey paper and in de Gasperin et al. (2020). The image resolution was 15″ with an rms noise of ∼1.6 mJy beam−1 and is shown in Fig. A.1a.
2.3. Apertif 1360 MHz
The APERture Tile In Focus (Apertif) phased-array feed (PAF) system is installed at the Westerbork Synthesis Radio Telescope (WSRT) and operates around 1400 MHz (see van Cappellen et al. 2022 for the full description). Because of the large field of view (the total area covered in one observation by such a mosaic is about 3.5° ×3°) combined with the sensitivity of the radio continuum images (about 30−40 μJy) and the spatial resolution of about 12″ × 12″/sin δ, Apertif is an ideal survey instrument and provides an ideal complement to the surveys done with LOFAR.
The Apertif imaging surveys started on 1 July, 2019, with the aim of covering up to ∼2500 square degrees of the northern sky. The first data from the Apertif surveys – providing radio continuum, HI, and polarisation images and cubes – have become publicly available in the first data-release, which was done in November 20203. The observations covering the position of Mrk 273 were taken on the 22 December 2020 (ObsID 201222016) with a total observing time of 11 h. 3C 196 was used as the flux calibrator. The data of the beams including Mrk 273 were flagged and cross calibrated by the automatic Apertif pipeline, Apercal (Adebahr et al. 2022). The self-calibration and the final image were made manually using MIRIAD (Sault et al. 1995), followed by the primary beam correction (Kutkin et al., in prep.4). The image is shown in Fig. A.1b and has an rms noise of 35 μJy beam−1 and a restoring beam of 12″ × 1″ (PA = 0°), see Table 1.
2.4. Archival data
2.4.1. Low resolution VLA 1.45 GHz and 4.7 GHz
We reprocessed archival VLA 1.45 GHz and 4.7 GHz observations of the target. 1.45 GHz observations were carried on 31 October 2010 with the VLA in C configuration (Project code: 10B-212). The target was observed for ∼50 min with an integration time of 1s and a 128 MHz bandwidth split into 64 channels of 2 MHz each. 3C 286 and J1349+5341 were used as the flux and phase calibrators respectively.
Mrk 273 was also observed at 4.7 GHz on 29 August 2010 with the VLA in D configuration (Project code: AL746). The target was observed for ∼6 min with an integration time of 1 s and a bandwidth of 1024 MHz. The entire bandwidth was split into 8 equal width sub-bands of 128 MHz, with a sub-band split into 64 channels of 2 MHz each. 3C 286 was used as the flux calibrator and J1339+6328 and J1349+5341 were used as the phase calibrators.
The data were manually flagged and cross-calibrated in Astronomical Image Processing System (AIPS, Greisen 2003). The flux scale was set according to Perley & Butler (2013). Flagging on the C band data was done using the task FLAGR in AIPS. This data was then self-calibrated and imaged using Common Astronomy Software Applications (CASA, McMullin et al. 2007). The final images were made using Briggs weighting with a robust = 0.5 for the 1.45 GHz and 4.7 GHz image. Self-calibration of the dataset was done initially with phase only self-calibration, followed by amplitude and phase self-calibration. The resolution obtained was 16.8″ × 13″ (PA = 15°) for the 1.45 GHz image and 18.6″ × 13.2″ (PA = 63°) for the 4.7 GHz image. The rms noise is 60 μJy beam−1 for the 1.45 GHz image and 28 μJy beam−1 for the 4.7 GHz image, see Fig. A.1 and Table 1. Primary beam correction was done for both images using the task impbcor in CASA.
2.4.2. High resolution VLA 4.7 GHz and 8.4 GHz
We have reprocessed VLA 4.7 GHz data to image the nucleus of Mrk 273 and match the spatial resolution of the LOFAR international stations data. Mrk 273 was observed with the VLA on 22 June 2011 in A configuration (Project code: AL746). The target was observed for ∼8 min with an integration time of 1 s and a bandwidth of 1024 MHz split into 8 equal width sub-bands. 3C 286 and J1349+5341 were used as the flux density calibrator and phase calibrator respectively. Similar to the low resolution VLA images, the data was flagged, cross-calibrated and bandpass calibrated in AIPS. This data was then self-calibrated and imaged in CASA with the Briggs weighting scheme and a robust = 0.5. For the self-calibration, phase only calibration was performed first, followed by amplitude and phase calibration. Primary beam correction was again done for the image using the task impbcor. The final resolution obtained was 0.45″ × 0.3″ and an rms noise of 30 μJy beam−1.
We have also used the 8.4 GHz image from Condon et al. (1991), made with the VLA in A configuration. The image has a resolution of 0.3″ × 0.3″ and an rms noise of 28 μJy beam−1. The images are shown in Fig. A.1 and image statistics are summarised in Table 1.
2.5. Flux density scale
An accurate flux density scale is necessary for studying the spectral properties of the emission. To check the accuracy of the flux scale of our low resolution images, we have plotted the integrated flux densities of Mrk 273 along with the values from NED. In our low resolution images, the integrated flux densities were measured using 3σrms contours as reference. We used a flux density scale error of 15% for LOFAR LBA, 10% for LOFAR HBA, and 5% for Apertif 1.36 GHz, VLA 1.45 GHz, 4.7 GHz and 8.4 GHz data. The plot is shown in Fig. 3 and shows that our flux scale is in agreement with the literature. Our 144 MHz flux density is lower, although in agreement within errors with the literature value at 150 MHz, which is taken from the 6C survey (Hales et al. 1990). This is likely due to the fact that 6C survey has a much lower resolution of ∼4.5′, and therefore collects all the low level flux also from the companions that we excluded.
![]() |
Fig. 3. Integrated flux densities of Mrk 273 from our low resolution images, measured using 3σ contours for reference and NED. The red points show measurements from this study and orange points mark the NED values. |
3. Results
We describe here the morphology of the large (∼300 kpc) and small (∼4 kpc) scale emission in Mrk 273, followed by the spectral properties and the estimation of the spectral age of the emission.
3.1. Morphology
Mrk 273 shows interesting morphology on both small and large-scales which we describe below.
3.1.1. Nuclear region
The radio morphology of the Mrk 273 nucleus has been observed before at GHz frequencies with sub-arcsecond resolution (Cole et al. 1999; Bondi et al. 2005) and multiple components were identified. As mentioned before, dual AGNs have been identified in the X-ray in the N and SW component, and another AGN has been proposed in the SE component. Our new 144 MHz sub-arcsecond image is shown in Fig. 2. Two bright components, N and SE, of ∼0.9 kpc size can be seen, separated by ∼0.7 kpc. Two bright knots of emission ∼0.6 kpc apart can also be seen for the first time in the SW component, along with some diffuse emission. We also observe diffuse emission of about ∼1 kpc and ∼0.4 kpc in the north-east (NE) and north-west (NW). In our long baseline image at 144 MHz with a lower minimum uv length and a 0.5″ resolution (Fig. A.1e), we also detect diffuse emission around the components. The interesting feature to note in the figure is the emission extending to ∼1.5 kpc southwards of the SE component, possibly linking it to the southern arc or ridge emission. Some more diffuse emission has been detected around the SW component. The N, SE and SW component had been detected before at 8.4 GHz by (Condon et al. 1991). However, at 144 MHz, the SE component is the brightest whereas the N component is brightest in the 8.4 GHz image. This suggests that the emission in the N component is absorbed at low frequencies. We measured the flux density of each component by fitting a 2D Gaussian over a region of the size of 3σrms level contours at every frequency, where σrms is the local rms noise in the image. The flux densities, luminosities and spectral indices of these components are listed in Table 2.
Flux densities and spectral indices for nuclear components.
3.1.2. Large-scale structure
Figure 2 shows the large-scale emission in Mrk 273 at 144 MHz. The 6″ map shows the southern arc with a remarkably smooth curvature. This arc has a surface brightness ranging from 0.15−0.30 mJy arcsec−2 and a total projected linear extent of ∼100 kpc (measured using 3σ contours as reference). The emission is slightly more diffuse at the end of the arc than near the nucleus. Although the giant radio continuum ‘plume’ had been reported before by Y99, this is the first time this large-scale structure has been resolved and other features identified. A new feature we observe is the large arc in the north, shown in the 30″ image in Fig. 2. This arc has a much lower surface brightness than its counterpart in the south, 0.04−0.08 mJy arcsec−2, and a larger projected linear extent of ∼190 kpc. We recover part of this emission in the Apertif image as well, however, we do not have enough sensitivity to recover this emission in the VLA 1455 MHz and 4700 MHz images.
Another new feature of interest is the bright ridge of emission, located at ∼25 kpc southwards from the nucleus, marked in Fig. 2. It is oriented almost perpendicular to the southern optical continuum tail, shown in Fig. 1, and is also coincident with the linear feature branching out of the optical tail to the east. It overlaps well with the ionised gas (RZ14,S16) and soft X-ray (Xia et al. 2001; L19) nebula in this region. The ridge has a surface brightness ranging from 0.5−0.7 mJy arcsec−2 and a projected linear extent of ∼45 kpc. No radio emission is seen to the north-east of the nucleus in the low resolution images, where extended soft X-ray emission can be seen in Fig. 1.
3.2. Spectral index maps
In order to trace the large-scale spectral properties in Mrk 273, we have constructed spectral index maps from 144−1455 MHz at ∼16″ and 20″, and 1455−4700 MHz at 20″. We do not construct spectral index maps of the nucleus with our sub-arcsecond images. This is due to the fact that the 1400 MHz image from Cole et al. (1999) is not publicly available in FITS format. This image is crucial to cover the frequency range from 144 MHz to 4700 MHz. Without the 1400 MHz image, making a spectral index map from 144–4700 MHz would not give an accurate picture of the spectral properties. For the low resolution spectral index maps, we have not used the 54 MHz image since it did not recover the extended emission as well as the higher frequency images. Spectral index maps require the images to recover emission on the same angular scales and have the same resolution. We have ensured that the data sets used for the low resolution images had short baselines to be sensitive to emission at an angular scale of ∼150″ which is the size of the southern arc (an interferometer with a shortest baseline of Dmin is sensitive to a largest angular scale of 0.6 λ/Dmin; Tamhane et al. 2015). To have the same resolution in both images for the spectral index map, the image with higher resolution was smoothed using a 2D Gaussian in the image plane with the task IMSMOOTH in CASA. To trace the spectrum of the northern arc, we also constructed a low resolution spectral index map from 144−1360 MHz after smoothing the Apertif 1360 MHz image to 20″ resolution.
Position offsets in images can be caused by phase calibration errors, and it is necessary to correct for such offsets and align the images before making the spectral index maps. For this purpose, we have fitted a 2D Gaussian over the nucleus of Mrk 273 in the images and derived the pixel position coordinates of the peak flux density. All images were then regridded and aligned using the position from one image as reference with the tasks IMHEAD and IMREGRID in CASA. The position offset after this procedure was ≤0.01 pixels, which is sufficient for our analysis. The spectral index maps were made with the task IMMATH in CASA. Throughout the paper, spectral index errors are calculated as
where ΔS1 and ΔS2 are errors in the flux densities which include statistical errors in the measurements as well as uncertainties in the overall scale.
To derive the spectral properties from 144−1400 MHz, we follow two approaches for the southern and northern arc regions. The analysis of the southern arc was done done using the spectral index image obtained at ∼16″. The main properties (that will be described in Sect. 3.3.2) were first derived using the Apertif image and then refined using the VLA image. The spectral index values of the southern arc and the trends in spectral properties with the two were found to be in agreement, and in the remainder of the paper the analysis will be done using the image obtained with the VLA (Fig. 4a) because it is less affected by calibration artefacts. For this map, emission within the 3σ contours was used.
![]() |
Fig. 4. Spectral index maps of large scale emission in Mrk 273. (a) Spectral index map from 144−1455 MHz made with the LOFAR and VLA images at 16.8″ × 13″ resolution. The local rms noise, σrms, in the 144 MHz image is 200 μJy beam−1 and in the 1455 MHz image is 60 μJy beam−1. The contours mark the 1455 MHz emission at levels (−3, 3, 5, 10, 15, 25)×σrms. The map shows a spectral gradient across the arc, with α ∼ 0.6 at the southern outer edge and ∼0.9 at the northern inner edge. The ridge has a steep spectral index of α ≈ 0.88 ± 0.07. (b) The 20″ spectral index map from 144−1360 MHz made with the LOFAR and Apertif images. The local rms noise, σrms, in the 144 MHz image is 100 μJy beam−1 and in the 1360 MHz image is 37 μJy beam−1. The contours mark the 1360 MHz emission at levels (−3, 3, 5, 10, 25, 100)×σrms′. The pixels outside the 3σ contours show the lower limit of the spectral index. The black cross marks the background galaxy present at the location of the northern arc. The inset maps show the spectral index errors. (a) Spectral index map |
In the northern arc region, the high sensitivity to low surface brightness emission of the Apertif image allows it to recover more emission than the VLA 1455 MHz image, and we use it to obtain the spectral indices in (at least part of) this region (see Sect. 3.3.2). We smoothed the Apertif map to match the LOFAR 20″ image and used > 3σ emission to construct a spectral index map from 144−1360 MHz (Fig. 4b). Since the LOFAR image recovered more emission in the northern arc than the Apertif image, we estimate a lower limit of the spectral index for pixels outside 3σ contours in the Apertif image. To do this, we have used a value of 3σ as the upper limit on the flux density for these pixels.
For the nucleus, we have used the integrated flux densities of the different components (see Fig. 2) from the sub-arcsecond images, to calculate the spectral indices. This is done because the high frequency sub-arcsecond images do not recover all the emission we see in the 144 MHz sub-arcsecond images.
3.3. Spectral properties
Our multi-frequency images allow us to probe the spectral properties of the radio emission from Mrk 273. Below, we first discuss the spectral properties of the nuclear components and then the large-scale spectral properties.
3.3.1. Nuclear region
To study the spectrum of the components in the nucleus of Mrk 273 (right panel in Fig. 2), we estimated their integrated flux densities and calculated their spectral indices. The flux densities for these components were measured from the sub-arcsecond images. For the 1400 MHz flux density, we used the measurements from Cole et al. (1999), who studied the 1400 MHz emission from the nuclear components with MERLIN at ∼0.2″ resolution. We assumed a flux scale error of 5% on their values.
The spectral indices of nuclear components are summarised in Table 2. In the N component, we find an inverted spectrum from 144−1400 MHz with = −0.26 ± 0.05, which steepens to
= 0.47 ± 0.06 and
= 0.60 ± 0.14. This tells us that the emission (at low frequencies) is absorbed. The SE component shows a spectral index of
= 0.75 ± 0.05, which steepens at higher frequencies to
= 1.14 ± 0.14. This steep spectral index suggests the presence of a radio AGN in this component. This would be in agreement with Vardoulaki et al. (2015), who constructed a low resolution spectral index map of the nuclear region and proposed the presence of a radio AGN in it. In the SW and NE components, we detect steep spectral indices over all frequency range. The SW component spectrum shows a straight steep spectrum over our frequency range, with
and
= 0.93 ± 0.14. The NE component has spectral indices of
and
= 1.08 ± 0.14. The NW component is detected in the 144 MHz and 4700 MHz image; however, it is not detected in the 1400 MHz image of Cole et al. (1999) and very faintly detected in the 8415 MHz image. We estimate a spectral index of
for this component.
3.3.2. Large-scale emission
Southern arc. In the integrated spectral index of the southern arc emission (excluding the nuclear region), we find no significant steepening of the spectrum up to 4700 MHz, as ,
and
(see Fig. 5 for spectral profiles at different locations in the southern arc). A relative steepening at the high frequency end would be expected in a radiatively ageing population, provided sufficient time has passed. To spatially resolve the spectral properties, we have constructed spectral index maps shown in Figs. 4 and 6 from 144−4700 MHz. The maps show several interesting features. The first property to note is that while no spectral gradient is observed along the arc, a gradient across the arc can be seen in the
map in Fig. 4a – from 0.8−0.9 (±0.13) at the inner edge of the arc to 0.5−0.6(±0.12) at the outer edge. The spectral index map made with LOFAR 144 MHz and Apertif 1360 MHz image at ∼16″ resolution also showed a similar gradient, thus increasing our confidence in the presence of this feature. In the
map, we detect a steep spectral index of ∼1.1 ± 0.2 at the tail end of the arc. The spectral index from 1455−4700 MHz at the tail end of the southern arc is steeper than the middle region, thus we observe a gradient along the arc at this frequency range. The second property to note is the spectral index of the ridge
, which is steeper than the typical spectral index measured in the arc. The average spectral index for the ridge in the
map is 0.96 ± 0.10. In these spectral index maps, the nuclear region shows very little curvature, with an average spectral index of
and
.
![]() |
Fig. 6. 20″ spectral index map from 1455−4700 MHz made with the VLA images. The local rms noise, σrms, in the 1455 MHz image is 80 μJy beam−1 and in the 4700 MHz image is 40 μJy beam−1. The contours mark the 4700 MHz emission at levels (3, 5, 10, 15, 20, 50)×σrms. The inset map shows the spectral index error. |
Northern arc. The sensitivity of the Apertif image allowed us to trace the spectral indices in the northern arc region. As can be seen in Fig. 4b, the northern arc has a typical spectral index of , within the 3σ 1360 MHz contours. This could also have contribution from the background source marked in the figure. Outside the contours,
. Overall, this suggests that the spectral index of the northern arc is likely steeper than its southern counterpart. We do not recover any emission in the northern arc in the 4700 MHz image, and do not discuss its properties further. We discuss the spectral properties of Mrk 273 further in Sect. 4.
3.4. Magnetic field and spectral age
To obtain the spectral age of the large-scale radio emission in Mrk 273, we need to first estimate the magnetic field strength. We have used the python version of the SYNCH code5 (Hardcastle et al. 1998) to estimate the equipartition magnetic field in different regions in the source. The regions are marked in Fig. 7. We assume that star formation has no contribution to the radio emission in these regions. Typically, while estimating the equipartition magnetic field, the ratio of energy in non-radiating (protons) to radiating particles is assumed to be zero (κ = 0), that is, no non-radiating particles are present. However, studies have shown that in more turbulent lobes of low power radio galaxies, with radio luminosities comparable to Mrk 273, internal lobe pressure is not sufficient to balance the pressure of the external medium (Croston et al. 2003; Dunn & Fabian 2004; Croston & Hardcastle 2014). The extra pressure needed is provided by either the higher magnetic field or the entrained material, that is, the non-radiating particles. Croston et al. (2018) found that a ratio of ∼10 between the non-radiating to radiating particles was required for the lobe pressure to balance the external pressure. Bicknell et al. (1998) also found that the energy density of radio jets in low power AGNs (like Seyferts) was dominated by thermal plasma. This is likely to be the case for the radio emission in Mrk 273, given its low luminosity (L1.4 GHz ≈ 5 × 1023 W Hz−1) and Seyfert 2 classification. Therefore, we have estimated the equipartition values, reported in Table 3, using both κ = 0 and κ = 10. Our estimates are in agreement with magnetic field strengths found in similar power radio galaxies (Heesen et al. 2018; Das et al. 2021).
![]() |
Fig. 7. VLA 1455 MHz image of Mrk 273 at a resolution of 20″ × 20″. The solid lines mark the positive contour levels at: 3 |
Flux densities and model fit results for southern arc and ridge regions
We have extracted the flux density from all the regions, and used an aged synchrotron spectrum that provided a reasonable fit to the flux densities. For this purpose, we have used all the images at 20″ resolution. We assumed a cylindrical geometry for all the regions, shown in Fig. 7. For regions S1 and S2 in the southern arc, we used a length of 30″ and radius of 15″. For the region in the ridge, R, we used a length of 50″ and radius of 8″. The injection index was varied over the typical range of 0.5−0.7, with steps of 0.01, using the Broadband Radio Astronomy Tools (BRATS6; Harwood et al. 2013, 2015) package to find the best fit injection index (αinj) for the integrated spectra of the regions. The best fit values for αinj were 0.55 for S1, 0.61 for S2 and 0.65 for R. These indices were then used to estimate the magnetic field strengths. These values are listed in Table 3.
The spectrum of particles emitting synchrotron radiation steepens due to preferential cooling of higher energy particles (Kellermann 1966; Pacholczyk 1970). Modelling the steepening of the radio spectra can therefore give us an estimate of the spectral age. If a break exists in the spectra shown in Fig. 5, it is above ∼5 GHz. To model these spectra and estimate the spectral ages, we fitted the JP model (Jaffe & Perola 1973) to the integrated spectra, using the task fitintegrated in BRATS. We note that there are discrepancies in spectral age between estimates using integrated spectra and a spatially resolved analysis, but our results should still be reliable within the regions. However, these estimates may change for an analysis at lower spatial resolution (Harwood 2017). The model fit results are summarised in Table 3 and shown in Fig. 5. We find ages of a few tens of Myr in all the regions. Despite the steeper spectra of the ridge compared to the southern arc, the ages are similar for both due to the higher magnetic field in the ridge. A gradient in spectral age can also be seen from S2 to S1, which is discussed further in Sect. 4.2.1.
In the lobes of radio galaxies, the equipartition magnetic field strength is related to the minimum energy. This can give us information about the properties of the external medium. In order to estimate these properties (and the possible impact of the emission) we use the relation given as
where umin is the minimum energy in erg cm−3 and Beq is in Gauss. The particle energy density, derived using the equation above, can be used to approximate the particle density of the surrounding medium. We estimated the relativistic gas pressure as Pr = (γ − 1)umin, where γ = 4/3 for relativistic gas (Pacholczyk 1970), and assumed that the regions are in pressure balance with their surrounding medium. The gas pressure of the surrounding medium can be approximated as Pm = η kT, where η is the particle density of the medium (Mack et al. 1998). Using the particle energy density obtained above, we roughly estimate a medium particle density of ηR = 1 × 10−3 cm−3 for T = 107 K for the ridge region. This is comparable to the electron number density estimated for the southern X-ray nebula (ηe = 3 × 10−3 cm−3) that overlaps with the ridge (Fig. 1), by L19. For the southern arc regions, we estimate a particle density of ∼3 × 10−4 cm−3, an order of magnitude lower than that found in the ridge. Therefore, if such a surrounding medium is present, the low particle density would explain why it is not detected in the soft X-ray images at their current sensitivity.
4. Discussion
We have traced a rare case of large-scale (∼300 kpc) radio continuum emission in a ULIRG. The LOFAR Dutch array images have allowed us to trace spectacular and smoothly curved arcs out to ∼100−200 kpc, in the north and south, while on kiloparsec scale a complex structure has been traced by the LOFAR International stations data. Using multi-frequency images, we derived the spectral properties from MHz to GHz frequencies, something not done in the very few cases known so far. These data now allow us to explore the properties of the merger system. We first discuss the radio properties of the nuclear components, and explore the possibility of the presence of a radio AGN in this system. Then we move outwards to the large-scale emission and discuss the different scenarios that could explain its origin, spectral properties and how they connect to the nuclear components. Finally, we discuss what this system tells us about presence of nuclear activity in mergers.
4.1. Nuclear components
The sub-arcsecond images show multiple components in the ∼4 kpc nuclear region of Mrk 273 (Fig. 2). The multi-wavelength properties of some of these components have been mentioned before in the introduction. The presence of absorbed X-ray AGNs has been confirmed in the N and SW component (Iwasawa et al. 2011; L19), and an AGN has also been proposed in the SE component in the optical emission line study of RZ14. In the radio, an AGN had also been proposed in the nucleus of Mrk 273 by Vardoulaki et al. (2015) but it could not be pinned down to one nuclear component.
-
N component: The radio emission of this component is understood to be dominated by a starburst (Bondi et al. 2005). However, Bondi et al. (2005), in their VLBI study of the nucleus, could not rule out the possibility that one of the compact components could be the radio counterpart of the AGN. As mentioned in Sect. 3.3.1, we detect broadband radio absorption in the N component spectrum. Although the presence of absorption in the nuclear region had been suggested by Clemens et al. (2010), they did not observe a turnover in the spectrum and could not pinpoint the location of the absorbed component. We note that these differences with respect to our results could be due to the low resolution (5″ and 13″) of their images, compared to our sub-arcsecond images of the nucleus. Using our low frequency sub-arcsecond image of the nucleus, we have detected the absorption to be in the N component. Investigating the absorption mechanism, either SSA (synchrotron self-absorption) or FFA (free–free absorption), can give us more information about this component. SSA requires high brightness temperatures, comparable to the kinetic temperatures of ∼1010 K of the relativistic electrons (Broderick & Condon 1975). Condon et al. (1991) estimated the brightness temperatures of the nuclear components using 0.3″ images of Mrk 273 and concluded that the brightness temperatures (∼6 × 103 K) were too low for SSA to be significant. On the other hand, a turnover in the spectrum is also expected in compact, nuclear starbursts due to optically thick FFA (Condon et al. 1991; Clemens et al. 2010; Leroy et al. 2011; Murphy 2013). We therefore consider FFA to be the dominant absorption mechanism for the N component, as also suggested by Clemens et al. (2010). We used an FFA model that assumes that the absorption is caused by an external homogeneous ionised gas screen, given as
where a and α are the amplitude and spectral index of the intrinsic synchrotron spectrum, and τν is the optical depth equal to (ν/νp)−2.1, where νp is the frequency at which the optical depth is unity. We fitted this model to the N component spectrum and found best fit parameters of a = 0.11 ± 0.01, α = 0.51 ± 0.03 and νp = 0.19 ± 0.01 GHz. The model fit is shown in Fig. 8. For an SSA model fit, we obtained similar model parameters and quality of fit to the data. Using the FFA parameters, we obtained an emission measure of 9.3 × 104 pc cm−6 (O’Dea 1998), ∼30 times lower than Clemens et al. (2010). This is due to their assumption of a turnover at ∼1 GHz whereas we estimate the turnover to be at ∼370 MHz. Using the deconvolved size from Condon et al. (1991) as the path length for the N component, we estimate an electron density of ∼25 cm−3 for the ionised absorbing medium, which corresponds to a column density of Ne ≈ 5 × 1021 cm−2. This is only a fraction of the column density of neutral hydrogen of NH = 1.7 × 1022 cm−2 (Cole et al. 1999) distributed in a ∼250 pc disc around the starburst. The presence of absorption is interesting because Nandi et al. (2021) found absorption and low spectral ages in several ULIRGs, concluding that these are GPS/CSS sources. However, our analysis shows that in the nuclear component of a ULIRG, where the radio emission is dominated by a starburst, the absorption can also be explained by FFA from the surrounding medium.
-
SE component: The SE component shows a “normal” spectral index at low frequencies,
∼ 0.75 ± 0.05, which steepens to
∼ 1.14 ± 0.14. Spectral index of star forming components are not expected to be steeper than ∼0.83 ± 0.13 (Niklas et al. 1997). In their study of spectral properties of ULIRGs, Vardoulaki et al. (2015) proposed the presence of a radio AGN in the nucleus of Mrk 273 on the basis of a spectral index steeper than that expected from star formation, from 1.4−8.4 GHz. The steep spectral index of the SE component leads us to conclude that the radio AGN proposed by Vardoulaki et al. (2015) is associated with this component. Its radio luminosity of L(144 MHz) = 4 × 1023 W Hz−1 would make it a low power radio AGN. The idea that this component is a separate nucleus with its own AGN was also suggested by RZ14 on the basis of the high optical ionisation state and compact continuum emission. The diffuse emission extending southwards from this component (see Fig. A.1e) could be the connection to the southern large-scale features i.e., the arc or the ridge.
-
SW and NE component: The LOFAR 144 MHz sub-arcsecond image shows for the first time, resolved structure in the SW component, in the form of two bright spots (see Fig. 2). Some diffuse emission can be seen around this component in the lower resolution image in Fig. A.1e. This component is aligned with the more diffuse emission in the NE component. The SW component is also aligned in the direction of the western outflow, detected by RZ14 (see their Fig. 1) and Leung et al. (2021). Both of these components show steep spectral indices varying from ∼0.8−1.0. The alignment of these two components and the overlap with optical outflows suggests that this emission could represent outflowing material. The outflows could carry the radio emitting electrons from the nuclear components, or could cause shock-acceleration of the electrons present in the nuclear region. At sub-arcsecond resolution with LOFAR, such steep-spectrum diffuse emission from radio outflows has been detected before in multiple objects, for example Arp 299 (Ramírez-Olivencia et al. 2018, 2021), Arp 220 (Varenius et al. 2016) and M 82 (Varenius et al. 2015).
![]() |
Fig. 8. Homogeneous FFA model best fit for the N component flux densities. A turnover can be seen at ∼370 MHz. |
Overall, we propose that the radio AGN in Mrk 273 is located in the SE component. This adds further evidence for the presence of three AGNs in the nucleus of Mrk 273. The low resolution image in Fig. A.1e shows diffuse emission around the nuclear components that had been missed by the previous sub-arcsecond images at GHz frequencies. The southward extension from the SE component, could be the connection between the radio AGN in the SE component and the southern arc or ridge.
4.2. Large-scale
The large-scale structure can be described as a combination of three components – the arcs in north and south, and the ridge. The direction of the southern arc, partly consistent with the optical tail at that location (see Fig. 1), and perpendicular to the ridge, suggests that the southern arc and the ridge are not physically related. L19 also proposed that the soft X-ray nebula, that overlaps with the ridge, is unrelated to the tidal tail that casts a shadow over the nebula emission. Therefore, we discuss these separately below.
4.2.1. Southern arc
Starbursts and AGNs are known to be the dominant power source in ULIRGs. In Mrk 273, the high SFR and the presence of at least two to possibly three AGNs in the nuclear region makes it necessary to explore both as the powering mechanisms of the radio emission. Below, we discuss the various scenarios for the origin of the southern arc emission.
Star formation. Radio emission from galaxies also has contribution from star formation, which emits both thermal (free-free) and non-thermal (synchrotron) radiation. The thermal component dominates above ∼30 GHz, and has a spectral index of ∼0.1. The non-thermal emission, from relativistic electrons accelerated by supernova remnants, dominates the spectrum up to ∼10 GHz. Assuming that all the radio emission in the arc at 1.4 GHz is from star formation, and using the calibration from Yun et al. (2001), we estimate an SFR of 24.2 ± 7.4 M⊙ yr−1. Although such high SFRs are seen in some ULIRGs, they are not expected at distances of ∼100 kpc from the nucleus, in the form of an arc. This SFR corresponds to an Hα surface brightness of ∼6 × 10−17 erg s−1 cm−2 arcsec−2 (Kennicutt 1998). S16 study the Hα emission in Mrk 273, and have a sensitivity of ∼1 × 10−17 erg s−1 cm−2 but recover no Hα emission coinciding with the southern arc, out to 100 kpc. There is also no optical continuum emission detected along the arc from the stars being formed. No emission from dust has been detected along the arc either, which could hide the optical emission. The optical continuum and Hα emission is only seen out to a radial extent of ∼45 kpc to the south of the nucleus (S16), which overlaps with the radio ridge. This tells us that star formation in the arc is not the origin of the radio emission.
Tidal origin. Another possibility is that the southern (and northern) arc have tidal origin. In this scenario, the radio emission would be associated with supernovae from the stars stripped from the merging galaxies. However, a lack of emission from stars (no optical continuum) and gas (no HI or ionised gas emission) over these regions makes it unlikely (RZ14,S16; Leung et al. 2021). The dynamical timescale of the merger is also much larger than the typical synchrotron lifetime of ∼40 Myr for the emitting electrons. Therefore, if the southern arc was made up of tidally stripped relativistic plasma left behind as the merger progressed, there would be more signs of spectral ageing, which we do not detect. This suggests that the southern arc does not have a tidal origin.
Galactic winds. Nuclear starbursts and AGNs in galaxies can also drive large-scale galactic winds, that can be detected at different wavelengths (Heckman et al. 1990; Veilleux et al. 2005). In the last few years, diffuse emission from starburst driven galactic winds has been detected at large radial distances from the nucleus, including at ∼60 kpc (Hodges-Kluck et al. 2020) and ∼50 kpc (Rupke et al. 2019), although the actual sizes could be even larger since these observations are limited by their field of view.
In the case of Mrk 273, recent wide-field integral field unit observations of optical emission lines have yielded no such detections (Leung et al. 2021) out to a radial extent of ∼40 kpc. Our high sensitivity low frequency radio observations could be more sensitive to such winds. Therefore, the southern arc could be a result of galactic winds originating from the nucleus, ploughing into the external medium. The curvature of the southern arc could be caused by the interaction of these winds with an external medium. We refer to this as Scenario A. Given the presence of powerful AGNs in the nucleus, we propose that these winds could be driven by AGNs, but some contribution from the nuclear starburst cannot be ruled out. However, compared to known cases of large-scale galactic winds, the southern arc has a much more collimated morphology and a significantly larger total spatial extent of ∼300 kpc (southern + northern arc). Although we cannot rule out this scenario, the morphological properties reduce our confidence in it.
AGN. The more likely case is that the southern arc is fuelled with relativistic electrons by an AGN in Mrk 273. The luminosity of the southern arc, L(144 MHz) = 2×1023 W Hz−1, is typical of radio lobes in AGNs, although it lies on the "radio quiet" end of the population (Panessa et al. 2019). Yun et al. (1999) have also suggested that the luminosity and high degree of polarisation of the extended emission in Mrk 273 points towards an AGN origin. In this case, the southern arc would be fuelled by a radio AGN in the nucleus, i.e. the SE component. However, the radio AGN could have either already been present in one of the merging galaxies − which we refer to as Scenario B − or could have been triggered by the merger − which we refer to as Scenario C.
-
In Scenario B, the smooth curvature of the arc would then trace the infall of the radio AGN. But such plasma left behind as the merger progresses would also be expected to show signs of strong (α ≥ 1.2 or Δα ≥ 0.5) spectral ageing, as discussed above, which we do not detect. We do find that the spectral index is steeper at the tail end (S2) than at the leading end (S1, Sect. 3.3.2). The spectral age of the arc emission is also higher, 36−42 Myr, at the tail end (S2) than at the leading end (S1), 20−22 Myr. If we assume that the age difference is due to the motion of the infalling galaxy as mentioned above, and use the difference in the spectral ages as the time taken by the infalling galaxy to cover the distance between the regions, we can approximate, to first order, an infall speed of the galaxy (v≈ d/tspec). Using 55 kpc for the distance between the two regions, we estimate an infall speed of ∼2500−3300 km s−1. During interaction of galaxies in non-cluster environments, the velocities are expected to be of the order of the typical gravitational motions in galaxies, i.e., several hundred km s−1. For instance, the collisions in Taffy galaxies are understood to occur at ∼600−800 km s−1 (see Condon et al. 1993; Vollmer et al. 2012). However, the velocity we estimate, albeit a first order approximation, is much higher than seen before for interacting galaxies. Therefore, if scenario B were true, it points to other processes, like in situ particle reacceleration, being present in the arc. Such processes have also been suggested to be present in the nearby ULIRG Mrk 231.
-
In Scenario C, the southern arc would be a radio lobe fuelled by a low power radio jet from AGN activity triggered by the merger, in the pre-coalescence stage. The low power means that the lobe could be more turbulent and less well collimated than in typical radio galaxies with powerful and highly collimated radio jets. A lobe fuelled by a radio jet can explain the large spatial extent, since the energy flux can be channelled better than in the case of winds, allowing the plasma to travel greater distances with relatively less energy loss. A radio jet would also produce a much more collimated morphology than AGN driven winds. The smooth curvature of the southern arc in this case would be due to interaction with the large-scale external medium bending the low power radio lobe. We note that we do not directly detect the presence of a radio jet in Mrk 273 in our sub-arcsecond images. This could be due to the jet getting decollimated, as it makes its way out of the nucleus, to a level where its surface brightness is below the sensitivity of our sub-arcsecond images. However, we do see an extension of the SE component towards the southern arc, in the smoothed LOFAR sub-arcsecond image (Fig. A.1e), as mentioned in Sect. 3.1.1.
A new feature observed in the southern arc is the presence of a spectral gradient across the arc, with α ∼ 0.6 at the outer edge and ∼0.9 at the inner edge (Fig. 4a). We also note an edge-brightened feature at the outer edge in the 6″ contour plot at 144 MHz. This suggests that the plasma is interacting with an external medium at this location, which rejuvenates the electrons at the outer edge. This would be expected if the plasma were moving radially outwards, like the wind in Scenario A or radio jet in Scenario C, rather than being formed locally as expected in Scenario B. Therefore, our spatially resolved spectral properties help us rule out Scenario B. Assuming that the southern arc plasma is in pressure balance with such a surrounding medium, we estimate the particle density of the medium to be lower (by an order of magnitude) than the southern X-ray nebula (L19). This could be a reason why such a large-scale medium has not been detected in the current soft X-ray images.
Using our spectral age estimate at the tail end of the southern arc (S2), we can also approximate, to first order, the speed of the galactic winds required in Scenario A or radio lobe in Scenario C. Assuming constant velocity, a distance of 100 kpc from the nucleus for S2, and a time range of 36−42 Myr gives a velocity range of ∼2300−2700 km s−1. For galactic winds in Scenario A, such a high speed at a distance of 100 kpc from the nucleus is highly unlikely. For instance, the speed of galactic wind at a distance of ∼50 kpc in SDSS J211824.06+001729.4 (Makani; Rupke et al. 2019) is only 120 km s−1, and for NGC 3079 at a distance of ∼60 kpc it is 500 km s−1 (Hodges-Kluck et al. 2020). However, for radio lobe in Scenario C, our estimate corresponds to 0.5−0.7% of the speed of light. This is in agreement with the current understanding of radio jets, which are understood to decelerate to sub-relativistic speeds at tens of kiloparsecs distance from the nucleus (Laing & Bridle 2002, 2014; Hardcastle & Croston 2020). Therefore, we propose Scenario C as a more favourable explanation for the origin of the southern arc.
As mentioned before, large-scale radio continuum emission has also been detected in the ULIRG Mrk 231, also known to host an AGN and a starburst. Morganti et al. (2016) had concluded that the large radio lobes could be fuelled by electrons from the central AGN by or from electrons accelerated in situ. Although Morganti et al. (2016) did not directly detect a radio jet from their data either, they detected a “bridge” of emission connecting the lobe to the nucleus, which supports the fuelling from central AGN hypothesis. Although the radio power of Mrk 273 is lower than Mrk 231 (L1.4 GHz = 1.4 × 1024 W Hz−1), the presence of multiple AGNs in the nucleus, including one radio AGN (in the SE component), suggests that processes similar to Mrk 231 could also be responsible for the radio emission in the southern arc of Mrk 273.
4.2.2. Northern arc
We detect, for the first time, a giant ∼190 kpc arc in the north of the nuclear region of Mrk 273, which is much larger than its southern counterpart (∼100 kpc). This component is also more diffuse and has lower surface brightness than the southern arc. It also does not show a smooth curvature. A background galaxy 2MFGC 11079 at z = 0.038 also contributes to the radio emission in this arc, marked in Fig. 2. In Scenario A, the northern arc would be formed by the counter-wind as it expands into the external medium. Scenario B cannot explain the origin of the northern arc, since the southern arc is created locally in that scenario. In Scenario C, the northern arc would be formed by the radio lobe fuelled by the counter jet.
In Sect. 3.3.2, we found that the spectral index in this arc is likely steeper than the southern arc. Given the low power of the radio AGN involved and the less collimated structure of the northern arc, it is unlikely that this is a result of relativistic beaming. Instead, the less collimated structure and larger size of the northern arc could be due to different conditions of the surrounding medium (i.e., halo material, circumgalactic medium) around it, compared to the south. A less dense medium in this region would mean that the wind (Scenario A) or radio lobe (Scenario C) could travel larger distances in the same time as the southern arc. The scarcity of interactions with a dense surrounding medium could cause a lack of rejuvenated electrons responsible for the radio emission, leading to the lower surface brightness and steeper spectral index (since the feature is not detected at 4700 MHz). The lack of curvature in this arc would also be explained by the lack of a medium in this location to bend the plasma. Such asymmetric radio emission has also been seen before in the ULIRG Mrk 231 (Morganti et al. 2016).
4.2.3. Ridge
A peculiar feature in Mrk 273 is the ridge of radio emission, at a distance of ∼25 kpc south of the nucleus. This overlaps very well with the soft X-ray emission, shown in Fig. 1, and also the ionised gas nebula observed in this region (RZ14, S16). The emission in the nebula is understood to be due to shocks in this region. RZ14 and S16 propose that the gas in the nebula is ionised by low velocity shocks (200−400 km s−1) and L19 conclude that the gas is heated to X-ray emitting temperatures by shocks that sweep through the medium. L19 also propose that the X-ray nebula in this region would have to be enriched by multiple outflows. We propose that the radio emission in this region is due to a reservoir of electrons deposited by multiple outflows, and then accelerated by the shocks that sweep the medium. The spectral age of the ridge, ∼20−32 Myr, is also in agreement with the timescale of enrichment of the nebula by multiple outflows, ≲0.1 Gyr, estimated by L19. The radio ridge also overlaps with the tidal tail seen in optical continuum, and therefore has some contribution from the star formation in the tail.
The systematically steeper spectral index of the ridge in comparison to the arc (≈2−3σ significance) suggests that these components are formed by two separate electron populations. The electron population that makes the ridge has a steeper energy distribution, i.e., less high energy electrons relative to the arc, leading to a steeper spectral index. This could be due to the different powering mechanisms for the ridge and the southern arc. The southern arc is fuelled by a radio AGN, which would be more efficient at accelerating electrons than the shocks present in the ridge.
We do not observe any radio emission to the north-east of the nucleus, where another ionised gas and soft X-ray nebula is seen (Fig. 1). RZ14 find quiescent gas kinematics in this region and conclude that the ionised gas represents tidal debris left over from a secondary merger event illuminated by one of the AGNs from the nucleus. L19 also propose that this nebula is not enriched by multiple outflows. The lack of radio emission then, would be in line with these observations, since there would be no pool of electrons deposited by outflows in this region and no shocks to accelerate them to synchrotron emitting energies.
4.3. Mrk 273 in the broader context of ULIRGs
Mrk 273 shows spectacular continuum and spectral properties on both large and small-scales. Although extended optical continuum and HI emission is known to occur in major mergers, radio continuum from such objects is still very rare (for a review, see Hibbard et al. 2001). Similarly, very few cases of extended radio emission have been identified in ULIRGs (Y99; Hayashi et al. 2021; Nandi et al. 2021), and even fewer have been spatially resolved. To the best of our knowledge, this is the first time the spectral properties of the resolved extended emission from an ULIRG have been studied (similar analysis was done for the LIRG Arp 299 by Ramírez-Olivencia et al. 2021). We propose that the southern arc is a formed by a radio AGN in the nuclear region of Mrk 273, and a star formation or tidal origin for the emission are unlikely. On the basis of the steep spectral index, we conclude that this radio AGN is associated with the SE component in the nucleus. In their high resolution study of the nuclear components, Bondi et al. (2005) also suggested that the morphology of the SE component resembled a radio AGN.
The presence of radio activity in Mrk 273 can be useful to link radio AGNs and ULIRGs. We find that fuelling from a radio AGN can explain the formation of both the arcs, where the radio activity has been triggered by the merger. This shows that activity can be triggered even before the nuclei have coalesced and strengthens the current understanding that activity can be triggered at different stages of the merging process (Tadhunter et al. 2011). Another link between ULIRGs and radio AGNs is that in the local universe, ULIRGs may be the early stage in the evolution of a merger into a radio galaxy (Nandi et al. 2021). Our results for Mrk 273 support this hypothesis. Although it does not appear very likely, another possibility is that the arcs are formed by AGN driven high speed winds. However, if this is the case, then this would be one of the largest known cases of galactic winds, with a total spatial extent of ∼300 kpc. In any case, our results, combined with the other ULIRGs mentioned in Sect. 1, give further evidence for the idea that large (∼300 kpc) radio structures in ULIRGs are signposts of AGNs in the system.
To further test this hypothesis, it is important that a large number of ULIRGs are inspected for the presence of large-scale low surface brightness radio emission. For this purpose, observations with high sensitivity to low-surface brightness features are needed. As illustrated by this work, low frequency surveys are particularly suited for this and, in particular, the ongoing LOFAR surveys LoLSS and LoTSS can provide these data. However, our work has also illustrated the importance of the spectral analysis for understanding the origin of the radio structure in ULIRG.
In the case of the present work, the addition of the image from the Apertif survey has allowed us to derive the (spatially resolved) spectral properties of the extended emission in Mrk 273. The subsequent confirmation (and refinement) provided by the pointed VLA observations is an important step. This illustrates the potential of the LOFAR-Apertif combination for studying a large sample of these objects. Other surveys, in particular EMU at ∼800 MHz, will provide comparable images for an even larger sky area. Studying the spectral properties is useful because it can help pin down the origin of the radio emission and provide an estimate for the timescale of merger triggered activity. This information would be crucial for our understanding of the relation between mergers and AGNs.
Acknowledgments
We would like to thank the referee, Loreto Barcos-Munoz, for their comments, especially on the VLA C band data reduction, which has significantly improved the quality of the paper. LOFAR, the Low Frequency Array designed and constructed by ASTRON, has facilities in several countries, that are owned by various parties (each with their own funding sources), and that are collectively operated by the International LOFAR Telescope (ILT) foundation under a joint scientific policy. This work makes use of data from the Apertif system installed at the Westerbork Synthesis Radio Telescope owned by ASTRON. ASTRON, the Netherlands Institute for Radio Astronomy, is an institute of the Dutch Science Organisation (De Nederlandse Organisatie voor Wetenschappelijk Onderzoek, NWO). This work was supported by the Medical Research Council [MR/T042842/1]. CT acknowledges support from STFC. AD acknowledges support by the BMBF Verbundforschung under the grant 05A20STA. RJvW acknowledges support from the ERC Starting Grant ClusterWeb 804208. FdG acknowledges support from the Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy – EXC 2121 “Quantum Universe” – 390833306. EAKA is supported by the WISE research programme, which is financed by the Dutch Research Council (NWO). JvL acknowledges funding from Vici research programme ‘ARGO’ with project number 639.043.815, financed by the Dutch Research Council (NWO). LCO acknowledges funding from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement No. 617199. DV acknowledges support from the Netherlands eScience Center (NLeSC) under grant ASDI.15.406. JMvdH acknowledges funding from the Europeaní Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement No. 291531 (‘HIStoryNU’). 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). KMH acknowledges funding from the State Agency for Research of the Spanish Ministry of Science, Innovation and Universities through the “Center of Excellence Severo Ochoa” awarded to the Instituto de Astrofísica de Andalucía (SEV-2017-0709); from grant RTI2018-096228-B-C31 (Ministry of Science, Innovation and Universities/State Agency for Research/European Regional Development Funds, European Union); and from the coordination of the participation in SKA-SPAIN, funded by the Ministry of Science and innovation (MICIN).
References
- Adebahr, B., Schulz, R., Dijkema, T. J., et al. 2022, Astron. Comput., 38, 100514 [NASA ADS] [CrossRef] [Google Scholar]
- Barcos-Muñoz, L., Leroy, A. K., Evans, A. S., et al. 2017, ApJ, 843, 117 [CrossRef] [Google Scholar]
- Baum, S. A., O’Dea, C. P., Dallacassa, D., et al. 1993, ApJ, 419, 553 [NASA ADS] [CrossRef] [Google Scholar]
- Bessiere, P. S., Tadhunter, C. N., Ramos Almeida, C., & Villar Martín, M. 2012, MNRAS, 426, 276 [NASA ADS] [CrossRef] [Google Scholar]
- Bicknell, G., Dopita, M., Tsvetanov, Z., & Sutherland, R. 1998, ApJ, 495, 680 [NASA ADS] [CrossRef] [Google Scholar]
- Bondi, M., Pérez-Torres, M.-A., Dallacasa, D., & Muxlow, T. W. B. 2005, MNRAS, 361, 748 [NASA ADS] [CrossRef] [Google Scholar]
- Broderick, J. J., & Condon, J. J. 1975, ApJ, 202, 596 [NASA ADS] [CrossRef] [Google Scholar]
- Cano-Díaz, M., Maiolino, R., Marconi, A., et al. 2012, A&A, 537, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Carilli, C. L., & Taylor, G. B. 2000, ApJ, 532, 95 [Google Scholar]
- Carniani, S., Marconi, A., Maiolino, R., et al. 2016, A&A, 591, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&A, 562, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Clemens, M. S., Scaife, A., Vega, O., & Bressan, A. 2010, MNRAS, 405, 887 [NASA ADS] [Google Scholar]
- Cole, G. H., Pedlar, A., Holloway, A. J., & Mundell, C. G. 1999, MNRAS, 310, 1033 [NASA ADS] [CrossRef] [Google Scholar]
- Condon, J. J., Huang, Z. P., Yin, Q. F., et al. 1991, ApJ, 378, 65 [NASA ADS] [CrossRef] [Google Scholar]
- Condon, J. J., Helou, G., Sanders, D. B., & Soifer, B. T. 1993, AJ, 105, 1730 [Google Scholar]
- Cresci, G., Mainieri, V., Brusa, M., et al. 2015, ApJ, 799, 82 [Google Scholar]
- Croston, J. H., & Hardcastle, M. J. 2014, MNRAS, 438, 3310 [Google Scholar]
- Croston, J. H., Hardcastle, M. J., Birkinshaw, M., & Worrall, D. M. 2003, MNRAS, 346, 1041 [NASA ADS] [CrossRef] [Google Scholar]
- Croston, J. H., Ineson, J., & Hardcastle, M. J. 2018, MNRAS, 476, 1614 [Google Scholar]
- Das, S., Kharb, P., Morganti, R., & Nandi, S. 2021, MNRAS, 504, 4416 [NASA ADS] [CrossRef] [Google Scholar]
- de Gasperin, F., Dijkema, T. J., Drabent, A., et al. 2019, A&A, 622, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Gasperin, F., Brunetti, G., Bruggen, M., et al. 2020, A&A, 642, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Gasperin, F., Williams, W. L., Best, P., et al. 2021, A&A, 648, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- De Lucia, G., Springel, V., White, S. D. M., et al. 2006, MNRAS, 366, 499 [NASA ADS] [CrossRef] [Google Scholar]
- Di Matteo, T., Springel, V., & Ilernquist, L. 2005, Nature, 433, 604 [NASA ADS] [CrossRef] [Google Scholar]
- Dunn, R. J. H., & Fabian, A. C. 2004, MNRAS, 355, 862 [NASA ADS] [CrossRef] [Google Scholar]
- Emonts, B. H. C., Morganti, R., Tadhunter, C. N., et al. 2006, A&A, 454, 125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Greisen, E. W. 2003, Information Handling in Astronomy - Historical Vistas (Netherlands: Springer), 109, 285 [Google Scholar]
- Hales, S. E. G., Masson, C. R., Warner, P. J., et al. 1990, MNRAS, 246, 256 [NASA ADS] [Google Scholar]
- Hardcastle, M. J., & Croston, J. H. 2020, New Astron. Rev., 88, 101539 [Google Scholar]
- Hardcastle, M. J., Birkinshaw, M., & Worrall, D. M. 1998, MNRAS, 294, 615 [CrossRef] [Google Scholar]
- Harwood, J. J. 2017, MNRAS, 466, 2888 [Google Scholar]
- Harwood, J. J., Hardcastle, M. J., Croston, J. H., & Goodger, J. L. 2013, MNRAS, 435, 3353 [NASA ADS] [CrossRef] [Google Scholar]
- Harwood, J. J., Hardcastle, M. J., & Croston, J. H. 2015, MNRAS, 454, 3403 [NASA ADS] [CrossRef] [Google Scholar]
- Hayashi, T. J., Hagiwara, Y., Imanishi, M., et al. 2021, MNRAS, 504, 2675 [NASA ADS] [CrossRef] [Google Scholar]
- Heald, G., Mckean, J., Pizzo, R., et al. 2010, in ISKAF, 57 [Google Scholar]
- Heckman, T. M., Armus, L., & Miley, G. K. 1990, ApJS, 74, 833 [Google Scholar]
- Heesen, V., Croston, J. H., Morganti, R., et al. 2018, MNRAS, 474, 5049 [Google Scholar]
- Hibbard, J. E., van Gorkom, J. H., Rupen, M. P., & Schiminovich, D. 2001, ASP Conf. Ser., 240, 659 [NASA ADS] [Google Scholar]
- Hodges-Kluck, E. J., Yukita, M., Tanner, R., et al. 2020, ApJ, 903, 35 [NASA ADS] [CrossRef] [Google Scholar]
- Hopkins, P. F., Hernquist, L., Cox, T. J., et al. 2006, ApJS, 163, 1 [Google Scholar]
- Hopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008, ApJS, 175, 356 [Google Scholar]
- Iwasawa, K., Mazzarella, J. M., Surace, J. A., et al. 2011, A&A, 528, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jackson, N., Tagore, A., Deller, A., et al. 2016, A&A, 595, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jackson, N., Badole, S., Morgan, J., et al. 2021, A&A, 658, A2 [Google Scholar]
- Jaffe, W. J., & Perola, G. C. 1973, A&A, 26, 423 [NASA ADS] [Google Scholar]
- Kellermann, K. I., & Kellermann, I. K., 1966, Aust. J. Phys., 19, 195 [NASA ADS] [Google Scholar]
- Kennicutt, R. C., Jr 1998, ApJ, 498, 541 [Google Scholar]
- Kim, D. C., & Sanders, D. B. 1998, ApJS, 119, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Knapen, J. H., Laine, S., Yates, J. A., et al. 1997, ApJ, 490, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Laing, R. A., & Bridle, A. H. 2002, MNRAS, 336, 328 [NASA ADS] [CrossRef] [Google Scholar]
- Laing, R. A., & Bridle, A. H. 2014, MNRAS, 437, 3405 [NASA ADS] [CrossRef] [Google Scholar]
- Leroy, A. K., Evans, A. S., Momjian, E., et al. 2011, ApJ, 739, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Leung, G. C. K., Coil, A. L., Rupke, D. S. N., et al. 2021, ApJ, 914, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, W., Veilleux, S., Iwasawa, K., et al. 2019, ApJ, 872, 39 [Google Scholar]
- Lonsdale, C., Farrah, D., & Smith, H. 2006, Astrophysics Update 2 (Springer Praxis Books), 285 [CrossRef] [Google Scholar]
- Mack, K. H., Klein, U., O’Dea, C. P., et al. 1998, A&A, 329, 431 [NASA ADS] [Google Scholar]
- McMullin, J. P., Waters, B., Schiebel, D., et al. 2007, ASP Conf. Ser., 376, 127 [NASA ADS] [Google Scholar]
- Mihos, J. C., & Hernquist, L. 1996, ApJ, 464, 641 [Google Scholar]
- Morabito, L. K., Jackson, N., Mooney, S., et al. 2022, A&A, 658, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Morganti, R., Veilleux, S., Oosterloo, T., et al. 2016, A&A, 593, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Murphy, E. J. 2013, ApJ, 777, 58 [Google Scholar]
- Nagar, N. M., Wilson, A. S., Falcke, H., et al. 2003, A&A, 409, 115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nandi, S., Das, M., & Dwarakanath, K. S. 2021, MNRAS, 503, 5746 [NASA ADS] [CrossRef] [Google Scholar]
- Nardini, E., Risaliti, G., Watabe, Y., et al. 2010, MNRAS, 405, 2505 [Google Scholar]
- Niklas, S., Klein, U., & Wielebinski, R. 1997, A&A, 322, 19 [NASA ADS] [Google Scholar]
- Norris, R. P., Hopkins, A. M., Afonso, J., et al. 2011, PASA, 28, 215 [Google Scholar]
- O’Dea, C. P. 1998, PASP, 110, 493 [Google Scholar]
- Offringa, A. R., & Smirnov, O. 2017, MNRAS, 471, 301 [Google Scholar]
- Offringa, A. R., de Bruyn, A. G., Biehl, M., et al. 2010, MNRAS, 405, 155 [NASA ADS] [Google Scholar]
- Offringa, A. R., Van De Gronde, J. J., & Roerdink, J. B. 2012, A&A, 539, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Offringa, A. R., McKinley, B., Hurley-Walker, N., et al. 2014, MNRAS, 444, 606 [Google Scholar]
- Pacholczyk, A. G. 1970, Radio Astrophysics. Nonthermal Processes in Galactic and Extragalactic Sources (San Francisco: Freeman) [Google Scholar]
- Panessa, F., Baldi, R. D., Laor, A., et al. 2019, Nat. Astron., 3 [Google Scholar]
- Perley, R. A., & Butler, B. J. 2013, ApJS, 204, 19 [Google Scholar]
- Pierce, J. C. S., Tadhunter, C. N., Gordon, Y., et al. 2022, MNRAS, 510, 1163 [Google Scholar]
- Ramírez-Olivencia, N., Varenius, E., Pérez-Torres, M., et al. 2018, A&A, 610, L18 [CrossRef] [EDP Sciences] [Google Scholar]
- Ramírez-Olivencia, N., Varenius, E., Pérez-Torres, M., et al. 2022, A&A, 658, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ramos Almeida, C., Tadhunter, C. N., Inskip, K. J., et al. 2011, MNRAS, 410, 1550 [NASA ADS] [Google Scholar]
- Rigopoulou, D., Spoon, H. W. W., Genzel, R., et al. 1999, AJ, 118, 2625 [Google Scholar]
- Rodríguez Zaurín, J., Tadhunter, C. N., & González Delgado, R. M. 2009, MNRAS, 400, 1139 [Google Scholar]
- Rodríguez Zaurín, J., Tadhunter, C. N., Rupke, D. S. N., et al. 2014, A&A, 571, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rupke, D. S. N., & Veilleux, S. 2013, ApJ, 768, 75 [Google Scholar]
- Rupke, D. S., Coil, A., Geach, J. E., et al. 2019, Nature, 574, 643 [NASA ADS] [CrossRef] [Google Scholar]
- Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 [Google Scholar]
- Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Data Analysis Software and Systems IV, eds. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, ASP Conf. Ser., 77, 433 [Google Scholar]
- Shimwell, T. W., Röttgering, H. J., Best, P. N., et al. 2017, A&A, 598 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shimwell, T. W., Tasse, C., Hardcastle, M. J., et al. 2019, A&A, 622, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shimwell, T. W., Hardcastle, M. J., Tasse, C., et al. 2022, A&A, 659, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shlosman, I., Begelman, M. C., & Frank, J. 1990, Nature, 345, 679 [Google Scholar]
- Silk, J., & Rees, M. J. 1998, A&A, 331, L1 [NASA ADS] [Google Scholar]
- Spence, R. A., Zaurín, J. R., Tadhunter, C. N., et al. 2016, MNRAS, 459, L16 [NASA ADS] [CrossRef] [Google Scholar]
- Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 776 [Google Scholar]
- Tadhunter, C., Holt, J., Delgado, R. G., et al. 2011, MNRAS, 412, 960 [NASA ADS] [Google Scholar]
- Tamhane, P., Wadadekar, Y., Basu, A., et al. 2015, MNRAS, 453, 2438 [Google Scholar]
- Tasse, C., Shimwell, T., Hardcastle, M. J., et al. 2020, A&A, 648, A1 [Google Scholar]
- Ulvestad, J. S., Wrobel, J. M., & Carilli, C. L. 1999, ApJ, 516, 127 [NASA ADS] [CrossRef] [Google Scholar]
- van Cappellen, W. A., Oosterloo, T. A., Verheijen, M. A. W., et al. 2022, A&A, 658, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Weeren, R. J., Williams, W. L., Hardcastle, M. J., et al. 2016, ApJS, 223, 2 [Google Scholar]
- van Weeren, R. J., Shimwell, T. W., Botteon, A., et al. 2021, A&A, 651, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vardoulaki, E., Charmandaris, V., Murphy, E. J., et al. 2015, A&A, 574, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Varenius, E., Conway, J. E., Martí-Vidal, I., et al. 2015, A&A, 574, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Varenius, E., Conway, J. E., Martí-Vidal, I., et al. 2016, A&A, 593, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 769 [NASA ADS] [CrossRef] [Google Scholar]
- Veilleux, S., Rupke, D. S., Kim, D. C., et al. 2009, ApJS, 182, 628 [NASA ADS] [CrossRef] [Google Scholar]
- Vollmer, B., Braine, J., & Soida, M. 2012, A&A, 547, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Williams, W. L., van Weeren, R. J., Röttgering, H. J., et al. 2016, MNRAS, 460, 2385 [Google Scholar]
- Xia, X. Y., Xue, S. J., Mao, S., et al. 2001, ApJ, 564, 196 [Google Scholar]
- Yuan, T. T., Kewley, L. J., & Sanders, D. B. 2010, ApJ, 709, 884 [NASA ADS] [CrossRef] [Google Scholar]
- Yun, M. S., Hibbard, J. E., Condon, J. J., & Reddy, N. 1999, A&AS, 266, 29 [Google Scholar]
- Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 803 [Google Scholar]
Appendix A: Continuum images
![]() |
Fig. A.1. Continuum images of Mrk 273. Images in the first and second row show the large-scale structure at (a) 54 MHz with LOFAR LBA, (b) 1360 MHz with Apertif, (c) 1455 MHz with VLA and (d) 4700 MHz with VLA. The third and fourth row show small-scale structure in the nuclear region with (e) LOFAR at 144 MHz, VLA at (f) 4700 MHz and (g) 8415 MHz. The contour levels in all images are 3 |
![]() |
Fig. A.1. continued. |
All Tables
All Figures
![]() |
Fig. 1. Multi-wavelength emission from Mrk 273. The blue contours map the 144 MHz emission from the LOFAR 144 MHz image, red contours map the soft X-ray emission (0.4−2 keV) from L19, overlaid on the optical continuum image from S16 at λ = 6572 Å. The optical continuum tail seen towards the south of the nucleus is marked. The extended soft X-ray nebulae seen in south-west (SW) and north-east (NE) can be seen as well. The Chandra X-ray image has been smoothed with a Gaussian kernel of σ = 10 pixels, using the CIAO package. |
In the text |
![]() |
Fig. 2. LOFAR 144 MHz images showing emission at different scales in Mrk 273. Left panel: LoTSS-DR2 image at 30″ resolution shows the low surface brightness large-scale emission, with the northern arc marked. Black cross in the northern arc marks the nearby galaxy 2MFGC 11079 at that location. Middle panel: LoTSS-DR2 image at 6″ resolution, showing the southern arc and ridge in Mrk 273. Right panel: sub-arcsecond image of the nucleus of Mrk 273 at a resolution of 0.3″ × 0.2″. The northern (N), south-eastern (SE), south-western (SW) and diffuse north-east (NE) components are marked. The contour levels in all images are 3 |
In the text |
![]() |
Fig. 3. Integrated flux densities of Mrk 273 from our low resolution images, measured using 3σ contours for reference and NED. The red points show measurements from this study and orange points mark the NED values. |
In the text |
![]() |
Fig. 4. Spectral index maps of large scale emission in Mrk 273. (a) Spectral index map from 144−1455 MHz made with the LOFAR and VLA images at 16.8″ × 13″ resolution. The local rms noise, σrms, in the 144 MHz image is 200 μJy beam−1 and in the 1455 MHz image is 60 μJy beam−1. The contours mark the 1455 MHz emission at levels (−3, 3, 5, 10, 15, 25)×σrms. The map shows a spectral gradient across the arc, with α ∼ 0.6 at the southern outer edge and ∼0.9 at the northern inner edge. The ridge has a steep spectral index of α ≈ 0.88 ± 0.07. (b) The 20″ spectral index map from 144−1360 MHz made with the LOFAR and Apertif images. The local rms noise, σrms, in the 144 MHz image is 100 μJy beam−1 and in the 1360 MHz image is 37 μJy beam−1. The contours mark the 1360 MHz emission at levels (−3, 3, 5, 10, 25, 100)×σrms′. The pixels outside the 3σ contours show the lower limit of the spectral index. The black cross marks the background galaxy present at the location of the northern arc. The inset maps show the spectral index errors. (a) Spectral index map |
In the text |
![]() |
Fig. 5. JP model fits to the spectrum of the regions marked in Fig. 7. |
In the text |
![]() |
Fig. 6. 20″ spectral index map from 1455−4700 MHz made with the VLA images. The local rms noise, σrms, in the 1455 MHz image is 80 μJy beam−1 and in the 4700 MHz image is 40 μJy beam−1. The contours mark the 4700 MHz emission at levels (3, 5, 10, 15, 20, 50)×σrms. The inset map shows the spectral index error. |
In the text |
![]() |
Fig. 7. VLA 1455 MHz image of Mrk 273 at a resolution of 20″ × 20″. The solid lines mark the positive contour levels at: 3 |
In the text |
![]() |
Fig. 8. Homogeneous FFA model best fit for the N component flux densities. A turnover can be seen at ∼370 MHz. |
In the text |
![]() |
Fig. A.1. Continuum images of Mrk 273. Images in the first and second row show the large-scale structure at (a) 54 MHz with LOFAR LBA, (b) 1360 MHz with Apertif, (c) 1455 MHz with VLA and (d) 4700 MHz with VLA. The third and fourth row show small-scale structure in the nuclear region with (e) LOFAR at 144 MHz, VLA at (f) 4700 MHz and (g) 8415 MHz. The contour levels in all images are 3 |
In the text |
![]() |
Fig. A.1. continued. |
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.