pith. sign in

arxiv: 2606.07838 · v1 · pith:ZFXS7SRYnew · submitted 2026-06-05 · ⚛️ physics.plasm-ph · astro-ph.SR· physics.space-ph

Hybrid simulations of the proton beam instabilities in the young solar wind. The formation of hammerhead-like distributions

R. A. L\'opez , Shaaban M. Shaaban , M. Lazar , L. Pezzini , S. Poedts , H. Fichtner , P. H. Yoon This is my paper

Pith reviewed 2026-06-27 20:06 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph astro-ph.SRphysics.space-ph
keywords solar windproton beamshammerhead distributionshybrid simulationsplasma instabilitiesright-handed wavesParker Solar Probequasilinear theory
0
0 comments X

The pith

Proton beam relaxation via right-handed wave instabilities produces hammerhead velocity distributions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents hybrid simulations of proton-beam-plasma systems with parameters typical of the young solar wind. It shows that beam relaxation proceeds through instabilities that generate right-handed polarized electromagnetic waves, and this process imprints hammerhead-type features onto the proton velocity distributions. Feature prominence scales with the magnetic power of the excited waves, which itself traces back to the free energy set by plasma beta and beam drift speed. The simulations follow the self-consistent nonlinear evolution, and their outcomes align with quasilinear theory predictions for the same conditions. A reader would care because the work supplies a concrete dynamical pathway connecting Parker Solar Probe observations of hammerhead distributions and associated waves.

Core claim

From the long-term evolution of these systems, it is found that beam relaxation is driven by instabilities and growing wave fluctuations, leading to HH-type features in the velocity distributions. The production of these features, as well as their prominence, depends on the magnetic power of the waves generated by the instabilities and, therefore, implicitly on the available free energy, quantified by the plasma beta parameter and the relative beam drift. The simulation results capture the self-consistent evolution of the instabilities and their nonlinear development. Linear theory, together with simulations, helps identify the nature of the unstable modes and the plasma conditions under whi

What carries the argument

Hybrid simulations that evolve proton-beam-plasma systems self-consistently, allowing instabilities of right-handed waves to drive beam relaxation and form hammerhead features.

If this is right

  • Hammerhead feature production and strength scale directly with wave magnetic power and therefore with plasma beta and beam drift.
  • Linear theory can identify the unstable modes and the plasma parameter ranges where they grow.
  • Quasilinear theory provides a fast complementary tool for interpreting the wave-particle interactions seen in the simulations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same relaxation pathway may operate in other space-plasma environments that host drifting proton beams and right-handed waves.
  • Varying the simulation beta and drift values across the observed PSP range could map which conditions produce the most prominent hammerhead shapes.
  • Combining hybrid runs with quasilinear calculations offers a practical route to model wave generation across larger volumes of the young solar wind.

Load-bearing premise

The chosen initial beam drifts, densities, and betas are representative of the young solar wind regions where PSP observes hammerhead distributions and right-handed waves.

What would settle it

Time-resolved Parker Solar Probe measurements that either show or fail to show hammerhead features forming specifically during intervals when right-handed wave power grows at the rates and wavelengths predicted by the simulated instability.

Figures

Figures reproduced from arXiv: 2606.07838 by H. Fichtner, L. Pezzini, M. Lazar, P. H. Yoon, R. A. L\'opez, Shaaban M. Shaaban, S. Poedts.

Figure 1
Figure 1. Figure 1: Unstable solutions from linear theory: normal￾ized wave-frequency (left) and growth-rates (right), as func￾tions of k∥ (upper) and both k∥ and k⊥ (lower). The cases correspond to Klein’s parameters in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (top panels) Normalized power spectra for the transverse magnetic field obtained from the FFT of By − iBz at ky = 0. RH waves are shown in ω > 0 and LH waves in ω < 0. Black dashed lines correspond to the linear theory solutions at θ = 0◦ . (bottom panels) Snapshot of the normalized power spectra for the transverse magnetic field. 3.2. RH instabilities Beyond the linear analysis presented in [PITH_FULL_IM… view at source ↗
Figure 4
Figure 4. Figure 4: Moments of the distribution for both species, proton core (top) and beam (bottom). The columns show the drift velocity, plasma beta, and temperature anisotropy. negligible for the core populations (middle panel of the second row). An alternative visualization of these cool￾ing and heating processes is obtained by considering the temperature anisotropy, defined as A = β⊥/β∥, for the beam and core population… view at source ↗
Figure 5
Figure 5. Figure 5: Contours of VDs in the vx-vy plane for various stages in the simulations for Case 1 (left) and Case 2 (right). −4 −2 0 2 4 vx/vA 10−4 10−3 10−2 10−1 100 f(vx) Ωpt = 0.0 Case 1 Total Ωpt = 286.72 Ωpt = 655.36 −5 0 5 vx/vA Ωpt = 0.0 Total Ωpt = 40.96 Case 2 Ωpt = 81.92 [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Parallel cuts of the VD function in the vx, for various stages of both simulations. Top panels for Case 1, bottom panels for Case 2. total (right panels) distributions, illustrating progressive deformations relative to their initial (Ωpt = 0) shapes in response to the enhanced right-handed fluctuations. Apart from the initial condition, we include snapshots during the linear growth phase, and once the satu… view at source ↗
Figure 7
Figure 7. Figure 7: Dynamical path of various simulations with dif￾ferent initial relative drift, beam temperature anisotropy vs. relative drift. All four simulations use the input parameters of [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
read the original abstract

Parker Solar Probe (PSP) observations in the young solar wind reveal new properties of both plasma particle velocity distributions (VDs) and associated electromagnetic (EM) wave fluctuations. The quasilinear (QL) kinetic theory of plasma wave instabilities has recently shown that new hammerhead (HH) proton distributions can be generated by the relaxation of proton beams through the instabilities of right-handed (RH) polarized waves. Such RH waves have indeed been reported in association with HH distributions. In this paper, new results from hybrid simulations of proton-beam-plasma systems with properties typical of those observed to excite EM-RH wave instabilities are presented. From the long-term evolution of these systems, it is found that beam relaxation is driven by instabilities and growing wave fluctuations, leading to HH-type features in the velocity distributions. The production of these features, as well as their prominence, depends on the magnetic power of the waves generated by the instabilities and, therefore, implicitly on the available free energy, quantified by the plasma beta parameter and the relative beam drift. The simulation results capture the self-consistent evolution of the instabilities and their nonlinear development. Linear theory, together with simulations, helps identify the nature of the unstable modes and the plasma conditions under which they arise. The good agreement with quasi-linear (QL) theory further indicates that it can serve as a computationally efficient complementary framework for interpreting the associated wave-particle interactions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents hybrid kinetic-ion simulations of proton beam-plasma systems initialized with parameters typical of the young solar wind. It reports that beam relaxation proceeds through right-handed electromagnetic instabilities, producing hammerhead-like features in the proton velocity distributions whose prominence scales with available free energy (plasma beta and relative beam drift). The results are stated to agree with quasi-linear theory and to capture the self-consistent nonlinear evolution.

Significance. If the central claim holds, the work supplies the first self-consistent numerical demonstration that RH-wave instabilities can generate the hammerhead distributions observed by PSP, thereby linking linear/QL theory to the long-term nonlinear regime and identifying the controlling role of beta and drift. This strengthens the interpretation of PSP wave-particle observations.

major comments (2)
  1. [Introduction / Simulation setup] The abstract and introduction assert that the chosen initial drifts, densities, and betas are 'typical of those observed to excite EM-RH wave instabilities,' yet no table or figure quantitatively compares the simulated parameter set (beta, n_b/n_0, v_d/v_A) against the specific PSP intervals that exhibit both HH distributions and RH waves. Because the formation and prominence of HH features are stated to depend directly on these free-energy parameters, this comparison is load-bearing for the applicability claim.
  2. [Numerical methods / Results] The central result rests on the long-term evolution producing HH features; however, the manuscript does not report resolution studies, particle-per-cell counts, or box-size convergence tests. Without these, it is unclear whether the reported HH morphology is robust against numerical diffusion or finite-particle noise, which directly affects the reliability of the nonlinear saturation stage.
minor comments (2)
  1. [Abstract] Notation for the beam density ratio and drift speed should be defined once at first use and used consistently; the abstract introduces 'relative beam drift' without an accompanying symbol.
  2. [Figure captions] Figure captions should explicitly state the time at which the velocity distributions are shown and whether they are averaged over the simulation domain.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments and positive assessment of the work's significance. We address each major comment below, indicating the revisions planned for the manuscript.

read point-by-point responses

  1. Referee: The abstract and introduction assert that the chosen initial drifts, densities, and betas are 'typical of those observed to excite EM-RH wave instabilities,' yet no table or figure quantitatively compares the simulated parameter set (beta, n_b/n_0, v_d/v_A) against the specific PSP intervals that exhibit both HH distributions and RH waves. Because the formation and prominence of HH features are stated to depend directly on these free-energy parameters, this comparison is load-bearing for the applicability claim.

    Authors: We agree that a quantitative comparison to specific PSP intervals would strengthen the applicability claim. In the revised manuscript we will add a table (or figure) that directly compares the simulation parameters (beta, n_b/n_0, v_d/v_A) to the observed values from the PSP intervals that exhibit both hammerhead distributions and right-handed waves. revision: yes

  2. Referee: The central result rests on the long-term evolution producing HH features; however, the manuscript does not report resolution studies, particle-per-cell counts, or box-size convergence tests. Without these, it is unclear whether the reported HH morphology is robust against numerical diffusion or finite-particle noise, which directly affects the reliability of the nonlinear saturation stage.

    Authors: We acknowledge that explicit convergence tests are necessary to confirm the robustness of the nonlinear saturation stage. In the revised manuscript we will include resolution studies, the adopted particle-per-cell counts, and box-size convergence tests demonstrating that the hammerhead morphology is insensitive to numerical diffusion and particle noise. revision: yes

Circularity Check

0 steps flagged

No significant circularity; simulations yield independent results

full rationale

The paper's core claim—that long-term hybrid simulations of proton-beam systems with observationally motivated initial drifts, densities, and betas produce HH-type velocity distribution features via RH-wave instabilities—is an output of the self-consistent particle-in-cell evolution, not a re-expression of the inputs. Linear and quasi-linear theory are invoked only for mode identification and as a complementary check, with the abstract explicitly stating that simulations capture the nonlinear development independently. No equations reduce the reported HH formation to a fitted parameter or prior self-citation by construction; the agreement with QL theory is presented as validation rather than the source of the result. Parameter selection is justified by external PSP observations, not by internal redefinition, satisfying the criteria for a self-contained derivation.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

Review performed on abstract only; free parameters are the input plasma beta and beam drift speed chosen to match observations; no new entities postulated.

free parameters (2)
  • plasma beta
    Input parameter varied to control available free energy; values chosen to represent young solar wind conditions.
  • relative beam drift
    Input parameter that sets the free energy available to drive the instability.
axioms (1)
  • domain assumption Hybrid approximation (fluid electrons, kinetic ions) captures the relevant wave-particle dynamics for right-handed electromagnetic modes.
    Implicit in the choice of hybrid simulation method; stated in the context of comparing to quasi-linear kinetic theory.

pith-pipeline@v0.9.1-grok · 5823 in / 1318 out tokens · 19366 ms · 2026-06-27T20:06:49.324748+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    A., Marsch, E., & F

    Araneda, J. A., Marsch, E., & F. -Vi˜ nas, A. 2008, PhRvL, 100, 125003, doi: 10.1103/PhysRevLett.100.125003 Bharati Das, S., Verniero, J. L., Badman, S. T., et al. 2026, The Astrophysical Journal Letters, 1001, L7, doi: 10.3847/2041-8213/ae52f6

    work page doi:10.1103/physrevlett.100.125003 2008

  2. [2]

    J., Albright, B

    Bowers, K. J., Albright, B. J., Yin, L., Bergen, B., & Kwan, T. J. T. 2008, Physics of Plasmas, 15, 055703, doi: 10.1063/1.2840133

    work page doi:10.1063/1.2840133 2008

  3. [3]

    2016, title Large Gradual Solar Energetic Particle Events , Living Reviews in Solar Physics, 13, 3, 10.1007/s41116-016-0002-5

    Desai, M., & Giacalone, J. 2016, Living Reviews in Solar Physics, 13, 3, doi: 10.1007/s41116-016-0002-5

    work page doi:10.1007/s41116-016-0002-5 2016

  4. [4]

    2026, The Astrophysical Journal, 997, 174, doi: 10.3847/1538-4357/ae2601

    Franci, L. 2026, The Astrophysical Journal, 997, 174, doi: 10.3847/1538-4357/ae2601

    work page doi:10.3847/1538-4357/ae2601 2026

  5. [5]

    P., & Lee, M

    Gary, S. P., & Lee, M. A. 1994, Journal of Geophysical Research: Space Physics, 99, 11297, doi: https://doi.org/10.1029/94JA00253

    work page doi:10.1029/94ja00253 1994

  6. [6]

    G., Verniero, J

    Klein, K. G., Verniero, J. L., Alterman, B., et al. 2021, ApJ, 909, 7, doi: 10.3847/1538-4357/abd7a0

    work page doi:10.3847/1538-4357/abd7a0 2021

  7. [7]

    B., et al

    Lario, D., Berger, L., Decker, R. B., et al. 2019, ApJ, 158, 12, doi: 10.3847/1538-3881/ab1e49

    work page doi:10.3847/1538-3881/ab1e49 2019

  8. [8]

    G., Wilson III, L

    Lario, D., Richardson, I. G., Wilson III, L. B., et al. 2022, ApJ, 925, 198, doi: 10.3847/1538-4357/ac3c47

    work page doi:10.3847/1538-4357/ac3c47 2022

  9. [9]

    2021, A&A, 656, A37, doi: 10.1051/0004-6361/202141149

    Lavraud, B., Kieokaew, R., Fargette, N., et al. 2021, A&A, 656, A37, doi: 10.1051/0004-6361/202141149

    work page doi:10.1051/0004-6361/202141149 2021

  10. [10]

    2023, Physics of Plasmas, 30, 063902, doi: 10.1063/5.0146529

    Le, A., Stanier, A., Yin, L., et al. 2023, Physics of Plasmas, 30, 063902, doi: 10.1063/5.0146529

    work page doi:10.1063/5.0146529 2023

  11. [11]

    2021, Journal of Geophysical Research: Space Physics, 126, 1, doi: 10.1029/2021JA029125 L´ opez, R

    Le, A., Winske, D., Stanier, A., et al. 2021, Journal of Geophysical Research: Space Physics, 126, 1, doi: 10.1029/2021JA029125 L´ opez, R. A., Shaaban, S. M., & Lazar, M. 2021, J. Plasma Phys., 87, 905870310, doi: 10.1017/S0022377821000593

    work page doi:10.1029/2021ja029125 2021

  12. [12]

    2021, A&A, 656, A36, doi: 10.1051/0004-6361/202141095 L´ opez, R

    Louarn, P., Fedorov, A., Prech, L., et al. 2021, A&A, 656, A36, doi: 10.1051/0004-6361/202141095 L´ opez, R. A. 2023, DIS-K: Dispersion Solver for Kappa Plasmas, 1.0, Zenodo, doi: 10.5281/zenodo.15991

    work page doi:10.1051/0004-6361/202141095 2021

  13. [13]

    Matteini, L., Landi, S., Velli, M., & Matthaeus, W. H. 2013, The Astrophysical Journal, 763, 142, doi: 10.1088/0004-637X/763/2/142

    work page doi:10.1088/0004-637x/763/2/142 2013

  14. [14]

    2022, ApJ, 926, 185, doi: 10.3847/1538-4357/ac402c

    Larson, D. 2022, ApJ, 926, 185, doi: 10.3847/1538-4357/ac402c

    work page doi:10.3847/1538-4357/ac402c 2022

  15. [15]

    N., Arr´ o, G., & L´ opez, R

    Pezzini, L., Bacchini, F., Zhukov, A. N., Arr´ o, G., & L´ opez, R. A. 2026, The Astrophysical Journal, 997, 158, doi: 10.3847/1538-4357/ae27a3

    work page doi:10.3847/1538-4357/ae27a3 2026

  16. [16]

    N., Bacchini, F., et al

    Pezzini, L., Zhukov, A. N., Bacchini, F., et al. 2024, ApJ, 975, 37, doi: 10.3847/1538-4357/ad7465

    work page doi:10.3847/1538-4357/ad7465 2024

  17. [17]

    D., Verniero, J

    Phan, T. D., Verniero, J. L., Larson, D., et al. 2022, Geophys. Res. Lett., 49, e96986, doi: 10.1029/2021GL096986

    work page doi:10.1029/2021gl096986 2022

  18. [18]

    2010, AIP Conf

    Pierrard, V., & Voitenko, Y. 2010, AIP Conf. Proc., 1216, 102, doi: 10.1063/1.3395812

    work page doi:10.1063/1.3395812 2010

  19. [19]

    H., Lazar, M., & Qureshi, M

    Saeed, S., Sarfraz, M., Yoon, P. H., Lazar, M., & Qureshi, M. N. S. 2017, Monthly Notices of the Royal Astronomical Society, 465, 1672, doi: 10.1093/mnras/stw2900

    work page doi:10.1093/mnras/stw2900 2017

  20. [20]

    2025, The Astrophysical Journal, 993, 1, doi: 10.3847/1538-4357/ae09ac

    Poedts, S. 2025, The Astrophysical Journal, 993, 1, doi: 10.3847/1538-4357/ae09ac

    work page doi:10.3847/1538-4357/ae09ac 2025

  21. [21]

    2024, A&A, 692, L6, doi: 10.1051/0004-6361/202452205

    Poedts, S. 2024, A&A, 692, L6, doi: 10.1051/0004-6361/202452205

    work page doi:10.1051/0004-6361/202452205 2024

  22. [22]

    M., Lazar, M., & Poedts, S

    Shaaban, S. M., Lazar, M., & Poedts, S. 2018, Monthly Notices of the Royal Astronomical Society, 480, 310, doi: 10.1093/mnras/sty1567

    work page doi:10.1093/mnras/sty1567 2018

  23. [23]

    2002, Journal of Geophysical Research: Space Physics, 107, SSH 8, doi: https://doi.org/10.1029/2002JA009264

    Tu, C.-Y., Wang, L.-H., & Marsch, E. 2002, Journal of Geophysical Research: Space Physics, 107, SSH 8, doi: https://doi.org/10.1029/2002JA009264

    work page doi:10.1029/2002ja009264 2002

  24. [24]

    L., Larson, D

    Verniero, J. L., Larson, D. E., Livi, R., et al. 2020, ApJS, 248, 5, doi: 10.3847/1538-4365/ab86af

    work page doi:10.3847/1538-4365/ab86af 2020

  25. [25]

    L., Chandran, B

    Verniero, J. L., Chandran, B. D. G., Larson, D. E., et al. 2022, ApJ, 924, 112, doi: 10.3847/1538-4357/ac36d5

    work page doi:10.3847/1538-4357/ac36d5 2022

  26. [26]

    2015, SoPh, 290, 1231, doi: 10.1007/s11207-015-0661-8

    Voitenko, Y., & Pierrard, V. 2015, SoPh, 290, 1231, doi: 10.1007/s11207-015-0661-8

    work page doi:10.1007/s11207-015-0661-8 2015