Gas, dust, and star formation in the positive AGN feedback candidate 4C 41.17 at z=3.8

We present new, spatially resolved [CI]1-0, [CI]2-1, CO(7-6), and dust continuum observations of 4C 41.17 at $z=3.8$ obtained with the IRAM NOEMA interferometer. This is one of the best-studied radio galaxies in this epoch and is arguably the best candidate of jet-triggered star formation at high redshift currently known in the literature. 4C 41.17 shows a narrow ridge of dust continuum extending over 15 kpc near the radio jet axis. Line emission is found within the galaxy in the region with signatures of positive feedback. Using the [CI]1-0 line as a molecular gas tracer, and multifrequency observations of the far-infrared dust heated by star formation, we find a total gas mass of $7.6\times 10^{10}$ M$_{\odot}$, which is somewhat greater than that previously found from CO(4-3). The gas mass surface density of $10^3$ M$_{\odot}$ yr$^{-1}$ pc$^{-2}$ and the star formation rate surface density of 10 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$ were derived over the 12 kpc$\times$8 kpc area, where signatures of positive feedback have previously been found. These densities are comparable to those in other populations of massive, dusty star-forming galaxies in this redshift range, suggesting that the jet does not currently enhance the efficiency with which stars form from the gas. This is consistent with expectations from simulations, whereby radio jets may facilitate the onset of star formation in galaxies without boosting its efficiency over longer timescales, in particular after the jet has broken out of the interstellar medium, as is the case in 4C 41.17.


Introduction
It is now widely recognized that powerful radio jets ejected by powerful active galactic nuclei (AGN) can have a major impact on the evolution of their host galaxies by regulating their gas content and interstellar medium conditions, thereby potentially facilitating the formation of new stars.Outflows of molecular, atomic, and warm ionized gas found in numerous galaxies are our prime observational evidence of feedback today.
Other potential mechanisms, by which the injection of AGN energy may regulate star formation in the host galaxies, are much less explored.This includes the enhancement of gas turbulence and pressure within the interstellar gas and star-forming clouds.A pressure increase on embedded gas clouds may temporarily facilitate star formation ("positive feedback"; e.g., Fragile et al. 2004Fragile et al. , 2017;;Wagner et al. 2012;Silk 2013).The existence of such positive feedback has been expected theoretically for several decades, but observational evidence is hard to collect, and is only found for very few individual targets in the nearby and distant universe (Croft et al. 2006;Salomé et al. 2015Salomé et al. , 2017)).Maiolino et al. (2017) suggest that velocity offsets found in a IRAM data are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/639/L13 Based on data obtained with the IRAM NOEMA interferometer through program S19DA.
nearby ultraluminous infrared galaxy (ULIRG) may also point to positive feedback associated with radio-quiet quasar activity.
An interesting question is whether positive feedback can boost the efficiency with which stars form from molecular gas beyond the levels found in intensely star-forming galaxies, as postulated by some models of galaxy evolution (Gaibler et al. 2012;Silk 2013;Mukherjee et al. 2018).Shocks driven into the gas by the radio jet may compress the denser regions of these clouds and lead to star formation from gas that was previously marginally Jeans-stable (e.g., Wagner et al. 2012;Fragile et al. 2017).The rare candidates at low redshift with signatures of positive feedback often seem to fall below the relationship between molecular gas mass and star formation rate density, suggesting these galaxies form stars at lower efficiency than those without AGN (Salomé et al. 2017).However, our clearest examples of positive feedback in the nearby Universe are either dwarf galaxies (e.g., Minkowski's Object, 3C 285, Salomé et al. 2015) or the outskirts of massive early-type galaxies (e.g., Cen A, Salomé et al. 2017).These examples have low gas-mass surface densities compared to actively star-forming galaxies, and therefore they do not allow us to infer whether or not AGN may boost star formation in environments dominated by dense molecular gas.Mukherjee et al. (2018) found from relativistic hydrodynamic simulations that the efficiency of star formation in such galaxies likely depends strongly on the detailed properties of the gas clouds, jet power, and interaction geometry, making it very difficult to predict theoretically whether the net effect of the radio jet on star formation is positive or negative.
To test whether positive feedback may enhance the star formation efficiency in massive galaxies in the early Universe, we observed the powerful radio galaxy 4C 41.17 at z = 3.792 (Chambers et al. 1990) with the NOEMA interferometer of the Institut de Radioastronomie Millimetrique (IRAM), probing the [CI] lines of atomic carbon and the dust continuum at ∼500 GHz rest-frame frequency.We observed [CI] rather than CO because this line is the more reliable tracer of molecular gas mass at high redshift, it is optically thin, and is less affected by gas excitation effects than the more commonly observed mid-J CO lines (e.g., Papadopoulos & Greve 2004;Valentino et al. 2018;Nesvadba et al. 2019).4C 41.17 is arguably our best example of positive AGN feedback in the early Universe, as shown by the alignment of extended, unpolarized Lyα emission with the radio jet axis (Dey et al. 1997), strongly kinematically perturbed gas (Steinbring 2014), and asymmetric absorption-emission spectra typical of very young stellar populations (P Cygni profiles; Dey et al. 1997).Bicknell et al. (2000) present a detailed analytical study of the potential gas compression by the radio jet, finding that at least parts of the extended Lyα emission associated with the radio jet axis and bright extended UV continuum may originate from star formation in gas compressed through the passage of the radio jet.A young starburst age of 30 Myr is also supported by multicomponent SED fitting (Rocca-Volmerange et al. 2013) and is consistent with the age estimate of the radio jet (Bicknell et al. 2000).
In the following we describe our observations of the gas, dust, and star formation in 4C 41.17, and discuss our results in the light of star formation in massive galaxies at high redshift with and without radio jets.We use the flat ΛCDM cosmology from Planck Collaboration VI ( 2020) with H 0 = 67.4km s −1 Mpc −1 , Ω M = 0.315, and Ω Λ = 1−Ω M .At z = 3.792 the luminosity distance is D L = 34.41Gpc, and 7.26 kpc are projected onto one arcsecond.
2. NOEMA interferometry of 4C 41.17 4C 41.17 was observed between June and December 2019 with the NOEMA interferometer in the 2 mm and 3 mm bands with either 9 or 10 antennae in the C configuration under average to good atmospheric conditions.The central frequencies were 102.705 and 168.895GHz, respectively, corresponding to the observed frequencies of [CI]1-0 and [CI]2-1 at z = 3.79198, respectively.The source was observed for a total of 7.0 h and 5.6 h in the 2 mm and 3 mm bands, respectively.The data were calibrated with the Continuum and Line Interferometer Calibration package CLIC 1 and with MWC349 and LKHA101 as flux 1 http://www.iram.fr/IRAMFR/GILDAS
calibrators.Combining both polarizations leads to a continuum root mean square (rms) of 7.3 µJy bm −1 and 23.9 µJy bm −1 in the 3 mm and 2 mm bands, respectively, measured over 3.8 GHz bandwidth.We rebinned the data to a spectral resolution of 20 MHz in both bands (corresponding to 58 km s −1 in band 3, and 36 km s −1 in band 2, respectively), and obtained a line rms of 0.1 and 0.3 mJy bm −1 and per spectral channel, respectively.With natural weighting, we obtain a beam size of 2.3 × 1.7 , and 1.4 × 1.0 at 3 mm and 2 mm, respectively.Figure 1 shows the resulting integrated spectra.
Figure 2 shows the [CI] and dust morphology.All line emission is concentrated within an area of full width at half maximum (FWHM) 2.8 × 1.7 around the radio core for [CI]1-0, and 1.8 × 1.1 for the [CI]2-1 and CO(7-6) lines, respectively.The major axis of the line emission exceeds the major axis size of the beam significantly in both bands, clearly showing that the emission is resolved along this axis, but not extended enough to neglect beam smearing.Gaussian deconvolution suggests major axis lengths of 1.6 and 1.3 for the [CI]1-0 and the [CI]2-1 and CO(7-6) lines, respectively.
The continuum, best seen in the 2 mm band, is strongly elongated and more extended along the major axis than the line emission observed in the same band, has an observed size of 2.4 × 1.1 , and a deconvolved size of 2.0 .For simplicity, and because our main goal is to rule out that the star formation rate surface density exceeds those measured in other galaxies, we adopt a common size of 1.6 × 1.1 for the area emitting [CI]1-0 and the dust continuum.
The continuum major axis is roughly aligned with the radio jet axis.Integrated continuum flux densities, which are measured after subtracting the line emission from the data cubes, are 488 ± 24 and 112 ± 7 µJy in the 2 mm and 3 mm band, respectively.Integrated line and continuum fluxes are listed in Table 1.As already noted by De Breuck et al. (2005), the nonthermal synchrotron emission from the radio source does not significantly change these measurements.
To extract line fluxes, we fit the [CI]1-0 line with a single Gaussian centered on a velocity of (228 ± 21) km s −1 relative to z = 3.97197, a FWHM line width of (372 ± 21) km s −1 and a flux of I CI10 = (0.26 ± 0.03) Jy km s −1 .

Gas and star formation properties
We used Eq.(4) of Alaghband-Zadeh et al. (2013), , where D L is the luminosity distance in units of Gpc, z the redshift, I CI the integrated line flux of [CI]1-0 in Jy km s −1 .The Einstein A coefficient, A 10 = 7.93 × 10 −8 is given in units of 10 −7 s −1 , and the carbon abundance, X CI = 3 × 10 −5 (Alaghband-Zadeh et al. 2013), in units of 10 −5 .We find a molecular gas mass of M H2 = (7.6 ± 1.0) × 10 10 M , which is somewhat greater than the 5.4 × 10 10 M previously found by De Breuck et al. (2005) from CO(4-3) with a 7 × 5 beam; this is a sufficient to cover most of the surrounding Lyα halo.This suggests that we have not resolved more extended line emission associated with 4C 41.17.
We estimated the star formation rates from the dust continuum in the NOEMA data at 2 mm and 3 mm and used previous Herschel/SPIRE and ground-based single-dish measurements taken from the compilation of Drouart et al. (2014) to determine the dust temperature, T D and the β parameter.We only used observations taken at ≥350 µm, where the AGN contamination is <5% in the spectral energy distribution (SED) of dust shown by Drouart et al. (2014).This suggests T D = (53 ± 3) K and β = 1.9, respectively, for a modified blackbody fit calculated in the same way as in Canameras et al. (2015).This SED overpredicts the observed flux densities by factors 2.5 and 2 in the 2 mm and 3 mm bands, respectively (Table 1).This may either be due to missing flux in the interferometry data or additional, fainter dust emission within the much larger beam size (10−20 ) of the single-dish observations.A potential candidate for additional far-infrared (FIR) emission would be diffuse emission from the large, 5 × 7 , Lyα halo previously found by van Breugel et al. (1999).In the following we use the fluxes obtained with NOEMA, which are best matched to the [CI] emission-line region rather than those derived from the singledish measurements.
We adopted a modified version of the relationship of Kennicutt (1998a) to estimate a star formation rate from the FIR luminosity.Using the Chabrier stellar initial mass function this implies SFR [M yr −1 ] = 2.5 × 10 −44 L FIR [erg s −1 ].We find a total star formation rate of SFR = 650 M yr −1 .This is a factor ∼2 greater than the SFR = 290 M yr −1 of star formation that van Breugel et al. (1999) found for the UV continuum in the same region, which is likely affected by dust extinction.Averaged over a region of 1.6 × 1.1 (11.6 kpc × 8 kpc at z = 3.79), this corresponds to star formation and gas mass surface densities of Σ SFR = 8.9 M yr −1 kpc −2 and 1040 M yr −1 pc −2 , respectively.Since we measured and resolved [CI] and the dust continuum with similar beam sizes, these results are robust, even if some of the discrepancy between the continuum measurements with NOEMA and single-dish telescopes should be due to missing flux instead of additional astrophysical sources.

Discussion and summary
Gas mass and star formation rate surface density in galaxies are closely related; they have a near-linear relationship that holds over more than six orders of magnitude in star formation and for molecular gas mass surface densities >10 M pc −2 (e.g., Kennicutt 1998b;Bigiel et al. 2008).Actively star-forming galaxies at high redshift tend to fall above that relationship, as first noted by Daddi et al. (2010) and Genzel et al. (2010), indicating shorter gas consumption times and potentially higher star formation efficiency.
Figure 3 shows where 4C 41.17 falls in this Schmidt-Kennicutt diagram between star formation and gas mass surface density, averaged over the 11.6 kpc × 8 kpc area previously defined in Sect.2, and the gas mass and star formation rate from Sect. 3. We also add the second galaxy with potential evidence of jet-triggered star formation in the early Universe, PKS 0529−549 at z = 2.6, to this analysis.This galaxy was previously discussed by Man et al. (2019), who point out that star formation in this source seems to be more efficient than in giant molecular clouds in the Milky Way.Both galaxies fall well above the typical relationship found for disk galaxies in this redshift range, which are on the main sequence of star formation, as already noted by Man et al. (2019).However, the star formation rates do not fall above those found in other intensely star-forming galaxies at the same cosmic epoch (Bothwell et al. 2010), which likely include the progenitors of powerful highredshift radio galaxies.
This result suggests that the star formation efficiency does not exceed that found in other strongly star-forming galaxies without AGN, at least when averaged over the ∼10 kpc regions where jettriggered star formation was found (Dey et al. 1997).We might L13, page 3 of 4 Fig. 3. Gas and star formation rate surface densities of 4C 41.17 (red diamond), PKS 0529−549 (blue dot), and several comparison samples: Giant molecular clouds in the Milky Way (Evans et al. 2014), local starburst galaxies (Kennicutt 1998a), star-forming galaxies at z = 1−2 (Tacconi et al. 2008;Valentino et al. 2018Valentino et al. , 2020)), and submillimeter galaxies (Bothwell et al. 2010).The solid purple and blue lines indicate the sequences for typical main-sequence and starburst galaxies of Daddi et al. (2010).
suspect that positive and negative feedback would both be present and partially cancel the effect on each other (Dugan et al. 2017;Mukherjee et al. 2018).Steinbring (2014) found outflow signatures in warm ionized gas in 4C 41.17, suggesting that negative feedback is at work, although the low line width of [CI]1-0 does not suggest that a large percentage of the cold neutral gas is currently being removed.However, in situ gas removal through winds would also lower the observed gas-mass surface densities.This would move the galaxy to the left in Fig. 3 and artificially enhance, not lower, the observed star formation efficiency.
A more likely explanation is that the phase of positive feedback is very short and only triggers the onset of star formation without boosting it over longer timescales.This is not unexpected from a theoretical point of view: initial hydrodynamic simulations suggested that jet-triggered star formation may be long-lived (Gaibler et al. 2012), but did not take into account energy and momentum losses through disk porosity.More recent work suggests that the impact of the jet on the star formation drops after a short initial "kick-off" phase (Fragile et al. 2017;Mukherjee et al. 2018).This would also be consistent with the very high star formation rates of 7300 M yr −1 required to form the 2.2 × 10 11 M of young stars found by Rocca-Volmerange et al. ( 2013) within 30 Myr, that is, an order of magnitude higher than today.Figure 2 also shows that the jet has already broken out of 4C 41.17, suggesting it may no longer be able to sustain a high overpressure within the gas.Jet-triggered star formation may also be limited by star formation itself once the starburst has evolved past the first few million years, as previously suggested by Fragile et al. (2017).At 30 Myr age (Rocca-Volmerange et al. 2013), the starburst in 4C 41.17 is clearly evolved enough to have produced the first generations of supernovae.These may make star formation in 4C 41.17 self-regulating, as has previously been suggested for other massive, intensely star-forming galaxies at z = 2−3 (e.g., Lehnert et al. 2013).
To summarize, we used new [CI]1-0 and 2-1, CO, and dust continuum interferometry obtained with IRAM/NOEMA to study star formation in 4C 41.17 at z = 3.8, which is the best example of positive feedback from radio jets at high redshift in the literature.We estimated a molecular gas mass of 7.6 × 10 10 M in an elongated region where signatures of jettriggered star formation have previously been found.The spatially resolved gas and star formation rate densities suggest that stars are formed at efficiencies that are comparable to submillimeter galaxies at the same epoch.The same is true for PKS 0529−549, the second galaxy with likely signatures of a very young, potentially jet-triggered, stellar population.
While our results do not call into question the previous good evidence for jet-triggered star formation in 4C 41.17 obtained in the UV and radio, they do suggest that this phase is very short and not sustained after the jet has broken out from the interstellar gas of the host galaxy.After few 10 Myr, jet-triggered star formation proceeds with similar efficiency as in other intensely star-forming galaxies, making the signatures of positive feedback difficult to identify observationally.This result is also consistent with recent observational and theoretical findings in low-redshift galaxies, which indicate that jet-triggered star formation does not necessarily lead to a boost in star formation efficiency, even at high gas-mass surface densities of up to 1000 M pc −2 .

Fig. 1 .
Fig. 1.Integrated spectra of [CI]1-0 (left),and the [CI]2-1 and CO(7-6) lines (right).The red and orange lines show the total fit to the line profile and the fits to the individual components, respectively.

Fig. 2 .
Fig.2.Dust (left) and [CI]1-0 morphology (right) of 4C41.17.Contours show the 1.4 GHz radio morphology.The radio core is between the two central components(Bicknell et al. 2000).The beam size is shown in the lower left corner of each panel.