Free Access
Issue
A&A
Volume 647, March 2021
Article Number A63
Number of page(s) 15
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202037677
Published online 11 March 2021

© ESO 2021

1. Introduction

Active galactic nuclei (AGNs) produce tremendous amounts of energy through the accretion of matter onto the super-massive black hole (SMBH). The release of this energy can heat and expel gas that would otherwise be available for star formation and for the growth of the SMBH. Thus the energy release can have a direct impact on the interstellar medium (ISM) of the host galaxy and this interplay, or feedback, links the evolution of the galaxy to the activity of the SMBH (e.g., Silk & Rees 1998; Di Matteo et al. 2005; Croton et al. 2006; McNamara & Nulsen 2007; Heckman & Best 2014; Harrison 2017).

Among the clear observational pieces of evidence of this interaction are outflows of gaseous matter observed in all phases of the ISM. Outflows of hot and warm ionized gas have been found in many AGNs over a range of redshifts suggesting that they could be a general feature (e.g., Veilleux et al. 2005; Tadhunter 2008; King & Pounds 2015; Morganti 2017).

In addition, several studies have shown outflows of matter from the cold ISM, which consists of molecular and atomic gas. Interestingly, the outflows of cold gas appear to carry higher mass outflow rates than their counterparts of ionized gas (e.g., Morganti et al. 2005a; Cicone et al. 2014; Morganti 2017; Fiore et al. 2017; Veilleux et al. 2020). There are different possible driving mechanisms capable of producing these outflows: radiative winds from the accretion disk, radiation pressure on dust, and radio jets. Observational evidence has been found for all of them with one being more likely than another depending on the type of host galaxy and nuclear activity mode (e.g., McNamara & Nulsen 2012; King & Pounds 2015; Morganti 2017).

In a number of cases, radio jets appear to be the most likely mechanism capable of accelerating the gas and producing outflows (e.g., Morganti et al. 1998, 2005a,b, 2013, 2016; Oosterloo et al. 2000; Mahony et al. 2013; Geréb et al. 2015; Aditya & Kanekar 2018; Aditya 2019). In particular, one of the findings supporting this is the higher fraction of outflows observed in young and restarted AGNs and, in some cases, the close morphological association of the region of the outflow and the radio jet (e.g., Holt et al. 2008).

Radio galaxies in their young phase are often characterised as compact steep spectrum (CSS) or GHz-peaked spectrum (GPS) sources based on the shape of the radio spectrum (e.g., O’Dea 1998; Orienti 2016). The radio continuum emission of these sources extends to a few kiloparsecs or less and it is still contained within the host galaxy. They are often hosted in galaxies with a richer ISM (Callingham et al. 2017; Holt 2009), which has also been confirmed by the study of H I gas (Geréb et al. 2015; Gupta et al. 2006; Chandola et al. 2011).

In restarted AGNs, the CSS or GPS source stems from the most recent or current cycle of AGN activity and is embedded in larger radio galaxies where the large-scale radio emission reflects previous cycles of AGN activity. The expected coupling of the jet with the ISM depends on the properties of the gas. In particular, in recent numerical simulations (Wagner & Bicknell 2011; Wagner et al. 2012; Mukherjee et al. 2016, 2018a; Cielo et al. 2018) the ISM is now modeled much more realistically (e.g., clumpy instead of smooth). An important result is that a clumpy ISM has a much stronger impact on the propagation of newly (re-)born jets: the jet interacts with dense clumps of the ISM, temporarily blocking its progress, which results in an over-pressured cocoon in the ISM that affects a much larger region of the host galaxy than the jet itself. Investigating whether this scenario is confirmed by the observations requires a high spatial resolution to resolve the distribution and kinematics of the outflowing gas and to compare it with the predictions of the simulations. Such a comparison has already been successful in some cases, for example, Mukherjee et al. (2018a) and Zovaro et al. (2019), but the number is still limited.

In this work, we focus on deriving the properties and structure of outflows of neutral atomic hydrogen (H I) gas observed in absorption in young or restarted radio galaxies. In order to resolve the outflows on sub-arcsecond angular resolution, we carried out Very Long Baseline Interferometry (VLBI) observations. This paper is a continuation of the work presented in Morganti et al. (2013) on 4C 12.50 and in Schulz et al. (2018) on 3C 236, both of which are restarted radio sources. Here, we present results from the remaining two sources of our initial sample, 4C 52.37 and 3C 293, and compare the properties of the H I outflows in all four sources in particular with respect to the evolution of the radio AGNs.

The paper is structured as follows: we first describe the sample and in particular the properties of 4C 52.37 and 3C 293. This is followed by a description of the VLBI observation and data reduction in Sect. 3. This is followed by the presentation of our results on the line and continuum in Sect. 4. We discuss in Sect. 5 the results of the two new objects and we also discuss some general results obtained by combining the findings for the full sample. We conclude this paper with a summary in Sect. 6 and provide additional plots in Appendix A. Throughout this paper, we use a ΛCDM cosmology of H0 = 70 km s−1 Mpc−1, Ωm = 0.3 and ΩΛ = 0.7.

2. Sample

The four sources 4C 12.50, 3C 236, 4C 52.37, and 3C 293 have been selected because they exhibit fast H I outflows, they are powerful radio sources, and they are likely to be young or restarted AGNs, but they are potentially at different stages in their AGN evolution. Based on their optical properties, 4C 52.37, 3C 236, and 3C 293 have been classified as low-excitation radio galaxies (Buttiglione et al. 2010; Best & Heckman 2012; de Gasperin et al. 2011), which makes the radio jet the most likely driving force behind the outflowing gas. Even though 4C 12.50 is classified as a high-excitation radio galaxy (Holt et al. 2011), Morganti et al. (2013) concluded that the H I outflow is jet-driven based on the location of the H I outflow. The VLBI observations revealed all of the H I outflow concentrated in a slightly extended cloud toward the hot spot of the southern radio jet. In 3C 236, part of the H I outflow was detected in the form of unresolved clouds towards the center of the AGN and as an extended cloud located towards the hot spot of the southeast radio lobe, but the VLBI data did not recover all of the outflow and suggested that part of the outflow might be in the form of a diffuse component.

Galaxy 4C 52.37 is a radio galaxy located at a redshift of 0.106 based on the Sloan Digital Sky Survey (SDSS; Abolfathi et al. 2018). The SDSS redshift corresponds to a systemic optical velocity of 31688 km s−1 and a linear scale of 1 mas ≈ 1.9 pc. The optical spectrum is characterized by strong broad permitted and forbidden emission lines. 4C 52.37 is included in the Compact Radio sources at Low Redshift (CORALZ) sample of compact and young radio sources (Snellen et al. 2004). The sub-arcsecond radio morphology has been classified as a compact symmetric object by previous VLBI observations (de Vries et al. 2009). So far, it is not clear whether there have been previous cycles of radio activity. The main H I absorption feature was first reported by Chandola et al. (2011) using the Giant Metrewave Radio Telescope (GMRT). The H I outflow was discovered by observations with the Westerbork Synthesis Radio Telescope (WSRT, Geréb et al. 2015; Maccagni et al. 2017) and covers more than 600 km s−1. For 4C 52.37, this study represents the first detailed investigation of the H I absorption with VLBI.

The radio galaxy 3C 293 is located at a redshift of 0.045 (Abolfathi et al. 2018). At this redshift an angular scale of 1 mas corresponds to 0.96 pc. Optical observations show a complex system of dust lanes and a companion galaxy (Heckman et al. 1986; Martel et al. 1999; Evans et al. 1999). The host galaxy of 3C 293 has undergone at least one merger event, but it is unclear whether the companion galaxy is connected to this (Evans et al. 1999; Capetti et al. 2002; Floyd et al. 2006). Star formation and radio activity are considered to be linked to the merger event (Tadhunter et al. 2005; Labiano et al. 2014). The large-scale radio emission extends over 2′ with a bright central region comprising two strong features extending about 2″ in size and oriented in an east-west direction (e.g., Haschick & Baan 1985; Beswick et al. 2002). Strong H I absorption has been detected from both features of the central region (Baan & Haschick 1981; Haschick & Baan 1985; Beswick et al. 2002) and a shallow, broad (∼1200 km s−1) blue wing corresponding to an H I outflow towards the western feature (Morganti et al. 2003, 2005a; Mahony et al. 2013). VLBI observations by Beswick et al. (2004) reveal a complex, extended radio morphology of the central region with H I absorption corresponding to the dust lane detected over most of the radio emission. However, the available bandwidth of the observations by Beswick et al. (2004) was insufficient to detect the outflow. Cold molecular gas in the form of CO has been detected in absorption and emission (Evans et al. 1999; Labiano et al. 2014) and associated with an asymmetric, warped disk rotating around the AGN. Labiano et al. (2014) did not detect a high-velocity molecular outflow associated with 3C 293 based on their CO observation.

3. Observation and data reduction

3.1. 4C 52.37

Galaxy 4C 52.37 was observed with a global VLBI array comprising the European VLBI Network (EVN), the Very Long Baseline Array (VLBA), and the Arecibo radio telescope on 2015 Oct 15 (project code: GN002C). The data were correlated at the Joint Institute for VLBI ERIC (JIVE) producing two data sets: a ‘continuum pass’ with four intermediate frequencies (IFs) each with 32 channels and 16 MHz bandwidth and a ‘spectral-line pass’ with one IF of 16 MHz bandwidth and 512 channels (see Table 1). We used J1545+5400 and J1638+5720 as the phase reference source and fringe finder (and bandpass) calibrator, respectively.

Table 1.

Properties of observations.

The data reduction process is similar to Schulz et al. (2018). We processed both the continuum and the spectral-line pass in the Astronomical Image Processing System (AIPS, version 31DEC16, Associated Universities, Inc. 1999), using the a-priori amplitude calibration and flag tables provided by the EVN calibration pipeline. We first started with the calibration of the continuum pass. We flagged the data if the station elevation was below 15°. Next, we corrected for the instrumental delay using the task FRING and J1638+5720. Due to the nature of the VLBI array, a single scan of J1638+5720 from all stations was not available. Therefore, we used two scans for the manual phase calibration. In the next step, we corrected for the phase, delay, and rate of both calibrators using a global fringe fit. The solutions for each calibrator were applied to themselves. In addition, the solution of the phase calibrator J1545+5400 was applied to our target source (4C 52.37). Finally, we performed the bandpass calibration using J1638+5720 and applied the solutions to all sources. Because the global fringe fit assumes a point source model, we imaged the calibrators and found them to be slightly extended. Therefore, we used the image of the calibrators as an input model for the global fringe fit and repeated the process up to the bandpass calibration.

For the spectral-line pass, we used the amplitude calibration provided by the EVN calibration pipeline and the calibration of the continuum pass. We first applied the solutions from the manual phase calibration of the continuum pass to the spectral-line pass. Afterwards, we took the solutions of the calibration from the continuum pass and applied them to the calibrator and target source of the spectral-line pass. In the last step, we performed the bandpass calibration. Afterwards, we produced a continuum and a spectral-line visibility data file of 4C 52.37. For the continuum data, we averaged all channels together, while the spectra-line data contains the full cube.

Flagging bad visibilities, imaging and self-calibration of the continuum data of 4C 52.37 was performed with DIFMAP (Shepherd et al. 1994; Shepherd 2011). We used natural weighting and set CLEAN windows to allow for better convergence of the CLEAN algorithm. Once a sufficiently good model was achieved, we performed phase self-calibration and imaged the data again. This process was repeated with smaller time intervals for the phase self-calibration as long as the solutions stayed smooth. Then, a constant gain correction factor over the entire observing time was determined with amplitude self-calibration and applied to the data. After flagging bad solutions, we repeated the smaller cycle of imaging and phase-self-calibration and performed another amplitude self-calibration with a smaller time interval. We repeated this process with decreasing time intervals for amplitude self-calibration. The noise level of the final continuum image is above the theoretical thermal noise level that is expected in phase referencing VLBI. Also, radio frequency interference in this frequency band can contribute to an increase in noise level.

For further processing of the spectral-line visibility data of 4C 52.37, we first inspected the data in detail and flagged bad visibilities. In order to keep the visibilities as uniform as possible across channels, we generally flagged entire scans or channels of a specific telescope instead of single visibilities. Next, we loaded the continuum image into AIPS and used it to perform a single phase self-calibration of the spectral-line visibility data. To improve sensitivity and consistency with the continuum data, we time-averaged data within 30 s. In the next step, we subtracted the continuum in the visibility domain using a fit to the first and last 100 channels. Because we are primarily interested in the broad absorption and in order to increase the sensitivity, we averaged the spectral-line data also in frequency. This reduced the velocity resolution by a factor of 3 to 20.7 km s−1. The resulting data cube was manually imaged in AIPS using robust weighting 1, a visibility taper of 0.3 at 10 Mλ, and a restoring beam of 20 mas. Like in Schulz et al. (2018) we found this combination of parameters to be sufficient. We only cleaned those channels of the cube that showed significant absorption in the area of the radio continuum. The imaged cube was Hanning-smoothed to further improve the sensitivity which reduces the effective velocity resolution by an additional factor of 2. The cube was de-redshifted in order to gain relative velocity information in the source rest frame. The full list of image parameters is given in Table 2.

Table 2.

Properties of images.

3.2. 3C 293

Galaxy 3C 293 was observed with a global VLBI array on 2015 March 2 (project code: GN002A). J1350+3034 and J1407+2827 served as phase reference source and bandpass calibrator, respectively. Correlation and data processing were done following the procedure outlined for 4C 52.37.

Unfortunately, an issue in the frequency setup of the observation resulted in a shift in the center frequency. As a result a significant part of the blueshifted wing observed in previous observations was shifted out of the band. Therefore, only a linear fit using the outer 200 channels on the high-velocity side of the spectrum could be used for the subtraction of the continuum in the visibility domain. We used a restoring beam of 25 mas for the continuum image and spectra-line cube instead of the highest angular resolution due to the large extent of the radio emission.

3.3. Archival radio data

We make use of archival Very Large Array (VLA) and WSRT data for all four sources in this paper. The VLA data of 3C 236 and 3C 293 were presented in Schulz et al. (2018) and Mahony et al. (2013), while the WSRT data on 4C 12.50 and 4C 52.37 were taken from Morganti et al. (2013) and Geréb et al. (2015). We also made use of VLBI data from 4C 12.50 (Morganti et al. 2013) and 3C 236 (Schulz et al. 2018).

In Appendix A we present the H I profiles of these objects, which have been fitted with Gaussian distribution (using the Python package LMFIT; Newville et al. 2014) to characterize the absorption spectrum. The fit parameters are listed in Table A.1 and plots of the fits are shown in Figs. A.2A.4. These parameters are slightly different from Maccagni et al. (2017) because that study used the busy function (Westmeier et al. 2014) to fit the spectrum. In order to be consistent, we use Gaussian distribution for all our objects for an easy comparison.

4. Results

4.1. 4C 52.37

As we will describe in the sections below, in 4C 52.37 we have succeeded in imaging the radio continuum on scales between ∼10 and 40 pc, and we have recovered a large fraction of the H I absorption and, in particular, most of the outflowing gas. The outflow appears to be located mostly in the central regions, specifically in the central ∼120 pc.

4.1.1. VLBI continuum

The top panel of Fig. 1 shows the continuum VLBI image of 4C 52.37 at a resolution of 20 mas and a zoom-in into the central 200 pc at 6 mas. At the higher angular resolution, the central region is clearly resolved into two distinct features, E1 and W1. They are located about 20 mas apart and are separated by a small drop in surface brightness. The western feature W1 represents the peak of the radio emission and is compact while the eastern one (E1) is marginally resolved at 6 mas resolution. Additionally, there are two prominent radio lobes, E2 and W2. The overall radio morphology is asymmetric and spans a projected linear size of ∼280 mas, which corresponds to ∼530 pc.

thumbnail Fig. 1.

Continuum images of 4C 52.37 (top) and 3C 293 (bottom) convolved with a synthesized beam of 20 mas and 25 mas, respectively, which is marked by the gray circle. The top panel also shows a higher resolution (6 mas) zoom-in of the central region of the radio source. The contour lines begin at three times the noise level of the corresponding image and increase logarithmically by a factor of 2. The solid, horizontal scale bar represents the spatial scales.

The majority of the integrated flux density stems from extended radio emission of W2 and E2 with W2 accounting for almost 50%. The central features E1 and W1 contain about 25% of the total flux density. In total, we measure a flux density of 429 mJy, which is about 20% below the value obtained, at lower resolution, by Geréb et al. (2015) and Maccagni et al. (2017) with the WSRT (577 mJy) and from the VLA survey FIRST at 1.4 GHz (575 mJy, Helfand et al. 2015). Even considering our conservative flux density uncertainty of 15%, this suggests that there is a small amount of extended low-surface brightness emission that is either below our sensitivity or resolved out or both. The shortest baseline of the VLBI array, which is about 460 Mλ, implies that extended emission larger than about 550 mas cannot be recovered.

The EVN observations at 1.659 GHz by de Vries et al. (2009) provide the best comparison on VLBI scales to our observations in both frequency and resolution. Our measured VLBI flux density is consistent, within the uncertainty, with the value obtained by de Vries et al. (2009) of 448.3 mJy ± 1.1 mJy. Also, the morphology of the source recovered by these two observations matches very well.

4.1.2. VLBI H I absorption

The H I absorption profile from the archival WSRT observations (Geréb et al. 2015) and the integrated spectrum from the new VLBI H I are shown in Fig. 2. Because of the lower spatial resolution of the WSRT observation, we can use this absorption spectrum as a reference of the total absorption that should have been recovered. The WSRT spectrum comprises of two main features: a narrow, deep absorption close to the systemic velocity and a broad, shallower feature blueshifted with respect to the systemic velocity. The latter has been interpreted in previous studies as outflowing H I gas (Geréb et al. 2015; Maccagni et al. 2017).

thumbnail Fig. 2.

H I absorption spectra of 4C 52.37. The unresolved WSRT spectrum is shown in red. The two spatially integrated VLBI profiles (in blue) were extracted from a region of 60 mas × 40 mas centered on E1. The dashed blue line corresponds to the VLBI spectrum that was obtained by integrating only pixels with a value of ≤ − 3σVLBI, cube, while for the solid blue line the integration was only limited to the region marked by the 3σVLBI, cont contour.

The spatially integrated VLBI H I absorption spectra have been extracted using the same method as in Schulz et al. (2018). First, we created a mask based on the continuum image using only those pixels with a value above three times the continuum noise level (σVLBI, cont). Second, for each channel we applied the mask and integrated either all remaining pixels of the channel (solid, blue line) or integrated only the remaining pixels with a value above three times the noise level of the cube (σVLBI, cube) leading to the dashed, blue line. The clipped and non-clipped spectra shown in Fig. 2 have been obtained by integrating the absorption in the central 60 mas × 40 mas around E1. This was done to prevent the spectrum from being dominated by the noise and because the majority of the absorption is concentrated in this region. This would make a comparison between the clipped and non-clipped spectrum meaningless because the latter would have been largely affected by noise. Therefore, we extracted the integrated VLBI spectrum around E1 (also see Sect. 4.1.3).

A comparison between the WSRT and the clipped VLBI spectrum shows (see Fig. 2) that our VLBI observation recovers a large fraction of the H I broad, blueshifted wing in this region. The narrow feature close to the systemic velocity shows a greater discrepancy between WSRT and VLBI. In the following sections, we argue that the non-detected absorption is largely due to low-column density gas across the source.

As mentioned in Sect. 3.3, we fitted Gaussian distributions using the Python package LMFIT (Newville et al. 2014). We find that both features of the WSRT absorption spectrum are well described by two Gaussian distributions. The fit parameters are listed in Table A.1 and plots of the fits are shown in Fig. A.1. These parameters are slightly different from Maccagni et al. (2017) because that study used the busy function (Westmeier et al. 2014) to fit the spectrum. However, the busy function is not suitable to fit the spectra of other sources in our sample because the spectra consist of multiple components. In the WSRT spectrum, the broad feature is blueshifted by ∼310 km s−1 and it has a full-width at half maximum (FWHM) of about 330 km s−1.

4.1.3. H I gas distribution

Most interestingly, in 4C 52.37 we have been able to trace the location of the outflowing H I gas. The central panel of Fig. 3 shows the integrated optical depth of the H I gas, ∫τdv, detected by our VLBI observation with respect to the continuum source. To calculate ∫τdv, we used the same clipping method as described in the previous section, which is a 3σ cut in the continuum and in the cube. The map shows a large concentration of H I gas toward the central region, over an elongated region of ∼140 pc, partly overlapping with the eastern component E1, but not W1. In addition, there are smaller clouds of H I gas distributed over the extent of the continuum emission.

thumbnail Fig. 3.

Top panel: position-velocity (PV) diagram of 4C 52.37 along a cut aligned with the largest extent of the central absorption feature. The solid, horizontal line corresponds to the systemic velocity. The dashed, vertical lines mark the size of the synthesized beam. Contour lines correspond to 3 (grey, solid), −2 (black, dashed), −3, and −5 (black, solid) times the noise level in the cube. Middle panel: contour lines show the continuum radio emission as in Fig. 1 with the color map corresponding to N H I T spin 1 $ N_{{\mathrm{H}\,\textsc{I}}} T^{-1}_{\mathrm{spin}} $. The dashed lines represent the cuts for the PV diagram. Bottom panel: same as top panel, but a cut along the position angle of W1 and E1. The top and bottom panels are centred on E1.

In order to characterize the H I gas further, we extracted position-velocity diagrams along different position angles. Figure 3 shows slices along the extent of the H I gas in the central region (panel a) and along the position angle of the radio emission in the central region (panel c). The position-velocity plots show that the gas in the central region covers the full range of velocities seen in the spatially integrated absorption spectrum of Fig. 2. Three features, labeled S1, S2, and S3, can be seen at increasing (blueshifted) velocity from the systemic velocity. All three features are located within the central 100 pc of the source. While S2 and S3 are unresolved, S1, which corresponds to the gas close to the systemic velocity, is extended (Fig. 3a). From Fig. 3 we can therefore derive the location of the outflowing H I. The gas does not cover the entire central region but it is concentrated towards E1. Since E1 and W1 are of comparable brightness, it would have been possible to detect against the latter a more extended distribution of gas of a similar column density if it were present in front of the continuum.

Figure 2 shows that there is a significant fraction of absorption missing close to the systemic velocity in the VLBI extracted around E1, in the velocity range of S1. The patches of absorption that are depicted in Fig. 3b cannot account for the undetected gas. Because they are weak absorption features, we are cautious in interpreting them, but we might be seeing the presence of regular rotating gas distributed across the radio source. However, we do not reach the optical depth across the radio source to fully recover the gas at the systemic velocity. In comparison, most of the absorption in the velocity range of S2 and S3, which is the blueshifted component, is recovered in the region around E1. Therefore, clouds of outflowing gas of a similar optical depth to S2 and S3 are unlikely to exist anywhere else toward the radio source and the undetected absorption at the systemic velocity must be in the form of compact clouds or diffuse gas of a lower optical depth. The only possible exception may be the small absorption feature (S4) at the edge of E2 (see Fig. 3). This absorption is also blueshifted with respect to the systemic velocity by about 300 km s−1. S4 is much fainter than S2 and S3 and located toward the edge of the radio continuum of E2 rather than the brighter parts. That is why we are cautious in interpreting this feature, but it could be a sign of outflowing clouds outside the central region.

The peak brightness of the continuum and the noise level of the cube imply an optical depth limit of about 0.035 at 3σVLBI, cube (∼0.011 at 1σVLBI, cube). The optical depth limit toward the peak of the western radio lobe (W2) is only slightly higher. W2 has a peak brightness of about 42 mJy beam−1, which corresponds to an optical depth limit of 0.042 at 3σVLBI, cube. At the location of E1 the peak optical depth of the absorption feature is about 0.06. Therefore, we would have been able to detect clouds such as S1 at the peak of the western radio lobe (W2) unless the gas clouds have a lower optical depth or the gas is diffuse. This is an indication that the missing H I absorption is in part due to low-column density gas that is distributed across the continuum source (see Sect. 4.1).

4.2. 3C 293

Given the much larger spatial scale, we managed to recover only a relatively small part of the continuum structure of 3C 293. Similarly, a small fraction of the H I absorption is also recovered. This is particularly the case for the H I outflow. However, the data show that at least a small fraction of the outflow can be seen against the central (∼50 pc) region.

4.2.1. VLBI continuum

The bottom panel of Fig. 1 shows the VLBI continuum emission recovered by our observations. The radio emission from the eastern and western structure spans over 2″ as seen in other high-resolution radio images (e.g., Beswick et al. 2004; Mahony et al. 2013). The total flux density from our image is about 539 mJy, which is about 17% of the total flux density recovered by Beswick et al. (2004) at a resolution of 30 mas. The difference can be best explained by the spatial filtering. The shortest baseline of our VLBI array is 400 Mλ, which means that any extended emission on scales larger than about 630 mas is resolved out. In contrast, Beswick et al. (2004) combined a global VLBI array similar to ours with the Multi-Element Radio Linked Interferometer Network (MERLIN) and VLA data. This has provided significantly shorter baselines, which improves the sensitivity for diffuse extended emission.

The brightest feature in our image matches the location of the VLBI core identified in Beswick et al. (2004). Interestingly, we find the peak brightness to be a factor of 4.1 brighter than in Beswick et al. (2004). This cannot be explained by the difference in resolution, because the synthesized beam is similar: 25 mas in our case, compared to 30 mas in Beswick et al. (2004). A frequency-dependent effect can also be excluded, because the frequency setup of both studies is chosen to investigate H I absorption. This suggests that the emission from the radio core is variable and the VLBI core is brighter in our observation compared to previous observations. Our data do not allow us to test whether the radio emission extending from the core also shows changes in brightness or distribution. However, this result warrants caution when comparing observations.

4.2.2. VLBI H I absorption

Figure 4 shows the H I absorption spectrum from the VLBI observations over-plotted to the spectrum obtained with the VLA with 1″ resolution by Mahony et al. (2013). At this lower resolution, the structure of the source is made of two – barely resolved – components, the eastern and western radio lobes. Mahony et al. (2013) extracted H I absorption spectra for the eastern and western part, separately. For comparison, we provide spatially integrated H I absorption spectra from our VLBI data for the eastern and western part.

thumbnail Fig. 4.

H I absorption spectra of 3C 293. The low-resolution VLA spectrum is shown in red. The solid blue and dashed violet lines correspond to the spatially integrated VLBI spectra of the eastern and western part of the source. The dotted, gray line represents the total integrated VLBI spectrum. The VLBI spectra were obtained by only integrating pixels with a value of ≤ − 3σVLBI, cube.

Overall, we recover a relatively small fraction of the VLA-detected absorption. In particular, we recover absorption from the narrow, deep absorption on both lobes. We also marginally detect, on the western side of the source, clouds blueshifted by up to −300 km s−1 with respect to the systemic velocity (see solid blue line in Fig. 4). We do not consider any features below a relative velocity of −450 km s−1. The issue with the frequency setting of our observation mentioned in Sect. 3.2 strongly affected our ability to perform an optimal continuum subtraction at the corresponding frequencies. As a result, our effectively usable bandwidth is comparable to that of Beswick et al. (2004).

The spatial distribution of the H I absorption obtained by our VLBI observation is shown in Fig. 5. On the eastern side, we detect H I gas distributed in three connected regions. The kinematic of the gas in this part of the source is mostly consistent with quiescent gas like the one originating from the narrow, deep absorption feature. On the western side, the H I distribution is more patchy. Most of the broad absorption is located toward the core and shows kinematics deviating from that of the regularly rotating gas (as described in Beswick et al. 2004).

thumbnail Fig. 5.

Left panel: position-velocity diagram of a representative location in the eastern part and column density normalized by spin temperature (bottom panel) for the same velocity range as in the right panel. Right panel: position-velocity diagram of 3C 293 at the location of the VLBI core (top panel) and column density normalized by spin temperature separated into systemic and outflowing H I gas (bottom two panels). The dashed lines in the top two panels correspond to the size of the synthesized beam. The dashed lines in the bottom panels correspond to the location of the slice for the position-velocity diagram.

We identify two kinematic components (S1 and S2). S1 is likely to be part of the narrow absorption feature similar to the gas recovered on the eastern side and also traced by the observations of Beswick et al. (2004). S2 is broader in velocity than S1 and blueshifted by up to −300 km s−1 with respect to the systemic velocity. Interestingly, S2 appears to be broader than what was detected by Beswick et al. (2004) at the location of the core (see their Fig. 8). Thus, this suggests that the gas originating the S2 absorption has a highly disturbed kinematics and could be part, as suggested by Haschick & Baan (1985), of the broad H I outflow observed by Morganti et al. (2005a) and Mahony et al. (2013). If this is the case, in 3C 293 also the outflow can already be observed in the inner nuclear region (≤40 pc).

We also tentatively detect gas clouds of a similar velocity to S2 at other locations in the western part of the source (see Fig. 5). They are at the edges of the recovered radio emission, which is why we are cautious in interpreting them. However, if confirmed they would indicate that the outflow is extended, as suggested by Mahony et al. (2013), and clumpy.

The gas related to S1 has a peak optical depth of about 0.1 while the gas of S2 peaks at an optical depth of about 0.024. Only at the location of the core are we able to reach the necessary optical depth of about 0.005 and 0.014 at 1σVLBI, cube and 3σVLBI, cube, respectively, to detect S2. We cannot detect clouds such as S2 toward the eastern part of the source. The tentatively detected gas cloud at the edge of the western lobe has a peak optical depth of about 0.22. If such clouds existed toward other locations of the western radio lobe where the continuum emission is at the same level or higher, we would have been able to detect them.

5. Discussion

In both radio galaxies, 4C 52.37 and 3C 293, we were able to recover (at least partly) and trace the H I outflows known from previous observations on scales of about 20 pc. We discuss below their structure as revealed on a few tens of parsec scale, their properties, and their location. We then compare them with already published H I VLBI results (see Sect. 5.4).

5.1. Properties of the H I gas and the outflow in 4C 52.37

As described in Sect. 3.1, in 4C 52.37 we recover most of the H I outflow and only part of the H I close to the systemic velocity. The outflowing gas is mainly concentrated in the central region (see Fig. 3) and largely contained within two clouds, S2 and S3, detected toward E1 (see Fig. 1), the eastern component of the central region. The location of the core in 4C 52.37 is still unclear and this has implications for understanding the origin of the outflow. This will be further discussed in Sect. 5.2. However, it is reasonable to assume that the active black hole is located either in E1 or W1, which are separated by about 40 mas (∼76 pc) in projected size. Table 3 lists the properties of H I gas in 4C 52.37 based on our observations and using the same clipping as for the absorption profile. We determine the spin temperature (Tspin) and normalized column density (NH I) of the outflow by integrating over S2 and S3. In order to calculate any subsequent parameters, we assume Tspin = 1000 K, which is reasonable considering the proximity of the gas to the AGN. This suggests a column density of ∼2.0 × 1022 cm−2. For a spherical geometry, this yields an H I mass of 1.8 × 105M. Following Heckman (2002), we estimate the mass outflow rate to be ∼4.3 M yr−1 with

M ˙ H I 30 r kpc N H I 10 21 cm 2 v 300 km s 1 Ω 4 π M yr 1 $$ \begin{aligned} \dot{M}_{{\text{ H}}{\small {{\text{ I}}}}} \sim 30 \frac{r_\star }{\mathrm{kpc}}\frac{N_{{\text{ H}}{\small {{\text{ I}}}}}}{10^{21}\,\mathrm{cm^{-2}}}\frac{v}{300\,\mathrm{km\,s^{-1}}}\frac{\Omega }{4\pi }\,M_\odot \,\mathrm{yr}^{-1} \end{aligned} $$(1)

Table 3.

Properties of the kinematically disturbed H I gas.

where r is the deprojected distance of the cloud, v its velocity, and Ω its solid angle, which is assumed to be π (see Table 3). We take the synthesis beam as an upper limit on the radius of the clouds because S2 and S3 are both unresolved (see Fig. 3).

The gas detected with velocities close to the systemic (S1) represents a fraction of what has been detected by the WSRT. Assuming also for this gas that Tspin = 1000 K, given the proximity to the AGN, we estimate an average column density of 1.0 × 1022 cm−2. Based on our VLBI data, we describe S1 by an ellipsoid that is resolved in one direction (95 pc × ≲38 pc) and estimate a lower limit on the mass of the H I gas of ∼2.3 × 105M. In the kinematics of this gas, we do not detect signatures of rotation. However, we cannot fully rule out the presence of a circumnuclear H I disk, for example, given the large fraction of missing absorbed flux. Circumnuclear disks of H I gas have been traced in the center of a number of radio galaxies such as 3C 236 (Schilizzi et al. 2001; Tremblay et al. 2010), Cygnus A (Struve & Conway 2010), NGC 4261 (Jones et al. 2000), and Centaurus A (van Gorkom et al. 1990; Morganti et al. 2008).

Our observations are limited by the optical depth that we can probe. While we are only able to probe an optical depth similar to E1 toward the peak of the western lobe (W2), S1 is extended and detected against fainter parts around E1. If there is a disk made of clouds like S1, we should have detected it against W1. Therefore, it is possible that the missing absorption comes from an extended, diffuse component of the systemic gas or lower column density clouds. The patches of absorption observed across the radio source could be related to this.

If an H I disk is present, the absorption recovered in the central region of the radio source would correspond to gas in an inner denser disk, while the undetected absorption would be related to the extended (and lower column density) outer part of the disk. It is possible that the radio jets are drilling into this gas disk. This could explain the asymmetric morphology of the radio emission. Radio jets that are expanding into the ISM have been observed for example in IC 5063 (Oosterloo et al. 2017) and NGC 1167 (Murthy et al. 2019). In both objects VLBI observations were only able to recover a small fraction of the absorption.

5.2. Origin of the H I outflow in 4C 52.37

The fact that most of the H I outflow in 4C 52.37 is concentrated in a very small region in the center of the radio source is important for understanding the driving mechanism of the outflow and the evolution of the radio source. As mentioned in Sect. 1, 4C 52.37 is classified as a low-excitation radio galaxy. As such, the most likely drivers of the outflow are the radio jets rather than radiative winds from the AGN. It is worth taking a closer look at the region of E1 and W1 where the AGN is likely to be located.

We compare the location and spectral index of W1 and E1 (Fig. 1) using the Gaussian model components fitted to the visibility data at 1.659 GHz and 4.99 GHz by de Vries et al. (2009). We find that the position of the model components at 1.659 GHz and 4.99 GHz match for W1 and E1. We also fit Gaussian components to W1 and E1 of our VLBI observation and find no difference in location with respect to the components from de Vries et al. (2009) at 1.659 GHz. In addition, we determine the spectral index α between the two frequencies for W1 and E1 to be 0.17 and −0.16, respectively. Both features seem to have a flat spectral index (|α| < 0.5). Due to the time span between our observation and the data from de Vries et al. (2009), we do not determine the spectral index between these data sets. As the radio jets of 4C 52.37 seem to be oriented in the plane of the sky and we find H I absorption toward the region of E1 marginally extending in the direction of W1, it is possible that the spectral index measurements at least of E1 are affected by an external absorber. This would be similar to the case of, for example, NGC 1052 (Kameno et al. 2001; Kadler et al. 2004; Sawada-Satoh et al. 2008; Baczko et al. 2019) where free-free absorption is observed toward the origin of the radio jets. An absorber might also explain the small emission gap between E1 and W1 visible in the 6 mas-resolution image (Fig. 1). Our VLBI data and de Vries et al. (2009) show W1 to be the brightest feature, but this may not be the case intrinsically if free-free absorption affected our measurements of E1. However, free-free absorption is frequency dependent and as such we should observe a change in position of the brightest feature in the components from de Vries et al. (2009) as is observed in NGC 1052. Simultaneous multi-frequency VLBI observations are necessary to characterize the nature of the central region of 4C 52.37 in greater detail.

Galaxy 4C 52.37 has been included in the CORALZ sample (Snellen et al. 2004) of candidate young radio sources based on its CSS-like spectrum and is considered to be a young radio source. We can use the same approach as O’Dea et al. (2001) to estimate a lower limit on its age based on the size of the radio source. As the source extends further to the west, we take the western lobe as a reference estimating its size to be about 130 mas (250 pc). This yields a minimum age of the radio source of tmin ≈ 2 × 104 yr. Thus, the jet would not have had sufficient time to disperse more of the ISM gas. This is consistent with the location of the H I outflow toward E1 close to the center of the AGN. However, the fact that we do not recover all the H I gas at the systemic velocity suggests that this gas is distributed over a larger area where it has not yet been affected by the jet. This and the cloud-like nature of the outflow would be consistent with numerical simulations by Wagner et al. (2012) and Mukherjee et al. (2017); Mukherjee et al. (2018a). These simulations have also shown that the propagation of the jet into the ISM and the amount of gas that is affected by it is determined by gas density, jet power, and inclination angle.

5.3. H I outflow in 3C 293

In Sect. 4.2 we have shown that, unlike for 4C 52.37, our VLBI observation of 3C 293 recovers a relatively small fraction of the H I absorption and, in particular, of the broad, blueshifted component. We strongly detect signatures of this component against the VLBI core and we tentatively detect clouds against other parts of the western radio lobe. We are cautious in interpreting the latter features, which is why we only derive physical values for the absorption against the VLBI core. At this location, we detect unresolved clouds of outflowing H I gas with velocities blueshifted and deviating up to 300 km s−1 with respect to the systemic velocity (S2). The clouds have column densities of ∼4 × 1021 cm−2. If we take the beam as an upper limit on the size and assuming a spherical geometry, we derive an H I mass of the clouds of 1.5 × 104M and a mass outflow rate of ∼0.36 M yr−1.

Mahony et al. (2013) investigated the H I outflow using high-resolution VLA observations. The study concluded that the outflow is largely extended across the western lobe. The authors state that it is unlikely that all of the outflow is concentrated within the VLBI core because of the implied high optical depth and the location of the H I outflow in the VLA image. Our findings support and are complementary to these results. The mass outflow rate of S2 is a fraction of the total rate estimated by Mahony et al. (2013) to be 8 − 50 M yr−1. We recover only a fraction of the radio emission from the lobe and, at the high spatial resolution of our observations, the core appears to be the brightest component in the western lobe. This suggests we are limited in our sensitivity to trace the outflow in the western lobe. However, our observations show that a fraction of the outflow exists on nuclear scales close to the AGN. Furthermore, we also find clouds at other locations that could be part of the outflow (see bottom right panel of Fig. 5) indicating that the outflowing gas might be extended and clumpy. The optical depth of these tentative detections suggests that we should have been able to find absorption against other parts of the western radio continuum emission, if distributed in clouds with similar properties. The lack thereof suggests that the outflow might have a diffuse component.

The cold ISM gas in 3C 293 is not the only ISM phase that has an outflowing gas component: Emonts et al. (2005) and Mahony et al. (2016) reported outflowing ionized gas. The latter study showed that disturbed ionized gas is also observed for a few kiloparsecs in the direction perpendicular to the radio jet. This has been interpreted as the result of a cocoon of shocked gas produced by the radio jet, expanding perpendicular to the jet propagation into the gas disk aligned with the dust-lane known to be present in the host galaxy of 3C 293. Such cocoons have been predicted in numerical simulations describing the interaction between a newly formed radio jet and the surrounding, clumpy medium (e.g., Wagner & Bicknell 2011; Mukherjee et al. 2016, 2017).

Our results are consistent with this scenario, which predicts that the interaction between jet and ISM begins as soon as the jet starts expanding into the gas and it occurs along and perpendicular to the jet propagation. Additionally, the longer the jet expands into the ISM the more the gas can get dispersed. The detection of clouds of outflowing H I toward the core, in addition to outflowing gas distributed across the western lobe, further supports this scenario.

5.4. Properties of H I outflows in our sample

We now compare the results on the properties of the H I outflows presented here for 4C 52.37 and 3C 293 with the other two objects (4C 12.50 and 3C 236) presented in Morganti et al. (2013) and Schulz et al. (2018). These four radio galaxies are the only objects for which fast outflows in H I absorption on milliarcsecond scales have been studied.

We can summarize the main results we have obtained as:

  1. We recover at least part of the H I outflow on the VLBI (20 pc) scale in all objects.

  2. This indicates that all the outflows include at least a component of relatively compact clouds (with masses in the range of 104 − 105M).

  3. We suggest that the clouds might be embedded in a more diffuse component. We detect this in 4C 12.50 (Morganti et al. 2013). Also in 3C 236 we concluded that at least part of the undetected outflow is made of a diffuse, lower column density component (Schulz et al. 2018). This could also be the case in 3C 293 since the outflow must be extended over the western lobe.

  4. The data alone provide only limited additional insight into the understanding of the driving mechanism. As already mentioned in the introduction, we remark that the radio galaxies are low luminosity optical AGNs (with the exception of 4C 12.50), and are all young or restarted radio galaxies. Therefore, the radio jet is a very likely mechanism for driving the outflows. The properties of the outflows we derive are consistent with what is predicted by the simulation of jet-driven outflow. The derived density of the clouds is also similar to that used in the simulations (e.g., Wagner & Bicknell 2011; Wagner et al. 2012; Mukherjee et al. 2016, 2018b).

In all four sources we find evidence for a clumpy distribution of the outflowing H I but at different levels. This is shown in the spatial and velocity distribution of the gas as well as the difference between the blueshifted H I absorption recovered at low and high angular resolution. Here, we attempt to explore whether our results show an evolutionary sequence as expected from numerical simulations by considering the age of the sources and the characteristics of the H I outflow both in the integrated spectrum as well as in the spatially resolved distribution of the outflowing gas. Numerical simulations (see, e.g., Wagner & Bicknell 2011; Wagner et al. 2012) have shown that radio jets will accelerate and disperse gas along their direction of propagation, but also transversal to it when they expand into a clumpy medium. With time, more and more of the ISM gas will be dispersed (Mukherjee et al. 2018b) leading to changes in the properties of the ISM and outflowing gas.

In Sect. 5.2, we have estimated a minimum age of the radio source of 4C 52.37 of tmin ≈ 2 × 104 yr based on the size of the radio continuum emission. This is of the same order as the minimum age estimated for 4C 12.50 (Morganti et al. 2013). In 4C 12.50, all of the H I outflow has been recovered with VLBI in the form of a slightly extended gas cloud toward the southern jet, while in 4C 52.37 most of the outflow is recovered and located toward the inner region of the radio source. It is still mostly compact in 4C 52.37, but the undetected absorption in this source suggests that a fraction of the outflowing gas is located elsewhere. Based on this, 4C 12.50 would be consistent with an early evolutionary stage of the jet-ISM interaction and 4C 52.37 would be at a slightly more advanced stage compared to 4C 12.50.

The minimum age of 4C 12.50 and 4C 52.37 is an order of magnitude lower than the age estimated for 3C 236 (O’Dea et al. 2001; Tremblay et al. 2010). In 3C 236, only a fraction of the H I outflow was recovered, mostly toward the nuclear region of the radio source, but also toward the hot spot of the southeast jet in the form of an extended cloud. This would be consistent with a later stage in the evolution of the radio source where the jet had significant time to disperse most of the gas. In agreement with this, also in the case of 3C 293 we recover only a fraction of the outflow, mostly at the location of the core.

In order to characterize H I spectra in more detail, we fitted between one and three Gaussian distributions to the low and high resolution spectra. The spectra and fit parameters are provided in Table A.1, respectively. While Gaussian functions do not represent a physical model, they provide a simple and consistent approach. The fits to the VLA spectra of 3C 236 and 3C 293 are consistent with the fits presented in Labiano et al. (2013) and Mahony et al. (2013). In the case of 4C 12.50 we only performed a fit of the WSRT spectrum, because the VLBI spectrum is consistent with it as shown in Morganti et al. (2013). For 4C 52.37 our fit parametrization is different from Maccagni et al. (2017), because this study used the busy function of Westmeier et al. (2014) to fit the spectrum.

The shape of the blue wing in the absorption spectra differs greatly among the four sources. In the case of 4C 12.50 and 4C 52.37 it is more detached from the systemic gas than in 3C 236 and 3C 293. In the latter two sources, the blue wing is shallow and smooth, including a smooth transition from the systemic H I. The difference in the blue wing is reflected by the properties of the Gaussian components that describe that part of the spectrum. The FWHM of 3C 236 (VLA component 3) is larger by factor of 2 than 4C 12.50 (WSRT component 2) and 4C2.37 (WSRT component 3). The component describing the wing of 3C 293 (VLA component 2) has a FWHM similar to that of the latter two sources, but the uncertainty of the fit result is very high most likely due to correlation with the parameters of component 3. It is therefore possible that the shape of the blue wing in the H I spectrum is a signature of the compactness of the H I outflow.

It is clear that the actual situation is more complicated, because the above arguments do not consider individual differences in the gas distribution and evolution of the system, such as merger activity (3C 293) and previous cycles of activities (4C 12.50, 3C 236, and 3C 293). The ages of the radio sources provide only limited information as these are lower limits. Also, 4C 12.50 is the only source in our sample where no absorption is detected toward the VLBI core. Interestingly, the density of the detected outflowing gas clouds is relatively similar among the four sources and the VLBI observations probe similar spatial scales in all of them. Therefore, it is intriguing that a connection can be drawn in the form of an evolutionary scheme for these sources and it shows the need for further studies of H I outflows in other systems in combination with numerical simulations.

6. Summary and conclusion

In this paper, we have presented results, using global spectral-line VLBI, for two additional sources from our investigation of the parsec-scale H I outflow in four radio galaxies. In the case of 4C 52.37, this is the first such VLBI study of the H I gas distribution. We recovered the majority of the H I absorption related to the outflowing gas and a significant fraction of the systemic H I gas with a resolution of 20 mas. We detected the outflow in the form of two unresolved clouds toward the central region of the radio source. They have a combined mass outflow rate of less than 4.3 M yr−1 and cover about 600 km s−1. The recovered systemic gas is extended and largely concentrated toward the central 100 pc, but it also shows signs of being clumpy and distributed across the radio source. We find signs that the radio jet of 4C 52.37 expands into the gas disk, but further observations are necessary.

In the case of 3C 293, our VLBI observation resolved out most of the extended radio emission and it suffered from technical difficulties limiting the frequency range over which we could probe the broad H I outflow. Nevertheless, we recovered a fraction of the H I gas over a velocity range of 300 km s−1 in the form of distinct gas clumps. Some of the clumps are located toward the VLBI core, close to the AGN, and are likely signs of the outflow in this source. As both 4C 52.37 and 3C 293 are low-excitation galaxies, it is reasonable to assume that only the radio jets are capable of driving the H I outflows.

All four sources (4C 12.50, 4C 52.37, 3C 236, 3C 293) in our sample have in common the fact that the outflowing H I gas is clumpy and that some of the gas is located toward the innermost region (< 50 pc) of the radio source. For 4C 12.50 and 4C 52.37 we find that almost all of the outflowing H I gas is concentrated in a small number of compact clouds, whereas for the other two sources our data implies that the H I outflow is more extended and might have a diffuse component. Taking into account our observational limitations, we can interpret the properties of the H I outflow as a sign of different stages of jet-ISM interaction in these sources. This would be consistent with numerical simulations that have shown that the radio jets would disperse more and more of the ISM gas while they are expanding over time.

Our results highlight the importance of high-angular-resolution observations to understand the gas properties and dynamics in connection to AGN activity. Upcoming facilities, such as the Square Kilometre Array and its pathfinders such as Apertif, MeerKAT, and ASKAP (Braun et al. 2015; Morganti et al. 2015; Oosterloo et al. 2010; Gupta et al. 2016; Johnston et al. 2007) and the Deep Synoptic Array 2000 (DSA-2000, Hallinan et al. 2019, will produce blind large-scale surveys of H I absorption that will increase the number of radio AGN with H I absorption, particularly at low radio powers. In addition to large-scale H I surveys to be conducted over the next decade, follow-up VLBI observations with new telescopes such as the next-generation Very Large Array (Murphy et al. 2018) and the Square Kilometre Array (SKA, Paragi et al. 2015), will be crucial for spatially resolving H I absorption against radio jets in order to quantify the impact of jet-ISM feedback on host galaxy properties (Nyland et al. 2018). The results from future H I absorption studies of much larger samples will ultimately inform numerical simulations of jet-ISM interactions (e.g., Wagner et al. 2012; Mukherjee et al. 2016, 2018a,b; Bicknell et al. 2018), thereby providing new insights into the role of jet-driven feedback in the broader context of galaxy evolution.

Acknowledgments

We thank the referee, Bjorn Emonts, for the insightful comments which helped to improve the manuscript. RS gratefully acknowledge support from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC Advanced Grant RADIOLIFE-320745. Basic research in radio astronomy at the US Naval Research Laboratory is supported by 6.1 Base Funding. The European VLBI Network is a joint facility of independent European, African, Asian, and North American radio astronomy institutes. Scientific results from data presented in this publication are derived from the following EVN project code: GN002. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The Long Baseline Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The Arecibo Observatory is a facility of the National Science Foundation (NSF) operated by SRI International in alliance with the Universities Space Research Association (USRA) and UMET under a cooperative agreement. The Arecibo Observatory Planetary Radar Program is funded through the National Aeronautics and Space Administration (NASA) Near-Earth Objects Observations program. Based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). This research has made use of NASA’s Astrophysics Data System Bibliographic Services. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013). This research made use of APLpy, an open-source plotting package for Python (Robitaille & Bressert 2012).

References

  1. Abolfathi, B., Aguado, D. S., Aguilar, G., et al. 2018, ApJS, 235, 42 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aditya, J. N. H. S. 2019, MNRAS, 482, 5597 [Google Scholar]
  3. Aditya, J. N. H. S., & Kanekar, N. 2018, MNRAS, 473, 59 [NASA ADS] [CrossRef] [Google Scholar]
  4. Associated Universities, Inc. 1999, Astrophysics Source Code Library [record ascl:9911.003] [Google Scholar]
  5. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Baan, W. A., & Haschick, A. D. 1981, ApJ, 243, L143 [NASA ADS] [CrossRef] [Google Scholar]
  7. Baczko, A. K., Schulz, R., Kadler, M., et al. 2019, A&A, 623, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Best, P. N., & Heckman, T. M. 2012, MNRAS, 421, 1569 [NASA ADS] [CrossRef] [Google Scholar]
  9. Beswick, R. J., Pedlar, A., & Holloway, A. J. 2002, MNRAS, 329, 620 [NASA ADS] [CrossRef] [Google Scholar]
  10. Beswick, R. J., Peck, A. B., Taylor, G. B., & Giovannini, G. 2004, MNRAS, 352, 49 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bicknell, G. V., Mukherjee, D., Wagner, A. Y., Sutherland, R. S., & Nesvadba, N. P. H. 2018, MNRAS, 475, 3493 [NASA ADS] [CrossRef] [Google Scholar]
  12. Braun, R., Bourke, T., Green, J. A., Keane, E., & Wagg, J. 2015, Advancing Astrophysics with the Square Kilometre Array (AASKA14), 174 [Google Scholar]
  13. Buttiglione, S., Capetti, A., Celotti, A., et al. 2010, A&A, 509, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Callingham, J. R., Ekers, R. D., Gaensler, B. M., et al. 2017, ApJ, 836, 174 [NASA ADS] [CrossRef] [Google Scholar]
  15. Capetti, A., Celotti, A., Chiaberge, M., et al. 2002, A&A, 383, 104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Chandola, Y., Sirothia, S. K., & Saikia, D. J. 2011, MNRAS, 418, 1787 [NASA ADS] [CrossRef] [Google Scholar]
  17. Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&A, 562, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Cielo, S., Bieri, R., Volonteri, M., Wagner, A. Y., & Dubois, Y. 2018, MNRAS, 477, 1336 [NASA ADS] [CrossRef] [Google Scholar]
  19. Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 11 [NASA ADS] [CrossRef] [Google Scholar]
  20. de Gasperin, F., Merloni, A., Sell, P., et al. 2011, MNRAS, 415, 2910 [NASA ADS] [CrossRef] [Google Scholar]
  21. de Vries, N., Snellen, I. A. G., Schilizzi, R. T., Mack, K.-H., & Kaiser, C. R. 2009, A&A, 498, 641 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  23. Emonts, B. H. C., Morganti, R., Tadhunter, C. N., et al. 2005, MNRAS, 362, 931 [NASA ADS] [CrossRef] [Google Scholar]
  24. Evans, A. S., Sanders, D. B., Surace, J. A., & Mazzarella, J. M. 1999, ApJ, 511, 730 [NASA ADS] [CrossRef] [Google Scholar]
  25. Fiore, F., Feruglio, C., Shankar, F., et al. 2017, A&A, 601, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Floyd, D. J. E., Perlman, E., Leahy, J. P., et al. 2006, ApJ, 639, 23 [NASA ADS] [CrossRef] [Google Scholar]
  27. Geréb, K., Maccagni, F. M., Morganti, R., & Oosterloo, T. A. 2015, A&A, 575, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Gupta, N., Salter, C. J., Saikia, D. J., Ghosh, T., & Jeyakumar, S. 2006, MNRAS, 373, 972 [NASA ADS] [CrossRef] [Google Scholar]
  29. Gupta, N., Srianand, R., Baan, W., et al. 2016, Proc. MeerKAT Science: On the Pathway to the SKA, 25-27 May, 2016 Stellenbosch, South Africa (MeerKAT2016). Online at https://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=277, 14 [Google Scholar]
  30. Hallinan, G., Ravi, V., Weinreb, S., et al. 2019, Bull. Am. Astron. Soc., 51, 255 [Google Scholar]
  31. Harrison, C. M. 2017, Nat. Astron., 1, 0165 [NASA ADS] [CrossRef] [Google Scholar]
  32. Haschick, A. D., & Baan, W. A. 1985, ApJ, 289, 574 [NASA ADS] [CrossRef] [Google Scholar]
  33. Heckman, T. M. 2002, in Extragalactic Gas at Low Redshift, eds. J. S. Mulchaey, & J. T. Stocke, ASP Conf. Ser., 254, 292 [NASA ADS] [Google Scholar]
  34. Heckman, T. M., & Best, P. N. 2014, ARA&A, 52, 589 [Google Scholar]
  35. Heckman, T. M., Smith, E. P., Baum, S. A., et al. 1986, ApJ, 311, 526 [NASA ADS] [CrossRef] [Google Scholar]
  36. Helfand, D. J., White, R. L., & Becker, R. H. 2015, ApJ, 801, 26 [NASA ADS] [CrossRef] [Google Scholar]
  37. Holt, J. 2009, Astron. Nachr., 330, 226 [NASA ADS] [CrossRef] [Google Scholar]
  38. Holt, J., Tadhunter, C. N., & Morganti, R. 2008, MNRAS, 387, 639 [NASA ADS] [CrossRef] [Google Scholar]
  39. Holt, J., Tadhunter, C. N., Morganti, R., & Emonts, B. H. C. 2011, MNRAS, 410, 1527 [NASA ADS] [Google Scholar]
  40. Johnston, S., Bailes, M., Bartel, N., et al. 2007, PASA, 24, 174 [Google Scholar]
  41. Jones, D. L., Wehrle, A. E., Meier, D. L., & Piner, B. G. 2000, ApJ, 534, 165 [NASA ADS] [CrossRef] [Google Scholar]
  42. Kadler, M., Ros, E., Lobanov, A. P., Falcke, H., & Zensus, J. A. 2004, A&A, 426, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Kameno, S., Sawada-Satoh, S., Inoue, M., Shen, Z.-Q., & Wajima, K. 2001, PASJ, 53, 169 [NASA ADS] [CrossRef] [Google Scholar]
  44. King, A., & Pounds, K. 2015, ARA&A, 53, 115 [NASA ADS] [CrossRef] [Google Scholar]
  45. Labiano, A., García-Burillo, S., Combes, F., et al. 2013, A&A, 549, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Labiano, A., García-Burillo, S., Combes, F., et al. 2014, A&A, 564, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Maccagni, F. M., Morganti, R., Oosterloo, T. A., Geréb, K., & Maddox, N. 2017, A&A, 604, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Mahony, E. K., Morganti, R., Emonts, B. H. C., Oosterloo, T. A., & Tadhunter, C. 2013, MNRAS, 435, L58 [NASA ADS] [CrossRef] [Google Scholar]
  49. Mahony, E. K., Oonk, J. B. R., Morganti, R., et al. 2016, MNRAS, 455, 2453 [NASA ADS] [CrossRef] [Google Scholar]
  50. Martel, A. R., Baum, S. A., Sparks, W. B., et al. 1999, ApJS, 122, 81 [NASA ADS] [CrossRef] [Google Scholar]
  51. McNamara, B. R., & Nulsen, P. E. J. 2007, ARA&A, 45, 117 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  52. McNamara, B. R., & Nulsen, P. E. J. 2012, New J. Phys., 14, 055023 [NASA ADS] [CrossRef] [Google Scholar]
  53. Morganti, R. 2017, Front. Astron. Space Sci., 4, 42 [NASA ADS] [CrossRef] [Google Scholar]
  54. Morganti, R., Oosterloo, T., & Tsvetanov, Z. 1998, AJ, 115, 915 [NASA ADS] [CrossRef] [Google Scholar]
  55. Morganti, R., Oosterloo, T. A., Emonts, B. H. C., van der Hulst, J. M., & Tadhunter, C. N. 2003, ApJ, 593, L69 [NASA ADS] [CrossRef] [Google Scholar]
  56. Morganti, R., Tadhunter, C. N., & Oosterloo, T. A. 2005a, A&A, 444, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Morganti, R., Oosterloo, T. A., Tadhunter, C. N., van Moorsel, G., & Emonts, B. 2005b, A&A, 439, 521 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Morganti, R., Oosterloo, T., Struve, C., & Saripalli, L. 2008, A&A, 485, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Morganti, R., Fogasy, J., Paragi, Z., Oosterloo, T., & Orienti, M. 2013, Science, 341, 1082 [NASA ADS] [CrossRef] [Google Scholar]
  60. Morganti, R., Sadler, E. M., & Curran, S. 2015, Advancing Astrophysics with the Square Kilometre Array (AASKA14), 134 [CrossRef] [Google Scholar]
  61. Morganti, R., Veilleux, S., Oosterloo, T., Teng, S. H., & Rupke, D. 2016, A&A, 593, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Mukherjee, D., Bicknell, G. V., Sutherland, R., & Wagner, A. 2016, MNRAS, 461, 967 [NASA ADS] [CrossRef] [Google Scholar]
  63. Mukherjee, D., Bicknell, G. V., Sutherland, R., & Wagner, A. 2017, MNRAS, 471, 2790 [NASA ADS] [CrossRef] [Google Scholar]
  64. Mukherjee, D., Wagner, A. Y., Bicknell, G. V., et al. 2018a, MNRAS, 476, 80 [NASA ADS] [CrossRef] [Google Scholar]
  65. Mukherjee, D., Bicknell, G. V., Wagner, A. E. Y., Sutherland, R. S., & Silk, J. 2018b, MNRAS, 479, 5544 [NASA ADS] [CrossRef] [Google Scholar]
  66. Murphy, E. J., Bolatto, A., & Chatterjee, S. 2018, in The ngVLA Science Case and Associated Science Requirements, ed. E. Murphy, ASP Conf. Ser., 517, 3 [NASA ADS] [Google Scholar]
  67. Murthy, S., Morganti, R., Oosterloo, T., et al. 2019, A&A, 629, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Newville, M., Stensitzki, T., Allen, D. B., & Ingargiola, A. 2014, https://doi.org/10.5281/zenodo.11813 [Google Scholar]
  69. Nyland, K., Harwood, J. J., Mukherjee, D., et al. 2018, ApJ, 859, 23 [NASA ADS] [CrossRef] [Google Scholar]
  70. O’Dea, C. P. 1998, PASP, 110, 493 [NASA ADS] [CrossRef] [Google Scholar]
  71. O’Dea, C. P., Koekemoer, A. M., Baum, S. A., et al. 2001, AJ, 121, 1915 [NASA ADS] [CrossRef] [Google Scholar]
  72. Oosterloo, T. A., Morganti, R., Tzioumis, A., et al. 2000, AJ, 119, 2085 [NASA ADS] [CrossRef] [Google Scholar]
  73. Oosterloo, T., Verheijen, M., & van Cappellen, W. 2010, Proceedings of the ISKAF2010 Science Meeting, 43 [CrossRef] [Google Scholar]
  74. Oosterloo, T., Raymond Oonk, J. B., Morganti, R., et al. 2017, A&A, 608, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Orienti, M. 2016, Astron. Nachr., 337, 9 [NASA ADS] [CrossRef] [Google Scholar]
  76. Paragi, Z., Godfrey, L., Reynolds, C., et al. 2015, Advancing Astrophysics with the Square Kilometre Array (AASKA14), 143 [Google Scholar]
  77. Robitaille, T., & Bressert, E. 2012, Astrophysics Source Code Library [record ascl:1208.017] [Google Scholar]
  78. Sawada-Satoh, S., Kameno, S., Nakamura, K., et al. 2008, ApJ, 680, 191 [NASA ADS] [CrossRef] [Google Scholar]
  79. Schilizzi, R. T., Tian, W. W., Conway, J. E., et al. 2001, A&A, 368, 398 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Schulz, R., Morganti, R., Nyland, K., et al. 2018, A&A, 617, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Shepherd, M. 2011, Astrophysics Source Code Library [record ascl:1103.001] [Google Scholar]
  82. Shepherd, M. C., Pearson, T. J., & Taylor, G. B. 1994, BAAS, 26, 987 [NASA ADS] [Google Scholar]
  83. Silk, J., & Rees, M. J. 1998, A&A, 331, L1 [NASA ADS] [Google Scholar]
  84. Snellen, I. A. G., Mack, K.-H., Schilizzi, R. T., & Tschager, W. 2004, MNRAS, 348, 227 [NASA ADS] [CrossRef] [Google Scholar]
  85. Struve, C., & Conway, J. E. 2010, A&A, 513, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Tadhunter, C. 2008, New Astron. Rev., 52, 227 [Google Scholar]
  87. Tadhunter, C., Robinson, T. G., González Delgado, R. M., Wills, K., & Morganti, R. 2005, MNRAS, 356, 480 [NASA ADS] [CrossRef] [Google Scholar]
  88. Tremblay, G. R., O’Dea, C. P., Baum, S. A., et al. 2010, ApJ, 715, 172 [NASA ADS] [CrossRef] [Google Scholar]
  89. van Gorkom, J. H., van der Hulst, J. M., Haschick, A. D., & Tubbs, A. D. 1990, AJ, 99, 1781 [NASA ADS] [CrossRef] [Google Scholar]
  90. Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 769 [NASA ADS] [CrossRef] [Google Scholar]
  91. Veilleux, S., Maiolino, R., Bolatto, A. D., & Aalto, S. 2020, A&ARv, 28, 2 [NASA ADS] [CrossRef] [Google Scholar]
  92. Wagner, A. Y., & Bicknell, G. V. 2011, ApJ, 728, 29 [NASA ADS] [CrossRef] [Google Scholar]
  93. Wagner, A. Y., Bicknell, G. V., & Umemura, M. 2012, ApJ, 757, 136 [NASA ADS] [CrossRef] [Google Scholar]
  94. Westmeier, T., Jurek, R., Obreschkow, D., Koribalski, B. S., & Staveley-Smith, L. 2014, MNRAS, 438, 1176 [NASA ADS] [CrossRef] [Google Scholar]
  95. Zovaro, H. R. M., Sharp, R., Nesvadba, N. P. H., et al. 2019, MNRAS, 484, 3393 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Fits to H I absorption sepctra

thumbnail Fig. A.1.

Fits of Gaussian functions to inverted H I absorption spectra of 4C 52.37. The fit paramaters are given in Table A.1. Top: WSRT spectum from Fig. 2 with two Gaussian functions fitted to it (red and yellow dotted lines). The combination of the two fit functions is shown as the black dashed line. The data is colored in gray. Bottom: clipped integrated VLBI spectrum from Fig. 2. Three Gaussian components were fitted to it (colored dotted lines) and their combination is shown as the black dashed line.

thumbnail Fig. A.2.

Inverted VLA H I absorption spectrum of 3C 293 of the western component from Fig. 4. Data from Mahony et al. (2013). Three Gaussian components were fitted to it (colored dotted lines) and their combination is shown as the black dashed line. Fit parameters are listed in Table A.1.

thumbnail Fig. A.3.

Fits of Gaussian functions to inverted H I absorption spectra of 3C 236. Data are taken from Schulz et al. (2018). Fit parameters are listed in Table A.1. Top: VLA spectrum with three Gaussian components fitted to it (coloured dotted lines). Their combination is shown as the black dashed line. Bottom: clipped integrated VLBI sepctra for the core and lobe region of 3C 236.

thumbnail Fig. A.4.

Inverted WSRT H I absorption spectrum of 4C 12.50 from Morganti et al. (2013). Three Gaussian components were fitted to it (colored dotted lines) and their combination is shown as the black dashed line. Fit parameters are listed in Table A.1.

Table A.1.

Parameters of the Gaussian fit to the H I absorption spectra.

All Tables

Table 1.

Properties of observations.

Table 2.

Properties of images.

Table 3.

Properties of the kinematically disturbed H I gas.

Table A.1.

Parameters of the Gaussian fit to the H I absorption spectra.

All Figures

thumbnail Fig. 1.

Continuum images of 4C 52.37 (top) and 3C 293 (bottom) convolved with a synthesized beam of 20 mas and 25 mas, respectively, which is marked by the gray circle. The top panel also shows a higher resolution (6 mas) zoom-in of the central region of the radio source. The contour lines begin at three times the noise level of the corresponding image and increase logarithmically by a factor of 2. The solid, horizontal scale bar represents the spatial scales.

In the text
thumbnail Fig. 2.

H I absorption spectra of 4C 52.37. The unresolved WSRT spectrum is shown in red. The two spatially integrated VLBI profiles (in blue) were extracted from a region of 60 mas × 40 mas centered on E1. The dashed blue line corresponds to the VLBI spectrum that was obtained by integrating only pixels with a value of ≤ − 3σVLBI, cube, while for the solid blue line the integration was only limited to the region marked by the 3σVLBI, cont contour.

In the text
thumbnail Fig. 3.

Top panel: position-velocity (PV) diagram of 4C 52.37 along a cut aligned with the largest extent of the central absorption feature. The solid, horizontal line corresponds to the systemic velocity. The dashed, vertical lines mark the size of the synthesized beam. Contour lines correspond to 3 (grey, solid), −2 (black, dashed), −3, and −5 (black, solid) times the noise level in the cube. Middle panel: contour lines show the continuum radio emission as in Fig. 1 with the color map corresponding to N H I T spin 1 $ N_{{\mathrm{H}\,\textsc{I}}} T^{-1}_{\mathrm{spin}} $. The dashed lines represent the cuts for the PV diagram. Bottom panel: same as top panel, but a cut along the position angle of W1 and E1. The top and bottom panels are centred on E1.

In the text
thumbnail Fig. 4.

H I absorption spectra of 3C 293. The low-resolution VLA spectrum is shown in red. The solid blue and dashed violet lines correspond to the spatially integrated VLBI spectra of the eastern and western part of the source. The dotted, gray line represents the total integrated VLBI spectrum. The VLBI spectra were obtained by only integrating pixels with a value of ≤ − 3σVLBI, cube.

In the text
thumbnail Fig. 5.

Left panel: position-velocity diagram of a representative location in the eastern part and column density normalized by spin temperature (bottom panel) for the same velocity range as in the right panel. Right panel: position-velocity diagram of 3C 293 at the location of the VLBI core (top panel) and column density normalized by spin temperature separated into systemic and outflowing H I gas (bottom two panels). The dashed lines in the top two panels correspond to the size of the synthesized beam. The dashed lines in the bottom panels correspond to the location of the slice for the position-velocity diagram.

In the text
thumbnail Fig. A.1.

Fits of Gaussian functions to inverted H I absorption spectra of 4C 52.37. The fit paramaters are given in Table A.1. Top: WSRT spectum from Fig. 2 with two Gaussian functions fitted to it (red and yellow dotted lines). The combination of the two fit functions is shown as the black dashed line. The data is colored in gray. Bottom: clipped integrated VLBI spectrum from Fig. 2. Three Gaussian components were fitted to it (colored dotted lines) and their combination is shown as the black dashed line.

In the text
thumbnail Fig. A.2.

Inverted VLA H I absorption spectrum of 3C 293 of the western component from Fig. 4. Data from Mahony et al. (2013). Three Gaussian components were fitted to it (colored dotted lines) and their combination is shown as the black dashed line. Fit parameters are listed in Table A.1.

In the text
thumbnail Fig. A.3.

Fits of Gaussian functions to inverted H I absorption spectra of 3C 236. Data are taken from Schulz et al. (2018). Fit parameters are listed in Table A.1. Top: VLA spectrum with three Gaussian components fitted to it (coloured dotted lines). Their combination is shown as the black dashed line. Bottom: clipped integrated VLBI sepctra for the core and lobe region of 3C 236.

In the text
thumbnail Fig. A.4.

Inverted WSRT H I absorption spectrum of 4C 12.50 from Morganti et al. (2013). Three Gaussian components were fitted to it (colored dotted lines) and their combination is shown as the black dashed line. Fit parameters are listed in Table A.1.

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.