{"id":"607bf2e6-c505-464e-a2ae-57f8b15780b3","arxiv_id":"2603.11718","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":13,"one_line_summary":"Finite ejecta thickness, not off-axis geometry, explains the achromatic afterglow peaks of XRF 080330 and GRB 080710, implying ~10^13 cm ejecta widths and longer central engine activity.","lead":"The authors re-analyze the early multi-wavelength afterglows of two gamma-ray events and argue that their achromatic peaks are produced by the finite thickness of the relativistic ejecta, not by viewing the jet off-axis. If correct, early afterglow light curves encode the duration of the central engine's activity and the density structure around the progenitor.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"GRB 080710 primary result depends on a posteriori 10× inflation of decay-phase errors; without it the model fails to reproduce the rise, so the joint Δ0~10^13 cm and finite-thickness claim for both bursts is conditional.","rationale":"The reader's weakest_assumption identifies exactly the issue I find most load-bearing: the primary inference setups are selected after looking at the fits, with data-weighting factors chosen to produce physically plausible or rising-phase-matching results. For GRB 080710, this selection is extreme (a factor-of-10 inflation of decay-phase uncertainties) and is the sole reason the model reproduces the rise that the central claim interprets. The paper is transparent about this, which is good, but transparency does not remove the conditionality. The central claim is stated for both events, and the GRB 080710 part is unsupported without a priori justification of the weight choice or a demonstration that conclusions are insensitive to it. I agree with the reader that the result is conditionally acceptable: the authors should justify the weighting a priori and test robustness. My stress-test does not change the verdict; it reinforces the conditional nature. I also considered a more physics-specific concern—that Δ0 may not be identifiable in the thin-shell regime for GRB 080710—but the paper's posterior does show a narrow Δ0, and the data-weighting selection is the more direct and damaging issue. The XRF 080330 Δ0 constraint appears robust across setups, so a full rejection would be too strong; CONDITIONAL (as the reader already concluded) is appropriate. The concrete test I propose would settle whether the GRB 080710 result is an artifact of the weighting: rerun with original uncertainties and check if the rise is still reproduced and the Δ0 posterior remains narrow.","tokens_in":32675,"tokens_out":6959,"duration_ms":65375,"concrete_test":"Re-run the GRB 080710 Bayesian inference using the original published 1σ uncertainties (no factor-of-10 inflation) with all optical/NIR and X-ray points included, the same priors, and the same model. If the posterior no longer produces a narrow Δ0~10^13 cm and β<1, or if the rising phase is not reproduced, the primary result is an artifact of the selected downweighting. As a secondary check, fix Δ0 to a negligible value (e.g., 10^10 cm) and see whether the rise can still be fitted, testing identifiability in the thin-shell regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that finite-thickness jet dynamics, not off-axis viewing, explains the achromatic peaks in both events rests most heavily on the GRB 080710 analysis. In §4.2, the authors state that 'in almost all cases we examined, the model failed to adequately reproduce the rising phase,' and they select as primary the setup obtained by 'conservatively inflating the uncertainties of the decaying-phase data by a factor of ten.' This is an explicitly a posteriori choice: the data weighting is adjusted until the model matches the feature (the rise) that the paper then interprets physically. With the original uncertainties, the model does not reproduce the rise; with the decay phase effectively removed from the likelihood, the rise is fitted and MAP values (E0~2×10^54 erg, Γ0~42, Δ0~1.3×10^13 cm, k~0.06, β~0.64) are quoted. The quoted χ²/d.o.f. for the full data set is 5.4, while the rise-only χ² is 2.2—a clear sign that the fit is only acceptable after discarding most of the data. Consequently, the inferred Δ0 and the conclusion that the peak is explained by finite-thickness dynamics rather than off-axis effects are conditional on this arbitrary weighting. The XRF 080330 analysis has the same structure, with X-ray uncertainties inflated by a factor of 2 to avoid energetically implausible E0~10^55 erg, though the Δ0 constraint there appears more robust because it is tied to the optical/NIR achromatic break. The paper's own admission that the model 'failed to adequately reproduce the rising phase' in other setups means that the central claim for GRB 080710 is not a robust, data-driven inference but a selection effect. This is the load-bearing weak point: if the weighting is not justified a priori, the central claim for both bursts collapses to a single-event claim for XRF 080330.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reanalyzes the early afterglows of XRF 080330 and GRB 080710 using the public numerical code Magglow, which incorporates finite-thickness ejecta and a generalized power-law circumburst density profile, within a Bayesian inference framework (MultiNest). The central claim is that the achromatic peaks/breaks at ~10^3–10^4 s in both events are best explained by jet dynamical evolution with finite shell thickness, with inferred initial radial widths Δ0 ~ 8.6×10^12 cm (XRF 080330) and ~1.3×10^13 cm (GRB 080710), rather than by off-axis viewing. The paper further claims that a free density slope k is strongly favored over canonical ISM (k=0) or wind (k=2) models, and that the implied central engine activity timescale Δ0/c is about an order of magnitude longer than the prompt T90/(1+z). The analysis includes posterior predictive light curves, model comparison via Bayesian evidence, and a falsifiable radio-prediction section.","tokens_in":33158,"tokens_out":3101,"duration_ms":31527,"significance":"If the central inference is robust, the paper would strengthen the case that finite-thickness ejecta dynamics are important for early afterglow interpretation and would provide a physical bridge between prompt and afterglow phases. The strengths of the manuscript are substantial: it uses a publicly available numerical code with self-consistent shell-thickness evolution, performs full Bayesian parameter estimation with nested sampling, reports posterior distributions and credible intervals, makes explicit model comparisons via evidence, and gives falsifiable radio predictions (§5.6) that are not used to set the model constants. The authors are also transparent about the main systematic limitation — the infinitesimally thin emission region (§5.1). These features make the paper valuable regardless of the outcome of the specific data-weighting choices, provided those choices are justified or their impact quantified.","major_comments":[{"comment":"The primary result for GRB 080710 is obtained by 'conservatively inflating the uncertainties of the decaying-phase data by a factor of ten.' The authors state that 'in almost all cases we examined, the model failed to adequately reproduce the rising phase,' so the setup is selected after inspecting the fits. Under the original uncertainties the model does not reproduce the rise; after the decay phase is effectively removed from the likelihood, the rise is fitted and the MAP values are quoted. The reported full χ²/d.o.f. = 5.4, versus 2.2 for the rise-only subset, shows that the fit is acceptable only after discarding most of the data. This makes the central claim — that the achromatic peak in GRB 080710 is explained by jet dynamics rather than off-axis effects — conditional on an a posteriori weighting choice. A pre-specified weighting rule, or a demonstration that the conclusion is inse","section":"§4.2 (GRB 080710, especially the paragraph beginning 'We performed Bayesian inference under several different conditions"},{"comment":"The X-ray uncertainties are inflated by a factor of two because with the original uncertainties the fit drives E0 to ~10^55 erg and θ_j > 0.1 rad, which the authors deem 'energetically implausible.' This is also a post hoc data‑weighting choice, though the authors argue that Δ0 and k are robust across setups. To make the primary result and the Bayesian model comparison (ΔlnZ > 20 favoring free k) credible, the paper should quantify how the posteriors of Δ0, k, and the evidence change when the original X-ray uncertainties are used, or when the X-ray data are treated as upper limits. Without this, the Δ0 ~ 10^13 cm conclusion for XRF 080330 is not yet shown to be independent of the selection.","section":"§4.1.1 (X-ray down-weighting for XRF 080330)"},{"comment":"The model treats the emission region as infinitesimally thin, and the authors point out that this systematically underestimates the cooling frequency and hence the X-ray flux. This is an honest caveat, but it undermines the use of the X-ray data in the same fit — the X-ray down-weighting is presented as a remedy, yet the optical/NIR fit itself has χ²/d.o.f. = 3.5 for XRF 080330 and 5.4 (full) for GRB 080710. Those values indicate systematic model–data discrepancies beyond the reported error bars. The claim that the finite-thickness treatment resolves the achromatic-peak question would be more convincing if the remaining optical/NIR discrepancies were addressed, or if the posterior predictive bands were shown to contain the data at the quoted χ² level. I recommend presenting the best-fit residuals and a χ² decomposition by band.","section":"§5.1 (thin emission region) and §4.1.2/§4.2.1 (χ²/d.o.f.)"},{"comment":"The one-dimensional posterior of log(Δ0/c) for GRB 080710 is bimodal, but the authors interpret this as 'an apparent effect arising from the geometry of the parameter space' without quantitative support. Since Δ0 is a central quantity, the bimodality should either be analyzed (e.g., by showing the two modes correspond to different physical branches of the model, or by presenting the posterior conditioned on p and β) or demonstrated to be a marginalization artifact through a likelihood-ratio or profile-likelihood calculation. As written, the claim that Δ0 = 1.3×10^13 cm is the representative value is not fully justified.","section":"§4.2.1 (bimodality of log(Δ0/c))"}],"minor_comments":[{"comment":"The jet-break time formula uses 'θ_obs ± θ_j' with a brief explanation, but the sign convention is not fully tied to the text. In §4.1.2 and §4.2.1, it would help to state explicitly which branch is used for each event and how the plus/minus choice affects the quoted T_jet values.","section":"Eq. (24)"},{"comment":"The captions say 'black dashed lines show the medians and the 95% credible intervals,' but the text says 'black dashed lines show the medians,' and the figure descriptions in the main text sometimes call them 'dashed lines.' Please make the caption consistent with the actual plot.","section":"Fig. 2 and Fig. 4 captions"},{"comment":"The column header reads 'T90/(1+z)' but the values (24 s, 65 s) are rest-frame durations; the comparison to Δ0/c is then only meaningful if the units are made explicit. Suggest writing 'cT90/(1+z) [s]' to match the text and to avoid implying that the gamma-ray duration itself is a length.","section":"Table 3"},{"comment":"The normalization A(n0,k) = n0 (3×10^35)^(k/2) and the statement 'n0 has units [cm^{k-3}]' are confusing; for k=2 the text says n0 is in cm^{-1}, for k=1 in cm^{-2}. Please clarify the units of n0 consistently, ideally by writing n0 as a density normalization at a reference radius (e.g., R = 10^17 cm).","section":"Equations (1) and footnote 1"},{"comment":"The 'local' prompt efficiency uses f_rad = cT90/Δ0, but the numerical values are not stated for the two events (only Δ0/c is in Table 3). Including f_rad and the resulting η_γ,local explicitly would help the reader check the claimed factor-of-~10 reduction.","section":"§5.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and uses modern Bayesian tools with a public code, and the authors are transparent about the thin-emission-region approximation. The central problem is statistical: the primary conclusions for both events depend on data-weighting factors that are introduced after seeing the fits. The GRB 080710 case is especially concerning because the model fails to reproduce the rise unless the decay-phase data are down-weighted by a factor of ten, and even then the full χ²/d.o.f. is 5.4. This is a load-bearing issue rather than a presentational one, so I recommend major revision. The authors could address it by: (i) pre-registering or systematically varying the down-weighting factors and showing the sensitivity of Δ0, k, and the off-axis/on-axis conclusions; (ii) reporting the fit with original uncertainties, even if it yields an energetically implausible E0, and discussing what physical constraint (e.g., E_K,jet limits) would rule that out; (iii) providing profile likelihoods or posterior conditionals to resolve the Δ0 bimodality in GRB 080710. If those are provided, the paper would be a solid contribution to the early-afterglow modeling literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look, but the two-event claim is really a one-event claim plus a conditional fit. The solid piece is XRF 080330: the achromatic optical/NIR break at ~1.5e3 s maps onto T_BM and pins Delta0 ~ 8.6e12 cm across all their inference setups, with k ~ 0.9 and strong Bayesian preference for free k over ISM/wind (Delta ln Z > 20). That is a genuine new result, and the transparency is a credit: they state the X-ray downweighting, the thin-emission-region caveat, and the radio predictions.\n\nThe soft spot is GRB 080710. The paper admits the model failed to reproduce the rise in almost all setups, and the 'primary' result is selected after inflating decay-phase uncertainties by a factor of 10. With full weights, chi2/dof is 5.4; the rise-only fit is 2.2. That is not an inference, that is a selection effect. So the joint conclusion that finite-thickness dynamics explain both events does not hold as stated. The Delta0 for GRB 080710 (~1.3e13 cm) is conditional on that weighting. The XRF 080330 result stands on its own.\n\nAlso, the off-axis disfavoring is presented strongly, but they never run a dedicated off-axis structured-jet model in the same Bayesian framework. Their posterior beta<1 is within their finite-thickness top-hat model, so it says the finite-thickness model prefers on-axis, not that off-axis is ruled out. The paper's own discussion admits structured jets could matter.\n\nThe Bayesian evidence comparison is only within the same model family (finite-thickness ejecta with k free vs fixed), so it doesn't validate the family itself.\n\nStill, the work is honest, the code is public, the limitations are spelled out, and the XRF 080330 constraint is a real step forward. If you work on early afterglows, you should know this. I'd send it to peer review, but the referee should push for a robustness analysis of the GRB 080710 weighting and a comparison against an off-axis model.","headline":"Worth a referee, but the two-event claim is really one solid constraint (XRF 080330) plus a conditional fit for GRB 080710 that leans on a post hoc 10x error inflation.","tokens_in":33736,"tokens_out":2464,"would_cite":true,"duration_ms":25192,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.70.Rz"],"model":"deepseek-v4-flash","headline":"Finite-thickness jet dynamics, not off-axis viewing, produce the late achromatic peaks in XRF 080330 and GRB 080710; the implied shell width of ~10^13 cm suggests the central engine ran about ten times longer than the gamma-ray burst.","keywords":["gamma-ray bursts","afterglow light curves","finite-thickness ejecta","achromatic peaks","X-ray flashes","circumburst density profile","Bayesian inference","relativistic jets"],"falsifier":"A direct check: rerun the Bayesian inference with the original, unaltered X-ray uncertainties for XRF 080330 and decay-phase uncertainties for GRB 080710, with the data weights fixed before any fitting. If the fit again drives the isotropic energy to about 10^55 erg (the solution the authors rejected as energetically implausible) or fails to reproduce the rising phase, then the quoted shell width of 10^13 cm and the density slopes are artifacts of the post-hoc weighting. An independent observable test: at about 10^5 s the model predicts a 1.4 GHz radio afterglow near 0.1 mJy for XRF 080330-lik","tokens_in":32559,"feed_emoji":"💥","tokens_out":8534,"duration_ms":70233,"temperature":0.7,"pith_summary":"This paper tries to establish that the achromatic peaks seen a few thousand seconds after two gamma-ray bursts, XRF 080330 and GRB 080710, are produced by the forward shock's own deceleration history once the ejecta's finite radial thickness is taken into account, rather than by the off-axis viewing geometry invoked in the original discovery papers. Using Bayesian inference over an afterglow model with a free circumburst density slope, the authors find a shell width of roughly 10^13 cm in both events, about an order of magnitude larger than the size implied by the prompt gamma-ray duration. If correct, this links the prompt and afterglow phases: the central engine kept working well after the gamma-rays stopped, and early afterglow light curves cannot be safely read with the thin-shell approximation. The analysis also favors a generalized density profile over the canonical uniform or steady-wind models, with XRF 080330 preferring a slope near k=1 and GRB 080710 an approximately uniform medium.","feed_headline":"Jet thickness, not viewing angle, shapes two GRB afterglow peaks","feed_subtitle":"A 10^13 cm shell width implies the central engine ran about ten times longer than the prompt burst.","key_machinery":"The central object is the finite-thickness ejecta shell with initial radial width Delta_0 and the three-phase dynamics it produces: free expansion at constant Lorentz factor, a transition phase in which the forward-shock Lorentz factor falls as R^-(2-k)/4, and the eventual approach to the Blandford-McKee self-similar deceleration. The load-bearing dimensionless parameter is xi_k = (l_S/Delta_0)^(1/2) Gamma_0^-(4-k)/(3-k), which conventionally separates thick-shell (xi<1) from thin-shell (xi>1) behavior but, as the paper stresses, does not by itself fix the ordering of the observer-frame timescales T_tr, T_gamma, and T_BM. The mechanism that carries the argument is the delay of the transition","core_discovery":"On the paper's own terms, the discovery is that the gradual rise and achromatic break in the multi-wavelength afterglows of XRF 080330 and GRB 080710 are fingerprints of the ejecta's transition from free expansion to the Blandford-McKee deceleration phase, with the timing set by the finite initial shell width. For XRF 080330 the event sits in the thick-shell regime (xi_k = 0.10), with the end of the transition phase at about 1.8x10^3 s matching the achromatic break; for GRB 080710 the comparable ordering of the transition, deceleration, and Blandford-McKee onset timescales around 2x10^3 s, produced by a lower Lorentz factor, shapes the peak in a thin-shell-like regime even though the shell i","pith_inferences":["Editorial inference: the fit statistics reported in the paper (chi^2 per degree of freedom near 4-5 for the full datasets) show that the forward-shock model does not fully capture the X-ray band and the late optical decay; the paper's conclusion that the shell width and density slope are unaffected by this tension assumes the excess is a separate component that does not correlate with the peak-tim","Editorial inference: a reanalysis with structured-jet or reverse-shock components could revise the GRB 080710 result, since the paper itself notes that the top-hat jet overproduces the late-time flux; in particular, the near-uniform density slope inferred for this burst may be an artifact of the top-hat assumption.","Editorial inference: applied to the growing sample of X-ray flashes now being discovered, the same analysis pipeline could turn early afterglow rise times into a population measurement of shell width, testing whether the roughly tenfold ratio between engine activity and prompt duration seen in these two events is universal.","Editorial inference: the near-on-axis solutions for both an X-ray flash and a classical burst suggest that achromatic late peaks are not a viewing-angle diagnostic; if confirmed, the circumburst density slope k, rather than the peak morphology, is the more informative discriminator of progenitor history."],"forward_implications":["If the central claim is right, the achromatic peaks in XRF 080330 and GRB 080710 are not evidence of off-axis jets; the viewing geometry in both events is nearly on-axis, so the earlier off-axis interpretations of these two bursts are not required.","The constraint on the shell width, about 10^13 cm for both bursts, turns the afterglow peak time into a direct measurement of central engine activity of roughly 300-470 s in the rest frame, about ten times longer than the prompt gamma-ray duration, implying a non-uniform radial structure of the ejecta.","Early afterglow light curves computed under the thin-shell approximation will misplace the deceleration onset and the achromatic break; finite-thickness dynamics must be included when interpreting early-time data.","Fixing the external density profile to a uniform medium (k=0) or a steady wind (k=2) is strongly disfavored for XRF 080330 by the Bayesian evidence, and the inferred slopes (k near 1 for XRF 080330, near 0 for GRB 080710) point to diversity in the progenitor's terminal mass-loss history.","The predicted 1.4 GHz radio afterglow differs sharply between the two events, roughly 0.1 mJy peak for XRF 080330 and 10^-3 mJy for GRB 080710, making radio follow-up a concrete way to break degeneracies left by optical and X-ray data."],"fun_headline_variants":["Thick ejecta, not viewing angle, sets GRB afterglow peaks","10^13 cm shell width drives early GRB afterglow peaks","GRB afterglows reveal engine ran longer than prompt burst","Finite shell thickness ties GRB prompt to afterglow","Bayesian fits favor thick ejecta over off-axis for GRB peaks"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central claim rests on a data-weighting decision made after looking at preliminary fits: Section 4.1.1 doubles the X-ray uncertainties for XRF 080330 and Section 4.2 inflates the decay-phase uncertainties tenfold for GRB 080710, so that the primary solutions are the ones that yield energetically plausible parameters and reproduce the rising phase, rather than the fits with unaltered weights.","fun_headline_variants_meta":{"raw":{"variants":["Thick ejecta, not viewing angle, sets GRB afterglow peaks","10^13 cm shell width drives early GRB afterglow peaks","GRB afterglows reveal engine ran longer than prompt burst","Finite shell thickness ties GRB prompt to afterglow","Bayesian fits favor thick ejecta over off-axis for GRB peaks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000281,"raw_usage":{"total_tokens":1537,"prompt_tokens":814,"completion_tokens":723,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":558,"completion_tokens_details":{"reasoning_tokens":629}},"tokens_in":558,"tokens_out":723,"duration_ms":5775,"temperature":1.0,"reasoning_tokens":629,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T18:19:38.110204+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check: rerun the Bayesian inference with the original, unaltered X-ray uncertainties for XRF 080330 and decay-phase uncertainties for GRB 080710, with the data weights fixed before any fitting. If the fit again drives the isotropic energy to about 10^55 erg (the solution the authors rejected as energetically implausible) or fails to reproduce the rising phase, then the quoted shell width of 10^13 cm and the density slopes are artifacts of the post-hoc weighting. An independent observable test: at about 10^5 s the model predicts a 1.4 GHz radio afterglow near 0.1 mJy for XRF 080330-lik","supporting_citations":[],"review_version":1}