## Feedback from low-luminosity radio galaxies: B2 0258+35

Suma Murthy 1 2 , Ra ; GLYPH&lt;11&gt; aella Morganti 2 1 ; , Tom Oosterloo 2 1 ; , Robert Schulz 2 , Dipanjan Mukherjee 3 , Alexander Y. Wagner 4 , Geo GLYPH&lt;11&gt; rey Bicknell 5 , Isabella Prandoni 6 , Aleksandar Shulevski 7

- 1 Kapteyn Astronomical Insititute, University of Groningen, Postbus 800, 9700 AV Groningen, The Netherlands e-mail: murthy@astro.rug.nl
- 2 ASTRON, the Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands
- 3 Dipartimento di Fisica Generale, Universita degli Studi di Torino , Via Pietro Giuria 1, 10125 Torino, Italy
- 4 Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577
- 5 Australian National University, Research School of Astronomy and Astrophysics, Cotter Rd., Weston, ACT 2611, Australia
- 6 INAF - Istituto di Radioastronomia, via P. Gobetti 101, 40129 Bologna, Italy
- 7 Anton Pannekoek Institute for Astronomy, University of Amsterdam, Postbus 94249, 1090 GE Amsterdam, The Netherlands

Received ..., accepted ...

## ABSTRACT

Low-luminosity radio-loud active galactic nuclei (AGN) are of importance in studies concerning feedback from radio AGN since a dominant fraction of AGN belong to this class. We report high-resolution Very Large Array (VLA) and European VLBI Network (EVN) observations of H i 21cm absorption from a young, compact steep-spectrum radio source, B2 0258 + 35, nested in the early-type galaxy NGC 1167, which contains a 160 kpc H i disc. Our VLA and EVN H i absorption observations, modelling, and comparison with molecular gas data suggest that the cold gas in the centre of NGC 1167 is very turbulent (with a velocity dispersion of GLYPH&lt;24&gt; 90 km s GLYPH&lt;0&gt; 1 ) and that this turbulence is induced by the interaction of the jets with the interstellar medium (ISM). Furthermore, the ionised gas in the galaxy shows evidence of shock heating at a few kpc from the radio source. These findings support the results from numerical simulations of radio jets expanding into a clumpy gas disc, which predict that the radio jets in this case percolate through the gas disc and drive shocks into the ISM at distances much larger than their physical extent. These results expand the number of low-luminosity radio sources found to impact the surrounding medium, thereby highlighting the possible relevance of these AGN for feedback.

Key words. galaxies: active - galaxies: individual: B2 0258 + 35 - radio lines: galaxies - galaxies: ISM

## 1. Introduction

Active galactic nuclei (AGN) have long been of importance in studies of galaxy evolution. The energy emitted by an active supermassive black hole (SMBH) is believed to a GLYPH&lt;11&gt; ect its host galaxy in a number of di GLYPH&lt;11&gt; erent ways, collectively referred to as AGN feedback. These e GLYPH&lt;11&gt; ects can range from regulating the growth of the SMBH (which can explain the relation between the SMBH mass and the velocity dispersion of the galactic bulge; e.g. Silk &amp; Rees 1998), and the possible quenching of the star formation via massive gas outflows (Harrison et al. 2018, and references therein), to preventing the cooling of the gas in hot halos, which in turn can prevent massive galaxies from forming (e.g. Croton et al. 2006; McNamara &amp; Nulsen 2012, and references therein).

(L1.4 GHz . 10 24 5 : WHz 1 ) are of interest to studies concerning GLYPH&lt;0&gt; this mechanical mode of feedback since they greatly outnumber their high-luminosity counterparts (Best et al. 2005).

It is thought that radio-loud AGN play an important role mostly in the last: by preventing the cooling of the circumgalactic and intergalactic gas through their large-scale radio lobes, they prevent the accretion of gas onto the galaxy. This e GLYPH&lt;11&gt; ect has been clearly observed in a number of relatively bright radio AGN that have been studied in X-rays (e.g. McNamara et al. 2009; Randall et al. 2011; Cavagnolo et al. 2011; Fabian 2012).

However, over the past years evidence has been building up to show that radio AGN also a GLYPH&lt;11&gt; ect the host galaxy directly through radio jet-ISM interactions (e.g. Oosterloo et al. 2017; Rodríguez-Ardila et al. 2017; Fabbiano et al. 2018; Maksym et al. 2018; Wang et al. 2011; Finlez et al. 2018; Alatalo et al. 2015; Croston et al. 2008). Low-luminosity radio AGN

Recent numerical simulations (e.g. Cielo et al. 2018; Mukherjee et al. 2018a) also support the importance of this mode of feedback from radio AGN. They predict that jets expanding into a clumpy gas disc a GLYPH&lt;11&gt; ect the gas kinematics significantly. The extent of this impact depends on the orientation of the jets with respect to the gas disc, among other parameters. If the jets expand directly into the gas disc, then the simulations predict that they will induce fast outflows of multi-phase gas as well as strong turbulence. These simulations further predict that the jetISM interaction will generate expanding cocoons that shock the gas in the disc and also a GLYPH&lt;11&gt; ect the star formation rate in the regions surrounding the AGN activity. The gas may also, in turn, a GLYPH&lt;11&gt; ect the evolution of the radio source itself (e.g. morphology and age; Mukherjee et al. 2018b).

This mechanical feedback is best studied by probing the kinematics of di GLYPH&lt;11&gt; erent phases of gas a GLYPH&lt;11&gt; ected by the radio jets. Such detailed studies of radio jet-ISM interactions can provide constraints to quantify the AGN feedback e GLYPH&lt;11&gt; ects, which can be incorporated into the simulations of galaxy evolution. In order to gain a better understanding of how this actually happens, we need to expand the number of cases for which detailed, highresolution observations of multiple phases of the gas are available so that a comparison with the predictions of numerical simulations is possible.

We present here one such study of a nearby ( z = 0.0165; Wegner et al. 1993), low-luminosity (L1 4GHz : = 2.1 GLYPH&lt;2&gt; 10 23 W Hz GLYPH&lt;0&gt; 1 ; Shulevski et al. 2012) radio source: B2 0258 + 35. This source is nested in a gas-rich galaxy, NGC 1167, and appears to be interacting strongly with the ISM. The galactic H i is well studied (Struve et al. 2010) and so are the properties of the radio source (Sanghera et al. 1995; Giroletti et al. 2005; Brienza et al. 2018). This radio galaxy is thus a good candidate to understand the nature of jet-ISM interactions and the possible impact the radio source has on the ambient medium.

B2 0258 + 35 is classified as a compact steep-spectrum (CSS) source (Sanghera et al. 1995). The sources in this class are compact with their synchrotron spectrum peaking at a few hundred MHz and are generally considered to be young (e.g. O'Dea 1998). Earlier radio continuum studies of B2 0258 + 35 (Giroletti et al. 2005) have found the radio source to be about 1 kpc in size, consisting of a core and two jets (see Figure 1, inset). The northern jet is comparatively faint and the southern jet is bent sharply. This asymmetry both in brightness and morphology suggests an interaction with the ISM, similar to that seen in IC 5063 (Morganti et al. 1998). The age of B2 0258 + 35 is estimated to be between 0.4 Myr (Brienza et al. 2018) and 0.9 Myr (Giroletti et al. 2005). Following the correlation between the age and the size of the source (Orienti &amp; Dallacasa 2014; Murgia 2003), it has also been found that the source is smaller than expected for this age, given an unconstrained expansion. This suggests that it may have been confined within the host galaxy for a protracted length of time (Brienza et al. 2018). However, we note that a hypothesis based on the above-mentioned correlations alone is quite uncertain, but as we show in the subsequent sections, the H i absorption and CO emission studies further strengthen this hypothesis.

The host galaxy, NGC 1167, was further found to host largescale ( GLYPH&lt;24&gt; 240 kpc) extremely low surface brightness radio lobes (Shulevski et al. 2012), about GLYPH&lt;24&gt; 110 Myr old (Brienza et al. 2018), which make up only 3% of the total source luminosity (4.75 mJy arcmin GLYPH&lt;0&gt; 2 at 145 MHz). Though such low surface brightness may imply that these structures are remnants of a previous activity, their spectral indices have been found to be 'normal' and not ultra-steep as expected for these sources (Brienza et al. 2018). Various possibilities have been put forth to explain this feature, including in situ particle reacceleration, multiple duty cycles, dense ISM smothering the large-scale jets that are still being fuelled at a low level, and magnetic draping enhancing the mixing within the lobes while suppressing the same between the radio lobes and the surrounding medium (Brienza et al. 2018; Adebahr et al. 2019).

NGC 1167 is a gas-rich (MH i = 1.5 GLYPH&lt;2&gt; 10 10 M ) early-type GLYPH&lt;12&gt; galaxy, with a regularly rotating 160 kpc H i disc (seen in emission; Figure 1). The disc has very regular kinematics within a radius of 65 kpc, and only in the very outer parts does it show signs of interactions, likely with a satellite galaxy. This indicates that the galaxy has not undergone a major merger in the last few billion years.

H i has also been detected in absorption against the CSS source with the WSRT (Struve et al. 2010) at a low spatial resolution ( GLYPH&lt;24&gt; 10 kpc; Figure 1). CO (1-0) emission has also been detected from NGC 1167 (Prandoni et al. 2007; O'Sullivan et al. 2015; Bolatto et al. 2017). The galaxy had originally been optically classified as a Seyfert 2 galaxy (Ho et al. 1997). However, more recent observations show that the central optical AGN has a LINER (low ionization nuclear emission-line region) spectrum (see Emonts 2006). NGC 1167 is also a part of the CALIFA

Fig. 1. Large-scale H i disc seen in emission along with the velocity field (in colour). The eastern portion of the rotating disc is the approaching side. H i is seen in absorption within a radius of GLYPH&lt;24&gt; 10 kpc (angular resolution of the WSRT observations GLYPH&lt;24&gt; 29 00 ; Struve et al. 2010). The largescale low surface brightness radio lobes (WSRT observations with a beam of 39 00 GLYPH&lt;2&gt; 33 00 ; Shulevski et al. 2012) are shown as black contours. The contours start at 400 GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 and increase by a factor of 2. The central CSS source (at 22 GHz with VLA A-array at an angular resolution of 0.12 00 Giroletti et al. 2005) is shown in the inset. The position angle of the H i disc is 75 GLYPH&lt;14&gt; (Struve et al. 2010), while that of the CSS source is 132 GLYPH&lt;14&gt; (Giroletti et al. 2005).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-4c6e0af7-f3d0-49f2-8f28-983300699bcf.jpg)

survey and their integral field unit (IFU) studies further confirm the LINER nature of the nucleus (Gomes et al. 2016a).

Here, we present high-resolution Very Large Array (VLA) and European VLBI Network (EVN) H i absorption studies of the central CSS source B2 0258 + 35. We detect resolved, unusually broad H i absorption that closely matches the CO emission profile, pointing towards a circumnuclear structure that is being disturbed as a result of strong interaction with the expanding radio jets.

In Section 2 we describe the VLA and EVN observations. We present the results in Section 3, discuss the implications of our results for studies seeking to quantify AGN feedback in Section 4, and provide a summary in Section 5.

Wherever required we have assumed a flat universe with H0 = 67.3 km s GLYPH&lt;0&gt; 1 Mpc GLYPH&lt;0&gt; 1 , GLYPH&lt;10&gt; GLYPH&lt;3&gt; = 0.685, and GLYPH&lt;10&gt; M = 0.315 (Planck Collaboration et al. 2014). At z = 0.0165, 1 00 corresponds to 0.349 kpc.

## 2. Observations and data reduction

## 2.1. Very Large Array observations

The H i observations with the VLA A-array were carried out in October 2008 (Proposal id: AS0955) for a total on-source time of 5.82 hours. A bandwidth of 6.25 MHz, centred at the redshifted H i line frequency of 1396.5 MHz, subdivided into 64 channels was used as the set-up. This gave us a raw velocity resolution of 21 km s GLYPH&lt;0&gt; 1 . 3C147 was used for flux calibration

Murthy et al.: Low-luminosity AGN: B2 0258 + 35

| Telescope    | GLYPH<23> obs (GHz)   | on-source (hours)   | BW (MHz)   | vel. resln (km s GLYPH<0> 1 )   | beam-size 00 GLYPH<2> 00 (6)   | PA ( GLYPH<14> ) (7)   | RMS map ( GLYPH<22> Jy beam GLYPH<0> 1 ) (8)   | RMS cube ( GLYPH<22> Jy beam GLYPH<0> 1 channel GLYPH<0> 1 )   |
|--------------|-----------------------|---------------------|------------|---------------------------------|--------------------------------|------------------------|------------------------------------------------|----------------------------------------------------------------|
| (1)          | (2)                   | (3)                 | (4)        | (5)                             | GLYPH<2>                       |                        |                                                | (9)                                                            |
| VLA          | 1.3965                | 5.85                | 6.25       | 21                              | 1.13 0.98                      | -75.11                 | 340                                            | 840                                                            |
| EVN GLYPH<3> | 1.4054                | 6.6                 | 16         | 6.7                             | 0.04923 GLYPH<2> 0.04136       | 61.32                  | 100                                            | 410                                                            |

Table 1. Observation details: The columns are: (1) Telescope used; (2) Central frequency of the band (GHz); (3) On-source time in hours; (4) Bandwidth (MHz); (5) Velocity resolution (km s); (6) Beam size; (7) Beam position angle ( ); (8) RMS noise on the continuum map ( / GLYPH&lt;14&gt; GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 ); and (9) RMS noise on the spectral cube at a velocity resolution mentioned in (5) ( GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 channel GLYPH&lt;0&gt; 1 ).

GLYPH&lt;3&gt; Participating antennas: E GLYPH&lt;11&gt; elsberg (Germany), Westerbork (the Netherlands), Jodrell-Bank (the United Kingdom), Medicina (Italy), Noto (Italy), Onsala (Sweden), Torun (Poland).

and 0319 + 415 for phase and bandpass calibration. The details of the observations are listed in Table 1. The observations were carried out when the VLA was being upgraded to the EVLA. The digital signal from the EVLA antennas had to be converted to analogue signals before being fed to the correlator, and this caused aliasing of power in the bottom 2 MHz of the band in all the EVLA-EVLA baselines. These baselines were flagged to avoid complications in solving for bandpasses and gains. Data were reduced in 'classic' AIPS (Astronomical Image Processing Software). After the initial flagging of EVLA-EVLA baselines and bad data, antenna-dependent gains, and bandpass solutions were determined using the observations on the calibrators. The gain solutions were further improved iteratively through the selfcalibration procedure (imaging and phase-only self-calibration cycles until the continuum map showed no further improvement) after which a round of amplitude and phase self-calibration and imaging was done. Then the continuum model was subtracted from the calibrated visibilities and the residual UV data affected by radio frequency interference (RFI) were flagged. The final continuum model was subtracted from the calibrated multichannel UV dataset. Residual continuum emission was removed by fitting a linear polynomial to each visibility spectrum. Finally, the continuum-subtracted data were shifted to the heliocentric frame. This dataset was then imaged to obtain the spectral cube.

The continuum map has an RMS noise of GLYPH&lt;24&gt; 340 GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 and a restoring beam of 1.11 00 GLYPH&lt;2&gt; 1.02 00 with a position angle of 75.11 . This was made using ROBUST GLYPH&lt;14&gt; = -2 weighting, averaging all the line-free channels together. The peak flux density of the target is 1.08 GLYPH&lt;6&gt; 0.05 Jy beam GLYPH&lt;0&gt; 1 . The integrated flux density is 1.97 GLYPH&lt;6&gt; 0.10 Jy. The uncertainty on the flux density scale at the observed frequency is assumed to be 5% (Perley &amp; Butler 2017). The continuum-subtracted spectral cube has the same restoring beam as the continuum map. The noise on the cube is GLYPH&lt;24&gt; 840 GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 channel GLYPH&lt;0&gt; 1 for a channel width of 21 km s GLYPH&lt;0&gt; 1 , without any spectral smoothing.

## 2.2. EVN observations

The phase-referenced H i observations with the European VLBI Network (EVN) were carried out in October 2012 (Project ID: ES070) for a total on-source time of 6.6 hours. We used a bandwidth of 16 MHz subdivided into 512 channels. Data reduction was carried out using 'classic' AIPS. The details of the observations are listed in Table 1.

We corrected for the instrumental delay, and then for delay and rate as a function of time after initial flagging of bad data, carried out bandpass calibration, and applied the solutions to the target. We then self-calibrated the target dataset to further improve the gain solutions. The self-calibration iterations are similar to those performed for the VLA data, but with the solution interval being shortened progressively from 30 minutes to 2 minutes for phase-only self-calibration. We applied these gain solutions to the UV data, subtracted the final continuum model from the calibrated UV data, and then flagged the residuals for RFI. We subtracted residual continuum emission by fitting a second-order polynomial to each visibility spectrum and imaged the residual UV data to obtain a spectral cube after shifting to the heliocentric frame.

The final continuum map (obtained with natural weighting ROBUST = 5 and averaging all the channels together) has an RMS noise of 100 GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 and a beam of 49.23 mas x 41.36 mas with a position angle of 61.32 . The peak GLYPH&lt;14&gt; flux of the target is 95 mJy beam GLYPH&lt;0&gt; 1 and the integrated flux is 193 mJy. This was also estimated from the AIPS task JMFIT with a single component Gaussian fit.

The spectral cube was also obtained with natural weighting (ROBUST = 5) and with the same restoring beam as the continuum map. The cube has an RMS noise of 0.4 mJy beam GLYPH&lt;0&gt; 1 channel GLYPH&lt;0&gt; 1 for a channel resolution of 6.7 km s GLYPH&lt;0&gt; 1 . We also tried other weighting schemes (uniform and robust) and found that the H absorption (see Section 3) was detectable only in the spectral i cube obtained with natural weighting.

## 3. Results

## 3.1. Continuum structure

Figure 2 shows the VLA continuum map of B2 0258 + 35 with the higher resolution 22 GHz map of Giroletti et al. (2005) overlaid. In the VLA map the source is slightly resolved and is about 1 kpc in size. Higher resolution observations (Sanghera et al. 1995; Giroletti et al. 2005) have shown that B2 0258 + 35 has plumelike lobes. In addition, there is a bright knot GLYPH&lt;24&gt; 0.1 00 south-east of the core. Further to the south-east the southern jet undergoes a sharp bend. There is no bright spot or bending in the northern jet. This asymmetry between the northern and the southern jets is reflected in our VLA map at a spatial resolution of GLYPH&lt;24&gt; 0.35 kpc. We note that the integrated flux density from our observations is consistent with that from the WSRT observations at a resolution of GLYPH&lt;24&gt; 10.1 kpc (1.8 Jy; Struve et al. 2010) within the flux measurement errors. Thus, there is no significant radio continuum outside the central GLYPH&lt;24&gt; 1 kpc region.

With the EVN we recover only three components (see Figure 3): the faint unresolved core (with an integrated flux density of GLYPH&lt;24&gt; 2.5 mJy), the bright extended knot ( GLYPH&lt;24&gt; 185 mJy) coincident with the peak of the VLA map and about 75 mas in size, and another knot GLYPH&lt;24&gt; 0.2 00 south-east of the bright knot ( GLYPH&lt;24&gt; 5 mJy), also resolved with an angular size of GLYPH&lt;24&gt; 55 mas. We do not detect any jet-like structures connecting the core and this bright spot. The total integrated flux density is GLYPH&lt;24&gt; 193 mJy, corresponding to about 10% of the flux density recovered with the VLA observations. Thus, the rest of the continuum flux is di GLYPH&lt;11&gt; use and distributed over length scales larger than the largest recoverable

Fig. 2. Top left: VLA continuum map of B2 0258 + 35 ( beam size: 1.13 00 GLYPH&lt;2&gt; 0.98 00 ; red contours). The contour levels start from 1 mJy beam GLYPH&lt;0&gt; 1 (3 GLYPH&lt;27&gt; ) and increase by a factor of two. The 3 GLYPH&lt;27&gt; negative contours are shown in grey and the beam is shown in the bottom left corner of the image. The regions from which the H i absorption spectra are extracted are labelled. The 22 GHz VLA map (by Giroletti et al. 2005, beam size: 0.12 00 GLYPH&lt;2&gt; 0.12 00 ) is overlaid in blue contours to show the underlying continuum structure. These contours start from 325 GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 (5 GLYPH&lt;27&gt; ) and increase by a factor of two. Bottom left, bottom right, and top right : The absorption spectra extracted from the regions are in black in the continuum map. The spectra have been Hanning-smoothed and resampled to a velocity resolution of GLYPH&lt;24&gt; 42 km s GLYPH&lt;0&gt; 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-62c95d19-5974-4581-b586-f7a51f7435c8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-e0e8573e-c2a6-45a8-8533-cee9001a859c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-21114199-0b0c-4c6f-b753-4354b273aa5e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-e07e9508-3dc1-417d-8f19-fb2da759f61e.jpg)

angular scale of our EVN array, GLYPH&lt;24&gt; 200 mas (which corresponds to 80 pc). B2 0258 + 35 has also been observed with the VLBA at 1.6 GHz by Giroletti et al. (2005). They report the detection of the bright knot and the core, but not the third component to the south of the bright knot. They measure a total flux density of 243 mJy over the entire radio source, consistent with our measurement within the flux measurement errors. However, the peak flux density we measure is an order of magnitude higher than their measurement. These di GLYPH&lt;11&gt; erences could be due to the much higher spatial resolution (a beam size of 10.4 mas GLYPH&lt;2&gt; 4.3 mas) of their observations.

## 3.2. H I absorption

We detect slightly resolved H i absorption against the radio continuum in our VLA data. The absorption profiles against di GLYPH&lt;11&gt; er- ent parts of the continuum are shown in Figure 2. Since the continuum is only a few beams across, we extracted the spectra from regions su GLYPH&lt;14&gt; ciently far apart so that each spectrum extracted is from a di GLYPH&lt;11&gt; erent beam element: against the northern jet, the region including the core and the bright knot, and the southern jet. We see strong absorption against the continuum in the beam elements covering the bright knot, and the southern jet. We also detect weak absorption against the northern jet. The full width at zero intensity (FWZI) of these H i absorption profiles is GLYPH&lt;24&gt; 545 km s GLYPH&lt;0&gt; 1 against the southern jet and the hotspot, while it is GLYPH&lt;24&gt; 500 km s GLYPH&lt;0&gt; 1 against the northern jet. The column density of the gas detected at the locations shown in Figure 2 are consistent with each other within the measurement errors; they range from (3 GLYPH&lt;6&gt; 0.6) GLYPH&lt;2&gt; 10 20 cm GLYPH&lt;0&gt; 2 against the northern jet to (2.2 GLYPH&lt;6&gt; 0.1) GLYPH&lt;2&gt; 10 20 cm GLYPH&lt;0&gt; 2 against the hotspot. This was estimated by assuming a spin temperature of 100 K and a covering factor of unity.

Fig. 3. The VLBI map of B2 0258 + 35 (beam size: 0.04923 00 GLYPH&lt;2&gt; 0.04136 00 ; red contours). The contours start at 400 GLYPH&lt;22&gt; Jy beam GLYPH&lt;0&gt; 1 (4 GLYPH&lt;27&gt; ) and increase by a factor of two at each step. The 4 GLYPH&lt;27&gt; negative contours are shown in grey and the beam is shown in the bottom right corner of the map. The black cross marks the VLBI core. The 22 GHz VLA image from Giroletti et al. (2005) is overlaid in blue contours for reference, as in Figure 2. Only 10% of the total VLA flux is recovered.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-28999833-1fe0-4cca-8324-dea5981e899e.jpg)

In addition, Figure 4 shows a comparison between the integrated VLA and the WSRT H i absorption profiles (Struve et al. 2010). Given the smaller bandwidth of the VLA observations and the higher noise level, the two profiles are in good agreement with each other. Thus, the higher spatial resolution of the VLA observations does not result in non-detection of a significant fraction of H i .

With the EVN observations, we do not detect any absorption at the location of the peak of the knot or the southern blob or against the core. However, we do find an absorption feature when integrated over the region of the knot (see Figure 5). The absorption profile has an FWZI of GLYPH&lt;24&gt; 322 km s GLYPH&lt;0&gt; 1 and a peak absorption of -0.85 mJy. Since a considerable amount of the continuum flux density is resolved out at the VLBI scale, the peak of the recovered absorption is GLYPH&lt;24&gt; 5 times lower than that of the VLA absorption spectrum. However, the H i column density derived is (1.70 GLYPH&lt;6&gt; 0.21) GLYPH&lt;2&gt; 10 20 cm GLYPH&lt;0&gt; 2 , which is consistent with the column density measured against the knot from the VLA spectrum. This detection thus implies a uniform coverage of the radio continuum by the absorbing gas and hence a covering factor of unity.

## 4. Discussion

With the information collected, we first wanted to investigate whether the H i seen in absorption is a part of the large-scale disc and is distributed over a few tens of kpc or is co-spatial with the jet and if so, whether there is any signature of jet-ISM interaction. We did this by comparing the absorption profile with a model absorption profile arising from the large-scale disc and then combining this with the results from molecular gas emission.

Fig. 4. Integrated VLA (blue) and WSRT (red) H i absorption profiles. The profiles are in excellent agreement with each other, implying there is no di GLYPH&lt;11&gt; use component lost in continuum and in absorption.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-83af349b-6987-408c-a1b3-a960f745145a.jpg)

Fig. 5. Integrated H i absorption spectrum against the knot (Figure 3) recovered with the EVN. The spectrum has been Hanning-smoothed to a velocity resolution of 54 km s GLYPH&lt;0&gt; 1 . The RMS noise on the Hanningsmoothed spectrum is 30 GLYPH&lt;22&gt; Jy and the FWZI of the absorption feature is GLYPH&lt;24&gt; 322 km s GLYPH&lt;0&gt; 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-8841ea98-9175-47ec-9e46-2bfa818691d0.jpg)

## 4.1. Origin of H I absorption

The shape of the absorption profile, in general, depends on the distribution of gas and on the morphology of the background radio source. As Figure 1 shows, the large-scale H i disc has a position angle of 75 GLYPH&lt;14&gt; and the approaching side of the disc is towards the knot and the southern jet of B2 0258 + 35 . The radio continuum in our case is asymmetric with the southern knot; the southern jet is considerably brighter than the northern jet (see Section 3.1) and has a position angle of 132 . If the absorption GLYPH&lt;14&gt; solely arose from H i belonging to this large disc, the southern lobe would give rise to much stronger absorption compared to the northern lobe. Thus, the integrated absorption profile would be asymmetric with stronger absorption towards the blueshifted velocities. Furthermore, the absorption profile towards the southern jet (for example, at position 'a' shown in Figure 2) should have a blueshifted wing instead of a redshifted one, unlike what is observed. To produce a symmetric absorption profile, the absorbing gas and hence also the radio jets should be along the

minor axis of the rotating disc, which is not the case here (see Figure 1).

This can be illustrated by modelling the H i absorption arising from an H i disc around the radio source, in this case the galactic H i disc. The parameters involved in this modelling are the rotation curve of the disc, inclination and position angles, and the velocity dispersion of the gas.

In the nuclear regions of the galaxy, the major contribution to the rotation curve comes from the stellar component of the galaxy and the SMBH. The SMBH mass of NGC 1167 is 4.4 GLYPH&lt;2&gt; 10 8 M (Kormendy &amp; Ho 2013). For a stellar velocity disGLYPH&lt;12&gt; persion of 219.6 km s GLYPH&lt;0&gt; 1 (Ho et al. 2009) the black hole sphere of influence extends up to GLYPH&lt;24&gt; 40 pc. To estimate the contribution to rotation by the stellar component, we used the Hernquist (1990) model with an e GLYPH&lt;11&gt; ective radius of 6.7 00 and a total stellar mass of 1.3 GLYPH&lt;2&gt; 10 13 M as parameters (Noordermeer &amp; van der Hulst GLYPH&lt;12&gt; 2007). The resulting rotation curve in the inner 1 kpc region of NGC 1167 is shown in Figure 6, left panel. We note that the rotation curve is essentially flat at 300 km s GLYPH&lt;0&gt; 1 even in the very central region of the galaxy.

The inclination (38 ) and the position angles (75 GLYPH&lt;14&gt; GLYPH&lt;14&gt; ) of the disc have been well constrained by Struve et al. (2010). In addition, we assume a velocity dispersion ( GLYPH&lt;27&gt; ) of 20 km s GLYPH&lt;0&gt; 1 , a typical value observed in the central regions of early-type galaxies. However, Struve et al. (2010) have found that a model with a velocity dispersion of 9 km s GLYPH&lt;0&gt; 1 explains the H i rotation curve of NGC 1167 very well. Hence, our assumption of 20 km s GLYPH&lt;0&gt; 1 for the velocity dispersion is quite conservative. The resulting absorption profile is shown in Figure. 6, right panel. The profile is asymmetric (as explained earlier), unlike the observed profile and does not match the observed location of the peak absorption. This supports our hypothesis that the absorbing H i has kinematics that is quite di GLYPH&lt;11&gt; erent from the large H i disc.

An extreme possibility that can give rise to an absorption profile centred at the systemic velocity is that the gas does not have any rotational motion, but is entirely chaotic. Another possibility that can lead to a symmetric H i absorption profile is that the disc has undergone a warp in the central few kiloparsecs of the galaxy with an orientation such that the kinematical minor axis in the inner region coincides with the jet axis. In this case all absorption should appear close to the systemic velocity and the large width of the absorption would also imply a large dispersion of H i .

However, the resolved H i absorption spectra shown in Figure 2 do show a sign of rotation with the absorption centroid shifting redwards as we move towards the northern jet. Thus, we argue that the absorption arises from a kiloparsec-scale gas disc with a velocity dispersion higher than that of the large H i disc (further supported by the nature of CO emission; see the following subsection) whose orientation is di GLYPH&lt;11&gt; erent from that of the large H i disc.

## 4.2. Comparison with molecular gas

In order to further strengthen this hypothesis about the morphology and kinematics of the absorbing H , we make use of the i available information on molecular gas. Observed in emission, the distribution of gas traced by the molecular component is independent of the size and asymmetry of the background continuum. Thus, when combined with the absorption studies, it can provide further insights into the kinematics and morphology of the gas under investigation.

CO(1-0) emission has been detected from NGC 1167 by many studies: Prandoni et al. (2007), O'Sullivan et al. (2015)

and Bolatto et al. (2017). The first two studies have made use of IRAM single-dish observations, while the last one is an interferometric study with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The single-dish observations detect CO(1-0) emission from the central GLYPH&lt;24&gt; 8 kpc (23 00 ; the resolution of the telescope) region of the galaxy. The interferometric observation of Bolatto et al. (2017) has a larger field of view, but is of poorer sensitivity. They report the detection of very faint CO emission from the region outside the IRAM field of view, but are not sensitive to the emission seen by IRAM in the central part of the galaxy. Furthermore, the CO detected by Bolatto et al. (2017) is coincident with the low-level star formation sites in NGC 1167 reported by Gomes et al. (2016b), and thus we only consider the deeper IRAM observations of CO in the central region for our discussion.

The IRAM emission spectra presented by Prandoni et al. (2007) and O'Sullivan et al. (2015) agree with each other in their line widths (FWHM of 261 GLYPH&lt;6&gt; 36 km s GLYPH&lt;0&gt; 1 and 233 GLYPH&lt;6&gt; 47 km s GLYPH&lt;0&gt; 1 respectively) and in the integrated emission: 0.93 GLYPH&lt;6&gt; 0.14 K km s GLYPH&lt;0&gt; 1 and 1.18 GLYPH&lt;6&gt; 0.1 K km s GLYPH&lt;0&gt; 1 respectively (shown in Figure 7, bottom panel).

The integrated CO(1-0) emission profile is shown in the bottom panel of Figure 7. The profile is inverted and plotted in arbitrary units for ease of comparison with the H i absorption profile. As mentioned earlier, the rotation curve for the galaxy is approximately flat at 300 km s GLYPH&lt;0&gt; 1 even within the inner kiloparsec region (Figure 6, left panel). Thus, we would expect the CO emission profile to be double-horned if it formed a regularly rotating disc. This is not the case, which suggests that the kinematics of CO is not dominated by rotation. This supports our earlier hypothesis of the presence of a kiloparsec-scale turbulent gas disc. Assuming the gas to be entirely chaotic, we find an upper limit to the velocity dispersion of CO to be GLYPH&lt;24&gt; 90 km s GLYPH&lt;0&gt; 1 , estimated via a single-component Gaussian fit to the emission profile of O'Sullivan et al. (2015), which has a higher signal-to-noise ratio compared to that of Prandoni et al. (2007).

Next, we find (as can be seen in Figure 7) that the line widths of the CO emission and the H i absorption profiles match, as do the natures of the two profiles. This suggests that these two components of gas (molecular and atomic) are spatially coincident, and hence CO is distributed within the central few kiloparsecs of the host galaxy. Thus, this spatial coincidence suggests that the H seen in absorption is also dominated by turbulence rather than i bulk rotation. In addition, the covering factor of the gas is unity since we detect H i absorption against the entire radio source. We thus posit that the gas seen in H i absorption and CO emission forms a kiloparsec-scale disc into which the radio jets are expanding and injecting turbulence, causing the observed line profiles.

The line ratio CO(2-1) CO(1-0) can also be used to investi-/ gate the presence of turbulence. This line ratio, as reported by O'Sullivan et al. (2015), is unity, which is twice the expected value (of 0.5) for quiescent molecular gas in normal ISM (e.g. Peñaloza et al. 2017) and similar to the value observed in other cases of strong jet-ISM interaction (for example IC 5063; Oosterloo et al. 2017). This further supports the presence of strong turbulence in the medium.

The presence of turbulence has been observed in other objects for example, in 3C326N (Guillard et al. 2012) and in NGC1266 (Alatalo et al. 2015). Following a similar approach to that used for these sources, we estimated the turbulent component of energy associated with the gas. The turbulent kinetic energy is given by Eturb = 3 2 Mmol / GLYPH&lt;27&gt; 2 (e.g. Guillard et al. 2012). The molecular gas is likely to be more massive than the atomic

Fig. 6. Left panel: Rotation curve for NGC 1167 in the inner 1 kpc. In the inner regions the major contributions to the rotation curve come from the stellar component and the supermassive black hole. Right panel: Observed absorption profile (blue) and the model absorption profile (in red) for the large-scale H i disc.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-042f24ec-c329-4235-bc61-7422a180de63.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-d99da788-16a8-4eb6-92a0-3ec28202a86d.jpg)

gas in the nuclear regions and hence will contain most of the turbulent energy. We estimate this energy associated with the molecular phase by assuming a velocity dispersion ( GLYPH&lt;27&gt; ) of 90 km s GLYPH&lt;0&gt; 1 . We would like to note here that the value of velocity dispersion was derived assuming the gas to be entirely chaotic with no systematic motion. The turbulent energy estimated in this way will thus be an upper limit. To estimate the molecular gas mass, we assume a range of values for the CO to H2 conversion factor (XCO): 0.2 GLYPH&lt;2&gt; 10 20 (K km s GLYPH&lt;0&gt; 1 ) GLYPH&lt;0&gt; 1 , a conservative value representing optically thin turbulent gas; 0.4 GLYPH&lt;2&gt; 10 20 (K km s GLYPH&lt;0&gt; 1 ) GLYPH&lt;0&gt; 1 representing dense gas in ultra-luminous infrared galaxies; and 2 GLYPH&lt;2&gt; 10 20 (K km s GLYPH&lt;0&gt; 1 ) GLYPH&lt;0&gt; 1 , the Galactic value (Bolatto et al. 2013). The corresponding molecular gas masses are 8.8 GLYPH&lt;2&gt; 10 6 M , 1.8 GLYPH&lt;12&gt; GLYPH&lt;2&gt; 10 7 M , and 8.8 GLYPH&lt;12&gt; GLYPH&lt;2&gt; 10 7 M . With these assumptions, we derive a turGLYPH&lt;12&gt; bulent kinetic energy (in the molecular phase) of 2.1 GLYPH&lt;2&gt; 10 54 erg, 4.3 GLYPH&lt;2&gt; 10 54 erg, and 2.1 GLYPH&lt;2&gt; 10 55 erg, respectively, corresponding to the XCO factors mentioned above.

spreads over the entire disc (e.g. Mukherjee et al. 2018a, section 3 (v)).

## 4.3. Can jets cause this turbulence?

In B2 0258 + 35, we find that both the H i absorption and the CO emission profiles suggest that the gas is turbulent. We suggest that radio jets expanding through the gas are the main cause of this turbulence (see Sections 4.1 and 4.2). The e GLYPH&lt;11&gt; ect of radiation or wind is unlikely to be significant since the optical AGN in B2 0258 + 35 is characterised by a LINER spectrum (Gomes et al. 2016a) with only weak emission lines.

For a scenario where jets are expanding in the presence of a gas disc, the numerical simulations (e.g. Mukherjee et al. 2018b) predict that when the jets are perpendicular to the plane of the disc, they break out of the ISM quickly ( . 100 kyr). They also give rise to energy bubbles, which expand laterally into the gas disc generating radial gas outflows at velocities as high as GLYPH&lt;24&gt; 500 km s GLYPH&lt;0&gt; 1 . The jets that are partially inclined to the gas disc or in the plane of the disc evolve much more slowly, and in the process also get deflected by the gas clouds. They too launch expanding bubbles in the gas disc, which evolve slowly. These simulations further predict that the mean velocity dispersion of the gas increases by several times its initial value, up to a few hundred kilometres per second, both at the point of direct interaction with the jet head and further out, as the energy bubble thus created

The jet power of B2 0258 + 35 is GLYPH&lt;24&gt; 10 44 erg s GLYPH&lt;0&gt; 1 , based on the scaling relation between cavity power (jet power) and radio luminosity (Cavagnolo et al. 2010). Over a time span of GLYPH&lt;24&gt; 0.4 Myr (the lower limit on the age of the AGN; Brienza et al. 2018), the jets are capable of depositing a total of GLYPH&lt;24&gt; 10 57 erg. Thus, the radio jets appear to be capable enough to provide the energy required to support the maximum level of turbulence even if only 1% of the total energy in the jet is converted to turbulent energy. However, we note that the jet power estimation is quite uncertain owing to the large scatter of about an order of magnitude in the scaling relation and several uncertainties in estimating the energy losses. For example, in the case of another low-luminosity radio AGN IC 5063 (Morganti et al. 2015), the jet power estimated using this prescription was shown to be an order of magnitude lower than that required to create the observed velocity dispersion in gas (Mukherjee et al. 2018b). If the jet power turns out to be higher, the e GLYPH&lt;14&gt; ciency required is even lower.

## 4.4. Impact of the radio continuum on the host galaxy

In order to understand the e GLYPH&lt;11&gt; ect of B2 0258 + 35, which is GLYPH&lt;24&gt; 1 kpc in size, on its host galaxy, we need to explore whether the impact of the radio jets is limited only to the central kpc region or extends to larger radii. Various multiwavelength studies have shown the e GLYPH&lt;11&gt; ects of radio jets can indeed span distances much larger than the size of the radio jets.

Mahony et al. (2016) found, using IFU spectroscopy of 3C 293, ionised gas outflows with linewidths &gt; 500 km s GLYPH&lt;0&gt; 1 at 12 kpc from the nucleus along the radio axis and outflows with linewidths as high as 300 km s GLYPH&lt;0&gt; 1 at distances as far as GLYPH&lt;24&gt; 3.5 kpc perpendicular to the radio axis. Similarly, in 4C 31.04 (Zovaro et al. 2019), it was found that the jet-induced shocks extend to regions four times the size of the observed radio lobes. RodríguezArdila et al. (2017) report the presence of two expanding ionised gas shells that are spatially coincident with the radio jets. Fabbiano et al. (2017, 2018) suggest, via deep Chandra imaging and spectroscopy, that the radio jets of ESO 428-G014 give rise to GLYPH&lt;24&gt; 2-3 kpc extended X-ray emission along the major axis of the radio source. In NGC 6764, the interaction between the parsec-

Fig. 7. Comparison between the H i emission profile (top panel, red: Struve et al. 2010); CO emission profiles, inverted and in arbitrary units for ease of comparison (bottom panel, blue dotted: Prandoni et al. 2007; black dashed: O'Sullivan et al. 2015); and the WSRT H i absorption profile (bottom panel, red). The CO emission profile does not match the H i emission profile, implying that the CO emission does not arise from the molecular counterpart of the large H i disc. The H i absorption and the CO emission profiles match, implying that CO and absorbing H i are spatially coincident.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-42e3e30a-6363-4e94-8722-b852d80da31e.jpg)

scale radio jets and the ISM is thought to be one of the causes of subkiloparsec-scale radio bubbles (Kharb et al. 2010). Maksym et al. (2017) and Wang et al. (2010) also find evidence, via X-ray studies, of the jet-ISM interaction shaping the inner kpc region of the host galaxies of NGC 3396 and NGC 4151, respectively.

Predictions from the numerical simulations of jets expanding at an inclination to a clumpy gas disc (e.g. Sutherland &amp; Bicknell 2007; Cielo et al. 2018; Mukherjee et al. 2018a,b) agree with these findings. They predict that the jet flow will be channelled through the disc in di GLYPH&lt;11&gt; erent directions, which can ionise the gas even in the regions far from the radio jets.

lar to the radio axis. The [O III] H / GLYPH&lt;12&gt; , [NII] H / GLYPH&lt;11&gt; , and [S II] H / GLYPH&lt;11&gt; line ratios are plotted on diagnostic diagrams introduced by Baldwin et al. (1981), here referred to as BPT diagrams in Figure 8. We see that the line ratios clearly fall in the low-ionisation nuclear emission-line region (LINER; Heckman 1980) part of the diagram. This is further confirmed by the IFU study of Gomes et al. (2015), who present a very detailed BPT diagram up to a diameter of GLYPH&lt;24&gt; 30 kpc of NGC 1167. Except for the low surface brightness star-forming spiral arms, the remaining regions clearly fall in the LINER portion of the diagram.

We can investigate the presence of such an e GLYPH&lt;11&gt; ect in B2 0258 + 35 using the long-slit spectroscopy presented by Emonts (2006). In this study, the slit is aligned with the major axis of the host galaxy (PA = 64 ), far removed from the radio GLYPH&lt;14&gt; axis (PA = 132 ). The spectra have been extracted at the location GLYPH&lt;14&gt; of the nucleus and at the o GLYPH&lt;11&gt; -nuclear regions 4 00 south-east and north-west of the nucleus. Thus the o GLYPH&lt;11&gt; -nuclear spectra are extracted at a distance of GLYPH&lt;24&gt; 1.3 kpc from the nucleus, perpendicu-

LINER spectra can be explained by various mechanisms such as photoionisation either by nonstellar UV continuum or shocks, or evolved stars. From the IFU analysis, Gomes et al. (2016a) identify two regions in the disc of NGC 1167 (see their Figure 1): the inner region ( GLYPH&lt;24&gt; 4 00 ; 1.4 kpc) is characterised by high equivalent width H GLYPH&lt;11&gt; emission, suggesting that evolved stars alone cannot be the source of ionisation, and the outer region with a lower equivalent width where the LINER spectrum can be explained by pAGB stars.

Fig. 8. BPT diagrams showing the line ratios (Emonts 2006) of the o GLYPH&lt;11&gt; -nuclear regions of B2 0258 + 35 . The blue stars represent the line ratios from the spectrum extracted from the south-western region, o GLYPH&lt;11&gt; set by GLYPH&lt;24&gt; 1.4 kpc from the nucleus, with PA = 64 . The red dots represent the line GLYPH&lt;14&gt; ratios from the spectrum extracted at a similar distance o GLYPH&lt;11&gt; the nucleus in the north-eastern region. The points clearly lie in the LINER portion of the diagram, suggesting the presence of shocks.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-43b30619-eab2-4e86-b7ab-8a88b2c7f90c.jpg)

Thus, wide-spread shocks are the most plausible cause of the LINER spectra in the inner region. A major merger cannot be the reason for these shocks because such a possibility has already been ruled out by Struve et al. (2010). We thus suggest that they are a result of the interaction between the radio jets and the ISM as described earlier.

It is more challenging to ascertain whether this interaction has an e GLYPH&lt;11&gt; ect on the rest of the disc, outside the inner few kpc. Interestingly, it is worth noting here that in a study of group dominant galaxies, O'Sullivan et al. (2015) have found that NGC 1167 lies in the red sequence of quenched galaxies. However, with the present data, it is not possible to ascertain whether this could have, at least partly, resulted from the AGN activity. Highresolution mapping of molecular and ionised gas may shed further light in this direction.

## 4.5. Impact of the host galaxy on the radio continuum

As discussed in the earlier sections, the radio jets in B2 0258 + 35 appear to be expanding into a gas disc, injecting a significant amount of turbulence into the disc. If this is the case, numerical simulations (e.g. Mukherjee et al. 2018a) predict that the ISM strongly influences the overall morphology of the radio continuum. A strong jet-ISM interaction, in addition to generating the shocks described above, would also launch subrelativistic winds perpendicular to the plane of the galaxy. These winds resemble the wide-angle outflows seen associated with AGN winds or nuclear starburst activities (see Figure 13 in Mukherjee et al. 2018b), but are driven by the pressure of escaping jet plasma. Buoyancy further contributes to the propagation and expansion of the bubbles to large radii. The parameters of the jet-disc simulations presented in Mukherjee et al. 2018a are not fine-tuned to be an accurate representation of B2 0258 + 35. Nevertheless, simulation D in that paper, representing 10 45 erg s GLYPH&lt;0&gt; 1 jets propagating at 45 GLYPH&lt;14&gt; into an inhomogeneous disc does reproduce several of the features of B2 0258 + 35.

In order to facilitate the comparison, we generate a synthetic radio surface brightness image from simulation D, using the ap- proach described in Bicknell et al. (2018); this approach was also used to simulate the jet-ISM interaction in IC 5063 (Mukherjee et al. 2018b). The main aim of this exercise is to illustrate the e GLYPH&lt;11&gt; ect of ISM on the radio continuum and investigate whether the low surface brightness radio lobes can originate from such jet-ISM interactions. The synthetic radio image at 1.4 GHz was obtained at a time of 9 8 : GLYPH&lt;2&gt; 10 5 yr, chosen so that the extent of the radio jets matches the observed size of the radio jets in B2 0258 + 35. The radio structure thus obtained is a result of radio plasma expanding into a turbulent, inhomogeneous disc of gas and is shown in the right panel of Figure 9.

There is a remarkable similarity between the morphology of the radio jets obtained from the simulations (see the yellow contours in Figure 9, left panel) and B2 0258 + 35. In the synthetic image, the northern jet is clearly disrupted by its interaction with the ISM, but maintains a generally straight path; the southern jet has a sharp bend, which is the result of its path being obstructed by a dense cloud. There is also a hot spot in the jet approximately halfway between the core and the deflection point. The di GLYPH&lt;11&gt; erence in the jet morphologies is the result of the fractal, lognormal distribution of density in the ISM, leading to an approximately straight path in one jet, but to a substantial obstruction in the other. There is, of course, a di GLYPH&lt;11&gt; erence between the observation and the synthetic image. In the source, the southern jet is deflected to the east; in the simulation, the southern jet is deflected to the west. This feature may be the result of the random distribution of dense clumps or the gravitational field of the galaxy and requires attention in more detailed simulations. Nevertheless, we conclude that the asymmetry in the radio continuum is a result of a strong jet interaction with an inhomogeneous disc.

As noted earlier, this interaction also results in kiloparsecscale bubbles whose overall direction is di GLYPH&lt;11&gt; erent from the jets and is determined by buoyancy of the jet plasma in the gravitational field of the galaxy. However, the kiloparsec-scale simulations do not follow the evolution of the jet lobes and the bubble beyond 4 kpc, so that a comparison with the large, low surface brightness radio structures on 100 kpc scales observed in B2 0258 + 35 (see Fig. 1) is not possible at this time. Moreover, as can be seen in Figure 10, the subrelativistic winds have only spread out to a diameter of 6 kpc in this duration, while the radio lobes in question are 240 kpc in size. Hence, the faint radio structures produced by the simulations with the present jet, galaxy, and ISM parameters do not explain the large radio structures. As shown in Brienza et al. (2018), these lobes do not show spectral steepening, and it has been suggested that they represent a previous episode of activity that switched o GLYPH&lt;11&gt; no more than a few tens of Myr ago, or that they are the result of large-scale jets that have been temporarily disrupted and smothered, but that are still fuelled, albeit at a low rate, by the nuclear engine (perhaps similar to what is happening in Centaurus A; e.g. McKinley et al. 2018; Morganti et al. 1999). In this scenario, the kiloparsec-scale bubbles in the model could be the connecting structure keeping the fuelling channel open. This is schematically represented in Figure 9, left panel.

## 5. Summary

We have presented high-resolution VLA-A array and EVN H i absorption observations of B2 0258 + 35, a low-power radio AGN in a gas-rich host galaxy NGC 1167, which also hosts large (240 kpc) radio lobes. With our VLA observations, the continuum source and the absorption, is slightly resolved. We detect absorption all against the source, and the profile also exhibits a slight asymmetry towards the blueshifted velocities. At VLBI scales

Fig. 9. Left panel: Cartoon of the possible scenario explaining the observed properties of B2 0258 + 35 and the host galaxy NGC 1167. The host galaxy is viewed edge-on. Right panel: Synthetic radio (1.4 GHz) surface brightness map from the simulation of the jet entering an edge-on gas disc at an angle of 45 . The regions in yellow contours indicate the bright radio jets, while the blue contours are the di GLYPH&lt;14&gt; GLYPH&lt;11&gt; use radio emission produced by the plasma percolating through the gas disc. The image has been obtained at an arbitrary instance in the simulation for qualitative comparison. The remarkable similarity should be noted between the radio jets in this map (yellow contours) and the morphology of B2 0258 + 35 (Figure 1 inset) in the absence of any fine-tuning. This further supports the proposed strong jet-ISM interaction in B2 0258 + 35.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-4c1cda2e-15bf-4fad-805f-d198917d5bd0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00374v1-figure-e48b4387-e1f3-43b8-9471-772224902a7a.jpg)

most of the continuum is resolved out, and we detect only GLYPH&lt;24&gt; 10% of the total continuum seen with the VLA. We also detect a weak absorption feature in the integrated VLBI spectrum of the bright knot seen along the southern jet and the estimated H i column density is in agreement with that estimated from the VLA observations. This implies a uniform screen of gas in front of the background radio source.

with an inhomogeneous disc are likely to be the cause of the widespread shocks, quenched star formation, and low surface brightness radio bubbles resulting from plasma leaking out in the direction perpendicular to the plane of the galaxy, possibly over multiple episodes of AGN activity. High-resolution mapping of ionised and molecular gas and targeted simulations can shed further light on the jet-ISM interaction in B2 0258 + 35.

We modelled the H i 21cm absorption and find that the kinematics of the absorbing gas should be di GLYPH&lt;11&gt; erent from the large H i disc. CO(1-0) emission was also detected from the host galaxy NGC 1167 and the emission profile agrees well with the integrated H i absorption profile suggesting that CO and H i are cospatial and are being a GLYPH&lt;11&gt; ected by the radio jets to a similar extent. The line ratio CO(2-1) CO(1-0) / GLYPH&lt;25&gt; 1 further supports this hypothesis. We conclude that the radio source is interacting with the ambient gas, causing the observed turbulence, which has also a GLYPH&lt;11&gt; ected the morphology of the radio source. Assuming the kinematics of molecular gas to be entirely chaotic, we find that the upper limit to the turbulent kinetic energy in the medium is GLYPH&lt;24&gt; 10 55 erg. We also find that the radio jets, in this case, are capable of depositing this amount of turbulence in the medium within 0.4 Myr, the estimated age of the radio source. The optical line ratios ([O III] H / GLYPH&lt;12&gt; , [NII] H / GLYPH&lt;11&gt; and [S II] H / GLYPH&lt;11&gt; ) in the regions a few kiloparsecs from the radio nucleus suggest that the emission there arises from shocks caused by the jet-ISM interaction. All these results support the fact that low-luminosity radio AGNs are capable of a GLYPH&lt;11&gt; ecting the ISM significantly and hence contribute to negative feedback in their host galaxies.

We further compare our results with the predictions from numerical simulations and conclude that the interactions of the jets

Acknowledgements. We would like to thank the referee for the useful comments. We also thank Bjorn Emonts for the help with the optical spectra of B2 0258 + 35, and Ewan O'Sullivan and Françoise Combes for providing their CO emission spectrum. SM would like to thank Aparajitha Ramesh, Anqi Li, and Francesco Santoro for the useful discussions. Part of the research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP 2007-2013) / / ERC Advanced Grant RADIOLIFE-320745.

## References

Adebahr, B., Brienza, M., &amp; Morganti, R. 2019, A&amp;A, 622, A209 Alatalo, K., Lacy, M., Lanz, L., et al. 2015, ApJ, 798, 31 Baldwin, J. A., Phillips, M. M., &amp; Terlevich, R. 1981, PASP, 93, 5 Best, P. N., Kau GLYPH&lt;11&gt; mann, G., Heckman, T. M., &amp; Ivezi´ c, Ž. 2005, MNRAS, 362, 9 Bicknell, G. V., Mukherjee, D., Wagner, A. Y., Sutherland, R. S., &amp; Nesvadba, N. P. H. 2018, MNRAS, 475, 3493 Bolatto, A. D., Wolfire, M., &amp; Leroy, A. K. 2013, ARA&amp;A, 51, 207 Bolatto, A. D., Wong, T., Utomo, D., et al. 2017, ApJ, 846, 159 Brienza, M., Morganti, R., Murgia, M., et al. 2018, A&amp;A, 618, A45 Cavagnolo, K. W., McNamara, B. R., Nulsen, P. E. J., et al. 2010, ApJ, 720, 1066 Cavagnolo, K. W., McNamara, B. R., Wise, M. W., et al. 2011, ApJ, 732, 71 Cielo, S., Bieri, R., V olonteri, M., Wagner, A. Y ., &amp; Dubois, Y . 2018, MNRAS, 477, 1336

Croston, J. H., Hardcastle, M. J., Kharb, P., Kraft, R. P., &amp; Hota, A. 2008, ApJ, 688, 190

Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 11 Emonts, B. H. C. 2006, PhD thesis, University of Groningen Fabbiano, G., Elvis, M., Paggi, A., et al. 2017, ApJ, 842, L4 Fabbiano, G., Paggi, A., Karovska, M., et al. 2018, ApJ, 865, 83 Fabian, A. C. 2012, ARA&amp;A, 50, 455 Finlez, C., Nagar, N. M., Storchi-Bergmann, T., et al. 2018, MNRAS, 479, 3892 Giroletti, M., Giovannini, G., &amp; Taylor, G. B. 2005, A&amp;A, 441, 89 Gomes, J. M., Papaderos, P., Kehrig, C., Vílchez, J. M., &amp; Lehnert, M. D. 2015, in IAU Symposium, Vol. 309, Galaxies in 3D across the Universe, ed. B. L. Ziegler, F. Combes, H. Dannerbauer, &amp; M. Verdugo, 105-108 Gomes, J. M., Papaderos, P., Kehrig, C., et al. 2016a, A&amp;A, 588, A68 Gomes, J. M., Papaderos, P., Vílchez, J. M., et al. 2016b, A&amp;A, 585, A92 Guillard, P., Boulanger, F., Pineau des Forêts, G., et al. 2012, ApJ, 749, 158 Harrison, C. M., Costa, T., Tadhunter, C. N., et al. 2018, Nature Astronomy, 2, 198 Heckman, T. M. 1980, A&amp;A, 87, 152 Hernquist, L. 1990, ApJ, 356, 359 Ho, L. C., Filippenko, A. V., &amp; Sargent, W. L. W. 1997, ApJS, 112, 315 Ho, L. C., Greene, J. E., Filippenko, A. V ., &amp; Sargent, W. L. W. 2009, ApJS, 183, 1

Kharb, P., Lister, M. L., &amp; Cooper, N. J. 2010, ApJ, 710, 764 Kormendy, J. &amp; Ho, L. C. 2013, ARA&amp;A, 51, 511 Mahony, E. K., Oonk, J. B. R., Morganti, R., et al. 2016, MNRAS, 455, 2453 Maksym, W. P., Fabbiano, G., Elvis, M., et al. 2017, ApJ, 844, 69 Maksym, W. P., Fabbiano, G., Elvis, M., et al. 2018, arXiv e-prints [ arXiv:1810.12926 ] McKinley, B., Tingay, S. J., Carretti, E., et al. 2018, MNRAS, 474, 4056 McNamara, B. R., Kazemzadeh, F., Ra GLYPH&lt;11&gt; erty, D. A., et al. 2009, ApJ, 698, 594 McNamara, B. R. &amp; Nulsen, P. E. J. 2012, New Journal of Physics, 14, 055023 Morganti, R., Killeen, N. E. B., Ekers, R. D., &amp; Oosterloo, T. A. 1999, MNRAS, 307, 750 Morganti, R., Oosterloo, T., Oonk, J. B. R., Frieswijk, W., &amp; Tadhunter, C. 2015, A&amp;A, 580, A1 Morganti, R., Oosterloo, T., &amp; Tsvetanov, Z. 1998, AJ, 115, 915 Mukherjee, D., Bicknell, G. V., Wagner, A. Y., Sutherland, R. S., &amp; Silk, J. 2018a, MNRAS, 479, 5544 Mukherjee, D., Wagner, A. Y., Bicknell, G. V., et al. 2018b, MNRAS, 476, 80 Murgia, M. 2003, PASA, 20, 19 Noordermeer, E. &amp; van der Hulst, J. M. 2007, MNRAS, 376, 1480 O'Dea, C. P. 1998, PASP, 110, 493 Oosterloo, T., Raymond Oonk, J. B., Morganti, R., et al. 2017, A&amp;A, 608, A38 Orienti, M. &amp; Dallacasa, D. 2014, MNRAS, 438, 463 O'Sullivan, E., Combes, F., Hamer, S., et al. 2015, A&amp;A, 573, A111 Peñaloza, C. H., Clark, P. C., Glover, S. C. O., Shetty, R., &amp; Klessen, R. S. 2017, MNRAS, 465, 2277 Perley, R. A. &amp; Butler, B. J. 2017, ApJS, 230, 7 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014, A&amp;A, 571, A16 Prandoni, I., Laing, R. A., Parma, P., et al. 2007, New A Rev., 51, 43 Randall, S. W., Forman, W. R., Giacintucci, S., et al. 2011, ApJ, 726, 86 Rodríguez-Ardila, A., Prieto, M. A., Mazzalay, X., et al. 2017, MNRAS, 470, 2845 Sanghera, H. S., Saikia, D. J., Luedke, E., et al. 1995, A&amp;A, 295, 629 Shulevski, A., Morganti, R., Oosterloo, T., &amp; Struve, C. 2012, A&amp;A, 545, A91 Silk, J. &amp; Rees, M. J. 1998, A&amp;A, 331, L1 Struve, C., Oosterloo, T., Sancisi, R., Morganti, R., &amp; Emonts, B. H. C. 2010, A&amp;A, 523, A75 Sutherland, R. S. &amp; Bicknell, G. V. 2007, ApJS, 173, 37 Wang, J., Fabbiano, G., Elvis, M., et al. 2011, ApJ, 736, 62 Wang, J., Fabbiano, G., Risaliti, G., et al. 2010, ApJ, 719, L208 Wegner, G., Haynes, M. P., &amp; Giovanelli, R. 1993, AJ, 105, 1251 Zovaro, H. R. M., Sharp, R., Nesvadba, N. P. H., et al. 2019, MNRAS, 484, 3393
Suma Murthy Raffaella Morganti Tom Oosterloo Robert Schulz Dipanjan Mukherjee Alexander Y. Wagner Geoffrey Bicknell Isabella Prandoni Aleksander Shulevski
Feedback from low-luminosity radio galaxies: B2 0258+35
Low-luminosity radio-loud active galactic nuclei (AGN) are of importance in
studies concerning feedback from radio AGN since a dominant fraction of AGN
belong to this class. We report high-resolution Very Large Array (VLA) and
European VLBI Network (EVN) observations of HI-21cm absorption from a young,
compact steep-spectrum radio source, B2 0258+35, nested in the early-type
galaxy NGC 1167, which contains a 160 kpc HI disc. Our VLA and EVN HI
absorption observations, modelling, and comparison with molecular gas data
suggest that the cold gas in the centre of NGC 1167 is very turbulent (with a
velocity dispersion of ~ 90 km/s) and that this turbulence is induced by the
interaction of the jets with the interstellar medium (ISM). Furthermore, the
ionised gas in the galaxy shows evidence of shock heating at a few kpc from the
radio source. These findings support the results from numerical simulations of
radio jets expanding into a clumpy gas disc, which predict that the radio jets
in this case percolate through the gas disc and drive shocks into the ISM at
distances much larger than their physical extent. These results expand the
number of low-luminosity radio sources found to impact the surrounding medium,
thereby highlighting the possible relevance of these AGN for feedback.