## Complex molecular gas kinematics in the inner 5 kpc of 4C12.50 as seen by ALMA

C. M. Fotopoulou 1 2 ; , K. M. Dasyra 1 3 ; , F. Combes 4 , P. Salomé 4 and M. Papachristou 1 3 ;

- 1 Department of Astrophysics, Astronomy &amp; Mechanics, Faculty of Physics, National and Kapodistrian University of Athens, Panepistimiopolis Zografos 15784, Greece
- 2 Max-Planck Institut für Astrophysik, Karl-Schwarzschild-Str. 1, D-85741 Garching, Germany
- 3 Institute for Astronomy, Astrophysics, Space Applications, and Remote Sensing, National Observatory of Athens, 15236 Penteli, Greece
- 4 Observatoire de Paris, LERMA, 61 Av. de l'Observatoire, 75014, Paris, France

## ABSTRACT

The nearby system 4C12.50, also known as IRAS 13451 + 1217 and PKS 1345 + 12, is a merger of gas-rich galaxies with infrared and radio activity. It has a perturbed interstellar medium (ISM) and a dense configuration of gas and dust around the nucleus. The radio emission at small ( GLYPH&lt;24&gt; 100 pc) and large ( GLYPH&lt;24&gt; 100 kpc) scales, as well as the large X-ray cavity in which the system is embedded, are indicative of a jet that could have a GLYPH&lt;11&gt; ected the ISM. We carried out observations of the CO(1-0), (3-2), and (4-3) lines with the Atacama Large Millimeter Array (ALMA) to determine basic properties (i.e., extent, mass, and excitation) of the cold molecular gas in this system, including its already-known wind. The CO emission reveals the presence of gaseous streams related to the merger, which result in a small ( GLYPH&lt;24&gt; 4kpc-wide) disk around the western nucleus. The disk reaches a rotational velocity of 200 kms 1 , and has GLYPH&lt;0&gt; a mass of 3.8( GLYPH&lt;6&gt; 0.4) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; . It is truncated at a gaseous ridge north of the nucleus that is bright in [O iii ]. Regions with high-velocity COemission are seen at signal-to-noise ratios of between 3 and 5 along filaments that radially extend from the nucleus to the ridge and that are bright in [O iii ] and stellar emission. A tentative wind detection is also reported in the nucleus and in the disk. The molecular gas speed could be as high as 2200 km s GLYPH&lt;0&gt; 1 and the total wind mass could be as high as 1.5( GLYPH&lt;6&gt; 0.1) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; . Energetically, it is possible that the jet, assisted by the radiation pressure of the active nucleus or the stars, accelerated clouds inside an expanding bubble.

Key words. ISM: jets and outflows - ISM: kinematics and dynamics - Line: profiles - Galaxies: active - Galaxies: nuclei Infrared: galaxies

## 1. Introduction

Galactic winds that are driven by the feedback of active galactic nuclei (AGNs) or young stars and that are detected in molecular gas tracers are now considered common (Sakamoto et al. 2006, 2009; Leon et al. 2007; Feruglio et al. 2010; Fischer et al. 2010; Fluetsch et al. 2019; Alatalo et al. 2011; Rangwala et al. 2011; Sturm et al. 2011; Krips et al. 2011; Dasyra &amp; Combes 2011, 2012; Aalto et al. 2012; Tsai et al. 2012; Morganti et al. 2013b; Combes et al. 2013; Spoon et al. 2013; Veilleux et al. 2013; Cicone et al. 2014; García-Burillo et al. 2014; George et al. 2014; Sakamoto et al. 2014; Alatalo et al. 2014; Tombesi et al. 2015; Aalto et al. 2016; Stone et al. 2016; González-Alfonso et al. 2017; Pereira-Santaella et al. 2018). In galaxies with active nuclei, the momentum rate of the molecular winds is often considerably higher, that is, about 20 times higher, than the pressure exerted by the AGN radiation (e.g., Cicone et al. 2014; Carniani et al. 2015). Multiple photon scatterings (Ishibashi &amp; Fabian 2015) and, more frequently, an energy-conserving expansion of an ionized gas bubble that leads to momentum boosting (King et al. 2011; Faucher-Giguère &amp; Quataert 2012; Zubovas &amp; King 2014) have been evoked to justify the high momentum rates of molecular winds. Radiation pressure can drive such winds during specific phases of the ionized medium expansion: when the expansion happens so rapidly compared to the radiative cooling that it is nearly adiabatic (as in the Sedov-Taylor phase of supernovae). Radio jets, when powerful or nearly relativistic, can e GLYPH&lt;14&gt; ciently drive adiabatically expanding bubbles as they rapidly deposit energy in the interstellar medium (ISM) for most expansion phases (Wagner et al. 2016). Indeed, several of the molecular winds in the above-mentioned studies were detected in galaxies with AGN jets.

Dasyra et al. (2016) and Oosterloo et al. (2017) examined what happens in the molecular gas when a radio jet impacts clouds. Focusing on the nearby galaxy IC5063, they found that the flux ratios of CO lines reveal the presence of highly excited and optically thin gas in the wind. This result indicated that jetISM interactions can leave traces in the kinematics and or in the / excitation of the impacted gas. In previous work, Dasyra et al. (2014) reported heating of the molecular gas in the wind of another radio galaxy, 4C12.50. This result emerged from the comparison of two molecular gas probes: the fraction of accelerated cold ( &lt; 25K) gas in CO data was less than one third of the total reservoir, which corresponded to the fraction of accelerated warm ( GLYPH&lt;24&gt; 400K) gas in H2 data (from the Spitzer Space Telescope; Dasyra et al. 2014).

The system 4C12.50, also known as IRAS 13451 + 1217 and PKS1345 + 12, is a good candidate for further studies of the impact of AGN feedback on the ISM. This system is an ultraluminous infrared galaxy (ULIRG) in the local Universe, which originated from the merger of other galaxies (Fig. 1). It has two nuclei 4.4 kpc away: a western nucleus, from which the radio emission emerges, and an eastern nucleus. We refer to these nuclei as main or primary, and secondary, respectively. In optical wavelengths,

the continuum emission of the primary nucleus is comparable to (i.e., lower by a factor of 1.3 than) that of the secondary nucleus within a radius of 1 00 . The K- and L-band images with Subaru show that the primary galaxy is roughly twice as massive (Imanishi &amp; Saito 2014). Contrarily, the primary nucleus has a significantly greater [O iii ] emission than the secondary nucleus (i.e., greater by a factor of three within the same radius), as derived from the subtraction of two Hubble Space Telescope images (Fig. 1). This indicates that the gas transfer toward the main nucleus has significantly progressed. The [O iii ] emission reveals nonregular structures, such as a filamentary ridge to the northwest of the main nucleus. Traces of this ridge can also be observed in the continuum emission, indicating the presence of stars. Several shells and tidal tails that are caused by the merger are seen at distances greater than 2 00 . Some of them are visualized in another optical, H GLYPH&lt;11&gt;; and continuum image, from which we subtracted a model of one bulge and two disks to bring up underlying structures (Dasyra et al. 2011). Further larger-scale tails are seen by Emonts et al. (2016). Stanghellini et al. (1993) suggested that the merging system is part of a poor cluster. A third, smaller galaxy is seen 18 kpc northwest of the main nucleus, which could potentially be part of a merging group.

The two unambiguously merging nuclei are embedded in an X-ray bubble of about 30 GLYPH&lt;2&gt; 15 kpc (Fig. 2; Siemiginowska et al. 2008). Another, discrete X-ray component is detected along the radio jet axis further to the south. Radio emission is detected in two very di GLYPH&lt;11&gt; erent scales. A superluminal jet is seen within 100 pc from the main nucleus at 2 to 6 cm wavelengths (Stanghellini et al. 1997; Lister et al. 2003). More di GLYPH&lt;11&gt; use emission, extending out to 100 kpc, is seen at 1.36 and 1.66 GHz (see Fig. 2; Stanghellini et al. 2005). This indicates the presence of a jet that restarted recently - as recently as a few thousand years (Lister et al. 2003). Tilts of the radio emission in all scales indicate that the jet is precessing and that it could impact the ISM at various locations. Indeed, an outflow of atomic gas had been seen in absorption by Morganti et al. (2004, 2013a) in clouds at the tip of the small-scale radio emission. The outflow had long been studied in its ionized phase (Holt et al. 2008, 2009, 2011). Dasyra &amp;Combes (2011), Guillard et al. (2012), and Dasyra &amp; Combes (2012) reported the outflow detection in rotational H2 lines and in CO(2 GLYPH&lt;0&gt; 3) absorption. The outflow had yet to be seen in CO emission.

To further investigate the physical conditions of the lowtemperature molecular gas in the wind of 4C12.50 (including its location, mass, and excitation), we acquired new millimeter interferometric data of CO (1 GLYPH&lt;0&gt; 0), (3 GLYPH&lt;0&gt; 2), and (4 GLYPH&lt;0&gt; 3). For simplicity, we adopted a GLYPH&lt;3&gt; CDMcosmology with H0 = 70 kms GLYPH&lt;0&gt; 1 Mpc GLYPH&lt;0&gt; 1 , GLYPH&lt;10&gt; M = 0.3, and GLYPH&lt;10&gt; GLYPH&lt;3&gt; = 0.7, which yields 2.2kpc per arcsecond and a luminosity distance of DL = 570 Mpc (Wright 2006) at the redshift of the source.

## 2. The ALMA data and their reduction

The Atacama Large Millimeter Array (ALMA) observed 4C12.50 for the cycle 2 program 2013.1.00180.S (PI Dasyra) in band 3 (July 2015), and in bands 7 and 8 (May 2015). At the time, the array had 37 operational antennas. In all bands, four spectral windows of 1.875 GHz bandwidth were employed. In band 3, we used one spectral window to cover the velocity range -3300 km s GLYPH&lt;0&gt; 1 to 2100 km s GLYPH&lt;0&gt; 1 around the CO(1-0) line, and one spectral window to cover the velocity range -3800 km s GLYPH&lt;0&gt; 1 to 1100 km s GLYPH&lt;0&gt; 1 around the HCS (3-2) line. + The other two spectral windows targeted the continuum at rest frame

2.641( GLYPH&lt;6&gt; 0.026) mm and 2.364( GLYPH&lt;6&gt; 0.021) mm. In band 7, we used two partially overlapping spectral windows to cover the velocity range -2600 km s GLYPH&lt;0&gt; 1 to 900 km s GLYPH&lt;0&gt; 1 around the CO(3-2) line, and two spectral windows to obtain the restframe 896( GLYPH&lt;6&gt; 3) GLYPH&lt;22&gt; m and 902( GLYPH&lt;6&gt; 3) GLYPH&lt;22&gt; m continuum. In band 8, we used two spectral windows to cover the velocity range -1350 km s GLYPH&lt;0&gt; 1 to 840 km s GLYPH&lt;0&gt; 1 around the CO(4-3) line, one spectral window to target the HCO (5-4) line, and one spectral window to obtain the con-+ tinuum at restframe 668.4( GLYPH&lt;6&gt; 1.7) GLYPH&lt;22&gt; m. For the line observations, we typically chose a velocity bin of 5( GLYPH&lt;6&gt; 1) km s GLYPH&lt;0&gt; 1 , which is well suited for the detection of lines originating from small cloud ensembles (either in emission or in absorption in front of a background continuum). For the continuum observations, we chose coarser velocity bins of GLYPH&lt;24&gt; 20 up to GLYPH&lt;24&gt; 90 km s GLYPH&lt;0&gt; 1 .

To create the cubes, we reran the Common Astronomy Software Applications (CASA) reduction routines delivered by ESO. The calibration routines were run with CASA version 4.3.1, 4.2.2, and 4.3.1 for bands 3, 7, and 8, respectively. The automated pipeline results were restored for bands 3 and 7. For band 8, the calibrations were fine tuned: the data from the antennas DA34 and DV02 were flagged due to high system temperatures caused by an error in the antenna position database. After a first basic flagging, the calibration pipeline computed the water vapor, system temperature, and antenna position calibration solutions and applied them to the measurement set (MS), that is, to the visibilities of all baselines, spectral windows, and target fields. It then extracted the pertinent spectral windows and continued with the masking of antenna-shadowed baselines and edge channels. The assignment of a flux model for the calibrator, and the computation of analytic expressions for the bandpass (phase and amplitude vs. frequency) and gain (phase and amplitude vs. time) calibrations were then done. The flux calibration solution was found using 3C273 in band 7, and using Titan in bands 3 and 8. Some band 8 baselines had to be excluded due to the extent of Titan. In all bands, the bandpass calibration solution was found using the QSO J1256-0547 (3C279). The phase was self-calibrated on 4C12.50 using additional observations dedicated for this purpose.

Once the calibrations were applied, the visibilities of the selected fields were converted into image cubes. For the imaging part of the pipeline, CASA version 5.4 was used. The image reconstruction was executed with the routine tclean that also performs the beam cleaning. The science-intended data were merged with the calibration-intended data of 4C12.50 during the image reconstruction phase for the sake of maximum uv plane coverage, even though the addition of the calibration-intended data had a small contribution to the total exposure time. As science target, 4C12.50 was observed for 24.4, 5.2, and 24 minutes, and as phase calibrator, 4C12.50 was observed for 2.0, 0.8, and 2.0 minutes in bands 3, 7, and 8, respectively. For the image reconstruction, we used Briggs weighting of the visibilities with robustness parameter of 2 for band 3, and 1 for bands 7 and 8. A robustness parameter of 2 corresponds to natural weighting, most appropriate for high-sensitivity detection experiments. Indeed, a robustness parameter below 2 led to considerable flux losses in band 3. In bands 7 and 8, a robustness parameter of either 1 or 2 led to identical flux measurements. In these cases, the lower value was selected, as it leads to a narrower beam with smaller sidelobes. During the image reconstruction, a primary beam correction was also applied. This procedure was repeated for continuum-free cubes, produced by the uvcontsub routine. First-order polynomials were employed for the modeling and the subtraction of the continuum in all spectral windows with

one exception: in band 8, we had to use a zero-order polynomial as the main spectral window extended from -400 km s GLYPH&lt;0&gt; 1 to 840 km s GLYPH&lt;0&gt; 1 of the line center, and only the right-hand side of the continuum could be sampled. Contrarily, in band 3 more than 1000 km s GLYPH&lt;0&gt; 1 were available on the left side of the line and a few hundred km s GLYPH&lt;0&gt; 1 on the right side of the line enabling the continuum determination. Access to the continuum flux on both sides of the line reduces the uncertainty in the continuum slope. The small number of continuum-free channels on the right side of the line and the dynamic range of the observations add to the uncertainty in the continuum level.

In all bands, the configuration of the antennas led to a subarcsecond beam. The beam was 0.62 00 GLYPH&lt;2&gt; 0.56 00 at a position angle (PA) of 14 GLYPH&lt;14&gt; in band 3, 0.61 00 GLYPH&lt;2&gt; 0.47 00 at a PA of 53 GLYPH&lt;14&gt; in band 7, and 0.43 00 GLYPH&lt;2&gt; 0.36 00 at a PA of -1 GLYPH&lt;14&gt; in band 8. At a common pixel scale of 0.08 00 and at a common velocity resolution of 20 km s GLYPH&lt;0&gt; 1 , the noise level is at 0.40 mJy beam for CO(1-0), at 1.1 mJy beam / / for CO(3-2), and at 1.2 mJy beam for CO(4-3). For the exam-/ ination of flux ratios, we produced a second set of datacubes, by convolving the data to a common beam of 0.63 00 GLYPH&lt;2&gt; 0.63 00 using the routine ia.convolve2d. In this fully homogenized dataset the noise level is at 0.41 mJy beam, 1.2 mJy beam, and 1.4 / / mJy beam for CO(1-0), CO(3-2), and CO(4-3), respectively. /

To create the deepest possible cubes for each line, we also looked into the archive for additional CO line data from observations of 4C12.50 as a calibrator for other programs. Indeed, the CO(1-0) line was observed with 38 antennas for 13.6 mins as a calibrator for the program 2013.1.00976.S. We calibrated the measurement set using the delivered calibration pipeline with CASA 4.7.0, and then we merged the calibrated measurement set with ours using the routine concat in CASA 5.4 to create a final cube with GLYPH&lt;24&gt; 40 mins of on-source exposure. As before, we used tclean to reconstruct this cube, which has a new (common) beam of 0.81 00 GLYPH&lt;2&gt; 0.74 00 at a PA of -22 . The common GLYPH&lt;14&gt; pixel scale of 0.08 00 was selected for comparison purposes. This CO(1-0) cube reached noise levels of 0.17 mJy at a spectral resolution of 100 km s GLYPH&lt;0&gt; 1 . Likewise, the CO(3-2) line was observed with 27 antennas for 5.8 mins as a calibrator for the program 2012.1.00797.S. We calibrated the data in CASA 4.1, and we merged the calibrated measurement set with ours in CASA 5.4. We again adopted the common pixel scale of 0.08 00 for comparison purposes. The beam was 0.63 00 GLYPH&lt;2&gt; 0.48 00 at a PA of 51 . GLYPH&lt;14&gt; This CO(3-2) cube reached noise levels of 1.1 mJy beam at a / spectral resolution of 20 km s GLYPH&lt;0&gt; 1 and 0.57 mJy at a spectral resolution of GLYPH&lt;24&gt; 100 kms GLYPH&lt;0&gt; 1 . No further data were found for CO(4-3). From here onwards, we refer to the concatenated low-spatialresolution cubes as deep cubes.

## 3. Results

## 3.1. Line fluxes

Spectra of all detected lines are shown in Figs. 3 and 4. The CO (1 GLYPH&lt;0&gt; 0), (3 GLYPH&lt;0&gt; 2), and (4 GLYPH&lt;0&gt; 3) lines that we detected provided an average redshift z of 0.122( GLYPH&lt;6&gt; 0.001) when fitted by Gaussian profiles. At that redshift, 1 00 corresponds to 2.2kpc. The integrated CO line fluxes from the high-spatial-resolution cubes are 13.6( GLYPH&lt;6&gt; 1.5) Jy km s GLYPH&lt;0&gt; 1 , 66( GLYPH&lt;6&gt; 2) Jy km s GLYPH&lt;0&gt; 1 and 70( GLYPH&lt;6&gt; 2) Jy km s GLYPH&lt;0&gt; 1 for CO (1 GLYPH&lt;0&gt; 0), (3 GLYPH&lt;0&gt; 2), and (4 GLYPH&lt;0&gt; 3), respectively. The integrated CO(1 GLYPH&lt;0&gt; 0) flux from the deep concatenated data, which have a high coverage of the uv-plane in short baselines, is 13.5( GLYPH&lt;6&gt; 1.2) Jy km s GLYPH&lt;0&gt; 1 . The results from the CO(1-0) data are in good agreement with the results of Evans et al. (1999) yielding a CO(1-0) flux of 14( GLYPH&lt;6&gt; 4) Jy km s GLYPH&lt;0&gt; 1 . Previous CO(1-0) interferometric ob- servations taken with the IRAM Plateau de Bure (PdB) array yielded a CO(1-0) flux of 18.2( GLYPH&lt;6&gt; 1.5) Jy km s GLYPH&lt;0&gt; 1 at a resolution of 4.0 00 GLYPH&lt;2&gt; 3.8 00 (Dasyra et al. 2014). The integrated CO(3 GLYPH&lt;0&gt; 2) flux from the deep concatenated cube is 67( GLYPH&lt;6&gt; 2) Jy km s GLYPH&lt;0&gt; 1 . For comparison, the CO(3-2) line flux is comparable to the value we had previously measured with the IRAM 30m single dish telescope: 50( GLYPH&lt;6&gt; 8) Jy km s GLYPH&lt;0&gt; 1 (Dasyra &amp; Combes 2012). This means that not only did the band 7 data not miss extended emission, but that they were more sensitive than the 30m telescope data. As in both bands the high-resolution data reach comparable flux levels to the deep and the ancillary data, we perform all of our calculations with them.

Additionally, HCO (5-4) + and a H2O line ( GLYPH&lt;23&gt; = 0 4[2,3] -3[3,0]) at restframe 448.00108 GHz were detected in band 8. The H2O emission, a shock tracer, is nearly unresolved, originating from the vicinity of the main nucleus of 4C12.50. The HCO (5-4) ion emission tracing the AGN radiation is also cir-+ cumnuclear with a low-signal extension towards the northwest. The overall fluxes of the HCO + (5-4) and H2O lines are 3.5( GLYPH&lt;6&gt; 0.2) and 2.1( GLYPH&lt;6&gt; 0.2) Jy km s GLYPH&lt;0&gt; 1 , respectively. HCS + (3 GLYPH&lt;0&gt; 2) was not detected in band 3. Unfortunately, we were unable to test the CO(2 GLYPH&lt;0&gt; 3) absorption previously reported by Dasyra &amp; Combes (2012). The concatenation of the two spectral windows of the band 7 observations at about -900 km s GLYPH&lt;0&gt; 1 (indicated in Fig. 4), together with bandpass calibration and continuum subtraction uncertainties, led to significant di GLYPH&lt;11&gt; erences in the continuum level and slope between the spectral window edges. These di GLYPH&lt;11&gt; erences can mimic absorption or emission features or remove real absorption or emission features near the concatenation region. In the nuclear spectrum of CO(3-2) for example, the continuum to the right of the line is higher than the continuum to the left of the line (Fig. 4).

## 3.2. Disk and ambient gas emission

The CO data indicate an overall gas extent that exceeds GLYPH&lt;24&gt; 12( GLYPH&lt;6&gt; 2) kpc along the east-west axis (Fig. 7). Its shape is determined by the superposition of spatially resolved structures and its total mass in H2, as inferred from the high-resolution CO(1-0) data, is 7.7( GLYPH&lt;6&gt; 0.8) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; (Fig. 3). The H2 mass is computed as in Solomon et al. (1997).

$$\colon \text{am at a } M _ { H _ { 2 } } = 3. 2 5 \times 1 0 ^ { 7 } \alpha S _ { C O ( 1 - 0 ) } \Delta V \nu _ { o b s } ^ { - 2 } D _ { L } ^ { 2 } / ( 1 + z ) \, M _ { \odot }, \text{ \quad \ } ( 1 )$$

where GLYPH&lt;11&gt; is the CO intensity to H2 mass conversion factor, SCO (1 GLYPH&lt;0&gt; 0) GLYPH&lt;1&gt; V is the integrated line flux in Jy km s GLYPH&lt;0&gt; 1 , GLYPH&lt;23&gt; obs is the observed CO(1-0) frequency in GHz, and DL is the source luminosity distance in megaparsec. Throughout our work we adopted an GLYPH&lt;11&gt; value of 0.8 M GLYPH&lt;12&gt; / (Kkms GLYPH&lt;0&gt; 1 pc 2 ) for consistency with the literature assuming a lower conversion factor for ULIRGs than for the Milky Way (Downes &amp; Solomon 1998).

A dynamically settled disk is detected around the main nucleus, from which the radio jet is launched. Its H2 mass within 0.8 00 (radius 1.8kpc) is 3.8( GLYPH&lt;6&gt; 0.4) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; . Compared to the CO(1-0) emission, the CO(3-2) or the CO(4-3) emission in the disk is more nucleated (Fig. 7). In all CO lines, the disk is centered at the main nucleus, and it follows a progression from the northeast to the southwest with increasing velocity (from negative to positive). Depending on the emission line and the galaxy side examined for the kinematics, the projected disk circular velocity is 150-200 km s GLYPH&lt;0&gt; 1 . The disk rotation pattern is best revealed by the fitting of a Gaussian line profile in every spatial pixel of the ALMA data (see Fig. 8 for the momenta maps). The rotation is clearly seen for the inner 2 kpc, and it agrees

with the results of Imanishi et al. (2016, 2018) for comparable or smaller scales. Counter-rotating blobs or high-velocity dispersion regions are also seen in the momenta maps, for example to the west of the nucleus.

A blueshifted stream from the primary nucleus reaches, in the plane of the sky, the secondary nucleus (see Fig. 7; -300 &lt; V -100kms &lt; GLYPH&lt;0&gt; 1 panel). This stream could either be a true bridge or a tail of the main nucleus in front of the secondary nucleus, similarly to the tail in projection between, for example, M51 and NGC5195 (Toomre &amp; Toomre 1972). In this scenario, the system originates from the major merger of a gas-rich galaxy with an elliptical galaxy that has a (potentially small) disk. The companion elliptical has gone through a series of turn-arounds and pericenters, creating loops and spiral extensions at each pericenter. Another scenario is that the system originates from the major merger of two gas-rich galaxies with bulges, and that gas stripping led to a ratio of secondary-disk flux over primary-disk flux that is low for the [O iii ] and even lower for the CO.

## 3.3. Galactic components (out of dynamical equilibrium) south of the primary nucleus

A tail-like structure south of the main nucleus is the brightest extra-nuclear region detected in our data (Figs. 9, 10). Its total flux, as deduced from the CO(1-0) high-resolution cube is 3.1( GLYPH&lt;6&gt; 0.6) Jy km s GLYPH&lt;0&gt; 1 , which corresponds to a mass of 1.7( GLYPH&lt;6&gt; 0.3) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; .

This structure is detected in all CO(1 GLYPH&lt;0&gt; 0), (3 GLYPH&lt;0&gt; 2), and (4 GLYPH&lt;0&gt; 3) lines and is brightest in CO(3-2). The detected emission is so bright that it dominates the kinematics of the region: its velocity dispersion is low, typical of that of a spiral arm (60-80 km s GLYPH&lt;0&gt; 1 ), while the velocity dispersion of the disk exceeds 120 km s GLYPH&lt;0&gt; 1 (Fig. 8). Its mean velocity is high (200-300 km s GLYPH&lt;0&gt; 1 ) compared to that of the disk, meaning that its redshifted component dominates.

The redshifted emission primarily originates from a structure that spirals out from the nucleus to the southwest. In its inner part, the gas reaches velocities of GLYPH&lt;24&gt; 300 km s GLYPH&lt;0&gt; 1 . Further away, towards the southwest, the gas velocity exceeds 500 kms GLYPH&lt;0&gt; 1 (Fig. 9). The structure is either a spiral arm that has conserved its initial angular momentum, or a tidal tail, or material that was accreted in a polar ring geometry during the merger. In any case, it can be related to the large-scale structure seen by Dasyra et al. (2014) at scales &gt; 4 kpc. Adding its mass to the above calculated mass of the tail-like structure, the mass of the south structure reaches 1.8( GLYPH&lt;6&gt; 0.4) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; . A secondary, fainter structure with di GLYPH&lt;11&gt; erent orientation and angular momentum vector, which spirals out from the nucleus to the southeast is also seen in the inner part of the emission. The same applies to the inner part of the blueshifted emission partly coinciding with the redshifted emission (Fig. 10). It might be yet another spiral arm or twisted tidal tail, or even part of a second bridge to the secondary nucleus. Due to the spatial overlap of the blueshifted and the redshifted components, we cannot reliably identify wind candidates or ascribe masses to individual components.

## 3.4. Other structures as wind components or candidates

## 3.4.1. Nuclear wind in CO

A tentative wind detection in the nuclear region is seen in the CO(1-0) data (Fig. 11). Components at multiple velocities are seen due to line-of-sight e GLYPH&lt;11&gt; ects and or to an intrinsic distribu-/

tion of velocities (potentially caused by the di GLYPH&lt;11&gt; erent e GLYPH&lt;14&gt; ciency of acceleration of the molecular gas at di GLYPH&lt;11&gt; erent densities). The molecular wind terminal velocity exceeds -1800 km s GLYPH&lt;0&gt; 1 for the blueshifted wing and 1400 km s GLYPH&lt;0&gt; 1 for the redshifted wing. Holt et al. (2003), in their optical spectroscopic study of 4C12.50, also found a very broad wind in the nucleus of the system. The [O iii ] line profile required three individual Gaussian components to be fit; one of them was centered at -1980 km s GLYPH&lt;0&gt; 1 , with a width of 1944 km s GLYPH&lt;0&gt; 1 . The terminal velocity of the [O iii ] was thus 2950 km s GLYPH&lt;0&gt; 1 . A Gaussian fit to the data provides the residual emission underlying the disk emission (Fig. 11), indicating a residual flux of 0.94( GLYPH&lt;6&gt; 0.12) Jy km s GLYPH&lt;0&gt; 1 . This corresponds to a nuclear wind mass of 5.3( GLYPH&lt;6&gt; 0.7) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; .

Despite the strong continuum and the high dynamic range (of order 1000; Fig. 4), the wind detection survives the test of altering pipeline parameter values (e.g., first- and second-order continuum fit, di GLYPH&lt;11&gt; erent spectral window ranges for the continuum fit) in both datasets, that is, in our high-resolution and in the deep cubes. In previous work carried out by the team using Plateau de Bure (PbD) data (Dasyra et al. 2014; Dasyra &amp; Combes 2012) this nuclear emission was not detected because the spectral window employed was not su GLYPH&lt;14&gt; cient (the PdB data reached -1500 km s GLYPH&lt;0&gt; 1 , whereas the detection in the current ALMA data exceeds -1800 km s GLYPH&lt;0&gt; 1 ; Fig. 11) and the PdB sensitivity was not adequate to detect it: the PdB data were sensitive down to 1.1mJy beam. /

## 3.4.2. Extended wind: high-velocity CO in radially extending filaments

Herschel observations indicated the possible existence of an extended molecular wind in 4C12.50: unlike in most other ULIRGs, the outflowing OH molecules in 4C12.50 are primarily seen in emission (Spoon et al. 2013). This result indicated that many OH molecules are located in lines of sight free of background emission or that the background emission is low compared to the OH emission. A region in which an extended outflow is likely to be detected in 4C12.50 is north of the main nucleus (Rodríguez Zaurín et al. 2007). There, a very distinct [O iii ] ridge is seen in optical Hubble Space Telescope data, referred to as an "arc" by Tadhunter et al. (2018). Our ALMA data indicate that the molecular gas disk is rather abruptly cut along this ridge at low velocities (see CO line contours near systemic velocity; Fig. 7). This result could be indicative of a turnover of a tidal tail or stream or of an interface where a gas phase transition from molecular to ionized occurs. The transition could be due to the deposition of energy by some mechanism, for example, stellar or AGN radiation, collision-related shocks, or jet. Filaments that radially extend from the nucleus towards the ridge are additionally seen in the [O iii ] image. Two of them are also distinct in the stellar continuum image. Holt et al. (2003) found kinematically distorted [O ii ] emission along a 160 -position-angle slit GLYPH&lt;14&gt; that contained such filaments.

High-velocity clouds are detected tentatively along these filaments in the ALMA data, with signal-to-noise ratios between 3 and 5 (Fig. 12). In a filament along a position angle of -10 GLYPH&lt;14&gt; , designated as F1 in Fig. 12 and prominently seen in [O iii ], we see CO(1-0) at about -2200 km s GLYPH&lt;0&gt; 1 . The total H2 gas mass in region F1 is 0.5( GLYPH&lt;6&gt; 0.1) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; .

In a second filament along a position angle of -30 GLYPH&lt;14&gt; that is designated as F2 in Fig. 12 and that is prominently seen in the stellar continuum (Fig. 13), we again detect molecular gas. This time, we see CO(1-0) emission with velocities of about -

1200 km s GLYPH&lt;0&gt; 1 , -2200 kms GLYPH&lt;0&gt; 1 , and with a corresponding H2 mass equal to 1.3( GLYPH&lt;6&gt; 0.2) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; . The gas in or near this filament is counter-rotating in all CO lines. Moreover, this is the only extranuclear region with CO(4-3) detection of gas out of dynamical equilibrium; this filament is very nicely outlined in the third panel of the last row of Fig. 7 for gas counter-rotating with respect to the local disk velocity field. The molecular gas mass in F1 and F2 is 1.8( GLYPH&lt;6&gt; 0.2) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; .

In both filaments F1 and F2, we also observe marginal emission at positive velocities. Focusing on velocities near 1370 km s GLYPH&lt;0&gt; 1 we see that the redshifted emission has the shape of a ring, which is centered at the nucleus but occupies the whole disk. The H2 mass inside it, excluding the nuclear component contribution, is calculated using the integrated CO(1-0) flux of the area between the ellipses of the left panel of the third row of Fig. 12: namely 1.8( GLYPH&lt;6&gt; 0.4) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; . To ensure that all extranuclear components are accounted for, we extracted the spectrum from a larger area around the nucleus as shown in the last panels of Fig. 12. Subtracting the spectrum of the nuclear wind and integrating the CO(1-0) flux in the velocity range: -980 km s GLYPH&lt;0&gt; 1 &lt; V -&lt; 400 km s GLYPH&lt;0&gt; 1 we calculate the mass of the blueshifted extended emission to be: 6.0( GLYPH&lt;6&gt; 1.0) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; .

Adding the mass probed by the redshifted emission diskwide and the blueshifted extended emission to the mass probed by F1 and F2, we find that the total mass of the extended wind reaches 1.0( GLYPH&lt;6&gt; 0.1) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; and tentatively spreads over an area of roughly 40kpc 2 .

Several arguments are in favor of the credibility of the wind detections in the filaments and in the extended disk-wide region despite their low S N ratios. At / GLYPH&lt;24&gt; 0.7 00 blueshifted and redshifted emission partially overlaps in the disk-wide region. The detections are located along a line that connects the main nucleus to the region of highest-velocity dispersion, as shown in Fig. 14. This is observed in both CO(1-0) and CO(3-2) maps. The shapes and the loci of the high-velocity blobs F1 and F2 coincide with those of their optical filamentary counterparts (in either the [O iii ] or stellar continuum image; Fig. 13). The CO(1-0) emission in F1 and F2 has components in two or more velocity ranges. Additionally, F2 is detected in both CO(1-0) and CO(4-3). Some of the blobs also trace the excess emission that we see after the subtraction of a bulge and two disks from the H GLYPH&lt;11&gt; + continuum image (Fig. 13). The detection of stars along some regions also indicates two potential origin mechanisms for the accelerated gas: either it is related to a locally generated stellar wind or to an AGN generated nuclear wind. We evaluate the most likely scenario in Section 4, noting that the conclusions we draw about the origin of localized winds hold even if some detections turn out not to be real.

## 3.5. Gas excitation

Because the accelerated gas can have di GLYPH&lt;11&gt; erent excitation from that of the ambient gas, CO flux ratios are examined as wind indicators. Observed flux ratios, computed from the data at the common beam resolution, are shown in Fig. 15. Clear di GLYPH&lt;11&gt; erences are seen between the excitation of the nucleus and that of the southern tidal structure. The gas in the southern structure is subthermally excited. A range of excitations is seen between this region and the disk. Contrarily, higher excitation is observed near the nucleus. At positive velocities, some gas displays a CO(4-3) = CO(3-2) ratio above the value of 1.78, which is the upper limit in case of optically thick gas emission. Excitation temperatures from these maps will be presented in a forthcoming paper (Paraschos et al. 2019, in prep.).

Potential excitation di GLYPH&lt;11&gt; erences can also be revealed by the comparison of multi-wavelength, multi-temperature gas probes. From previous Spitzer data, we know that in the warm ( GLYPH&lt;24&gt; 400 K ) H2 phase, the outflow has a mass of 5.2 GLYPH&lt;2&gt; 10 7 M GLYPH&lt;12&gt; (Dasyra et al. 2011), whereas the disk has a mass of 1.4 GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; (Dasyra et al. 2011). Therefore, the fraction of warm H2 in the wind exceeds 30%. From our new ALMA data, we find that the H2 mass of the wind is as high as 1.5 GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; whereas that of the ambient gas is 7.7 GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; . Therefore, the fraction of cold H2 in the wind is less than 20%. This confirms our past findings of heating of the accelerated gas. It is noteworthy that some of the cold gas mass in the wind may be unaccounted for due to the spatial overlap of multiple structures out of dynamical equilibrium in our ALMA data. Still, the mass loss cannot be as high as required to alter our initial conclusion. The e GLYPH&lt;11&gt; ects of gas acceleration and heating could possibly delay star formation temporarily.

Table 1. Properties of tentative extranuclear wind detection

| Region               | M ( GLYPH<2> 10 8 M GLYPH<12> )   |   d ( 00 ) |   j V mean j ( kms GLYPH<0> 1 ) |   ˙ M ( M GLYPH<12> yr GLYPH<0> 1 ) |   L kin ( GLYPH<2> 10 44 ergs GLYPH<0> 1 ) |   ˙ MV ( GLYPH<2> 10 36 ergcm GLYPH<0> 1 ) |
|----------------------|-----------------------------------|------------|---------------------------------|-------------------------------------|--------------------------------------------|--------------------------------------------|
| F1                   | 0.5( GLYPH<6> 0.1)                |        0.6 |                            2200 |                                  90 |                                        1.3 |                                        1.2 |
| F2                   | 0.6( GLYPH<6> 0.2)                |        0.8 |                            2180 |                                  70 |                                        1.1 |                                        1   |
| F2                   | 0.7( GLYPH<6> 0.2)                |        0.8 |                            1120 |                                  50 |                                        0.2 |                                        0.3 |
| redshifted disk-wide | 1.8( GLYPH<6> 0.4)                |        0.8 |                            1370 |                                 140 |                                        0.9 |                                        1.2 |
| blueshifted extended | 6.0( GLYPH<6> 1.0)                |        1.2 |                             700 |                                 160 |                                        0.3 |                                        0.7 |

Notes. Using 0.1 00 as a typical error for d and 10 km s GLYPH&lt;0&gt; 1 as a typical error for Vmean , the errors in ˙ , M Lkin , and ˙ MV are roughly 20% of their values.

## 4. Discussion

## 4.1. Energy output of the galaxy at radio and infrared wavelengths

To evaluate whether the jet, the AGN radiation pressure, or the starburst can sustain winds, we first need to calculate the energy output of the galaxy at radio and infrared wavelengths. To obtain the jet power, Pj , we fitted the SED at radio, IR, optical, UV , and X-ray wavelengths. Because the jet and the dust can emit at common frequencies, we simultaneously fitted their SEDs using a standard GLYPH&lt;31&gt; 2 minimization method (Newville et al. 2014). For the dust emission, we used a modified black-body law at long wavelengths, coupled with a power law at short wavelengths (Casey 2012). We found the dust temperature to be 55 K and L IR to be 2.5 GLYPH&lt;2&gt; 10 12 L GLYPH&lt;12&gt; (in good agreement with Dasyra et al. 2014). For Pj , we fitted the synchrotron radiation using a broken power law with an exponential cut o GLYPH&lt;11&gt; (Rybicki &amp; Lightman 1979). At radio wavelengths, up to 4 GLYPH&lt;2&gt; 10 11 Hz, this model is well determined from our data (Fig. 16). The integral of the fit provides the lower limit for the synchrotron power, which is 9.1 GLYPH&lt;2&gt; 10 43 erg = s. At higher frequencies however, the dust, stellar, and AGN emission outshine the jet emission. To find an upper limit for the synchrotron power, we fitted the data under the extreme assumption that 0.6% of the 6 GLYPH&lt;2&gt; 10 14 Hz emission originates from the jet. The fraction 0.6% corresponds to the part of the 5092 Å flux that is enclosed in a radius of 150 pc in the Hubble Space Telescope data; 150 pc is the distance at which the radio jet interacts with ISM clouds per Morganti et al. (2013a). In this model (Fig. 16), Pj is 6.2 GLYPH&lt;2&gt; 10 44 erg = s. The X-rays may also be bolstered by a contribution from various Compton components related to the jet or

magnetic fields (Finke 2016). To obtain the maximum possible contribution of the jet to the X-ray emission, we again adopted a power law with a cut-o GLYPH&lt;11&gt; that was steeper than before due to the absence of pertinent data above 10 19 Hz. The maximum power output in the X-rays is 8.4 GLYPH&lt;2&gt; 10 43 erg = s. In summary, we find the jet power to be in the range 10 44 GLYPH&lt;0&gt; 8 GLYPH&lt;2&gt; 10 44 erg = s. For comparison, Guillard et al. (2012) calculated the power of the radio-jet of 4C12.50 from the monochromatic 178 MHz flux, following Punsly (2005). They adopted a calibration of the monochromatic flux to the bolometric output that takes into account X-rays and the plasma thermal energy in the lobes. Their approach yielded Pj = 3 GLYPH&lt;2&gt; 10 45 erg s GLYPH&lt;0&gt; 1 . Likewise, using the calibration between the monochromatic 1.4 GHz luminosity and the bolometric jet luminosity following Sulentic et al. (2010), Pj can be as high as 10 46 erg s GLYPH&lt;0&gt; 1 .

To obtain the force exerted on the gas due to AGN radiation pressure we assume that the AGN luminosity is well described by its L IR value, as the flux that was absorbed and reemitted by the dust is higher than that seen in the optical. Following Veilleux et al. (2009), we ascribe half the L IR to the AGN. Based on our previous analysis, the AGN-related part of the L IR is 4.7 GLYPH&lt;2&gt; 10 45 erg s GLYPH&lt;0&gt; 1 . The force exerted on the gas due to the AGN radiation pressure is LAGN = c = 1.6 GLYPH&lt;2&gt; 10 35 erg cm GLYPH&lt;0&gt; 1 .

To evaluate whether or not the starburst can sustain winds, we need to calculate the energy released by SNe and the force exerted on the gas due to stellar radiation pressure. We assume that for every 100 M GLYPH&lt;12&gt; formed, there is one SN. We also assume that each SN ejects material with a kinetic energy of 10 51 erg . Given a star formation rate (SFR) of GLYPH&lt;24&gt; 200 M GLYPH&lt;12&gt; yr GLYPH&lt;0&gt; 1 , computed from the Kennicutt (1998) formula for L IR / 2, the power released by the SN is 6 GLYPH&lt;2&gt; 10 43 ergs GLYPH&lt;0&gt; 1 . Also taking into account the fraction of the total star formation that can be ascribed to any given area (as explained in Sections 4.2 and 4.3) we calculate the power released locally by the SN. Furthermore, using the fraction of the total L IR = 2 that can be attributed to a local starburst, we calculate the force exerted on the gas due to stellar radiation pressure.

## 4.2. Extended wind energetics: gas acceleration not sustained by a local starburst

We begin our energetics study of the extended wind by ruling out that the low-S N CO detections in filaments are accelerated / by the local starburst. For this purpose, we attribute a part of the SFR or of L IR to each ellipse of F1 and F2 of Fig. 12 using the ratio of the locally enclosed CO(1-0) emission over the total CO(1-0) emission in our high-resolution data. The filament F1, for example, comprises roughly 8% of the total CO emission. Based on this fraction, the local L IR is 3.8 GLYPH&lt;2&gt; 10 44 erg s GLYPH&lt;0&gt; 1 and the local SFR is GLYPH&lt;24&gt; 15 M GLYPH&lt;12&gt; yr GLYPH&lt;0&gt; 1 . In this case, the mechanical energy from the SN is 4.8 GLYPH&lt;2&gt; 10 42 erg s GLYPH&lt;0&gt; 1 and the force due to stellar radiation pressure is 1.3 GLYPH&lt;2&gt; 10 34 erg cm GLYPH&lt;0&gt; 1 . To calculate the energetics of the gas in F1, we need to use the distance d of the accelerated gas from the spot of the wind generation. For this purpose, we use the mean value of the ellipse semi-axes (0.4 00 ) as the distance from the star-forming region. Using the mass and the mean velocity presented in Section 3.4.2, we find that the wind kinetic luminosity, Lkin , computed as (1 = 2) MV d 3 GLYPH&lt;0&gt; 1 , is 2.0 GLYPH&lt;2&gt; 10 44 erg s GLYPH&lt;0&gt; 1 . The wind momentum rate, ˙ MV (Combes et al. 2013), is 1.8 GLYPH&lt;2&gt; 10 36 erg cm GLYPH&lt;0&gt; 1 . In this formula, ˙ M is the mass-flow rate of the accelerated gas, equal to MVd GLYPH&lt;0&gt; 1 . We find that neither the SN nor the stellar radiation pressure in the area can drive the wind along the filament, each being at least two orders of magnitude short of the required levels.

The same calculation for the filament F2 shows that it comprises roughly 9% of the total CO emission. Therefore, the local L IR is 4.2 GLYPH&lt;2&gt; 10 44 erg s GLYPH&lt;0&gt; 1 and the local SFR is GLYPH&lt;24&gt; 20 M GLYPH&lt;12&gt; yr GLYPH&lt;0&gt; 1 . The mechanical energy released by the SN is 5.4 GLYPH&lt;2&gt; 10 42 erg s GLYPH&lt;0&gt; 1 and the force due to stellar radiation pressure is 1.4 GLYPH&lt;2&gt; 10 34 erg cm GLYPH&lt;0&gt; 1 . The distance d is 0.7 00 and calculating the wind energetics as in the case of F1 we find that Lkin is 1.4 GLYPH&lt;2&gt; 10 44 erg s GLYPH&lt;0&gt; 1 and ˙ MV is 1.5 GLYPH&lt;2&gt; 10 36 erg cm GLYPH&lt;0&gt; 1 . In this case neither the SN nor the stellar radiation pressure in the area are adequate drivers of the wind.

Inversely, how high does the local SFR need to be to sustain the wind? If the SN were a local generation mechanism, the SFR would need to be higher than 630 M GLYPH&lt;12&gt; yr GLYPH&lt;0&gt; 1 and 440 M GLYPH&lt;12&gt; yr GLYPH&lt;0&gt; 1 in the areas of F1 and F2, respectively, in order to sustain the wind. Therefore, the local SNe cannot drive the flow by themselves in the areas of the filaments. However, given that the velocity in the filaments is the highest that we have observed, the local starburst is likely to assist in their acceleration. Still, the main driver of all wind candidates needs to be sought in a central mechanism.

## 4.3. Overall wind energetics: sustainability by a central mechanism

To evaluate whether or not a power source at the nucleus (i.e., within our beam) could sustain the wind, we evaluate whether the wind can be sustained via the nuclear starburst or via the jet and the AGN. For the nuclear starburst, we compare the energy deposited by SNe to the wind kinetic luminosity and, as above, the force exerted on the gas due to stellar radiation pressure to the wind momentum rate. As explained in Section 4.1 the power released by the SN is 6 GLYPH&lt;2&gt; 10 43 ergs GLYPH&lt;0&gt; 1 . Taking further into account that the nuclear SN can only be ascribed 10-20% of the total star formation, then the total SN power is 6-12 GLYPH&lt;2&gt; 10 42 ergs GLYPH&lt;0&gt; 1 . The fraction ascribed to nuclear star formation is computed from the fraction of CO(1-0) emission within the radius of a beam (0.3 00 ). Furthermore, taking into account that 10-20% of L IR / 2 can be ascribed to the nuclear starburst, that is, 4.7-9.4 GLYPH&lt;2&gt; 10 44 erg s GLYPH&lt;0&gt; 1 , the force exerted on the gas due to stellar radiation pressure, LSB = c , is then 1.6-3 GLYPH&lt;2&gt; 10 34 erg cm GLYPH&lt;0&gt; 1 . Alternatively, the force exerted on the gas due to the AGN radiation pressure is LAGN = c = 1.6 GLYPH&lt;2&gt; 10 35 erg cm GLYPH&lt;0&gt; 1 .

To find the central mechanism that could sustain the extended wind, we need to compute the respective wind kinetic luminosity and momentum rate assuming the distance of the accelerated clouds from the main nucleus. For this purpose, we compute Lkin and ˙ MV for each individual region of Fig. 12 using its distance from the nucleus. We sum the results for the extended regions with the F1, the F2 (using the two velocity ranges of F2 as two di GLYPH&lt;11&gt; erent regions), the redshifted disk-wide and the extended blueshifted detections from Table 1. We find that the mass-flow rate ˙ M carried by the extended wind is as high as GLYPH&lt;24&gt; 500 M GLYPH&lt;12&gt; yr GLYPH&lt;0&gt; 1 . Then, Lkin (extended) is 3.8 GLYPH&lt;2&gt; 10 44 erg s GLYPH&lt;0&gt; 1 and ˙ MV (extended) is 4.4 GLYPH&lt;2&gt; 10 36 erg cm GLYPH&lt;0&gt; 1 . Therefore, the mechanical power of the SN is insu GLYPH&lt;14&gt; cient to drive the flow, considering that it is low and that only a part of it is deposited in the ISM. Likewise, the stellar radiation pressure within the central 0.3 00 cannot drive the wind, since the momentum rate of the wind is two orders of magnitude higher than the force exerted on gas by the stellar radiation pressure. As this momentum rate is also an order of magnitude higher than LAGN = c , the AGN radiation pressure is unable to drive this wind, unless an energy-conserving expansion has significantly boosted the momentum rate of the gas (i.e., by a factor of 30). Such a boost is high, but not impossible, as shown by Cicone et al. (2014).

To add to the above numbers those for the circumnuclear wind, we set d equal to 0.3 00 , which is approximately half of the semi-major axis of the ellipse in Fig. 11. The mass of the circumnuclear wind is 5.3( GLYPH&lt;6&gt; 0.7) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; , and the mean velocity of either the blueshifted or the redshifted component is GLYPH&lt;24&gt; 800 km s GLYPH&lt;0&gt; 1 . Using these numbers and assuming that the H2 clouds and the COclouds have similar spatial and velocity distributions, we find that the kinetic luminosity of the nuclear wind is 1.3 GLYPH&lt;2&gt; 10 44 ergs GLYPH&lt;0&gt; 1 and the momentum rate is 3.3 GLYPH&lt;2&gt; 10 36 erg cm GLYPH&lt;0&gt; 1 . Adding these numbers to the values we calculated for the extended wind, the total kinetic luminosity and the total momentum rate of the outflow reach 5.1 GLYPH&lt;2&gt; 10 44 ergs GLYPH&lt;0&gt; 1 and 7.7 GLYPH&lt;2&gt; 10 36 erg cm GLYPH&lt;0&gt; 1 , respectively. These numbers may yet increase, considering that part of the circumnuclear wind is unresolved and thus the used radius may be overestimated.

Overall, a combination of the jet, radiation pressure, and SN mechanical power needs to be invoked. The jet is the most likely generation mechanism given that it carries the most power. The jet alone could have driven the flow if a past event, which has been recorded for this system, were found to have deposited enough energy onto the ISM. In this case we are observing a fossil outflow which expands into an already carved cavity (Fluetsch et al. 2019). This could explain the radial symmetry of the wind even though the radiation pressure is equally wellsuited to explain it. However, the force due to the radiation pressure cannot drive the outflow by itself, but could assist it. The same applies to the local generation mechanism along the filaments: a local starburst could further accelerate the gas that has already been accelerated.

## 5. Conclusions

We obtained ALMA data of the radio galaxy 4C12.50 in order to determine the millimeter properties of its outflow, which was previously known to exist from large- and small-angularresolution observations at other wavelengths. We mapped the COdistribution, kinematics, and excitation with ALMA at a resolution of GLYPH&lt;24&gt; 0.5 00 and found the following results.

- -The main gaseous disk is rather compact, extending to radii of GLYPH&lt;24&gt; 6 kpc in CO (1 GLYPH&lt;0&gt; 0) and GLYPH&lt;24&gt; 2 kpc in CO (4 GLYPH&lt;0&gt; 3). Some corotating CO emission is also seen in the secondary nucleus and in a bridge connecting the nuclei.
- -A shock-tracer, H2O, and a hard-ionization-field tracer, HCO (5-4) are seen in the main nucleus. +
- -Several extranuclear structures are seen. These include a prominent tidal tail south of the main nucleus with gas seen in negative and positive velocities.
- -The CO emission is abruptly cut along an [O iii ] ridge north of the main nucleus, where no CO is detected. This indicates a phase transition of the molecular gas.
- -Extended wind components are tentatively seen in the CO(1-0) line. The detections are also seen in regions with [O iii ] and stellar continuum emission. They include highvelocity (-2000 kms GLYPH&lt;0&gt; 1 ) filaments that extend radially from the nucleus to the [O iii ] ridge, a redshifted disk-wide emission (V GLYPH&lt;24&gt; 1350 kms GLYPH&lt;0&gt; 1 ) and an extended blueshifted emission (V GLYPH&lt;24&gt; -700 kms GLYPH&lt;0&gt; 1 ). The total mass is 1.0( GLYPH&lt;6&gt; 0.1) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; .
- -A circumnuclear wind candidate is detected in emission, in CO(1-0), through broad blueshifted and redshifted line wings. The emission peaks within 200 pc from the radio core, and it coincides with a region of high CO excitation (within the velocity range of the disk), as indicated by its CO(4-3) / CO(3-2) line ratio. The mass of the circumnuclear wind is equal to 5.3( GLYPH&lt;6&gt; 0.7) GLYPH&lt;2&gt; 10 8 M GLYPH&lt;12&gt; .
- -The total mass of the wind is as high as 1.5( GLYPH&lt;6&gt; 0.1) GLYPH&lt;2&gt; 10 9 M GLYPH&lt;12&gt; .
- -Both at the nucleus and at the extra-nuclear regions, the wind can be sustained mainly by the jet. The radiation pressure of the AGN can help, in particular near the nucleus. It is plausible that the starburst also contributes, in particular for any extended wind components. However, the starburst cannot drive the wind alone based on its radiation pressure or its energy release by SN remnants.

Acknowledgements. This paper makes use of the ALMA data ADS JAO.ALMA / 2013.1.00180.S, 2013.1.00976.S, 2012.1.00797.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan) and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI NRAO and NAOJ. K. M. Dasyra acknowledges financial / support by the Hellenic Foundation for Research and Innovation (HFRI) and the General Secretariat For Research and Technology, under grant number 1882. We would like to thank the referee, B. Emonts, for detailed comments which led to a significant improvement of the paper.

## References

Aalto, S., Costagliola, F., Muller, S., et al. 2016, A&amp;A, 590, A73 Aalto, S., Garcia-Burillo, S., Muller, S., et al. 2012, A&amp;A, 537, A44 Alatalo, K., Blitz, L., Young, L. M., et al. 2011, ApJ, 735, 88 Alatalo, K., Nyland, K., Graves, G., et al. 2014, ApJ, 780, 186 Batcheldor, D., Tadhunter, C., Holt, J., et al. 2007, ApJ, 661, 70 Carniani, S., Marconi, A., Maiolino, R., et al. 2015, A&amp;A, 580, A102 Casey, C. M. 2012, MNRAS, 425, 3094 Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&amp;A, 562, A21 Combes, F., García-Burillo, S., Casasola, V., et al. 2013, A&amp;A, 558, A124 Dasyra, K. M. &amp; Combes, F. 2011, A&amp;A, 533, L10 Dasyra, K. M. &amp; Combes, F. 2012, A&amp;A, 541, L7 Dasyra, K. M., Combes, F., Novak, G. S., et al. 2014, A&amp;A, 565, A46 Dasyra, K. M., Combes, F., Oosterloo, T., et al. 2016, A&amp;A, 595, L7 Dasyra, K. M., Ho, L. C., Netzer, H., et al. 2011, ApJ, 740, 94 Downes, D. &amp; Solomon, P. M. 1998, ApJ, 507, 615 Emonts, B. H. C., Morganti, R., Villar-Martín, M., et al. 2016, A&amp;A, 596, A19 Evans, A. S., Kim, D. C., Mazzarella, J. M., Scoville, N. Z., &amp; Sanders, D. B.

- 1999, ApJ, 521, L107

Faucher-Giguère, C.-A. &amp; Quataert, E. 2012, MNRAS, 425, 605

Feruglio, C., Maiolino, R., Piconcelli, E., et al. 2010, A&amp;A, 518, L155

Finke, J. D. 2016, ApJ, 830, 94

- Fischer, J., Sturm, E., González-Alfonso, E., et al. 2010, A&amp;A, 518, L41

Fluetsch, A., Maiolino, R., Carniani, S., et al. 2019, MNRAS, 483, 4586

García-Burillo, S., Combes, F., Usero, A., et al. 2014, A&amp;A, 567, A125

George, R. D., Ivison, R. J., Smail, I., et al. 2014, MNRAS, 442, 1877

González-Alfonso, E., Fischer, J., Spoon, H. W. W., et al. 2017, ApJ, 836, 11

- Guillard, P., Ogle, P. M., Emonts, B. H. C., et al. 2012, ApJ, 747, 95
- Holt, J., Tadhunter, C. N., &amp; Morganti, R. 2003, MNRAS, 342, 227
- Holt, J., Tadhunter, C. N., &amp; Morganti, R. 2008, MNRAS, 387, 639
- Holt, J., Tadhunter, C. N., &amp; Morganti, R. 2009, MNRAS, 400, 589
- Holt, J., Tadhunter, C. N., Morganti, R., &amp; Emonts, B. H. C. 2011, MNRAS, 410, 1527
- Imanishi, M., Nakanishi, K., &amp; Izumi, T. 2016, AJ, 152, 218
- Imanishi, M., Nakanishi, K., &amp; Izumi, T. 2018, ApJ, 856, 143
- Imanishi, M. &amp; Saito, Y. 2014, ApJ, 780, 106
- Ishibashi, W. &amp; Fabian, A. C. 2015, MNRAS, 451, 93
- Kennicutt, Robert C., J. 1998, ApJ, 498, 541
- King, A. R., Zubovas, K., &amp; Power, C. 2011, MNRAS, 415, L6
- Krips, M., Martín, S., Eckart, A., et al. 2011, ApJ, 736, 37
- Leon, S., Eckart, A., Laine, S., et al. 2007, A&amp;A, 473, 747
- Lister, M. L., Kellermann, K. I., Vermeulen, R. C., et al. 2003, ApJ, 584, 135
- Morganti, R., Fogasy, J., Paragi, Z., Oosterloo, T., &amp; Orienti, M. 2013a, Science, 341, 1082
- Morganti, R., Frieswijk, W., Oonk, R. J. B., Oosterloo, T., &amp; Tadhunter, C. 2013b, A&amp;A, 552, L4
- Morganti, R., Oosterloo, T. A., Tadhunter, C. N., et al. 2004, A&amp;A, 424, 119

Newville, M., Stensitzki, T., Allen, D. B., &amp; Ingargiola, A. 2014, LMFIT: Non-

Linear Least-Square Minimization and Curve-Fitting for Python

- Oosterloo, T., Raymond Oonk, J. B., Morganti, R., et al. 2017, A&amp;A, 608, A38 Peng, C. Y., Ho, L. C., Impey, C. D., &amp; Rix, H.-W. 2002, AJ, 124, 266
- Pereira-Santaella, M., Colina, L., García-Burillo, S., et al. 2018, A&amp;A, 616, A171
- Punsly, B. 2005, ApJ, 623, L9

Rangwala, N., Maloney, P. R., Glenn, J., et al. 2011, ApJ, 743, 94

Article number, page 7 of 19

Rodríguez Zaurín, J., Holt, J., Tadhunter, C. N., &amp; González Delgado, R. M. 2007, MNRAS, 375, 1133

Rybicki, G. B. &amp; Lightman, A. P. 1979, Radiative processes in astrophysics

Sakamoto, K., Aalto, S., Combes, F., Evans, A., &amp; Peck, A. 2014, ApJ, 797, 90

Sakamoto, K., Aalto, S., Wilner, D. J., et al. 2009, ApJ, 700, L104

Sakamoto, K., Ho, P. T. P., &amp; Peck, A. B. 2006, ApJ, 644, 862

Siemiginowska, A., LaMassa, S., Aldcroft, T. L., Bechtold, J., &amp; Elvis, M. 2008, ApJ, 684, 811

Solomon, P. M., Downes, D., Radford, S. J. E., &amp; Barrett, J. W. 1997, ApJ, 478, 144

Spoon, H. W. W., Farrah, D., Lebouteiller, V., et al. 2013, ApJ, 775, 127 Stanghellini, C., O'Dea, C. P., Baum, S. A., et al. 1997, A&amp;A, 325, 943 Stanghellini, C., O'Dea, C. P., Baum, S. A., &amp; Laurikainen, E. 1993, ApJS, 88, 1

Stanghellini, C., O'Dea, C. P., Dallacasa, D., et al. 2005, A&amp;A, 443, 891 Stone, M., Veilleux, S., Meléndez, M., et al. 2016, ApJ, 826, 111 Sturm, E., González-Alfonso, E., Veilleux, S., et al. 2011, ApJ, 733, L16 Sulentic, J., Zamfir, S., &amp; Marziani, P. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 427, Accretion and Ejection in AGN: a Global

View, ed. L. Maraschi, G. Ghisellini, R. Della Ceca, &amp; F. Tavecchio, 371 Tadhunter, C., Rodríguez Zaurín, J., Rose, M., et al. 2018, MNRAS, 478, 1558 Tombesi, F., Meléndez, M., Veilleux, S., et al. 2015, Nature, 519, 436 Toomre, A. &amp; Toomre, J. 1972, ApJ, 178, 623

- Tsai, A.-L., Matsushita, S., Kong, A. K. H., Matsumoto, H., &amp; Kohno, K. 2012, ApJ, 752, 38

Veilleux, S., Meléndez, M., Sturm, E., et al. 2013, ApJ, 776, 27

Veilleux, S., Rupke, D. S. N., Kim, D. C., et al. 2009, ApJS, 182, 628

Wagner, A. Y., Bicknell, G. V., Umemura, M., Sutherland , R. S., &amp; Silk, J. 2016,

Astronomische Nachrichten, 337, 167

Wright, E. L. 2006, PASP, 118, 1711

Zubovas, K. &amp; King, A. R. 2014, MNRAS, 439, 400

Fig. 1. Optical images of 4C12.50 from the Hubble Space Telescope , previously presented by Batcheldor et al. (2007) and Dasyra et al. (2011). The crosses mark the location of the peak continuum emission in the ALMA band 3 data. The same applies to all following figures with the crosses. Upper left : Gas emission comprising primarily [O iii ] emission with potential contribution from H GLYPH&lt;12&gt; . The image was created from the subtraction of a Hubble Space Telescope ACS image at 5092 Å (filter FR459M) from an ACS image at 5580 Å (filter F550M). Upper right : Continuum emission at 5092 Å (filter FR459M). Lower left : H GLYPH&lt;11&gt; and continuum emission (filters FR647M at 6616 Å). Lower right : Residual image originating from the subtraction of galactic component models (i.e., a bulge and two disks) from an optical Hubble Space Telescope image of H GLYPH&lt;11&gt; and 5900 Å continuum emission (Batcheldor et al. 2007) using GALFIT (Peng et al. 2002). The subtraction method was described in detail in Dasyra et al. (2011). Its purpose was to bring up any structures underlying the bright regular emission.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-40a72ce4-45df-48eb-9d60-f2cd0b3723d0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-f128f082-fbf0-4a6e-9781-746b36694ee3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-9a1d154d-15af-4da3-94ee-1f008fd9f175.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-3ef20b6a-538e-4b69-b415-3c503d7441ef.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-5638b32b-4a30-460b-985c-1e3f602e174c.jpg)

Fig. 2. Left: X-ray emission (Siemiginowska et al. 2008) plotted over an IRAM Plateau de Bure CO(1-0) image (Dasyra et al. 2014). Middle, right: Large and small scale radio emission at GLYPH&lt;24&gt; 1.4 and 5 GHz (Stanghellini et al. 2005) compared to the same CO(1-0) image.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-a0e33646-c6e3-42b7-bbd7-7413a8a759ec.jpg)

Fig. 3. CO(1-0) , CO(3-2) and CO(4-3) galaxy spectra derived from our high-resolution data cube (upper panels) and the deep data (lower panels) that originated from the merging of the measurements of two programs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-949cb9b8-c33f-42da-a995-8712afb1e74a.jpg)

Fig. 4. Spectra of all molecular gas lines and continuum emission detected in the ALMA data of the program 2013.1.00180.S. Left : Spectra of the entire galaxy. Right : Spectra of the nuclear pixel. The dashed line in CO(3-2) spectrum corresponds to the frequency at which the stitching of the data of the two spectral windows, in which the continuum level and slope are not known with su GLYPH&lt;14&gt; cient accuracy for the determination of the presence of a wind, took place.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-d50150a2-9ac2-4b60-a9bb-037cb4a10677.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-15c40131-f54b-4fb1-b9d0-5bd8e31dc3a2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-63859ed0-3da3-47fe-bb6b-81d0bf37fa42.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-49663b3d-1f2e-4a18-b2fa-ab9b07404d09.jpg)

3

Fig. 5. H2O(4 GLYPH&lt;0&gt; 3) in 4C12.50. The left panel shows the ALMA image, collapsed in the indicated velocity range. The beam is shown with a shaded ellipse. The open ellipse corresponds to the area used for the extraction of the spectrum in the middle panel. In the right panel we plot the CO contours (levels: 2 GLYPH&lt;27&gt; ,3 GLYPH&lt;27&gt; and 5 GLYPH&lt;27&gt; ) over the Hubble Space Telescope [O iii ] image of Fig.1. Crosses are as in Fig 1. The same applies to all following figures with similar layout.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-49642012-ef0f-450c-bb26-24a14b9fcfa7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-304498b5-f6ef-4add-9c24-e772096c974b.jpg)

3

Fig. 6. HCO (5-4) in 4C12.50. Contour levels are at 2 + GLYPH&lt;27&gt; ,3 GLYPH&lt;27&gt; ,5 GLYPH&lt;27&gt; and 10 GLYPH&lt;27&gt; .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-c729649e-bdf2-4081-9db1-04ce08f0d409.jpg)

Jy

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-419fbe55-a385-476f-8e31-10fcd1030ae0.jpg)

J

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-915657c8-0de5-4f30-a454-71ff95b4deac.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-d2c9d237-7fb8-47f8-9dbe-94fb1b5fce09.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-09e69c22-9e42-4f8d-84c6-09db1cb61d25.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-9708261f-ab5f-4014-a85e-38a1e9d821b3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-7f82bc31-45b7-4ffe-aa50-259f70293a1b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-12215b4e-cbf1-4cae-be07-1ed26651f598.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-9cec77d5-4784-400a-98c0-f3ef50d606f3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-b1eb7317-715c-47a6-ae3e-692b24c4176a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-f7ecb98e-febd-4566-bd46-51134c649045.jpg)

3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-96cc6293-b1aa-4886-b0af-83a1a2a9bada.jpg)

3

Fig. 7. Upper three panels: CO(1-0), CO(3-2) and CO(4-3) emission in di GLYPH&lt;11&gt; erent velocity ranges. Filled ellipses show the beam size at each frequency. The same applies to all of the following figures. Lower panel: CO contours with levels 2 GLYPH&lt;27&gt; and 3 GLYPH&lt;27&gt; over the [O iii ] image in the optical showing the extent of the molecular gas near systemic velocity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-f24d9bde-7890-433a-b32e-6a7c8d46658a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-dc2f1489-576d-450f-a8b4-763d203325e5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-66b9f3ca-0b15-426e-abd3-085adb07a140.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-9bb9e704-4918-4b1c-a260-10ae0e5116ae.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-2282fcf3-0fac-436d-8db4-264300311952.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-40383360-eeb2-4e4c-9711-6bb5733a80e0.jpg)

Jy

Fig. 8. CO(1-0), CO(3-2) and CO(4-3) momenta maps for pixels with flux exceeding three times the noise root mean square. Left: Flux. Middle: Velocity. Right: Velocity dispersion.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-b56b0cb3-1742-4d9e-9dbb-ad958b8116bb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-08312e13-10ce-4424-9133-a244e90e4894.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-a152baa4-0653-4085-86ae-27fb2f23dbb8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-eda3f9ec-2649-418a-964e-19f24bfce95f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-a64cfc4c-aab5-454c-b069-fdc283b8b1b0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-9e112d74-d8c7-492a-ac56-bc7d112a7722.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-a445e2cb-c261-4406-b0bb-49c678e5e6f8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-f1c98779-04e0-4dc3-b779-1d0846226e95.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-d3fa7040-f3fc-4e26-aad4-a356fd325848.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-de8bf92f-c0aa-42c7-bda4-e08496639f6b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-4ee26330-0528-4270-a9d0-4202110680e0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-7044fbca-c6e4-46e2-ac80-3af4d4b94304.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-5433241f-0d88-417b-9dc4-c1a3096b6fa9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-bb2b4664-521a-4a51-a267-71e7d84c7c1a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-d3b0dd54-2e70-461c-b8ae-44564f155f63.jpg)

3

Fig. 9. Redshifted part of the emission south of the main nucleus. The contour levels in the right panels are: 2 GLYPH&lt;27&gt; ,3 GLYPH&lt;27&gt; ,5 GLYPH&lt;27&gt; for the first three rows and 2 GLYPH&lt;27&gt; ,3 GLYPH&lt;27&gt; ,4 GLYPH&lt;27&gt; for the last row.

Fig. 10. Blueshifted part of the emission south of the main nucleus. The contour levels in the right panels are: 2 GLYPH&lt;27&gt; ,3 GLYPH&lt;27&gt; ,5 GLYPH&lt;27&gt; ,10 GLYPH&lt;27&gt; ,11 GLYPH&lt;27&gt; ,13 GLYPH&lt;27&gt; for CO(1-0) and 2 GLYPH&lt;27&gt; ,3 GLYPH&lt;27&gt; ,5 GLYPH&lt;27&gt; for CO(3-2) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-294a5eed-fb94-48e8-a7e4-21b755f708d8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-0e2719b9-c53a-4b42-a90a-824b58db57f0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-cb3e74dc-d273-4b0b-84b2-50d1d0776203.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-ec558ad4-a982-4e7a-89f9-068c7f4fec3e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-c9461a7a-bd66-473f-ae80-f0bcfd26fdfb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-b787bc72-b7dc-4484-a333-5a6baacbbad8.jpg)

Fig. 11. Nuclear wind candidate. Average image of the blueshifted and the redshifted emission of the nuclear wind (left). A Gaussian fit to the nuclear spectrum is performed and fit residuals are presented in the middle and right panels, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-bca7ee6a-d598-4e7a-a385-474f314c634c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-812f1cc1-fdb1-4308-a4d3-4859dff21bf5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-773787ce-be22-4c9c-9d93-d8913e268053.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-da24ecc2-6cd1-49f0-b933-aa912029b689.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-69624056-f38c-4ddd-bff8-cc80eed9d40a.jpg)

Jy

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-042def5c-bd8d-4bab-833d-b1ea070cbb4c.jpg)

Jy

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-53006b79-293f-404f-b16c-cb6dbde3b291.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-c1e92429-1d3f-4e81-be17-ac6bb6c7f7b6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-8bb7c7a2-04b0-46d2-a4c9-cfdf4e038a23.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-1f014f02-9c72-45dc-b9d8-2d5ff6f9b893.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-1a3041e3-eef8-4e72-a68a-351eeea26b44.jpg)

3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-6e524b10-beb9-4fc4-a183-7fddf9e7b059.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-4e1183ce-9722-4808-8f35-8907bde359da.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-c9f5e3d0-52c3-427b-92c4-83b88db3f993.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-7e03e31f-eb29-4f9c-b773-2b35974559ce.jpg)

3

Fig. 12. Tentative CO detections of high-velocity clouds in directions radially extending from the main nucleus to outer parts of the galaxy. The contour levels in the right panels are: 2 GLYPH&lt;27&gt; ,3 GLYPH&lt;27&gt; and 4 GLYPH&lt;27&gt; . In the last two panels, the solid-line and the dotted spectra (middle) correspond to the emission of the solid-line and the dotted ellipses (left), respectively. The latter marks the extent of the nuclear wind which is shown in Fig. 11.

1′′  = 2.2 kpc Fig. 13. Loci of tentative detections of high-velocity molecular clouds compared to the stellar and warm, ionized gas emission. The contours are as in Fig. 12 for regions that exceed a S N of 3. Contours of comparable velocity ranges are plotted with the same colors. /

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-93e02615-815a-435d-b3ba-17674161c5e2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-c1d81b58-1a0a-4cb2-b9b6-1439760f054c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-ce6f7747-c7b4-42bb-83a4-97eefa5134a2.jpg)

## A &amp; A proofs: manuscript no. AA\_34416

Fig. 14. Loci of tentative detections of high-velocity molecular clouds compared to the velocity and the velocity dispersion maps of CO(1-0) and CO(3-2) of Fig. 8. The contours are as in Fig. 13.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-18282901-af30-476b-bc2f-3822e11c1ca1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-99870b3b-6207-4a42-be67-efaeb8da789a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-f6c9f0c5-e2c7-451e-ace1-6fa3295c169e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-61e546ea-3898-4d39-b4df-5b9767bf9395.jpg)

Fig. 15. CO(3-2) CO(1-0), = CO(4-3) CO(1-0) and CO(4-3) CO(3-2) flux = = ratios. Values exceeding 9, 16 and 1.78 for CO(3-2) CO(1-0) , = CO(4-3) CO(1-0) and CO(4-3) CO(3-2) respectively, indicate the presence of optically thin gas. First row: Flux ratios in the velocity range = = of the entire disk. Second row: Flux ratios in disk velocities that show high excitation close to the nucleus. Third row: Flux ratios in a sub-range of the velocities that shows excitation in the southern tail-like structure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-e76c2fd8-f37a-4a15-bbe8-3f9da90d0237.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-d3bf34d1-d992-432d-b4c3-60886867f332.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-e26fc474-6be0-4852-8337-0267a028938c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-adf14839-465a-4c19-bf43-8a366fa5e509.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-0c238d34-9823-4474-8a48-15b91091db73.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-ac8e638c-657b-44a5-8317-1ece06886731.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-de55f732-0e7b-4eab-9084-157f13fbea83.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-f7041e38-18f2-4e79-b2f7-b31ea3e6a607.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-50743b47-ccb7-4015-b590-c29727ed7759.jpg)

Fig. 16. Fit of the spectral energy distribution of 4C12.50. The infrared part of the spectrum was fit with a modified black body at long wavelengths (red dashed-dotted line) and with a power-law representing a series of modified black bodies at short wavelengths (red thick line). The total infrared model is plotted with a solid orange line. A broken power-law approximation with exponential cuto GLYPH&lt;11&gt; is used for the radio emission. Minimum and maximum synchrotron flux models are given for the radio to optical range (solid cyan curve and dotted cyan curve, respectively; see the text for more details). A maximum model for the inverse Compton emission in the X-rays is shown with a dashed line.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01011v1-figure-917c4996-c158-47a2-8bc5-d65926a66f5c.jpg)
C. M. Fotopoulou K. M. Dasyra F. Combes P. Salomé M. Papachristou
Complex molecular gas kinematics in the inner 5 kpc of 4C12.50 as seen by ALMA
The nearby system 4C12.50, also known as IRAS 13451+1217 and PKS 1345+12, is
a merger of gas-rich galaxies with infrared and radio activity. It has a
perturbed interstellar medium (ISM) and a dense configuration of gas and dust
around the nucleus. The radio emission at small ($\sim$100 pc) and large
($\sim$100 kpc) scales, as well as the large X-ray cavity in which the system
is embedded, are indicative of a jet that could have affected the ISM. We
carried out observations of the CO(1-0), (3-2), and (4-3) lines with the
Atacama Large Millimeter Array (ALMA) to determine basic properties (i.e.,
extent, mass, and excitation) of the cold molecular gas in this system,
including its already-known wind. The CO emission reveals the presence of
gaseous streams related to the merger, which result in a small
($\sim$4kpc-wide) disk around the western nucleus. The disk reaches a
rotational velocity of 200 $kms^{-1}$ , and has a mass of
3.8($\pm$0.4)$\times$10${^9}M_{\odot}$. It is truncated at a gaseous ridge
north of the nucleus that is bright in [O III]. Regions with high-velocity CO
emission are seen at signal-to-noise ratios of between 3 and 5 along filaments
that radially extend from the nucleus to the ridge and that are bright in [O
III] and stellar emission. A tentative wind detection is also reported in the
nucleus and in the disk. The molecular gas speed could be as high as 2200
$kms^{-1}$ and the total wind mass could be as high as
1.5($\pm$0.1)$\times$10$^9M_{\odot}$. Energetically, it is possible that the
jet, assisted by the radiation pressure of the active nucleus or the stars,
accelerated clouds inside an expanding bubble.