{"id":"0b1cc9b7-7f65-47da-93e0-d3b97f92270b","arxiv_id":"2505.21625","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Only three sodium events appear in 17 years of HARPS spectra of Beta Pictoris, all linked to deep calcium absorption, and two 2019 exocomets show non-Keplerian acceleration interpreted as destruction.","lead":"After correcting thousands of archived spectra for Earth's atmosphere, this study finds that sodium vapor from exocomets around Beta Pictoris is extremely rare and appears only during the strongest calcium events. It also reports two comets in 2019 that suddenly accelerated away from normal orbital motion, which the authors interpret as the comets breaking apart near the star.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fragmentation claim rests on single-Gaussian centroid fits; an unresolved second absorption component could produce the apparent non-Keplerian acceleration and line-width jumps without dynamical disruption.","rationale":"Read in good faith, the paper does what it claims: it uses a large homogeneous HARPS sample, applies explicit molecfit telluric correction, and documents Bayesian fits with stated priors. The Na I detections are few but individually credible, and the disk-variability caution in Sec. 3.3 is honest. The fragmentation claim is the most speculative and the most novel part of the paper, and it is explicitly labelled a hypothesis. Its evidential basis is only two six-hour sequences, one of which is blended. The centroid-fitting ambiguity is concrete and testable, so the appropriate outcome remains CONDITIONAL/UNCHANGED rather than ACCEPT or REJECT. The reader's reference-spectrum concern is real but applies mainly to the long-lived -5 km/s feature and to absolute depth calibration; it does not directly threaten the -30 km/s 2019 kinematics. That is why I mark agreement as partial.","tokens_in":18954,"tokens_out":6173,"duration_ms":71940,"concrete_test":"Re-fit the full 2019-12-11 time series with an additional Gaussian component whose priors mirror the second component used on night 2 (Table A.2, tau4/lambda4), and recompute the main component's centroid acceleration curve. If, after including this component, the main feature is consistent with a constant Keplerian acceleration over the full six hours, the reported non-Keplerian departure—and the fragmentation inference built on it—is not supported; if the non-Keplerian shape persists, the concern is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2's central inference—that the 2019 blue-shifted features depart from Keplerian motion because the nuclei fragmented—depends on the fitted centroid v0 of a single exocomet Gaussian component tracing the absorbing cloud's radial velocity. The most anomalous part of the time series is precisely where that assumption is least secure: on night 2 the authors introduce a second static component near -39 km/s (Fig. 9), omit unstable fits at the end of the night, and report a simultaneous, dramatic broadening of the line (Fig. 10). A second unresolved or emerging component can shift and broaden the centroid of a single-Gaussian fit without any dynamical acceleration, and the paper argues against tail-egress but does not test this alternative for night 1. The reference-spectrum instability identified by the reader is credible for the year-long -5 km/s feature and for disk-line variability (Sec. 3.3, Fig. A.5), but the 2019 features are well separated from the disk line; the unresolved-blend degeneracy is the more load-bearing risk for the paper's most novel conclusion. The year inconsistency (2019 in the abstract, 2018 in Sec. 4) is a smaller but related presentational problem.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a systematic analysis of 9071 archival HARPS spectra of Beta Pictoris (198 nights, 2003-2020) aimed at finding new exocomet phenomenology in the Ca II H&K lines and the Na I D doublet. The authors apply molecfit-based telluric correction, construct an 80th-percentile reference spectrum with Tikhonov smoothing, and fit the exocomet and disk absorption with a Bayesian multi-component Gaussian model (JAX/NumPyro) with explicit, tabulated priors. The main results are: (i) exocometary Na I absorption is generally absent, with only two ~2%-deep events (March 2008 and March 2009) and a weaker (~1%) feature persisting over 13 nights in February 2004, all coincident in time and radial velocity with exceptionally deep Ca II absorption; (ii) a long-lived blue-shifted Ca II feature near -5 km/s that appears continuously for about a year (2017-2018), which the authors note may be difficult to explain with the classical exocomet picture; and (iii) two blue-shifted, strongly accelerating Ca II features observed on consecutive nights in December 2019 that depart from Keplerian motion, which the authors hypothesize to be evidence for the destruction of the comet nuclei shortly after periastron passage. The paper also documents significant variability in the Na I and Ca II disk lines that complicates the reference-spectrum approach.","tokens_in":19179,"tokens_out":8606,"duration_ms":92465,"significance":"If the results hold, this is a valuable contribution to exocomet studies. The Na I rarity result constrains photoionization and sublimation models, and the simultaneous Na I + Ca II detections provide a physical link between extreme outgassing and detectable neutral sodium that goes beyond the classical Falling Evaporating Bodies framework. The reported disk-line variability is an important caution for the common percentile-based reference-spectrum methodology, and the long-lived 2017-2018 feature challenges existing models. The fragmentation scenario for the 2019 events, with its analogy to Kreutz-family sungrazers such as comet Lovejoy, is an intriguing and falsifiable hypothesis. The paper's strengths include its unusually large and uniform dataset (9071 spectra over 17 years), the systematic telluric correction, transparent Bayesian fits with explicit priors (Tables A.1, A.2), and the authors' own candid statements about the limitations of the reference-spectrum assumption.","major_comments":[{"comment":"The central novel claim - that the 2019 features depart from Keplerian motion because the comet nuclei fragmented - rests entirely on the fitted centroid v0 of single-Gaussian fits to a blended, time-varying absorption region. The paper itself introduces a second static component near -39 km/s on night 2 (Fig. 9), omits unstable fits at the end of that night, and reports a simultaneous dramatic broadening of the line (Fig. 10). A second unresolved or gradually emerging absorption component could shift and broaden the centroid of a single-Gaussian fit without any dynamical acceleration of the absorbing body. The paper's argument against the tail-egress scenario (Section 3.2) does not test this unresolved-blend degeneracy, particularly for night 1 where only one component was fitted. I request a direct test: refit night 1 with a two-component model (or a time-varying blend) and compare the model evidence, and/or show that single-component residuals remain clean throughout the sequence. Without such a test, the attribution of the non-Keplerian departures to nucleus destruction (abstract and Section 4) is not uniquely supported.","section":"Section 3.2, Figs. 8-10"},{"comment":"The reference spectrum is the baseline against which all time-variable absorption is measured, and the paper itself documents that the disk lines vary substantially: the Na I D2 line depth varies by more than a factor of two (from ~5% to ~13%, Fig. A.4), and the Ca II disk line can be narrower and shifted on a single night (Fig. A.5). The long-lived -5 km/s feature of 2017-2018 (Fig. A.6) sits close to the disk line, exactly where the reference spectrum is least certain, and the authors themselves ask whether it might be 'due to some other component of the circumstellar material' (Section 3.3). The impact of the documented disk-line variability on the parameters of the long-lived feature and on the weaker Na I fits is not quantified. I ask that the reference spectrum be perturbed within its stated 0.5-2% uncertainty (or reconstructed with different percentile thresholds) and that the affected claims be re-derived; this is particularly important because the 2019 features are well separated from the disk line, so the reference-spectrum concern is most acute precisely for the long-lived feature claim, not for the fragmentation claim.","section":"Sections 2.3, 3.3, Figs. A.4-A.6"},{"comment":"The claim of a 'general absence' of exocometary Na I absorption, with only two clear (~2%) events out of 198 nights, is established by visual inspection of night-averaged spectra with an ad hoc '>~2%' threshold; no per-night detection limit is given. Given that telluric Na I is present at the 1-5% level in 71 of 198 nights (Fig. A.3), that the Na I disk lines vary in depth by more than a factor of two (Fig. A.4), and that the weakest fitted exocomet feature has a constrained line depth of only 0.89 +/- 0.13% (Fig. 7), the sensitivity of the search should be quantified before concluding that Na I is generally absent. A map of the minimum detectable line depth per night (or per epoch) would make the central sodium rarity claim falsifiable and would also justify the statement that the Earth's atmospheric sodium is not being mistaken for exocometary absorption in the remaining nights.","section":"Section 3.1, Figs. A.2-A.3"}],"minor_comments":[{"comment":"Section 4 states that the two anomalous events were 'observed in 2018', whereas Section 3.2 and the abstract date them to 11 and 12 December 2019; the year should be corrected to 2019.","section":"Section 4"},{"comment":"The caption of Fig. 6 dates the strongest sodium event to 'March 17 2008', while Section 3.1 and Fig. 5 give 'March 16 2008'; the dates should be reconciled.","section":"Fig. 6 vs. Section 3.1"},{"comment":"The abstract contains a typo, 'blue-shifted exocomes', which should read 'exocomets'.","section":"Abstract"},{"comment":"The final paragraph of the conclusions contains a grammatically broken sentence beginning 'We systematically corrected these spectra using Molecfit and that exocometary sodium absorption is relatively rare'; this should be rephrased.","section":"Section 4, last paragraph"},{"comment":"The companion paper is referred to inconsistently as both 'Jaworska et al. in prep.' and 'Jarworska et al. in prep.' within the same section; the spelling should be unified.","section":"Section 3.2"},{"comment":"The vertical axis of Fig. 10 is labeled 'Centroid velocity v0 (km/s)' with values ranging from about 24 to 34 km/s, while the text quotes the features at approximately -32 km/s; please clarify whether the plotted quantity is |v0| or else correct the sign convention.","section":"Fig. 10"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies on a companion paper (Jaworska et al., in prep.) for the dynamical interpretation of blue-shifted exocomets and the icy-progenitor scenario; since that paper is not available for review, the testable content here is limited to the spectroscopic analysis itself, which is a substantial and well-documented dataset. The fragmentation claim is the most novel result and is also the point most in need of the additional test requested in major comment 1. I would also gently encourage the authors to include a statement of code and data availability, since the JAX/NumPyro workflow and the public HARPS archive would otherwise make the analysis straightforwardly reproducible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read this one. First, it's a genuinely useful observational paper: the first systematic, telluric-corrected Na I exocomet survey over 198 nights of HARPS data, and it finds that sodium absorption is rare and occurs only during the deepest Ca II events. Second, the eye-catching claim—two 2019 comets showing non-Keplerian acceleration, interpreted as post-periastron nucleus destruction—is plausible but far from proven; it rests on single-Gaussian centroid fits that an unresolved second component could mimic.\n\nThe paper does a lot right. It grinds through 9071 spectra, corrects tellurics with molecfit, builds a reference spectrum with documented truncation and Tikhonov smoothing, and fits lines with explicit Bayesian priors (Tables A.1, A.2). The Na I detections at ~1–2% depth are supported by fits with decent constraints. The authors also honestly document disk line variability (Figs A.4, A.5), which is relevant because that variability is the main threat to their reference-spectrum assumption.\n\nWhere it gets soft: the 'general absence' of Na I is asserted from visual inspection of night-averaged spectra, with no formal detection-limit map. That makes the rarity claim a statement about what's visible, not a quantitative upper limit. The bigger issue is the fragmentation story. The stress-test gets it right: the centroid of a single Gaussian can shift and broaden if a second, weaker component emerges or changes, without any dynamical acceleration. The authors handled the known blend on night 2 by adding a component, but night 1 is still vulnerable, and the line-width jump in Fig 10 is precisely the signature of a blend. Their argument against tail egress doesn't rule out an unresolved blend. Minor point: the conclusions say '2018' while the abstract and results say '2019'—sloppy.\n\nWho's this for? Exocomet and debris-disk researchers will want to see the Na I survey and the 2019 events. It deserves a serious referee. I'd send it out, with a request for a detection-limit analysis and a more careful treatment of unresolved blends before the destruction claim is accepted. That's a revision, not a reject.","headline":"A careful HARPS study that establishes rare, Ca II-linked Na I exocomet absorption, but the 2019 non-Keplerian events are not yet proof of nucleus destruction.","tokens_in":19758,"tokens_out":3386,"would_cite":true,"duration_ms":36630,"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":"Sodium absorption from beta Pictoris exocomets is almost never observed, and two 2019 Ca II events show sudden non-Keplerian acceleration that the paper interprets as comet nuclei breaking apart after periastron.","keywords":["exocomets","beta Pictoris","HARPS spectroscopy","calcium H&K lines","sodium D lines","circumstellar disk variability","telluric correction","non-Keplerian acceleration"],"falsifier":"A re-analysis of the 2019 nights that models the Ca II disk line as time-variable, allowing its depth, width, and velocity to drift on hour timescales, would settle the destruction claim: if the apparent acceleration reversal and line broadening disappear once the disk line is free to vary, they are disk variability rather than comet fragmentation. Separately, detecting sodium absorption in an exocomet event with only moderate Ca II depth would directly refute the claim that Na I becomes detectable only during the most extreme Ca II outgassing.","tokens_in":18727,"feed_emoji":"☄️","tokens_out":7695,"duration_ms":76812,"temperature":0.7,"pith_summary":"The paper tries to establish two things from nearly two decades of HARPS spectra of beta Pictoris. First, exocometary sodium absorption is rare: in 198 nights only two clear ~2% events plus one weaker ~1% feature lasting 13 nights appear, and they coincide with the deepest Ca II absorptions at the same redshift, implying sodium becomes visible only when the same bodies evaporate calcium extremely strongly, because most sodium is quickly photo-ionized close to the star. Second, two blue-shifted accelerating Ca II features observed on 11 and 12 December 2019 deviate sharply from Keplerian motion within hours, which the authors interpret as destruction of the comet nuclei shortly after periastron passage, analogous to sungrazing comets like Lovejoy. The paper also argues that the Ca II and Na I disk lines vary on observable timescales, so static reference spectra can misattribute disk variability to exocomets. If right, these results sharpen the use of sodium as a tracer of extreme outgassing and add a concrete, observable fate for stargrazing exocomets.","feed_headline":"Two beta Pictoris exocomets appear to shatter near the star","feed_subtitle":"17 years of HARPS data also show sodium only in the deepest calcium events.","key_machinery":"The analysis rests on three working parts. An unbiased reference spectrum is built by taking the 80th percentile of flux values per wavelength channel and correcting the truncation bias with a truncated-normal formula, so that time-variable exocomet absorption is not baked into the baseline; telluric water lines near the sodium doublet are removed before any sodium analysis. Exocomet lines are modeled as homogeneous Gaussian clouds with central optical depth $\\tau_0$, a fill factor $f$ for stellar-disk coverage, and doublet oscillator-strength ratio fixed at 2, fitted with an auto-differentiable Bayesian sampler that yields posterior velocities, widths, depths, and accelerations. For the 2019 features, the paper computes families of Keplerian orbits consistent with the early linear acceleration and then compares the later velocity curves against those predictions, so the departure from Keplerian motion is the observable that carries the destruction claim.","core_discovery":"After systematically telluric-correcting 9071 HARPS spectra and calibrating them against an unbiased reference spectrum, the paper finds that exocometary Na I is generally absent: only two nights (March 2008 and March 2009) show clear ~2% absorption in both sodium D lines, plus a ~1% feature persisting ~13 nights in February 2004. In all cases the sodium appears red-shifted at the same radial velocity as exceptionally deep Ca II H&K absorption, so the paper concludes that the Na I is produced by the same evaporating bodies, and that detectable sodium requires the most extreme Ca II outgassing because sodium is otherwise photo-ionized. For the 2019 events, the paper shows two blue-shifted Ca II features on consecutive nights that accelerate roughly linearly at first, then suddenly depart from Keplerian acceleration; one reverses its acceleration within about an hour, and line widths grow simultaneously. The paper hypothesizes this is the final fragmentation of the comet nuclei and blow-out of the Ca II tail shortly after periastron, rather than an egress artifact. It further reports a ~5 km/s blue-shifted Ca II feature persisting from early 2017 to mid-2018 that the classical exocomet model may not explain, and documents Na I disk-line depth variations by more than a factor of two as well as a night where the Ca II disk line narrows and shifts, concluding that the circumstellar disk lines themselves are observably variable.","pith_inferences":["If sodium traces only extreme outgassing, comparing Ca II and Na I line depths across a large sample could provide a quantitative threshold for photo-ionization shielding, effectively mapping where in the coma neutral sodium can survive.","The destruction scenario predicts that some stargrazing exocomets should show a characteristic sequence within a single night: linear acceleration, then a sudden acceleration reversal or jump with line broadening; this sequence could be searched for automatically in other archival time-series of beta Pictoris and of other debris-disk stars.","The year-long 2017-2018 feature, if it is disk gas, would connect exocomet activity to the slowly varying circumstellar disk and might change how the Falling Evaporating Bodies population is separated from the disk in future analyses.","A testable extension: re-fit the 2019 nights with a model that lets the disk line parameters vary with time; if the non-Keplerian acceleration survives, the destruction interpretation is strengthened, and if it does not, the event is disk variability."],"forward_implications":["Sodium D-line absorption can be used as a specific diagnostic of the most extreme Ca II-producing exocomets; most exocomet events will be invisible to sodium surveys.","Multi-species observations of rare deep events can probe the evaporation and ionization structure of exocomet clouds, though abundance ratios will require physical cloud models rather than homogeneous ones.","The 2019 events imply that stargrazing comets can be destroyed within hours near or after periastron, and that sudden non-Keplerian acceleration with line broadening is a signature of that destruction.","The long-lived 2017-2018 Ca II feature and the variable disk lines mean that some 'exocomet' absorption may actually be circumstellar disk gas, so future surveys need time-variable disk models.","If the 2019 pair are on the same orbit, their different behavior after two hours suggests fragmentation is stochastic, not a predictable orbital effect."],"supporting_citations":[{"why":"Supplies the previous statistical analysis of HARPS exocomet Ca II populations that this study extends, and the earlier percentile-based reference-spectrum approach.","marker":"Kiefer et al. (2014)"},{"why":"Provides the truncated-normal formula used to unbias the reference spectrum from the 80th-percentile flux samples.","marker":"Tallis (1961)"},{"why":"Provides the telluric-correction code used to remove water lines near the sodium doublet.","marker":"Smette et al. (2015)"},{"why":"Provides the telluric-correction implementation applied to the HARPS spectra.","marker":"Kausch et al. (2015)"},{"why":"Establishes the framework for constraining exocomet orbits from observed acceleration, which the paper uses to define the Keplerian baseline.","marker":"Kennedy (2018)"},{"why":"Documents the post-periastron destruction of sungrazing comet Lovejoy, the solar-system analogue for the proposed nucleus fragmentation.","marker":"Sekanina & Chodas (2012)"},{"why":"Provides the mean-motion-resonance formation model and the prediction that blue-shifted exocomets are rare, framing the dynamical interpretation of the 2019 events.","marker":"Beust et al. (2024)"},{"why":"Supplies the modern multi-species exocomet line-fitting model with doublet optical depths and fill factors that the Na I and Ca II fits build on.","marker":"Vrignaud et al. (2024)"},{"why":"Provides the Fe I disk-line measurement of the systemic velocity, which sets the zero-point for the radial velocities in this paper.","marker":"Kiefer et al. (2019)"},{"why":"First reported accelerating Ca II absorption events in beta Pictoris, the phenomenon the 2019 features are compared with.","marker":"Ferlet et al. (1987)"}],"fun_headline_variants":["Two beta Pic exocomets appear to shatter near star: HARPS","Sodium in beta Pic exocomets only during deepest calcium dips","17 yrs of HARPS: exocomets break up, sodium mostly absent","Beta Pictoris exocomets: sodium rare, only when Ca II is extreme","Long-lived Ca II absorption challenges exocomet model in beta Pic"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the 80th-percentile reference spectrum is a faithful static baseline for the stellar and circumstellar disk lines, so that all residual time-variable absorption can be attributed to exocomets; since the paper itself shows the Na I disk lines vary in depth by more than a factor of two and the Ca II disk line can narrow and shift on a single night, disk variability could contaminate the long-lived 2017-2018 feature and the fitted parameters of the 2019 events.","fun_headline_variants_meta":{"raw":{"variants":["Two beta Pic exocomets appear to shatter near star: HARPS","Sodium in beta Pic exocomets only during deepest calcium dips","17 yrs of HARPS: exocomets break up, sodium mostly absent","Beta Pictoris exocomets: sodium rare, only when Ca II is extreme","Long-lived Ca II absorption challenges exocomet model in beta Pic"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000388,"raw_usage":{"total_tokens":2194,"prompt_tokens":1241,"completion_tokens":953,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":857,"completion_tokens_details":{"reasoning_tokens":855}},"tokens_in":857,"tokens_out":953,"duration_ms":9899,"temperature":1.0,"reasoning_tokens":855,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:26:18.288382+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A re-analysis of the 2019 nights that models the Ca II disk line as time-variable, allowing its depth, width, and velocity to drift on hour timescales, would settle the destruction claim: if the apparent acceleration reversal and line broadening disappear once the disk line is free to vary, they are disk variability rather than comet fragmentation. Separately, detecting sodium absorption in an exocomet event with only moderate Ca II depth would directly refute the claim that Na I becomes detectable only during the most extreme Ca II outgassing.","supporting_citations":[{"cited_title":"2014, , 514, 462","cited_arxiv_id":null,"evidence_quote":"Supplies the previous statistical analysis of HARPS exocomet Ca II populations that this study extends, and the earlier percentile-based reference-spectrum approach."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the truncated-normal formula used to unbias the reference spectrum from the 80th-percentile flux samples."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the framework for constraining exocomet orbits from observed acceleration, which the paper uses to define the Keplerian baseline."},{"cited_title":"& Chodas , P","cited_arxiv_id":null,"evidence_quote":"Documents the post-periastron destruction of sungrazing comet Lovejoy, the solar-system analogue for the proposed nucleus fragmentation."},{"cited_title":"2024, , 683, A89","cited_arxiv_id":null,"evidence_quote":"Provides the mean-motion-resonance formation model and the prediction that blue-shifted exocomets are rare, framing the dynamical interpretation of the 2019 events."},{"cited_title":"2024, , 684, A210","cited_arxiv_id":null,"evidence_quote":"Supplies the modern multi-species exocomet line-fitting model with doublet optical depths and fill factors that the Na I and Ca II fits build on."},{"cited_title":"2019, , 621, A58","cited_arxiv_id":null,"evidence_quote":"Provides the Fe I disk-line measurement of the systemic velocity, which sets the zero-point for the radial velocities in this paper."},{"cited_title":"M., & Vidal-Madjar , A","cited_arxiv_id":null,"evidence_quote":"First reported accelerating Ca II absorption events in beta Pictoris, the phenomenon the 2019 features are compared with."}],"review_version":1}