{"id":"1594d09a-af2f-4f63-b87b-fb0a5d0699f0","arxiv_id":"2411.11234","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":15,"one_line_summary":"Prompt GRB gamma-ray spectra can be reproduced by synchrotron and SSC emission from shear-accelerated electrons in a mixed jet-cocoon plus internal-shock electrons in the jet core.","lead":"A model combining shear-accelerated electrons in a gamma-ray burst's jet-cocoon boundary with shock-accelerated electrons in the jet core reproduces the Band-shaped, bimodal, and Band-cut gamma-ray spectra seen in three bright bursts. The work offers a physical explanation for the long-standing Band function and ties prompt emission to structured jets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sec. 2.1 Eq. (1) anchors the exponential velocity profile at r=0 rather than the stated inner boundary r0, biasing the shear-acceleration solution and the fitted B_cn and beta_cn,0.","rationale":"The reader's verdict is CONDITIONAL, and I concur. The most load-bearing issue is not the general flexibility of the two-zone model (an evidence-quality concern) but a concrete internal inconsistency in the defining equation of the model: the exponential velocity profile does not satisfy its stated boundary condition at r0. Because the derived magnetic field B_cn and inner-edge velocity beta_cn,0 are headline results and depend sensitively on the shear profile, this inconsistency must be resolved before the quantitative claims can be accepted. The proposed test—recompute with the corrected boundary condition and refit the three bursts—is decisive and feasible with the existing code. The paper's later power-law comparison (Fig. 6) does not address this issue, since both profiles appear to be anchored at r=0 rather than r0. I therefore do not request a different verdict; I ask that the authors fix this before acceptance. The reader's weakest_assumption identified the velocity profile and particle transport as fragile, and noted a small boundary inconsistency in the rationale, so my concern is a sharper version of the same point rather than a new objection.","tokens_in":15216,"tokens_out":14551,"duration_ms":140096,"concrete_test":"Recompute the model with the corrected profile u_cn(r)=beta_cn,0 exp[-(r-r0) ln(beta_cn,0/beta_cn,2)/(r2-r0)] for r0<r<r2, using the same transport solution (Eqs. 6-7) and the same fitting procedure, and refit the three GRB spectra. Compare the best-fit B_cn and beta_cn,0 with Table 1; if either shifts by more than ~10-20%, the headline parameters are unreliable. The authors should also report the assumed r0/r2 and r1/r2 values so the profile is fully specified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The model's quantitative predictions rest on the exponential velocity profile in Eq. (1), u_cn(r)=beta_cn,0 exp[-(r/r2) ln(beta_cn,0/beta_cn,2)]. The text explicitly states that beta_cn,0 and beta_cn,2 are the fluid velocities at r0 and r2, respectively (Sec. 2.1). However, this formula gives u_cn(0)=beta_cn,0 and u_cn(r2)=beta_cn,2, not u_cn(r0)=beta_cn,0. For the adopted geometry (theta_jet=0.07 rad, theta_cn=0.7 rad, R=10^15 cm), r0/r2 = theta_jet/(theta_cn/2) = 0.2. With beta_cn,0=0.9 and beta_cn,2=1/sqrt(3)~0.577, Eq. (1) yields u_cn(r0)=0.9 exp(-0.2 ln(0.9/0.577))~0.82, about 9% lower than the nominal inner-edge velocity. The shear-acceleration solution (Eqs. 6-7) and the spectral index mu_infty (Eq. 13) depend on xi0-xi2 = artanh(u(r0))-artanh(u(r2)); using 0.82 instead of 0.9 reduces this rapidity difference from ~0.81 to ~0.72, steepening the accelerated electron distribution and shifting the SSC peak. Since the fits in Sec. 3 report beta_cn,0 as the inner-edge velocity but the actual profile is anchored at r=0 (or the inner boundary is implicitly redefined), the derived B_cn=54-450 G and beta_cn,0=0.83-0.91c are not uniquely defined by the stated model. The paper never lists the values of r0, r1, r2 used in the fits, leaving the ambiguity unresolved. This is an internal inconsistency, not a disagreement with external consensus, and it directly affects the headline quantitative claims.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a two-zone structured jet model for GRB prompt emission: an ultra-relativistic jet core with internal-shock-accelerated electrons, plus a surrounding sub-relativistic mixed jet-cocoon (MJC) layer where shear acceleration produces a power-law electron distribution. The model adds synchrotron and SSC emission from both zones and claims that the resulting SEDs reproduce Band-like, bimodal, and Band-Cut spectra. The manuscript applies the model to Fermi GBM+LAT spectra of GRBs 090926A, 131108A, and 160509A, reporting good fits and derived MJC parameters B_cn = 54-450 G and beta_cn,0 = 0.83-0.91c. The authors also discuss an X-ray excess produced by the SSC component and a gamma-gamma opacity check for GRB 160509A.","tokens_in":2158,"tokens_out":2395,"duration_ms":75858,"significance":"If the quantitative claims hold, the model is a significant alternative to one-zone emission scenarios: it offers a physical mechanism for the Band function through shear acceleration, explains bimodal and Band-Cut spectra without invoking separate ad hoc components, and makes falsifiable predictions such as an SSC-driven X-ray excess, a prompt IR/optical flash from MJC synchrotron, and a sub-GeV to GeV spectral tail. The use of public Fermi data and the comparison with three well-known bursts are strengths, as is the explicit discussion of the gamma-gamma transparency condition. However, the evidence is currently conditional: the transport solution is adopted from the literature rather than re-derived, the fitted SEDs are presented without statistical measures, and many parameters that set the spectral shape are not reported. The significance of the derived B_cn and beta_cn,0 ranges therefore depends on the fixes requested below.","major_comments":[{"comment":"The velocity profile u_cn(r) = beta_cn,0 exp[-(r/r2) ln(beta_cn,0/beta_cn,2)] satisfies u_cn(0)=beta_cn,0 and u_cn(r2)=beta_cn,2, yet the text states that beta_cn,0 and beta_cn,2 are the fluid velocities at r0 and r2, respectively. Since r0 > 0 in the adopted jet-cocoon geometry (r0 is set by the jet half-opening angle and r2 by the cocoon half-opening angle), the inner-edge velocity actually used in the transport solution differs from the nominal beta_cn,0; for the stated angles the difference is of order ten percent. The rapidity difference xi0-xi2 in Eq. (8), which controls the spectral index mu_infinity in Eq. (13) and hence the SSC peak, is therefore not the one implied by the quoted beta_cn,0. The manuscript should either modify Eq. (1) to anchor the profile at r0 (e.g., replacing r by r-r0 and normalizing over r0 to r2) or explicitly set r0=0 and define beta_cn,0 accordingly. The values of r0, r1, and r2 used in the fits are never given, so this ambiguity directly affects the headline claims B_cn=54-450 G and beta_cn,0=0.83-0.91c and prevents the fits from being reproducible.","section":"§2.1, Eq. (1)"},{"comment":"The paper states that the three observed spectra are 'well fit' by the model, but no goodness-of-fit statistics are provided: there are no chi-square or likelihood values, no degrees of freedom, no residual plots, and no parameter uncertainties. Given that the model has a large number of free parameters (at least the 16 parameters listed in the text, several of which strongly affect the peak positions and normalizations), a visual comparison alone is insufficient to support the derived parameter ranges. Please provide a quantitative comparison, including at least reduced chi-square or Cash statistic values, parameter confidence intervals, and ideally a model-selection test (e.g., delta-AIC or delta-BIC) against a simple Band-function fit to the same data.","section":"§3, Table 1 and Figure 4"},{"comment":"The reported fits are not reproducible from the manuscript because Table 1 lists only a subset of the model parameters. Missing are the radii r0, r1, r2, the emission radius R, the opening angles theta_cn and theta_jet, the electron number normalization N0, the turbulence parameters (eta, q, kb, kd, delta_B/B, chi), and the cooling/Compton parameter Y_cn used for each burst. Several of these parameters set the absolute flux and the peak energies of the SSC_cn and Syn_jet components, so their omission leaves the fitted B_cn and beta_cn,0 values underdetermined. Please include a complete parameter table or provide the model spectra in electronic form.","section":"§3, Table 1 and §2.2"},{"comment":"The model's SSC_cn peak and the resulting Band-like shape depend sensitively on the assumed injection of a mono-energetic electron population at gamma_e,inject approximately Gamma_jet in the SBL. This assumption is taken from PIC simulations, but no sensitivity study is given: the injection distribution could have a finite width, and the value of gamma_e,inject is varied burst by burst (6.11e2, 5.02e2, 3.31e2) based on external estimates of Gamma_jet that carry their own uncertainties. Because gamma_e,inject is degenerate with beta_cn,0 and B_cn in setting the spectral shape, the paper should demonstrate that the conclusions are robust to reasonable variations in the injection Lorentz factor and to a non-monoenergetic injection spectrum.","section":"§2.2 and Table 1"},{"comment":"The manuscript adopts the Webb et al. (2018) steady-state, incompressible, strong-scattering solution without verifying that these conditions are satisfied for the MJC region, and it does not quantify the effect of relaxing them. Section 4 notes that a power-law velocity profile leaves the 'primary conclusions' unchanged, but Figure 6 is only a qualitative comparison of the distribution shape; it does not show how the fitted B_cn and beta_cn,0 would shift, nor does it address time dependence or compressibility. Since the transport solution is the physical core of the model, a direct integration of Eq. (6) for the same parameters, or a parameter scan over profile shape and equation-of-state assumptions, is needed before the quoted parameter ranges can be considered robust.","section":"§2.2, Eqs. (6)-(14); §4, Figure 6"}],"minor_comments":[{"comment":"The phrase 'is potentially explained the spectral characteristics' should be 'may explain the spectral characteristics', and 'on-borad' should be 'on-board'.","section":"Abstract and §1"},{"comment":"The text says 'as usually observed with BTASE'; this should be 'BATSE'.","section":"§2.3"},{"comment":"The phrase 'the tip of an ice-burger' should be 'the tip of an iceberg'.","section":"§4"},{"comment":"The text refers to a 'broken power-low function'; this should be 'broken power-law function'.","section":"Eq. (17) and surrounding text"},{"comment":"The legend for GRB 090926A labels the combined curve 'SSC_cn + Syn_cn', whereas the text and the other panels identify the high-energy hump as the Syn_jet component; please clarify which components are actually included in each plotted curve.","section":"Figure 4"},{"comment":"The legends showing only '0.99', '0.9', and '0.8' should explicitly state that these values are beta_cn,0, and the line styles should be identified consistently across panels.","section":"Figures 2 and 3"},{"comment":"The sentence 'The SSCcn-component almost dominates the observed in the keV-MeV-GeV band' is missing the word 'spectrum' and should be rephrased.","section":"§3, GRB 160509A paragraph"},{"comment":"The first sentence of the fourth paragraph uses 'Combing'; this should be 'Combining'.","section":"§4"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important question and the proposed mechanism is physically interesting, but the quantitative claims currently outrun the evidence. In particular, the anchoring inconsistency of the velocity profile in Eq. (1) and the complete absence of fit statistics and parameter uncertainties make the derived B_cn and beta_cn,0 ranges difficult to evaluate. I recommend major revision rather than rejection because the issues are fixable within the manuscript's scope. The authors should also state clearly that the three bursts are a proof-of-concept sample and not a systematic population study. I saw no indication of inappropriate citation behavior; the reliance on Webb et al. (2018) and on co-authored PIC work is disclosed in the text and is a normal scientific dependence, though it should be tested numerically before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper you're asking about is a two-zone model for GRB prompt spectra: shear acceleration in the mixed jet-cocoon (MJC) region plus internal shocks in the jet core, with synchrotron and SSC from both. The new part is the specific spectral synthesis and its application to three Fermi GRBs (090926A, 131108A, 160509A). The underlying transport solution is from Webb et al. 2018, and the structured jet-cocoon picture is not new, so the novelty is moderate but real.\n\nWhat I like: the physical motivation is legitimate, and the idea that the Band function could be a blend of a sub-relativistic cocoon SSC peak and a jet-core synchrotron component is worth taking seriously. The authors are also candid about limitations—they test the exponential vs. power-law velocity profile and show the qualitative electron distribution is robust. The derived B_cn (54-450 G) is in a plausible range.\n\nThe soft spots are more than cosmetic. Most importantly, Eq. (1) is internally inconsistent: the text says β_cn,0 and β_cn,2 are the velocities at r0 and r2, but the formula anchors the exponential at r=0, so u_cn(r0) is about 9% lower than β_cn,0 for their adopted geometry. Since the shear profile sets the acceleration efficiency and the spectral index, this changes the predicted spectra and the fitted B_cn and β_cn,0. This is a load-bearing bug, not a typo; the paper never lists r0, r1, r2, so the ambiguity is unresolved.\n\nSecond, the fitting is not statistically grounded. There are roughly a dozen free parameters per burst, no error bars, no goodness-of-fit, no residuals, and no comparison to a simpler Band-function fit. With that many knobs, \"well fit\" is a weak claim. The bursts are selected for the very shapes the model is built to produce, so the evidence is suggestive, not confirmatory.\n\nThat said, the central claim is not circular: the Band-like shape follows from the shear-acceleration solution, which is independent of the data. The issues are fixable with a clearer geometry and a proper likelihood analysis.\n\nThis paper is for GRB prompt-emission specialists and particle-acceleration theorists. It deserves peer review, but the referee should demand a corrected Eq. (1), a full list of the geometric parameters, and a quantitative fit statistic. My own verdict is conditional: I wouldn't quote the derived B_cn or β_cn,0 until those points are addressed.\n\nRecommendation: send it to review, with the expectation of major revision.","headline":"A plausible but over-flexible two-zone shear-acceleration model for GRB prompt spectra; the quantitative fits are undercut by an inconsistency in the velocity profile and a lack of statistical rigor.","tokens_in":16397,"tokens_out":6725,"would_cite":true,"duration_ms":56277,"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 prompt gamma-ray burst spectra, including the canonical Band function, bimodal spectra, and Band-Cut spectra, can be produced by a structured jet in which shear-accelerated electrons in a sub-relativistic mixed…","keywords":["gamma-ray bursts","structured jets","shear acceleration","synchrotron self-Compton","Band function","jet-cocoon structure","GRB 090926A","non-thermal radiation"],"falsifier":"A single detection of a prompt photon above about 6 GeV from GRB 160509A would violate the paper's transparency calculation for the cocoon region, since the model requires the gamma-gamma optical depth to drop below unity only below that energy; more broadly, a stacked search of GRB spectra that finds no keV X-ray excess in bursts fitted as Band-Cut would undercut the cocoon SSC component.","tokens_in":14952,"feed_emoji":"💥","tokens_out":6094,"duration_ms":54293,"temperature":0.7,"pith_summary":"The paper argues that the observed gamma-ray spectra of gamma-ray bursts can be explained by a structured jet with two emission zones: an ultra-relativistic core where internal shocks accelerate electrons, and a sub-relativistic mixed jet-cocoon layer where shear acceleration energizes electrons. In the cocoon layer, synchrotron self-Compton cooling dominates and its keV peak anchors the low-energy part of the Band function, while synchrotron emission from the jet core adds the higher-energy component. The authors fit time-integrated Fermi GBM and LAT spectra of GRBs 090926A, 131108A, and 160509A and derive cocoon magnetic fields from 54 to 450 G and inner-edge velocities from 0.83c to 0.91c. If correct, the prompt spectrum is not evidence for a single radiation mechanism but for a two-zone structured outflow.","feed_headline":"Shear acceleration can shape GRB Band spectra","feed_subtitle":"Two-zone jet model fits three bright bursts and yields cocoon fields of 54–450 G.","key_machinery":"The engine of the argument is the steady-state cosmic-ray transport equation for shear acceleration in the strong-scattering limit, solved analytically for an incompressible flow with an exponential velocity profile $u_{\\rm cn}(r)=\\beta_{\\rm cn,0} e^{-kr}$. The solution provides a shear-accelerated electron distribution $f_0$ that rises as $p^{-\\mu_0}$ below the injection momentum and falls as $p^{-\\mu_\\infty}$ above it, with the indices set by the velocity contrast across the layer and by Kolmogorov turbulence. Equating the shear acceleration time with the synchrotron plus synchrotron self-Compton cooling time sets the maximum electron energy, and the SSC cooling then dominates. This cocoon electron distribution, combined with a broken power-law electron distribution for internal-shock electrons in the jet core, is the mechanism that produces the Band-like and Band-Cut spectral energy distributions.","core_discovery":"The central claim is that the Band function and related prompt GRB spectra arise from the superposition of two physically distinct electron populations. Shear acceleration in the mixed jet-cocoon region, with a steadily decaying velocity profile and strong scattering, produces electrons up to Lorentz factor roughly $10^4$; their synchrotron and synchrotron self-Compton emission peaks at keV energies and supplies the low-energy spectral component. Internal-shock-accelerated electrons in the jet core, with Lorentz factors $10^4$ to $10^5$, emit synchrotron radiation peaking around the keV to MeV band, supplying the high-energy component. Adding the two reproduces the observed bimodal and Band-Cut shapes. For the three bursts studied, the model gives cocoon magnetic fields of 54 to 450 G and inner-edge velocities of 0.83c to 0.91c, with the jet core highly relativistic and strongly magnetized.","pith_inferences":["A testable extension beyond the paper is that time-resolved spectra should show the keV SSC peak and the MeV-to-GeV synchrotron hump varying with different temporal lags if they originate in two distinct zones.","The same two-zone picture could be extended to explain the low-energy X-ray excess reported in some BATSE bursts; a stacked search over many GRBs for the predicted excess at a few keV would provide a statistical test.","The paper fixes the cocoon velocity profile to an exponential decay and only checks a power-law profile in a comparison figure; the quoted $B_{\\rm cn}$ and $\\beta_{\\rm cn,0}$ ranges should be read as profile-dependent, and constraining the actual breakout dynamics would sharpen the fits."],"forward_implications":["If the model is right, the canonical Band function is a composite, so fitting GRB spectra with a single Band function may be averaging two physically separate emission components.","The keV X-ray excess seen in some bursts can be identified with the peak of the cocoon SSC component, giving a concrete physical origin to that excess.","The high-energy MeV-to-GeV hump in bimodal GRBs is attributed to the jet core's synchrotron emission, so its presence traces the core parameters rather than the cocoon physics.","Derived cocoon magnetic fields of 54 to 450 G and inner-edge velocities of 0.83c to 0.91c become observational constraints on jet-cocoon structure that simulations and afterglow modeling should reproduce.","The model provides a natural way to produce Band-Cut spectra, in which the high-energy cutoff is not intrinsic to the emission process but reflects the relative strengths of the cocoon SSC and core synchrotron components."],"supporting_citations":[{"why":"Supplies the analytical steady-state solution for the shear-accelerated particle distribution that the model uses as its core cocoon component.","marker":"Webb et al. 2018"},{"why":"Establishes the strong-scattering limit and the shear acceleration transport regime adopted for the mixed jet-cocoon region.","marker":"Rieger & Duffy 2005"},{"why":"Defines the empirical Band function that the model aims to reproduce from physical emission components.","marker":"Band et al. 1993"},{"why":"Provides the bimodal spectrum of GRB 090926A used as one of the three case-study targets.","marker":"Ackermann et al. 2011"},{"why":"Reports Band-Cut spectra that the model claims to explain with the two-zone emission structure.","marker":"Ackermann et al. 2013"},{"why":"Supplies the scattering-time and mean-free-path relations linking turbulence properties to particle transport in the shear layer.","marker":"Liu et al. 2017"},{"why":"Particle-in-cell simulation evidence for particle acceleration at the shear boundary layer, which motivates the electron injection layer.","marker":"Alves et al. 2014"},{"why":"Source for the jet-cocoon structure and for the constraint that the outer cocoon velocity stays below the local sound speed.","marker":"Ramirez-Ruiz et al. 2002"},{"why":"Reports the low-energy X-ray excess in BATSE bursts, which the model reproduces as the cocoon SSC peak.","marker":"Preece et al. 2000"}],"fun_headline_variants":["Shear acceleration in GRB jets yields Band spectra","Two-zone jet model reproduces GRB Band function","Cocoon shear plus core shocks explain GRB spectra","Shear-driven electrons set GRB Band shape","Mixed jet-cocoon shear shapes GRB spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on the assumption that the mixed jet-cocoon layer is a steady, incompressible shear flow with an exponential velocity profile, strong scattering, and a mono-energetic electron injection at a Lorentz factor roughly equal to the jet Lorentz factor; if the real layer is time-dependent, compressible, differently shaped, or injects electrons differently, the predicted spectra and the derived field and velocity values would change.","fun_headline_variants_meta":{"raw":{"variants":["Shear acceleration in GRB jets yields Band spectra","Two-zone jet model reproduces GRB Band function","Cocoon shear plus core shocks explain GRB spectra","Shear-driven electrons set GRB Band shape","Mixed jet-cocoon shear shapes GRB spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000461,"raw_usage":{"total_tokens":2408,"prompt_tokens":1150,"completion_tokens":1258,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":766,"completion_tokens_details":{"reasoning_tokens":1181}},"tokens_in":766,"tokens_out":1258,"duration_ms":9544,"temperature":1.0,"reasoning_tokens":1181,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:46:47.009275+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single detection of a prompt photon above about 6 GeV from GRB 160509A would violate the paper's transparency calculation for the cocoon region, since the model requires the gamma-gamma optical depth to drop below unity only below that energy; more broadly, a stacked search of GRB spectra that finds no keV X-ray excess in bursts fitted as Band-Cut would undercut the cocoon SSC component.","supporting_citations":[],"review_version":1}