No Arabic abstract
We present waveform observations of electromagnetic lower hybrid and whistler waves with f_ci << f < f_ce downstream of four supercritical interplanetary (IP) shocks using the Wind search coil magnetometer. The whistler waves were observed to have a weak positive correlation between partialB and normalized heat flux magnitude and an inverse correlation with T_eh/T_ec. All were observed simultaneous with electron distributions satisfying the whistler heat flux instability threshold and most with T_{perp,h}/T_{para,h} > 1.01. Thus, the whistler mode waves appear to be driven by a heat flux instability and cause perpendicular heating of the halo electrons. The lower hybrid waves show a much weaker correlation between partialB and normalized heat flux magnitude and are often observed near magnetic field gradients. A third type of event shows fluctuations consistent with a mixture of both lower hybrid and whistler mode waves. These results suggest that whistler waves may indeed be regulating the electron heat flux and the halo temperature anisotropy, which is important for theories and simulations of electron distribution evolution from the sun to the earth.
A statistical analysis of 15,210 electron velocity distribution function (VDF) fits, observed within $pm$2 hours of 52 interplanetary (IP) shocks by the $Wind$ spacecraft near 1 AU, is presented. This is the second in a three-part series on electron VDFs near IP shocks. The electron velocity moment statistics for the dense, low energy core, tenuous, hot halo, and field-aligned beam/strahl are a statistically significant list of values illustrated with both histograms and tabular lists for reference and baselines in future work. The beam/strahl fit results in the upstream are currently the closest thing to a proper parameterization of the beam/strahl electron velocity moments in the ambient solar wind. This work will also serve as a 1 AU baseline and reference for missions like $Parker Solar Probe$ and $Solar Orbiter$. The median density, temperature, beta, and temperature anisotropy values for the core(halo)[beam/strahl] components, with subscripts $ec$($eh$)[$eb$], of all fit results respectively are $n{scriptstyle_{ec(h)[b]}}$ $sim$ 11.3(0.36)[0.17] $cm^{-3}$, $T{scriptstyle_{ec(h)[b], tot}}$ $sim$ 14.6(48.4)[40.2] $eV$, $beta{scriptstyle_{ec(h)[b], tot}}$ $sim$ 0.93(0.11)[0.05], and $mathcal{A}{scriptstyle_{ec(h)[b]}}$ $sim$ 0.98(1.03)[0.93]. The nuanced details of the fitting method and data product description were published in Paper I and the detailed analysis of the results will be shown in Paper III.
Analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within $pm$2 hours of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 AU, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partial-to-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions) while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature, even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.
Analysis of 15314 electron velocity distribution functions (VDFs) within $pm$2 hours of 52 interplanetary (IP) shocks observed by the emph{Wind} spacecraft near 1 AU are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar velocity distribution function, while both the halo and beam/strahl components were best fit to bi-kappa velocity distribution function. This is the first statistical study to show that the core electron distribution is better fit to a self-similar velocity distribution function than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The range of values defined by the lower and upper quartiles for the kappa exponents are $kappa{scriptstyle_{ec}}$ $sim$ 5.40--10.2 for the core, $kappa{scriptstyle_{eh}}$ $sim$ 3.58--5.34 for the halo, and $kappa{scriptstyle_{eb}}$ $sim$ 3.40--5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents are $s{scriptstyle_{ec}}$ $sim$ 2.00--2.04, and asymmetric bi-self-similar core exponents are $p{scriptstyle_{ec}}$ $sim$ 2.20--4.00 for the parallel exponent, and $q{scriptstyle_{ec}}$ $sim$ 2.00--2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.
Shock accelerated electrons are found in many astrophysical environments, and the mechanisms by which they are accelerated to high energies are still not completely clear. For relatively high Mach numbers, the shock is supercritical, and its front exhibit broadband fluctuations, or ripples. Shock surface fluctuations have been object of many observational and theoretical studies, and are known to be important for electron acceleration. We employ a combination of hybrid Particle-In-Cell and test-particle methods to study how shock surface fluctuations influence the acceleration of suprathermal electrons in fully three dimensional simulations, and we give a complete comparison for the 2D and 3D cases. A range of different quasi-perpendicular shocks in 2D and 3D is examined, over a range of parameters compatible with the ones observed in the solar wind. Initial electron velocity distributions are taken as kappa functions, consistent with solar wind emph{in-situ} measurements. Electron acceleration is found to be enhanced in the supercritical regime compared to subcritical. When the fully three-dimensional structure of the shock front is resolved, slightly larger energisation for the electrons is observed, and we suggest that this is due to the possibility for the electrons to interact with more than one surface fluctuation per interaction. In the supecritical regime, efficient electron energisation is found also at shock geometries departing from $theta_{Bn}$ very close to 90$^circ$. Two dimensional simulations show indications of unrealistic electron trapping, leading to slightly higher energisation in the subcritical cases.
Upstream of shocks, the foreshock is filled with hot ions. When these ions are concentrated and thermalized around a discontinuity, a diamagnetic cavity bounded by compressional boundaries, referred to as a foreshock transient, forms. Sometimes, the upstream compressional boundary can further steepen into a secondary shock, which has been observed to accelerate particles and contribute to the primary shock acceleration. However, secondary shock formation conditions and processes are not fully understood. Using particle-in-cell simulations, we reveal how secondary shocks are formed. From 1D simulations, we show that electric fields play a critical role in shaping the shocks magnetic field structure, as well as in coupling the energy of hot ions to that of the shock. We demonstrate that larger thermal speed and concentration ratio of hot ions favors the formation of a secondary shock. From a more realistic 2D simulation, we examine how a discontinuity interacts with foreshock ions leading to the formation of a foreshock transient and a secondary shock. Our results imply that secondary shocks are more likely to occur at primary shocks with higher Mach number. With the secondary shocks previously proven ability to accelerate particles in cooperation with a planetary bow shock, it is even more appealing to consider them in particle acceleration of high Mach number astrophysical shocks.