{"id":"195b6254-8a48-4377-ae79-d1391cef9a35","arxiv_id":"2607.24214","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"Quasistatic electromagnetic theory and QS-PIC simulations show spatiotemporal current filamentation dominates near an ultrarelativistic beam front while OTSI dominates downstream, enabling blazar-regime modeling.","lead":"A quasistatic theory shows that magnetic filamentation, not the oblique two-stream mode, dominates right at the front of an ultrarelativistic beam entering fresh plasma. The same approximation lets particle simulations reach blazar-jet densities that ordinary kinetic codes cannot touch.","discovery_kind":"new_method","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The linear core (Eq. 5, its asymptotes, and the PIC/QS-PIC cross-check) is solid; the least-secured element of the strongest claim is the asserted αb/γb-independence of the CFI–OTSI transition locus, which rests on unshown QuaSSis sweeps anchored to full PIC at a single parameter point.","rationale":"I agree with the reader's ACCEPT/HIGH and with the assessment that the linear central claim is well supported by analytics plus two independent numerical methods, with the one input parameter (Γ via measured ky) honestly accounted for. My concern differs in locus from the reader's weakest_assumption: the reader flagged cold-QSA validity into the deep nonlinear blazar stage, which is real but attaches to the demonstration run (Fig. 3), not to the strongest claim, which is explicitly linear. Within the strongest claim itself, the only piece whose evidence base is materially thinner than the rest is the αb/γb-independence of the transition locus: shown at one parameter point, asserted (without displayed data) across eight decades in αb, and dependent on a QS-PIC code benchmarked at that same single point. This is a low-risk gap rather than an inconsistency, because the QSA ordering strictly improves in the direction of extrapolation and the normalized-locus claim is parameter-free given Eq. (5). Hence UNCHANGED: the verdict stays ACCEPT, and the proposed second-anchor test is a strengthening check rather than a condition. No code is shipped for QuaSSis, which the reader already noted; the proposed test is executable with the described calder setup and does not require the new code to verify the physics claim itself.","tokens_in":13505,"tokens_out":11611,"duration_ms":360690,"concrete_test":"Run one additional full-PIC simulation at a second parameter point still within full-PIC reach, e.g. αb=10⁻³, γb=10³ (ωp/ωb≈700, below the stated ≲1000 limit), plus a matched QuaSSis run. Measure the Bz/Ey transition locus in (kpξ, Γτ) variables and the filament extent at fixed e-fold number Ne. If the transition sits at kpξ≈2Γτ/√27 (extent ≈0.4Ne kp⁻¹) and QuaSSis reproduces full PIC, the universality claim gains a second independent anchor spanning ~2 decades in ωp/ωb; a systematic drift with αb/γb would qualify the claimed refutation of Ref. [32].","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation chain Eqs. (1)–(5) checks out: the QSA drops only τ-derivatives (the ξ-derivatives cancel exactly in the d'Alembertian in co-moving coordinates), the neglect of γ_b^{-2}∂²ξ against k_y² is satisfied by many orders of magnitude in both runs, and the saddle analysis of Eq. (5) self-consistently yields the ST-OTSI pole at K≈ikp and the purely growing ST-CFI saddle at K=√(kpΓτ/ξ). The full-solution comparison in Fig. 1(d,e) — amplitude and phase across the whole ξ–τ window, with Γ=1.7×10⁻³ ωp computed (not fitted) from αb, γb and the measured ky/kp≈5.5 — is genuinely strong validation of the growth-law claim, so I do not press the usual \"τ^{1/2} vs τ^{2/3} exponent discrimination\" worry.\n\nThe soft spot is the last clause of the strongest claim: the transition locus kpξ ≲ 2Γτ/√27 being \"independent of αb and γb.\" Mathematically this follows trivially once Eq. (5) is written in (kpξ, Γτ) variables, so the physical content is entirely empirical: that real beam-plasma systems at widely different αb, γb actually follow Eq. (5) with the same normalized transition. Direct full-PIC evidence exists at exactly one point (αb=0.03, γb=2×10⁴). The wide-range verification (αb∈[10⁻¹⁰,10⁻²], γb∈[10²,10⁶]) is stated in one sentence — \"We have further verified it through additional simulations\" — with no data shown, and those runs necessarily use QuaSSis, which is benchmarked against full PIC only at that same single point. The claim is also what carries the paper's refutation of the SVEA-based confinement scaling of Ref. [32], so it bears real weight. Mitigation: the extrapolation direction (smaller αb, larger γb) is precisely where the QSA ordering ωbs≪ωp improves, so the risk is low — but it is presently an extrapolation from one anchor point, not a demonstrated universality.","agreement_with_reader":"partial"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript develops a linear, fully electromagnetic, two-dimensional cold-fluid theory of ultrarelativistic beam-plasma instabilities using the quasistatic approximation without the slowly varying envelope approximation. The central result, Eq. (5), admits two spatiotemporal regimes: an OTSI-like mode downstream, Eq. (6), and a purely growing current-filamentation mode near the beam front, Eq. (7), with a normalized transition near \\(k_p\\xi\\simeq2\\Gamma\\tau/\\sqrt{27}\\). A full-PIC calculation at \\(\\alpha_b=0.03\\), \\(\\gamma_b=2\\times10^4\\), an independent QS-PIC calculation, and the numerical solution of Eq. (5) are compared in amplitude, phase, and field polarization. The authors then use the QuaSSis code to follow an \\(\\alpha_b=10^{-10}\\), \\(\\gamma_b=10^6\\) pair beam into the nonlinear regime, interpreting saturation through beam trapping and the subsequent front-filament wakefields.","tokens_in":14015,"tokens_out":8776,"duration_ms":268275,"significance":"If supported, this is a substantial and useful result: it provides a unified spatiotemporal description of CFI and OTSI at the beam front, directly challenging an earlier SVEA-based conclusion that ultrarelativistic spatiotemporal CFI is unimportant. The paper has several concrete strengths: a compact derivation avoiding the SVEA, closed-form asymptotic solutions, a one-to-one comparison against full PIC and an independently implemented QS-PIC code, and a computationally efficient route to a blazar-relevant regime inaccessible to conventional PIC. The value of the \\(\\Gamma\\) used in the main comparison is calculated rather than freely fitted, although it depends on the measured \\(k_y\\). The normalized transition prediction is falsifiable and practically useful. The wide-parameter evidence for that prediction and the reproducibility of the extreme nonlinear QS-PIC run should be documented more fully.","major_comments":[{"comment":"Analytical-solutions section, following Eq. (7): the claimed independence of the normalized CFI extent from \\(\\alpha_b\\) and \\(\\gamma_b\\), and the resulting invalidation of Ref. [32], are supported in the manuscript only by the sentence “We have further verified it through additional simulations.” Once the data are plotted against \\(k_p\\xi\\) and \\(\\Gamma\\tau\\), Eq. (5) makes the normalized locus mathematically expected; the empirical issue is whether simulations across \\(\\alpha_b=10^{-10}\\)–\\(10^{-2}\\), \\(\\gamma_b=10^2\\)–\\(10^6\\) follow it. Please show a compact verification—e.g. rescaled transition/filament extent versus e-fold number for representative points—identify which runs use full PIC or QuaSSis, and state convergence. Otherwise this specific universality claim should be narrowed.","section":"Analytical solutions"},{"comment":"Extending the QSA section and Fig. 3: the key nonlinear demonstration at \\(\\alpha_b=10^{-10}\\), \\(\\gamma_b=10^6\\) gives almost no numerical information. Mesh size, time step, domain, macroparticle numbers, noise initialization, moving-window treatment, and resolution of the \\(\\sim k_p^{-1}\\) filaments are not reported. Because direct full-PIC verification is impossible at these parameters, the conclusion that QuaSSis captures filament saturation, coalescence, and the \\(E_x\\) jump near \\(\\omega_p\\tau\\simeq2.4\\times10^{10}\\) rests heavily on numerical trustworthiness. Please add a parameter table and at least a concise resolution/particle-count convergence check for the saturation fields and the coalescence event, and state the QSA ordering after transverse heating.","section":"Extending the QSA to extreme blazar jet regimes"}],"minor_comments":[{"comment":"PIC-simulations section: explain how the dominant value \\(k_y/k_p\\simeq5.5\\) was extracted (for example, from the early-time Fourier spectrum) and how sensitive the Fig. 1(e) comparison is to this choice. Although \\(\\Gamma\\) is not freely fitted, it is conditional on this measured wavenumber.","section":"PIC simulations"},{"comment":"Following Eq. (7): define precisely what is meant by a “temporal e-fold” \\(N_e\\) in the estimate that CFI extends over \\(\\sim0.4N_e k_p^{-1}\\). The relevant local spatiotemporal growth differs from both \\(\\Gamma\\) and the downstream ST-OTSI exponent.","section":"Analytical solutions"},{"comment":"Fig. 1(d,e): specify the normalization and color scale of \\(|F_\\perp|\\), and state quantitatively how the dashed nonlinear-regime boundary in panel (d) was identified.","section":"Figure 1"},{"comment":"Fig. 3 and surrounding text: define whether \\(n_b\\) in the saturation-energy estimates denotes the density per pair species or the total pair density, and clarify why \\(k_y\\simeq k_p\\) is appropriate for the extreme run given the different dominant wavenumber in the Fig. 1 benchmark.","section":"Extending the QSA to extreme blazar jet regimes"},{"comment":"The phrases “all spatiotemporal scales” in the Introduction and “identical field evolution” in the QuaSSis benchmark are broader than the evidence shown. Qualifying them by the QSA ordering and the demonstrated parameter point would make the claims more precise.","section":"Introduction / QuaSSis benchmark"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript appears appropriate in scope for a concise physics Letter. Some technical ingredients—including details of the saddle analysis, generalized beam profiles, and robustness to distributed seeds—are deferred to companion Ref. [41]. It would be useful for the editorial process to ensure that this companion manuscript is available to reviewers and that no evidence essential to the Letter's central claims exists only there."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The core result is real: under a fully electromagnetic quasistatic model that drops the slowly-varying-envelope approximation, spatiotemporal CFI owns the first few skin depths of the beam front (growth ~ exp(2√(Γτ kp ξ))), while spatiotemporal OTSI takes over farther back. That hierarchy, and the transition locus kp ξ ≲ 2Γτ/√27, is the new piece. It directly contradicts the inductive-SVEA conclusion in Pathak et al. that CFI becomes irrelevant at large γb.\n\nWhat they did well is the cross-check. Equation (5) plus the two saddle-point asymptotes are clean. Figure 1 matches full PIC in both amplitude and phase across the whole ξ–τ window, with Γ taken from the measured ky, not fitted. The Ey-vs-Bz polarization split in Figure 2 is exactly what the theory predicts. The new QS-PIC code then reproduces the same linear window at a tiny fraction of the cost and lets them reach αb = 10^{-10}, γb = 10^6 — a regime ordinary PIC cannot touch. The nonlinear saturation estimates (bounce frequency matching the local spatiotemporal growth rate) line up with the runs they show.\n\nThe soft spot is proportionate, not fatal. The claim that the transition locus is independent of αb and γb follows mathematically once you nondimensionalize Eq. (5), but the physical verification is one full-PIC anchor point plus a one-sentence statement that QuaSSis sweeps cover the rest. Those sweeps are not shown, and QuaSSis itself is only benchmarked against full PIC at that same point. The extrapolation direction (smaller αb, larger γb) is where the QSA ordering gets better, so the risk is low, but the refutation of Pathak rests partly on that unshown sweep. Cold-fluid QSA into deep nonlinear filament coalescence is also an assumption, not a proof; it looks plausible from the scalings they recover, but it is not rigorously closed here. Code is not shipped.\n\nThis is for people who actually run or analyze relativistic beam-plasma systems — accelerator beams with finite length, or blazar pair cascades. The linear theory is tight enough that I would cite it; the extreme-parameter demonstration is the part that changes what is computable. Send it to referees. It deserves the time.","headline":"Solid linear EM-QSA theory that cleanly puts CFI at the beam front and OTSI downstream; the extreme blazar runs are the real payoff, with one under-shown universality claim.","tokens_in":14788,"tokens_out":593,"would_cite":true,"duration_ms":16154,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Near an ultrarelativistic beam front, spatiotemporal current filamentation wins over oblique two-stream growth, and a quasistatic model captures both.","keywords":["streaming instabilities","current filamentation","oblique two-stream","quasistatic approximation","ultrarelativistic beams","blazar jets","particle-in-cell","spatiotemporal growth"],"falsifier":"A controlled PIC or laboratory measurement in which, for a cold ultrarelativistic beam, the magnetic filaments fail to appear within a few skin depths of the front, or the front-to-downstream transition distance grows or shrinks with beam density or Lorentz factor instead of staying proportional only to the number of e-folds.","tokens_in":14486,"feed_emoji":"⚡","tokens_out":1024,"duration_ms":21749,"temperature":0.7,"pith_summary":"Relativistic particle beams streaming through plasma drive instabilities that amplify fields and reshape the beam, but when the beam is tenuous and ultrarelativistic the timescales of beam and plasma differ by many orders of magnitude, making full kinetic simulation impractical. This paper shows that the quasistatic approximation yields a single linear electromagnetic equation that tracks the full unstable spectrum without assuming a slowly varying envelope. Solving that equation reveals a previously unreported hierarchy: magnetic filamentation dominates within a few plasma skin depths of the beam front, while the oblique two-stream mode takes over farther downstream. The same approximation powers a particle-in-cell code that reaches the extreme density ratios of blazar pair beams and follows the instabilities into the nonlinear filament-coalescence regime. Agreement among the analytic scalings, ordinary PIC runs, and the quasistatic code supports the claim that the front-side filaments set the wakefields that later heat the beam bulk.","feed_headline":"Filamentation wins at the beam front; quasistatics make it computable","feed_subtitle":"A single electromagnetic equation and a fast PIC code reach blazar-pair densities once thought unreachable","key_machinery":"The quasistatic equation for the pseudo-potential Psi (Eq. 5), obtained by changing to co-moving coordinates and dropping slow laboratory-time derivatives of plasma and fields; it unifies electrostatic and inductive modes without a slowly-varying-envelope assumption and supplies both the analytic saddle-point solutions and the advance used by the quasistatic PIC code.","core_discovery":"In the linear, fully electromagnetic quasistatic description of a cold ultrarelativistic beam that continuously meets fresh plasma, spatiotemporal current filamentation grows as exp(2 sqrt(Gamma tau k_p xi)) and dominates within a few skin depths of the front, while spatiotemporal oblique two-stream growth dominates farther back; the transition locus scales as k_p xi less than or similar to 2 Gamma tau / sqrt(27) and is independent of beam density ratio and Lorentz factor.","pith_inferences":["If front-side filaments control the later wakefield spectrum, blazar-pair-beam models that omit the beam head may systematically mis-estimate the residual GeV cascade and the residual intergalactic magnetic field.","The density- and gamma-independent transition length suggests a simple experimental diagnostic: measure the axial extent of magnetic filaments after a fixed number of e-folds; departure from ~0.4 N_e skin depths would signal warm-beam or ion-motion corrections.","Once ions are mobile or an external guide field is present, the same quasistatic machinery should reveal whether the front-side CFI window shrinks or expands, offering a direct test of suppression scenarios proposed for blazar beams."],"forward_implications":["Laboratory beams with finite length will show magnetic filamentation confined to the head and electrostatic chevron patterns farther back, even when temporal theory predicts pure OTSI dominance.","Quasistatic PIC can now reach blazar-relevant density ratios (alpha_b ~ 10^{-10}, gamma_b ~ 10^6) and follow the system past saturation into filament coalescence and wakefield-driven longitudinal heating.","Downstream bulk heating after OTSI saturation is set by wakefields driven by the short, CFI-generated filaments at the front rather than by bulk temporal CFI alone.","The same framework extends immediately to longitudinal modes (two-stream, self-modulation, hosing) treated in the companion paper."],"fun_headline_variants":["Filamentation dominates near beam front in quasistatic model","Spatiotemporal CFI overtakes oblique mode within skin depths of front","QSA reveals CFI growth exp(2 sqrt(Gamma tau k_p xi)) at beam edge","Filamentation leads at front; oblique two-stream farther back","Quasistatic EM model captures full unstable spectrum for dilute beams"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The cold-fluid closure and the ordering that plasma and field profiles change much faster along the beam than in laboratory time must remain valid from linear growth all the way through deep nonlinear filament coalescence.","fun_headline_variants_meta":{"raw":{"variants":["Filamentation dominates near beam front in quasistatic model","Spatiotemporal CFI overtakes oblique mode within skin depths of front","QSA reveals CFI growth exp(2 sqrt(Gamma tau k_p xi)) at beam edge","Filamentation leads at front; oblique two-stream farther back","Quasistatic EM model captures full unstable spectrum for dilute beams"]},"model":"grok-4.5","effort":"low","cost_usd":0.003549,"raw_usage":{"total_tokens":1190,"prompt_tokens":789,"num_sources_used":0,"completion_tokens":81,"cost_in_usd_ticks":35488000,"prompt_tokens_details":{"text_tokens":789,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":320,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":789,"tokens_out":81,"duration_ms":6183,"temperature":1.0,"reasoning_tokens":320,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T20:47:12.988530+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A controlled PIC or laboratory measurement in which, for a cold ultrarelativistic beam, the magnetic filaments fail to appear within a few skin depths of the front, or the front-to-downstream transition distance grows or shrinks with beam density or Lorentz factor instead of staying proportional only to the number of e-folds.","supporting_citations":[],"review_version":1}