{"id":"4e5ad090-97ec-4170-9648-351390f83dea","arxiv_id":"2506.15286","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"JWST/MIRI MRS data show a fast ionized gas outflow in the central ~6 pc of Centaurus A, with velocities up to about 1000-1400 km/s and a mass outflow rate near 2 solar masses per year.","lead":"This paper maps gas motions in the center of the nearest giant radio galaxy, Centaurus A, with JWST. It finds a fast outflow of hot ionized gas in the inner six parsecs, while warm molecular gas shows only weak signs of outflow.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Outflow rate depends on an unmeasured 50% broad Pfα fraction; the quoted 1.6-2.9 Msun/yr does not include this uncertainty.","rationale":"The paper is a careful, multi-line MIRI-MRS study. The detection of fast ionized gas in the unresolved nuclear aperture is supported by consistent broad components in [Ar II] and [Ne III], p-v diagrams reaching roughly +1000 and -1400 km/s, enhanced velocity dispersions perpendicular to the radio jet, and shock-sensitive line ratios; this is real evidence. My concern is specifically the quantitative outflow rate: it is built on the 50% Pfα broad-fraction assumption, which the authors explicitly acknowledge is unmeasurable in Pfα itself, and the current 1.6-2.9 Msun/yr range is narrower than the true systematic uncertainty. However, this does not overturn the qualitative outflow detection, and the energy-budget conclusion is robust to a factor-of-several change in dM/dt because the jet power exceeds L_kin by roughly an order of magnitude. The reader's CONDITIONAL verdict already captures this quantitative uncertainty, so I recommend no change to the verdict. The proposed test would either validate the 50% assumption or provide a corrected, more honest range.","tokens_in":25638,"tokens_out":16143,"duration_ms":176253,"concrete_test":"Fit the nuclear Pfα line with the broad Gaussian shape (centroid and width) fixed to the [Ar II] three-Gaussian broad component, leaving only its amplitude free, and measure the broad flux fraction with a 68% confidence interval. Then recompute M_out and dM/dt using the recovered fraction. For an independent cross-check, perform a free three-Gaussian fit on [Ne II] at 12.814 um, which is bright and has no strong PAH confusion in the nuclear aperture, and compare its broad-component flux fraction and velocity width with the [Ar II]/[Ne III] values. If the recovered fraction lies outside 0.35-0.65, the paper should report a propagated dM/dt range that reflects this systematic uncertainty rather than only the two v_out,max choices.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The mass outflow rate in Sect. 3.5 rests on an unmeasured assumption: that about 50% of the Pfα flux belongs to the broad outflow component. Pfα is described as 'broad and relatively faint against the bright continuum, which complicates an accurate kinematic decomposition', so its broad fraction is taken by analogy from two-Gaussian fits to [Ar II] and [Ne III]. That analogy is not secure because recombination lines (Pfα) and collisionally excited lines ([Ar II], [Ne III]) have different density dependences, and the two-Gaussian versus three-Gaussian decompositions yield different broad-component widths and centroids (e.g., [Ar II] sigma_broad = 567 km/s in the two-Gaussian fit vs 702 km/s in the three-Gaussian fit). The quoted dM/dt = 1.6-2.9 Msun/yr interval only reflects two choices of v_out,max and not the uncertainty in the 50% fraction, nor the unresolved-nature ambiguity between a coherent outflow and a turbulent NLR or jet-cloud interaction. Because this rate is a headline quantitative result used in the energy budget comparison, the assumption should be tested directly before the number is used as a benchmark.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents JWST/MIRI MRS 5–28 μm observations of the central ~100–200 pc of Centaurus A and analyzes the kinematics of ionized fine-structure lines and warm H2 rotational lines. The ionized gas is found to be extended along the radio jet direction, with broad nuclear components (σ ≈ 600 km s−1 for [Ar II] and [Ne III]), high velocity wings reaching about +1000 and −1400 km s−1 in position-velocity diagrams, and enhanced velocity dispersions perpendicular to the jet. These features are interpreted as a jet/AGN-driven ionized outflow within the central unresolved ~6 pc region. Using the Pfα line and assuming that half of its flux belongs to a broad outflow component, the authors derive an ionized gas mass outflow rate dM/dt ≈ 1.6–2.9 M⊙ yr−1 and conclude that both the AGN radiation field and the radio jet power are sufficient to drive it. The warm H2 gas is mostly rotational with mild noncircular motions, and no fast H2 outflow is detected. The central claim is that the MRS data provide clear evidence of an ionized gas outflow with |v| > 500 km s−1 originating from the central 6 pc.","tokens_in":25799,"tokens_out":3986,"duration_ms":48669,"significance":"If the qualitative outflow detection holds, this paper is a valuable addition to the small set of spatially resolved, multi-phase studies of jet-ISM interaction in a nearby radio galaxy, and it strengthens the observational case that low-to-intermediate-power radio jets can drive ionized outflows on tens-of-parsec scales. The explicit comparison with simulations of jet-inflated bubbles, the use of multiple independent fine-structure lines, and the resolved p-v diagrams are strengths. The quantitative mass outflow rate, however, depends on an unmeasured broad-flux fraction and on a standard one-zone formula, and the unresolved nuclear geometry leaves room for alternative interpretations. The paper would be strengthened by a direct propagation of the assumed flux-fraction uncertainty and by a more cautious framing of the outflow interpretation.","major_comments":[{"comment":"The mass outflow rate rests on the assumption that approximately half of the Pfα flux belongs to the broad outflow component. This fraction is not measured: the text states that Pfα is 'broad and relatively faint against the bright continuum, which complicates an accurate kinematic decomposition,' and the 50% value is transferred from two-Gaussian fits to [Ar II] and [Ne III]. The transfer is not secure because recombination lines and collisionally excited lines have different density and temperature dependences, and the two-Gaussian versus three-Gaussian decompositions yield different broad-component parameters (e.g., [Ar II] σ_broad = 567 km s−1 in two Gaussians versus 702 km s−1 in three Gaussians; [Ne III] 306.6 versus 605.6 km s−1 in Table A.1). Since M_out and dM/dt scale linearly with this fraction, a plausible range of 20–80% would change the quoted rate by a factor of about 2.5 in each direction. The quoted interval 1.6–2.9 M⊙ yr−1 reflects only the two choices of v_out,max, not the flux-fraction uncertainty, even though the text claims that uncertainties include the fraction. I request a sensitivity analysis over the assumed broad fraction and, if possible, a direct decomposition of Pfα using the same fitting methodology as for [Ar II] and [Ne III].","section":"Section 3.5"},{"comment":"The interpretation of the unresolved broad nuclear components as a coherent outflow is not unique. The broad profiles and the high-velocity p-v features are all confined to the central unresolved region (≲0.35″ or 6 pc), so a radially expanding outflow, a turbulent AGN narrow-line region, a jet-cloud interaction, or scattered light from the AGN continuum could in principle produce similar unresolved kinematics. The paper's statement that the features are 'consistent with simulation predictions' is correct, but consistency with one class of simulations does not exclude these alternatives. Because the final summary claims 'clear evidence' of an outflow, the authors should either add diagnostics that distinguish a coherent outflow from a turbulent/broad-line region (e.g., line-ratio variations across the broad component, or a resolved velocity gradient in the broad component traced by [Ne III] where it is slightly more extended) or temper the claim to 'consistent with an outflow' throughout the abstract and summary.","section":"Section 3.4.2 and Section 4"},{"comment":"The statement that the mass outflow rate is a lower limit because the nuclear region is unresolved is not self-evident from the adopted formula dM/dt = 3 × M_out × v_out,max / R_out. If R_out is taken as the beam FWHM (6 pc) but the true emitting region is smaller, the rate would be higher, as stated; however, if part of the broad component is not outflowing gas, or if the assumed 50% flux fraction is an overestimate, the rate would be lower. The lower-limit claim is therefore contingent on the flux-fraction and geometry assumptions. Please specify the assumed physical radius explicitly and discuss how the result depends on it, or remove the unqualified 'lower limit' statement.","section":"Section 3.5"}],"minor_comments":[{"comment":"'The latter might be to related to gas streamers' contains a typo ('to related' should be 'related').","section":"Abstract"},{"comment":"'we switched offthe background corrections' is missing a space; also the phrase 'switched off' is colloquial for a pipeline setting and could be rephrased as 'disabled'.","section":"Section 2.1"},{"comment":"In the sentence 'as infered from the neon line ratios,' 'infered' should be 'inferred'.","section":"Section 3.3"},{"comment":"The label 'Pf' in the ch1 spectrum is ambiguous; using 'Pfα' in the figure would match the text and table.","section":"Figure 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of A&A and addresses a timely topic. The main quantitative result (dM/dt ≈ 1.6–2.9 M⊙ yr−1) is currently presented with an unquantified and possibly dominant systematic uncertainty in the broad-flux fraction. I would advise the editor that a sensitivity analysis over this fraction and a softening of the 'clear evidence' claim should be required before publication; these are fixable in revision and do not undermine the value of the MRS data set."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives the sharpest look yet at the inner 100–200 pc of Cen A with MIRI/MRS, and the central detection is real: broad line components with σ ∼ 600 km/s in [Ar II] and [Ne III], velocities reaching +1000/−1400 km/s in position–velocity diagrams, and a spatially resolved picture that ties the ionized gas to the jet direction. That is a solid, multi-line kinematic case for fast nuclear gas. The H2 analysis is careful too — the warped-disk versus radial-motion modeling is honest about degeneracies, and the absence of a fast warm-H2 outflow is a useful contrast with systems like IC 5063.\n\nThe soft spot is exactly where the reader's stress test lands: the mass outflow rate in Sect. 3.5 depends on the assumption that roughly half the Pfα flux belongs to the broad component, even though the line is too faint for its own kinematic decomposition. That fraction is borrowed from two-Gaussian fits to [Ar II] and [Ne III]. The problem is not just that it is an assumption — all such estimates have assumptions — but that the two-Gaussian and three-Gaussian fits give materially different broad-component widths and centroids, and recombination lines like Pfα do not have the same density dependence as collisionally excited lines. The quoted 1.6–2.9 M⊙/yr range only spans two velocity prescriptions, not the uncertainty in the flux fraction or the unresolved radius. That makes the headline number less firm than the text implies.\n\nA few lesser points: the sigma maps for H2 S(5) and [Ar II] use an artificially lowered instrumental resolution (footnote 1), which is a bit awkward but does not affect the main detection. And the summary's phrase \"clear evidence\" of an outflow outruns the unresolved, model-dependent decomposition — the broad component could in principle include an NLR or jet-cloud interaction. The authors do not fully rule those out. None of this kills the qualitative result, but it should temper the claim.\n\nThis is a single-object study, but it is the nearest radio galaxy and the MRS data are new and well reduced. The analysis is detailed, the modeling is standard, and the authors flag the main assumption themselves. I would send it to a serious referee, mainly to push for a quantitative treatment of the Pfα fraction and a slightly more cautious summary. For anyone working on jet-driven feedback in nearby AGN, it is worth reading and citing.","headline":"Genuine MIRI-MRS detection of fast ionized gas in Cen A's nucleus, but the mass outflow rate leans on an assumed Pfα broad-flux fraction that should be tested before the number becomes a benchmark.","tokens_in":26569,"tokens_out":2025,"would_cite":true,"duration_ms":23442,"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":"JWST/MIRI MRS observations of Centaurus A show a fast ionized gas outflow whose fastest gas, above 500 km/s, comes from the central 6 parsecs, with an ionized mass outflow rate near 1.6-2.9 solar masses per year and no fast molecular…","keywords":["Centaurus A","NGC 5128","radio galaxy","active galactic nuclei","ionized gas outflow","jet-ISM interaction","JWST MIRI MRS","mid-infrared spectroscopy"],"falsifier":"Resolve the hydrogen recombination line with an independent, higher-signal observation of the central 0.3 arcsec: if a kinematic decomposition shows that the broad Pf alpha fraction is far below half, or if the broad component is not seen in a recombination line, the outflow rate is not supported. In addition, a position-velocity map along the jet at sub-arcsecond resolution showing the +1000/-1400 km/s gas to be localized clouds rather than a coherent expanding shell would count against the jet-driven bubble interpretation.","tokens_in":25369,"feed_emoji":"🌌","tokens_out":6938,"duration_ms":73770,"temperature":0.7,"pith_summary":"This paper uses JWST/MIRI MRS mid-infrared spectroscopy to map the inner 100-200 parsecs of Centaurus A, the nearest radio galaxy, in both ionized and warm molecular gas. It claims that ionized gas is being expelled from the unresolved central 6 parsecs at speeds above 500 km/s, with line wings reaching roughly +1000 and -1400 km/s. The inferred ionized mass outflow rate is 1.6-2.9 solar masses per year, and the kinematics resemble simulation predictions of a jet-inflated bubble expanding into the galaxy's interstellar medium. In contrast, the warm molecular hydrogen shows rotation and noncircular motions but no fast outflow, suggesting that the jet couples only weakly to the molecular gas on these scales.","feed_headline":"Fast ionized outflow erupts from Cen A's central 6 parsecs","feed_subtitle":"MIRI spectroscopy maps jet-driven gas racing at up to 1,400 km/s, with no fast molecular counterpart","key_machinery":"The load-bearing observational object is the unresolved nuclear spectrum of the [Ar II] 6.99 micron and [Ne III] 15.56 micron lines, decomposed into one, two, and three Gaussian components; the broadest component, with sigma around 600 km/s and confined to a 6 pc aperture, is the outflow signature. Tilted-ring kinematic modeling separates the rotating disk component and exposes residual noncircular motions, while line-ratio diagnostics such as [Ne V]/[Ne II], [Ne III]/[Ne II], and [Fe II]/Pf alpha provide the shock-excitation evidence. The mass outflow rate is carried by the Pf alpha flux, the electron density from the two [Ne V] lines, and the standard expression dM/dt = 3 M v/R, evaluated for a compact unresolved source, which makes the rate a lower limit.","core_discovery":"The central claim is that Cen A's nuclear region contains a fast ionized gas outflow confined to the inner 6 pc, detected through broad components in the mid-infrared [Ar II] and [Ne III] line profiles with velocity dispersion around 600 km/s and extreme line-of-sight velocities near +1000 and -1400 km/s in position-velocity diagrams. Using the Pf alpha hydrogen recombination line, an electron density from the [Ne V] line ratio, and the assumption that about half of the Pf alpha flux is in the broad outflow component, the paper derives an ionized gas mass outflow rate of 1.6-2.9 solar masses per year. The paper also shows that both AGN radiation pressure and the radio jet power are energetically sufficient to launch the outflow, and that mid-infrared line ratios favor a contribution from shocks in addition to AGN photoionization. Unlike the ionized gas, the warm H2 lines show no clear fast nuclear outflow, only a weak blueshifted wing, which the authors attribute to low nuclear warm-H2 column density and the jet being launched nearly perpendicular to the circumnuclear disk.","pith_inferences":["A testable extension the paper leaves implicit: if the outflow is energy-driven by the jet, the momentum boost relative to AGN radiation (a factor of 16-64) should grow with distance from the nucleus, so mapping the outflow at 10-50 pc would distinguish a momentum-conserving from an energy-conserving bubble.","The same mid-infrared diagnostics could classify other optically obscured nearby AGNs: a compact broad component in low-ionization lines plus an elevated [Fe II]/Pf alpha ratio may be a generic jet-inflated bubble signature even when the narrow-line region is hidden by dust.","A deeper observation of the nuclear Pf alpha line, or of a brighter hydrogen recombination line at comparable spectral resolution, would directly test the assumed 50 percent broad-flux fraction and could revise the outflow rate by up to a factor of two.","Comparing the line shapes of low-ionization [Ar II] with high-ionization coronal lines at the same spectral resolution would settle whether the broad wings are truly outflowing gas or scattered AGN light, since scattering and photoionization would imprint different ionization-dependent profiles."],"forward_implications":["Cen A's nucleus is currently ejecting roughly 1.6-2.9 solar masses per year of ionized gas from a region only 6 pc across, a significant channel for removing gas from the very center of the galaxy.","The fastest ionized gas is confined to the unresolved nucleus, so the outflow is compact and the derived mass outflow rate is a lower limit that could rise with better spatial resolution.","The absence of a fast warm-H2 outflow implies that the jet's mechanical coupling to the molecular circumnuclear disk is weak on these scales, consistent with a jet launched almost perpendicular to the disk.","Both AGN radiation pressure and the radio jet can energetically drive the outflow, so energy arguments alone do not identify the driver; the jet is favored by its relativistic apparent speed.","Enhanced velocity dispersions perpendicular to the jet and the shock-sensitive line ratios align the observed kinematics with simulation predictions for a jet-inflated bubble expanding into the interstellar medium."],"supporting_citations":[{"why":"Supplies the adopted distance to Cen A (3.5 Mpc, 1 arcsec = 17 pc), the earlier near-infrared kinematics of the central arcsecond, and the warped-disk and jet backflow interpretations this work builds on.","marker":"Neumayer et al. 2007"},{"why":"Maps the cold molecular gas components including the circumnuclear disk, nuclear ring, nuclear disk, and streamers used for comparing H2 morphology and noncircular motions.","marker":"Espada et al. 2017"},{"why":"Provides the jet-inflated bubble and outflow simulations whose predicted kinematic signatures, including fast velocities and enhanced dispersions, match the observed features.","marker":"Mukherjee et al. 2016"},{"why":"Supplies the methodology and equation used to convert the Pf alpha flux and electron density into ionized gas mass and mass outflow rate.","marker":"Hermosa Muñoz et al. 2025"},{"why":"Provides the comparison radio galaxy IC 5063, showing broad ionized components and a jet-driven outflow, as well as the nuclear and extended line ratios used in the diagnostic diagram.","marker":"Dasyra et al. 2024"},{"why":"Established with Spitzer/IRS that high-ionization mid-infrared lines extend along the jet direction over tens of arcseconds, giving the larger ionized cone that the MRS field covers only in its inner part.","marker":"Quillen et al. 2008"},{"why":"Gives the jet position angle and apparent subluminal motion used to argue that the jet is likely the main influence on the nuclear scales.","marker":"Hardcastle et al. 2003"},{"why":"Provides an estimate of the radio jet power used in the energetics comparison showing the jet can launch the outflow.","marker":"Croston et al. 2009"},{"why":"Supplies the AGN plus fast-shock model tracks used to interpret the nuclear neon line ratios and support the shock-excitation conclusion.","marker":"Feltre et al. 2023"},{"why":"Concurrent work reported broad ionized lines in Cen A, cited alongside this paper as the first such detections.","marker":"Marconcini et al. 2025"}],"fun_headline_variants":["Cen A's ionized gas races at 1,400 km/s in nucleus","Fast ionized outflow confined to Cen A's inner 6 pc","JWST spots jet-driven outflow in Cen A's core","Ionized gas outflow in Cen A reaches 1,400 km/s","Cen A's fast outflow: ionized gas, not molecular"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the broad, faint component of the nuclear gas emission is truly gas moving outward in an organized outflow, and that roughly half of the hydrogen recombination flux sits in that component, rather than the broad wings being scattered light, gas lit up near the black hole, or a chance superposition of clouds.","fun_headline_variants_meta":{"raw":{"variants":["Cen A's ionized gas races at 1,400 km/s in nucleus","Fast ionized outflow confined to Cen A's inner 6 pc","JWST spots jet-driven outflow in Cen A's core","Ionized gas outflow in Cen A reaches 1,400 km/s","Cen A's fast outflow: ionized gas, not molecular"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001327,"raw_usage":{"total_tokens":5517,"prompt_tokens":1178,"completion_tokens":4339,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":794,"completion_tokens_details":{"reasoning_tokens":4244}},"tokens_in":794,"tokens_out":4339,"duration_ms":26826,"temperature":1.0,"reasoning_tokens":4244,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:38:07.764844+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Resolve the hydrogen recombination line with an independent, higher-signal observation of the central 0.3 arcsec: if a kinematic decomposition shows that the broad Pf alpha fraction is far below half, or if the broad component is not seen in a recombination line, the outflow rate is not supported. In addition, a position-velocity map along the jet at sub-arcsecond resolution showing the +1000/-1400 km/s gas to be localized clouds rather than a coherent expanding shell would count against the jet-driven bubble interpretation.","supporting_citations":[{"cited_title":"2007, , 671, 1329","cited_arxiv_id":null,"evidence_quote":"Supplies the adopted distance to Cen A (3.5 Mpc, 1 arcsec = 17 pc), the earlier near-infrared kinematics of the central arcsecond, and the warped-disk and jet backflow interpretations this work builds on."},{"cited_title":"E., et al","cited_arxiv_id":null,"evidence_quote":"Maps the cold molecular gas components including the circumnuclear disk, nuclear ring, nuclear disk, and streamers used for comparing H2 morphology and noncircular motions."},{"cited_title":"V., Sutherland , R., & Wagner , A","cited_arxiv_id":null,"evidence_quote":"Provides the jet-inflated bubble and outflow simulations whose predicted kinematic signatures, including fast velocities and enhanced dispersions, match the observed features."},{"cited_title":"M., Paraschos , G","cited_arxiv_id":null,"evidence_quote":"Provides the comparison radio galaxy IC 5063, showing broad ionized components and a jet-driven outflow, as well as the nuclear and extended line ratios used in the diagnostic diagram."},{"cited_title":"C., Bland-Hawthorn , J., Green , J","cited_arxiv_id":null,"evidence_quote":"Established with Spitzer/IRS that high-ionization mid-infrared lines extend along the jet direction over tens of arcseconds, giving the larger ionized cone that the MRS field covers only in its inner part."},{"cited_title":"J., Worrall , D","cited_arxiv_id":null,"evidence_quote":"Gives the jet position angle and apparent subluminal motion used to argue that the jet is likely the main influence on the nuclear scales."},{"cited_title":"Unveiling the Fast Acceleration of AGN-Driven Winds at Kiloparsec Scales","cited_arxiv_id":"2503.24359","evidence_quote":"Concurrent work reported broad ionized lines in Cen A, cited alongside this paper as the first such detections."}],"review_version":2}