{"id":"e59d2dd0-b2a0-427f-bf2e-f9a6bc75c4d2","arxiv_id":"2505.08912","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The CO(2-1) to CO(1-0) line ratio drops by about 20 percent within 40 percent of the optical radius in barred galaxies while remaining roughly flat in unbarred galaxies.","lead":"This paper maps the ratio of two carbon monoxide emission lines across 11 nearby galaxies and finds that the ratio declines with radius in barred galaxies but stays nearly constant in unbarred ones. If the trend holds, surveys that use only the higher-frequency CO line to estimate molecular gas mass would need to correct for galaxy morphology.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The bimodal r21 gradient claim is a post hoc visual split of 11 galaxies without a slope fit or significance test; the 20% drop and SA/SB separation are not demonstrated statistically.","rationale":"The reader's strongest claim matches the abstract. The paper has genuine strengths: ALMA TP calibration (2.4%), consistent pipeline, honest discussion of arbitrary thresholds, and the outlier-vs-disk r21 offset is consistent across all 11 galaxies. However, the new headline result—the bimodal radial gradient separating SA from SB—rests on visual inspection of normalized profiles with no formal statistics. I consider this the most load-bearing concern because, even if the tilted-ring decomposition were perfect, a post hoc split of 11 galaxies into two visually distinct groups is not itself evidence of bimodality. The paper contains no test that the two groups' gradients differ, no uncertainty on the 20% figure, and no sensitivity analysis for the reclassification of NGC 3521 or exclusion of NGC 628. The tilted-ring concern raised by the reader is real but secondary: mis-assigned non-circular gas could bias the disk profiles, but the direction is not obviously aligned with the claimed effect, and the paper's own threshold sensitivity discussion shows the decomposition is fragile for at least one galaxy (NGC 4579). The concrete test I propose would settle the statistical issue directly on the existing data. Hence I agree with the reader's CONDITIONAL verdict; the claim should be quantified or softened before the paper is accepted as is.","tokens_in":19871,"tokens_out":10366,"duration_ms":105489,"concrete_test":"Obtain or digitize the normalized disk r21 profiles for the 11 galaxies in Figure 4 (top-right). For each galaxy, fit a linear slope over 0-0.4 R25 using bootstrap resampling of the underlying PV pixels to get slope uncertainties. Then run a permutation test comparing the slopes of the shallow group (SA plus NGC 3627/NGC 4321) against the steep group (SB plus NGC 4536/NGC 4579). As a secondary diagnostic, apply a Hartigan dip test to the distribution of the 11 slopes to assess bimodality. Rerun both tests after reclassifying NGC 3521 as SAB (per RC3) and, if possible, after measuring NGC 628. If the group separation p-value is >0.05 or the dip test is non-significant, the bimodality claim should be softened to a tentative trend.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (abstract; Section 4.3; Figure 4) is that the radial gradient of the disk r21 is bimodal: all SA galaxies (plus two SABs) have shallow gradients out to 0.4 R25 while all SB galaxies (plus two other SABs) drop by ~20%. This grouping is presented visually, with no slope fitted to the normalized profiles, no uncertainties on the 20% drop, and no statistical test comparing the two groups. The only significance test in the paper (Welch's t-test in Section 4.2) concerns the disk vs outlier r21 distributions, not the radial gradients. With n=11, a post hoc split into 6+5 after inspecting the profiles can arise by chance; the probability is not evaluated. The grouping is also sensitive to classification decisions: NGC 3521 is reclassified from SAB to SA (Section 2); NGC 628 is excluded (Section 2); and Section 4.1 states that lowering the threshold for NGC 4579 from 10 sigma to 5 sigma reclassifies the entire PV diagram as disk, which would move a steep-group member. The normalization to the central 28 arcsec value (Figure 4) introduces correlated, unpropagated errors. Thus the headline bimodality is not established at the claimed confidence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents ALMA Total Power observations of 12CO(J=1-0) for 12 nearby galaxies from the FACTS survey and compares them with archival ALMA TP 12CO(J=2-1) data at approximately 56 arcsecond resolution. After constructing ratio maps and major-axis position-velocity diagrams, the authors fit tilted-ring models to the CO(1-0) cubes with 3DBarolo, separate each galaxy into a rotating 'disk' component and a kinematic 'outlier' component, and measure the line ratio r21 in each component. They find that outliers have systematically higher r21 than the disk (Welch's t-test, Section 4.2), that the average r21 does not depend significantly on bar morphology, and that the radial gradient of the disk r21 is bimodal: galaxies classified as SA plus two SABs show constant or shallow gradients out to 0.4 R25, while all SB galaxies plus the remaining two SABs decline by about 20% within 0.4 R25, with the decline plausibly associated with the stellar bar radius. The paper also gives a crude estimate that a radially varying r21 would raise the Schmidt-Kennicutt index by about 20% for barred galaxies.","tokens_in":20092,"tokens_out":4710,"duration_ms":49063,"significance":"If the bimodal gradient is real, it is an important result: it would demonstrate that CO(2-1)/CO(1-0) varies systematically with bar morphology on kiloparsec scales, with direct implications for CO(2-1)-only molecular gas mass estimates and for the slope of the Schmidt-Kennicutt relation. The paper's strengths are its carefully documented ALMA TP calibration (execution-block scaling, ozone-line removal, propagation of flux uncertainties), its explicit separation of disk and outlier kinematic components, and the consistency of the global r21 values with earlier work. The data handling is transparent enough that the central claim can be tested, but that test is not currently supplied in the manuscript.","major_comments":[{"comment":"The headline bimodality is established only by visual inspection. No slope is fitted to the normalized radial profiles, the quoted ~20% decline is not accompanied by an uncertainty, and no statistical test compares the shallow-gradient group with the steep-gradient group. With n=11 and a grouping decided after inspecting the profiles, a chance separation is not excluded. Please add a quantitative definition of the two groups (for example, a fitted linear slope over 0-0.4 R25 with its uncertainty, or a likelihood-ratio/bimodality test) and report the resulting significance.","section":"Section 4.3, Figure 4"},{"comment":"The morphology grouping is sensitive to classification decisions that are not tested. NGC 3521 is moved from SAB to SA (Section 2), NGC 628 is excluded from the sample (Section 2), and Section 4.1 reports that lowering the disk/outlier threshold for NGC 4579 from 10 sigma to 5 sigma reclassifies the entire PV diagram as disk; NGC 4579 is in the steep-gradient group. Please show that the bimodality and its grouping survive these choices (for example, using RC3 classifications unchanged, or excluding NGC 4579), or soften the claim accordingly.","section":"Section 4.1; Section 2; Table 1"},{"comment":"The disk component itself may be contaminated by non-circular gas in barred galaxies because the tilted-ring model fixes radial inflow/outflow to zero and uses a sigma contour to separate disk from outliers. Since the bimodality claim is based solely on the disk r21 profiles, systematic mis-assignment of non-circular gas in bars could create or enhance the steep decline. I ask for a robustness check, for example refitting with radial motions allowed or with a different contour level, and showing that the disk r21 profiles are unchanged.","section":"Section 4.1"},{"comment":"Normalizing each profile by the value averaged over the central 28 arcseconds introduces correlated errors that are not propagated, and adjacent radial bins spaced by 5.6 arcseconds (one tenth of the beam) are not independent. The quoted ~20% decline needs an uncertainty that accounts for this covariance, and the central normalization should be checked against an alternative reference, such as a fitted intercept rather than a single-beam average.","section":"Figure 4, Section 4.3"}],"minor_comments":[{"comment":"The beam size notation uses '56.6' without an arcsecond unit in several places; please standardize this as 56.6 arcseconds (or 56\".6) throughout.","section":"Section 2.2"},{"comment":"The Schmidt-Kennicutt slope estimate contains garbled mathematics in the printed text: '~log 0.050.7 0.63' and '~log 0.10.63 0.5' appear to be missing fraction formatting. Please typeset these as log10(0.7/0.63) and log10(0.63/0.5), or equivalent.","section":"Section 4.4"},{"comment":"The caption says the right-hand panels show the radial change of r21, but it does not define the plotted error bars beyond 'standard deviation of the values on the PV diagram,' nor does it explain whether the disk+outlier values (blue diamonds) are area-weighted in the same way as the disk and outlier values. Please clarify.","section":"Figure 2 caption"},{"comment":"The footnote for NGC 3521 contains a typo, 'treat is as a nonbarred SA,' which should read 'treat it as a nonbarred SA.'","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"This is a solid observational paper with valuable ALMA TP data and careful calibration work. The main issue is that the central bimodality claim is presented with visual grouping and no quantitative significance test; this is fixable within the scope of the manuscript. I would be satisfied if the authors add a slope-based or likelihood-based test of the two groups and robustness checks for the classification and threshold choices. I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi — short take on Komugi et al., ALMA FACTS II.\n\nWorth knowing upfront: this is a careful, well-calibrated data paper, but the headline claim — that the r21 radial gradient splits cleanly between barred and unbarred galaxies — is not backed by any statistical test. Treat it as a suggestion, not a result.\n\nWhat's genuinely good: both CO lines come from ALMA, with the TP array data reduced carefully and the relative flux scales cross-checked. The global average r21 = 0.61 ± 0.06 matches earlier work, so the measurement passes sanity checks. The disk/outlier separation with 3DBarolo is a reasonable idea, and the finding that outlier components have higher r21 in every single galaxy is new and plausible. The SK slope discussion is explicitly crude, and the authors say so — that's honest.\n\nThe weak part is the bimodality claim. Figure 4 normalizes each galaxy to its central value and then the authors visually split the 11 galaxies into two groups: SA plus two SABs flat, SB plus two other SABs declining ~20% by 0.4 R25. No slope fits, no error bars on the decline, no significance test between groups. With n=11 and a split chosen after looking at the data, that separation isn't established. The classification is fragile: NGC 3521 is reclassified from SAB to SA; NGC 628 is excluded; NGC 4579's disk/outlier division changes entirely if the threshold drops from 10σ to 5σ. The paper acknowledges some of this, which is good, but the abstract still sells the bimodality as the main result.\n\nThe outlier-disk difference also rests on correlated pixels, so the p=0.009 from Welch's t-test overstates the evidence; still, the consistency across galaxies makes it worth reporting.\n\nBottom line: send to peer review if you're the editor. The data and outlier result deserve publication. But ask the authors to fit the gradients properly, put error bars on the 20%, run a permutation test on the group separation, and soften the language. After that, it's a solid paper even if the bimodality doesn't survive.","headline":"A careful data paper with a new and plausible outlier-disk r21 effect, but the headline barred/unbarred bimodality is a post hoc visual split with no statistical support.","tokens_in":20697,"tokens_out":3333,"would_cite":true,"duration_ms":32748,"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":"Barred spiral galaxies show a roughly 20% decline in the CO line ratio inside the bar radius.","keywords":["CO line ratio","molecular gas","galactic bars","tilted-ring modeling","ALMA Total Power","galaxy kinematics","Schmidt-Kennicutt relation","nearby galaxies"],"falsifier":"A decisive test is to measure $r_{21}$ at about 100 parsec resolution along the bars of several SB galaxies using ALMA's 12 m array together with the Total Power array, and to model the gas kinematics with a bar-potential fit instead of a pure tilted ring. The bimodality claim predicts low $r_{21}$ throughout the bar and a rise near the bar ends; if the roughly 20% decline before 0.4 $R_{25}$ disappears or does not track the bar radius under that model, the claim fails.","tokens_in":19647,"feed_emoji":"📡","tokens_out":11237,"duration_ms":98350,"temperature":0.7,"pith_summary":"This paper asks whether the ratio between the two carbon monoxide emission lines $^{12}$CO($J=2{-}1$) and $^{12}$CO($J=1{-}0$), called $r_{21}$, varies systematically across galaxy disks and with galaxy morphology. Using ALMA Total Power observations of 11 nearby galaxies, the authors find an average luminosity-weighted $r_{21}$ of 0.61, but the radial behavior splits into two groups. A group containing all SA galaxies and two weakly barred SAB galaxies keeps $r_{21}$ nearly constant out to 40% of the optical radius, while all strongly barred SB galaxies and the remaining SABs show a steep decline of roughly 20% over the same range, at about the radius of the stellar bar. The result matters because many surveys observe only the $J=2{-}1$ line and convert to molecular gas mass with a constant ratio; a systematic radial and morphology-dependent variation would bias those conversions and affect inferred star formation scaling relations.","feed_headline":"CO line ratio falls 20% inside the bars of spiral galaxies","feed_subtitle":"ALMA mapping of 11 galaxies shows the ratio is flat in unbarred disks but drops within the bar radius.","key_machinery":"The central object is the line ratio $r_{21}=I_{\\mathrm{CO}(2{-}1)}/I_{\\mathrm{CO}(1{-}0)}$, measured from ALMA Total Power maps of CO(1–0) and matched CO(2–1) cubes, both smoothed to the same roughly 57 arcsecond beam. The key tool is the position–velocity diagram along the major axis, compared with a tilted-ring model of pure circular rotation fit to the CO(1–0) cube; the model contour at 30$\\sigma$ (10$\\sigma$ for two galaxies) defines the disk, and everything outside is a kinematic outlier. Radial profiles of the disk component, normalized to the central-beam value and to the optical radius $R_{25}$ or to the bar radius from the $m=2$ Fourier amplitude, are what expose the bimodal gradient. The load-bearing step is the assumption that this decomposition cleanly separates disk gas from non-circular components, so the radial profiles reflect the disk rather than contamination by nuclear or bar-driven gas.","core_discovery":"The central claim is that the radial gradient of $r_{21}$ is bimodal in nearby galaxies. After fitting a tilted-ring model to the CO(1–0) cube with the 3DBarolo code and using a 30$\\sigma$ (or 10$\\sigma$, for two galaxies) model contour to separate the rotating disk from kinematic outliers, the authors find two clearly separated classes: disks that keep $r_{21}$ constant or nearly flat out to 0.4 of the optical radius $R_{25}$, and disks that fall by about 20% before reaching 0.4 $R_{25}$, with the break located at the characteristic radius of the stellar bar. All SA galaxies fall in the flat class, all SB galaxies in the declining class, and the two SAB galaxies split between the classes. The same decomposition shows that kinematic outliers, gas at velocities that deviate from pure circular rotation, have systematically higher $r_{21}$ than the disk in every galaxy, with average values of 0.69 versus 0.62. The paper concludes that large-scale $r_{21}$ variation is real and morphology-dependent, not a calibration artifact, and is large enough to matter for molecular mass estimates made from CO(2–1) alone.","pith_inferences":["If confirmed with a larger sample, the bimodality predicts that barred galaxies will appear to have systematically shallower radial CO(2–1) extent than unbarred galaxies of similar mass, because low-$r_{21}$ outer gas is undercounted relative to the center.","A clean cross-check is available in existing data: the Schmidt–Kennicutt index for barred galaxies should be higher when molecular masses come from CO(1–0) than when they come from CO(2–1) at fixed ratio.","The two SAB galaxies in the steep-gradient group, NGC 4536 and NGC 4579, are natural targets for high-resolution $r_{21}$ mapping: the bimodality picture predicts sub-kpc structure resembling SB galaxies, while the shallow-gradient SABs should resemble SA disks.","Because barred and unbarred galaxies have different metallicity gradients, the net morphology dependence of inferred molecular gas fractions could partly cancel; separating the two requires abundance measurements along the same radii."],"forward_implications":["Surveys that observe only CO(2–1) and assume a constant $r_{21}$ will overestimate molecular gas surface density in the centers of barred galaxies and underestimate it at larger radii, flattening the apparent range of gas densities.","The Schmidt–Kennicutt power-law index $N$ inferred from CO(2–1) at fixed $r_{21}$ will be systematically too low for barred galaxies; the paper's simple estimate turns $N=1$ into $N\\simeq1.17$ when a radially varying $r_{21}$ is used.","The absence of a difference in average $r_{21}$ between SA, SAB, and SB galaxies implies that the bar effect is not a global change in molecular gas excitation but a redistribution of the ratio with radius.","Beyond the bar radius, even the steep-decline galaxies return to a shallow or constant $r_{21}$ gradient, so the morphology dependence is localized to the bar region and the central disk."],"supporting_citations":[{"why":"Supplies the PHANGS CO(2–1) Total Power data cubes that the paper matches to its new CO(1–0) maps.","marker":"A. K. Leroy et al. (2021)"},{"why":"Documents the global $r_{21}$ measurements and the roughly 30% inter-telescope flux scatter that this study's ALMA calibration is designed to overcome.","marker":"A. K. Leroy et al. (2022)"},{"why":"Provides the 3DBarolo code used to fit the tilted-ring model that separates disk from outlier components.","marker":"E. M. Di Teodoro & F. Fraternali (2015)"},{"why":"Gives the large velocity gradient calculations that map $r_{21}$ to gas density and kinetic temperature, the physical interpretation used throughout.","marker":"S. Sakamoto et al. (1994)"},{"why":"Shows that replacing CO(1–0) with CO(2–1) at fixed $r_{21}$ changes the Schmidt–Kennicutt index $N$, the direct consequence this paper extends.","marker":"Y. Yajima et al. (2021)"},{"why":"Provides the high-resolution $r_{21}$ map of the barred galaxy NGC 1365 with low ratios inside the bar, the comparison case for the new large-scale barred gradients.","marker":"F. Egusa et al. (2022)"},{"why":"Supplies the characteristic bar radii from the $m=2$ Fourier amplitude used to normalize the radial profiles.","marker":"S. Díaz-García et al. (2016)"},{"why":"Establishes the ALMA Total Power reduction and ratio method in M83 that the FACTS survey extends to a larger sample.","marker":"J. Koda et al. (2020)"}],"fun_headline_variants":["CO ratio drops 20% inside barred galaxies' bars","Barred disks show steep CO ratio drop","CO line ratio falls sharply at bar radius in spirals","In barred galaxies, CO ratio drops 20%","Barred spirals show 20% CO ratio decline in inner disk"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the tilted-ring model fit to CO(1–0), with radial inflow and outflow forced to zero and a 10–30$\\sigma$ contour threshold, cleanly separates the rotating disk from kinematic outliers, so the radial $r_{21}$ profiles of the disk are not biased by non-circular bar gas.","fun_headline_variants_meta":{"raw":{"variants":["CO ratio drops 20% inside barred galaxies' bars","Barred disks show steep CO ratio drop","CO line ratio falls sharply at bar radius in spirals","In barred galaxies, CO ratio drops 20%","Barred spirals show 20% CO ratio decline in inner disk"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000749,"raw_usage":{"total_tokens":3419,"prompt_tokens":1113,"completion_tokens":2306,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":729,"completion_tokens_details":{"reasoning_tokens":2226}},"tokens_in":729,"tokens_out":2306,"duration_ms":16490,"temperature":1.0,"reasoning_tokens":2226,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:45:18.189248+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test is to measure $r_{21}$ at about 100 parsec resolution along the bars of several SB galaxies using ALMA's 12 m array together with the Total Power array, and to model the gas kinematics with a bar-potential fit instead of a pure tilted ring. The bimodality claim predicts low $r_{21}$ throughout the bar and a rise near the bar ends; if the roughly 20% decline before 0.4 $R_{25}$ disappears or does not track the bar radius under that model, the claim fails.","supporting_citations":[{"cited_title":"2022, ApJ, 935, 64","cited_arxiv_id":null,"evidence_quote":"Provides the high-resolution $r_{21}$ map of the barred galaxy NGC 1365 with low ratios inside the bar, the comparison case for the new large-scale barred gradients."}],"review_version":1}