No Arabic abstract
The evolution of the electron heat flux in the solar wind is regulated by the interplay between several effects: solar wind expansion, that can potentially drive velocity-space instabilties, turbulence and wave-particle interactions, and, possibly, collisions. Here we address the respective role played by the solar wind expansion and the electron firehose instability, developing in the presence of multiple electron populations, in regulating the heat flux. We carry out fully kinetic, Expanding Box Model simulations and separately analyze the enthalpy, bulk and velocity distribution function skewness contributions for each of the electron species. We observe that the key factor determining electron energy flux evolution is the reduction of the drift velocity of the electron populations in the rest frame of the solar wind. In our simulations, redistribution of the electron thermal energy from the parallel to the perpendicular direction after the onset of the electron firehose instability is observed. However, this process seems to impact energy flux evolution only minimally. Hence, reduction of the electron species drift velocity in the solar wind frame appears to directly correlate with efficiency for heat flux instabilities
Using hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvenic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in the linear frequency and residual energy of the Alfv{e}nic fluctuations. In response, these fluctuations modify their interactions and spatial anisotropy to maintain a scale-by-scale critical balance between their characteristic linear and nonlinear frequencies. Once the temperature anisotropy is sufficiently negative, the plasma becomes unstable to kinetic firehose instabilities, which excite rapidly growing magnetic fluctuations at ion-Larmor scales. The consequent pitch-angle scattering of particles maintains the temperature anisotropy near marginal stability, even as the turbulent plasma continues to expand. The resulting evolution of parallel and perpendicular temperatures does not satisfy double-adiabatic conservation laws, but is described accurately by a simple model that includes anomalous scattering. Our results have implications for understanding the complex interplay between macro- and micro-scale physics in various hot, dilute, astrophysical plasmas, and offer predictions concerning power spectra, residual energy, ion-Larmor-scale spectral breaks, and non-Maxwellian features in ion distribution functions that may be tested by measurements taken in high-beta regions of the solar wind.
The role of solar wind expansion in generating whistler waves is investigated using the EB-iPic3D code, which models solar wind expansion self-consistently within a fully kinetic semi-implicit approach. The simulation is initialized with an electron velocity distribution function modeled after Parker Solar Probe observations during its first perihelion at 0.166 au, consisting of a dense core and an anti-sunward strahl. This distribution function is initially stable with respect to kinetic instabilities. Expansion drives the solar wind into successive regimes where whistler heat flux instabilities are triggered. These instabilities produce sunward whistler waves initially characterized by predominantly oblique propagation with respect to the interplanetary magnetic field. The excited waves interact with the electrons via resonant scattering processes. As a consequence, the strahl pitch angle distribution broadens and its drift velocity reduces. Strahl electrons are scattered in the direction perpendicular to the magnetic field, and an electron halo is formed. At a later stage, resonant electron firehose instability is triggered and further affects the electron temperature anisotropy as the solar wind expands. Wave-particle interaction processes are accompanied by a substantial reduction of the solar wind heat flux. The simulated whistler waves are in qualitative agreement with observations in terms of wave frequencies, amplitudes and propagation angles. Our work proposes an explanation for the observations of oblique and parallel whistler waves in the solar wind. We conclude that solar wind expansion has to be factored in when trying to explain kinetic processes at different heliocentric distances.
We present results of two-dimensional fully kinetic Particle-In-Cell simulation in order to shed light on the role of whistler waves in the scattering of strahl electrons and in the heat flux regulation in the solar wind. We model the electron velocity distribution function as initially composed of core and strahl populations as typically encountered in the near-Sun solar wind as observed by Parker Solar Probe. We demonstrate that, as a consequence of the evolution of the electron velocity distribution function, two branches of the whistler heat flux instability can be excited, which can drive whistler waves propagating in the direction oblique or parallel to the background magnetic field. First, oblique whistler waves induce pitch-angle scattering of strahl electrons, towards higher perpendicular velocities. This leads to the broadening of the strahl pitch angle distribution and hence to the formation of a halo-like population at the expense of the strahl. Later on, the electron velocity distribution function experiences the effect of parallel whistler waves, which contributes to the redistribution of the particles scattered in the perpendicular direction into a more symmetric halo, in agreement with observations. Simulation results show a remarkable agreement with the linear theory of the oblique whistler heat flux instability. The process is accompanied by a significant decrease of the heat flux carried by the strahl population.
Understanding the nature of the turbulent fluctuations below the ion gyroradius in solar-wind turbulence is a great challenge. Recent studies have been mostly in favor of kinetic Alfven wave (KAW) type of fluctuations, but other kinds of fluctuations with characteristics typical of magnetosonic, whistler and ion Bernstein modes, could also play a role depending on the plasma parameters. Here we investigate the properties of the sub-proton-scale cascade with high-resolution hybrid-kinetic simulations of freely-decaying turbulence in 3D3V phase space, including electron inertia effects. Two proton plasma beta are explored: the intermediate $beta_p=1$ and low $beta_p=0.2$ regimes, both typically observed in solar wind and corona. The magnetic energy spectum exhibits $k_perp^{-8/3}$ and $k_|^{-7/2}$ power laws at $beta_p=1$, while they are slightly steeper at $beta_p=0.2$. Nevertheless, both regimes develop a spectral anisotropy consistent with $k_|sim k_perp^{2/3}$ at $k_perprho_p>1$, and pronounced small-scale intermittency. In this context, we find that the kinetic-scale cascade is dominated by KAW-like fluctuations at $beta_p=1$, whereas the low-$beta$ case presents a more complex scenario suggesting the simultaneous presence of different types of fluctuations. In both regimes, however, a non-negligible role of ion Bernstein type of fluctuations at the smallest scales seems to emerge.
To properly describe heating in weakly collisional turbulent plasmas such as the solar wind, inter-particle collisions should be taken into account. Collisions can convert ordered energy into heat by means of irreversible relaxation towards the thermal equilibrium. Recently, Pezzi et al. (Phys. Rev. Lett., vol. 116, 2016, p. 145001) showed that the plasma collisionality is enhanced by the presence of fine structures in velocity space. Here, the analysis is extended by directly comparing the effects of the fully nonlinear Landau operator and a linearized Landau operator. By focusing on the relaxation towards the equilibrium of an out of equilibrium distribution function in a homogeneous force-free plasma, here it is pointed out that it is significant to retain nonlinearities in the collisional operator to quantify the importance of collisional effects. Although the presence of several characteristic times associated with the dissipation of different phase space structures is recovered in both the cases of the nonlinear and the linearized operators, the influence of these times is different in the two cases. In the linearized operator case, the recovered characteristic times are systematically larger than in the fully nonlinear operator case, this suggesting that fine velocity structures are dissipated slower if nonlinearities are neglected in the collisional operator.