{"id":"22a4d3e7-bac3-4a99-b2e4-3c594b8ee465","arxiv_id":"2608.12606","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A near-Sun interplanetary shock's upstream foreshock resolves into four wave families, showing a proton-beam-driven scattering field and a weak compressive fast-mode component that changes resonance energies by up to 13%.","lead":"Parker Solar Probe crossed a very fast, near-parallel shock in March 2023, and the authors separated its upstream wave field into four distinct families for the first time. The data show the shock-accelerated protons drive the very waves that scatter them, with a weak compressive component shifting the resonance energies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The precursor scale L_EP rests on a stationarity assumption that the paper asserts but does not test; if L_EP is a temporal rise rather than a spatial diffusion scale, the central self-regulation comparison loses its quantitative anchor.","rationale":"The reader's verdict is CONDITIONAL, anchored on the LP-OB density-proxy identification. That is a fair concern for the compressive component and the 13% energy-shift sub-claim, but it is outside the self-regulation loop identified as the paper's central claim. The loop is supported by the match between the quasi-linear mean free path from the resonant wave amplitude and the empirical mean free path inferred from L_EP. That inference assumes the upstream energetic-particle profile is stationary in the shock frame. A single-spacecraft time series cannot distinguish a stationary spatial precursor from a time-dependent source, and the fit e-folding time (~26 min) is comparable to the interval over which the shock-accelerated population may still be building. The split-half and energy-resolved re-fit proposed above is a direct, computational check; if the profile is stationary, the central claim holds as stated and the CONDITIONAL verdict stands. If not, the empirical half of the comparison fails, though the quasi-linear estimate itself remains; the paper would need an alternative demonstration that the beam remains anisotropic. In either case, this does not warrant rejection at this stage, so the reader's CONDITIONAL verdict should remain unchanged.","tokens_in":30206,"tokens_out":17329,"duration_ms":189696,"concrete_test":"Split the inner-foreshock interval 06:15-07:13 UT into two independent halves and re-fit PEP(x) in each; if the two e-folding lengths differ by more than ~30%, the exponential time profile is not a stationary spatial precursor. In the same pass, fit L_EP separately for EPI-Lo energy channels near 0.3, 0.5, 1, 2, and 4 MeV: a stationary diffusive precursor should show L(E) = kappa(E)/u1 increasing with energy according to the observed resonant wave spectrum, whereas a temporal rise would produce speed-ordered profiles that break the common x mapping. If both checks are consistent with a stationary, energy-dependent diffusion length, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative step is the comparison in Section 5 between the quasi-linear mean free path, lambda_parallel ~ r_g (B0/deltaB)^2 (Eq. 6), and an empirical lambda_parallel = 3 kappa_eff / v derived from kappa_eff = u1 L_EP, with L_EP ~ 5.48 R_sun obtained by fitting PEP(t) against x = u1 (t_shock - t) (Appendix B). This time-to-space conversion is valid only if the precursor is stationary in the shock frame. The text asserts stationarity ('Fully developed means stationary'), but the only support offered is the constancy of the same time series that defines the exponential profile; a spatially uniform but temporally rising energetic-particle population at the shock also produces an exponential time profile with a constant e-folding time. On the ~26 min e-folding timescale of the fit, time-dependent injection in a CME-driven shock is not excluded. If L_EP is a temporal growth scale rather than a spatial diffusion scale, the empirical mean free path is not measured, and the claimed 'half the precursor scale' agreement with Eq. (6) does not test the self-regulation loop. The reader's LP-OB density-proxy concern, while valid, affects the peripheral compressive component and the 13% shift, not the core loop.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents Parker Solar Probe observations of the fast, near-parallel interplanetary shock of 2023 March 13 at 0.24 AU. Using a wavelet decomposition of the upstream magnetic field, the authors separate the foreshock wave field into four families: right-hand and left-hand circularly polarized waves, a field-aligned linear family, and an oblique linear family. They identify the first three as cyclotron-resonant with backstreaming suprathermal-to-MeV protons and the oblique family as a non-resonant, compressive fast-magnetosonic component. The central quantitative step is a comparison in Section 5 between a quasi-linear parallel mean free path lambda_par ~ r_g (B0/deltaB)^2 and an empirical lambda_par = 3 kappa_eff / v with kappa_eff = u1 L_EP from an exponential fit to the energetic-particle pressure precursor. The agreement (lambda_par ~ 0.7-1.7 R_sun versus 1.9-2.9 R_sun) is interpreted as evidence that the beam drives the waves that scatter it, leaving the beam anisotropic enough to sustain the drive.","tokens_in":30385,"tokens_out":13941,"duration_ms":126219,"significance":"If the interpretation holds, this is the first in situ resolution of the self-regulated foreshock of a fast near-parallel shock close to the Sun, and it provides a concrete, quantitative test of the quasi-linear resonant-scattering loop that underlies diffusive shock acceleration. The paper's strengths include the use of independent measurements for the polarization/wavenumber classification, the hodogram and compressibility checks, the explicit derivation of the quasi-linear estimate from measured B0 and deltaB_res rather than from the precursor fit, and the detailed appendices documenting the single-spacecraft estimators. The qualitative picture of beam-driven RH and LH waves over a common band, with both helicities present, is credible and internally consistent. The quantitative self-regulation claim, however, rests on the time-to-space conversion x = u1(tshock - t), which requires the precursor to be stationary in the shock frame; this assumption is asserted but not independently tested. The compressive fast-mode identification also depends on the spacecraft-potential density proxy, for which a direct phase calibration at the LP-OB frequencies is not provided.","major_comments":[{"comment":"The central comparison in Section 5, between the quasi-linear lambda_par ~ 0.7-1.7 R_sun from Eq. (6) and the empirical lambda_par = 3 kappa_eff / v ~ 1.9-2.9 R_sun, uses kappa_eff = u1 L_EP with L_EP obtained by fitting PEP(t) against x = u1(tshock - t) in Appendix B. This time-to-space conversion is valid only if the upstream precursor is stationary in the shock frame. The paper asserts stationarity in Section 2 ('Fully developed means stationary: the family statistics and the pressure gradient hold steady'), but the support offered is the constancy of the same time series that defines the exponential profile; a spatially uniform but temporally rising energetic-particle population at the shock would produce an exponential time profile with the same constant e-folding time. The single-spacecraft geometry cannot distinguish these cases, and the 'family statistics hold steady' statement does not break the degeneracy. Because the self-regulation claim is quantitatively anchored by this comparison, the manuscript should either provide an independent stationarity test (for example, energy-dependent e-folding scales from EPI-Lo channels, which should scale with kappa(E)/u for a spatial diffusion profile) or explicitly reframe the comparison as conditional on stationarity. Without such a test, the statement 'the measured mean free path, half the precursor scale' is not supported.","section":"Appendix A.4 and Appendix C"},{"comment":"The identification of the LP-OB family as fast magnetosonic rests on the density-field cross-phase Phi(delta_n, delta_|B|) ~ 0, where delta_n is derived from the spacecraft floating potential (Appendix A.4). The paper states that the potential method fails above tau_c^{-1} >~ 10 kHz, well above the 0.003-0.03 Hz LP-OB band, but no measured tau_c for this interval is reported and no direct cross-calibration of the potential-density fluctuations against quasi-thermal-noise density within the LP-OB band is shown. The measured amplitude ratio |delta_n/n| / |delta_B_par/B0| ~ 1.4 (Appendix C) is attributed to the 'multiplicative uncertainty' of the density calibration, which indicates that amplitude fidelity is not independently established. If the potential response has a frequency-dependent phase or contamination, the in-phase signature and the derived 3-13% resonance-energy shift would not be uniquely attributable to fast magnetosonic compression. Please provide a quantitative phase calibration of the spacecraft-potential density at these frequencies or an independent density fluctuation measurement over the LP-OB band.","section":"Appendix A.4 and Appendix C"}],"minor_comments":[{"comment":"The multi-panel Figures 2 and 3 are extremely dense; the small panel labels and overlaid distributions make it difficult to verify the family separations visually. Consider publishing full-resolution versions or splitting the spectrograms and the occurrence/statistics panels into separate figures.","section":"Figures 2 and 3"},{"comment":"The split of the linearly polarized class at theta_kB = 45 degrees is presented without a physical justification or a sensitivity test. Because the LP-OB family is defined by this threshold, it would strengthen the classification to show how the power fractions and the cross-phase separation vary when the threshold is moved within, say, 35-55 degrees.","section":"Section 3"},{"comment":"The paper acknowledges that the LP-FA family could be a coherent superposition of RH and LH packets rather than an independent linear mode, but the branch assignment that follows 'rests on the drive' and is not independently verified. Please state explicitly whether the physical conclusions require LP-FA to be an independent resonant family, or clarify how the drive argument distinguishes a directly excited shear-Alfven mode from a superposition artifact.","section":"Section 4"},{"comment":"Several key references are to works that are submitted, in press, or arXiv-only (for example, Kouloumvakos et al. 2026, Capanema et al. 2026, Giacalone et al. 2026, Raptis et al. 2026). The claims that rely on these should be clearly marked as preliminary, and the companion paper should be cited with its current status.","section":"References"},{"comment":"The abstract states that the measured mean free path is 'half the precursor scale,' while Section 5 gives ranges of 0.7-1.7 R_sun versus 1.9-2.9 R_sun, i.e., a factor between roughly 1.1 and 4. Please ensure the wording reflects the range and the associated uncertainties.","section":"Abstract and Section 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the observational analysis is careful, with a credible qualitative case for a beam-driven, self-regulated foreshock. The main risk is the stationarity assumption behind L_EP; because PSP is a single spacecraft, this cannot be fully resolved with the present data, but the authors can either provide an energy-dependent e-folding analysis or soften the quantitative claim. The density-proxy issue is secondary but deserves a quantitative response. I would support publication after a major revision that addresses the stationarity point and clarifies the LP-FA classification status."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first in-situ decomposition of a fast sub-AU shock foreshock into four wave families—RH, LH, LP-FA, and a resolved oblique compressive LP-OB component. That part is real and worth refereeing. The resonant-scattering loop is credible in outline, but the central quantitative comparison, the mean free path against L_EP, is weaker than the paper admits because stationarity is asserted rather than tested.\n\nWhat is actually new: the polarization-resolved four-family classification, with near-equal RH and LH power over a common wavenumber band, is not in the ISEE-3 or bow-shock literature. The tools are mostly standard (wavelet decomposition, SVD polarization, cross-phase), but the observational result is new. The LP-OB identification is the most delicate and the authors are honest about its limits: they check coherence against a phase-randomized control, show the in-phase density–field signature, and quantify why the known local drivers fall short. Their exclusion analysis is quantitative rather than hand-waving. The near-equal RH and LH power is also physically interesting, since it is what full pitch-angle scattering needs.\n\nSoft spots, in proportion: the stationarity assumption for L_EP is the load-bearing one. The fit of PEP(t) against x = u1(tshock − t) is valid only if the precursor is stationary in the shock frame. The paper says “fully developed means stationary,” but that is supported only by the constancy of the same time series that defines the exponential profile. A temporally rising injection in a CME-driven shock would also produce an exponential time profile with a constant e-folding time. If L_EP is a temporal growth scale rather than a spatial diffusion scale, then the empirical mean free path is not measured, and the “half the precursor scale” agreement with the quasi-linear estimate does not test the self-regulation loop. This is not a peripheral concern; it is the quantitative anchor of the headline claim. I would want either a sub-interval stability test of the fit, or a direct statement that the comparison is order-of-magnitude only.\n\nThe LP-OB density-proxy concern is real but secondary. A frequency-dependent phase lag in the spacecraft-potential density would blur the fast-mode identity, and the authors themselves note that a second spacecraft would settle it. The 13% resonance-energy shift inherits this uncertainty, but it is a minor part of the paper. Error bars are missing from Table 1 and the L_EP fit, and no code or data version identifiers are given; these are addressable.\n\nThe paper deserves a serious referee. It is a single-event study with a strong claim, and the stress-test concern is the one a referee should push on. The wave-family classification and the resonant loop do not collapse if L_EP is uncertain, but the quantitative headline needs defense or softening.\n\nRecommendation: send to peer review, with referees who understand single-spacecraft time-series ambiguities.","headline":"A genuinely new four-family decomposition of a fast PSP shock foreshock, with a credible resonant-scattering loop whose quantitative anchor (L_EP) rests on an untested stationarity assumption.","tokens_in":31042,"tokens_out":2300,"would_cite":true,"duration_ms":24599,"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":"The upstream foreshock of a near-parallel shock at 0.24 AU is self-regulated: the beam from the shock drives the waves that scatter it, and the measured mean free path is about half the precursor scale.","keywords":["interplanetary shocks","solar wind","space plasmas","plasma astrophysics","Alfvén waves","foreshock","diffusive shock acceleration","Parker Solar Probe"],"falsifier":"Recompute the density–field cross-phase in the LP-OB band using the quasi-thermal-noise density estimate from the same spacecraft where it is available at comparable resolution; if the near-zero phase difference between density and magnetic-field magnitude is not reproduced, the fast-magnetosonic identification fails. For the self-regulation loop, run a kinetic simulation of a θ_Bn ≈ 8°, M_A ≈ 7.5 shock with a self-consistently driven foreshock and check that the resonant-band transverse amplitude yields a quasi-linear mean free path within a factor of a few of the measured 0.7–2.9 R_sun and that the RH and LH power stays within a factor of two over the shared wavenumber band.","tokens_in":29947,"feed_emoji":"☀️","tokens_out":11733,"duration_ms":98085,"temperature":0.7,"pith_summary":"Parker Solar Probe's crossing of a fast (~2800 km/s), near-parallel interplanetary shock at 0.24 AU on 2023 March 13 provides an in situ look at the self-regulation that diffusive shock acceleration assumes. The paper separates the upstream wave field into four families — right-hand and left-hand circular, field-aligned linear, and oblique linear — and argues that the first three are excited by the backstreaming suprathermal-to-MeV proton beam through cyclotron resonance, while scattering that same beam. The measured parallel mean free path, roughly 0.7–2.9 solar radii against a precursor scale of about 5.5 solar radii, leaves the beam anisotropic enough to keep driving the waves. The weak oblique family is compressive and fast magnetosonic, and its field-magnitude modulation shifts the resonance energies of the scattering families by up to 13% along the precursor. If the interpretation is right, acceleration at shocks inside 0.3 AU is set in a foreshock the shock builds for itself.","feed_headline":"The protons a shock flings back create the waves that scatter them","feed_subtitle":"In situ data from 0.24 AU show the foreshock is self-regulated, with mean free path about half the precursor scale.","key_machinery":"The central identity is the cyclotron-resonance condition ω − k∥v∥ = ±Ωci, which assigns each measured wavenumber to a proton energy E_res ≈ ½m_p(Ωci/k∥)² and, together with the quasi-linear mean-free-path estimate λ∥ ∼ r_g (B0/δB)², closes the loop between the beam and the scattering field. The family separation itself is carried by a Morlet-wavelet spectral-matrix analysis that returns signed ellipticity, propagation angle, wavenumber, and the density–field cross-phase for every time–frequency bin. The in-phase density–field cross-phase is the identifying test that places the compressive part of the oblique linear family on the fast magnetosonic branch.","core_discovery":"On the paper's own terms, the discovery is that the foreshock of the 2023 March 13 shock is a closed loop. The suprathermal-to-MeV protons streaming back from the ramp drive three cyclotron-resonant wave families — right-hand circular (R-mode, ~36% of power), left-hand circular (L-mode, ~37%), and field-aligned linear (~21%) — and those same families scatter the beam with a parallel mean free path λ∥ ≈ 0.7–1.7 R_sun from the quasi-linear estimate (up to 2.9 R_sun from the empirical effective diffusion coefficient), comparable to the measured precursor scale L_EP ≈ 5.48 R_sun. Because only two to three scattering lengths span the precursor, the beam isotropizes only partially and retains the field-aligned drift that feeds the instability. Outside this loop sits the remaining ~6%: an oblique, linearly polarized, compressive population that is fast magnetosonic by its in-phase density–field cross-phase, is not produced by any local beam-driven channel, and modulates the resonant energies of the kinetic families by 3–13% through its δB∥/B0 fluctuations.","pith_inferences":["A testable extension: apply the same four-family decomposition to other Parker Solar Probe shock crossings with different Mach numbers and check whether the ratio λ∥/L_EP stays near 0.1–0.5; if it does, the self-regulation condition is a universal attractor rather than a property of this event.","If the LP-OB family is indeed pre-existing fast-mode turbulence advected from the solar wind, its amplitude just upstream of a shock is a probe of the ambient seed-fluctuation level rather than a shock product; comparing LP-OB amplitude to quiet-time solar-wind fast-mode levels would separate the two.","The 3–13% resonance-energy modulation implies the injection energy for diffusive shock acceleration varies with distance upstream; using a fixed resonance energy in transport models would misplace the spectral cutoff.","The near-equal RH and LH power implies the backstreaming beam's backward pitch-angle hemisphere is substantially populated; measured pitch-angle distributions of 0.3–1 MeV protons should show a filled counter-streaming component rather than a narrow beam."],"forward_implications":["The foreshock of a near-parallel shock inside 0.3 AU is self-regulating: the backstreaming proton beam supplies the scattering that controls its own diffusion, so no external wave-amplitude prescription is needed to model upstream transport.","The near-equal power of the two circularly polarized families over a common wavenumber band means both pitch-angle hemispheres of the beam are scattered, which is what allows the distribution to stay partly anisotropic while still feeding the instability.","The measured mean free path of roughly 0.7–2.9 R_sun against a precursor scale of 5.48 R_sun puts the beam in a diffusive but not fully trapped regime; two or three scattering lengths span the precursor, consistent with the observed exponential pressure profile.","The compressive fast-magnetosonic component, though only a few percent of the wave power, shifts the resonant energies of the scattering families by up to 13% along the precursor, so the energy that resonates with a given wave changes as the shock approaches.","At stronger or longer-driven shocks, where the precursor pressure ratio approaches order unity, the resonant-band amplitude should approach δB/B0 ~ 1 and the mean free path collapse toward the gyroradius, pushing the system toward the nonlinear, high-rigidity regime invoked for galactic cosmic rays."],"supporting_citations":[{"why":"Supplies the quasi-linear scattering theory that gives the parallel mean free path λ∥ from the resonant wave spectrum, the core of the self-regulation argument.","marker":"Lee (1983)"},{"why":"Establishes the right-hand resonant ion/ion instability and the hot-beam condition for comparable left-hand growth, the excitation mechanism invoked for the RH and LH families.","marker":"Gary (1985)"},{"why":"Predicts the steady cosmic-ray-pressure-gradient compression used as the baseline forced response that Section 4 rules out for the LP-OB amplitude.","marker":"Drury & Voelk (1981)"},{"why":"Gives the magnetic-pressure (ponderomotive) compression by Alfvénic waves, another candidate source excluded for the compressive family.","marker":"Hollweg (1971)"},{"why":"Supplies the cold-plasma R- and L-mode dispersion branches overlaid on the measured wavenumber–frequency distributions to assign the families.","marker":"Stix (1992)"},{"why":"Calibrates the spacecraft-potential electron density relation n∝e^{−Vsc/Vpe}+C_b, the proxy whose in-phase response underpins the LP-OB fast-mode identification.","marker":"Chen et al. (2013)"},{"why":"The earlier ISEE-3 test of quasi-linear foreshock theory at 1 AU that this study extends by resolving polarization and wavenumber at a near-Sun shock.","marker":"Kennel et al. (1986)"},{"why":"Models the beam-driven shear-Alfvén response used to place the LP-FA family on the resonant Alfvén branch.","marker":"Vainio (2003)"}],"fun_headline_variants":["Shock builds its own scattering foreshock, PSP shows","Proton-driven waves scatter the beam that drives them","Foreshock self-regulation resolved at 0.24 AU","Near-parallel shock: observed self-driven wave loop"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's identification of the oblique compressive family as fast magnetosonic depends on the spacecraft-potential electron density being an in-phase, amplitude-faithful proxy for the true plasma density in the 0.003–0.03 Hz band; if that proxy has a frequency-dependent phase lag or contamination, the near-zero density–field cross-phase that pins the LP-OB family to the fast branch loses its identifying power.","fun_headline_variants_meta":{"raw":{"variants":["Shock builds its own scattering foreshock, PSP shows","Proton-driven waves scatter the beam that drives them","Foreshock self-regulation resolved at 0.24 AU","Near-parallel shock: observed self-driven wave loop"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000369,"raw_usage":{"total_tokens":2037,"prompt_tokens":1060,"completion_tokens":977,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":676,"completion_tokens_details":{"reasoning_tokens":921}},"tokens_in":676,"tokens_out":977,"duration_ms":8884,"temperature":1.0,"reasoning_tokens":921,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:04:53.076946+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the density–field cross-phase in the LP-OB band using the quasi-thermal-noise density estimate from the same spacecraft where it is available at comparable resolution; if the near-zero phase difference between density and magnetic-field magnitude is not reproduced, the fast-magnetosonic identification fails. For the self-regulation loop, run a kinetic simulation of a θ_Bn ≈ 8°, M_A ≈ 7.5 shock with a self-consistently driven foreshock and check that the resonant-band transverse amplitude yields a quasi-linear mean free path within a factor of a few of the measured 0.7–2.9 R_sun and that the RH and LH power stays within a factor of two over the shared wavenumber band.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the right-hand resonant ion/ion instability and the hot-beam condition for comparable left-hand growth, the excitation mechanism invoked for the RH and LH families."},{"cited_title":"F., Coroniti, F","cited_arxiv_id":null,"evidence_quote":"The earlier ISEE-3 test of quasi-linear foreshock theory at 1 AU that this study extends by resolving polarization and wavenumber at a near-Sun shock."}],"review_version":1}