{"id":"bf5e6e77-00e2-4c4d-a327-0c586b2302e8","arxiv_id":"2411.12398","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"JWST and ALMA reveal stratified molecular gas in MCG-05-23-16: cold CO and H2 trace a nuclear spiral and ring, while warmer H2 shows higher velocity dispersion in between, plus a possible localized star-formation-driven outflow.","lead":"Astronomers mapped warm and cold molecular gas in the Seyfert galaxy MCG-05-23-16 using JWST and ALMA. They find that cold gas follows a nuclear spiral and ring while warmer gas is more turbulent, including a clump that may be a small star-formation-driven outflow.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"R5 outflow interpretation rests on a single near/far-side geometry and an unpropagated 3DBAROLO model-parameter degeneracy; a cross-check of both would settle it.","rationale":"I agree with the reader that the morphological stratification is well-supported but the kinematical claims are model-dependent. The most load-bearing uncertainty is the near/far-side assignment and the model-parameter degeneracy in 3DBAROLO, exactly as identified in the reader's weakest_assumption. I recommend keeping the CONDITIONAL verdict: the outflow and inflow interpretations require systematic checks; the stratification conclusion does not. No ad hominem, and no manufactured concern: the concern is directly tied to Sect. 3.2.3, Fig. B.2/B.3, and Sect. 4. I have not found a separate critical flaw beyond this; the paper's internal reasoning is consistent, and the morphology claims are supported by the data.","tokens_in":27255,"tokens_out":1299,"duration_ms":11794,"concrete_test":"Re-fit the S(3) and CO(2-1) datacubes with 3DBAROLO (or a comparable tilted-ring code) while (a) fixing i=75 deg for S(3) instead of the fitted 37 deg, and (b) swapping the near/far-side geometry; if the R5 residual still appears as a spatially coherent redshifted knot in the blueshifted side in both cases, the outflow interpretation is robust, but if it vanishes under either change, the outflow claim becomes model-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the stratified cold/warm molecular gas picture, which is robustly supported by morphology alone. The load-bearing dynamical sub-claims are (1) inflow along the minor axis and (2) the R5 feature as an outflow, both derived from 3DBAROLO residuals in Sect. 3.2.3 and interpretation in Sect. 4. The weakest link is the assumed geometry: the near/far-side assignment (SE near, NW far, from extinction) converts identical residual patterns into inflow vs outflow. If the SE is actually far, the R5 sign and the minor-axis residuals flip meaning. Moreover, the 3DBAROLO fits return discordant inclinations for CO (i=75 deg) vs S(1) (i=33 deg) and S(3) (i=37 deg), with no systematic error budget or validation that the fitted tilted-ring model is unique. A warped disc or a two-component model could absorb the residuals now attributed to non-circular motions; the paper notes BAROLO cannot reproduce central dispersions but does not propagate this into the outflow claim. The R5 outflow (PA~270, velocities to ~350 km/s in PVD) is also based on a single pseudo-slit and residual-moment maps without an independent check against the full datacube. The morphological stratification survives independently; the outflow and inflow claims do not.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents JWST/MIRI MRS observations of the pure-rotational H2 lines S(1) through S(5) together with ALMA CO(2-1) observations of the Seyfert galaxy MCG-05-23-16, and combines them with archival HST and VLA images. The authors map fluxes, velocities, and velocity dispersions of the warm molecular gas, compare them with the CO morphology and kinematics, fit a single tilted-ring rotating-disc model with 3DBAROLO to CO(2-1), H2 S(1), and S(3), extract position-velocity diagrams along selected axes, and derive rotational-diagram temperatures, column densities, and masses for six regions. The central claim is a stratification of the molecular gas: the colder gas traced by CO(2-1), S(1), and S(2) lies in the nuclear spiral, ring, and connecting arms, whereas warmer gas traced by S(3) through S(5) fills the inter-arm regions and shows high velocity dispersions, with one knot (R5) at about 350 pc west of the nucleus interpreted as a possible star-formation-driven molecular outflow. The paper also reports non-circular motions along the kinematic minor axis, attributed to bar-driven inflow, and concludes that the compact VLA jet does not significantly affect the molecular gas.","tokens_in":27612,"tokens_out":5283,"duration_ms":46079,"significance":"If the claims hold, the paper provides one of the few spatially resolved comparisons of cold and warm molecular gas kinematics in a Seyfert galaxy, exploiting the complementary resolutions of JWST/MIRI and ALMA. The morphological stratification - cold gas in the nuclear spiral, ring, and connecting arms, and warmer, more turbulent gas in the inter-arm regions - is directly supported by the maps and the S(1)/S(3) excitation map, and it is a testable and potentially important result for barred galaxies and AGN fueling. The paper also demonstrates a useful methodological approach by combining 3DBAROLO residual analysis with PVDs for warm H2 lines, and the authors are appropriately cautious in presenting the outflow mass as an upper limit. However, the dynamical sub-claims (inflow and the R5 outflow) rest on assumptions about the disc model and the near/far-side geometry that are not yet cross-checked; these need to be strengthened before the kinematic interpretation can be considered secure.","major_comments":[{"comment":"The interpretation of the minor-axis residuals as inflow and of the R5 knot as outflow in Sect. 4 rests on the near/far-side assignment (southeast near, northwest far, based on dust extinction). Since the residual velocities are computed relative to a single 3DBAROLO rotating-disc model, flipping this geometric assignment would reverse the sign of the radial residual pattern and would change the R5 interpretation from outflow to a motion on the blueshifted side of the galaxy. Please test the alternative near/far-side geometry explicitly, for example by re-deriving the residual maps and PVDs under the opposite assignment or by using a kinematic model with the two possible sides, and show how the inflow and outflow claims fare in each case.","section":"Sect. 4, near/far-side geometry"},{"comment":"The 3DBAROLO fits return very different inclinations for CO (i=75 degrees) and for H2 S(1) (i=33 degrees) and S(3) (i=37 degrees), which the paper attributes to the smaller MIRI FOV. Because the residuals that define the non-circular motions are computed with respect to these fitted models, the amplitude and sign of the residual velocities are sensitive to the adopted inclination and position angle, which are not given error bars or systematically explored. Please provide a quantitative exploration of the allowed range of (i, PA) for each tracer (for example, fixed versus free fits or a grid of models) and show how the residual maps and the R5 PVD change within that range.","section":"Sect. 3.2.3"},{"comment":"The R5 outflow claim is based on a single pseudo-slit PVD along PA=87 degrees and on the residual moment maps of S(3). The paper notes that BAROLO cannot reproduce the high central velocity dispersions but does not propagate this limitation into the outflow interpretation, and no independent check on the full datacube is presented. Please verify the R5 red wing with a two-Gaussian decomposition of the S(3) line in that aperture and/or with channel maps covering the velocities around +350 km/s, so that the non-circular component is confirmed independently of the model subtraction.","section":"Sect. 3.2.3, Fig. 7"},{"comment":"In Sect. 4 the mass outflow rate is computed as Mdot_out = M_R5 times v_out divided by r_out, using the total mass of R5 (about 4000 solar masses), a representative velocity of about 200 km/s, and r=88 pc, with the caveat that only part of the gas is actually outflowing. As stated, this makes the upper limit (<0.01 solar masses per year) not directly comparable to the SFR of 0.013 solar masses per year, since the numerator and denominator refer to different gas components. A two-component line fit that isolates the outflowing part of R5 would provide a self-consistent mass-outflow rate and a firmer basis for the mass-loading comparison.","section":"Sect. 4, paragraph beginning 'The mass of the H2 knot'"}],"minor_comments":[{"comment":"The paper switches between '3DBAROLO' and 'BAROLO' (for example in Sect. 3.2.3 and in Fig. B.3); please use the full name consistently or define the abbreviation once.","section":"Throughout"},{"comment":"In the caption of Fig. 3, the beam sizes are quoted as '0.14 x 0.30 and 0.68 x 0.83' without units; please add arcseconds.","section":"Fig. 3 caption"},{"comment":"The statement that 'if we subtract the S(5) from the S(1) velocity map' is not illustrated; consider adding a residual map to make this comparison quantitative.","section":"Sect. 3.2.1"},{"comment":"The central positions used for the 3DBAROLO fits are given in the text for S(1) and S(3), but the adopted coordinate for the CO fit is only described as the peak of the 200 GHz continuum; please list the actual coordinates for all three tracers for reproducibility.","section":"Sect. 3.2.3"},{"comment":"In Table 2, several sigma entries are quoted as '< 41', '< 84', '< 29', and '< 131' km/s, but the table caption does not explain whether these are upper limits from non-detections or from the fitting; please clarify the notation.","section":"Table 2 caption"},{"comment":"The abstract carefully says the R5 kinematics 'are consistent with outflowing gas', while the conclusions repeat the claim without the 'consistent with' caveat; consider harmonizing the strength of the wording.","section":"Abstract and Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is generally well presented and the data are valuable. The main concern is that the dynamical sub-claims, particularly the minor-axis inflow and the R5 outflow, are presented with more confidence than the current model robustness warrants. The morphological stratification conclusion is strong and will likely survive the requested tests, so I recommend major revision rather than rejection. Please ensure the requested geometry tests and model-degeneracy exploration are performed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the paper on MCG-05-23-16. The main claim holds: cold molecular gas (CO(2-1), H2 S(1)/S(2)) traces the nuclear spiral, ring, and arms, while warmer gas (S(3)-S(5)) fills the inter-arm regions with high velocity dispersion. That stratification is directly visible in the maps and does not depend on any model. The paper is the first joint JWST/MIRI + ALMA kinematic study of this galaxy, and it does a clean job with standard reduction and line fitting. The rotational diagrams and two-temperature fits are careful, and the R5 knot at ~350 pc west, with its S(3) dispersion ~160 km/s, PAH 11.3 micron emission, and temperature ~260 K, is a genuinely interesting candidate for a localized, star-formation-driven molecular outflow. The conclusion that the compact VLA jet does not significantly disturb the molecular gas is argued from spatial anti-correlation; it is reasonable but weaker.\n\nWhere the paper gets soft is the dynamical interpretation. The minor-axis inflow and the R5 outflow are read from 3DBAROLO residuals, and that reading depends on two choices the paper does not fully stress-test. First, the near/far side is assumed (SE near, NW far) from extinction; flipping it flips the sign of the radial motions and the outflow direction. Second, the 3DBAROLO fits give very different inclinations for CO (75 deg) vs S(1) (33 deg) and S(3) (37 deg), with no error budget or uniqueness check; the paper's explanation (smaller MIRI FOV) is plausible but not demonstrated. The outflow also rests on a single pseudo-slit PVD, not an independent datacube cross-check. The authors do note that BAROLO cannot reproduce the central dispersions, but they don't propagate that limitation into the outflow mass/rate. The good news is the paper is appropriately hedged — 'consistent with' and 'possibly' — and the mass outflow rate is presented as an upper limit with mass loading <0.8. So the soft spot is real but not enough to sink the main result.\n\nThe citation pattern is clean; the GATOS reduction references are the actual sources for the methods. I'd send this to a referee. The observational material is valuable, the stratification is robust, and the outflow candidate is worth pursuing with a more careful dynamical model (e.g., two-component fits, warped disc, or a full datacube residual analysis). If I were the editor, I'd ask for the geometry/disc-model robustness check before accepting the outflow as more than a candidate. For a reading group, it's a decent example of the JWST+ALMA synergy but not a methods landmark.","headline":"A solid JWST+ALMA single-galaxy study whose stratification result is robust, but the inflow/outflow sub-claims rest on one 3DBAROLO geometry and need a cross-check before being quoted.","tokens_in":28286,"tokens_out":4026,"would_cite":false,"duration_ms":38591,"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":"Combining JWST mid-infrared H$_2$ lines with ALMA CO(2-1), the paper shows that molecular gas in the Seyfert galaxy MCG-05-23-16 is stratified: cold gas rotates in a bar-driven spiral and ring while warmer, turbulent gas fills the…","keywords":["molecular gas stratification","warm molecular gas","H2 rotational lines","Seyfert galaxy","AGN fueling","bar-driven inflow","star-formation-driven outflow","CO(2-1) kinematics"],"falsifier":"Observe the R5 region at higher spatial and spectral resolution in H$_2$ S(3) and in a cold-gas tracer such as CO: if the positive-velocity component disappears once a warped or two-component disk model replaces the single tilted ring, or if CO is detected there, the outflow interpretation is refuted; detecting whether the 11.3 micron PAH emission is really co-spatial with the outflowing gas would test the star-formation driver.","tokens_in":27029,"feed_emoji":"🌀","tokens_out":11202,"duration_ms":101653,"temperature":0.7,"pith_summary":"The paper combines JWST mid-infrared spectroscopy with ALMA millimeter-wave imaging of the Seyfert galaxy MCG-05-23-16 to show that its molecular gas is stratified by temperature. The cold gas traced by CO(2-1) and the lower pure-rotational H$_2$ lines sits in a nuclear trailing spiral, a surrounding ring, and connecting arms, while warmer and more turbulent gas traced by H$_2$ S(3), S(4), and S(5) fills the spaces between those structures. The authors interpret the cold-gas morphology and kinematics, including S-shaped zero-velocity contours and non-circular residuals along the kinematic minor axis, as the signature of a nuclear bar funneling gas inward. One warm-gas clump roughly 350 pc west of the nucleus shows velocities that do not fit the rotating-disk model and is read as a localized molecular outflow, probably driven by recent star formation traced by 11.3 micron PAH emission rather than by the AGN's compact radio jet. If the interpretation is right, the galaxy is simultaneously a clear case of bar-driven AGN fueling and a case where a compact jet passes through a molecular disk without strongly disrupting it.","feed_headline":"Stratified gas reveals how a Seyfert galaxy feeds its black hole","feed_subtitle":"JWST and ALMA show cool gas funneling inward as warmer, turbulent gas fills the spaces between spiral arms.","key_machinery":"The argument is carried by comparing molecular gas tracers at different temperatures: CO(2-1) with the lower H$_2$ pure-rotational lines S(1) and S(2) for cold gas at tens of Kelvin, and S(3), S(4), and S(5) for warmer gas at hundreds of Kelvin. The quantitative tool is 3DBAROLO, a named 3D tilted-ring model that fits a rotating disk to each datacube; subtracting the model residual identifies gas that is not in circular rotation. The 'forbidden velocity' criterion, gas appearing redshifted on the blueshifted side of the rotation pattern, is what converts the R5 residual into an outflow candidate. Supporting measurements are the S(1)/S(3) excitation map, rotational-diagram fits that separate warm and hot components and give their temperatures and masses, and the 11.3 micron PAH feature used as a star-formation tracer.","core_discovery":"The central claim is a temperature stratification of the molecular interstellar medium in MCG-05-23-16, seen by mapping five H$_2$ pure-rotational lines together with CO(2-1). The cold gas forms a ~350 pc nuclear trailing spiral connected to a ~1.4 kpc ring, with roughly regular rotation reaching $\\pm 250$ km s$^{-1}$ and residual velocities along the minor axis that the authors attribute to inflow on elliptical orbits in a barred potential. The warmer H$_2$ transitions show more disturbed kinematics, including clumps with velocity dispersion up to about 160 km s$^{-1}$, located precisely where CO(2-1) is absent. One such clump, R5, about 350 pc west of the nucleus, has redshifted velocities on the blueshifted side of the galaxy, a 'forbidden' position in the rotation model, and the paper interprets it as outflowing warm gas most plausibly driven by localized star formation rather than by the ~200 pc VLA jet, which shows no spatial association with the molecular gas. The paper concludes that cold molecular gas rotates and flows inward while warmer, higher-dispersion molecular gas occupies the inter-arm volume, and that the AGN jet does not significantly affect the molecular gas disk.","pith_inferences":["If this stratification is generic, other barred Seyferts observed with JWST and ALMA should show high-dispersion warm H$_2$ filling the gaps between CO spiral arms and rings; that is a testable prediction for a small survey.","The R5 outflow may belong to a broader class of star-formation-driven molecular outflows located near the ends of nuclear spirals or bars, well away from the jet axis; similar knots should appear in other galaxies where PAH emission coincides with high-dispersion H$_2$.","The apparent non-interaction between the jet and the molecular gas could be a resolution or projection effect; deeper ALMA observations at sub-100 pc scales might reveal jet-induced turbulence that the current 70 pc beam washes out.","Because S(5) was outside the field at R5, the clump's temperature rests on a single-component fit and is uncertain; obtaining S(5) coverage there would test whether the 'outflowing' gas is truly hot and excited or merely blended kinematic components."],"forward_implications":["Molecular gas is being delivered to the AGN through a bar-driven trailing spiral and ring, so AGN fueling in this S0 galaxy appears to work by secular inflow rather than by a major merger.","Warm and cold molecular gas are spatially anti-correlated in the central kiloparsec, so studies that trace only CO can miss the turbulent warm phase that fills the inter-arm regions.","If the R5 knot is an outflow, star formation rather than the AGN jet can drive warm molecular outflows in a Seyfert galaxy, with a mass loading factor below 0.8 relative to the local star formation rate.","The warm molecular gas mass is less than 5 percent of the cold gas mass in the central 630 pc under the adopted Milky-Way conversion factor, so the warm phase is a small but kinematically distinct reservoir.","The compact jet's lack of a detectable effect on the molecular gas implies that jet-disk orientation is a decisive factor for whether AGN feedback disturbs the cold ISM."],"supporting_citations":[{"why":"Supplies the 3DBAROLO 3D tilted-ring code used to fit the rotating disk and define non-circular residuals for CO(2-1), S(1), and S(3).","marker":"Di Teodoro & Fraternali 2015"},{"why":"Provides the VLA 8.4 GHz detection of the compact ~200 pc jet whose spatial association with the gas the paper tests.","marker":"Orienti & Prieto 2010"},{"why":"Reports the dust-lane morphology and the A_V=1 extinction used for line corrections and for the near-side/far-side geometry.","marker":"Prieto et al. 2014"},{"why":"Gives the previous optical IFU kinematics (major-axis PA and absence of clear ionized outflows) as the baseline the molecular kinematics are compared with.","marker":"Ruschel-Dutra et al. 2021"},{"why":"Calibration used to convert the 11.3 micron PAH luminosity into the star formation rate that supports the star-formation-driven outflow interpretation.","marker":"Shipley et al. 2016"},{"why":"Provides the equations used to estimate the outflow mass rate and kinetic energy upper limit for R5.","marker":"Speranza et al. 2024"},{"why":"Supplies the Milky-Way CO-to-H$_2$ conversion factor used to compute the cold molecular gas mass.","marker":"Bolatto et al. 2013"},{"why":"Describes the rotational-diagram method and provides comparison warm and hot gas masses for infrared-bright Seyferts.","marker":"Pereira-Santaella et al. 2014"}],"fun_headline_variants":["Cold gas spirals in, warm gas churns out in Seyfert galaxy","JWST and ALMA reveal molecular gas strata in AGN","Warm gas turbulence and cold gas inflow mapped in MCG-05-23-16","Seyfert's stratified gas: cool disk, turbulent warm inter-arm","Molecular gas stratification exposes Seyfert feeding and outflow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All the dynamical interpretations assume that one 3DBAROLO tilted-ring rotating disk is the correct baseline for the CO and H$_2$ kinematics, and that the southeast side of the galaxy is the near side because it is more extincted; if those assumptions fail, the inferred inflows and the R5 outflow weaken, although the observed temperature stratification would survive.","fun_headline_variants_meta":{"raw":{"variants":["Cold gas spirals in, warm gas churns out in Seyfert galaxy","JWST and ALMA reveal molecular gas strata in AGN","Warm gas turbulence and cold gas inflow mapped in MCG-05-23-16","Seyfert's stratified gas: cool disk, turbulent warm inter-arm","Molecular gas stratification exposes Seyfert feeding and outflow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000807,"raw_usage":{"total_tokens":3682,"prompt_tokens":1220,"completion_tokens":2462,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":836,"completion_tokens_details":{"reasoning_tokens":2368}},"tokens_in":836,"tokens_out":2462,"duration_ms":17834,"temperature":1.0,"reasoning_tokens":2368,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:34:47.255731+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the R5 region at higher spatial and spectral resolution in H$_2$ S(3) and in a cold-gas tracer such as CO: if the positive-velocity component disappears once a warped or two-component disk model replaces the single tilted ring, or if CO is detected there, the outflow interpretation is refuted; detecting whether the 11.3 micron PAH emission is really co-spatial with the outflowing gas would test the star-formation driver.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the 3DBAROLO 3D tilted-ring code used to fit the rotating disk and define non-circular residuals for CO(2-1), S(1), and S(3)."},{"cited_title":"2021, MN- RAS, 507,","cited_arxiv_id":null,"evidence_quote":"Gives the previous optical IFU kinematics (major-axis PA and absence of clear ionized outflows) as the baseline the molecular kinematics are compared with."},{"cited_title":"V ., Papovich, C., Rieke, G","cited_arxiv_id":null,"evidence_quote":"Calibration used to convert the 11.3 micron PAH luminosity into the star formation rate that supports the star-formation-driven outflow interpretation."},{"cited_title":"A., et al","cited_arxiv_id":null,"evidence_quote":"Provides the equations used to estimate the outflow mass rate and kinetic energy upper limit for R5."},{"cited_title":"D., Wolfire, M., & Leroy, A","cited_arxiv_id":null,"evidence_quote":"Supplies the Milky-Way CO-to-H$_2$ conversion factor used to compute the cold molecular gas mass."},{"cited_title":"P., et al","cited_arxiv_id":null,"evidence_quote":"Describes the rotational-diagram method and provides comparison warm and hot gas masses for infrared-bright Seyferts."}],"review_version":1}