{"id":"94b10801-58b0-4409-acd8-4c0459270cad","arxiv_id":"2412.01878","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Lithium enhancements in several polluted white dwarfs are most consistent with Big Bang and Galactic nucleosynthesis rather than continental crust or spalled-ring exomoons, with one object defying all three explanations.","lead":"This paper reports new spectra and uniform atmospheric models for six polluted white dwarfs to test three explanations for their unusually high lithium levels. The data support a Big Bang and Galactic nucleosynthesis origin for three of them, rule out spalled-ring exomoons, and leave one object unexplained.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Steady-state abundances and the Section 5 hypothesis rankings depend on diffusion timescales extrapolated below the Koester et al. (2020) grid; the ad hoc 0.2 dex uncertainty may underestimate the factor-of-3 corrections implied by Heinonen et al. (2020).","rationale":"The paper is a careful, well-hedged evaluation, and the exomoon spallation timescale argument is independently compelling. However, the 'most plausible' ranking for the three positively classified white dwarfs depends on steady-state abundance ratios whose diffusion corrections are roughly 0.5 dex and whose uncertainty of 0.2 dex is assumed rather than derived from a comparison with a valid grid. The authors explicitly flag in Section 6.3 that updated diffusion coefficients can change relative timescales by a factor of three, which would shift the inferred log(Li/Ca)_SSP by about 0.5 dex and could alter the compatibility assignments for WD J1644 and WD J2317. This does not overturn the paper's central claim, but it does justify the CONDITIONAL verdict and makes a specific recomputation necessary before the hypothesis ranking can be considered secure.","tokens_in":33480,"tokens_out":7148,"duration_ms":74979,"concrete_test":"Compute diffusion timescales for Li, Na, Mg, K, Cr, and Fe at each object's Teff and log g using the Heinonen et al. (2020) method, or full envelope diffusion calculations, instead of the linear extrapolation of Koester et al. (2020); compare log(tau_el/tau_Ca) to Table 5. If any value differs by more than 0.2 dex, recompute the Table 4 SSP abundances and rerun the Section 5 compatibility analysis, checking whether any object moves across the BBN, CI-chondrite, or continental-crust matching boundaries.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central conclusion that Big Bang and Galactic nucleosynthesis best explains WD J1644, WD J1824, and WD J2317 is reached using steady-state abundances from Equation 2, whose corrections are the relative diffusion timescales of Section 3.2.4. These timescales are not computed at the stars' parameters: every Teff (3540-5020 K) lies below the Koester et al. (2020) grid, so the authors linearly extrapolate log(tau_el/tau_Ca) from Teff < 10,000 K and add a 0.2 dex uncertainty by a single random draw. The correction itself is about 0.5 dex for Li/Ca (Table 5), so a factor-of-3 error in relative timescales changes inferred log(Li/Ca)_SSP by about 0.5 dex. The paper itself notes in Section 6.3 that Heinonen et al. (2020) found factor >= 3 differences in relative diffusion timescales using updated coefficients, yet the adopted grid still uses Paquette et al. (1986) coefficients. If, for example, true tau_Li/tau_Ca is smaller than extrapolated, WD J1644's SSP Li/Ca would shift from the BBN/CI-matching region toward CI only, and WD J2317's match to the [Fe/H] = -3 bin could move by roughly 1.5 sigma. Since Section 5's compatibility matrix and Table 8 hinge on these shifts, the ranking of hypotheses is not secure until the diffusion coefficients are recomputed at the actual Teff and log g values.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents new SOAR/Goodman and Gemini/GMOS-N spectroscopy of three previously known lithium-polluted white dwarfs (WD J1824+1213, WD J2317+1830, LHS 2534) and one non-lithium DZ (SDSS J1636+1619), and reanalyzes two additional Li-polluted systems from the literature using a homogeneous set of low-Te_eff white dwarf atmospheric models. It infers accreted-planetesimal abundances under increasing-phase, steady-state, and decreasing-phase accretion assumptions, uses SAGA main-sequence abundances together with a BBN lithium assumption to construct Galactic nucleosynthetic reference ratios, evaluates the three published explanations for white-dwarf Li enhancement, and concludes that Big Bang plus Galactic nucleosynthesis is most plausible for WD J1644-0449, WD J1824+1213, and WD J2317+1830; that SDSS J1330+6435 is too weakly constrained; that LHS 2534, as presently modeled, is inconsistent with all three hypotheses; and that the spalled exomoon scenario is highly unlikely in all cases.","tokens_in":33785,"tokens_out":8791,"duration_ms":97472,"significance":"The paper is a useful contribution because it applies one consistent modeling framework to all known Li-polluted cool DZs and confronts three competing hypotheses with a common abundance dataset. The new observations, the homogeneous model fits, the Monte Carlo treatment of abundance and diffusion uncertainties, and the use of independent age and kinematic diagnostics are genuine strengths. If the diffusion-timescale systematics can be made robust, the conclusion that BBN/Galactic nucleosynthesis is the most plausible explanation for at least three of the systems would be a meaningful step beyond the earlier discovery papers. The paper also provides a concrete quantitative argument against the spallation-exomoon hypothesis and appropriately hedges the LHS 2534 result via the phrase 'as presently modeled.'","major_comments":[{"comment":"The steady-state abundances that drive the hypothesis rankings are corrected using relative diffusion timescales that are linearly extrapolated below the Koester et al. (2020) grid for every object in Table 2. The added 0.2 dex uncertainty is assigned by a single random draw, and Section 6.3 itself reports that the updated Heinonen et al. (2020) coefficients change relative Ca/Si diffusion timescales by a factor of at least three, i.e. by about 0.5 dex. Since the SSP corrections in Table 5 are of order 0.5 dex for Li/Ca (log tau_Li/tau_Ca between 0.51 and 0.65 dex), a factor-of-three change in tau_Li/tau_Ca would shift the inferred log(Li/Ca)_SSP values in Table 4 by roughly 0.5 dex and would directly alter the compatibility statements in Section 5 and Table 8. This is the load-bearing step for the central claim that BBN/Galactic nucleosynthesis is most plausible for WD J1644-0449, WD J1824+1213, and WD J2317+1830, so the conclusion is not secure until the diffusion coefficients are recomputed at the actual Teff and log g values with updated physics, or until a conservatively propagated systematic error is incorporated into the Section 5 verdicts.","section":"§3.2.4, Eq. (2), Table 5, and §5/Table 8"},{"comment":"The log(Li/Ca) predictor bins used as Galactic nucleosynthesis references are constructed by taking a fixed BBN lithium abundance and each SAGA star's measured Ca abundance, thereby intentionally neglecting post-BBN lithium production. The paper acknowledges that this underestimates the true protostellar Li/Ca by up to about 0.5 dex at solar metallicity, but this systematic offset is not included in the bin dispersions quoted in Table 4 and is not propagated into the sigma-based compatibility statements in Section 5. The effect is likely small for the [Fe/H] = -2 and -3 bins that are central to WD J2317+1830, but it is not negligible for the [Fe/H] = 0 and -1 comparisons used for WD J1644-0449 and WD J1824+1213; the assumed Li evolution should be stated explicitly and the systematic should be included in the comparison.","section":"§3.5, Eq. (4), Table 4"},{"comment":"For WD J1644-0449, the continental-crust scenario is retained as compatible only if the accretion is in the decreasing phase (Table 8), while the BBN/Galactic nucleosynthesis scenario is assessed using the steady-state values from Equation 2. The decreasing-phase analysis is explicitly labeled in Section 3.2.3 as 'a general indicator rather than a firm value,' and no uncertainty is attached to the time since accretion ceased. Because the accretion phase is not independently known for WD J1644-0449, the relative ranking of the two remaining hypotheses for this object is partly a statement about the assumed phase; the paper should either add a quantitative treatment of the phase uncertainty or soften the conclusion for this object accordingly.","section":"§3.2.3 and §5.1.2"}],"minor_comments":[{"comment":"The sentence 'The new spectra spectra we present' contains a duplicated word, and the text 'WD J1824+1214' should read 'WD J1824+1213.'","section":"§2.1"},{"comment":"The conclusion text 'WD J644–0449' is a typo for 'WD J1644–0449.'","section":"§7"},{"comment":"The caption contains 'are shown are shown' and should be corrected to 'are shown.'","section":"Figure 6 caption"},{"comment":"The spallation timescale argument assumes that the hypothetical exomoon-forming rings are located near Europa's orbital distance where the Jovian radiation environment is modeled; a brief sentence justifying this geometric choice would strengthen the timescale argument against the exomoon hypothesis.","section":"§4"},{"comment":"Equation (5) introduces MineLi and MineCa without explicitly restating their definitions immediately after the equation; the surrounding text defines them only indirectly, so a one-sentence definition would improve readability.","section":"Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the authors have been appropriately careful with several conclusions, especially LHS 2534. My main concern is the diffusion-timescale extrapolation below the Koester et al. (2020) grid; because the central ranking in Table 8 depends on SSP abundances, this needs to be addressed with updated coefficients or a defensible systematic error treatment before acceptance. There are no concerns about novelty or citation practices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a genuinely useful paper. It brings new SOAR/Gemini spectra for three Li-polluted white dwarfs plus a metal-polluted control, re-derives all abundances with one atmospheric model grid, and evaluates the three proposed explanations—Big Bang plus Galactic nucleosynthesis, continental crust, and spalled-ring exomoons—on equal footing. The strongest new result is the quantitative case against the exomoon hypothesis: the threshold-energy and proton-flux argument, plus the progenitor-lifetime mismatch for WD J2317+1830, makes that channel look very unlikely. The K/Ca and Cr/Ca constraints that rule out unaltered continental crust for WD J1644 and LHS 2534 are also solid. The paper is refreshingly honest about its own limits: it flags the diffusion-timescale extrapolation, the questionable LHS 2534 magnetic modeling, and the exclusion of two recently published Li-polluted white dwarfs whose authors favored continental crust.\n\nThe load-bearing soft spot is Section 3.2.4. Every Teff is below the Koester et al. (2020) grid, so the relative diffusion timescales are linearly extrapolated and assigned a single ad hoc 0.2 dex uncertainty. The paper itself cites Heinonen et al. (2020) as finding factor >=3 differences with the older Paquette coefficients, then keeps using Paquette. For Li/Ca the correction is about 0.5 dex, so a factor-of-3 change in tau_Li/tau_Ca shifts the inferred SSP Li/Ca by roughly that amount—enough to move WD J1644 and WD J2317 between compatibility bins. The stress-test estimate of a ~1.5 sigma shift for WD J2317 is plausible. The 0.2 dex uncertainty thus understates the systematic risk. The qualitative conclusion that spallation is unlikely survives this; the ranking between BBN and CI-chondrite-like compositions for individual objects is less secure. The paper would be substantially stronger with a sensitivity analysis using Heinonen et al. coefficients or a direct calculation at the stars' Teff and log g.\n\nLHS 2534 defies all three hypotheses as modeled; the authors attribute this to magnetism and spots, which is plausible but means that object contributes little until better models exist. Excluding the two recent continental-crust-favoring objects is defensible on homogeneity grounds, but it leaves a gap in the sample.\n\nBottom line: solid, careful work, worth serious refereeing. A referee should ask for the diffusion-timescale sensitivity analysis and a more principled uncertainty estimate. The central argument holds up in its strongest form—exomoons are not the answer—while the BBN-versus-crust ranking needs that extra work before it is secure.","headline":"A careful, honest evaluation of the three lithium-excess hypotheses, with a solid case against spallation but a diffusion-timescale extrapolation that leaves the BBN-vs-crust ranking less secure than the paper claims.","tokens_in":34358,"tokens_out":2046,"would_cite":true,"duration_ms":22742,"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 paper argues that the lithium excess in three polluted white dwarfs comes from Big Bang and Galactic nucleosynthesis rather than continental crust or spallation-formed exomoons.","keywords":["white dwarfs","polluted white dwarfs","lithium enhancement","exoplanetesimal accretion","Big Bang nucleosynthesis","Galactic nucleosynthesis","continental crust","spallation"],"falsifier":"Compute the actual diffusion timescales for these five white dwarfs at their effective temperatures using updated diffusion coefficients; if any relative timescale (for example $\\log(\\tau_{\\rm Li}/\\tau_{\\rm Ca})$) differs from the extrapolated grid value by more than about 0.2 dex, the inferred planetesimal ratios change enough to move stars between the Big Bang bins and the continental crust match. A simpler check: obtain a K abundance or tighter K limit for LHS 2534, since its K/Ca is the ratio that currently rules out every hypothesis.","tokens_in":33274,"feed_emoji":"🔭","tokens_out":7069,"duration_ms":71402,"temperature":0.7,"pith_summary":"The paper asks where the extra lithium seen in some polluted white dwarfs comes from. It tests the three live hypotheses—primordial (Big Bang plus Galactic nucleosynthesis), geological (accreted continental crust), and an exomoon built from proton-spalled ring material—against fresh spectra and a homogeneous set of atmospheric models for five lithium-polluted white dwarfs. The authors conclude that the first hypothesis is the most plausible for WD J1644−0449, WD J1824+1213, and WD J2317+1830, that the spalled-exomoon route is highly unlikely for every case, and that LHS 2534, as presently modeled, matches none of the three. If right, the result turns lithium abundance in cool white dwarfs into a probe of how early a planetary system formed, while demoting the more exotic formation stories.","feed_headline":"Big Bang lithium explains three white dwarfs' excess","feed_subtitle":"Three candidate sources are tested against five lithium-polluted white dwarfs; one fits none.","key_machinery":"The load-bearing machinery is the abundance-recovery chain for polluted white dwarfs. Photospheric metal abundances come from model atmospheres built for cool, dense, helium-rich white dwarfs; those are then converted to the accreted body's abundances using phase-dependent corrections: in the increasing phase no correction is needed, in the steady state $\\log(\\mathrm{el_1/el_2})_{\\rm SSP} = \\log(\\mathrm{el_1/el_2})_{\\rm phot} + \\log(\\tau_{\\rm el_2}/\\tau_{\\rm el_1})$, and in the decreasing phase an extra time-decay term is added. The relative diffusion timescales $\\tau$ are linearly extrapolated from published grids down to these stars' effective temperatures, with a 0.2 dex uncertainty added in quadrature to cover the extrapolation. The corrected ratios are then matched against three reference sets—CI chondrites, continental crust, and binned main-sequence star abundances that represent Galactic nucleosynthetic evolution—with system ages and Galactic kinematics used to decide which matches are physically credible.","core_discovery":"The central claim is that the lithium excesses are inherited from the natal gas of old, metal-poor planetary systems rather than from any special planetary or ring process. Using a homogeneous set of cool white dwarf atmosphere models, the paper converts measured photospheric abundances into parent planetesimal abundances via diffusion-phase corrections. Those corrected ratios for WD J1644−0449, WD J1824+1213, and WD J2317+1830 land on the curves that Big Bang plus Galactic nucleosynthesis predicts for metallicities between roughly $[\\mathrm{Fe/H}] = -1$ and $-3$. The continental crust scenario fails where K/Ca can be measured—the crust is too potassium-rich—and for LHS 2534 the Fe/Ca and Cr/Ca are far too high for crust; the spallation scenario fails because Jupiter-like radiation belts do not supply enough protons above the threshold energy within the progenitor's lifetime. The paper therefore ranks Big Bang and Galactic nucleosynthesis first for three stars, leaves SDSS J1330+6435 undecided, and reports that LHS 2534, as modeled, matches none of the three.","pith_inferences":["A testable prediction of the paper's ranking is that future samples of cool polluted white dwarfs will show a high-Li/Ca tail concentrated in kinematically old, metal-poor populations rather than a smooth spread caused by crustal composition.","The same proton-threshold argument that kills spallation for lithium also bears on the beryllium-rich white dwarfs that motivated the exomoon model; recomputing those with the lower high-energy proton flux would test whether the original Be interpretation survives.","If magnetic surface spots are skewing LHS 2534's abundance fit, then the other magnetic Li-polluted white dwarfs found subsequently may also be misclassified; time-resolved spectroscopy of their line strengths could test whether spotted surfaces are the fourth explanation.","A sharper test for WD J2317+1830 would be a high signal-to-noise K measurement: detection near the current upper limit would confirm the very metal-poor bin, while a low K/Ca would indicate thermal depletion on top of primordial lithium."],"forward_implications":["A lithium excess in a polluted white dwarf becomes usable as a rough age and metallicity indicator for its planetary system, since the inferred Li/Ca should track how little Galactic enrichment occurred before the planets formed.","The spalled-exomoon explanation can be deprioritized for lithium; future lithium detections should first be compared with the Big Bang/Galactic nucleosynthesis predictions rather than with ring-spallation models.","The continental crust explanation survives only in heavily altered form: for most of these stars the potassium ratio rules out unprocessed crust, and where it does not, large systematic uncertainties still prevent a clean identification.","A null-detection control (SDSS J1636+1619) shows that even without lithium, upper limits on Na, K, Mg, and Fe can constrain whether an accreted body was primitive or thermally altered."],"supporting_citations":[{"why":"Proposed the Big Bang and Galactic nucleosynthesis hypothesis and supplied the first lithium detections and fits for WD J1644−0449 and SDSS J1330+6435.","marker":"Kaiser et al. (2021)"},{"why":"Provided the high-resolution spectra for WD J1824+1213, WD J2317+1830, and LHS 2534 and proposed the continental crust hypothesis.","marker":"Hollands et al. (2021)"},{"why":"Introduced the spalled icy exomoon hypothesis whose beam flux and timescales the paper re-examines.","marker":"Doyle et al. (2021)"},{"why":"Supplied the diffusion timescale grids that are extrapolated to low effective temperatures for the steady-state abundance corrections.","marker":"Koester et al. (2020)"},{"why":"Derived the increasing-, steady-state-, and decreasing-phase equations used to convert photospheric abundances to parent body abundances.","marker":"Harrison et al. (2018)"},{"why":"Contributed the aggregated main-sequence stellar abundances from which the Galactic nucleosynthetic evolution curves are built.","marker":"Suda et al. (2017)"},{"why":"Gives the CI Chondrite and solar system reference abundances against which all lithium enhancements are measured.","marker":"Lodders (2019)"},{"why":"Provides the continental crust reference composition used for the geologic differentiation comparison.","marker":"Rumble et al. (2019)"}],"fun_headline_variants":["Big Bang lithium wins over crust and exomoons","Excess lithium in white dwarfs points to Big Bang","White dwarf lithium: cosmic origin beats crust and moons","Three white dwarfs' lithium points to Big Bang, not moons","Lithium-polluted white dwarfs favor Big Bang origins"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The inferred planetesimal abundances rest on relative diffusion timescales that are linearly extrapolated below the published grid's temperature range; if the true timescales differ by more than the ad hoc 0.2 dex uncertainty, the abundance ratios that decide which hypothesis wins could shift.","fun_headline_variants_meta":{"raw":{"variants":["Big Bang lithium wins over crust and exomoons","Excess lithium in white dwarfs points to Big Bang","White dwarf lithium: cosmic origin beats crust and moons","Three white dwarfs' lithium points to Big Bang, not moons","Lithium-polluted white dwarfs favor Big Bang origins"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001053,"raw_usage":{"total_tokens":4487,"prompt_tokens":1076,"completion_tokens":3411,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":692,"completion_tokens_details":{"reasoning_tokens":3330}},"tokens_in":692,"tokens_out":3411,"duration_ms":26059,"temperature":1.0,"reasoning_tokens":3330,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:51:54.565576+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the actual diffusion timescales for these five white dwarfs at their effective temperatures using updated diffusion coefficients; if any relative timescale (for example $\\log(\\tau_{\\rm Li}/\\tau_{\\rm Ca})$) differs from the extrapolated grid value by more than about 0.2 dex, the inferred planetesimal ratios change enough to move stars between the Big Bang bins and the continental crust match. A simpler check: obtain a K abundance or tighter K limit for LHS 2534, since its K/Ca is the ratio that currently rules out every hypothesis.","supporting_citations":[],"review_version":1}