On the competition between radial expansion and Coulomb collisions in shaping the electron velocity distribution function: Kinetic simulations

We present numerical simulations of the solar wind using a fully kinetic model which takes into account the effects of particle's binary collisions in a quasi-neutral plasma in spherical expansion. Starting from an isotropic Maxwellian velocity distribution function for the electrons, we show that the combined effect of expansion and Coulomb collisions leads to the formation of two populations: a collision-dominated cold and dense population almost isotropic in velocity space and a weakly collisional, tenuous field-aligned and antisunward drifting population generated by mirror force focusing in the radially decreasing magnetic field. The relative weights and drift velocities for the two populations observed in our simulations are in excellent agreement with the relative weights and drift velocities for both core and strahl populations observed in the real solar wind. The radial evolution of the main moments of the electron velocity distribution function is in the range observed in the solar wind. The electron temperature anisotropy with respect to the magnetic field direction is found to be related to the ratio between the collisional time and the solar wind expansion time. Even though collisions are found to shape the electron velocity distributions and regulate the properties of the strahl, it is found that the heat flux is conveniently described by a collisionless model where a fraction of the electron thermal energy is advected at the solar wind speed. This reinforces the currently largely admitted fact that collisions in the solar wind are clearly insufficient to force the electron heat flux obey the classical Spitzer-Härm expression where heat flux and temperature gradient are proportional to each other. The presented results show that the electron dynamics in the solar wind cannot be understood without considering the role of collisions.