{"id":"25c13307-78bd-45c2-bf6e-402e75db97f8","arxiv_id":"2507.23100","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"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.","lead":"A deep Chandra look at the quasar PKS 0023-26 shows hot gas hugging the radio jet in the galaxy's center and no sign of a surrounding galaxy cluster. The finding is a data point for how young radio jets heat and reshape their host galaxies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 0.9 keV thermal component is identified without testing an alternative soft power-law component, and the hot-cocoon and Mach-number claims would not follow if such a component fit equally well.","rationale":"The stress-test confirms the reader's conditional verdict, with one refinement. The most load-bearing assumption is not the PSF modeling but the spectral identification of the soft excess as a thermal plasma. The morphology evidence is relatively robust: the surface brightness profile exceeds the simulated PSF at radii 0.4-5 arcsec, and the 2D Gaussian fit gives a deconvolved size much larger than the 0.07 arcsec blur, so even a modest PSF error would not erase the extended emission. In contrast, the distinction between a 0.9 keV apec component and a steep power-law soft excess is a classic degeneracy in low-count AGN spectra, especially when the AGN continuum is heavily absorbed and no strong lines are individually detected. The paper's own model list omits this alternative, so the reported p-value, while correctly computed against the absorbed power-law model, does not establish the apec interpretation. The proposed refit is a small, decisive test using existing data. If the apec model remains strongly preferred, the paper's central claim stands; if not, the interpretation should be reduced to an extended non-thermal or photoionized soft X-ray component. The paper is otherwise careful and the data are valuable, so the conditional recommendation is appropriate.","tokens_in":23924,"tokens_out":8160,"duration_ms":91239,"concrete_test":"Refit the combined 160 ks spectrum with the alternative model phabs*(zphabs*(powlaw1d + zgauss) + powlaw1d2), with the second power law free of intrinsic absorption or with its own small column, and compare Cstat to the apec model (model 5 in Table 2) using the same energy range and the same simulation-based significance procedure. Specifically, compute Delta Cstat between the second-power-law model and the apec model when both have the same number of free parameters; also inspect residuals around 0.56 keV (O VII) and 0.7-1.1 keV (Fe-L) for line structure that would favor a thermal plasma. If the power-law soft component matches the apec within statistical noise, the thermal interpretation should be treated as unproven and the central claim downscoped.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the soft excess in the central r < 7 kpc spectrum is genuinely thermal (apec, kT ~ 0.9 keV). The model set in Section 2.2 includes one absorbed power law for the AGN continuum, Galactic and intrinsic absorption, a Fe-Kalpha Gaussian, and then the apec component. It does not include an alternative soft component such as an unabsorbed or lightly absorbed second power law, which in a Type 2 quasar could represent scattered nuclear emission, a jet, or photoionized gas. Table 2 shows that adding the apec improves Cstat from 458.6/441 (model 3) to 449.2/439 (model 5), and the paper reports a p-value < 0.0002 from simulations against the no-apec model. That test does not discriminate between a thermal plasma and a steep power law (e.g., Gamma ~ 3-4) over 0.3-2 keV with only ~690 net counts and no individually detected thermal lines. If a second power law fits the data as well as the apec, then the hot-gas cocoon, jet-ISM shock heating, and Mach number M ~ 1.75-2 (Section 4.4) are not established; the extended morphology alone cannot prove thermal emission. The paper notes a two-temperature model is unconstrained (Section 3.3) but does not report the single-power-law soft-excess comparison that addresses the relevant degeneracy.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a deep (160 ks) Chandra ACIS-S observation of the compact steep-spectrum radio source PKS 0023-26, a Type 2 quasar at z=0.322. The authors report a non-detection of large-scale diffuse X-ray emission out to ~60 kpc, placing an upper limit of L(0.5-2 keV) < 3.2e42 erg/s and concluding that the source is not in a rich cluster but at most in a poor, low-temperature group. Within the central r<7 kpc region, the X-ray spectrum is modeled with an absorbed power law, a Fe-Kalpha line, and a thermal apec component with kT ~ 0.9 keV. The X-ray morphology is found to be elongated, aligned with the radio source and extending beyond it, with enhanced emission near the northern radio lobe coincident with the CO(3-2)/CO(2-1) line-ratio peak. The authors interpret this as a jet-driven shock heating the ISM, estimate a Mach number M ~ 1.75-2, and discuss implications for jet-ISM feedback and the molecular gas environment.","tokens_in":24343,"tokens_out":4996,"duration_ms":67855,"significance":"If the central claims hold, the paper provides a valuable multi-wavelength case study of a young radio source interacting with its ISM, connecting X-ray hot gas, radio morphology, and ALMA molecular gas data. The non-detection of a cluster-scale atmosphere is an important environmental constraint for a source previously suggested to reside in a rich cluster. The analysis is careful in several respects: it uses CHART/MARX PSF simulations with EDSER repositioning and BLURR=0.07 arcsec, reports p-value simulations for the significance of the Fe-Kalpha line and the apec component, and includes MCMC posterior checks. However, the central interpretive step - that the soft excess is genuinely thermal - is not fully established, and the quoted Mach number is not measured independently of that assumption. The observational material is strong; the interpretation needs additional testing or appropriate caveats.","major_comments":[{"comment":"The claim that the central r<7 kpc emission contains a thermal (apec) component is not tested against the most relevant alternative: a second, soft power-law component. The model set in Sec. 2.2 includes only one power law for the AGN, and Table 2 shows that adding apec improves Cstat from 458.6/441 (model 3) to 449.2/439 (model 5). With only ~690 net counts and no individually detected thermal lines, a steep power law (e.g., Gamma~3-4) over 0.3-2 keV could plausibly produce a similar improvement. The reported p-value<0.0002 tests apec versus no apec, not apec versus a soft power law. I ask the authors to fit an alternative model with a second unabsorbed or lightly absorbed power law replacing the apec component, report the Cstat, the best-fit photon index, and a simulation-based p-value for the apec versus power-law comparison. Without this test, the hot-cocoon, shock-heating, and Mach-number conclusions do not follow uniquely from the data.","section":"Sec. 3.3, Table 2"},{"comment":"The Mach number M ~ 1.75-2 is not an independent measurement but follows almost entirely from the assumption that kT=0.9 keV is the post-shock temperature. Using v_shock = sqrt(16 kT / 3 mu m_p) for the shock velocity and c_s = sqrt(gamma kT / mu m_p) for the post-shock sound speed with gamma=5/3 gives M = sqrt(16/5) ~ 1.79 by construction, regardless of the actual shock dynamics. The paper itself notes that the shocks are not resolved in X-rays. This should be presented as a consistency estimate under the strong-shock, fully-ionized assumption, not as a measured Mach number; the abstract and conclusions should be worded accordingly. A direct measurement would require a spatial temperature or density discontinuity or an independent shock velocity measurement.","section":"Sec. 4.4"},{"comment":"The extended-emission and cocoon conclusions rely on the accuracy of the CHART/MARX PSF simulations, specifically on the use of BLURR=0.07 arcsec. The profile comparison in Fig. 3 and the 2D Gaussian fit in Fig. 4 are suggestive, but I ask for a robustness test: vary the intrinsic PSF broadening over a plausible range and verify that the excess at radii >0.32 arcsec and the deconvolved Gaussian widths remain significant. This would strengthen confidence that the elongation in Figs. 5 and 6 is not an artifact of PSF modeling. If such tests are already available, they should be reported explicitly.","section":"Sec. 3.2 and Sec. 2.1.1"}],"minor_comments":[{"comment":"The text says \"difusse\" emission; this should be \"diffuse.\"","section":"Sec. 3.1"},{"comment":"The uncertainties on the apec temperature are given as 90% in Table 2 but as 1-sigma in the text (kT=0.91+0.19-0.37 keV). Please make the confidence levels consistent.","section":"Sec. 3.3"},{"comment":"For the sector fits, the table provides parameter values but no fit statistics (Cstat/d.o.f.) or p-values for the differences between sectors; adding these would help the reader assess the significance of the claimed north-south asymmetry and the weak Fe-Kalpha line in sector 3.","section":"Sec. 3.4, Table 3"},{"comment":"The emission-line identifications for the second X-ray source assume z=0.322 without a spectroscopic redshift for that object; this should be stated more prominently so that the line identifications are not over-interpreted.","section":"Appendix B"},{"comment":"The hot-gas mass estimate depends on the assumed cocoon volume and on the assumed density distribution. The text does note the dependence, but it would be helpful to give the electron-density or emission-measure value that corresponds to the quoted mass and luminosity.","section":"Sec. 4.3"}],"recommendation":"major_revision","confidential_remarks":"The Chandra data and the environmental non-detection are solid and valuable. The main issue is interpretive: the thermal nature of the soft excess and the Mach number derivation need to be either strengthened with additional spectral fitting or substantially caveated. I do not think the paper should be rejected, but the central claims need to be made robust against the soft-power-law degeneracy before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is a solid, careful observational paper, not a breakthrough. The new Chandra data on PKS 0023-26 are well reduced, and the two headline results hold up: there is no bright cluster-scale X-ray atmosphere around this CSS quasar, and there is a real extended soft X-ray component within ~7 kpc that aligns with the radio source and the ALMA CO(3-2)/CO(2-1) peak. The non-detection, the PSF simulations, and the p-value/MCMC checks are done properly.\n\nThe main soft spot is the spectral interpretation. The central claim is that the soft excess is thermal, kT~0.9 keV, and they quote M~1.75-2 for a shock. But with ~690 net counts and no resolved lines, an additional steep power-law (e.g., scattered nuclear emission) could fit the 0.3-2 keV band equally well. They compare apec against no apec, but not apec against an extra power law. That's a real omission given the source is a Type 2 quasar. If a second power law fits, the hot cocoon and Mach number don't follow. The extended morphology is suggestive but not sufficient to prove thermal emission.\n\nAlso, the Mach number is not a measurement. It comes from assuming kT=0.9 keV is the post-shock temperature and applying strong-shock jump conditions. That's a reasonable estimate, but the abstract presents it as a result. Similarly, the abstract says the environment is 'consistent only with' a poor group; the data actually just put an upper limit that also allows no group at all. These are wording issues, not fatal flaws.\n\nI disagree with the reader's verdict somewhat: I think the paper is more solid than a 'conditional' verdict implies for the observational parts. The environment result is important for the CSS/cluster connection, and the alignment with CO line ratios is a nice multi-wavelength step. But the thermal interpretation needs the alternative model test before the paper is fully convincing. I'd send it to a good referee with a request to add that test and tone down the Mach-number language. The authors are experienced; they'll know exactly what's being asked.","headline":"Careful new Chandra data on PKS 0023-26, with a robust environment non-detection and a plausible but under-tested thermal interpretation.","tokens_in":24839,"tokens_out":2775,"would_cite":true,"duration_ms":33417,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"PKS 0023-26's X-rays trace a jet-shocked hot cocoon, not a cluster-scale atmosphere.","keywords":["active galaxies","radio galaxies","radio jets","compact steep spectrum sources","X-ray astronomy","jet-ISM interaction","molecular gas","galaxy groups"],"falsifier":"A re-analysis of the same Chandra data that, after PSF subtraction, finds the central soft excess is consistent with a steep power law rather than a 0.9 keV thermal plasma, or that removes the extended surface-brightness excess when a different PSF blur is used, would overturn the cocoon interpretation. Likewise, detecting diffuse X-ray emission above 3e42 erg/s within a 60 kpc radius would overturn the poor-group conclusion.","tokens_in":23764,"feed_emoji":"🔭","tokens_out":9874,"duration_ms":108242,"temperature":0.7,"pith_summary":"PKS 0023-26 is a compact, powerful radio source inside a quasar, and earlier optical work left open whether it sits in a galaxy cluster or a group. This paper used a 160-kilosecond Chandra exposure to answer that question, and the answer is neither the simple cluster picture nor the simple group picture: there is no bright cluster-scale hot atmosphere out to about 60 kpc, but there is genuine hot gas within the central 7 kpc. That gas has a temperature near 0.9 keV, its shape is elongated along the radio source and extends beyond the lobes, and its brightest patch falls on the northern radio lobe where the molecular gas shows its highest CO(3-2)/CO(2-1) ratio. The paper reads these facts together as evidence that the young jet is driving a supersonic shock, with Mach number roughly 1.75-2, into the interstellar medium and heating a cocoon of gas around the radio source. If this is right, PKS 0023-26 becomes a resolved example of jet-ISM feedback inside a modest environment, and the same interaction is what raises the CO excitation.","feed_headline":"X-ray cocoon, not cluster, surrounds quasar PKS 0023-26","feed_subtitle":"A 160-kilosecond Chandra exposure finds hot 0.9 keV gas aligned with the radio lobes and a Mach 2 shock.","key_machinery":"The load-bearing object is the resolved soft X-ray cocoon around the radio source, isolated by subtracting the simulated Chandra point-spread function from the image and modeled spectrally as an apec component added to the absorbed AGN. Three measurements work together: the excess surface brightness over the simulated PSF at radii of roughly 0.4-5 arcsec; the 0.9 keV thermal component in the central spectrum; and the spatial alignment of the soft X-ray residual with the 87 GHz radio contours and the ALMA CO ratio map. The cocoon geometry gives a volume of about 52 $kpc^{3}$ and a hot-gas mass near 3.8e8 solar masses, and comparing the 0.9 keV post-shock temperature with the sound speed yields the supersonic Mach number of 1.75-2.","core_discovery":"The central discovery is that PKS 0023-26's kiloparsec-scale X-ray emission is dominated by a hot, collisionally ionized plasma associated with the radio source rather than by an extended cluster atmosphere. In the r<7 kpc spectrum, an absorbed power-law AGN plus a Fe-Kalpha line leaves a soft excess that requires an apec thermal component with kT=0.$91^{{+0.19}}$_{-0.37} keV and a 0.5-2 keV luminosity near 3.8e42 erg/s. The surface brightness profile is broader than the simulated point-spread function beyond about 0.32 arcsec, and the residual emission is elongated with a deconvolved 1-$\\sigma$ width of roughly 0.36 by 0.24 arcsec along the radio axis, matching the separation of the radio lobes. Soft-band X-rays extend beyond the radio source and perpendicular to it, with a clear enhancement at the northern lobe that coincides with the peak of the CO(3-2)/CO(2-1) line ratio and with the region where CO(1-0) is suppressed. The paper concludes that this is a cocoon of gas shock-heated by the expanding jet, with a Mach number of 1.75-2, and that the absence of large-scale diffuse emission rules out a rich cluster, allowing at most a poor group with kT<0.5 keV.","pith_inferences":["The authors leave implicit that if the cocoon is partly pre-existing hot gas from a poor group, the same observation strategy applied to other compact steep-spectrum sources could map jet-driven heating on 1-10 kpc scales even where no cluster is present.","The combination of the X-ray peak and the CO line-ratio peak predicts that even higher-J CO transitions and warm H2 emission should be brightest at the northern lobe; pointed ALMA observations of CO(6-5) or near-infrared H2 could test this directly.","If part of the 0.9 keV soft excess is photoionized rather than shock-heated, the true Mach number would be lower than 1.75-2; deeper, sector-resolved Chandra spectra of the northern lobe could separate the two emission mechanisms.","A hard X-ray observation above 10 keV would measure the intrinsic AGN power law and the reflection responsible for Fe-Kalpha, tightening the estimate of the thermal excess and its shock origin."],"forward_implications":["If the claim is correct, the optical companions of PKS 0023-26 are not evidence of a rich cluster; the source lives in a poor group or less, and its abundant cold molecular gas was probably accreted from merging companions rather than cooled from a cluster atmosphere.","The aligned X-ray morphology makes PKS 0023-26 a resolved example of a young compact radio source shocking its host ISM, alongside systems such as 3C171 and 3C305, extending the jet-feedback picture to a quasar with ALMA-resolved molecular gas.","The coincidence of the X-ray peak with the CO(3-2)/CO(2-1) peak and the CO(1-0) lacuna supports the idea that jet-driven X-ray irradiation raises CO to higher rotational levels and suppresses the lowest transition near the nucleus.","The Mach 1.75-2 shock implies the jet compresses the ISM by roughly a factor of two, which is the mechanism invoked to explain the wrapped morphology of the molecular gas around the radio lobes.","The upper limit on large-scale thermal emission, below about 3e42 erg/s, sharpens the environment census: powerful quasars in luminous X-ray clusters are rare, and PKS 0023-26 fits that statistical pattern."],"supporting_citations":[{"why":"Supplies the ALMA CO morphologies, the CO(3-2)/CO(2-1) line-ratio map, and the CO(1-0) lacuna that the X-ray peak is compared against.","marker":"T. Oosterloo et al. 2025, Paper I"},{"why":"Provides the ALMA detection of the radio core, the molecular gas geometry, and the jet power used in the outflow discussion.","marker":"R. Morganti et al. 2021"},{"why":"Gives the earlier XMM-Newton luminosity that the paper's central luminosity measurement and environment reassessment must match.","marker":"B. Mingo et al. 2014"},{"why":"Provides the X-ray luminosity profile of a poor cluster used to argue that a cluster atmosphere around PKS 0023-26 would have been detected.","marker":"E. O'Sullivan et al. 2010"},{"why":"Provides a rich, cool-core group comparison showing that such a group would also have been visible in the Chandra spectrum.","marker":"E. O'Sullivan et al. 2017"},{"why":"Supplies the 3C171 aligned X-ray/radio shock-heating result that motivates the same interpretation for PKS 0023-26.","marker":"M. J. Hardcastle et al. 2010"},{"why":"Supplies the 3C305 example of jet-driven X-ray emission used as a direct comparison for the cocoon interpretation.","marker":"M. J. Hardcastle et al. 2012"},{"why":"Provides the theoretical basis for X-ray irradiation altering CO excitation and explaining the CO(1-0) lacuna.","marker":"P. R. Maloney et al. 1996"},{"why":"Supports the adopted 0.07 arcsec intrinsic blur in the Chandra PSF simulations used to establish the extended emission.","marker":"J. Ma et al. 2023"},{"why":"Provides the redshift, the warm outflow velocity, and the ionized-gas properties used in the shock and mass-outflow estimates.","marker":"F. Santoro et al. 2020"}],"fun_headline_variants":["Chandra shows shock-heated cocoon around quasar PKS 0023-26","Quasar PKS 0023-26 cocooned by hot gas, not a cluster","Jet-driven X-ray cocoon explains quasar's hot gas","PKS 0023-26: X-rays trace jet shock, not cluster halo","Hot X-ray cocoon from jet shock in PKS 0023-26"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Chandra point-spread function, simulated with an intrinsic blur of 0.07 arcsec, truly reproduces the telescope's response, and that the soft X-ray excess inside 7 kpc is genuinely a hot plasma rather than an additional power-law component; if either part gives way, the extended cocoon and the Mach-number estimate are not established.","fun_headline_variants_meta":{"raw":{"variants":["Chandra shows shock-heated cocoon around quasar PKS 0023-26","Quasar PKS 0023-26 cocooned by hot gas, not a cluster","Jet-driven X-ray cocoon explains quasar's hot gas","PKS 0023-26: X-rays trace jet shock, not cluster halo","Hot X-ray cocoon from jet shock in PKS 0023-26"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001052,"raw_usage":{"total_tokens":4548,"prompt_tokens":1206,"completion_tokens":3342,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":822,"completion_tokens_details":{"reasoning_tokens":3234}},"tokens_in":822,"tokens_out":3342,"duration_ms":26384,"temperature":1.0,"reasoning_tokens":3234,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:04:07.150985+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A re-analysis of the same Chandra data that, after PSF subtraction, finds the central soft excess is consistent with a steep power law rather than a 0.9 keV thermal plasma, or that removes the extended surface-brightness excess when a different PSF blur is used, would overturn the cocoon interpretation. Likewise, detecting diffuse X-ray emission above 3e42 erg/s within a 60 kpc radius would overturn the poor-group conclusion.","supporting_citations":[{"cited_title":"The changing impact of radio jets as they evolve: The view from the cold gas","cited_arxiv_id":"2506.20448","evidence_quote":"Supplies the ALMA CO morphologies, the CO(3-2)/CO(2-1) line-ratio map, and the CO(1-0) lacuna that the X-ray peak is compared against."}],"review_version":1}