pith. machine review for the scientific record. sign in

arxiv: 2605.01895 · v1 · submitted 2026-05-03 · 🌌 astro-ph.SR · physics.plasm-ph

Recognition: unknown

Formation of Suprathermal Electron Populations in the Expanding, Turbulent Solar Wind

Maximilien P\'eters de Bonhome , Fabio Bacchini , Luca Pezzini , Viviane Pierrard
Authors on Pith no claims yet

Pith reviewed 2026-05-09 16:35 UTC · model grok-4.3

classification 🌌 astro-ph.SR physics.plasm-ph
keywords solar windsuprathermal electronskinetic simulationparticle-in-cellturbulenceexpansionfirehose instabilitynonthermal distributions
0
0 comments X

The pith

Solar wind simulations reveal how expansion and turbulence generate suprathermal electron tails.

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

The paper deploys the first fully kinetic particle-in-cell simulation of an expanding, turbulent plasma matching heliospheric conditions to trace the origin of nonthermal electron features routinely seen in solar wind observations. Expansion weakens the magnetic field and cools the plasma perpendicular to it, pushing the system toward the firehose instability threshold, while anisotropic turbulence simultaneously heats perpendicularly and seeds nonthermal distributions. Suprathermal tails then form preferentially along the parallel direction and survive even after the instability begins to regulate the anisotropy. The work supplies a single self-consistent kinetic picture connecting expansion, turbulence, and instabilities to the observed nonthermal populations that influence energy transport throughout the heliosphere.

Core claim

In the first fully kinetic particle-in-cell simulation of an expanding turbulent plasma under heliospheric conditions, expansion-driven weakening of the magnetic field adiabatically cools the plasma perpendicularly to the mean field while leaving the parallel temperature largely unchanged, driving the system toward the firehose-instability threshold; concurrently, strongly anisotropic turbulence leads to perpendicular heating and the development of nonthermal features, after which suprathermal electron populations preferentially develop in the parallel direction, forming pronounced power-law tails even under weakly compressive, highly Alfvénic conditions, and persist despite anisotropy being

What carries the argument

The coupled kinetic evolution of expansion-induced adiabatic perpendicular cooling, turbulence-driven perpendicular heating, and firehose-instability regulation that channels energy into parallel suprathermal electron tails.

Load-bearing premise

The specific expansion rate, turbulence amplitude, and initial conditions chosen in the simulation accurately capture real solar-wind behavior without important missing physics such as additional wave modes or boundary effects.

What would settle it

Comparison of the simulated parallel and perpendicular electron velocity distributions against in-situ data from spacecraft traversing regions of measured solar-wind expansion and Alfvénic turbulence levels, checking whether parallel power-law tails appear at the predicted amplitudes.

Figures

Figures reproduced from arXiv: 2605.01895 by Fabio Bacchini, Luca Pezzini, Maximilien P\'eters de Bonhome, Viviane Pierrard.

Figure 1
Figure 1. Figure 1: Volume rendering of the z-component of the mag￾netic-field fluctuations, δBz/Bg, shown at the initial expan￾sion time (t = 0.5τexp), mid-expansion time (t = 1.0τexp), and final expansion time (t = 1.5τexp). et al. 2007). This results in an ion thermal speed vthi0 = p 3kBTi0/mi ≈ 0.026c. Alfv´enic turbulence in our simulation is obtained through the aforementioned charge-independent forc￾ing applied on part… view at source ↗
Figure 2
Figure 2. Figure 2: (a) Evolution of the box-averaged magnetic-field strength normalized by the initial background magnetic field. (b) Evolution of the root-mean-square (rms) speed of electrons (in blue), ions (in red), and the magnetic-field fluctuations normalized by the background magnetic field in black. (c) Evolution of the average temperature ratio, T⊥/T∥, for electrons (in blue) and ions (in red), with the associated f… view at source ↗
Figure 3
Figure 3. Figure 3: (a) Distribution of (T⊥e/T∥e, β∥e) for electrons at t = 0.5τexp (left), t = 1.0τexp (middle), and t = 1.5τexp (right). The distribution is scaled as the amount of cells (nc) normalized by the total amount of cells (Nc). Isocontours of the theo￾retical oblique EFI are plotted as blue lines with T⊥e/T∥e = 1 − 1.29/β0.97 ∥e (dotted), T⊥e/T∥e = 1 − 1.32/β0.61 ∥e (dashed), and T⊥e/T∥e = 1 − 1.36/β0.47 ∥e (solid… view at source ↗
Figure 4
Figure 4. Figure 4: Time evolution of the perpendicular (fe(v⊥e)) and parallel (fe(v∥e)) electron VDFs. The green line corresponds to a Maxwellian fit at t = 0, while the dashed black line corresponds to a Kappa distribution fit (e.g., C. S. Salem et al. 2023) at t = 1.5τexp with parameter κ = 6.2 for fe(v⊥e) and κ = 5 for fe(v∥e). The solid blue line indicates the electron VDF at t = 0.5τexp (onset time of expansion). The in… view at source ↗
read the original abstract

Nonthermal features are ubiquitously observed in electron velocity distribution functions in the solar wind, yet their origin in the collisionless, turbulent, expanding solar-wind plasma remains unclear. We investigate how solar-wind expansion and Alfv\'enic turbulence jointly generate and regulate these features using the first fully kinetic particle-in-cell simulation of an expanding turbulent plasma under heliospheric conditions. In our setup, expansion-driven weakening of the magnetic field adiabatically cools the plasma perpendicularly to the mean field while leaving the parallel temperature largely unchanged, driving the system toward the firehose-instability threshold. Concurrently, strongly anisotropic turbulence leads to perpendicular heating and the development of nonthermal features. Subsequently, we find that suprathermal electron populations preferentially develop in the parallel direction, forming pronounced power-law tails even under weakly compressive, highly Alfv\'enic conditions, and persist despite anisotropy regulation by the firehose instability. The preferentially parallel energization suggests the involvement of parallel electric fields or resonant wave--particle interactions, rather than simple velocity-space redistribution. These results provide the first direct evidence of the emergence of nonthermal-electron features in a unified kinetic framework linking expansion, turbulence, and instabilities in the solar wind.

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 / 1 minor

Summary. The paper reports results from the first fully kinetic particle-in-cell simulation of an expanding, Alfvénic turbulent plasma under heliospheric conditions. It claims that expansion-driven perpendicular cooling combined with anisotropic turbulence drives the plasma toward the firehose threshold while generating suprathermal electron populations that form pronounced parallel power-law tails; these tails persist even after firehose regulation of anisotropy, implying involvement of parallel electric fields or resonant interactions rather than simple redistribution. The work positions this as direct evidence for a unified kinetic mechanism linking expansion, turbulence, and instabilities in the solar wind.

Significance. If the simulation results prove robust across parameter space, the work would be significant for solar-wind physics: it provides a self-consistent kinetic demonstration that nonthermal electron features can emerge from the interplay of expansion cooling, Alfvénic turbulence, and firehose regulation without invoking separate mechanisms. This addresses an open question about the origin of observed suprathermal electrons. The use of a single unified PIC framework is a strength, though the absence of quantitative measures (e.g., exact spectral indices, error bars) and sensitivity tests in the presented description reduces immediate applicability to observations.

major comments (2)
  1. [Abstract] Abstract: The central claim that parallel power-law tails arise from the joint action of expansion, turbulence, and firehose regulation is load-bearing for the 'unified framework' attribution, yet the description is limited to a single run under unspecified 'heliospheric conditions' with no mention of control simulations at varied expansion rates or turbulence amplitudes. If the tails weaken or vanish under modest changes to these parameters, the generality of the mechanism is not established.
  2. [Abstract] Abstract: The manuscript asserts 'first direct evidence' and 'pronounced power-law tails' but supplies no quantitative details such as measured spectral indices, their uncertainties, or comparison to observed solar-wind statistics. Without these, it is difficult to assess whether the simulated tails match heliospheric data or are an artifact of the chosen initial beta, anisotropy, or resolution.
minor comments (1)
  1. [Abstract] The abstract would benefit from a brief statement of the simulation domain size, particle number, and time-stepping scheme to allow readers to gauge numerical convergence.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their detailed and insightful report. Their comments have helped us clarify the scope and quantitative aspects of our simulation results. We address each major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim that parallel power-law tails arise from the joint action of expansion, turbulence, and firehose regulation is load-bearing for the 'unified framework' attribution, yet the description is limited to a single run under unspecified 'heliospheric conditions' with no mention of control simulations at varied expansion rates or turbulence amplitudes. If the tails weaken or vanish under modest changes to these parameters, the generality of the mechanism is not established.

    Authors: We recognize the importance of establishing the robustness of the observed power-law tails across parameter space. Our study focuses on a single, computationally demanding simulation chosen to represent typical heliospheric conditions (as detailed in Section 2 of the manuscript). Performing additional control simulations would require substantial additional resources and is beyond the scope of the current work. However, we have revised the abstract to remove the implication of broad generality and instead highlight the results for the simulated conditions. We have also added a paragraph in the discussion section addressing how the mechanism might depend on expansion rate and turbulence amplitude, drawing on supporting evidence from reduced models and linear analysis. revision: partial

  2. Referee: [Abstract] Abstract: The manuscript asserts 'first direct evidence' and 'pronounced power-law tails' but supplies no quantitative details such as measured spectral indices, their uncertainties, or comparison to observed solar-wind statistics. Without these, it is difficult to assess whether the simulated tails match heliospheric data or are an artifact of the chosen initial beta, anisotropy, or resolution.

    Authors: We agree that including quantitative details strengthens the manuscript. In the revised version, we have added measurements of the spectral indices of the parallel suprathermal tails, along with estimates of their uncertainties from the fitting procedure. We have also added a comparison to typical observed spectral indices in the solar wind from spacecraft data. These additions are now in the results section and referenced in the abstract. The 'first direct evidence' phrasing has been revised to 'provides direct evidence' to accurately reflect the scope of our single simulation study. revision: yes

standing simulated objections not resolved
  • The lack of additional control simulations to demonstrate robustness across parameter variations, which we cannot perform due to computational limitations.

Circularity Check

0 steps flagged

No circularity: results emerge from explicit numerical dynamics

full rationale

The paper reports outcomes of a single fully kinetic PIC simulation with specified expansion rate, turbulence amplitude, and initial conditions. Nonthermal parallel tails and firehose regulation are observed consequences of the time-stepped Vlasov-Maxwell evolution under those inputs; no equation is rearranged to equal its own fitted parameter, no prediction is statistically forced by a prior fit, and no uniqueness theorem or ansatz is imported via self-citation to close the argument. The framework is self-contained as a forward simulation whose outputs are not definitionally identical to the chosen setup.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the numerical model and initial conditions rather than new physical axioms or entities. No new particles or forces are invented; it uses standard plasma physics.

free parameters (2)
  • expansion rate
    The rate at which the plasma expands, chosen to mimic solar wind but value not given in abstract.
  • turbulence amplitude
    Strength of the Alfvénic turbulence imposed in the simulation.
axioms (2)
  • domain assumption The solar wind plasma is collisionless on the scales of interest
    Standard assumption for kinetic simulations of solar wind, invoked implicitly in the PIC setup.
  • domain assumption The turbulence is highly Alfvénic and weakly compressive
    Stated in the abstract as the condition under which the features develop.

pith-pipeline@v0.9.0 · 5527 in / 1311 out tokens · 49579 ms · 2026-05-09T16:35:54.296430+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

49 extracted references · 47 canonical work pages

  1. [1]

    JGR , author =

    A kinetic exospheric model of the solar wind with a nonmonotonic potential energy for the protons , volume =. JGR , author =. 2003 , pages =. doi:10.1029/2002JA009487 , language =

    work page doi:10.1029/2002ja009487 2003

  2. [2]

    A&A , author =

    A kinetic model of the solar wind with. A&A , author =. 1997 , note =

    1997

  3. [3]

    ApJ , author =

    A. ApJ , author =. 2004 , pages =. doi:10.1086/382866 , language =

    work page doi:10.1086/382866 2004

  4. [4]

    A&A , author =

    A kinetic model of solar wind acceleration driven by ambipolar electric potential and velocity-space diffusion , volume =. A&A , author =. 2025 , pages =. doi:10.1051/0004-6361/202554250 , abstract =

    work page doi:10.1051/0004-6361/202554250 2025

  5. [5]

    ApJL , author =

    Adaptive. ApJL , author =. 2021 , pages =. doi:10.3847/2041-8213/ac37c2 , abstract =

    work page doi:10.3847/2041-8213/ac37c2 2021

  6. [6]

    ApJ , author =

    Hybrid-kinetic. ApJ , author =. 2019 , pages =. doi:10.3847/1538-4357/ab20cc , abstract =

    work page doi:10.3847/1538-4357/ab20cc 2019

  7. [7]

    ApJ , author =

    On the. ApJ , author =. 2021 , pages =. doi:10.3847/1538-4357/ac1067 , abstract =

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

  8. [8]

    ApJ , author =

    A. ApJ , author =. 2019 , pages =. doi:10.3847/1538-4357/aaf1be , abstract =

    work page doi:10.3847/1538-4357/aaf1be 2019

  9. [9]

    2015 , pages =

    ApJ , author =. 2015 , pages =. doi:10.1088/0004-637X/800/2/88 , number =

    work page doi:10.1088/0004-637x/800/2/88 2015

  10. [10]

    JPlPh , author =

    Proton fire hose instabilities in the expanding solar wind , volume =. JPlPh , author =. 2017 , pages =. doi:10.1017/S0022377817000071 , abstract =

    work page doi:10.1017/s0022377817000071 2017

  11. [11]

    FrASS , author =

    Electron-. FrASS , author =. 2022 , pages =. doi:10.3389/fspas.2022.951628 , abstract =

    work page doi:10.3389/fspas.2022.951628 2022

  12. [12]

    R., Uzdensky, D

    ApJ , author =. 2013 , pages =. doi:10.1088/0004-637X/770/2/147 , number =

    work page doi:10.1088/0004-637x/770/2/147 2013

  13. [13]

    ApJ , author =

    Fully. ApJ , author =. 2022 , pages =. doi:10.3847/1538-4357/ac8a94 , abstract =

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

  14. [14]

    ApJ , author =

    Alpha. ApJ , author =. 2026 , pages =. doi:10.3847/1538-4357/ae27a3 , abstract =

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

  15. [15]

    and Péters de Bonhome, Maximilien and Els, Paul and Argyropoulos, Konstantinos-Xanthos and Ly, Minh Nhat and Grošelj, Daniel , year =

    Bacchini, Fabio and Gorbunov, Evgeny A. and Péters de Bonhome, Maximilien and Els, Paul and Argyropoulos, Konstantinos-Xanthos and Ly, Minh Nhat and Grošelj, Daniel , year =. Particle-in-. doi:10.48550/ARXIV.2602.15939 , abstract =

    work page doi:10.48550/arxiv.2602.15939

  16. [16]

    PhPl , author =

    Resonant electron firehose instability:. PhPl , author =. 2003 , pages =. doi:10.1063/1.1590982 , abstract =

    work page doi:10.1063/1.1590982 2003

  17. [17]

    A&A , author =

    Precision electron measurements in the solar wind at 1 au from. A&A , author =. 2023 , pages =. doi:10.1051/0004-6361/202141816 , abstract =

    work page doi:10.1051/0004-6361/202141816 2023

  18. [18]

    ApJL , author =

    Particle-in-cell. ApJL , author =. 2019 , pages =. doi:10.3847/2041-8213/ab0c95 , abstract =

    work page doi:10.3847/2041-8213/ab0c95 2019

  19. [19]

    JGRA , author =

    A model of turbulence in magnetized plasmas:. JGRA , author =. 2008 , pages =. doi:10.1029/2007JA012665 , abstract =

    work page doi:10.1029/2007ja012665 2008

  20. [20]

    Servidio, F

    Local. PhRvL , author =. 2012 , pages =. doi:10.1103/PhysRevLett.108.045001 , language =

    work page doi:10.1103/physrevlett.108.045001 2012

  21. [21]

    Nature Communications , author =

    Evidence for electron. Nature Communications , author =. 2019 , pages =. doi:10.1038/s41467-019-08435-3 , abstract =

    work page doi:10.1038/s41467-019-08435-3 2019

  22. [22]

    LRSP , author =

    Kinetic. LRSP , author =. 2006 , file =. doi:10.12942/lrsp-2006-1 , language =

    work page doi:10.12942/lrsp-2006-1 2006

  23. [23]

    A&A , author =

    Characteristics of solar wind suprathermal halo electrons , volume =. A&A , author =. 2020 , pages =. doi:10.1051/0004-6361/202038830 , abstract =

    work page doi:10.1051/0004-6361/202038830 2020

  24. [24]

    ApJ , author =

    Coronal. ApJ , author =. 2020 , pages =. doi:10.3847/1538-4357/ab7b7a , abstract =

    work page doi:10.3847/1538-4357/ab7b7a 2020

  25. [25]

    ApJS , author =

    Electrons in the. ApJS , author =. 2020 , pages =. doi:10.3847/1538-4365/ab4cec , abstract =

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

  26. [26]

    ApJ , author =

    Radial. ApJ , author =. 2022 , pages =. doi:10.3847/1538-4357/ac6605 , abstract =

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

  27. [27]

    ApJ , author =

    Radial. ApJ , author =. 2024 , pages =. doi:10.3847/1538-4357/ad7d05 , abstract =

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

  28. [28]

    Kappa distributions: theory and applications in plasmas , isbn =

    Livadiotis, George , year =. Kappa distributions: theory and applications in plasmas , isbn =

  29. [29]

    & Lazar, M

    Kappa. SoPh , author =. 2010 , pages =. doi:10.1007/s11207-010-9640-2 , language =

    work page doi:10.1007/s11207-010-9640-2 2010

  30. [30]

    ApJ , author =

    Self-induced. ApJ , author =. 2019 , pages =. doi:10.3847/1538-4357/ab4c30 , abstract =

    work page doi:10.3847/1538-4357/ab4c30 2019

  31. [31]

    ApJL , author =

    Collisionless. ApJL , author =. 2020 , pages =. doi:10.3847/2041-8213/aba591 , abstract =

    work page doi:10.3847/2041-8213/aba591 2020

  32. [32]

    ApJ , author =

    The. ApJ , author =. 2024 , pages =. doi:10.3847/1538-4357/ad1329 , abstract =

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

  33. [33]

    F., Goldberger, M

    The. RSPSA , author =. 1956 , pages =. doi:10.1098/rspa.1956.0116 , abstract =

    work page doi:10.1098/rspa.1956.0116 1956

  34. [34]

    JGR , author =

    Kinetic models of the solar wind , volume =. JGR , author =. 1971 , pages =. doi:10.1029/JA076i031p07479 , language =

    work page doi:10.1029/ja076i031p07479 1971

  35. [35]

    ApJ , author =

    The. ApJ , author =. 2022 , pages =. doi:10.3847/1538-4357/ac85b8 , abstract =

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

  36. [36]

    JGRA , author =

    Electron kappa distribution and quasi‐thermal noise , volume =. JGRA , author =. 2014 , pages =. doi:10.1002/2014JA020353 , abstract =

    work page doi:10.1002/2014ja020353 2014

  37. [37]

    2018, Phys

    Particle. PhRvL , author =. 2018 , pages =. doi:10.1103/PhysRevLett.121.255101 , language =

    work page doi:10.1103/physrevlett.121.255101 2018

  38. [38]

    ApJ , author =

    Turbulence versus. ApJ , author =. 2019 , pages =. doi:10.3847/1538-4357/ab3e01 , abstract =

    work page doi:10.3847/1538-4357/ab3e01 2019

  39. [39]

    E&SS , author =

    Intermittency in. E&SS , author =. 2019 , pages =. doi:10.1029/2018EA000535 , abstract =

    work page doi:10.1029/2018ea000535 2019

  40. [40]

    MNRAS , author =

    Dependence of kinetic plasma waves on ion-to-electron mass ratio and light-to-. MNRAS , author =. 2020 , pages =. doi:10.1093/mnras/staa977 , abstract =

    work page doi:10.1093/mnras/staa977 2020

  41. [41]

    JGR , author =

    Electron parameter correlations in high-speed streams and heat flux instabilities , volume =. JGR , author =. 1976 , pages =. doi:10.1029/JA081i013p02377 , language =

    work page doi:10.1029/ja081i013p02377 1976

  42. [42]

    JGRA , author =

    Regulation of the solar wind electron heat flux from 1 to 5. JGRA , author =. 1994 , pages =. doi:10.1029/94JA02068 , abstract =

    work page doi:10.1029/94ja02068 1994

  43. [43]

    Journal of Geophysical Research , author =

    Solar wind electrons , volume =. JGR , author =. 1975 , pages =. doi:10.1029/JA080i031p04181 , language =

    work page doi:10.1029/ja080i031p04181 1975

  44. [44]

    MNRAS , author =

    Firehose instabilities triggered by the solar wind suprathermal electrons , volume =. MNRAS , author =. 2019 , pages =. doi:10.1093/mnras/sty3377 , language =

    work page doi:10.1093/mnras/sty3377 2019

  45. [45]

    ApJ , author =

    Quantifying the. ApJ , author =. 2023 , pages =. doi:10.3847/1538-4357/acd769 , abstract =

    work page doi:10.3847/1538-4357/acd769 2023

  46. [46]

    ApJ , author =

    Mirror. ApJ , author =. 2017 , pages =. doi:10.3847/1538-4357/aa67e0 , abstract =

    work page doi:10.3847/1538-4357/aa67e0 2017

  47. [47]

    JGRA , author =

    Oblique proton fire hose instability in the expanding solar wind:. JGRA , author =. 2008 , pages =. doi:10.1029/2008JA013416 , abstract =

    work page doi:10.1029/2008ja013416 2008

  48. [48]

    JGRA , author =

    Parallel proton fire hose instability in the expanding solar wind:. JGRA , author =. 2006 , pages =. doi:10.1029/2006JA011667 , abstract =

    work page doi:10.1029/2006ja011667 2006

  49. [49]

    GeoRL , author =

    Evolution of the solar wind proton temperature anisotropy from 0.3 to 2.5. GeoRL , author =. 2007 , pages =. doi:10.1029/2007GL030920 , abstract =

    work page doi:10.1029/2007gl030920 2007