{"id":"1ade4c71-a7c7-4fe5-9c37-686024d9c0d5","arxiv_id":"1908.05652","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"ALMA plus IRAM-30m observations of SiO, SiS, and CS toward IRC+10216 show clumpy shells and molecular abundance profiles that support episodic mass loss over the past few hundred years.","lead":"This paper maps three molecules (SiO, SiS, and CS) around the nearby evolved star IRC+10216 using ALMA and IRAM-30m telescopes, revealing clumpy shells and measuring how far each molecule extends before being destroyed. The observations indicate the star's mass loss has been episodic, with dense shells ejected every few hundred years.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The episodic-mass-loss conclusion depends on the azimuthal-averaged brightness bumps being over-dense shells; the radiative-transfer analysis itself is degenerate between density and abundance enhancements, and the chemical model cannot rule out abundance variations.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the shells could be abundance variations rather than density enhancements. The paper itself acknowledges this degeneracy in Section 4 and relies on dust-shell coincidence to break it. My independent check confirms that no quantitative test in the paper actually discriminates between the two scenarios. The chemical model does not help because it is incomplete and inherits several unverified parameters. Thus the central claim of episodic mass loss is plausible but not proven; the word 'proves' in the abstract is too strong. This does not change the overall verdict from CONDITIONAL, since the observational data are valuable and the episodic interpretation is well-supported by prior dust and CO work. The conditional accept with the stipulation that the density interpretation be treated as tentative seems appropriate. I agree with the reader's assessment.","tokens_in":22163,"tokens_out":771,"duration_ms":9333,"concrete_test":"Recompute the azimuthal-averaged brightness profiles of SiO J=2-1, SiS J=5-4 and 6-5, and CS J=2-1 using radiative-transfer models with (a) a constant fractional abundance and a gas density that is enhanced by a factor of 2-3 in a thin shell at 420 R*, and (b) a constant gas density and a fractional abundance enhanced by the same factor at the same radius. If scenario (b) fits the observed line ratios (e.g., SiS 5-4/6-5) as well as scenario (a), then the density interpretation is not uniquely constrained by the present single-line (or two-line) data. A stronger test would be to observe a second rotational line of SiO and CS at comparable resolution (e.g., SiO J=3-2 and CS J=3-2) and check whether the shell-to-inter-shell brightness ratio changes in a way that only density enhancement can reproduce.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that IRC+10216 has undergone episodic mass loss rests on interpreting the enhanced emission shells (e.g., the common shell near 420 R*) as over-dense gas shells. The paper's own Section 4 states that a brightness enhancement can be produced either by increasing the fractional abundance locally or by increasing the total gas density, and that 'we cannot conclude which of the two scenarios is more realistic' until chemical predictions are tested. The chemical model (Section 5) fails to reproduce any abundance oscillations, but this failure is attributed to model simplicity (no gas-dust interaction, unknown SiS photodissociation rate, ad hoc N_H/A_v). Therefore, the argument that the shells are density rather than chemical features relies mainly on spatial coincidence with dust arcs (Mauron & Huggins 1999) and on the absence of an obvious physical process at ~420 R* that would change chemistry. This is an external, morphological correlation, not an internal proof. If the brightness bumps were caused by abundance variations tied to, e.g., selective photochemistry or grain-surface processes at certain radii, the episodic mass-loss conclusion would not follow from these data alone. The abstract's wording 'proves' overstates the logical support. This is a genuine load-bearing weakness, but it does not invalidate the strong observational evidence for shells; it only weakens the density interpretation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents ALMA Band 3 observations at sub-arcsecond resolution combined with IRAM-30m OTF data of several lines of SiO, SiS, and CS (including isotopologues) towards the carbon-rich AGB star IRC+10216. The authors map the spatial distribution, identify arcs and shells in the azimuthal averages, and use an LVG radiative transfer code to derive radial fractional abundance profiles. They then run a chemical model to predict photo-dissociation radii. The main scientific conclusion is that the spatial distribution of the gas proves episodic, variable mass loss on timescales of hundreds of years, with the shell near 420 R* having a kinematic age of ~365 years.","tokens_in":22424,"tokens_out":4492,"duration_ms":39817,"significance":"If confirmed, the result that the observed shells trace over-dense gas episodes would be an important constraint on mass-loss mechanisms in AGB stars. The paper's strengths include the combination of ALMA and IRAM-30m data to recover short-spacing flux, the detection of multiple isotopologues, the careful use of updated collisional rate coefficients, tests of IR pumping, and the public availability of reduced data cubes (stated in the footnote). The spatial coincidence of gas and dust shells (Mauron & Huggins 1999) is a valuable morphological comparison. The paper is honest in acknowledging degeneracies, but the central claim as stated in the abstract overreaches the evidence.","major_comments":[{"comment":"The episodic mass-loss conclusion depends on the brightness-enhanced shells being over-dense. Section 4 explicitly states that the same emission enhancement can be produced by increasing the total gas density or the local fractional abundance, and that \"we cannot conclude which of the two scenarios is more realistic\" until chemical predictions are tested. The chemical model in Section 5 does not reproduce the abundance oscillations, but the authors attribute this to model simplicity (missing gas-dust interaction, uncertain N_H/A_v, unknown SiS photo-dissociation rate) rather than treating it as evidence against abundance variations. The density interpretation therefore rests on morphological coincidence with dust arcs (Mauron & Huggins 1999) and the absence of an obvious chemical mechanism at ~420 R*. This is an external correlation, not an internal proof, and the abstract's \"proves\" and Section 7's \"pose strong evidences\" overstate the support. Please reframe the conclusion as \"consistent with\" episodic mass loss and discuss how future multi-line observations or gas-grain chemical models could break the degeneracy.","section":"Section 4, Section 6.1"},{"comment":"The chemical model's parent abundances for CS, SiO, and SiS are taken from the radiative transfer fits to the same observations (Table 3, flagged \"This work\"), and the SiS photo-dissociation rate is assumed equal to that of SiO (Section 5). The resulting photo-dissociation radius ordering (SiS < SiO < CS) is therefore at least partly inherited from fitted inputs and an assumed rate, not an independent prediction. The conclusions (Section 7) present the chemical model as \"able to predict a photo-dissociation radii spatial sequence that is consistent with the observations\", which overstates the model's independence. I recommend re-running the chemical model with literature-based parent abundances and a range of SiS photo-dissociation rates to assess robustness, and clearly separating the fitted inputs from the predicted outputs in the presentation.","section":"Section 5, Table 3, Section 7"},{"comment":"The abundance derivation for SiO and CS relies on a single rotational line each, and on two lines for SiS, under the LVG approximation with fixed mass-loss rate, temperature profile, and gas-to-dust ratio taken from the literature. As the paper itself notes, \"a perfect fitting could be achieved by modifying the density, the input abundances, and the excitation conditions as a function of radius in different combinations\". This degeneracy means the reported abundances (f(SiO)~1e-7, f(SiS)~1e-6, f(CS)~1e-6) and especially the local abundance bumps are not uniquely constrained. Adding an exploration of the density/excitation parameter space or quantitative uncertainty estimates would strengthen the presentation of the radial abundance profiles.","section":"Section 4"}],"minor_comments":[{"comment":"The word \"subarsecond\" appears in the abstract and should read \"subarcsecond\".","section":"Abstract and Section 1"},{"comment":"The adjective \"inhomogeneus\" appears twice and should read \"inhomogeneous\".","section":"Sections 3.1 and 3.2"},{"comment":"In the concluding paragraph, \"SiO photo-dissociates slightly farther away but very close to SiO\" is tautological; the intended meaning is presumably \"very close to SiS\". Please correct.","section":"Section 7"},{"comment":"The statement that N_H/A_v is \"1.5 times lower than the classical value\" should state explicitly that this lower value is an ad hoc adjustment chosen to reproduce earlier observations, so that its impact on the predicted photo-dissociation radii is transparent.","section":"Section 5"},{"comment":"The text says the chemical model predicts fall-off distances ~1.5 times larger than observed; this is a non-negligible discrepancy that sits somewhat uncomfortably with the abstract's \"compatible\" characterization. A quantitative statement of the comparison, including calibration and model uncertainties, would help.","section":"Section 6.2"},{"comment":"The figure marks a \"5-sigma level detection limit\" in each panel, but the text does not explain how this limit was estimated from the azimuthal averages. A sentence in Section 3 describing the procedure would be helpful.","section":"Figure 7"}],"recommendation":"major_revision","confidential_remarks":"The observational dataset and imaging are strong, and the paper is commendably transparent about its limitations. However, the central claim of episodic mass loss is presented too strongly relative to the density/abundance degeneracy, and the chemical-model agreement is partly circular because the parent abundances and the SiS photo-dissociation rate are adopted from the same radiative-transfer fits. The requested robustness tests and rephrased conclusions are within the scope of a revision; this is not a rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the data: these are the best maps yet of SiO, SiS, and CS in IRC+10216, ALMA at ~0.6 arcsec with IRAM-30m short spacings merged in, so the large-scale flux is actually recovered. The shell radii and abundance profiles from the LVG fits are a solid step up from the older PdBI/BIMA/SMA work. I trust the observational sections; the careful calibration, the uv-plane joining, and the comparison of cleaning methods show experience. The dust/gas shell correlation in Figure 10 is convincing, and checking the 16.5-year expansion shift between the new gas maps and Mauron & Huggins optical data is exactly the right thing to do.\n\nThe weak link is not the data but the inference chain from brightness bumps to episodic mass loss. The paper says itself in Section 4 that an enhancement can be fit by raising the fractional abundance or the gas density, and that they cannot decide between the two without chemistry. The chemical model then inherits the fitted CS, SiO, and SiS abundances as parent values and assumes the SiS photodissociation rate equals that of SiO, so the predicted photodissociation ordering is partly circular. I also note the paper is honest that the chemical model fails to reproduce the observed abundance oscillations. That means the density-shell interpretation rests mainly on the morphological match to dust arcs in scattered light plus the absence of an obvious chemical trigger at ~420 R*. That is decent evidence, but it is not proof. Section 7's wording, 'strong evidences', is right; the abstract's 'proves' overstates.\n\nThis is not a fatal flaw. The shells appear at near-coincident radii in three molecules and line up with independent dust structures, so over-dense gas from past enhanced mass loss is the simplest reading. The abundance values are model-dependent and lack error bars, and the SiS J=5-4 maser complicates the two-line fit, but the ordering of photodissociation radii (SiS < SiO < CS) is consistent with prior work and with the adopted rates.\n\nWho gets value: AGB/CSE modelers, astrochemists, and observers wanting reference maps of the archetypal carbon-rich envelope. It deserves a serious referee. I would recommend acceptance with minor revision: soften 'proves', state the density/abundance degeneracy explicitly in the abstract, and give realistic error bars on the abundance profiles if feasible. I'd bring it to reading group.","headline":"Strong new maps of SiO/SiS/CS in IRC+10216 that will become a reference; the episodic mass-loss conclusion is probably right, but the abstract's 'proves' outruns the model evidence.","tokens_in":23070,"tokens_out":2506,"would_cite":true,"duration_ms":25680,"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":"High-resolution millimetre maps of SiO, SiS, and CS around the carbon star IRC+10216 reveal concentric gas shells that the paper argues are direct evidence of episodic mass loss on timescales of hundreds of years.","keywords":["IRC+10216","AGB circumstellar envelopes","episodic mass loss","SiO","SiS","CS","photo-dissociation","millimetre astronomy"],"falsifier":"Map low-J CO or another total-gas tracer at the same angular resolution as these data and check whether the brightness bumps at roughly 2.8, 6, 8.2, and 10.6 arcseconds appear at the same radii. If a total-gas tracer shows no corresponding enhancement, the shells are abundance or excitation features and the episodic mass-loss conclusion does not follow. A second check would be to resolve one shell in at least two rotational lines of the same molecule with different critical densities: if the bump is purely excitation, its apparent radius should shift between lines.","tokens_in":21963,"feed_emoji":"⭐","tokens_out":9197,"duration_ms":83320,"temperature":0.7,"pith_summary":"This paper maps three gas-phase molecules, SiO, SiS, and CS, across the circumstellar envelope of the carbon star IRC+10216 at sub-arcsecond resolution and combines the interferometric maps with single-dish data to recover the large-scale emission. The maps reveal a series of bright arcs and roughly concentric shells in all three species, with shells appearing at similar radii in different molecules, most notably near 420 stellar radii. The paper argues that these features are over-dense shells produced by past episodes of enhanced mass loss, not structure from a steady wind, and estimates a kinematic age of about 365 years for the shell at 420 $R_*$. If correct, the mass-loss history of this well-studied star is episodic on timescales of centuries, and models of such envelopes must incorporate time-variable mass loss. The paper also derives mean fractional abundances ($f(\\mathrm{SiO})\\sim 10^{-7}$, $f(\\mathrm{SiS})\\sim 10^{-6}$, $f(\\mathrm{CS})\\sim 10^{-6}$) and a photo-dissociation sequence with SiS destroyed closest to the star and CS surviving farthest out.","feed_headline":"Gas shells show IRC+10216's mass loss comes in bursts","feed_subtitle":"Three molecules trace concentric shells, one dating back about 365 years, pointing to episodic stellar winds.","key_machinery":"The central device is the azimuthally averaged radial brightness profile of each line in the channel at the source's systemic velocity, the channel that selects gas moving in the plane of the sky. Bumps in these profiles define the shells, and they are compared with predictions of a spherically symmetric large-velocity-gradient radiative-transfer model whose input is a radial fractional-abundance profile per molecule. The model reproduces the average radial decline and the order of photo-dissociation fall-off (SiS first, then SiO, then CS), but it cannot by itself distinguish an over-density from an abundance enhancement; the paper therefore leans on the spatial coincidence of the gas shells with dust arcs seen in earlier scattered-light images to tie the bumps to density. A complementary chemical model with roughly 8,000 reactions predicts the photo-dissociation sequence but not the shell oscillations, which the authors attribute to missing gas-dust surface chemistry.","core_discovery":"The central claim is that the spatial distribution of the molecular gas around IRC+10216 is inconsistent with a smooth, constant mass-loss wind and instead shows the imprint of episodic mass loss. The evidence is the set of brightness-enhanced shells seen in the azimuthally averaged emission of SiO, SiS, and CS at the systemic velocity, with several shared radii of roughly 150, 315, 420, and 560 $R_*$. Because the shells appear in several molecules at the same locations and line up with dust arcs seen in scattered light, the paper interprets them as over-dense shells formed by enhanced mass-loss events. Under a constant expansion speed of 14.5 km s$^{-1}$, the shell at about 420 $R_*$ has a kinematic age of about 365 years. The paper also measures the radial extent of each species, with CS reaching about 1000 $R_*$, and uses radiative-transfer and chemical models to infer that the envelope is clumpy enough that interstellar ultraviolet radiation can penetrate deeper than a homogeneous wind would allow, shifting photo-dissociation inward.","pith_inferences":["If the shells are over-dense, a consistent episodic history could be reconstructed from the full set of shell radii and the expansion speed; the paper does not attempt such a reconstruction, but the data would support dating individual mass-loss surges.","A decisive test would be to map a tracer of total gas column density, such as low-J CO, at comparable resolution: if the same shell radii appear in CO, the density interpretation is confirmed, whereas if only the parent molecules show bumps, chemistry or excitation is the cause.","The inferred clumpy envelope has a general implication: photo-dissociation calculations for AGB envelopes that assume smooth density will overestimate molecular survival radii whenever mass loss is episodic.","Because the radiative-transfer model cannot distinguish density from abundance enhancements, observations of multiple rotational lines of the same molecules at the same positions would break the degeneracy and strengthen or weaken the episodic conclusion."],"forward_implications":["The mass-loss rate of IRC+10216 has not been constant over the past few centuries; reproducing the shells requires episodes of enhanced mass loss separated by timescales of hundreds of years.","The shared shell radii across SiO, SiS, and CS imply that the enhanced episodes affected the envelope roughly isotropically, rather than being limited to one molecule or one direction.","The photo-dissociation order SiS then SiO then CS means the three molecules probe different depths of the envelope, with CS emission tracing the outer wind out to about 1000 $R_*$.","Clumpy, low-density regions between shells allow interstellar ultraviolet radiation to penetrate deeper, so photo-dissociation radii inferred from emission sizes are smaller than a smooth-wind model would predict.","Radiative-transfer fits that reproduce the shells require local abundance enhancements of a factor of 2-3, which the paper notes could equally be read as local density enhancements of similar size."],"supporting_citations":[{"why":"Provides the dust-arc images in scattered light whose spatial coincidence with the gas shells supports the over-density interpretation.","marker":"Mauron & Huggins (1999)"},{"why":"Supplies the radiative-transfer code, the physical envelope model, and earlier abundance determinations that the new fits revise.","marker":"Agúndez et al. (2012)"},{"why":"Supplies the updated mass-loss rate and kinetic-temperature profile adopted as fixed parameters in the model.","marker":"Guélin et al. (2018)"},{"why":"Determines the systemic and terminal expansion velocities used to convert observed radii into kinematic shell ages.","marker":"Cernicharo et al. (2000)"},{"why":"Gives earlier interferometric size measurements of the same molecules against which the new, short-spacing-corrected maps are compared.","marker":"Bieging & Tafalla (1993)"},{"why":"Provides the SiO photo-dissociation rate used in the chemical model.","marker":"Heays et al. (2017)"},{"why":"Provides the CS photo-dissociation rate used in the chemical model.","marker":"Pattillo et al. (2018)"},{"why":"Updates the SiS destruction reaction with atomic oxygen and the Si + SO formation reaction in the chemical network.","marker":"Zanchet et al. (2018)"}],"fun_headline_variants":["Envelope shells reveal episodic mass loss in IRC+10216","Molecule shells time IRC+10216's bursty stellar wind","ALMA maps shells that expose star's erratic mass loss","IRC+10216's gas shells betray centuries-long mass-loss cycles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The bright shells in the maps are really denser gas, not just regions where chemistry or excitation makes the molecules shine brighter; the models cannot tell these apart, so the episodic mass-loss conclusion rests on that identification.","fun_headline_variants_meta":{"raw":{"variants":["Envelope shells reveal episodic mass loss in IRC+10216","Molecule shells time IRC+10216's bursty stellar wind","ALMA maps shells that expose star's erratic mass loss","IRC+10216's gas shells betray centuries-long mass-loss cycles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"total_tokens":1429,"prompt_tokens":1120,"completion_tokens":309,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":736,"completion_tokens_details":{"reasoning_tokens":235}},"tokens_in":736,"tokens_out":309,"duration_ms":3490,"temperature":1.0,"reasoning_tokens":235,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:07:33.045973+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Map low-J CO or another total-gas tracer at the same angular resolution as these data and check whether the brightness bumps at roughly 2.8, 6, 8.2, and 10.6 arcseconds appear at the same radii. If a total-gas tracer shows no corresponding enhancement, the shells are abundance or excitation features and the episodic mass-loss conclusion does not follow. A second check would be to resolve one shell in at least two rotational lines of the same molecule with different critical densities: if the bump is purely excitation, its apparent radius should shift between lines.","supporting_citations":[{"cited_title":"J.\\ 1999, , 349, 203","cited_arxiv_id":null,"evidence_quote":"Provides the dust-arc images in scattered light whose spatial coincidence with the gas shells supports the over-density interpretation."},{"cited_title":"H., & Tafalla, M.\\ 1993, , 105, 576","cited_arxiv_id":null,"evidence_quote":"Gives earlier interferometric size measurements of the same molecules against which the new, short-spacing-corrected maps are compared."},{"cited_title":"N., Bosman, A","cited_arxiv_id":null,"evidence_quote":"Provides the SiO photo-dissociation rate used in the chemical model."},{"cited_title":"J., Cieszewski, R., Stancil, P","cited_arxiv_id":null,"evidence_quote":"Provides the CS photo-dissociation rate used in the chemical model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Updates the SiS destruction reaction with atomic oxygen and the Si + SO formation reaction in the chemical network."}],"review_version":1}