REVIEW 2 major objections 6 minor 1 cited by
The changing impact of radio jets as they evolve: The view from the cold gas
T0 review · 2 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Extreme CO line ratios reveal where a radio jet is reshaping cold gas in PKS 0023-26.
desk verdict Solid, careful ALMA line-ratio work on one radio galaxy; the central observational results hold up, but the multi-kpc cocoon interpretation is honestly flagged by the authors as not unique because of the ongoing merger. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central diagnostic is the set of brightness-temperature line ratios R21 and R32 built from three CO transitions observed at matched spatial and spectral resolution. These ratios trace molecular gas excitation and optical depth; extreme values mark gas heated or stirred by the jet and its expanding cocoon, letting the authors map jet influence independently of kinematics. The cocoon of shocked gas created by the jet-ISM interaction, predicted by simulations to expand perpendicular to the jet, provides the mechanism that carries the interpretation of high ratios several kiloparsecs away from the radio axis.
What would settle it
Spatially resolve CO(1-0), CO(2-1), and CO(3-2) in gas-rich galaxy mergers that lack active nuclei; if regions with R32 > 0.6 extending about 2 to 5 kpc are found in such mergers, the cocoon interpretation for the outer high-excitation region of PKS 0023-26 would need to be abandoned in favour of merger turbulence.
Extended reading notes
Core claim
The paper establishes that in PKS 0023-26, a powerful young radio galaxy, the radio jet modifies the cold molecular interstellar medium in ways that are clearly visible in CO excitation and morphology. The brightness-temperature ratio R32 (CO(3-2)/CO(2-1)) peaks at values of about 1.0 to 1.2 in the region between the core and the northern lobe, and the region with R32 > 0.6 extends roughly 2 to 5 kpc beyond the radio source, with R21 > 1.9 in the central beam where CO(1-0) is undetected. These ratios are far above those of normal star-forming galaxies (R21 around 0.6 and R32 between 0.2 and 0.6). The authors interpret the coincidence of the extreme ratios with the radio emission as evidence of jet-ISM interaction: a destructive interaction at the northern lobe depletes molecular gas there, while a slowly expanding cocoon of shocked gas drives turbulence and high excitation over a larger volume, including perpendicular to the jet. They conclude that AGN impact can be traced through line ratios and gas distribution as well as kinematics, and that most of the jet energy escapes the galaxy.
Load-bearing premise
The load-bearing premise is that the jet cocoon, not the ongoing galaxy merger, produces the high CO excitation beyond the northern lobe and perpendicular to the jet; if merger-driven turbulence alone can produce the same line ratios, the cocoon interpretation for the outer high-excitation region would be unsupported.
Editorial extensions
If this is right
- Surveys of AGN feedback that rely only on gas kinematics will miss jet influence; spatially resolved mapping of multiple CO transitions is a complementary and sometimes more sensitive tracer.
- In young radio galaxies one can expect a spatial stratification of feedback: ejective in the central kiloparsec, destructive at the radio lobes, and a milder, cocoon-driven enhancement of turbulence and excitation on multi-kiloparsec scales.
- The small fraction of jet power transferred to the cold gas (kinetic power below about 0.05 percent of the jet power) suggests most of the energy escapes to the circumgalactic medium, so jets may affect galaxy halos more than host ISMs.
- High CO excitation in regions oriented perpendicular to radio jets can serve as a signpost for a growing jet cocoon in gas-rich galaxies, extending the diagnostic beyond individual objects.
- The coexistence of different interaction modes within one source implies that spatially resolving the full extent of the radio source is necessary to quantify AGN feedback reliably.
Reading between the lines
- If the cocoon interpretation is correct, other young radio galaxies of similar power should show high R32 regions elongated perpendicular to their jet axes; a statistical study of resolved line ratios could test this pattern.
- The paper's own merger caveat suggests a clean control: spatially resolved CO line-ratio imaging of gas-rich galaxy mergers without active nuclei would show whether high excitation regions beyond the radio source are unique to jet-ISM coupling or a generic merger phenomenon.
- The non-detection of CO(1-0) at the core could point to X-ray or cosmic-ray heating suppressing low-J CO emission; deep observations of higher-J CO or isotopologues at high angular resolution could separate the heating mechanisms.
- The measured low coupling efficiency implies that a full energy budget requires multi-phase measurements—hot X-ray gas, warm ionised gas, and cold molecular gas—rather than a single tracer, to quantify where jet energy is ultimately deposited.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents new ALMA CO(1-0) and CO(3-2) observations of the young radio galaxy PKS 0023−26 (z=0.32188), combined with the existing CO(2-1) data of Morganti et al. (2021). The data are matched in resolution and velocity binning, absorption is removed from the CO(1-0) cube via clean components, and the same source mask derived from the highest-S/N transition is applied to all three transitions to avoid S/N-dependent biases in line ratios. The main observational results are: (1) CO(1-0) is not detected at the radio core while CO(2-1) and CO(3-2) are bright there, implying R21 ≳ 1.9 in the central beam; (2) the CO(3-2) emission shows a broad central profile (FWZI ~600 km/s) and high velocity dispersion between the core and the northern lobe; (3) molecular gas is depleted at the northern lobe and appears wrapped around it; and (4) the R32 = I(3-2)/I(2-1) ratio peaks near 1.0–1.2 between the core and the N lobe and remains above 0.6 over a region extending up to ~5 kpc from the source, including an extension beyond the N lobe. The authors interpret the inner region as an AGN-driven outflow and the outer high-ratio gas as evidence of a jet cocoon inducing turbulence and enhanced excitation, while explicitly acknowledging that turbulence from the ongoing galaxy merger is a viable alternative for part of the outer region.
Significance. If the interpretation holds, the paper is a valuable demonstration that AGN feedback can be traced through molecular excitation and gas distribution, not only kinematics, on scales from the central kiloparsec to several kiloparsec, and that different feedback modes (ejective and preventative) can coexist in one source. The observational work is careful: the use of a common mask derived from the highest-S/N transition, the removal of absorption before imaging, and the quoting of extreme ratios against 3-sigma upper limits are all sound practices. The line ratios are measured rather than fitted, and the energetic estimates are explicitly parameterized with literature conversion factors and stated geometric assumptions. The paper also openly identifies the main degeneracy—merger-induced turbulence versus jet-cocoon excitation—and calls for non-AGN merger control observations. The central-kpc line-ratio anomalies and the gas depletion around the N lobe are robust; the multi-kpc cocoon interpretation is plausible but not uniquely established.
major comments (2)
- [§4.1, §5, Fig. 8] The most novel interpretation—that the extended R32 > 0.6 region reaching ~2–5 kpc beyond the N lobe and perpendicular to the radio axis is produced by the expanding jet cocoon—is not uniquely supported. The manuscript itself states in §4.1 that turbulence from accreted gas settling in the merger 'might explain the region of enhanced excitation beyond the N lobe' and in §5 that this region is 'difficult to reconcile' with the cocoon picture. Because the galaxy is demonstrably interacting (tidal tails, companions at similar redshift), spatial coincidence with the radio axis is not a decisive discriminator. If merger-driven turbulence alone produces R32 > 0.6 over several kpc, the multi-kpc cocoon claim falls, although the central-kpc anomalies and N-lobe depletion remain valid. Please either provide a quantitative test separating the mechanisms (e.g., comparison with tidal-tail morphology, modeling of merger-driven excitation, or a non-AGN merger control sample) or consistently downgrade the abstract and §3.3 to present the cocoon as one of several viable explanations.
- [§3.3, Fig. 8] The R32 map is presented without an uncertainty map or an explicit significance threshold. The statements that R32 reaches 1.0–1.2 between the core and the N lobe and that R32 > 0.6 extends several kpc require a quantitative definition of which pixels are considered detected in both transitions. Please add noise-based confidence contours, or state the S/N threshold used to construct the ratio map and to define the R32 > 0.6 region; without this, the spatial extent of the high-excitation region cannot be assessed.
minor comments (6)
- [§3.3] The notation 'R21 /greaterorsimilar1.9' should be typeset as R21 ≳ 1.9.
- [§4.3] The object 'IC 5053' appears to be a typo for 'IC 5063' (compare §4.1).
- [§4.3] The phrase 'in a ways that can go beyond' should read 'in ways that can go beyond'.
- [Table 1, §2] The three transitions have different maximum recoverable scales (4.7, 3.8, and 2.5 arcsec for CO(1-0), CO(2-1), and CO(3-2)); please comment on whether the shorter maximum recoverable scale of the CO(3-2) data affects the line-ratio measurements in the outer low-surface-brightness regions, noting that any missing CO(3-2) flux would bias R32 low, not high.
- [§4.1, References] The companion paper (Siemiginowska et al. 2026) is cited as providing X-ray evidence for enhanced hot gas at the N lobe and the R32 peak; since it is listed as submitted, the supporting X-ray interpretation should be marked as preliminary.
- [§4.1] The phrase 'the latterly expanding cocoon' is unusual; consider 'the laterally expanding cocoon'.
Circularity Check
No significant circularity: the key quantities are measured line ratios and gas distributions, and the jet-impact interpretation is an empirical correlation with an explicitly considered merger alternative.
full rationale
The paper's central claims rest on directly measured ALMA quantities: CO(1-0), CO(2-1), and CO(3-2) integrated intensities, the R32 = I(3-2)/I(2-1) and R21 brightness-temperature ratios, velocity dispersions, and the spatial correspondence of these with the radio continuum. None of these is derived from the conclusion; no parameter is fitted to a subset of the data and then renamed as a prediction. The statement that extreme line ratios in a region aligned with the radio axis indicate conditions different from the undisturbed gas is an observational contrast, not a definitional identity: the ratios are measured independently of the radio-axis alignment, and the undisturbed-gas comparison uses values in the same galaxy at large radii. The non-detection of CO(1-0) at the core is an observational upper limit; calling it 'extreme' via R21 > 1.9 is a calibrated statement of the same detection, not a fitted input. Energetic estimates use literature conversion factors (e.g., alpha_CO = 0.89 or 4.3 M_sun/(K km/s pc^2)) and explicitly stated assumptions about v_out and r_out; these do not feed back into the detection claims. The paper's most interpretive step, the attribution of the outer high-R32 region to a jet cocoon, is not circular because the manuscript itself presents the merger-turbulence alternative in Sect. 4.1 and notes in the Conclusions that the region beyond the N lobe is 'difficult to reconcile' with the cocoon explanation. Self-citations, notably Morganti et al. (2021), which supplies the CO(2-1) data and the earlier scenario, are cumulative: the new CO(1-0) and CO(3-2) observations test that scenario independently. A minor reference to 'Morganti et al. in prep.' for similar scales in other powerful AGN is ancillary and not load-bearing. No uniqueness theorem or ansatz is imported from the authors' prior work. Overall, the derivation chain is self-contained with respect to the measurements; the main vulnerability, if any, is the causal attribution of the outer excitation to the cocoon rather than the merger, which is a scientific alternative the authors explicitly acknowledge, not a circularity.
Assumptions & free parameters
free parameters (4)
- CO-to-H2 conversion factor alpha_CO =
0.89 M_sun/(K km/s pc^2) (ULIRG-type) or 4.3 (Galactic)
- Outflow velocity v_out =
200 to 300 km/s
- Outflow radius r_out =
~1 kpc
- Excitation temperature Tex from absorption =
5 +/- 1 K
assumptions (3)
- domain assumption CO line ratios in undisturbed gas of normal star-forming galaxies are R21 about 0.6 +/- 0.2 and R32 about 0.2 to 0.6
- domain assumption The radio source age is about 1e5 yr from the low-frequency spectral turnover
- domain assumption Hydrodynamic simulations of jets in clumpy media (Mukherjee et al. 2016, 2018a) describe the cocoon and cloud-shredding processes
Cite this review
Pith. "Pith review of The changing impact of radio jets as they evolve: The view from the cold gas." pith.science (2026). https://pith.science/paper/FAL6TSVZ
@misc{pith2026250620448,
author = {Pith},
title = {Pith review of: The changing impact of radio jets as they evolve: The view from the cold gas},
year = {2026},
howpublished = {\url{https://pith.science/paper/FAL6TSVZ}},
note = {Machine review of arXiv:2506.20448}
}
read the original abstract
We present ALMA CO(1-0) and CO(3-2) observations of a powerful young radio galaxy, PKS 0023-26, hosted by a far-infrared bright galaxy. The galaxy has a luminous optical AGN and a very extended distribution of molecular gas. We used these observations (together with available CO(2-1) data) to trace the impact of the AGN across the extent of the radio emission and beyond on scales of a few kpc. Despite the strength of the optical AGN, the kinematics of the cold molecular gas is strongly affected only in the central kpc, and is more weakly affected around the northern lobe. We found other signatures of the substantial impact of the radio AGN, however. Most notably, extreme line ratios of the CO transitions in a region aligned with the radio axis indicate conditions very different from those observed in the undisturbed gas at large radii. The non-detection of CO(1-0) at the location of the core of the radio source implies extreme conditions at this location. Furthermore, on the scale of a few kpc, the cold molecular gas appears to be wrapped around the northern radio lobe. This suggests that a strong jet-cloud interaction has depleted the northern lobe of molecular gas, perhaps as a result of the hot wind behind the jet-induced shock that shreds the clouds via hydrodynamic instabilities. The higher gas velocity dispersion and molecular excitation that we observed close to this location may then be the result of a milder interaction in which the expanding jet cocoon induces turbulence in the surrounding interstellar medium. These results highlight that the impact of an AGN can manifest itself not only in the kinematics of the gas, but also in molecular line ratios and in the distribution of the gas. Although the radio plasma and the cold molecular gas are clearly coupled, the kinetic energy that is transferred to the ISM is only a small fraction of the energy available from the AGN.
Figures
Figures from the paper (5 more)
Forward citations
Cited by 1 Pith paper
-
Chandra X-ray Observatory study of the X-ray emission of PKS 0023-26 and comparison with recent ALMA results
PKS 0023-26 lacks an X-ray bright galaxy cluster but shows a small, hot gas cocoon aligned with its radio jet, pointing to jet-driven shock heating of the host ISM.
Reference graph
Works this paper leans on
-
[1]
Molecular gas excitation and outflow properties of obscured quasars at z$\sim$0.1
Audibert, A., Combes, F., García-Burillo, S., et al. 2019, A &A, 632, A33 Audibert, A., Ramos Almeida, C., García-Burillo, S., et al. 2023, A&A, 671, L12 Audibert, A., Ramos Almeida, C., García-Burillo, S., et al. 2025, arXiv:2505.02759 Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207 . Callingham, J. R., Ekers, R. D., Gaensler, B. M., et a...
work page Pith review arXiv 2019
-
[11]
Figure A.1 shows the two absorption profiles plotted in optical depth τ. The profiles are extracted from the high spatial- and velocity resolution CO(1-0) and CO(2-1) cubes, given th e small velocity width of the absorption. −40 −20 0 20 40 V (km s−1) −0.02 0.00 0.02 0.04 0.06 0.08 0.10 0.12 τ PKS 0023−26 CO(1-0) CO(2-1) Fig. A.1. Profile of the absorption i...
work page 2021
-
[13]
all show the same structure of a core and two lobes (se e Fig. 3 and Morganti et al. (2021)). Table B.1 and Fig. B.1 give the integrated flux for the three components and their spectr al indices α (defined through S ∝ να) as derived from these fluxes. The spectral indices of the lobes are quite steep and steeper than seen at lower frequencies (Morganti et al...
work page 2021
- [21]
-
[24]
2011, ApJ, 736, 62 Westmeier, T., Kitaeff, S., Pallot, D., et al
Wang, J., Fabbiano, G., Elvis, M., et al. 2011, ApJ, 736, 62 Westmeier, T., Kitaeff, S., Pallot, D., et al. 2021, MNRAS, 506, 3962 Y ajima, Y ., Sorai, K., Miyamoto, Y ., et al. 2021, PASJ, 73, 257 Y oung, L. M., Meier, D. S., Crocker, A., et al. 2022, ApJ, 933, 90 Zakamska, N. L., Strauss, M. A., Krolik, J. H., et al. 2006, AJ , 132, 1496 Zovaro, H. R. M....
work page 2021
-
[79]
2002, A&A, 395, 13 Mingo, B., Hardcastle, M
Mellema, G., Kurk, J.D., Röttgering, H.J.A. 2002, A&A, 395, 13 Mingo, B., Hardcastle, M. J., Croston, J. H., et al. 2011, ApJ , 731,
work page 2002
-
[95]
Mingo, B., Hardcastle, M. J., Croston, J. H., et al. 2014, MNR AS, 440, 269 Miura, R. E., Espada, D., Hirota, A., et al. 2021, MNRAS, 504,
work page 2014
-
[174]
Carilli, C. L. & Walter, F. 2013, ARA&A, 51, 105 Cicone, C., Severgnini, P ., Papadopoulos, P . P ., et al. 2018, ApJ, 863, 143 Costa-Souza, J. H., Ri ffel, R. A., Souza-Oliveira, G. L., et al. 2024, ApJ, 974, 127 Cresci, G., Tozzi, G., Perna, M., et al. 2023, A&A, 672, A128 Croston, J. H., Kraft, R. P ., & Hardcastle, M. J. 2007, ApJ, 66 0, 191 Dasyra, K....
work page 2013
Show all 13 references
-
[293]
1988, MNRAS, 235, 403 Tadhunter, C., Holt, J., González Delgado, R., et al
Tadhunter, C., Fosbury, R., de Serego Alighieri, S., Bland, J., Danziger, I., Goss, W., Snijders, M. 1988, MNRAS, 235, 403 Tadhunter, C., Holt, J., González Delgado, R., et al. 2011, M NRAS, 412, 960 Tadhunter, C., Morganti, R., Rose, M., et al. 2014, Nature, 5 11, 440 Tamhane...
1988
-
[705]
Sutherland, R. S. & Bicknell, G. V . 2007, Ap&SS, 311,
2007
-
[1550]
S., Tadhunter, C
Ramos Almeida, C., Bessiere, P . S., Tadhunter, C. N., et al. 2013, MNRAS, 436, 997 Rose, T., McNamara, B. R., Combes, F., et al. 2024, MNRAS, 533 , 771 Roy, N., Heckman, T., Overzier, R., et al. 2024, ApJ, 970, 69 Ruffa, I., Prandoni, I., Davis, T. A., et al. 2022, MNRAS, 510,...
2013
-
[2021]
observed against the small, compact component in the N lobe which has been de- tected in VLBI and that has a size of about 100 pc Tzioumis et al. (2002). The results from Gaussian fits to the CO(1-0) an d CO(2-1) profiles are given in Table A.1. From the detection of absorption ...
2002
-
[6198]
2023, A&A, 673, A13 Morganti, R., Oosterloo, T., Oonk, J
Montoya Arroyave, I., Cicone, C., Makroleivaditi, E., et al . 2023, A&A, 673, A13 Morganti, R., Oosterloo, T., Oonk, J. B. R., Frieswijk, W., T adhunter, C. 2015, A&A, 580, A1 Morganti, R., Oosterloo, T., Tadhunter, C., et al. 2021, A&A , 656, A55 Mukherjee, D., Bicknell, G. V...
2023
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.