{"id":"e16f4234-fbad-4e2b-b2dc-923611e12a21","arxiv_id":"2411.19489","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"This paper decomposes SiO emission in 118 massive star-forming regions and finds that broad, outflow-driven emission dominates the SiO luminosity, with nine sources likely affected by expanding ultra-compact H II regions.","lead":"This paper maps silicon monoxide (SiO) emission in 146 massive star-forming regions and separates it into broad and narrow velocity components. It finds that over 66% of the SiO luminosity comes from broad emission tied to outflows, while a handful of sources show SiO possibly produced by expanding ultra-compact H II regions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Headline 66% broad-SiO fraction is an unweighted mean of per-source contributions, not the luminosity-weighted global fraction asserted in the abstract.","rationale":"After reading the paper's methods and results, the load-bearing step for the central claim is the computation of the 66% statistic. Section 4.4 defines the narrow contribution per source (Table A1) and reports the mean of these values (33.78%); the abstract and conclusions convert this into 'above 66% of L_SiO across the whole sample can be attributed to broad SiO.' However, the mean of per-source ratios is not the aggregate luminosity fraction. A source with very low L_SiO and 50% narrow contribution contributes as much to the mean as a source with 100 times higher L_SiO and 25% narrow contribution, but contributes almost nothing to the total luminosity. The correct quantity for the claim is sum(L_broad)/sum(L_total), which is straightforward to compute from Table A2. Without it, the headline number is not established; the qualitative conclusion that broad components dominate may still hold, but the specific 'above 66%' is unsupported and the paper's quantitative comparison to previous work (e.g., Liu et al. 2022's 60%) is on shaky ground. I also considered the reader's weakest assumption about the H13CO+ line-width threshold. That concern is real but its likely effect is conservative: H13CO+ drawn from a shocked region or multiple velocity components would broaden the threshold, causing some genuinely broad SiO to be classified as narrow and thus reducing the measured broad fraction. The reported 66% would then underestimate the true broad contribution, which does not threaten the 'majority from outflows' conclusion. The luminosity-aggregation issue, by contrast, can push the headline number in either direction and is not yet tested. Hence I identify it as the single most load-bearing concern. The reader's verdict of CONDITIONAL remains appropriate; if the global fraction falls below 50%, the verdict should be reconsidered. I recommend no change to the reader's verdict, but request the additional computation.","tokens_in":45389,"tokens_out":7926,"duration_ms":66714,"concrete_test":"Recompute the luminosity-weighted global fraction from Table A2: sum over the 118 decomposed sources of L_SiO,broad divided by sum of (L_SiO,broad+L_SiO,narrow). Compare with 66.22% = 100 - 33.78%. Also bootstrap over sources (10^4 resamples) to place an uncertainty on the global fraction. If the global fraction is below 66%, revise the abstract and conclusions; if it is below 50%, the claim that most SiO luminosity traces outflows is not supported by this statistic. Report both the per-source mean and the global fraction explicitly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim, that 'above 66%' of the sample's SiO luminosity arises from broad components, is computed in Section 4.4 as the unweighted mean of per-source narrow-component contributions (the 'Contribution N' column of Table A1; mean 33.78%). This is not the same as the global luminosity fraction sum(L_broad)/sum(L_broad+L_narrow) across the 118 decomposed sources, which is the quantity the abstract's 'majority of L_SiO across the whole sample' requires. Low-luminosity sources with narrow fractions near 50% (I08076-3556, I16297-4757) receive equal weight with luminous broad-dominated sources (e.g., I18056-1952 at 28% narrow, I15290-5546 at 23%), so the two statistics can differ materially. Since the 66% figure anchors the abstract, the conclusions, and the comparison with previous work, the paper must report the luminosity-weighted fraction and its uncertainty. The H13CO+ line-width threshold concern raised by the reader is secondary: if H13CO+ at position B is broadened by outflows or multiple velocity components, the threshold is overestimated, so some genuine broad SiO is misclassified as narrow, biasing the broad fraction downward; this would make the 66% a conservative lower limit rather than threaten the outflow conclusion.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analysis of ALMA Band 3 observations of SiO(2-1) and H13CO+(1-0) toward 146 massive star-forming regions from the ATOMS survey. SiO is detected in 136 regions, and 118 are decomposed into broad and narrow Gaussian components using scousepy, with the broad/narrow classification based on comparison with the H13CO+ line width measured at the position of broadest SiO. The authors report that the majority of the SiO luminosity (above 66%) comes from broad components, which they attribute to outflows; they further identify nine 'surrounded' UC HII regions where SiO may trace expanding HII regions. They also examine correlations of SiO luminosity with bolometric luminosity and with the presence of UC HII regions.","tokens_in":45638,"tokens_out":7282,"duration_ms":60291,"significance":"If correct, the paper would provide the first large-sample, high-resolution census of the origin of SiO emission in massive star-forming regions, supporting the use of SiO as an outflow tracer while identifying a minority population in which SiO traces HII-region expansion. The sample size (118 decomposed sources) and the use of public ALMA data and public spectral-fitting tools (scousepy, MADCUBA) are strengths. However, the headline 66% statistic is computed as an unweighted mean of per-source narrow-component fractions rather than as a luminosity-weighted global fraction, and the threshold-based decomposition rests on an untested assumption about H13CO+ line widths. These issues affect the quantitative central claim and require revision before the result can be accepted as stated.","major_comments":[{"comment":"The statement that 'the majority of L_SiO (above 66%) can be attributed to broad SiO' is computed as the unweighted mean of the per-source narrow-component contributions (the 'Contribution N' column in Table A1; mean 33.78%), not as the luminosity-weighted global fraction sum(L_broad)/sum(L_broad + L_narrow) across the 118 decomposed sources. These two statistics can differ materially: sources with narrow fractions near 50% (I08076-3556, I16297-4757) receive the same weight as luminous broad-dominated sources (I18056-1952, 28.3% narrow; I15290-5546, 22.6% narrow). Since the 66% figure anchors the abstract, Section 5.1, and the conclusions, the authors must report the luminosity-weighted fraction and its uncertainty (including fitting, distance, and flux-calibration contributions), and revise the text if that fraction differs materially.","section":"Section 4.4 / Abstract"},{"comment":"The broad/narrow classification uses the H13CO+ line width measured at position B as an upper limit for the quiescent line width across each source. This assumption is load-bearing: if H13CO+ at B is broadened by outflow entrainment or by multiple velocity components along the line of sight, the threshold is overestimated; if it is narrower than the ambient line width elsewhere, the threshold is underestimated. In neither case is the direction of the bias in the broad luminosity fraction established. I request a sensitivity analysis (e.g., repeating the decomposition with thresholds of 0.5x, 1x, and 2x the adopted FWHM_B, or using FWHM_N from Table A1) and reporting the resulting range of the broad fraction. Without this, the quantitative claim is not robust.","section":"Section 4.1"},{"comment":"The identification of nine 'surrounded' UC HII sources as having SiO produced by expanding HII regions is based on spatial morphology (SiO surrounding the 3 mm continuum and associated with 8 um and 1.3 GHz emission). No quantitative kinematic or chemical diagnostics are presented to distinguish HII-expansion shocks from outflows driven by embedded low-mass protostars in the surrounding gas, which is the standard interpretation of SiO in such environments. The paper would be considerably strengthened by a test such as comparing the SiO line widths, peak velocities, and abundances of these sources with predictions for expanding HII shells, or by showing the absence of 4.5 um outflow signatures toward them. As it stands, this claim is a plausible speculation rather than a demonstrated result, and the abstract and conclusions should reflect that distinction.","section":"Section 5.4 / Section 4.3.3"}],"minor_comments":[{"comment":"The section title 'H13CO+mission in active and quiescent regions' should read 'H13CO+ emission in active and quiescent regions'.","section":"Section 4.1"},{"comment":"The caption describes 'red filled circles' while the text refers to 'red empty circles'; please make the symbol description consistent.","section":"Figure 2"},{"comment":"The typesetting of the prefactors in Equations (3) and (4) is garbled (e.g., '2×10^11 2kT_ex...'); please correct the math display.","section":"Equations (3)-(4)"},{"comment":"There is a typo 'catergory' in the first paragraph; also, units such as 'Mo' should be written as 'M_sun' in Section 2.1 and Table A1.","section":"Section 4.2"},{"comment":"The KS tests compare sources classified as having outflows (based on Baug et al., in prep) against those without; since that classification uses SiO and HCO+ maps, the statistical comparison is not fully independent. Please note this limitation explicitly.","section":"Section 4.3.1"}],"recommendation":"major_revision","confidential_remarks":"The unweighted-mean issue with the 66% statistic is likely to be noticed by readers and should be fixed before publication; the luminosity-weighted global fraction should be computed and quoted in the abstract. The 'surrounded' HII region claim is speculative and should be toned down or supported with quantitative diagnostics. Also, the analysis relies on an unpublished paper (Baug et al., in prep) for the outflow classification; please ensure that this dependency is clearly flagged and that the relevant data are made available."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's real contribution is the systematic decomposition of SiO (2-1) profiles into broad and narrow components across 118 massive star-forming regions at 2.5\". That is new. So is the morphological classification of UC HII sources into coincident, offset, and surrounded, with the surrounded group plausibly tracing expanding HII regions. The data work is careful: public tools (scousepy, MADCUBA), tabulated per-source fits, ALMA archive availability, and a transparent description of the 18 sources excluded for absorption. No circular derivation. The qualitative case that most SiO arises from outflows is well supported by line-width distributions, spatial association with IR sources and filament skeletons, and the general agreement with previous work by Cosentino and others. I would not hesitate to cite this paper for the sample statistics and the new morphological classes.\n\nNow the soft spots, in proportion. The stress-test concern is correct: Section 4.4 computes the mean of per-source narrow contributions (33.78%) and then the abstract states that \"above 66%\" of L_SiO across the whole sample is broad. Those are different statistics. The correct global figure is sum(L_broad)/sum(L_broad+L_narrow) over the 118 sources. Looking at Table A2, low-luminosity sources with near-50% narrow fractions (I08076-3556, I16297-4757) get equal weight with luminous broad-dominated sources (I18056-1952 at 28% narrow). The unweighted mean is not what the abstract claims. That said, because the luminous sources are mostly broad-dominated, the luminosity-weighted broad fraction is likely at least as high as 66%—so the conclusion is probably robust, but the paper must report the proper weighted fraction and its uncertainty. This is an easy fix and an important one for the abstract, conclusions, and comparison with previous work.\n\nThe H13CO+ threshold assumption is secondary. If H13CO+ at position B is broadened by outflow entrainment, the threshold is overestimated and real broad SiO gets misclassified as narrow, biasing the broad fraction down. That makes the 66% a conservative lower limit rather than a threat to the outflow conclusion. Worth a sensitivity test, but not a fatal flaw.\n\nTwo other concerns are minor but worth noting. First, the outflow catalog is from Baug et al. (in prep), unpublished. The KS tests and the \"with/without outflow\" labels lean on that work; the authors should either present enough of it in the appendix or mark the conclusion as provisional. Second, the evolutionary-stage interpretation is not supported by the lack of correlation between L_SiO and Lbol/M, which the paper itself acknowledges. The proposed age ordering of the three UC HII classes is speculative, and the dismissal of the null result is a bit too quick.\n\nOverall, this is a solid survey paper that deserves a serious referee. I would send it to peer review with a request to replace the 66% statistic with the luminosity-weighted value and uncertainty, and to tone down the evolutionary claims. The main conclusion—that SiO is a robust outflow tracer over a large sample, with a minority population tracing HII expansion—will survive that revision.","headline":"The core result—most SiO traces outflows—is solid, but the headline 66% statistic is an unweighted per-source mean, not the luminosity-weighted fraction the abstract claims; fix that before accepting.","tokens_in":46331,"tokens_out":2670,"would_cite":true,"duration_ms":23682,"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":"Across 118 massive star-forming regions, more than two-thirds of silicon monoxide emission comes from broad velocity components tracing protostellar outflows.","keywords":["SiO emission","molecular outflows","UC HII regions","massive star formation","shock tracer","line profile decomposition","H13CO+","ATOMS survey"],"falsifier":"Measure H$^{13}$CO$^{+}$ (1-0) line widths at several quiescent positions away from outflows in each of the 118 sources and compare them with the single value taken at the broadest-SiO position; if many quiescent positions are appreciably narrower, the broad/narrow classification threshold moves and the 66 percent broad-luminosity fraction would change.","tokens_in":45164,"feed_emoji":"🌌","tokens_out":6770,"duration_ms":51595,"temperature":0.7,"pith_summary":"This paper asks where silicon monoxide (SiO) emission in massive star-forming regions actually comes from, using ALMA maps of 146 regions from the ATOMS survey. The authors decompose SiO (2-1) line profiles into broad and narrow Gaussian components and use the H$^{13}$CO$^{+}$ (1-0) line width as the dividing line between shock-broadened and quiescent gas. They find that broad SiO accounts for more than 66 percent of the total SiO luminosity across 118 decomposed sources, tying most SiO to protostellar outflows, while narrow SiO contributes about 34 percent and tracks gentler shocks. For nine sources whose SiO surrounds the ionized gas of an ultra-compact H II region, the paper argues the SiO is generated by the expansion of the H II region rather than by outflows. If correct, this makes SiO a reliable outflow tracer over large samples while flagging a minority population where SiO traces H II-region feedback.","feed_headline":"Over two-thirds of SiO emission across 118 regions traces outflows","feed_subtitle":"Broad SiO components carry 66% of the light; nine compact H II regions show SiO from expansion instead.","key_machinery":"The load-bearing object is the broad/narrow decomposition of SiO (2-1) line profiles carried out pixel-by-pixel with scousepy, a spectral-line fitting tool. Each SiO component is classified by comparing its Gaussian FWHM with the H$^{13}$CO$^{+}$ (1-0) line width measured at the position of broadest SiO; components wider than that threshold count as broad (shock-broadened) and narrower components count as narrow. A second classification, applied to the 58 sources with UC H II regions, sorts them by the spatial relation between SiO emission and the half-peak contour of 3 mm continuum: coincident, offset, and surrounded. The surrounded class, where SiO rings the ionized gas and aligns with 8 micron and 1.3 GHz emission, is the evidence for H II-expansion-driven SiO.","core_discovery":"The central claim is that the majority of SiO (2-1) luminosity in the ATOMS sample, above 66 percent, belongs to broad velocity components produced by strong protostellar outflows. Narrow components, contributing roughly one third of the luminosity, peak near ambient velocities and arise from lower-velocity shocks or young magnetic-precursor shocks. The paper additionally identifies nine 'surrounded' ultra-compact H II sources in which SiO emission encircles the 3 mm continuum and aligns with 8 micron and 1.3 GHz emission, interpreting this as SiO produced by the expanding H II region compressing molecular gas. Sources hosting UC H II regions show a weaker $L_{\\rm bol}$-$L_{\\rm SiO}$ correlation and higher SiO luminosity than sources without them, which the paper attributes to UV photochemistry from the H II regions. The result establishes SiO as a strong outflow tracer across a 118-source sample while delimiting a small subset where SiO traces H II expansion.","pith_inferences":["If the 66 percent broad fraction holds under a stricter quiescent-line-width calibration, SiO luminosity could become a quantitative proxy for outflow power in unresolved or extragalactic observations.","The nine surrounded sources predict that expanding H II regions produce SiO shells at the interface between ionized gas and molecular cloud; a targeted search for such shells in other UC H II samples would test this minority channel directly.","Because only two sources show narrow SiO exceeding half of the total luminosity, the magnetic-precursor-dominated phase implied by the narrow component may be short; a larger sample with evolutionary indicators could quantify that phase's duration.","Combining the kinematic decomposition with H40alpha emission measures and PDR tracers could separate UV-photochemistry effects from pure shock effects in the UC H II subsample."],"forward_implications":["SiO (2-1) can be used as a statistically reliable outflow tracer in massive star-forming samples, since broad components carry the majority of its luminosity.","Narrow SiO emission, about one third of the luminosity, traces slower or younger shocks near ambient velocities, so it can complement broad SiO in evolutionary studies.","In the nine 'surrounded' UC H II sources, SiO traces expanding H II regions, so those sources should be set aside when SiO is used to measure outflow activity.","The weaker $L_{\\rm bol}$-$L_{\\rm SiO}$ correlation in UC H II sources indicates that UV photochemistry modifies SiO emission at later evolutionary stages.","The group A/group B spatial classification (broad more extended in 94 sources, narrow more extended in 24) provides a statistical baseline for testing outflow evolution with larger samples."],"supporting_citations":[{"why":"Establishes SiO as a shock tracer with strongly enhanced abundance in outflow regions.","marker":"Martin-Pintado et al. 1992"},{"why":"Provides the grain sputtering mechanism that releases silicon into the gas phase in shocks.","marker":"Schilke et al. 1997"},{"why":"Models the narrow SiO component as arising from magnetic-precursor shocks in young outflows.","marker":"Jiménez-Serra et al. 2004, 2005"},{"why":"Supplies the SiO-versus-H13CO+ line-width comparison method and links narrow SiO to H II region expansion.","marker":"Cosentino et al. 2018, 2020"},{"why":"Connects broad SiO components to high-velocity shocks from powerful outflows in high-mass star-forming regions.","marker":"Duarte-Cabral et al. 2014"},{"why":"Earlier ACA observations found two SiO velocity components; this work extends that analysis at higher angular resolution.","marker":"Liu et al. 2022"},{"why":"Provides the H13CO+ filament skeletons used to test the spatial association of SiO emission.","marker":"Zhou et al. 2022"},{"why":"Supplies the H40alpha UC H II identifications and emission measures used to classify the UC H II sources.","marker":"Zhang et al. 2023"},{"why":"Provides the scousepy spectral decomposition tool used for the pixel-by-pixel Gaussian fitting.","marker":"Henshaw et al. 2016"}],"fun_headline_variants":["Outflows produce 66% of SiO light in ATOMS massive clumps","SiO mostly from outflows, except nine expanding HII regions","ALMA survey: SiO traces outflows, with a few HII-region surprises","Two-thirds of SiO emission in massive clumps comes from outflows","Most SiO in star-forming regions arises from outflows, not ambient gas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that the H$^{13}$CO$^{+}$ line width measured at the position of the broadest SiO emission is an upper limit to the quiescent gas line width across each source, so any SiO component wider than that value is necessarily shock-broadened.","fun_headline_variants_meta":{"raw":{"variants":["Outflows produce 66% of SiO light in ATOMS massive clumps","SiO mostly from outflows, except nine expanding HII regions","ALMA survey: SiO traces outflows, with a few HII-region surprises","Two-thirds of SiO emission in massive clumps comes from outflows","Most SiO in star-forming regions arises from outflows, not ambient gas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000594,"raw_usage":{"total_tokens":2852,"prompt_tokens":1088,"completion_tokens":1764,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":704,"completion_tokens_details":{"reasoning_tokens":1666}},"tokens_in":704,"tokens_out":1764,"duration_ms":12549,"temperature":1.0,"reasoning_tokens":1666,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:07:31.550128+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure H$^{13}$CO$^{+}$ (1-0) line widths at several quiescent positions away from outflows in each of the 118 sources and compare them with the single value taken at the broadest-SiO position; if many quiescent positions are appreciably narrower, the broad/narrow classification threshold moves and the 66 percent broad-luminosity fraction would change.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Models the narrow SiO component as arising from magnetic-precursor shocks in young outflows."}],"review_version":1}