{"id":"421e91c3-4085-42ee-935f-79e2489eb7d3","arxiv_id":"2507.05051","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"2I/Borisov is confirmed as carbon-depleted and NH2-rich among measured comets, with a significant rise in NH2 production after its 2020 splitting event.","lead":"Astronomers mapped the interstellar comet 2I/Borisov for 126 days around its closest approach to the Sun using the VLT's MUSE spectrograph. The comet's dust and gas stayed smooth and symmetric, and it released a statistically significant extra burst of NH2 gas after its nucleus split in March 2020.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Post-split NH2 excess rests on an unverified Haser scalelength scaling; the raw flux rise is far smaller than the reported Q(NH2) rise.","rationale":"The reader identified the Haser-model radial-profile assumption as the weakest point, and I agree that this is the right area. My pass sharpens it: the truly load-bearing sub-assumption is the r_h^2 growth of the Haser scalelengths combined with the 5,000 km aperture, because it differentially amplifies late-epoch NH2 production rates and therefore controls both the new post-split enhancement claim and the NH2/CN classification. The paper is otherwise careful: the MUSE dataset is unique, the starkiller processing and Molecfit telluric corrections are sensible, the comparison to literature values is mostly consistent, and the authors state the unresolved Prodan et al. discrepancy and the aperture caveat openly. However, the headline results are presented as >95% significant even though the conversion from raw flux to production rate carries an unverified, strongly r_h-dependent model factor. A re-analysis under alternative scalelength prescriptions is straightforward and would settle whether the 30% raw-flux increase really implies a 125% production-rate increase. This does not require changing the reader's CONDITIONAL verdict: the classification and the post-split excess are plausible but not yet demonstrated to be model-independent, so the paper should either add the robustness check or soften the significance claim.","tokens_in":32160,"tokens_out":11482,"duration_ms":133098,"concrete_test":"Using the same MUSE datacubes, recompute Q(NH2) for all epochs under three alternative prescriptions: (i) scalelengths held fixed at their 1-au values; (ii) scalelengths scaled linearly with r_h; (iii) a vectorial steady-state model with the NH2 parent scalelength fitted to the azimuthally averaged radial profile where SNR permits. For each prescription, refit the pre-outburst trend used in Section 5.1 and test whether the 2020 Mar 16 and 19 values fall above the 95% prediction interval, and recompute the average Q(NH2)/Q(CN). If the excess persists under all three alternatives, the post-split NH2 claim is robust to this concern; if it disappears under any plausible alternative, the claim in Section 5.1 must be downgraded to a model-dependent result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is not spherical symmetry per se, but the specific Haser scalelength, outflow velocity, and r_h^2 scaling used to convert small-aperture fluxes into production rates, at exactly the epochs where the new post-split claim is made. In Table 4, the NH2 line flux rises from 0.56e-15 erg/s/cm2 on Feb 28 (r_h=2.70 au) to 0.73e-15 on Mar 19 (r_h=3.00 au), a 30% increase, whereas Q(NH2) rises from 2.4 to 5.4e24 mol/s, a 125% increase. The amplification comes from the model: the g-factor scales as r_h^-2 and the parent scalelength of NH2 is assumed to grow as r_h^2, reaching 36,900 km at 3 au, seven times the 5,000 km aperture radius. Thus the derived Q at 3 au is a steep function of the assumed scalelength growth and of the fixed v=0.5 km/s. If the true scalelength grows more slowly with heliocentric distance, or if the post-split coma is not in steady-state Haser equilibrium (an outburst is the least steady-state situation), the March Q(NH2) values are over-corrected and the >95% excess over the pre-outburst trend reported in Section 5.1 may vanish. The same model dependence affects Q(NH2)/Q(CN), and therefore the 'NH2-rich' classification, because NH2 and CN have very different scalelengths (parent lp 4.1e3 vs 1.3e4 km at 1 au in Table 3). The authors explicitly note in Section 3.4 that the radial profile cannot be verified at this aperture due to low SNR; the unresolved factor-several discrepancy with Prodan et al. (2024) in Q(NH2) is a warning that model and aperture choices matter at this level.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the full VLT/MUSE campaign on 2I/Borisov: 16 epochs from 2019 November 14 to 2020 March 19, covering dust colour and morphology, gas maps of C2, NH2, and CN, and Haser-model production rates extracted from a 5,000 km aperture. The main conclusions are that the coma was uniformly active with a persistent north-west jet and no seasonal evolution; that 2I is carbon-depleted (Q(C2)/Q(CN) ~ 0.1-0.3) and relatively NH2-rich compared with Solar System comets; and that after the March 2020 outburst/splitting event the NH2 production rate increased above the 95% prediction interval, with smaller rises in CN and C2 and a slight reddening of the dust. These results are interpreted as evidence of a volatile-rich, minimally processed interstellar comet whose nucleus is not strongly heterogeneous.","tokens_in":32594,"tokens_out":10625,"duration_ms":114910,"significance":"If the production-rate claims hold, this is the first compositional characterisation of an interstellar comet throughout its perihelion passage and a benchmark dataset for LSST-era ISO follow-up. The strengths are substantial: 16 epochs of a challenging faint target, careful handling of Galactic-plane stellar contamination with starkiller, consistency of the C2 and CN rates with independent literature measurements, standard external model parameters rather than fitted target values, and public data products and analysis scripts. The NH2-rich and post-split-increase claims are, however, not yet at the same level of robustness, because they depend on unverified Haser scalelength/velocity assumptions and on an unresolved factor-several discrepancy with Prodan et al. (2024).","major_comments":[{"comment":"The post-split NH2 excess is amplified by the assumed Haser scalelength scaling, so the central claim is not yet robust. In Table 4, the NH2 flux rises from 0.56e-15 erg/s/cm2 on 2020 Feb 28 to 0.73e-15 on 2020 Mar 19 (about 30%), while Q(NH2) rises from 2.4e24 to 5.4e24 mol/s (about 125%). The amplification comes from adopting g proportional to r_h^-2 and l_p(NH2) proportional to r_h^2, which makes l_p(NH2) = 36,900 km at r_h = 3 au, about 7.4 times the 5,000 km aperture radius. If the true l_p grows more slowly with heliocentric distance, or if the post-outburst coma is not in Haser steady state (an outburst is the least steady-state situation), the March Q(NH2) values are over-corrected and the reported >95% excess over the pre-outburst trend may not survive. The authors themselves note in §3.4 that the radial profile cannot be verified at this aperture. I request a quantitative sensitivity analysis: recompute the March epochs with alternative outflow velocities (e.g., 0.3 and 1 km/s), with l_p scalings of r_h^1 and r_h^1.5, and with a non-steady-state model, and report the resulting range of Q(NH2) and of the post-split significance.","section":"§3.4, Table 3, §5.1, Table 4"},{"comment":"The relatively NH2-rich classification is not yet robust to the unresolved factor-several discrepancy with Prodan et al. (2024). The text lists possible causes (different NH2 band, different aperture) but performs no test. Because the classification rests on Q(NH2)/Q(CN), and because NH2 and CN have very different Haser parent scalelengths (4.1e3 km versus 1.3e4 km at 1 au, Table 3), a systematic error of the size suggested by the Prodan et al. comparison could move 2I from the high-NH2 tail to the normal range in Figure 10. I ask for a direct cross-calibration: compute Q(NH2) from the same MUSE spectra using Prodan et al.'s (0,8,0) band, g-factor, scalelengths, and aperture, and compare the resulting Q(NH2)/Q(CN) with the published values.","section":"§4.4, Fig. 10"},{"comment":"The significance test underlying the claim that NH2 increased after splitting uses a weighted linear-regression extrapolation with only two post-outburst epochs (2020 Mar 16 and 19). The quoted uncertainties in Table 4 reflect observed-versus-modelled spectral scatter and do not include the Haser model parameters or the CN 1.4 missing-flux factor. I recommend reporting the pre-outburst sample size and fit parameters, and adding a sensitivity test that includes the model systematics (e.g., a Monte Carlo over g-factors, scalelengths, and velocities) before the >95% language is used.","section":"§5.1, Fig. 11"}],"minor_comments":[{"comment":"Please state unambiguously whether the extraction aperture for production rates is 5,000 km radius or diameter; §3.2.1 refers to a 5,000 km diameter aperture for dust colour, while Figure 3 says within a 5,000 km radius. Since the aperture size enters the Haser conversion, this mismatch must be resolved.","section":"§3.4, Fig. 3 caption"},{"comment":"The data DOI is currently a placeholder; the final DOI should be included at acceptance.","section":"Footnote 3"},{"comment":"The role of the PSG synthetic spectra is hard to follow; clarify that synthetic spectra were used to measure band fluxes, while production rates inferred directly from PSG were inconsistent and not used.","section":"§3.4"},{"comment":"The phrase 95% prediction confidence interval is non-standard; use prediction interval for extrapolated point comparisons.","section":"§3.2.1 and Fig. 11"},{"comment":"The column headed l is explained in the table note as item 10 but is not labelled; please define l in the header or note.","section":"Table 1"},{"comment":"Comet designations such as C/2016R2 should follow the IAU style C/2016 R2 for consistency.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and the dataset is valuable. I do not see circularity: production rates are measured, and the comparison samples come from the literature. The path to acceptance is a focused revision: a Haser-parameter sensitivity analysis for the post-split epochs, a cross-calibration against Prodan et al. (2024), and a resolution of the aperture-radius ambiguity. I would not require new observations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious referee: this is the only full-perihelion compositional record we have for an interstellar comet, 16 MUSE epochs across 126 days. The paper does three things well. It documents the reduction carefully, including the starkiller star-subtraction work that made the crowded-field post-perihelion epochs usable. It places the C2 and CN production rates in context with independent measurements, and those agree. And it makes a clear, falsifiable claim about post-splitting behavior: NH2 rose relative to the pre-outburst trend, while C2 and CN rose with larger uncertainties and the dust reddened slightly.\n\nThe soft spot is exactly where the stress-test points: the post-split NH2 excess is amplified by the Haser model. The March 19 NH2 line flux is only ~30% above Feb 28, but Q(NH2) more than doubles, because the assumed parent scalelength grows as r_h^2 and reaches ~37,000 km at 3 au, seven times the 5,000 km aperture radius. With v fixed at 0.5 km/s and the radial profile unverified at that aperture (the authors say so in Section 3.4), the magnitude of the NH2 increase is not robust. That does not kill the paper, but it means the headline post-split claim should be stated as model-dependent, and the authors should show how Q(NH2) varies under alternative scalelength growth laws or outflow velocities. The unresolved factor-several discrepancy with Prodan et al. (2024) in Q(NH2) is a warning that model choices matter at exactly this level.\n\nThe NH2-rich classification itself is less fragile, because it was already indicated in Bannister et al. (2020) and the ratio Q(NH2)/Q(CN) is high across the whole campaign, not just post-split. The carbon-depleted status is confirmed by C2/CN ratios consistent with independent work. So the central classification holds; the new temporal detail is the part that needs scrutiny.\n\nData availability is a minor issue: the DOI is promised but not yet in the manuscript. For a reference dataset, that should be resolved before publication.\n\nIn sum: this is a solid, careful observational paper with one load-bearing but addressable model-dependence issue. It deserves peer review and, after revision, publication. I would bring it to reading group and cite it for the dataset.","headline":"First full-perihelion compositional record of an interstellar comet; the post-split NH2 claim is model-dependent but the dataset deserves a serious referee.","tokens_in":33263,"tokens_out":2393,"would_cite":true,"duration_ms":23788,"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":"The 16-epoch spectroscopic campaign establishes that interstellar comet 2I/Borisov was carbon-depleted and NH2-rich, and that its early-2020 nucleus split produced a statistically significant surge in NH2 gas.","keywords":["interstellar comet","2I/Borisov","comet composition","NH2 radical","gas production rates","Haser model","carbon depletion","nucleus splitting"],"falsifier":"Re-measure the NH2, C2, and CN production rates in apertures of 5,000 km, 10,000 km, and 20,000 km on the same data: if the derived production rates change systematically with aperture beyond the stated uncertainties, the Haser profile assumption fails and the claimed NH2 enhancement after the split is not robust.","tokens_in":31939,"feed_emoji":"☄️","tokens_out":5403,"duration_ms":50640,"temperature":0.7,"pith_summary":"The paper tries to establish, from the only dataset that tracked an interstellar object's composition through an entire perihelion passage, what 2I/Borisov was made of and how it changed when its nucleus split. It reports that the comet's C2, NH2, and CN production rates declined gently toward and past perihelion, then all increased after the 2020 March outburst and splitting event, with NH2 rising above the 95% prediction interval of pre-outburst trends while C2 and CN rose with larger uncertainty and the dust reddened slightly. The result, if correct, makes 2I/Borisov a carbon-depleted but relatively NH2-rich comet compared with measured Solar System comets, and a volatile-rich, uniformly active object that experienced minimal processing before leaving its parent system. A sympathetic reader would care because this is the first direct compositional window into the interior of an object from another planetary system.","feed_headline":"Interstellar comet's split unleashed a surge of NH2 gas","feed_subtitle":"A 126-day portrait of 2I/Borisov shows a carbon-poor but NH2-rich coma that changed after its nucleus broke apart.","key_machinery":"The load-bearing machinery is the Haser-model conversion of band fluxes into molecular production rates: measured fluxes are divided by fluorescence efficiencies (g-factors) and matched to parent and daughter scalelengths at a fixed outflow velocity of 0.5 km/s, with the model's radial profile presumed to hold inside the 5,000 km aperture used to avoid stars. Complementing this, the paper builds spatially resolved gas and dust maps from an integral-field spectrograph and tests post-splitting changes against 95% prediction intervals extrapolated from pre-outburst linear regressions.","core_discovery":"On the paper's own terms, 2I/Borisov's coma was compositionally and morphologically smooth: dust emission was mostly featureless apart from a persistent jet-like structure toward the north-west, and C2, NH2, and CN were all emitted symmetrically around the photocentre. The gas production rates, derived with a Haser model inside a 5,000 km aperture, declined gently through perihelion until the March 2020 outburst and splitting of the nucleus, after which NH2 production rose significantly above the pre-outburst trend while C2 and CN also increased though with larger uncertainties, and the dust colour reddened slightly. From abundance ratios Q(C2)/Q(CN) between 0.1 and 0.3 and Q(NH2)/Q(CN) averaging 1.7, the paper concludes that 2I/Borisov belongs to the carbon-depleted class of comets and is relatively NH2-rich compared with Solar System comets, implying a volatile-rich nucleus with minimal surface processing.","pith_inferences":["One extension the authors leave implicit: if NH2's parent is NH3, the post-split NH2 surge would make ammonia a major interior ice in 2I/Borisov, and future interstellar-object spectroscopy could test whether NH3-rich interiors are common.","The tentative decline of Q(C2)/Q(CN) with heliocentric distance, combined with the NH2 and C2 decline tracking H2O rather than CO, suggests a shared volatile source; a multi-wavelength campaign on a future ISO could test whether C2 and NH2 parents are trapped in the same ice phase as water.","The dust reddening after the split, if connected to fresh larger particles, predicts that the fragment's debris should show a steeper size distribution; searching archival space-based images for the fragment's color could confirm this."],"forward_implications":["If 2I/Borisov is genuinely NH2-rich and carbon-depleted, then interstellar comets can have volatile inventories that differ from typical Solar System comets in the same diagnostic ratios that separate comet classes locally.","The significant NH2 increase after the split implies the freshly exposed interior released more NH2 than the pre-split surface, which would make 2I's outer layer slightly volatile-depleted relative to its bulk.","The smooth, symmetric gas coma and persistent dust jet with no seasonal variation indicate a uniformly active nucleus rather than discrete, rotationally modulated source regions.","The 126-day campaign demonstrates that integral-field spectroscopy can recover both gas and dust information even when an interstellar object crosses the Galactic plane, which will matter for planning observations of future interstellar objects."],"supporting_citations":[{"why":"Supplies the initial gas production rates from the opening epochs of this campaign and the first NH2-rich assessment that the full dataset extends.","marker":"Bannister et al. (2020)"},{"why":"Provides the model that converts measured band fluxes into molecular production rates through assumed parent-daughter scalelengths.","marker":"Haser (1957)"},{"why":"Defines the carbon-depleted comet class and supplies the C2 g-factor and scalelengths used in the production-rate calculation.","marker":"A'Hearn et al. (1995)"},{"why":"Reports the Hubble Space Telescope resolution of the split fragment that brackets the March 2020 outburst and anchors the post-splitting analysis.","marker":"Jewitt et al. (2020)"},{"why":"Documents the March 4-9 outburst that the paper uses to define the pre- and post-splitting epochs.","marker":"Drahus et al. (2020)"},{"why":"Provides the high CO abundance context and the declining H2O trend that the paper compares with its NH2 and C2 behavior.","marker":"Bodewits et al. (2020)"},{"why":"Supplies the high CO measurement supporting formation beyond the CO snowline, which the paper's volatile-rich conclusion extends.","marker":"Cordiner et al. (2020)"},{"why":"Provides the comparison sample of Solar System comet NH2/CN ratios and the 1.4 flux correction factor for the CN band.","marker":"Fink (2009)"},{"why":"Supplies the NH2 scalelengths and a comparison sample of NH2/CN ratios used in the NH2-rich classification.","marker":"Cochran et al. (2012)"},{"why":"Describes the star-subtraction tool that allowed the Galactic-plane observations to be included in the dust and gas maps.","marker":"Ridden-Harper et al. (2025)"}],"fun_headline_variants":["Borisov's split unleashed NH2 surge","Interstellar comet's burst flaunts NH2 richness","Carbon-poor 2I/Borisov gets NH2-rich after split","Comet 2I/Borisov: NH2 rises after nucleus break","126-day portrait: Borisov's split boosts NH2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All production rates assume the coma's gas follows a spherical Haser model with a constant outflow speed of 0.5 km/s and literature scalelengths, but the data are too faint to verify the radial profile at the 5,000 km aperture, so a mismatch would shift the derived Q(NH2), Q(C2), and Q(CN) and could change the NH2-rich classification.","fun_headline_variants_meta":{"raw":{"variants":["Borisov's split unleashed NH2 surge","Interstellar comet's burst flaunts NH2 richness","Carbon-poor 2I/Borisov gets NH2-rich after split","Comet 2I/Borisov: NH2 rises after nucleus break","126-day portrait: Borisov's split boosts NH2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000258,"raw_usage":{"total_tokens":1595,"prompt_tokens":968,"completion_tokens":627,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":584,"completion_tokens_details":{"reasoning_tokens":540}},"tokens_in":584,"tokens_out":627,"duration_ms":6138,"temperature":1.0,"reasoning_tokens":540,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:33:25.292238+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-measure the NH2, C2, and CN production rates in apertures of 5,000 km, 10,000 km, and 20,000 km on the same data: if the derived production rates change systematically with aperture beyond the stated uncertainties, the Haser profile assumption fails and the claimed NH2 enhancement after the split is not robust.","supporting_citations":[{"cited_title":"T., Deam , S","cited_arxiv_id":null,"evidence_quote":"Describes the star-subtraction tool that allowed the Galactic-plane observations to be included in the dust and gas maps."}],"review_version":1}