Formation of Suprathermal Electron Tails in an Expanding, Turbulent Solar Wind: Insights from Fully Kinetic Particle-in-Cell Simulations

As the solar wind propagates through interplanetary space, adiabatic expansion preferentially cools the plasma in the direction perpendicular to the mean magnetic field, while leaving the temperature parallel to the field largely unaffected. The combined effect of the growing temperature anisotropy and the more rapid decrease of magnetic energy relative to the parallel pressure naturally drives the plasma toward the firehose instability threshold. Concurrently, the turbulent cascade from large to small scales leads to kinetic-scale dissipation, resulting in plasma heating and the potential development of suprathermal tails in velocity distribution functions. A central open question is how turbulence-driven heating competes with expansion-induced temperature anisotropies to regulate the onset and nonlinear evolution of kinetic instabilities. In this work, we present the first fully kinetic three-dimensional particle-in-cell (PIC) simulations of an expanding-box system that includes large-scale turbulent forcing, mimicking Alfvénic fluctuations. Our simulations reveal the emergence of suprathermal tails in the electron velocity distribution functions driven by expansion, suggesting an origin in the interplay between turbulence and the firehose instability. This work aims to bridge solar wind observations and theoretical models by providing a unified, fully kinetic framework that captures the coupled effects of expansion, turbulence-driven heating, and kinetic instabilities at electron scales.