{"id":"b57d1a92-acc9-40d3-9076-430658d89818","arxiv_id":"2505.04164","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"HC3N emission reveals 45 outflows in 146 massive star-forming clumps and appears to complement SiO, favoring low-velocity outflow gas.","lead":"This ALMA survey identifies 45 outflows, 44 of them bipolar, in 146 massive star-forming clumps using emission from the molecule HC3N. The study suggests HC3N can complement traditional outflow tracers such as SiO, especially for catching slow-moving outflow gas.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The low-velocity outflow claim hinges on the assumption that H13CO+ (J=1-0) is wing-free; if H13CO+ itself contains outflow emission, the VLSR±FWHM cut in Sec. 3.1 removes the slow HC3N gas the paper says HC3N recovers.","rationale":"The paper is a genuinely useful catalog: the 45 HC3N outflows are documented with moment maps, PV diagrams, opening angles, and quantitative correlations with clump properties. The concern is not whether HC3N outflows exist; it is whether the paper's signature claim about low-velocity recovery is supported by the reduction procedure. The reader's weakest_assumption identifies the same H13CO+ masking issue, and the recommended CONDITIONAL verdict already reflects the need for a check. I therefore see no reason to move the verdict. The separate weakness in the SiO comparison (KS p = 0.225; SiO catalog from a private communication) is real and should temper the wording about 'slightly better' performance, but it is secondary to the masking issue because the mask affects both the HC3N catalog and the comparison sample. If the proposed test shows that H13CO+ profiles are single-component and wing-free across the sample, the low-velocity claim would be substantially strengthened.","tokens_in":28740,"tokens_out":4953,"duration_ms":55194,"concrete_test":"For each of the 45 outflow fields, fit the H13CO+ (J=1-0) spectrum toward the driving core with a single Gaussian and with a two-component model (narrow core plus broad wing). If a broad component is detected above, say, 3σ in a substantial fraction of fields, refit the cloud velocity using only the narrow component and rerun the wing-identification and parameter pipeline from Secs. 3.1 and 3.4 with the narrower VLSR ± FWHM range. Then compare the re-derived vlobe distribution and cumulative low-velocity fractions with the published values; if the re-derived distribution contains additional low-velocity channels or sources, the current mask excludes the very emission used to support the low-velocity claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Secs. 3.1 and 3.4, the systemic/cloud velocity is fixed by a single Gaussian fit to H13CO+ (J=1-0), and all HC3N emission within VLSR ± FWHM is excluded from the outflow wings. The paper asserts that H13CO+ 'typically does not have contribution from outflowing material' but provides no diagnostic test for wings in the H13CO+ profiles. If H13CO+ has broad non-Gaussian emission (which can occur in shocked regions), the excluded interval is too wide: the slowest HC3N wing channels are discarded, terminal velocities are biased high, outflow mass/momentum/energy are underestimated, and the comparison with SiO uses differently masked populations. This cuts directly against the central claim that HC3N is especially useful for low-velocity components, because those are precisely the channels the procedure removes. The paper acknowledges the difficulty in Sec. 3.4 ('low-velocity outflowing gas always mixes with the cloud component'), but the key assumption that H13CO+ is uncontaminated is not tested. A quantitative test is therefore needed before the low-velocity claim can be regarded as secure.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a systematic search for molecular outflows in the HC3N (J=11-10) line toward 146 massive star-forming clumps from the ALMA-ATOMS survey. The authors identify 45 outflows (44 bipolar, 1 explosive) using wing emission outside a VLSR±FWHM interval set by a Gaussian fit to H13CO+ (J=1-0), with confirmation from moment maps and position-velocity diagrams. They derive outflow masses, momenta, energies, dynamical times, and related quantities under LTE, optically thin, and constant abundance assumptions, and compare the resulting catalog against an SiO (J=2-1) outflow catalog from the same survey. The paper argues that HC3N is a useful complementary outflow tracer, particularly for detecting low-velocity outflow components, and that outflow detection does not correlate with clump evolutionary stage.","tokens_in":28989,"tokens_out":3020,"duration_ms":33079,"significance":"If the central claims hold, this would be the first large-sample demonstration that HC3N can serve as an outflow tracer in massive star-forming regions, complementing SiO and CO by tracing slower, less contaminated gas. The paper provides a substantial catalog (Table 2) with derived parameters and identifies 10 sources with separate fast-jet and slow-wide components, which is of interest to models of jet- and wind-driven outflows. The analysis is systematic and uses public survey data, with clear figures and a reproducible methodology. However, the headline claim that HC3N is 'slightly better' at detecting low-velocity components rests on a Kolmogorov-Smirnov test that does not reach statistical significance (p=0.225), and the main comparison relies on an unpublished SiO catalog. The outflow parameter derivation also carries a large systematic uncertainty from the assumed HC3N/H2 abundance (spanning about a factor of 280 between literature values), which is acknowledged but not propagated. The paper's value is primarily as a catalog; its interpretive claims currently outrun the statistical evidence.","major_comments":[{"comment":"The outflow definition uses a single Gaussian fit to H13CO+ (J=1-0) to set VLSR±FWHM as the 'cloud component' and removes all HC3N emission within that range. The paper states that H13CO+ 'typically does not have contribution from outflowing material' but provides no diagnostic test of this assumption. If H13CO+ itself has non-Gaussian wings in shocked regions, the excluded interval is overestimated, and the slowest HC3N wing channels are discarded, which would bias the terminal velocities high and directly undermine the paper's specific claim that HC3N recovers low-velocity outflow components. Since the central claim of the paper depends on this masking choice, a quantitative test is needed; for example, comparing the H13CO+ profiles against single-Gaussian models in the outflow regions, or re-deriving outflow properties with an alternative cloud-velocity definition.","section":"§3.1, §3.4"},{"comment":"The claim that HC3N traces lower-velocity components than SiO is supported by a KS test on the vlobe distributions with p=0.225, which does not reject the null hypothesis that the two samples are drawn from the same distribution at any conventional significance level. The cumulative fraction plot in Figure 9 (right) shows a visual offset, but the stated test does not establish a statistically significant difference. The authors should either apply a more sensitive test (e.g., a two-sample Anderson-Darling test or a comparison of the full velocity distributions with the masks taken into account) or, if the data support it, moderate the claim to 'no statistically significant difference, with a possible trend at low velocities'.","section":"§4.2.2, Figure 9"},{"comment":"The derived outflow masses, momenta, and energies assume a single [HC3N/H2] abundance of 5×10^-9, while the cited literature values range from 1.4×10^-8 (Mendoza et al. 2018) to 5.1×10^-11 (Taniguchi et al. 2018b), a factor of about 280. This systematic uncertainty is not propagated into the catalog values in Table 2 or into the correlation analyses. Because the abundance factor is roughly constant across sources, the Spearman correlations with clump mass and luminosity may be insensitive to it, but the absolute values (and the direct comparison with Towner et al. 2024 in Figure 5) are not. The paper should at least provide a quantitative statement of how the assumed abundance and Tex affect the reported ranges in Table 1 and the slopes in Equations (10)-(15).","section":"§4.1, Equations (10)-(15)"},{"comment":"The comparison with SiO outflows relies on a 'Baug et al. (private communication)' catalog that is not publicly available and is not described beyond the number of detected outflows (153). This makes the central comparison non-reproducible for other researchers. The SiO catalog should be published or otherwise made accessible, or the authors should re-run the comparison using a publicly available SiO outflow sample. At minimum, the paper must state the detection threshold and masking procedure used for the SiO catalog so that differences in sensitivity do not drive the reported comparison.","section":"§4.2.2"}],"minor_comments":[{"comment":"There is a typo: 'adpoted' should be 'adopted' in the sentence preceding Equation (1).","section":"§3.3"},{"comment":"The formatting of Table 2 is dense and some entries are unclear (e.g., the superscript 'a' marks and the missing PA values for the explosive outflow lobes). Please clarify the meaning of blank cells and ensure the table is readable in the published version.","section":"Table 2"},{"comment":"The wording 'the KS-test do not provide enough evidence to claim the similarity' is imprecise; for p=0.087 the test does not reject the null hypothesis of identical distributions, so one should say 'no significant difference was detected' rather than 'a moderate difference is present'.","section":"§4.2.1"},{"comment":"The paper describes the explosive outflow in IRAS 15520-5234 as one outflow but Table 2 lists seven components O1-O7. Please clarify whether these are seven lobes of one explosive outflow or seven separately identified outflows, and how the opening angle and PV analysis treat them.","section":"§3.1"},{"comment":"The text says 'we assumed a mean inclination angle, θ, of 53.7°' and uses it to correct parameters with powers of sinθ and cosθ; please state explicitly that the quoted vlobe values in Table 2 are the observed (not inclination-corrected) terminal velocities, since the comparison in Figure 9 appears to use these values.","section":"§3.4"},{"comment":"The notation 'HC 3N' and 'HC3N' is used inconsistently. Please use a single convention throughout the text, tables, and figures.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper's main scientific claim (low-velocity advantage) is currently supported mainly by a non-significant KS test and by a comparison with an unpublished SiO catalog from the same group. I would encourage the editor to require that the SiO catalog be released or at least fully described, and to ask the authors to add the H13CO+ wing test before the paper can be considered for publication. The catalog itself is useful and likely sound, but the interpretive framing needs strengthening."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful contribution is the catalog: 45 HC3N outflows among 146 massive clumps, with moment maps, PV diagrams, and derived parameters. That is the first systematic sample of this sort, and the figures look honest. The detection work is careful, with confirmed/probable flags, and the authors are transparent about HC3N's weaker sensitivity relative to SiO. The comparison with Towner et al. shows that their mass, momentum, and energy ranges are plausible.\n\nThe soft spot is the central claim. The KS test on terminal velocities gives p=0.225, so the distributions are not significantly different. Calling HC3N 'slightly better' at low velocities is an interpretation the data do not back at the 95% level. The cumulative fraction plot does show a modest excess at low velocities, so the claim is not baseless, but the statistics should be presented as suggestive, not as a result.\n\nThe masking assumption deserves a closer look. The cloud velocity range is set by a single Gaussian fit to H13CO+; everything inside VLSR±FWHM is cut. If H13CO+ has non-Gaussian wings in any of these regions, the slowest HC3N wing gas is removed, which would directly suppress the low-velocity signal the paper wants to measure. The authors acknowledge the mixing problem in Sec 3.4 but offer no test. This is fixable: check the H13CO+ profiles for residuals, or compare with an independent cloud velocity tracer.\n\nAlso, the SiO comparison depends on an unpublished private communication. That is a real problem for reproducibility. The comparison should be based on a public catalog or the SiO data should be released.\n\nOutflow parameters include large systematic uncertainties (abundance alone spans a factor of ~280), and there are no per-source errors. The authors disclose this, and it is standard for this kind of work, so I treat it as a limitation, not a flaw.\n\nOverall this is a useful catalog paper with an overstated abstract. A referee can reasonably ask for a wing test on H13CO+, public access to the SiO sample, and softer language about the low-velocity advantage. I would send it to peer review.","headline":"The HC3N outflow catalog is a genuinely useful new sample, but the headline claim that HC3N beats SiO at low velocities rests on a non-significant KS test and an untested masking assumption.","tokens_in":29652,"tokens_out":3440,"would_cite":true,"duration_ms":31679,"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":"A dense-gas molecule usually used for cores, HC3N, also works as an outflow tracer and catches the low-velocity gas that SiO misses.","keywords":["protostellar outflows","HC3N","cyanoacetylene","massive star-forming regions","outflow tracers","low-velocity outflow gas","SiO","shock chemistry"],"falsifier":"Take the 45 HC3N outflow hosts and inspect the H13CO+ (J=1–0) spectra for non-Gaussian wings or asymmetries; if any host shows outflow wings in H13CO+ within the adopted VLSR ± FWHM, then the claimed low-velocity HC3N component is partly an artifact of choosing too wide a cloud window.","tokens_in":28562,"feed_emoji":"🔭","tokens_out":6144,"duration_ms":54210,"temperature":0.7,"pith_summary":"This paper argues that cyanoacetylene (HC3N), a molecule usually used to trace dense gas, can also serve as a practical tracer of protostellar outflows, and that it is especially good at detecting the slow-moving outflow gas that standard shock tracers such as SiO tend to miss. Working with ALMA data for 146 massive star-forming clumps from the ATOMS survey, the authors identify 45 outflows in HC3N (J=11–10), 44 of them bipolar and one explosive. They find that outflow detection does not correlate with the evolutionary stage of the host clump, which they interpret as HC3N being an unbiased outflow tracer. The study also shows that HC3N outflows appear at low velocity more often than SiO outflows in the same fields, and that HC3N can additionally pick up fast, collimated jet-like flows in a subset of sources.","feed_headline":"HC3N catches the slow outflow gas other tracers miss","feed_subtitle":"Survey of 146 massive clumps finds 45 HC3N outflows, including slow lobes that SiO misses.","key_machinery":"The machinery has three parts. First, the velocity range of the ambient cloud is defined by a single Gaussian fit to the optically thin H13CO+ (J=1–0) line toward each dust core; HC3N emission outside $V_{\\rm LSR} \\pm {\\rm FWHM}$ is treated as outflow wing emission. Second, outflow candidates are screened through moment-0/1/2 maps and position–velocity diagrams that must show elongated structure and a Hubble-law wedge. Third, HC3N column densities are computed pixel-by-pixel in LTE with $T_{\\rm ex}=50$ K and a fixed HC3N/H2 abundance of $5\\times10^{-9}$, from which masses, momenta, energies, dynamical times, and mechanical luminosities are derived for each lobe. The comparison with SiO outflows relies on the same ATOMS fields and includes Kolmogorov–Smirnov tests on lobe position angle, extent, and terminal velocity.","core_discovery":"The central claim is that HC3N (J=11–10) emission can be used as a reliable complementary tracer of protostellar outflows in high-mass star-forming regions, with an advantage for the low-velocity components that carry much of the outflow mass and energy. In the ALMA-ATOMS sample of 146 clumps, the authors identify 45 outflow candidates (31 confirmed, 14 probable), derive their physical parameters from LTE, optically thin column densities, and compare them with SiO (J=2–1) outflows. The HC3N outflow lobes show nearly the same orientation, extent, and terminal-velocity distribution as SiO lobes, but a slightly larger fraction of HC3N lobes is found at low velocities; ten sources even show separate fast narrow and slow wide components in position-velocity diagrams. The derived outflow mass, momentum, and energy correlate moderately (Spearman ρ ≈ 0.4–0.6) with clump mass and bolometric luminosity, but show no correlation with the $L_{\\rm bol}/M_{\\rm clump}$ evolutionary-stage proxy. The authors conclude that HC3N can detect outflows at all evolutionary stages and is an unbiased tracer that complements SiO and other traditional outflow tracers.","pith_inferences":["If HC3N indeed traces slow outflow components that SiO misses, combining the two tracers should raise the total measured outflow mass and momentum per source compared to SiO-only estimates; this is testable by re-deriving outflow parameters with combined masks.","The lack of correlation with evolutionary stage may reflect that HC3N is produced in shocks from grain-surface sputtering across a broad range of conditions; a chemical model prediction would be that HC3N abundance jumps by an order of magnitude in C-shocks regardless of protostellar age.","Because HC3N is rarer and more confined than CO, HC3N outflow surveys at moderate resolution could complement CO-based surveys in crowded or heavily contaminated regions, provided sensitivity is sufficient (the paper notes its own detection rate is sensitivity-limited).","The ten sources with both fast narrow and slow wide components suggest a two-wind geometry; high-angular-resolution follow-up of these ten could test whether the two components trace the jet and the entrained outflow as distinct kinematic structures."],"forward_implications":["Outflow surveys that include HC3N (J=11–10) can recover a larger share of the slow, massive gas that dominates outflow energetics.","HC3N can be used alongside SiO to separate slow wide-angle outflow components from fast collimated jets in the same source.","Because HC3N outflow detection does not correlate with clump evolutionary stage, surveys can use it to find outflows across a wide range of massive clumps without preselection.","Outflow mass, momentum, and energy scale with clump mass and luminosity, so HC3N-derived parameters can feed empirical calibrations of outflow feedback in high-mass star formation.","The 45-outflow catalog provides a new sample for follow-up in other tracers and for estimating outflow-driven feedback in the ATOMS clumps."],"supporting_citations":[{"why":"First observation of HC3N abundance enhancement in shocked regions, motivating HC3N as an outflow tracer.","marker":"Bachiller & Gutiérrez 1997"},{"why":"Proposed the shock chemistry mechanism (grain sputtering plus C2H2 + CN reaction) for HC3N enhancement and gave an abundance estimate adopted here.","marker":"Mendoza et al. 2018"},{"why":"Previous detections of HC3N outflows in massive clouds of the Central Molecular Zone, providing precedent and a comparison for detection rates.","marker":"Lu et al. 2021"},{"why":"Comparison of 13CO, SiO, and HC3N line widths and spatial distributions in one massive young stellar object, showing similar spatial structure.","marker":"Zinchenko et al. 2021"},{"why":"Detected extended HC3N emission wings toward star-forming cores, interpreted as outflowing gas, and provided an abundance value.","marker":"Taniguchi et al. 2018b"},{"why":"Found a strong correlation between SiO and HC3N column densities, supporting a common shocked-gas origin.","marker":"Wang et al. 2022"},{"why":"SiO outflow parameters in massive protoclusters used as the comparison benchmark for HC3N outflow parameter ranges.","marker":"Towner et al. 2024"},{"why":"Defines the ATOMS survey sample, clump mass, and luminosity measurements used for all correlations.","marker":"Liu et al. 2020"},{"why":"Provides the dust-core catalog and hot-core/ultra-compact H II region classifications used to identify outflow hosts and evolutionary stages.","marker":"Liu et al. 2021"}],"fun_headline_variants":["HC3N reveals slow outflow gas that SiO misses","45 HC3N outflows found in 146 massive star-forming clumps","HC3N: unbiased tracer of protostellar outflows at all stages","Slow outflows: HC3N outperforms SiO in survey of 146 clumps","HC3N catches low-velocity outflow lobes other tracers miss"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The cloud velocity window is fixed by a single Gaussian fit to H13CO+ (J=1–0), and all HC3N inside VLSR ± FWHM is excluded as cloud emission; if H13CO+ itself carries outflow wings, the window is too wide and the very slow outflow gas that the paper claims to detect is the first thing cut away.","fun_headline_variants_meta":{"raw":{"variants":["HC3N reveals slow outflow gas that SiO misses","45 HC3N outflows found in 146 massive star-forming clumps","HC3N: unbiased tracer of protostellar outflows at all stages","Slow outflows: HC3N outperforms SiO in survey of 146 clumps","HC3N catches low-velocity outflow lobes other tracers miss"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000639,"raw_usage":{"total_tokens":3038,"prompt_tokens":1137,"completion_tokens":1901,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":753,"completion_tokens_details":{"reasoning_tokens":1805}},"tokens_in":753,"tokens_out":1901,"duration_ms":13026,"temperature":1.0,"reasoning_tokens":1805,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:35:17.438345+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the 45 HC3N outflow hosts and inspect the H13CO+ (J=1–0) spectra for non-Gaussian wings or asymmetries; if any host shows outflow wings in H13CO+ within the adopted VLSR ± FWHM, then the claimed low-velocity HC3N component is partly an artifact of choosing too wide a cloud window.","supporting_citations":[{"cited_title":"J., et al","cited_arxiv_id":null,"evidence_quote":"Provides the dust-core catalog and hot-core/ultra-compact H II region classifications used to identify outflow hosts and evolutionary stages."}],"review_version":1}