{"id":"94fb0a30-1f85-43b7-974c-34d21ae43f79","arxiv_id":"2507.13979","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Simultaneous VLT/MUSE maps of dust, [OI], C2, NH2, and CN in comet 67P's 2021 coma reveal that NH2 and CN are linked to dust fans, with NH2 scale lengths 1.5-1.9 times longer in one region, hinting at extended sources.","lead":"Using the MUSE spectrograph on the VLT, the authors mapped dust and four gas species in the coma of comet 67P across twelve nights around its 2021 perihelion. The maps show that NH2 and CN emission are linked to specific dust fans, suggesting some of these gases are released from dust grains far from the nucleus.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"NH2 scale-length contrast and CN–Fan B correlation are degenerate with anisotropic nucleus outgassing; extended-source claim lacks a jet-geometry null test.","rationale":"The reader's weakest assumption correctly identifies dust-continuum subtraction as a risk to the gas maps, and the 1–5% residual at the red/blue extremes is plausible. However, the CN maps used a local polynomial subtraction with a dedicated residual test (Appendix B), and NH2 is in the central spectral region where the global method is stated to be robust. The more decisive gap for the paper's central claim is interpretive: the enhanced scale-length contrast and the CN-dust spatial correlation are both consistent with a simple nucleus-source jet sharing the same active region as the dust fan. The paper does not run any jet or anisotropic-outflow model, nor does it provide a CN radial-scale-length comparison, so the extended-source interpretation is not uniquely supported. This does not invalidate the valuable morphological dataset or the cautiously worded conclusions, but it means the strongest claim as summarized by the reader overstates the evidence. The verdict should remain conditional: the paper would be strengthened by a jet-geometry null test or by softening the extended-source statements to 'not inconsistent with.' The suggested concrete test (3D jet+Haser fitting) would settle whether the NH2 scale-length contrast survives a geometric null hypothesis; if it does not, the extended-source conclusion for NH2, and analogously for CN, cannot be sustained.","tokens_in":37689,"tokens_out":6758,"duration_ms":93887,"concrete_test":"Run a 3D model combining a spherical Haser coma with a collimated jet cone in the NE quadrant (standard non-extended parent and daughter scale lengths) and fit the resulting synthetic radial profiles with the same quadrant Haser procedure used in Sec. 3.3.2. If the apparent NE parent scale length is inflated by ~1.5–1.9x, then the observed NH2 enhancement is explained by jet geometry alone, removing the need for an extended source; optionally, apply the same test to the CN jet profile if SNR permits, comparing against prompt-production models using HCN rates from Biver et al. (2023).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim (Secs. 3.3.2, 3.1.2, 4) is that enhanced NH2 scale lengths (1.5–1.9x, Table 3) and a CN structure coincident with dust Fan B indicate extended dust sources. Both diagnostics are degenerate with an alternative the paper never models: a collimated, jet-like outflow from a localized nucleus source. A jet cone produces a radial column-density profile in its quadrant that is shallower than the spherical Haser assumption (roughly ~1/rho versus ~1/rho^2), so fitting the spherically symmetric Haser model (Eq. 3) biases the fitted parent scale length upward. The NE quadrant is exactly where the large curved NH2 structure sits, so the factor 1.5–1.9 may be a geometric artifact of fitting a jet with a spherical model, not evidence of a physically extended source. For CN, Appendix B only rules out residual dust continuum; it does not rule out the simpler possibility that the same active region emits both dust (Fan B) and HCN (or another CN parent), with CN produced promptly by photolysis. HCN is known to sublimate from southern mid-latitudes (Läuter et al. 2022), and no radial profile or scale-length comparison for CN is presented to distinguish prompt from extended production. The paper's cautious wording ('consistent with', 'tentative') reduces risk, but the reader's strong claim states that NH2 and CN 'are produced from dust grains away from the nucleus,' which the presented evidence does not uniquely support.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents VLT/MUSE integral-field observations of comet 67P from May 2021 to March 2022, together with archival 2016 data, and uses them to build simultaneous maps of dust, [OI], C2, NH2, and CN. After telluric correction and dust-continuum subtraction, the authors characterize dust fans A-F, spectral slope maps, green/red [OI] ratios, and the evolving morphology of each gas species. The two principal claims are: (1) in the pre-perihelion northern coma, the fitted NH2 parent scale length in the NE quadrant is 1.5-1.9 times larger than in the other quadrants (Table 3), interpreted as possible evidence for an extended NH2 source; and (2) a collimated CN structure on 30 September 2021 is spatially correlated with dust Fan B and survives a dust-residual test (Appendix B), tentatively indicating an extended CN source carried by larger dust particles. The paper is explicitly descriptive and does not attempt detailed modeling.","tokens_in":38017,"tokens_out":17246,"duration_ms":182510,"significance":"The dataset is a valuable contribution: twelve epochs of simultaneous optical gas and dust mapping of a Rosetta target, with careful processing steps including Molecfit telluric correction, per-spaxel dust subtraction, azimuthal enhancement, and cross-checked position-angle measurements via Voigt fits and continuous wavelet transforms. The explicit dust-residual test for the 30 September CN map (Appendix B) is a good example of due diligence. If the extended-source interpretation for NH2 and CN survives the tests requested below, the result would be significant for understanding nitrogen chemistry and dust-gas coupling in 67P, and for linking ground-based and Rosetta observations. At present, however, the central inference is not uniquely supported: the NH2 scale-length contrast is degenerate with anisotropic jet geometry, and the CN correlation does not distinguish extended dust production from prompt HCN sublimation from the same active region.","major_comments":[{"comment":"Equation (3) as printed is the volume density of the daughter species, not the column density: the denominator contains rho^2, giving units of molecules cm^-3, whereas Ncoma(rho) from Eq. (2) is a column density with units of molecules cm^-2. The standard Haser column-density formula has rho to the first power (or requires a line-of-sight integration). Fitting Eq. (3) directly to Eq. (2) is dimensionally inconsistent unless the code actually uses the 1/rho form or performs the projection. Please verify the implemented formula and correct Eq. (3) accordingly; the fitted scale lengths in Table 3 depend directly on this.","section":"Section 3.3.2, Eq. (3)"},{"comment":"The claimed 1.5-1.9x enhancement of the NE NH2 parent scale length rests on fitting the spherically symmetric Haser model to quadrant profiles of a strongly anisotropic coma. A localized, jet-like outflow produces a shallower radial column-density profile in its own quadrant than the spherical average, which biases the fitted parent scale length upward. Since the NE quadrant is exactly where the large curved NH2 structure sits, the enhancement may be a geometric artifact of the model rather than evidence for a physically extended source. Please add a null test: fit a conical-jet or otherwise anisotropic Haser variant to the same quadrant profiles, or compare the NE profile against the spherical-Haser expectation after masking the structure, and show that the NE-vs-average contrast survives. Table 3 should also list the fitting uncertainties, since Figure 11 currently shows large error ranges.","section":"Section 3.3.2, Table 3"},{"comment":"All gas maps are produced by subtracting a scaled, shape-corrected version of the 2016 reference dust spectrum, with acknowledged 1-5% over- and under-subtraction at the blue and red extremes. Appendix B tests residual dust only for the 30 September CN map. The NH2 maps and scale-length fits in Section 3.3.2, as well as the C2 and [OI] morphology maps, are produced by the same pipeline but are not tested for sensitivity to residual dust. Please run an injection test: add residual-dust maps at the 1-5% level with spatial structure matched to Fans A and B, re-fit the NH2 scale lengths, and show that the NE enhancement is robust. Also consider presenting residual maps for the NH2 and C2 bands.","section":"Section 2.2, Appendix B"},{"comment":"The CN-Fan B correlation is consistent with an extended dust source, but it is also consistent with the same southern active region releasing both dust and HCN, with CN produced promptly by photolysis. The authors acknowledge this and state that they could not fit CN scale lengths, yet the abstract and conclusions say the correlation 'potentially reveals' an extended source origin. Please present either a CN radial profile or a Haser scale-length comparison for the high-SNR 30 September map, or revise the wording so that prompt HCN and extended dust sources are both described as consistent with the data.","section":"Sections 3.1.2, 4, and Appendix B"},{"comment":"The spectral-slope conversion factor in Eq. (1) appears to be too small by a factor of 10. For example, a linear flux increase of 10% across 2000 Angstrom corresponds to 5%/1000 Angstrom, but Eq. (1) with F6000 normalized to unity gives 0.5%/1000 Angstrom. Please correct the conversion factor or show the derivation. This affects all spectral slope values in Section 3.2 and Figure 9, including the spatially resolved slopes used to argue for larger dust particles in the Fan B/CN region.","section":"Section 2.2, Eq. (1)"}],"minor_comments":[{"comment":"The table reports no uncertainties for the fitted scale lengths, although Figure 11 shows large error regions; please add numerical uncertainties so that the significance of the NE-vs-average contrast can be assessed.","section":"Table 3"},{"comment":"The spectral slope maps are presented without error maps or per-pixel uncertainties; the text states that an uncertainty of 1-2%/1000 Angstrom should be assumed. Please state this in the figure caption and ideally provide error maps.","section":"Section 3.2, Figure 9"},{"comment":"The text says the campaign comprises 12 epochs, while Table 1 lists 13 observation dates including the rejected 10 January 2022 observation; please clarify the count in the text or table.","section":"Section 2.1, Table 1"},{"comment":"The frequent citation of 'Ivanova et al., in prep.' is not ideal for a journal publication; please replace it with a published reference or remove the specific quantitative comparisons until the work is available.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of Planetary and Space Science and the dataset is genuinely useful. The main concern is that the central extended-source claim is currently under-supported by the combination of a spherically symmetric Haser fit to anisotropic quadrants and a dust-subtraction pipeline whose residuals are tested only for one CN map. These issues are addressable, but they are load-bearing for the paper's headline result, so I recommend major revision rather than minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The dataset is the real product, and it deserves to be published. This is the first MUSE simultaneous mapping of dust, [OI], C2, NH2, and CN across 12 epochs of 67P's 2021 apparition, and the processing is careful: Molecfit telluric correction, iterative 2016-reference dust continuum subtraction, a dedicated gas-dust disentangling test for CN (Appendix B), and CWT/Voigt methods for dust fan position angles. The morphology is consistent with prior 2021 work (Boehnhardt et al. 2024, Biver et al. 2023, Bonev et al. 2023), and the raw data sit in the ESO archive. As a reference dataset for linking ground-based coma work to Rosetta, it has real value.\n\nThe interpretive layer is where I'd push back. The stress-test concern lands: the 1.5-1.9x NH2 scale-length enhancement in the NE quadrant comes from fitting a spherically symmetric Haser model. A collimated jet from a localized source produces a shallower radial column profile, which biases the fitted parent scale length upward, and the NE quadrant is exactly where the big curved NH2 structure sits. The paper never models or tests that geometry. Three epochs and large covariance uncertainties make the contrast weaker still. That said, the authors' language is honestly hedged — 'consistent with,' 'we can only report' — so this is a need-a-null-test-or-softer-claim problem, not a fatal one. Same for CN-Fan B: Appendix B rules out residual dust continuum, but not prompt HCN photolysis from the same active region that feeds the dust fan, and the authors themselves cite Läuter et al. (2022) for southern HCN sources.\n\nMinor soft spots: the spectral slope maps have no formal error bars (the reader is told to assume ±1-2%/1000 Å), which is thin for the 'larger redder particles at Fan B' argument; and the acknowledged 1-5% over/under-subtraction at the blue and red spectral extremes is only partly mitigated by the separate 1D polynomial treatment for C2 and CN. Circularity is a non-issue — the Haser inputs are standard literature values and the comparisons are internal to their own quadrant fits.\n\nWho gets value: coma observers, the Rosetta-link community, anyone planning MUSE comet work. Send it to peer review. A serious referee should ask for a jet-geometry null test, formal slope uncertainties, and a radial CN profile; with those, it becomes a stronger paper.","headline":"Valuable MUSE dataset of 67P's 2021 apparition, with an extended-source interpretation that is suggestive but not uniquely demonstrated — the Haser fits need a jet-geometry null test.","tokens_in":38575,"tokens_out":6337,"would_cite":true,"duration_ms":69645,"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":"VLT/MUSE maps of comet 67P show NH2 and CN are released from dust, not only the nucleus.","keywords":["comet 67P","coma morphology","NH2 radical","CN radical","extended source","Haser scale length","dust spectral slope","integral field spectroscopy"],"falsifier":"Recompute all gas maps with an independent dust model built from each epoch's own gas-free windows, then refit the Haser profiles: if the NE quadrant's 1.5-1.9x NH2 scale-length enhancement and the CN-FanB correlation disappear, the extended-source conclusion collapses.","tokens_in":37531,"feed_emoji":"☄️","tokens_out":6613,"duration_ms":69661,"temperature":0.7,"pith_summary":"This paper uses 12 epochs of VLT/MUSE integral-field spectra of comet 67P/Churyumov-Gerasimenko across its 2021 perihelion to build simultaneous maps of dust and four gas species. It argues that the radical NH2 in the northern pre-perihelion coma is partly produced from dust grains far from the nucleus, because fitted parent scale lengths in that quadrant are 1.5-1.9 times longer than elsewhere. It also reports a CN jet spatially coincident with a dust fan, verified to be free of dust contamination, implying that larger redder dust particles can carry CN parent material outward. The maps further separate the coma into an evolving water and C2 component tied to nucleus sublimation and a stable NH2 and CN component tied to seasonal and distributed sources. If correct, the work means remote optical observations of comae must account for gas released from dust, not only from the nucleus.","feed_headline":"Comet gas comes partly from dust grains, new 67P maps show","feed_subtitle":"Fitted NH2 scale lengths run 1.5-1.9x longer in one region, and a CN jet matches a dust fan with no residual dust contamination.","key_machinery":"The argument is carried by radial Haser model fitting applied to quadrant wedge intensity profiles, comparing fitted parent scale lengths against known values, with a normalized chi-squared map test that verifies the CN signal at Fan B is not residual dust continuum. The dust subtraction itself uses a shape-corrected reference dust spectrum, scaled to each epoch, and the radial profiles are built from azimuthal-median-enhanced maps that reveal the narrow gas structures.","core_discovery":"The paper's central claim is that two radicals in the coma of 67P, NH2 and CN, have spatially localized structures that cannot be explained by direct nucleus release alone: in the northern pre-perihelion coma the fitted NH2 parent scale length is 1.5-1.9 times larger than in other quadrants, and a narrow CN jet coincides with dust Fan B after a residual-dust test found no contamination. The authors interpret both as evidence that some NH2 and CN are produced from dust grains away from the nucleus, possibly through ammoniated salts and larger dust particles that preserve parent species. The paper also shows that H2O ([OI]) and C2 track the subsolar point and nucleus sublimation while NH2 and CN remain morphologically stable, defining two evolutionary regimes in the coma.","pith_inferences":["If the dust-borne NH2 comes from ammonium salts, the NH2 enhancement should track the largest, reddest dust particles; a direct pixel-by-pixel correlation of NH2 maps with spectral-slope maps over the same epochs would test that.","The same 1.5-1.9x scale-length signature might appear in other Jupiter-family comets with similar viewing geometry; a survey of one or two apparitions would show whether the mechanism is common.","The CN-FanB association predicts enhanced HCN or other CN parents along the same position angle in radio maps, which could be searched for with submillimetre interferometry.","If the reference-dust subtraction method is the limiting step, validating it against a synthetic coma model with a known gas distribution would put bounds on the 1-5% residual problem."],"forward_implications":["Total NH2 and CN production rates in 67P cannot be derived from nucleus outgassing alone; a dust-borne fraction must be included.","Single Haser scale lengths fitted to whole-coma radial profiles will bias production rates whenever a localized extended source is present.","Spatial coincidence of gas jets with dust fans, plus redder spectral slopes, gives a remote-sensing diagnostic for dust particles that preserve volatile parent species.","The stability of NH2 and CN structures across epochs suggests seasonal illumination, not short-term activity, controls their morphology.","Ground-based maps like these can bridge to in-situ measurements by locating distributed sources that close-range instruments sample only partially."],"supporting_citations":[{"why":"Supplies the 2016 reference dust spectrum and the continuum-subtraction method used to extract all gas maps.","marker":"Opitom et al. (2020)"},{"why":"Supplies the MUSE-tested Haser model parameters (vgas, parent and daughter scale lengths) used in the NH2 fits.","marker":"Bannister et al. (2020)"},{"why":"Supplies NH2 production-rate estimates used as starting values for the Haser fitting.","marker":"Lara et al. (2011)"},{"why":"Provides the NH2 fluorescence efficiency and the scale-length framework used to interpret the fits.","marker":"A'Hearn et al. (1995)"},{"why":"Establishes NH2 as a product of NH3 and provides the scale-length context for the extended-source interpretation.","marker":"Fink et al. (1991)"},{"why":"Reports in-situ detection of CN and evidence for a distributed source, the key comparative evidence for the CN extended-source claim.","marker":"Hanni et al. (2020)"},{"why":"Proposes ammonium salts as a nitrogen reservoir, the mechanism proposed for the NH2 extended source.","marker":"Altwegg et al. (2020)"},{"why":"Provides complementary evidence for ammonium salts on the nucleus, supporting the NH2 extended-source mechanism.","marker":"Poch et al. (2020)"},{"why":"Established the southern CN enhancement after equinox, the seasonal framework used to interpret the CN maps.","marker":"Opitom et al. (2017)"},{"why":"Independently characterized dust fans in the 2021 apparition, used to align Fans B through F and support the CN-FanB correlation.","marker":"Boehnhardt et al. (2024)"}],"fun_headline_variants":["67P's NH2 and CN come partly from dust, not just nucleus","Dust grains feed 67P's NH2 and CN coma","67P maps: NH2 and CN trace dust fans, hinting at extended sources","Comet 67P: NH2 and CN born on dust, not just nucleus","Dust holds parent molecules for 67P's NH2 and CN"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The dust subtraction that produces every gas map assumes that the 2016 reference dust spectrum, rescaled to 2021 with a polynomial fit, matches the 2021 dust coma well enough that whatever remains is real gas; if the subtracted dust leaves spatially varying residuals (the paper allows 1-5% at the blue and red ends), the NH2 scale-length contrast and CN jet could appear where none exist.","fun_headline_variants_meta":{"raw":{"variants":["67P's NH2 and CN come partly from dust, not just nucleus","Dust grains feed 67P's NH2 and CN coma","67P maps: NH2 and CN trace dust fans, hinting at extended sources","Comet 67P: NH2 and CN born on dust, not just nucleus","Dust holds parent molecules for 67P's NH2 and CN"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000476,"raw_usage":{"total_tokens":2402,"prompt_tokens":1028,"completion_tokens":1374,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":644,"completion_tokens_details":{"reasoning_tokens":1272}},"tokens_in":644,"tokens_out":1374,"duration_ms":10531,"temperature":1.0,"reasoning_tokens":1272,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:11:11.107771+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute all gas maps with an independent dust model built from each epoch's own gas-free windows, then refit the Haser profiles: if the NE quadrant's 1.5-1.9x NH2 scale-length enhancement and the CN-FanB correlation disappear, the extended-source conclusion collapses.","supporting_citations":[{"cited_title":"MUSE observations of comet 67P/Churyumov-Gerasimenko: A reference for future comet observations with MUSE","cited_arxiv_id":"2010.05064","evidence_quote":"Supplies the 2016 reference dust spectrum and the continuum-subtraction method used to extract all gas maps."},{"cited_title":", author Lin , Z.Y","cited_arxiv_id":null,"evidence_quote":"Supplies NH2 production-rate estimates used as starting values for the Haser fitting."},{"cited_title":", author Combi , M.R","cited_arxiv_id":null,"evidence_quote":"Establishes NH2 as a product of NH3 and provides the scale-length context for the extended-source interpretation."},{"cited_title":"Evidence of ammonium salts in comet 67P as explanation for the nitrogen depletion in cometary comae","cited_arxiv_id":"1911.13005","evidence_quote":"Proposes ammonium salts as a nitrogen reservoir, the mechanism proposed for the NH2 extended source."},{"cited_title":"Ammonium salts are a reservoir of nitrogen on a cometary nucleus and possibly on some asteroids","cited_arxiv_id":"2003.06034","evidence_quote":"Provides complementary evidence for ammonium salts on the nucleus, supporting the NH2 extended-source mechanism."},{"cited_title":", author Snodgrass , C","cited_arxiv_id":null,"evidence_quote":"Established the southern CN enhancement after equinox, the seasonal framework used to interpret the CN maps."}],"review_version":1}