{"id":"a1002128-e04b-4330-9dae-508bb91da3e5","arxiv_id":"2505.08876","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The CO 2-1/1-0 ratio across M83 varies systematically with bar, spiral arms, and interarm regions, indicating that galactic dynamics rather than stellar feedback primarily drives the physical state of molecular gas.","lead":"Using ALMA data, this paper maps the ratio of two carbon monoxide emission lines across the barred spiral galaxy M83 at 46 parsec resolution, showing that the ratio rises in spiral arms and bars and falls between them. The pattern suggests molecular gas is heated or compressed by galactic structure and rotation more than by star formation feedback.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Beam-filling by diffuse interarm CO(1-0) is not quantified; if interarm 46-pc beams mix compact dense gas with diffuse low-excitation gas, the R21 arm/interarm contrast and the inferred 2-3x density/temperature change could be artifacts.","rationale":"The reader and I converge on the same weakest point. The paper's calibration work, deep CLEANing, inclusion of total-power data, and the R21 <= 1 histogram are real strengths; the large-scale pattern is likely real. However, converting that pattern into a claim about the physical conditions of the bulk molecular gas inside clouds requires the beam to isolate clouds, and Section 4's assertion is not backed by a quantitative isolation statistic. The diffuse extended CO(1-0) seen in Figure 1 is direct evidence that beam mixing is plausible. The proposed cloud-catalog test would settle it with already-published data. Other concerns, such as the snapshot-to-evolution inference and HII-mask completeness, are secondary: the map itself is the strong product, but the causal interpretation remains conditional. Hence the reader's CONDITIONAL verdict and moderate confidence remain appropriate, and no adjustment to the verdict is needed.","tokens_in":25920,"tokens_out":7323,"duration_ms":84248,"concrete_test":"For every 46-pc beam, use the CO(1-0) cloud catalog (Hirota et al. 2024) to compute the fraction of CO(1-0) flux in that beam associated with the dominant cataloged cloud. Repeat the radial and azimuthal R21 analysis for beams with this fraction larger than 0.8 (isolated) and smaller than 0.5 (mixed). If the arm-interarm R21 contrast persists in the isolated-beam subset with similar amplitude, the beam-filling objection is refuted. If the contrast disappears or weakens substantially, the low interarm R21 is dominated by diffuse gas, and the claim of a factor 2-3 change in cloud physical conditions would need to be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is stated, not established, in Section 4: 'The resolution of 46 pc is close to the typical cloud diameter of 40 pc... When the clouds are isolated in each beam, the line ratio does not suffer from the beam filling factor.' The paper never quantifies the fraction of beams with isolated clouds. Section 4.1 explicitly notes that CO(1-0) is more spatially extended than CO(2-1), especially toward interarm regions. That is exactly the situation in which the 46-pc beam contains diffuse, low-excitation CO(1-0) that contributes little to CO(2-1), plus compact denser clumps. Then R21_obs ~ f_dense * R21_cloud, where f_dense = I10,dense/I10,total. A change in f_dense from ~0.5 (interarm) to ~1 (arms) gives the observed jump from ~0.5-0.6 to ~0.7-0.8 with no change in cloud density or temperature. Figure 4's R21(T) versus R21(I) comparison does not remove this degeneracy because both quantities use the same spatial beam and the same diffuse component. Without an isolation or diffuse-flux quantification, the central inference that the bulk molecular gas within clouds changes its n_H2/T_k by a factor of 2-3, and hence the rotation-synchronized evolutionary interpretation, is not uniquely supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a full-disk map of the CO J=2-1/1-0 line ratio (R21/10) in the barred spiral galaxy M83, constructed from ALMA 12m+7m+Total Power observations that are jointly imaged to a common 46 pc resolution. The authors report systematic large-scale variations: low R21/10 (≲0.7) in interarm regions and higher values (≳0.7) in the bar and spiral arms, with even higher values (0.8-1.0) localized around HII regions, and a radial sequence of high values at the center, lower values in the bar, a rise near the bar ends, and a gradual decline outward. They interpret the pattern as evidence that the density and/or temperature of the bulk molecular gas changes by a factor of ~2-3 in response to galactic structures and dynamics rather than stellar feedback, and they argue that this implies molecular gas and clouds survive and evolve over at least a substantial fraction of the galactic rotation timescale. The paper compares the results with Milky Way HRG/LRG classifications and with other barred galaxies, and discusses the limited, localized impact of star formation feedback.","tokens_in":26188,"tokens_out":6812,"duration_ms":67909,"significance":"If the central interpretation holds, this is a valuable and possibly general result: high-resolution, full-disk R21 imaging of a nearby barred spiral can directly trace where and when molecular gas physical conditions change, and the data appear to support dynamically driven evolution over feedback-dominated evolution. The observational product is strong: the paper combines ALMA 12m, 7m, and Total Power data, carefully matches beams to 46 pc, discusses calibration and sensitivity cuts transparently, and provides a detailed appendix on imaging and HII-region masking. The comparison with the Milky Way and other barred galaxies is constructive, and the authors are appropriately cautious in several places (the causal role of feedback, the crudeness of the LVG conversion, the sensitivity bias). The main risk is that the quantitative claim of a factor-of-2-3 change in gas density/temperature rests on untested beam-filling and on external, model-dependent calibrations, so the significance of the physical conclusion is somewhat uncertain until those issues are quantified.","major_comments":[{"comment":"The central interpretation assumes that the 46 pc beam is mostly filled by isolated molecular clouds, so that the observed R21 variation reflects changes in gas density and temperature rather than beam filling. The manuscript states that 'when the clouds are isolated in each beam, the line ratio does not suffer from the beam filling factor' but does not quantify the fraction of beams that satisfy this condition, nor the contribution of diffuse, low-excitation CO(1-0) emission to the interarm flux. Section 4.1 explicitly notes that CO(1-0) is more spatially extended than CO(2-1), especially toward interarm regions, which is exactly the situation in which the low interarm R21 could result from a larger diffuse-to-clump flux ratio rather than from lower n_H2/T_k in the bulk gas. The comparison in Figure 4 between R21/10(I) and R21/10(T) does not remove this degeneracy because both quantities use the same beam and the same diffuse component. This is load-bearing for the abstract and Section 5.4 claim that the bulk molecular gas changes its density/temperature by a factor of 2-3; please quantify the diffuse flux fraction and the beam-isolation fraction, for example by comparing matched-beam maps at coarser resolution or by modeling a two-component (diffuse plus compact) beam-filling decomposition.","section":"Section 4 (paragraph beginning 'The R21/10 values reported here...') and Section 5.4"},{"comment":"The quantitative statement that R21 increases by about a factor of 2-3 in density and/or kinetic temperature between interarm and arm/bar regions is imported from external LVG calculations and from Milky Way LRG/HRG/VHRG classifications rather than derived from an inversion of the M83 data. R21 alone has well-known degeneracies among n_H2, T_k, N_CO, line opacity, and beam filling, and Section 3.1 describes the HRG/LRG boundary as a 'crude guideline.' The abstract and conclusions nevertheless present the factor of 2-3 as a quantitative result. Please either provide an explicit LVG-based inversion with assumed M83 parameters and show how the inferred n_H2/T_k ranges depend on those assumptions, or systematically soften the quantitative claim to state that the observations are consistent with such changes under the adopted conversion. Without this, the strongest quantitative conclusion in the abstract is not uniquely supported by the data.","section":"Section 3.1 and Section 5.4"},{"comment":"The inference that the large-scale R21 pattern demonstrates evolutionary synchronization with galactic rotation is based on a single snapshot. The observed azimuthal sequence (interarm low, arm high, interarm low again) could in principle be maintained by a quasi-static population of molecular clouds whose physical conditions depend on local environment, rather than by individual clouds that travel across the disk and evolve along the flow. The argument that short cloud lifetimes would destroy the large-scale pattern assumes that clouds must be advected between structures; a stationary pattern could also be produced by short-lived clouds forming in situ with different properties in different environments. If the authors wish to claim that the evolution is governed by galactic rotation timescales, they should provide additional support, for example a kinematic test, an azimuthal-phase-resolved analysis that accounts for the pattern speed and possible corotation, or an explicit statistical model showing that a transient-cloud population cannot reproduce the observed azimuthal coherence.","section":"Section 5.1 and Section 5.4"}],"minor_comments":[{"comment":"The notation is inconsistent: the paper mostly uses R21/10 but Table 1 and some figure captions use R21; please unify the notation.","section":"Throughout"},{"comment":"There is a typo in 'sidebloes' (should be 'sidelobes').","section":"Appendix A.1"},{"comment":"The CO(2-1) field is smaller than the CO(1-0) field, and the TP filling uses CO(1-0) scaled by 0.55 outside the observed CO(2-1) field. The paper states that this transition zone lies outside the final field of view, but it would be helpful to state explicitly in Section 2.1 whether any pixels in the R21 map near the outer disk and interarm regions could be affected by this procedure.","section":"Section 2.1 and Appendix A.2.2"},{"comment":"The statement that the radii of R21/10 > 0.8 regions are ≲100 pc around prominent HII regions is based on visual inspection; please either quantify this (e.g., with radial profiles or median sizes) or explicitly label it as a qualitative impression.","section":"Section 4.4"},{"comment":"The value of trot,p = 6 Gyr at R = 3 kpc is very sensitive to the adopted pattern speed because Ω and Ωp are nearly equal near corotation; it would be helpful to state this sensitivity and to show how the conclusion would change for a small variation in Ωp.","section":"Section 5.1"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a genuinely impressive observational dataset and a clear, well-illustrated R21 map. The main issue is interpretive: the quantitative physical conclusions depend on beam-filling and on external LVG classifications, and the dynamical-timescale argument is somewhat stronger than a snapshot can support. I would like the revision to focus on quantifying the beam-filling/diffuse-flux contribution and on properly caveating or reframing the factor-of-2-3 claim; with those changes the paper could be a strong contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The full-disk R21 map at 46 pc is the real contribution here, and it is good. The imaging is careful, the calibration discussion is honest, and the large-scale pattern is visually and statistically clear: low ratios in interarms, high in arms and bar, high in the center, low again through the bar region. That pattern is not something I doubt. The paper also does a decent job showing the high-ratio gas is localized near HII regions, and the low-ratio gas dominates both area and a large share of the flux.\n\nThe soft spot is exactly where the stress-test note points. The paper asserts that the 46 pc beam isolates molecular clouds and therefore the line ratio is not beam-filling limited, but it never quantifies how many beams actually contain isolated clouds. CO(1-0) is more extended than CO(2-1), especially in the interarm regions. That is precisely the situation where the observed arm/interarm R21 contrast could be driven by changing fraction of diffuse low-excitation gas in the beam rather than by a 2-3x change in the density and temperature of the gas within clouds. The R21(T) versus R21(I) comparison does not break this degeneracy. So the claim that the bulk molecular gas within clouds is evolving dynamically is plausible but not uniquely supported by the data as presented.\n\nThat said, the concern is not fatal to the paper as a whole. The map is real, the large-scale pattern is real, and the feedback-localization claim (R21 above 0.8 mostly within ~100 pc of HII regions) is less sensitive to the beam-filling issue because high ratios require genuinely excited gas. The main overreach is in the section that converts R21 variations directly into a factor-of-2-3 density/temperature change and uses that to argue against short cloud lifetimes and feedback-driven evolution. The paper itself is honest about some limitations, including the causality issue around HII regions and the cursory arm/interarm separation, but the beam-isolation step is treated as established rather than tested.\n\nFor readers in molecular gas and galaxy evolution, this is worth engaging with seriously. It deserves peer review. A good referee should ask for a quantification of the isolated-cloud fraction, an accounting of diffuse CO(1-0) flux, and ideally a cloud-by-cloud R21 analysis using the existing catalog. With those additions, the dynamical-evolution argument would be much stronger. Without them, the paper still stands as a valuable data product and a useful observationally grounded challenge to the feedback-dominated picture, just with a weaker interpretation than the abstract implies.","headline":"A valuable, carefully made 46 pc R21 map of M83 with a robust large-scale pattern, but the paper's central dynamical-evolution interpretation leans on an unquantified beam-filling assumption.","tokens_in":26763,"tokens_out":2524,"would_cite":true,"duration_ms":27970,"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":"Molecular gas in M83 changes state with galactic rotation, not with star formation feedback.","keywords":["CO(2-1)/CO(1-0) line ratio","molecular gas","barred spiral galaxies","galactic dynamics","stellar feedback","M83","ALMA","interstellar medium"],"falsifier":"Reconstruct the M83 CO maps at about 10 pc resolution in one interarm and one arm field: if the interarm CO(1-0) breaks into diffuse low-level gas that CO(2-1) misses while the arm emission stays in compact bright clumps, then the low interarm R21 is in part a beam-filling artifact rather than a true change in the bulk gas state.","tokens_in":25722,"feed_emoji":"🌌","tokens_out":6217,"duration_ms":58461,"temperature":0.7,"pith_summary":"This paper maps the CO J=2-1/1-0 line ratio (R21) across the barred spiral galaxy M83 at 46 parsec resolution, resolving individual molecular clouds. The map shows that R21 is low (<0.7) in interarm regions and high (>0.7) in the bar and spiral arms, implying that the gas density and kinetic temperature rise by a factor of roughly 2-3 as gas flows into the arms. This rise occurs even in arm regions with no star formation, while star-forming HII regions only push R21 to 0.8-1.0 within about 100 pc. The authors conclude that galactic structures and rotation timescales, rather than stellar feedback, govern the evolution of the bulk molecular gas, and that molecular clouds survive for at least ~100 Myr.","feed_headline":"Galactic rotation, not star feedback, drives M83's gas evolution","feed_subtitle":"A 46-pc CO line-ratio map shows gas density and temperature rising with spiral-arm flow, while star-formation effects stay within ~100 pc.","key_machinery":"The central tool is the CO J=2-1/1-0 brightness ratio R21, a diagnostic of molecular gas excitation: higher R21 indicates higher H2 volume density and/or kinetic temperature, calibrated by Large Velocity Gradient modeling that ties R21~0.6-0.7 to typical clouds of density ~300 $cm^{-3}$ and temperature ~10 K. The analysis relies on combining interferometer and total-power data to reach 46 pc resolution, close to the ~40 pc typical cloud diameter, and adopts the Milky Way classification of Low Ratio Gas (<0.7), High Ratio Gas (0.7-1.0), and Very High Ratio Gas (>1.0). Multi-wavelength comparisons (H-$\\alpha$, infrared, and HI images) separate regions with and without star formation to isolate dynamical compression from feedback heating.","core_discovery":"Using high-resolution CO(1-0) and CO(2-1) observations of M83 with a common 46 pc round beam, this paper constructs the line ratio R21 = I21/I10 and finds systematic large-scale variations tied to galactic structures. Azimuthally, R21 rises from below 0.7 in the interarm gas to at least 0.7 in the bar and spiral arms; radially, it is very high (~1.0) in the central 500 pc, drops to about 0.6 through the bar, rises again around the bar end, and declines toward the outer disk. Since R21 is sensitive to H2 density and kinetic temperature, this pattern evidences that the bulk molecular gas is compressed and heated by about a factor of 2-3 as it rotates through spiral-arm shocks, and that this evolution is synchronized with galactic rotation rather than triggered by star formation. The influence of stellar feedback appears localized to within about 100 pc of HII regions, and massive, non-star-forming molecular structures in interarm regions show low R21. The paper argues that molecular gas/cloud evolution in barred spirals is dynamically driven, with cloud lifetimes lasting a substantial fraction of the rotation timescale.","pith_inferences":["If the pattern holds, star formation prescriptions in galaxy simulations should couple molecular gas state to large-scale kinematics (arm crossings) rather than only to local feedback.","A testable extension is to compare R21 at matched physical resolution in a strongly flocculent or non-barred spiral; a weaker arm/interarm contrast would confirm that the bar and density waves are the driver rather than rotation alone.","The authors' assumption of resolved clouds could be checked statistically: the fraction of 46 pc beams containing a single velocity-coherent cloud can be measured from the CO(1-0) data cube, and R21 recomputed after excluding blended beams."],"forward_implications":["Molecular gas and clouds in barred spirals evolve on galactic rotation timescales (order 100 Myr or more), not the short 5-30 Myr cloud lifetimes proposed in recent literature.","Stellar feedback modifies molecular gas conditions only in localized patches within about 100 pc; feedback is not the main driver of the bulk gas evolution across a disk.","CO(2-1) cannot be treated as an equivalent mass tracer to CO(1-0) when galactic structures are resolved, since the ratio varies systematically with environment.","The observed interarm-to-arm R21 contrast provides a direct, resolved test of cloud assembly and disruption along galactic flows.","Similar R21 patterns in other barred spirals indicate that dynamically driven gas evolution is general among local barred spiral galaxies."],"supporting_citations":[{"why":"Provides the LRG/HRG classification and Milky Way arm/interarm R21 variations that set the 0.7 dividing line and the interpretation of galactic-structure-driven gas evolution.","marker":"Sakamoto et al. 1997b"},{"why":"Establishes the LVG radiative-transfer framework that connects the line ratio to gas density and kinetic temperature.","marker":"Goldreich & Kwan 1974"},{"why":"Supplies calibrated external-galaxy R21 mapping and LVG models showing that a factor 2-3 change in density/temperature spans the low-to-high ratio transition.","marker":"Koda et al. 2012"},{"why":"Provides the M83 CO(1-0) data, structure definitions (center, bar, disk), rotation curve, and mass measurements that this paper builds on for imaging and interpretation.","marker":"Koda et al. 2023"},{"why":"Shows a similar radial R21 pattern in the barred spiral NGC 1300, supporting the generality of the dynamically driven evolution scenario.","marker":"Maeda et al. 2022"},{"why":"Delivers Milky Way galactic center and disk R21 measurements that define the expected HRG/VHRG behavior and calibration reference.","marker":"Sawada et al. 2001"},{"why":"Provides an ALMA R21 map of NGC 3627 with similar variations, and raises the sensitivity-bias caveat that this paper discusses.","marker":"den Brok et al. 2023"}],"fun_headline_variants":["Spiral flow, not star birth, governs M83's gas state","Galactic rotation sets the pace of M83's gas evolution","M83's gas state follows spiral structure, not star feedback","Rotation, not star formation, drives M83's gas heating"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the 46 pc beam mostly isolates individual molecular clouds, so that the observed R21 variations reflect genuine changes in gas density and temperature rather than beam filling by diffuse, low-excitation gas.","fun_headline_variants_meta":{"raw":{"variants":["Spiral flow, not star birth, governs M83's gas state","Galactic rotation sets the pace of M83's gas evolution","M83's gas state follows spiral structure, not star feedback","Rotation, not star formation, drives M83's gas heating"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000561,"raw_usage":{"total_tokens":2760,"prompt_tokens":1133,"completion_tokens":1627,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":749,"completion_tokens_details":{"reasoning_tokens":1554}},"tokens_in":749,"tokens_out":1627,"duration_ms":10913,"temperature":1.0,"reasoning_tokens":1554,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:46:11.124418+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Reconstruct the M83 CO maps at about 10 pc resolution in one interarm and one arm field: if the interarm CO(1-0) breaks into diffuse low-level gas that CO(2-1) misses while the arm emission stays in compact bright clumps, then the low interarm R21 is in part a beam-filling artifact rather than a true change in the bulk gas state.","supporting_citations":[{"cited_title":"S., Leroy , A","cited_arxiv_id":null,"evidence_quote":"Provides an ALMA R21 map of NGC 3627 with similar variations, and raises the sensitivity-bias caveat that this paper discusses."}],"review_version":1}