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.

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 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.

Foreshock disturbances -- large-scale (~1000 km to >30,000 km), transient (~5-10 per day - lasting ~10s of seconds to several minutes) structures [1,2] - generated by suprathermal (>100 eV to 100s of keV) ions [3,4] arise upstream of Earth's bow shock formed by the solar wind colliding with the Earth's magnetosphere. They have recently been found to accelerate ions to energies of several keV [5,6]. Although electrons in Saturn's high Mach number (M > 40) bow shock can be accelerated to relativistic energies (nearly 1000 keV) [7], it has hitherto been thought impossible to accelerate electrons at the much weaker (M < 20) Earth's bow shock beyond a few 10s of keV [8]. Here we report observations of electrons energized by foreshock disturbances to energies up to at least ~300 keV. Although such energetic electrons have been previously reported, their presence has been attributed to escaping magnetospheric particles [9,10] or solar events [11]. These relativistic electrons are not associated with any solar activity nor are they of magnetospheric origin. Further, current theories of ion acceleration in foreshock disturbances cannot account for electrons accelerated to the observed relativistic energies [12-17]. These electrons are clearly coming from the disturbances, leaving us with no explanation as to their origin.

Here we present a commentary on the role of small-scale (e.g., few tens of meters) wave-particle interactions in large-scale (e.g., few tens of kilometers) processes and their capacity to accelerate particles from thermal to suprathermal or even to cosmic-ray energies. The purpose of this commentary is to provoke thought and further investigation into the relative importance of electromagnetic waves in the global dynamics of many macroscopic systems. We specifically focus on the terrestrial radiation belts and collisionless shock waves as examples.

In this report, we present a dynamical systems' approach to study the exact nonlinear wave-particle interaction in relativistic regime. We give a particular attention to the effect of wave obliquity on the dynamics of the orbits by studying the specific cases of parallel ($\theta=0$) and perpendicular ($\theta=-\pi/2$) propagations in comparison to the general case of oblique propagation $\theta=]-\pi/2, 0[$. We found that the fixed points of the system correspond to Landau resonance, and that the dynamics can evolve from trapping to surfatron acceleration for propagation angles obeying a Hopf bifurcations condition. Cyclotron-resonant particles are also studied by the construction of a pseudo-potential structure in the Lorentz factor $\gamma$. We derived a condition for which Arnold diffusion results in relativistic stochastic acceleration. Hence, two general conclusions are drawn : 1) The propagation angle $\theta$ can significantly alter the dynamics of the orbits at both Landau and cyclotron-resonances. 2) Considering the short-time scales upon which the particles are accelerated, these two mechanisms for Landau and cyclotron resonant orbits could become potential candidates for problems of particle energization in weakly collisional space and cosmic plasmas.

Using bifurcation theory on a dynamical system simulating the interaction of a particle with an obliquely propagating wave in relativistic regimes, we demonstrate that uniform acceleration arises as a consequence of Hopf bifurcations of Landau resonant particles. The acceleration process arises as a form of surfatron established through the locking in pitch angle, gyrophase, and physical trapping along the wave-vector direction. Integrating the dynamical system for large amplitudes ($\delta B/B_0\sim0.1$) obliquely propagating waves, we find that electrons with initial energies in the keV range can be accelerated to MeV energies on timescales of the order of milliseconds. The Hopf condition of Landau resonant particles could underlie some of the most efficient energization of particles in space and astrophysical plasmas.