Non-Maxwellianity of ion velocity distributions in the Earth's magnetosheath

Collisions are nearly negligible in many space and astrophysical plasmas, allowing charged-particle velocity distribution functions (VDFs) to depart from local thermodynamic equilibrium (LTE). How collisionless plasmas relax these non-LTE distributions and convert turbulent energy into particle heating remains an open question. We investigate deviations from LTE in ion velocity distribution functions (iVDFs) within collisionless plasma turbulence using high-resolution measurements from the Magnetospheric Multiscale (MMS) mission. We find that the iVDFs' non-bi-Maxwellian features are widespread and can be significant. Their complexity increases with ion plasma beta and turbulence intensity, with pronounced high-order non-LTE features emerging during intervals of large-amplitude magnetic field fluctuations. In addition, we show that turbulence-induced magnetic curvature plays a significant role in ion scattering and contributes to the isotropization of the iVDF. These results highlight the complex interaction between turbulence and the velocity distribution of charged particles, providing new insight into the kinetic processes responsible for energy conversion in collisionless plasmas.