{"id":"64f8c807-b92a-47a7-a06e-87d299f2862c","arxiv_id":"2504.15011","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Off-axis observers of structured gamma-ray-burst jets should see quasi-thermal photosphere spectra that are fainter, softer, and slower to evolve than on-axis spectra, and the Einstein Probe and SVOM should detect such emission from 170817A-like short bursts within 200 Mpc.","lead":"This paper computes the spectrum of the hot photosphere of a structured gamma-ray-burst jet seen from different viewing angles, including the off-axis geometry of GRB 170817A. It predicts that the Einstein Probe and SVOM telescopes can detect such off-axis thermal emission from short bursts at up to 200 Mpc, providing a new probe of jet structure.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The finite-boundary early-time spectra, including the claimed hard component, are computed at radii below the saturation radius where Eq. (6) does not apply; the paper's own §4 caveat flags this but its quantitative impact on Figures 3-5 is untested.","rationale":"The reader identified the saturation-regime assumption as the weakest premise, and I agree. My stress-test focuses that concern on the regime where the paper's finite-boundary novelty lives. The infinite-boundary spectra and the later-time finite-boundary spectra are less exposed, because by the time the photosphere approaches the saturation radius the saturated formula is adequate. But the early-time curves in Figures 3-5, especially the hard component, rely on radii r < rs, where Eq. (6) is not valid. This is not an internal algebraic inconsistency; it is a domain-of-application issue that can be addressed by including the acceleration phase. Even though T_obs = T0/2 may be the correct pre-saturation observed temperature, the Lorentz factor entering the density and Doppler factor is not the coasting value, so the early optical depth and emitted flux are misestimated. The paper's own Section 4 caveat is honest but does not quantify the effect. If unsaturated dynamics suppress the hard component, the paper's claim about the finite-boundary treatment and its observational implications weakens; if the hard component survives, the central claim is strengthened. I do not elevate the photon-number-conservation normalization in Eq. (18) to the primary concern, because a recalibration of A would shift absolute flux levels but would not address the saturation-domain problem; nevertheless, it remains a secondary validation task. Since the detectability claim for EP-WXT and SVOM-ECLAIRs depends on absolute flux as well as spectral shape, the concrete test above should be run before the central claims are accepted. The appropriate verdict is therefore CONDITIONAL, matching the reader's assessment.","tokens_in":19456,"tokens_out":18397,"duration_ms":188701,"concrete_test":"Recompute the finite-boundary spectra of Figure 5 at tobs = 10^-5, 10^-3, and 10^-2 s using unsaturated dynamics: set Γ(r,θ) = min(r/r0, Γ(θj)), use the corresponding density n'(r,θ) = L/(4πmpβc^3Γ^2 r^2), Doppler factor D = [Γ(1−β cosθ)]^-1, and the pre-saturation temperature T_obs = T0/2 constant for r < rs, then the r^-2/3 law for r > rs. Compare the resulting early flux and the 10-100 keV spectral shape with Figures 5(a)-(d). If the early flux changes by more than a factor of 2 or the hard-component spectral index changes by more than 0.3, the finite-boundary claim requires revision; if the curves are essentially unchanged, the saturation caveat is not fatal.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central novelty of the paper is the finite-boundary result of Section 3.2: early radiation escapes from smaller radii and produces an enhanced flux plus a hard high-frequency component that persists off-axis. This result inherits the temperature prescription of Eq. (6), which is derived under the saturated-regime assumption rs(θj) < rph(θj) and η(θj) = Γ(θj). For the fiducial parameters used in Figures 3-5, however, the earliest plotted observer times correspond to engine times at which the finite outer boundary rout = β(θj)cˆt lies below the saturation radius: taking r0 = 10^7 cm (the scale used in Eq. 30) and Γc = 300, rs ≈ 3×10^9 cm, while rout ≈ 3×10^5 cm at ˆt = 10^-5 s and ≈ 3×10^7 cm at ˆt = 10^-3 s. Photons escaping at these times originate from the acceleration phase, where Γ(r) ≈ r/r0 < Γ(θj), so the density used in Eq. (16), the Doppler factor in Eq. (5), and the optical depth in Eq. (26) are all evaluated with the wrong local Lorentz factor. Applying the saturated temperature law in this regime can overestimate the early flux and distort the spectrum. The paper acknowledges in Section 4 that an unsaturated situation may hold for parts of the jet, but does not identify that the finite-boundary early-time curves in Figures 3-5 sit squarely in that regime. Since the hard component is one of the two headline spectral signatures, the claim is not yet supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents numerical calculations of instantaneous photosphere emission spectra from a structured gamma-ray burst jet, with an angularly dependent luminosity and Lorentz factor, for observers at different viewing angles. The authors compute spectra for constant and variable central-engine luminosity, compare an infinite outer boundary with a finite outflow boundary, and estimate the detectability of quasi-thermal photosphere emission by EP-WXT, SVOM-ECLAIRs, and Swift-BAT. The central claims are that off-axis spectra are fainter, peak at lower energies, and evolve more slowly than on-axis spectra; that treating the outflow boundary as finite enhances the early flux and produces a hard high-frequency component; and that short GRBs similar to GRB 170817A should be detectable to about 200 Mpc for viewing angles below 10 degrees.","tokens_in":19894,"tokens_out":8549,"duration_ms":81564,"significance":"If the finite-boundary result holds, the paper offers a new observational handle on GRB jet structure through quasi-thermal photosphere spectra, linking prompt-emission modeling to the structured-jet picture favored by afterglow observations. The forward-modeling approach is a strength: the spectra are derived self-consistently from the assumed jet structure, and the short-GRB jet parameters are taken from independent afterglow fits (Li et al. 2019) rather than fitted to the predicted outputs. The predicted differences between on-axis and off-axis temporal and spectral evolution are falsifiable with current and upcoming wide-field X-ray instruments. However, the quantitative early-time results and the detectability estimates rest on two assumptions that are not adequately validated: the application of the saturated-regime temperature law to radii below the saturation radius, and the photon-number-conservation normalization of the last-scattering probability.","major_comments":[{"comment":"The finite-boundary early-time spectra are computed at emission radii far below the saturation radius. For the fiducial parameters used in the figures (r0=10^7 cm, Γc=300, Eq. 30), rs=Γr0≈3×10^9 cm, whereas rout=βct is about 3×10^5 cm at t=10^-5 s and 3×10^7 cm at t=10^-3 s. In this regime the local Lorentz factor is still Γ(r)≈r/r0, not Γ(θj), so the density in Eq. (16), the Doppler factor in Eq. (5), and the optical depth in Eq. (26) are evaluated with an incorrect Lorentz factor. Since the hard high-frequency component and the early flux enhancement are the paper's headline new signatures, their quantitative support is missing unless the acceleration phase is modeled. The Section 4 caveat about an unsaturated situation acknowledges the general issue but does not identify that the earliest finite-boundary curves in Figs. 3-5 lie in the sub-saturation-radius regime; the authors should either implement a proper acceleration-phase treatment or restrict the finite-boundary claims to times when rout exceeds rs.","section":"Section 3.2, Figs. 5-6, Eqs. (6), (8), (24)"},{"comment":"The normalization of the last-scattering probability P(r,Ω) is replaced by global photon-number conservation over all observer directions. However, Eq. (11) uses P as a probability density weighting the contribution of individual injected photons to a specific observer, for which the natural normalization is ∫∫P dr dΩ=1 per injection direction. Because the Doppler factor D and the optical depth τ depend on the observer direction θv, the constant A determined by Eq. (19) is an observer-averaged quantity. The manuscript does not show that this global normalization is equivalent to the per-observer probability normalization, and the issue affects the absolute flux level in every figure, including the detectability estimates. A derivation or a numerical validation against the spherically symmetric limit (e.g., Pe'er 2008) is needed.","section":"Section 3.1, Eqs. (18)-(19)"},{"comment":"The detectability conclusion that EP-WXT and SVOM-ECLAIRs can detect 170817A-like bursts within a viewing angle of 10 degrees out to 200 Mpc depends on the correctness of the early-time finite-boundary spectra, which are affected by the sub-saturation-radius problem noted above. In addition, the flux-angle dependence in Fig. 9 is presented only for a few parameter variations, with fixed r0 and fixed luminosity-history indices; given that r0 enters the peak-energy scaling in Eq. (30) and that the short-GRB case uses θv=27.6 degrees with Lc=10^51 erg/s, the 200-Mpc reach should be framed as conditional on the acceleration-phase treatment and on the parameter choices, not as a robust prediction.","section":"Section 3.4, Fig. 9, Eq. (30)"}],"minor_comments":[{"comment":"The text refers to a \"Plank distribution\"; this should be \"Planck distribution\".","section":"Section 2.2, Eq. (13)"},{"comment":"Equation (22) cites Eq. (26) for the optical depth, but Eq. (26) is introduced later in Section 3.2; please reorder the equations or adjust the cross-reference.","section":"Section 3.1, Eq. (22)"},{"comment":"The luminosity history is written as 10^{ar log t+br}; this is equivalent to a broken power law, and writing Lc(t)=Lcp(t/tp)^ar and Lc(t)=Lcp(t/tp)^ad for the two branches would be clearer.","section":"Section 3.3, Eq. (28)"},{"comment":"The numerical method is described only in words; a reproducibility statement with the grid resolution, integration tolerances, and convergence checks would strengthen the paper, especially because the finite-boundary case requires solving for r2 in Eq. (25) numerically.","section":"Section 3.2 and Appendix C"},{"comment":"The peak-energy scaling in Eq. (30) is stated without intermediate steps; a short derivation in an appendix would help readers verify the exponents 5/12 ke - 8/3 kΓ.","section":"Eq. (30)"}],"recommendation":"major_revision","confidential_remarks":"The topic is timely and the forward-modeling approach is a step forward, but the two load-bearing issues—sub-saturation-radius emission in the finite-boundary early-time spectra and the photon-number-conservation normalization of P(r,Ω)—need to be addressed before publication. If the authors can quantify the acceleration-phase correction or show that it is negligible for the plotted times, and can justify or fix the normalization, the paper could become a solid contribution to the structured-jet photosphere literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. The paper gives the first off-axis instantaneous photosphere spectra for structured jets, and a concrete EP-WXT/SVOM-ECLAIRs detection forecast. The framework is useful. But the headline new result—the hard high-frequency component from the finite-boundary treatment—is computed at engine times when the outflow is still in the acceleration phase, below the saturation radius, where the temperature law used in the calculation does not apply.\n\nWhat is actually new and good: the off-axis viewing-angle dependence of instantaneous photosphere spectra for structured jets; the finite-boundary catch-up geometry; the spectral evolution from Fν ∝ ν^2 to Fν ∝ ν^1.5; and a quantitative detectability estimate for 170817A-like short GRBs. The calculation is internally consistent, no parameter is fitted to the predicted outputs, and the short-GRB jet parameters are taken from an earlier afterglow fit, so the circularity burden is low. The citation practice looks normal.\n\nSoft spots, in order of severity.\n\nFirst, and this is load-bearing: the stress-test note is right. With Γc=300 and r0=10^7 cm, rs≈3×10^9 cm. At tobs=10^-5–10^-3 s, the finite outer boundary is at 3×10^5–3×10^7 cm. Those early spectra in Figures 3–5 come from radii where Γ(r)≈r/r0 < Γ(θj), so the density, Doppler factor, and optical depth are all evaluated with the wrong local Lorentz factor. The hard component is consequently not supported. The paper's own §4 caveat about unsaturated jets does not connect the dots to its own figures. The fix is to either restrict the finite-boundary analysis to times when rout > rs, or to implement the acceleration-phase temperature/density profile.\n\nSecond: the photon-number-conservation normalization in Eqs. 18–19 is nonstandard. It might be right, but it needs explicit validation—a Monte Carlo scattering test or a derivation showing it is equivalent to the usual probability normalization. The absolute flux, and therefore the detectability claim, rests on it.\n\nThird, a minor one: the detector forecast compares band-integrated flux to sensitivity limits without a background or signal-to-noise treatment. Fine as an order-of-magnitude statement, not as a detection-rate prediction.\n\nWho is this for: GRB prompt-emission theorists and mission planners. It deserves a serious referee. I would send it to review, with a request to fix the acceleration-phase problem and validate the normalization. I would cite the off-axis structured-jet framework, but not the hard component, in its current form.","headline":"Useful off-axis structured-jet photosphere spectra, but the headline hard component rests on applying the saturated temperature law in the acceleration phase; conditionally accept after major revision.","tokens_in":20401,"tokens_out":4659,"would_cite":true,"duration_ms":43837,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Structured GRB jets emit photosphere spectra that depend strongly on viewing angle: off-axis emission is fainter, softer, and slower to evolve, and the usual infinite-boundary treatment hides an early hard high-frequency component.","keywords":["gamma-ray bursts","photosphere emission","structured jets","off-axis viewing angle","quasi-thermal spectra","GRB 170817A","finite outflow boundary","SVOM-ECLAIRs"],"falsifier":"Detect the early X-ray spectrum of a 170817A-like short burst at about 200 Mpc viewed 10–30 degrees off-axis: the finite-boundary, saturated-jet model predicts a quasi-thermal spectrum with a $\\nu^2 \\to \\nu^{1.5}$ low-energy slope, a peak migrating into the 0.5–4 keV band, and a hard tail that persists past 100 s at $\\theta_v = 20^\\circ$; observing instead a pure blackbody cutoff, a fast on-axis-like evolution, or no thermal component at all would rule out the model.","tokens_in":19266,"feed_emoji":"🔭","tokens_out":13402,"duration_ms":106260,"temperature":0.7,"pith_summary":"This paper asks what the quasi-thermal photosphere emission of a gamma-ray burst looks like when the jet is structured — a fast inner core with power-law wings in luminosity and Lorentz factor — and the observer sits off-axis. It claims that the instantaneous spectrum depends strongly on viewing angle $\\theta_v$: off-axis emission is fainter, peaks at lower energy, and evolves far more slowly than on-axis emission, because the photosphere radius and temperature vary with latitude. It further argues that the common infinite-boundary treatment of the outflow is inadequate at early times; with a finite boundary, early flux is several orders of magnitude higher, and a hard high-frequency tail appears that persists longest for large viewing angles. For 170817A-like short bursts, the model predicts that EP-WXT and SVOM-ECLAIRs should detect this thermal emission out to about 200 Mpc for $\\theta_v \\lesssim 10^\\circ$, giving observers a new handle on jet structure.","feed_headline":"Two X-ray missions can catch off-axis GRB glow to 200 Mpc","feed_subtitle":"Structured jets viewed from the side look fainter and softer, and the thermal signal reveals the burst's geometry.","key_machinery":"The load-bearing object is the last-scattering photosphere of the structured jet: the shell where the optical depth to the observer reaches unity, whose radius and temperature vary with polar angle through the power-law profiles $L(\\theta_j)$ and $\\Gamma(\\theta_j)$. The argument runs through the comoving temperature profile $T'(r,\\theta_j)$ with its saturation radius $r_s = \\Gamma(\\theta_j) r_0$ and photosphere radius $r_{\\rm ph}$ (set by $\\tau = 1$), the Doppler-boosted probability density $P(r,\\theta,\\phi)$ for the last scattering, and the time-delay geometry $t_{\\rm obs}/(1+z) = \\hat{t} + r u/(\\beta c)$ that maps emission latitude onto observed time. The decisive element is the treatment of the outflow's outer boundary: an infinite boundary integrates the optical depth to infinity, whereas the finite-boundary treatment lets photons catch up with the expanding edge at $r_{\\rm out} = \\beta(\\theta_j) c \\hat{t}$, which brightens the early flux and produces the hard high-frequency component. The viewing-angle dependence enters through the Doppler factor $D = [\\Gamma(\\theta_j)(1 - \\beta(\\theta_j)\\cos\\theta)]^{-1}$ and the geometric relation $\\theta_j(\\theta,\\phi,\\theta_v)$.","core_discovery":"The central claim is that the observed photosphere spectrum of a structured GRB jet is a strong function of the viewing angle. Because luminosity and Lorentz factor fall off as power laws away from the jet core, an off-axis line of sight samples a larger photosphere radius at lower temperature: the flux density and peak energy drop, and the spectrum takes much longer to settle into its quasi-saturated shape — about $10^{-4}$ s on-axis versus $\\sim 10^2$ s at $\\theta_v = 20^\\circ$. The paper also claims that the standard infinite-boundary approximation is wrong at early times: when the outflow's finite outer boundary is treated properly, early photons escape before accumulating an artificial optical depth, so the flux is orders of magnitude higher and a power-law-like hard component sits above the thermal peak, disappearing gradually as the photosphere radius converges to the infinite-boundary value. As a corollary, the spectral peak energy tracks inversely with the luminosity history (hard-to-soft in the rise, soft-to-hard in the decay), with all evolution delayed for off-axis observers. The paper closes by predicting that EP-WXT and SVOM-ECLAIRs can detect quasi-thermal photosphere emission from 170817A-like short bursts out to roughly 200 Mpc when the viewing angle is below about $10^\\circ$.","pith_inferences":["The predicted drift of the low-energy spectral index from $\\nu^2$ toward $\\nu^{1.5}$, together with the persistent hard tail, offers a discriminating test against synchrotron emission in joint fits of the same off-axis events, since the two mechanisms predict different index trajectories over time.","If the saturation assumption fails at some latitudes ($r_s > r_{\\rm ph}$), the clean scaling of peak energy with viewing angle breaks down; measuring the peak energy versus viewing angle across a sample of off-axis short bursts would directly probe which latitudes of the jet are actually saturated.","Wide, soft X-ray surveys such as EP-WXT may be systematically biased toward off-axis thermal events, because the redshifted peak energy of off-axis emission lands squarely in their 0.5–4 keV band; this selection effect could be checked by comparing the inferred viewing-angle distribution of soft-band and hard-band detected GRBs.","Extending the calculation to include sub-photospheric Comptonization, which the paper itself lists as needed future work, would show whether the hard high-frequency component survives spectral processing by scattering before the photons escape."],"forward_implications":["Off-axis photosphere spectra are fainter, peak at lower energy, and evolve more slowly than on-axis spectra, so the saturation timescale itself — from about $10^{-4}$ s on-axis to about $10^2$ s at $20^\\circ$ — becomes a viewing-angle diagnostic.","The finite-boundary treatment raises the early flux by orders of magnitude and reveals a hard high-frequency component that persists for tens of seconds at large viewing angles; the infinite-boundary approximation cannot reproduce these features.","With a variable central engine, the peak flux tracks the luminosity history while the peak energy anti-correlates with it (hard-to-soft during the rise, soft-to-hard during the decay), with off-axis evolution delayed relative to $(1+z) t_p$.","EP-WXT and SVOM-ECLAIRs should detect quasi-thermal emission from 170817A-like short bursts out to about 200 Mpc as long as the viewing angle stays below roughly $10^\\circ$.","The detection flux falls with viewing angle more steeply than the luminosity profile, because the optical depth also grows off-axis; for EP-WXT, the most favorable viewing angle is structure-dependent rather than on-axis, since the off-axis peak energy moves into its 0.5–4 keV band."],"supporting_citations":[{"why":"Supplies the last-scattering probability density function that the paper generalizes to azimuth-dependent structured jets in Equation (14).","marker":"Pe'er 2008"},{"why":"Provides the instantaneous specific-flux integral (Equation 11) and the uniform-jet spectral evolution that the structured-jet calculation extends.","marker":"Deng & Zhang 2014"},{"why":"Source of the comoving temperature profile in Equation (6), the thermal structure underpinning all the computed spectra.","marker":"Mészáros & Rees 2000"},{"why":"Basis of the infinite-boundary optical depth integral (Equation 15) that the finite-boundary treatment is designed to replace.","marker":"Abramowicz et al. 1991"},{"why":"Defines the saturation radius and photosphere radius that anchor the temperature profile and the optical-depth-unity condition.","marker":"Daigne & Mochkovitch 2002"},{"why":"Establishes the saturated photosphere regime ($r_s < r_{\\rm ph}$) and the photosphere temperature entering the peak-energy estimates.","marker":"Pe'er et al. 2007"},{"why":"Evidence that photosphere emission can interpret the spectral feature of GRB 170817A, motivating the short-GRB detectability calculation.","marker":"Meng et al. 2018"},{"why":"Provides the jet structure and viewing angle of GRB 170817A used in Section 3.4 to generate the predicted off-axis spectra.","marker":"Li et al. 2019"}],"fun_headline_variants":["Off-axis GRB jets show cooler, softer thermal glow","Structured jet viewing angle shapes photosphere spectrum","Finite jet edge boosts early light from GRB photosphere","EP and SVOM can spot off-axis short GRBs out to 200 Mpc","Viewing angle sets thermal peak and luminosity in GRBs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the jet is already fully accelerated at every latitude before its light escapes (the photosphere lies beyond the saturation radius), so the simple temperature profile of Equation (6) holds everywhere — a condition the paper itself says may fail for dim parts of the jet.","fun_headline_variants_meta":{"raw":{"variants":["Off-axis GRB jets show cooler, softer thermal glow","Structured jet viewing angle shapes photosphere spectrum","Finite jet edge boosts early light from GRB photosphere","EP and SVOM can spot off-axis short GRBs out to 200 Mpc","Viewing angle sets thermal peak and luminosity in GRBs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000533,"raw_usage":{"total_tokens":2602,"prompt_tokens":1024,"completion_tokens":1578,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":640,"completion_tokens_details":{"reasoning_tokens":1492}},"tokens_in":640,"tokens_out":1578,"duration_ms":10122,"temperature":1.0,"reasoning_tokens":1492,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:36:50.252322+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Detect the early X-ray spectrum of a 170817A-like short burst at about 200 Mpc viewed 10–30 degrees off-axis: the finite-boundary, saturated-jet model predicts a quasi-thermal spectrum with a $\\nu^2 \\to \\nu^{1.5}$ low-energy slope, a peak migrating into the 0.5–4 keV band, and a hard tail that persists past 100 s at $\\theta_v = 20^\\circ$; observing instead a pure blackbody cutoff, a fast on-axis-like evolution, or no thermal component at all would rule out the model.","supporting_citations":[],"review_version":1}