{"id":"80bcb7e2-a2a7-4d8c-b55f-ff9c83911a5e","arxiv_id":"2607.24234","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Without the slowly-varying-envelope approximation, spatiotemporal current-filamentation instability dominates near an ultrarelativistic beam front and is overtaken by oblique two-stream instability only farther back.","lead":"A unified quasistatic model shows that magnetic filamentation, not the usual oblique two-stream mode, dominates right behind a relativistic beam front, with the switchover fixed by a simple ratio of distance to growth time. The result matters for plasma accelerators, laser-plasma experiments, and astrophysical jets where beams always have a leading edge.","discovery_kind":"unification","skeptic_critique":null,"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript develops a unified, fully electromagnetic, cold-fluid theory of the linear stage of streaming instabilities driven by bounded ultrarelativistic beams in unmagnetized plasma. Using the quasistatic approximation (QSA) in comoving coordinates and deliberately avoiding the slowly varying envelope approximation (SVEA), the authors derive a single master equation [Eq. (15)] for the spatiotemporal evolution of perturbations. For transversely wide beams they reduce it to Eq. (18), solve it by double Laplace transform and steepest descent (Appendices A–C, including an exact Bessel-function solution in the ∂²_ξ ≫ 1 limit), and find that spatiotemporal CFI — with a scaling distinct from the earlier SVEA result of Ref. [56] — dominates a front region ξ ≲ ξ_tr ≃ 0.385 Γ_CFI τ, while spatiotemporal OTSI takes over deeper in the beam; the transition corresponds to a coalescence of saddle points. The model is then extended to Gaussian longitudinal profiles [Eq. (32)], to longitudinal TSI for wide beams, and to SMI/HI for narrow beams, which are shown to share the OTSI spatiotemporal dynamics in the k_y → ∞ limit. 1D and 2D PIC simulations (CALDER) are used throughout to validate the asymptotic predictions. The central claims are (i) the previously unreported front-region spatiotemporal CFI under conditions where temporal theory predicts OTSI dominance, (ii) the parameter-free transition boundary, and (iii) the unification of SMI/HI with OTSI. The derivation is standard and透明ly","tokens_in":35405,"tokens_out":8544,"duration_ms":278720,"significance":"If it holds, this is a useful and genuinely new contribution to the theory of relativistic beam-plasma instabilities. The strengths are substantial: (i) a single parameter-free master equation spanning CFI, OTSI, TSI, SMI and HI; (ii) a quantitative, falsifiable prediction — the CFI/OTSI transition at ξ_tr = (2/√27)Γ_CFI τ and the ~0.4 N_e skin-depth extent of front filaments, independent of beam density and Lorentz factor; (iii) an exact Bessel-function Green's function in the sharp-front limit and fully worked steepest-descent asymptotics (including the second-order coalescing saddle at C = C_cr), cross-checked against direct numerical inversion (Figs. 10–12); (iv) independent PIC validation for flat-top and Gaussian beams with stated parameters, plus a consistent SMI/HI/OTSI unification in the k_y → ∞ limit. The results bear directly on AWAKE/FACET-II-class experiments, on pair-beam propagation in blazar/IGM debates, and they motivate quasistatic PIC as the practical tool for dilute ultrarelativistic beams. The departure from the earlier SVEA-based CFI picture [56] is derived, not asserted, and the conditions of validity (QSA, ξ ≲ τ/3, cold fluid, γ_b ≫ 1) are stated plainly.","major_comments":[{"comment":"The comparison in Fig. 2(c) of fixed-position simulation data (a window spanning 0 ≤ kpξ ≲ 4) against the critical-point solution Eq. (30) is not strictly justified as written. Eq. (30) describes an observer moving along the ray ξ = C_cr Γ_CFI τ; at fixed ξ, the asymptotics are Eq. (28) (with a time-dependent local rate d/dτ[2√(Γξτ)] = √(Γξ/τ)) or Eq. (29), depending on the sign of ξ − ξ_tr(τ). Since ξ_tr sweeps through the diagnostic window during the run, the observed constant slope ≈ 1.09 Γ_CFI plausibly reflects the window straddling the moving transition layer, but this should be demonstrated — e.g., by overplotting the local slope predicted by Eq. (28)/(29) evaluated at the windowed ξ, or by extracting the amplitude along the moving ray ξ = ξ_tr(τ) from the simulation. This point is load-bearing because Fig. 2(c) is the evidence for the transition region, one of the paper's central","section":"§III.C, Fig. 2(c)"},{"comment":"The analytical results are Green's functions for a front-localized impulse (Eqs. 24–30) or for an extended disturbance imposed at τ = 0 (Appendix C), whereas the PIC modes grow from particle-discreteness noise that is generated continuously and throughout the beam volume. Appendix C partially addresses this by showing robustness to an extended initial seed (including the neat cancellation of the e^{Γτ} terms between K and the J±_3 loop integrals), but a continuously generated volume source is a different convolution. Given that the theory–simulation agreement in Figs. 2–4 is the principal validation of the central claim, a short argument in §III.C that convolution of the Green's function with broadband volume noise is dominated by the same saddle contributions (or a direct numerical check along the lines of Figs. 10–12) would close this gap.","section":"§III.B–C and Appendix C: seeding of the instability"}],"minor_comments":[{"comment":"The caption of Fig. 9 states θ = 0.12 (C = 0.6 > C_cr), but the text of §B.2 says the saddle points 'move off the imaginary axis as illustrated in Fig. 9 for θ = 0.6'. One of the two values is wrong (presumably C = 0.6 with θ = 0.12); please reconcile.","section":"Appendix B.2 / Fig. 9"},{"comment":"The text states that 'the crossed derivative term ∂ξ∂τ in Eq. (18) does not vanish for large γ_b'. Eq. (18) contains the mixed operator ∂²_ξ∂²_τ, not ∂ξ∂τ; a genuine ∂ξ∂τ term appears only after the e^{−iξ} envelope substitution leading to Eq. (21). Please reword for precision.","section":"§III.A, paragraph after Eq. (18)"},{"comment":"The sentence 'which is far more restrictive than the condition ξ/τ ≪ √(n_b/γ_b) associated with Eq. (20)' appears to misreference: Eq. (20) is the SVEA result of Ref. [56] whose validity condition is ξ/τ ≪ Q; the condition ξ/τ ≪ Γ_CFI ∼ √(n_b/γ_b) belongs to the solution of Eq. (19).","section":"§III.A"},{"comment":"Please state how the absolute normalization (vertical offset) of the theoretical curves relative to the PIC data is chosen in Figs. 2(b–d), 4(a), and 7. Also in Fig. 2 the caption symbols 'fEy²' and 'fBz²' appear to be formatting artifacts.","section":"Figs. 2, 4, 7"},{"comment":"The caption refers to 'asymptotic analytical predictions from ... (c) Eq. (18)' for the OTSI panel; presumably the asymptotic solution Eq. (29) is meant. Similarly, Fig. 4(a) labels the ordinate |F_y| while the text defines and uses F_⊥; please unify the notation.","section":"Fig. 7 caption; Fig. 4(a)"},{"comment":"The analysis is restricted to the x–y plane, with a claim of ready generalization to 3D. One or two sentences on whether the CFI–OTSI hierarchy and the ξ_tr ≃ 0.385 Γ_CFI τ boundary are unchanged when the transverse wavevector can take any azimuthal direction (and whether the k_y → ∞ SMI/HI connection acquires geometric factors in 3D) would strengthen §V.","section":"§II and §V: 2D restriction"},{"comment":"The TSI benchmark uses n_b/n_p = 0.6, which stretches the 'dilute beam' ordering underlying the fluid reduction. The good agreement with Eq. (37) is reassuring, but a brief remark on why the dilute-beam asymptotics remain quantitatively accurate at this density ratio would be useful.","section":"§IV.A, Fig. 5"}],"recommendation":"minor_revision","confidential_remarks":"Reference [65] is a companion paper by essentially the same group, \"submitted to Phys. Rev. Lett.\" and cited here for supporting material (e.g., the magnetically dominated character of the front-region modes). The editor may wish to verify that the division of results between the two submissions is clean and that neither depends on the other for its central claims. I found the self-citation pattern otherwise appropriate; the relation to the earlier SVEA-based CFI treatment [56] is discussed substantively rather than dismissively."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real news here is simple: drop the SVEA, keep full electromagnetics under the quasistatic approximation, and you get one master equation that covers CFI, OTSI, TSI, SMI, and HI. Near the beam front they find a spatiotemporal CFI that grows like exp(2√(Γ_CFI τ ξ)), with a sharp transition at ξ_tr ≃ 0.385 Γ_CFI τ into spatiotemporal OTSI. That front solution is not in Pathak et al. or their own earlier OTSI paper, and it is exactly where SVEA fails. That is the new result.\n\nThey do the work properly. Linearized cold-fluid + Maxwell to the PDE is transparent. Double Laplace + steepest descent is written out in the appendices, including the exact Bessel form in the ∂_ξ² ≫ 1 limit and the coalescence at the critical saddle. PIC runs (flat-top and Gaussian pair beams, plus narrow-beam SMI/HI) match growth laws and field polarization in the linear window. Growth rates enter only as coefficients; the exponents and the 0.385 factor are derived, not fitted.\n\nSoft spots are the usual ones and they own them. QSA forces ∂_τ ≪ ∂_ξ, so purely temporal OTSI is out and the domain is restricted to ξ ≲ τ/3. Cold, dilute, ultrarelativistic closures are standard for this class of problem; if those break (denser or mildly relativistic beams), the hierarchy can change. No public code is a minor annoyance, not a soundness issue. None of this undercuts the linear claim they actually make.\n\nThis is for people who model bounded relativistic beams in accelerators, laser-plasma transport, or dilute astrophysical pair beams and who still reach for temporal growth rates by default. The citation pattern is clean: they build on their OTSI work and on Pathak, then go past both. I would bring it to reading group, I would cite the front CFI solution and the master equation, and a serious editor should send it to referees. Engage with it.","headline":"Solid unified quasistatic theory that actually finds a new front-region spatiotemporal CFI and a clean CFI–OTSI handoff; math and PIC checks hold up.","tokens_in":35749,"tokens_out":540,"would_cite":true,"duration_ms":18100,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.35.Qz","52.40.Mj","52.65.Rr","52.27.Ny"],"model":"grok-4.5","headline":"A unified quasistatic model shows spatiotemporal current filamentation prevails near a relativistic beam front, and only later yields to oblique two-stream growth.","keywords":["beam-plasma instabilities","current filamentation","oblique two-stream","quasistatic approximation","spatiotemporal growth","self-modulation","hosing","ultrarelativistic beams"],"falsifier":"A high-resolution particle-in-cell run of a dilute ultrarelativistic pair beam should show magnetically dominated, non-oscillatory growth matching the spatiotemporal CFI formula near the front and electric-field-dominated oscillatory OTSI growth only past the predicted transition ξ_tr ≃ 0.385 Γ_CFI τ; failure of that spatial hierarchy would falsify the claim.","tokens_in":36010,"feed_emoji":"⚡","tokens_out":990,"duration_ms":25415,"temperature":0.7,"pith_summary":"When a bounded ultrarelativistic particle beam plows into fresh plasma, streaming instabilities do not grow at a single fixed rate everywhere. This paper builds one fully electromagnetic quasistatic equation that tracks every unstable mode along the beam without the slowly varying envelope approximation. It finds that a previously unreported spatiotemporal form of the current filamentation instability (CFI) dominates right at the front, while the oblique two-stream instability (OTSI) takes over only farther back. The same equation also recovers the self-modulation and hosing instabilities of long, narrow beams, tying those accelerator modes to the classical streaming family. Particle-in-cell simulations confirm the predicted hierarchy, giving a practical route to model dilute relativistic beams that standard temporal theory misrepresents.","feed_headline":"Beam-front filamentation wins before oblique two-stream takes over","feed_subtitle":"One quasistatic equation maps the full CFI–OTSI handoff and recovers accelerator hosing modes","key_machinery":"The quasistatic master equation for the perturbed plasma density (Eq. 15 / reduced form Eq. 18), obtained by decoupling slow beam evolution from fast plasma response in comoving coordinates and retaining full electromagnetic coupling without the slowly varying envelope approximation; it supplies the Green’s function whose asymptotics separate front CFI from rear OTSI and recover narrow-beam modes.","core_discovery":"Under conditions where standard temporal theory expects OTSI to dominate, a spatiotemporal CFI actually prevails in a front region ξ ≲ 0.385 Γ_CFI τ (precisely where the slowly varying envelope approximation fails) and is superseded by spatiotemporal OTSI only deeper in the beam; the same cold-fluid quasistatic master equation unifies this competition with longitudinal two-stream, self-modulation, and hosing modes.","pith_inferences":["If the front CFI region is robust, laboratory diagnostics that sample only the beam head may systematically report magnetic rather than electrostatic signatures even when bulk theory predicts OTSI.","The same master equation could be extended to mildly magnetized or warm beams to test whether the CFI–OTSI transition surface survives once thermal or cyclotron terms are restored.","Astrophysical pair cascades from TeV blazars may filament over a larger fraction of the beam head than pure temporal CFI estimates suggested, altering the competition with inverse-Compton cooling."],"forward_implications":["Near-front magnetic filaments of length ~0.4 N_e skin depths are generic for dilute ultrarelativistic beams, independent of small density ratio or large Lorentz factor.","Short Gaussian beams (σ_x ~ plasma skin depth) can be treated with the same equation once Γ_CFI is allowed to vary with local density.","Self-modulation and hosing share the identical spatiotemporal envelope equation as OTSI in the large-transverse-wavenumber limit, so narrow-beam accelerator instabilities are continuous with classical streaming modes.","Quasistatic particle-in-cell codes become a viable large-scale tool for dilute astrophysical pair beams whose dynamical timescales are inaccessible to standard PIC."],"fun_headline_variants":["Spatiotemporal CFI leads at beam front before OTSI overtakes","Quasistatic model shows CFI first, OTSI deeper in relativistic beams","Front-region filamentation precedes oblique two-stream takeover","One master equation unifies CFI–OTSI handoff and hosing modes","SVEA fails where spatiotemporal CFI dominates beam front"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"Plasma and field quantities must vary much faster along the beam than in time, so the analysis is restricted to a limited region behind the front and cannot capture purely temporal oblique growth farther back.","fun_headline_variants_meta":{"raw":{"variants":["Spatiotemporal CFI leads at beam front before OTSI overtakes","Quasistatic model shows CFI first, OTSI deeper in relativistic beams","Front-region filamentation precedes oblique two-stream takeover","One master equation unifies CFI–OTSI handoff and hosing modes","SVEA fails where spatiotemporal CFI dominates beam front"]},"model":"grok-4.5","effort":"low","cost_usd":0.003796,"raw_usage":{"total_tokens":1253,"prompt_tokens":876,"num_sources_used":0,"completion_tokens":96,"cost_in_usd_ticks":37964000,"prompt_tokens_details":{"text_tokens":876,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":281,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":876,"tokens_out":96,"duration_ms":5715,"temperature":1.0,"reasoning_tokens":281,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T20:21:47.697307+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A high-resolution particle-in-cell run of a dilute ultrarelativistic pair beam should show magnetically dominated, non-oscillatory growth matching the spatiotemporal CFI formula near the front and electric-field-dominated oscillatory OTSI growth only past the predicted transition ξ_tr ≃ 0.385 Γ_CFI τ; failure of that spatial hierarchy would falsify the claim.","supporting_citations":[],"review_version":1}