{"id":"84e9b6f8-e018-437f-b31f-5a00a5904dbd","arxiv_id":"2607.24531","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"Peculiar long GRBs arise from magnetar–massive-WD mergers whose accretion-driven giant-flare forest powers the main burst and ejects crustal r-process material for the kilonova.","lead":"The paper argues that peculiar long gamma-ray bursts like GRB 211211A and 230307A come from neutron-star–white-dwarf mergers that repeatedly fire magnetar giant flares. Those flares both make the main burst and eject crust that builds r-process elements for the kilonova, tying prompt emission, afterglow, and heavy-element production into one engine.","discovery_kind":"unification","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The kilonova fit demands Mej ~10^-2 M_sun, but the proposed mechanism delivers 10^-5–10^-3 M_sun; the bridge — a relativistic beaming factor applied to non-relativistic ejecta plus ~10-1000x more unseen flares — is internally inconsistent with the spike census the paper itself derives.","rationale":"The reader's flagged concern (cold, non-accreting crust under sustained hyperaccretion and repeated flares) is real and the authors themselves concede it in Section 4, so the reader's weakest_assumption is a fair choice. However, I find a sharper, internally checkable problem one step downstream: even granting the Cehula et al. per-flare ejecta numbers wholesale, the arithmetic connecting the mechanism (10^-5–10^-3 M_sun from 800–2000 flares) to the fitted kilonova ejecta (10^-2 M_sun for both bursts) does not close. The beaming-factor argument misapplies a relativistic gamma-ray beaming correction to non-relativistic quasi-spherical ejecta, and the fallback of ~10-1000x unseen flares is in tension with the paper's own interpretation of tau_min as the single-spike duration, and pushes toward unbinding a large fraction of the entire crust. This is an internal-consistency problem rather than a disagreement with consensus, and it is directly quantifiable with existing published inputs (Patel et al. yields, Cehula et al. ejecta scalings, the observed light curves). I do not move the verdict off CONDITIONAL: the reader already assigned high correctness risk and flagged that the Mej claim is not secure; my concern sharpens where that risk lives (the multiplicity/beaming closure, not only the crust-temperature caveat) and yields a falsifiable numerical check, but it does not independently demonstrate the framework is wrong — a generous reading (higher per-flare ejecta at 100x GF energies, modest geometric correction, magnetar boost dominating the light curve) could narrow the gap. If the concrete test shows a persistent >=10x shortfall, the r-process pillar of the \"unified picture\" should be downgraded to a hypothesis and the verdict revisited toward REJECT for the nucleosynthesis claim specifically, while the prompt/plateau unification could still stand.","tokens_in":18084,"tokens_out":2854,"duration_ms":99026,"concrete_test":"Do the self-consistent budget the paper skips. Take the Patel et al. (2025b) / Cehula et al. (2024) ejecta-mass scaling as a function of flare energy, evaluated at this paper's own inferred spike energies (E_iso ~ 6.6x10^48 and 1.6x10^49 erg), and compute the maximum flare multiplicity N_max consistent with the observed light-curve morphology (spikes of width tau_min must not overlap into a continuum, and far-side-only flares contribute at most a modest geometric factor). Then compare N_max x m_ej(E_flare) against the fitted 1.2x10^-2 M_sun and against the available crustal mass above the shock-coupling depth. If N_max x m_ej falls short by more than ~a factor of a few — or if reaching 10^-2 M_sun requires ejecting most of the crust or cadences below tau_min — the r-process/kilonova leg of the unification fails quantitatively even granting the cold-crust assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires the GF-forest mechanism to supply the ejecta that the joint fit actually measures. Section 3.3/Table 2 fit log10(Mej/M_sun) = -1.91 for BOTH bursts, i.e. ~1.2x10^-2 M_sun. The mechanism-side accounting is: m_ej ~ 10^-8–10^-6 M_sun per flare (Cehula et al. 2024, Table 1) times ~800 spikes (211211A) or ~2000 spikes (230307A), where the spike counts come from dividing T90 by tau_min (Section 3.1). That gives 10^-5–10^-3 M_sun — one to three orders of magnitude short of the fitted value. The paper closes this gap in two steps, and both are problematic. (1) It invokes the monopolar beaming factor f_b ~ 10^-5–10^-3 to argue Mej is a \"lower limit\" (end of Section 2, and point 4 of Section 3.3). But f_b is the relativistic beaming of the gamma-ray SPIKES; the ejecta parcels are non-relativistic (v ~ 0.1–0.4c) and quasi-isotropically launched from the crust. The number of flares missed because their beamed gamma-ray emission points away is a geometric covering-fraction correction bounded by the fraction of the closed-field-line region invisible to us — plausibly a factor of a few, not 1/f_b ~ 10^3–10^5. (2) Point 4 concedes the need for \"at least ten times more misaligned GF events,\" i.e. N >= 10^4–10^6 flares over 13–18 s, an inter-flare cadence of ~10 us–1 ms. That contradicts the paper's own spike census: tau_min ~ 9–16 ms is interpreted as the duration of a single spike, and at 10^4+ flares the light curve would be a smooth continuum of overlapping spikes, not a \"forest\" with 10-ms resolved variability — unless one assumes the extra flares are all far-side events that eject mass but contribute no flux, in which case the geometric correction is again only modest. There is also a crustal mass-budget tension: ~10^-2 M_sun is a substantial fraction (roughly 20–100%) of the entire NS crust, to be unbound by surface-localized fireball shocks in ~15 s while the star continues to operate as a coherent magnetar engine, with the crust needing to reform (","agreement_with_reader":"partial"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The paper proposes a unified model for peculiar long GRBs (211211A, 230307A) based on tidally disrupted NS–WD mergers with a pre-merger magnetar. Hyperaccretion amplifies the toroidal field to B_phi ~ 1e15–5e16 G, triggering repeated magnetar giant flares whose initial spikes form the main burst as a \"forest\" of ~800–2000 spikes (N = T90/tau_min). The accretion-to-propeller transition produces the EE and the MB–EE trough; post-merger magnetar spin-down powers the X-ray plateau and boosts the kilonova. The key new ingredient is r-process production: each flare's e±–γ fireball shocks the crust and unbinds mej ~ 1e-8–1e-6 Msun (per Cehula et al. 2024), synthesized via α-rich freeze-out; the ensemble yields Mej ~ 1e-5–1e-3 Msun, argued to suffice for the kilonova when boosted by spin-down and beaming corrections. Joint afterglow+kilonova MCMC fits reproduce both events' light curves, but return Mej ~ 1.2e-2 Msun for both bursts, which the authors attribute to unobserved misaligned flares.","tokens_in":18658,"tokens_out":6986,"duration_ms":224062,"significance":"If the scenario holds, it would provide the first unified explanation of peculiar LGRBs — prompt three-episode structure, plateau, and kilonova — while identifying a concrete, observationally accessible r-process site (giant-flare crustal ejecta), which would be of broad interest. The manuscript has real strengths: joint MCMC fits to the full multi-wavelength datasets of both bursts, explicit per-spike energetics tied to an inferred toroidal field, a specific MB–EE trough interpretation as the accretion-to-propeller transition, falsifiable elements (spike-width statistics tied to Galactic GFs; predicted flare multiplicity), and an unusually candid caveat section. These give the model testable content beyond a phenomenological fit. However, the headline nucleosynthesis claim currently rests on unquantified extrapolations and a beaming argument that is internally inconsistent with the paper's own spike census, so the significance is conditional on resolving the major comments below.","major_comments":[{"comment":"The central r-process claim fails its own accounting. The mechanism supplies mej~1e-8–1e-6 Msun per flare x 800–2000 spikes = Mej~1e-5–1e-3 Msun (§2, Panel c), but Table 2 fits log Mej = -1.91 (~1.2e-2 Msun) for both bursts. The bridge offered — 'at least ten times more misaligned GF events' and the f_b~1e-5–1e-3 beaming factor — misapplies f_b, which is the relativistic beaming of the gamma-ray spikes; the ejecta parcels are non-relativistic (v~0.1–0.4c) and launched quasi-isotropically from the crust, so the unseen-flare correction is bounded by the covering fraction of the closed-field-line region (a factor of a few, not 1/f_b~1e3–1e5). Moreover, >=10x more flares means N>~1e4 over 13–18 s, i.e. ~1 ms cadence, below the tau_min~9–16 ms that §3.1 identifies with a single spike; the forest would overlap into a smooth continuum, contradicting the paper's own spike census. The summary sta","section":"§3.3, point 4; §2 (penultimate paragraph); Table 2"},{"comment":"The per-spike field estimate (E_spike ~ B_phi^2 DeltaR^3/8pi, giving B_phi ~ few x 1e16 G) ignores the integrated budget. The true total gamma-ray energy of the flare forest is ~5.3e51 erg (211211A) and of order 1e52 erg (230307A) — independent of f_b, since N_true x (f_b E_iso,spike) = E_iso,MB. A single toroidal reservoir at B~5e16 G contains only ~7e50 erg within the NS volume, so the field must be dissipated and re-amplified ~10–40 times within ~10–20 s, requiring a mean re-amplification power ~1e51–1e52 erg/s. This should be compared explicitly with the accretion-driven dynamo of Zhong et al. (2023); without it, the per-spike B estimate and the 'thousands of flares' picture are not shown to be energetically self-consistent.","section":"§3.1"},{"comment":"The entire per-flare yield (mej~1e-8–1e-6 Msun, Ye>~0.4, alpha-rich freeze-out) rests on the Cehula et al. (2024) cold, non-accreting, beta-equilibrated crust. Here the NS hyperaccretes at >~1e-2 Msun/s during the flaring epoch and fires ~1e3 flares in ~10 s; the paper itself concedes (§4) that polar shock heating and repeated ejection outpacing beta-equilibrium recovery are unmodeled. Because the ejecta mass and composition scale directly with these assumptions, this is load-bearing for the r-process claim rather than a peripheral caveat. The authors should quantify the sensitivity (e.g., how mej and the yield pattern degrade for a heated or depleted crust), or explicitly make the abstract's 'collectively yields Mej >~ 1e-5–1e-3 Msun' conditional on this assumption.","section":"§2, Panel (c); §4 caveats"},{"comment":"The plateau analysis (Zhong et al. 2023, 2024) and the kilonova boost yield different (P0, B_s,post) for the same post-merger magnetar, reconciled here by a newly posited angular structure of the spin-down radiation — a collimated component powering the plateau and a quasi-isotropic component boosting the kilonova. This is introduced to absorb a tension, not derived, and no fraction f_iso of the spin-down power is specified or checked against the total rotational energy for the fitted P0 and B_s,post. Since the plateau and the magnetar-boosted kilonova are both pillars of the 'unified' claim, the authors should quantify f_iso, demonstrate energetic consistency with both components, and state a testable consequence (e.g., off-axis plateau behavior).","section":"§3.3, point 2"}],"minor_comments":[{"comment":"The expression f_b ≃ 0.5 x (1/2) theta_j^2 contains an unexplained extra factor 0.5 relative to the standard monopolar beaming factor f_b = (1 - cos theta_j) ≈ theta_j^2/2. With the fitted theta_j = 0.028–0.072 rad (Table 2), the standard expression gives 4e-4–2.6e-3, not the quoted 1e-5–1e-3; please clarify the definition and its role in the ejecta-mass argument.","section":"§2, penultimate paragraph"},{"comment":"The spike counts N = T90/tau_min (~800 and ~2000) assume a 100% duty cycle of back-to-back spikes; since tau_min is strictly a lower bound on the spike width, these N values are upper limits and the resulting Mej range is likewise an upper envelope, which should be stated. Please also specify how tau_min was measured for these bursts.","section":"§3.1"},{"comment":"Roming et al. (2005) is cited for the Swift XRT light curves of GRB 211211A; the XRT repository reference is Evans et al. (2009) (already in the list), while Roming et al. is the UVOT instrument paper.","section":"§3.3, first paragraph"},{"comment":"With ~14 free parameters, the paper shows only best-fit curves and a parameter table. Please add MCMC corner plots (or at least prior ranges) and a residuals panel, and comment on why both bursts return identical best-fit log Mej = -1.91 — coincidence, or a prior/boundary effect.","section":"Fig. 2 and Table 2"},{"comment":"The fitted opacity for GRB 211211A, kappa = 0.88 cm^2/g, lies below the 3–20 cm^2/g range quoted in §3.2 and implies lanthanide-poor ejecta. This sits awkwardly with the r-process interpretation advocated here (versus the earlier 56Ni interpretation of Zhong et al. 2023) and deserves explicit discussion, as is done for GRB 230307A in point 3.","section":"Table 2"},{"comment":"'...will be addressed in forthcoming work by Cehula et al. (2024)' — presumably a forthcoming paper by those authors is meant; the reference and tense should be fixed.","section":"§4, final caveat"},{"comment":"The photospheric radius is set to vt with no thermalization-depth or electron-scattering correction; this is standard but should be stated as an assumption, since it affects the late-time color evolution shown in Fig. 2.","section":"Eq. (7)"},{"comment":"Check the author-name convention for 'Qiumu et al. (2026)' — the surname appears to be listed first.","section":"§1"}],"recommendation":"major_revision","confidential_remarks":"The engine architecture, accretion profiles, toroidal amplification, and propeller interpretation are all imported from the authors' own Zhong et al. (2023, 2024) and re-applied to the same two bursts, with discrepancies absorbed by new, untested hypotheses (structured spin-down radiation, unseen flare populations). This makes the 'unified picture' considerably less prior-independent than the abstract suggests. The paper is publishable in principle if the ejecta-mass gap is confronted honestly rather than bridged by the beaming argument; I encourage the editor to insist on that specific revision."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new move here is welding Zhong’s 2023 NS–WD accretion/propeller engine to the Cehula/Patel GF fireball–crust ejecta channel so the same forest of accretion-triggered giant-flare spikes both builds the main burst and supplies the missing r-process. That is a real compositional step, not just a re-label. The prompt structure (MB forest, propeller EE, trough as phase switch), the X-ray plateau, and the joint afterglow+kilonova MCMC fits are cleanly executed and reproduce the two flagship light curves with parameters that sit in the expected ballpark.\n\nEnergetics check out at order of magnitude: spike energies imply Bϕ a few ×10^16 G, consistent with their earlier amplification curves, and τ_min ~10 ms matches known GF initial-spike widths once the long tails are cut. The authors also flag the two biggest physics caveats themselves (hot accreting crust, repeated ejections outrunning β-equilibrium), which is honest.\n\nThe soft spot that actually matters is the mass budget. Mechanism side gives Mej ~10^{-5}–10^{-3} M⊙ from ~800–2000 parcels. Fit side wants ~10^{-2} M⊙ for both bursts. Invoking the relativistic beaming factor of the gamma spikes to multiply the non-relativistic crustal ejecta is not valid geometry; a modest covering-fraction correction (factor of a few) is all you can claim, and “ten times more far-side flares” already starts to erase the resolved 10 ms variability the paper itself uses as the spike census. Crustal mass fraction also gets uncomfortable at 10^{-2} M⊙. Low fitted κ for 211211A is a lesser tension if one still wants third-peak material.\n\nCircularity is moderate, not fatal: same group’s engine, same two events, new structured-spin-down hypothesis to reconcile magnetar parameters. No code, but the semi-analytic setup is re-implementable.\n\nThis is for people who already care about peculiar LGRB progenitors and alternative r-process sites. It is a serious theoretical framework, not a finished hydro+nucleosynthesis calculation. I would send it to referees; the load-bearing assumptions are clear enough that a good referee can force the needed follow-up. Worth a reading-group slot if the group works multi-messenger or magnetar engines. I would cite it as the current best single-engine NS–WD attempt while noting the Mej gap.","headline":"Coherent single-engine NS–WD picture that fixes the r-process hole with GF crustal ejecta, but the fitted Mej sits an order of magnitude above what the spike census can deliver and the beaming patch does not cleanly close it.","tokens_in":19689,"tokens_out":641,"would_cite":true,"duration_ms":20409,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Peculiar long GRBs arise from NS–WD mergers that fire repeated magnetar giant flares, making both the burst and the r-process kilonova.","keywords":["gamma-ray bursts","magnetars","neutron star–white dwarf mergers","r-process nucleosynthesis","kilonovae","giant flares","accretion disks","X-ray plateaus"],"falsifier":"A high-resolution spectroscopic campaign on a future peculiar long GRB that either fails to detect the expected second- and third-peak r-process lines once the afterglow has faded, or finds an ejecta mass and velocity structure incompatible with thousands of discrete crustal parcels.","tokens_in":19038,"feed_emoji":"💥","tokens_out":1051,"duration_ms":21616,"temperature":0.7,"pith_summary":"Peculiar long gamma-ray bursts such as GRBs 211211A and 230307A show a hard main burst, a soft extended emission separated by a trough, an X-ray plateau, and a kilonova with r-process signatures. The paper argues that these events are tidally disrupted mergers of a pre-merger magnetar with a massive white dwarf. Hyperaccretion from the constant-entropy disk amplifies the neutron star’s toroidal field until repeated giant flares erupt; their initial spikes form the main burst, while the later propeller phase produces the extended emission. Each flare’s electron–positron–photon fireball shocks the crust, ejecting neutron-rich material that undergoes α-rich freeze-out and builds enough r-process ejecta to power the observed kilonova once boosted by magnetar spin-down, which also lights the X-ray plateau. The picture therefore unifies prompt emission, afterglow, kilonova, and heavy-element production inside a single NS–WD engine.","feed_headline":"Magnetar flares in NS–WD mergers make peculiar long GRBs","feed_subtitle":"One engine supplies the burst forest, the trough, the X-ray plateau, and the r-process kilonova","key_machinery":"Accretion-amplified magnetar giant flares: hyperaccretion drives the toroidal field to ~10^15–5×10^16 G, triggering repeated starquakes whose e±–γ fireballs both radiate the main-burst spikes and unbind successive parcels of crust (each ~10^{-8}–10^{-6} M⊙) that freeze out into r-process nuclei.","core_discovery":"A tidally disrupted neutron-star–white-dwarf merger with a pre-merger magnetar (surface field ≳10^14 G) and a massive white dwarf produces the entire observed phenomenology of peculiar long GRBs: the main burst is a forest of accretion-amplified giant-flare initial spikes, the extended emission and the trough arise from the accretion-to-propeller transition, crustal ejecta from the associated fireballs synthesize r-process nuclei via α-rich freeze-out, and post-merger spin-down powers both the X-ray plateau and an extra boost to the kilonova.","pith_inferences":["If repeated flares deplete or heat the crust faster than β-equilibrium can be restored, the r-process yield per spike should decline with time, imprinting a measurable spectral evolution across the main burst.","The model predicts that MeV gamma-ray lines from freshly synthesized r-process nuclei could appear during or shortly after the main burst, offering a prompt nucleosynthesis diagnostic.","Similar accretion-amplified giant-flare forests may operate in other high-accretion magnetar systems, potentially linking peculiar GRBs to a broader class of magnetar-driven transients."],"forward_implications":["The main-burst minimum variability timescale directly counts the number of giant-flare spikes and therefore lower-limits the cumulative crustal ejecta mass.","The MB–EE trough is a clean observational signature of the accretion-to-propeller transition and can be used to diagnose the accretion-rate history.","Post-merger magnetar parameters inferred from the X-ray plateau and from the kilonova boost need not be identical if the spin-down luminosity is angularly structured.","Jet beaming implies the true r-process yield can exceed the line-of-sight estimate by an order of magnitude or more.","The same engine can explain both events that show spectroscopic r-process features and those that so far show only photometric kilonova candidates."],"fun_headline_variants":["NS–WD mergers amplify magnetar flares into peculiar long GRBs","Giant-flare forests from NS–WD disks unify odd long GRBs","Magnetar GFs in tidally disrupted NS–WD mergers make long GRBs","Accretion-driven magnetar flares explain peculiar long GRB signals","NS–WD merger flares forge r-process ejecta and long GRB parts"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The cold, non-accreting crust calculation still holds under sustained hyperaccretion and thousands of successive flares, so each spike continues to unbind enough neutron-rich crust for α-rich freeze-out.","fun_headline_variants_meta":{"raw":{"variants":["NS–WD mergers amplify magnetar flares into peculiar long GRBs","Giant-flare forests from NS–WD disks unify odd long GRBs","Magnetar GFs in tidally disrupted NS–WD mergers make long GRBs","Accretion-driven magnetar flares explain peculiar long GRB signals","NS–WD merger flares forge r-process ejecta and long GRB parts"]},"model":"grok-4.5","effort":"low","cost_usd":0.004291,"raw_usage":{"total_tokens":1437,"prompt_tokens":974,"num_sources_used":0,"completion_tokens":81,"cost_in_usd_ticks":42908000,"prompt_tokens_details":{"text_tokens":974,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":382,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":974,"tokens_out":81,"duration_ms":7081,"temperature":1.0,"reasoning_tokens":382,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T12:18:36.082747+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A high-resolution spectroscopic campaign on a future peculiar long GRB that either fails to detect the expected second- and third-peak r-process lines once the afterglow has faded, or finds an ejecta mass and velocity structure incompatible with thousands of discrete crustal parcels.","supporting_citations":[],"review_version":1}