{"id":"ff771728-9b13-4683-a407-735e30ed66cc","arxiv_id":"2508.06707","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"JWST/NIRSpec IFU observations of six z~3.5-4 radio galaxies show jet-aligned ionized outflows with kinetic powers of 10^43.2-10^45 erg/s, implying significant jet-driven feedback at early epochs.","lead":"This paper maps the kinematics and energetics of ionized gas outflows in six powerful radio galaxies at z=3.5-4 using JWST/NIRSpec. It finds that radio jets drive large-scale outflows carrying 0.15-2% of the AGN's bolometric luminosity, but only about 1% of jet mechanical power couples to the warm gas.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tabulated data contradict the headline 0.15–2% outflow efficiency; recomputed values are 0.005–0.2%, undermining the 'sufficient feedback' claim.","rationale":"The reader's weakest_assumption focused on the L_jet–L_radio calibration and noted an internal inconsistency between the 0.15–8% and 0.15–2% ranges, but did not identify that the paper's own tables yield 0.005–0.2%, more than an order of magnitude lower. This concern is more fundamental because it does not depend on any external calibration or assumed scaling; the headline efficiency and the 'sufficient to impact galaxy evolution' conclusion are directly contradicted by the manuscript's reported measurements. If the recomputation confirms these ratios, the central quantitative claim is unsupported and the paper's interpretation must change from 'sufficient' to 'energetically sub-dominant in the ionized phase,' even though the morphological and kinematic evidence for jet-driven gas disturbance remains valuable. I therefore recommend rejection as currently written, with the possibility of resubmission after correcting the efficiency calculation and revising the conclusions accordingly.","tokens_in":35958,"tokens_out":7361,"duration_ms":81987,"concrete_test":"Recompute epsilon_kin = KE_outflow / L_AGN for all six sources using Table 2 col. 7 and Table 1 col. 8, and recompute KE_outflow/L_jet using Table 1 col. 9. If the ratios are 0.005–0.2% and up to 2% respectively, the abstract/conclusion efficiency numbers are wrong and the 'sufficient to impact' conclusion must be revised downward. Also trace the 0.15–2% and 0.15–8% values to their derivation (check units of L_AGN, possible factor-100 normalization) to see which table or text needs correction.","verdict_should_be":"REJECT","load_bearing_attack":"The central quantitative claim — outflow kinetic power = 0.15–2% L_bol — is contradicted by the paper's own tables. §4.2 states epsilon_kin ~0.15–8%; abstract/conclusions say 0.15–2%. But dividing Table 2 KE_outflow by Table 1 L_AGN gives: TGSS 0.025%, TNJ0205 0.017%, 4C03 0.20%, 4C19 0.20%, TNJ0121 0.015%, TNJ1338 0.005% — i.e., 0.005–0.2%, one to two orders of magnitude lower than all reported ranges. Since the 'sufficient to impact galaxy evolution' assertion depends on exceeding the ~0.5% feedback threshold, none of the tabulated sources meet it. The <1% jet-coupling result is also affected: 4C03 has KE_outflow/L_jet ≈ 2% from the same tables, not <1%. This is not an external calibration issue; the text and tables are internally irreconcilable. The claim as stated is unsupported by the presented measurements.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents JWST/NIRSpec IFU observations of six powerful radio galaxies at z ≈ 3.5–4, mapping the morphology, kinematics, and energetics of warm ionized gas. It derives W80 line widths, electron densities from [SII], ionized gas masses, outflow rates, kinetic powers, and compares outflow energetics with AGN bolometric luminosity and with a newly calibrated L_jet–L_1.4GHz relation. The main claims are that all sources show large-scale jet-aligned outflows, that outflow kinetic power reaches 0.15–2% of L_bol (sufficient for feedback), and that only ≲1% of jet mechanical power is transferred to warm ionized gas. The kinematic data and spatial correlations are valuable, but the quantitative headline numbers are internally inconsistent with the paper's own tables.","tokens_in":36253,"tokens_out":5641,"duration_ms":65980,"significance":"If the quantitative results are correct, this would be the largest uniform IFU sample of jet-driven feedback at z > 3.5, with direct implications for models of AGN feedback and galaxy quenching. Strengths include the publicly available data, the consistent methodology across six targets, the spatially resolved nature of the maps, and the comparison with radio-quiet quasars. However, the central quantitative claims are not supported by the paper's own tabulated values, and key systematic assumptions are not propagated into the energetics. The current version requires major revision before the conclusions can be relied upon.","major_comments":[{"comment":"The quoted kinetic efficiencies do not match the tabulated numbers. Dividing Table 2 KE_outflow by Table 1 L_AGN gives log10(KE/L_bol) = -3.60, -3.77, -2.69, -2.70, -3.82, -4.33 for TGSS, TNJ0205, 4C03, 4C19, TNJ0121, TNJ1338 respectively, i.e. 0.025%, 0.017%, 0.20%, 0.20%, 0.015%, 0.005%. The abstract and conclusions state 0.15%–2%, while §4.2 states 0.15%–8%. The maximum value from the tables is 0.2%, below the ~0.5% feedback threshold discussed in §4.2. This internal inconsistency directly undermines the 'sufficient to impact galaxy evolution' claim.","section":"§4.2, Tables 1 and 2"},{"comment":"The claim that 'only 0.1–1%' of jet mechanical power couples to the warm ionized gas is contradicted by 4C03: from Table 1 log Ljet = 46.65 and Table 2 log KE = 44.96, giving KE/Ljet ≈ 2%, not <1%. The abstract's '≲1%' is thus not uniformly supported. Since the coupling-efficiency headline depends on this ratio, the paper should present the per-source ratio and revise the stated range accordingly, or explain why 4C03 is treated differently.","section":"§4.4, Figure 23, Table 1 and 2"},{"comment":"The new L_jet–L_1.4GHz scaling relation is calibrated on low-redshift X-ray cavity samples and then applied to z ≈ 3.5–4 radio galaxies. The HzRGs lie near the upper end of the calibration range, and the relation implicitly assumes no redshift evolution in jet composition, environment, or cavity energetics. The quoted scatter in the calibration (intercept uncertainty ±1.23 dex) is not propagated into the KE/Ljet ratios or into the N17 re-analysis. Because the main coupling conclusion is a ratio involving Ljet, the authors should show how the result changes under plausible systematic shifts (e.g., ±0.5 dex) or with independent jet-power estimators.","section":"§4.3, Eq. (5)"},{"comment":"The outflow masses and kinetic powers rest on several unquantified assumptions: (i) electron densities from single-Gaussian fits to [SII], which likely trace the narrow, non-outflow gas rather than the high-velocity component; (ii) clumping factor Ce = 1; (iii) solar oxygen abundance; (iv) for TGSS1530, Av = 0.8 borrowed from another source. These enter the mass and power linearly or, for extinction, exponentially. A sensitivity analysis is needed to show whether the 0.15–2% and <1% conclusions are robust to reasonable variations in these assumptions.","section":"§3.1.3, Eqs. (2) and (4)"}],"minor_comments":[{"comment":"W50 is used in Eq. (1) and later in §3.1.3 for v_out, but the paper only defines W80. Please define W50 or replace with W80 for consistency.","section":"§2.2, Eq. (1)"},{"comment":"Typo: 'TGSS150' should be 'TGSS1530'.","section":"§2.2"},{"comment":"The text describes an enhancement 'around ~25 kpc', but the listed outflow size for TGSS is 21 kpc. Please check the radial distance scale and the axis origin.","section":"Figure 6, bottom-left panel"},{"comment":"The text says 'only a small fraction (∼0.3%)' of jet energy is in the emisson-line outflows, while the abstract and §4.4 quote ≲1%. These numbers should be stated consistently.","section":"§4.5"},{"comment":"Table 1 states L_AGN is from [OIII]5007, but TGSS1530 lacks [OIII] coverage. Please specify what estimator is used for TGSS (e.g., Hα-based) and whether this affects the reported L_AGN.","section":"Table 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The kinematic analysis is a strong observational contribution, and the data appear to be of good quality. However, the internal arithmetic inconsistency between the text and Tables 1–2 is serious and must be resolved before the paper is suitable for publication. The authors should either correct the reported efficiencies and revise the feedback interpretation, or provide a detailed reconciliation of the numbers. A sensitivity analysis for the mass/energy assumptions would also substantially strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth reading for the data: spatially resolved W80 maps, the new detailed analysis of TGSS1530, and the systematic comparison of six z~4 radio galaxies. The alignment between broad-line regions and radio lobes is convincingly shown, and the reduced cubes are public. That part is solid.\n\nThe problem is the energetics. The abstract and conclusions say the outflow kinetic power is 0.15–2% of L_bol, but Table 2 divided by Table 1 gives 0.005–0.2%. The stress-test note is correct—I checked the arithmetic. Only 4C03 reaches 0.2% and everything else is an order of magnitude lower. Section 4.2 also quotes 0.15–8%, so the paper is internally inconsistent at three different levels. Since the “sufficient to impact galaxy evolution” argument depends on crossing the ~0.5% threshold, the claim is unsupported by their own measurements. This is not a calibration nuisance; it's a load-bearing arithmetic error.\n\nOther soft spots: the L_jet–L_radio relation (Eq 5) is fitted to low-z cavity samples and extrapolated without showing scatter or validating at high-z. The outflow mass estimates rely on solar abundance, clumping factor unity, and a single extinction for TGSS1530—all standard, but none propagated into error bars. The shortfalls are worth noting, but they are secondary to the internal inconsistency.\n\nWho benefits: anyone working on jet-driven feedback at high redshift will want the kinematic maps and the TGSS1530 analysis. But cite it for the data, not for the efficiencies.\n\nRecommendation: this deserves peer review, but it needs a major revision to reconcile the text with the tables and to temper the conclusions accordingly. The observational core is valuable; the quantitative summary as written is not correct.","headline":"The new IFU kinematics for six z>3.5 radio galaxies are valuable and the TGSS1530 analysis is genuinely new, but the headline outflow-efficiency range is contradicted by the paper's own tables, so the main feedback claim needs major revision before it can be taken seriously.","tokens_in":36798,"tokens_out":3279,"would_cite":false,"duration_ms":34510,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Radio jets at z~4 drive galaxy-scale outflows of warm ionized gas, but couple less than ~1% of their mechanical energy into that gas.","keywords":["AGN feedback","radio jets","high-redshift galaxies","ionized gas outflows","JWST NIRSpec IFU","jet-ISM coupling","W80 kinematics","Ljet-Lradio scaling"],"falsifier":"Measure the jet power of these six galaxies independently of Equation (5)—for example from X-ray cavity enthalpy, inverse-Compton lobe emission, or lobe minimum-energy arguments—and compare with the radio-derived values; a discrepancy of more than the calibration scatter would overturn the $<1\\%$ coupling claim. Separately, deep X-ray observations could test whether a hot $T\\sim10^7$ K phase with $E_{\\rm th}\\sim100\\times$ the ionized outflow kinetic energy actually exists; its absence would falsify the paper's proposed energy budget.","tokens_in":35811,"feed_emoji":"🔭","tokens_out":9505,"duration_ms":98840,"temperature":0.7,"pith_summary":"The paper uses JWST/NIRSpec integral-field spectroscopy of six powerful radio galaxies at $z \\approx 3.5$--$4$ to ask whether radio jets can regulate galaxy growth during the peak epoch of galaxy formation. It finds that warm ionized gas is turbulent on kiloparsec scales in every source, with non-parametric line widths $W_{80} \\approx 950$--$2500$ km/s, and that the broadest, fastest gas lies along the radio jet axis. From resolved maps it derives ionized gas masses of $\\sim 1$--$8\\times10^9\\,M_\\odot$, mass outflow rates of $80$--$950\\,M_\\odot\\,\\mathrm{yr}^{-1}$, and kinetic powers of $10^{43.2}$--$10^{45.0}$ erg/s, corresponding to $0.15\\%$--$2\\%$ of the AGN bolometric luminosity—enough, by common feedback thresholds, to affect the host galaxy. It then compares these outflows to jet powers from a new radio-luminosity–jet-power calibration and concludes that less than $\\sim 1\\%$ of the jet's mechanical energy couples to the warm ionized gas; the rest is argued to reside in shock-heated hot gas or in the shock ionization that powers the emission-line nebulae.","feed_headline":"At z~4, radio jets drive outflows but couple under 1% of energy","feed_subtitle":"Six JWST galaxies show jet-aligned turbulence strong enough to affect galaxy evolution.","key_machinery":"The analysis is carried by spatially resolved $W_{80}$ maps—a non-parametric line-width measure enclosing 80% of the line flux—built from the NIRSpec IFU cubes of $[{\\rm O\\,III}]\\(\\lambda5007\\) (or H$\\$\\alpha$$), together with radial profiles computed in 1-kpc boxes along the projected radio axis. These maps are converted into resolved outflow properties through a standard prescription that uses $[{\\rm S\\,II}]$ doublet ratios for electron density and the $|v|>500$ km/s high-velocity channel to isolate the outflow component, yielding gas mass, $\\dot{M}_{\\rm out}$, momentum flux, and kinetic power $\\dot{E} = \\frac12 \\dot{M} $v^{2}$$. The coupling-efficiency conclusion additionally leans on a newly fi","core_discovery":"On the paper's own terms, the discovery is a dual one. All six galaxies show broad $[{\\rm O\\,III}]\\(\\lambda5007\\) (or H$\\$\\alpha$$ for TGSS1530) emission with $W_{80}$ from about 950 to more than 2500 km/s extending over tens of kiloparsecs, and the highest line widths, highest outflow rates, and peak kinetic powers are spatially coincident with radio cores or lobes. This spatial coincidence is the direct evidence for jet-driven feedback: the gas is not merely disturbed, it is disturbed along the jet's direction. Quantitatively, the outflows carry $0.15\\%$--$2\\%$ of $L_{\\rm bol}$, which the paper argues is sufficient to impact galaxy evolution, yet only $\\approx0.1$--$1\\%$ of the inferred jet m","pith_inferences":["If the $L_{\\rm jet}$--$L_{1.4\\,\\rm GHz}$ relation carries large intrinsic scatter or evolves with redshift, the quoted $<1\\%$ coupling efficiency will shift; the qualitative conclusion that jet energy is not dominated by warm-ionized bulk motion is more robust, since even order-of-magnitude changes in $L_{\\rm jet}$ leave $\\dot{E}_{\\rm out}\\ll L_{\\rm jet}$.","A direct observational test is to measure the hot gas thermal energy of these same galaxies with deep X-ray observations: if no $T\\sim10^7$ K phase with $E_{\\rm th}\\sim100\\times$ the ionized outflow kinetic energy is present, the paper's proposed energy sink must be replaced.","Applying the same $W_{80}$ and outflow methodology to a matched sample of radio-quiet AGNs at $z\\sim4$ would isolate the jet contribution; the paper's comparison already suggests jets produce a directional turbulence component absent in radio-quiet systems, but a controlled sample would make that quantitative.","If the compact-source efficiency trend is confirmed with larger samples, it would imply feedback is most effective before jets break out of the dense interstellar medium, shifting attention to young compact steep-spectrum sources as the principal feedback agents at early epochs."],"forward_implications":["Radio jets at $z>3.5$ are capable of redistributing hundreds of solar masses per year and depositing kinetic power at or above the $\\sim0.5\\%$ $L_{\\rm bol}$ threshold that feedback models associate with quenching, so jet-driven feedback is a viable mechanism during the peak epoch of galaxy formation.","Because only $\\lesssim1\\%$ of the jet mechanical energy appears in warm ionized bulk motion, ionized-outflow studies alone systematically underestimate the feedback energy budget; the dominant reservoir should be a hot, X-ray-emitting phase with $E_{\\rm th}\\sim 100\\times$ the ionized outflow kinetic energy.","The apparent anti-correlation between $\\dot{E}_{\\rm out}/L_{\\rm jet}$ and radio source size implies compact, young jets couple most efficiently to the dense ISM, so early-stage jet feedback may be the most important for stirring the host galaxy.","Shock models with $v_{\\rm shock}\\approx600$--$1000$ km/s can produce the observed emission-line luminosities with shock front areas of 20--80 kpc$^2$, comparable to the outflow areas, so the line-emitting gas may itself be a calorimeter of jet energy loss rather than a passive tracer."],"supporting_citations":[{"why":"Establishes the spatially resolved outflow methodology applied here and provides the TNJ1338 kinematics and energetics that anchor the sample; the present analysis extends its per-source approach to six galaxies.","marker":"Roy et al. (2024)"},{"why":"Provides the JWST/NIRSpec IFU data and initial detection of extended nebular gas for four of the six high-redshift radio galaxies.","marker":"Wang et al. (2024)"},{"why":"Supplies adopted systemic redshifts, radio core positions, and the companion-galaxy interpretation that this paper argues against for some sources.","marker":"Wang et al. (2025)"},{"why":"Hydrodynamic simulations of jets in clumpy gas-rich media that predict cocoon-driven asymmetric outflows and the less-than-1% coupling efficiency the paper claims to match.","marker":"Mukherjee et al. (2018)"},{"why":"Jet-ISM interaction simulations used to interpret the low coupling efficiency and the high thermal/kinetic energy partition in the ISM.","marker":"Meenakshi et al. (2022)"},{"why":"One of the low-redshift X-ray cavity samples whose jet-power errors enter the least-squares fit of the new jet-power–radio-luminosity relation.","marker":"Rafferty et al. (2006)"},{"why":"Supplies additional X-ray cavity jet-power measurements spanning the radio luminosity range used to calibrate the scaling relation.","marker":"Cavagnolo et al. (2010)"},{"why":"Provides the [O III] 5007 luminosity-based ionized gas mass prescription adopted for the outflow mass and rate estimates.","marker":"Cano-Díaz et al. (2012)"},{"why":"Shock model grid used to translate emission-line luminosities into shock energy fluxes and required shock front areas, supporting the discussion of where jet energy goes.","marker":"Allen et al. (2008)"}],"fun_headline_variants":["Jet-driven outflows at z~4: powerful but only ~1% energy coupled","JWST reveals jet-aligned gas outflows: 0.15-2% of AGN bolometric power","High-z radio jets drive massive outflows, yet couple <1% of energy","z~4 jets stir gas: outflows carry 0.15-2% of AGN power, <1% jet coupling"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The jet mechanical powers are inferred from a 1.4 GHz radio luminosity–jet power relation calibrated on low-redshift X-ray cavity samples; if that calibration is biased or has large scatter when applied at $z\\sim4$, the headline coupling-efficiency result of $\\lesssim1\\%$ is not secure.","fun_headline_variants_meta":{"raw":{"variants":["Jet-driven outflows at z~4: powerful but only ~1% energy coupled","JWST reveals jet-aligned gas outflows: 0.15-2% of AGN bolometric power","High-z radio jets drive massive outflows, yet couple <1% of energy","z~4 jets stir gas: outflows carry 0.15-2% of AGN power, <1% jet coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000223,"raw_usage":{"total_tokens":1369,"prompt_tokens":897,"completion_tokens":472,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":641,"completion_tokens_details":{"reasoning_tokens":366}},"tokens_in":641,"tokens_out":472,"duration_ms":6199,"temperature":1.0,"reasoning_tokens":366,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T22:35:57.936476+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the jet power of these six galaxies independently of Equation (5)—for example from X-ray cavity enthalpy, inverse-Compton lobe emission, or lobe minimum-energy arguments—and compare with the radio-derived values; a discrepancy of more than the calibration scatter would overturn the $<1\\%$ coupling claim. Separately, deep X-ray observations could test whether a hot $T\\sim10^7$ K phase with $E_{\\rm th}\\sim100\\times$ the ionized outflow kinetic energy actually exists; its absence would falsify the paper's proposed energy budget.","supporting_citations":[{"cited_title":"2025, A&A, 696, A88, doi: 10.1051/0004-6361/202553668","cited_arxiv_id":null,"evidence_quote":"Supplies adopted systemic redshifts, radio core positions, and the companion-galaxy interpretation that this paper argues against for some sources."}],"review_version":1}