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.

We investigate the scattering of strahl electrons by microinstabilities as a mechanism for creating the electron halo in the solar wind. We develop a mathematical framework for the description of electron-driven microinstabilities and discuss the associated physical mechanisms. We find that an instability of the oblique fast-magnetosonic/whistler (FM/W) mode is the best candidate for a microinstability that scatters strahl electrons into the halo. We derive approximate analytic expressions for the FM/W instability threshold in two different $\beta_{\mathrm c}$ regimes, where $\beta_{\mathrm c}$ is the ratio of the core electrons' thermal pressure to the magnetic pressure, and confirm the accuracy of these thresholds through comparison with numerical solutions to the hot-plasma dispersion relation. We find that the strahl-driven oblique FM/W instability creates copious FM/W waves under low-$\beta_{\mathrm c}$ conditions when $U_{0\mathrm s}\gtrsim 3w_{\mathrm c}$, where $U_{0\mathrm s}$ is the strahl speed and $w_{\mathrm c}$ is the thermal speed of the core electrons. These waves have a frequency of about half the local electron gyrofrequency. We also derive an analytic expression for the oblique FM/W instability for $\beta_{\mathrm c}\sim 1$. The comparison of our theoretical results with data from the \emph{Wind} spacecraft confirms the relevance of the oblique FM/W instability for the solar wind. The whistler heat-flux, ion-acoustic heat-flux, kinetic-Alfv\'en-wave heat-flux, and electrostatic electron-beam instabilities cannot fulfill the requirements for self-induced scattering of strahl electrons into the halo. We make predictions for the electron strahl close to the Sun, which will be tested by measurements from \emph{Parker Solar Probe} and \emph{Solar Orbiter}.

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.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

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.

Remarkable progress has been made in understanding turbulent astrophysical plasmas in past decades including, notably, the solar wind and the interstellar medium. In the case of the solar wind, much of this progress has relied on in situ measurements from space-borne instruments. However, ground-based radio observations also have played a significant role and have the potential to play an even bigger role. In particular, using distant background sources (quasars, pulsars, satellite beacons) to transilluminate the foreground corona and solar wind, a variety of radio propagation phenomena can be used to map plasma properties of the solar corona and heliosphere, as well as the warm interstellar medium. These include angular broadening, interplanetary and interstellar scintillations, and differential Faraday rotation. These observations are highly complementary to in situ observations of the solar wind, and could be a mainstay of investigations into turbulence of the ISM. We point out that the Next Generation Very Large Array (ngVLA) fulfills all the requirements necessary to exploit radio observations of astrophysical turbulence fully.

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

Plasma processes occurring in the corona and solar wind can be probed by studying the thermodynamic properties of different ion species. However, most in-situ observations of positive ions in the solar wind are taken at 1 AU, where information on their solar source properties may have been irreversibly erased. In this study we aimed to use the properties of alpha particles at heliocentric distances between 0.3 and 1 AU to study plasma processes occurring at the points of observation, and to infer processes occurring inside 0.3 AU by comparing our results to previous remote sensing observations of the plasma closer to the Sun. We reprocessed the original Helios positive ion distribution functions, isolated the alpha particle population, and computed the alpha particle number density, velocity, and magnetic field perpendicular and parallel temperatures. We then investigated the radial variation of alpha particle temperatures in fast solar wind observed between 0.3 and 1 AU. Between 0.3 and 1 AU alpha particles are heated in the magnetic field perpendicular direction, and cooled in the magnetic field parallel direction. Alpha particle evolution is bounded by the alpha firehose instability threshold, which provides one possible mechanism to explain the observed parallel cooling and perpendicular heating. Closer to the Sun our observations suggest that the alpha particles undergo heating the perpendicular direction, whilst the large magnetic field parallel temperatures observed at 0.3 AU may be due the combined effect of double adiabatic expansion and alpha particle deceleration inside 0.3 AU.

In the near future, Parker Solar Probe and Solar Orbiter will provide the first comprehensive in-situ measurements of the solar wind in the inner heliosphere since the Helios mission in the 1970s. We describe a reprocessing of the original Helios ion distribution functions to provide reliable and reproducible data to characterise the proton core population of the solar wind in the inner heliosphere. A systematic fitting of bi-Maxwellian distribution functions was performed to the raw Helios ion distribution function data to extract the proton core number density, velocity, and temperatures parallel and perpendicular to the magnetic field. We present radial trends of these derived proton parameters, forming a benchmark from which new measurements in the inner heliosphere will be compared to. The new dataset has been made openly available for other researchers to use, along with the source code used to generate it.

The Arbitrary Linear Plasma Solver (ALPS) is a parallelised numerical code that solves the dispersion relation in a hot (even relativistic) magnetised plasma with an arbitrary number of particle species with arbitrary gyrotropic equilibrium distribution functions for any direction of wave propagation with respect to the background field. ALPS reads the background momentum distributions as tables of values on a $(p_{\perp},p_{\parallel})$ grid, where $p_{\perp}$ and $p_{\parallel }$ are the momentum coordinates in the directions perpendicular and parallel to the background magnetic field, respectively. We present the mathematical and numerical approach used by ALPS and introduce our algorithms for the handling of poles and the analytic continuation for the Landau contour integral. We then show test calculations of dispersion relations for a selection of stable and unstable configurations in Maxwellian, bi-Maxwellian, $\kappa$-distributed, and J\"uttner-distributed plasmas. These tests demonstrate that ALPS derives reliable plasma dispersion relations. ALPS will make it possible to determine the properties of waves and instabilities in the non-equilibrium plasmas that are frequently found in space, laboratory experiments, and numerical simulations.

We present a long-duration ($\sim$10 years) statistical analysis of the temperatures, plasma betas, and temperature ratios for the electron, proton, and alpha-particle populations observed by the \emph{Wind} spacecraft near 1 AU. The mean(median) scalar temperatures are $T{\scriptstyle_{e, tot}}$ $=$ 12.2(11.9) eV, $T{\scriptstyle_{p, tot}}$ $=$ 12.7(8.6) eV, and $T{\scriptstyle_{\alpha, tot}}$ $=$ 23.9(10.8) eV. The mean(median) total plasma betas are $\beta{\scriptstyle_{e, tot}}$ $=$ 2.31(1.09), $\beta{\scriptstyle_{p, tot}}$ $=$ 1.79(1.05), and $\beta{\scriptstyle_{\alpha, tot}}$ $=$ 0.17(0.05). The mean(median) temperature ratios are $\left(T{\scriptstyle_{e}}/T{\scriptstyle_{p}}\right){\scriptstyle_{tot}}$ $=$ 1.64(1.27), $\left(T{\scriptstyle_{e}}/T{\scriptstyle_{\alpha}}\right){\scriptstyle_{tot}}$ $=$ 1.24(0.82), and $\left(T{\scriptstyle_{\alpha}}/T{\scriptstyle_{p}}\right){\scriptstyle_{tot}}$ $=$ 2.50(1.94). We also examined these parameters during time intervals that exclude interplanetary (IP) shocks, times within the magnetic obstacles (MOs) of interplanetary coronal mass ejections (ICMEs), and times that exclude MOs. The only times that show significant alterations to any of the parameters examined are those during MOs. In fact, the only parameter that does not show a significant change during MOs is the electron temperature. Although each parameter shows a broad range of values, the vast majority are near the median. We also compute particle-particle collision rates and compare to effective wave-particle collision rates. We find that, for reasonable assumptions of wave amplitude and occurrence rates, the effect of wave-particle interactions on the plasma is equal to or greater than the effect of Coulomb collisions. Thus, wave-particle interactions should not be neglected when modeling the solar wind.

The firehose and mirror instabilities are thought to arise in a variety of space and astrophysical plasmas, constraining the pressure anisotropies and drifts between particle species. The plasma stability depends on all species simultaneously, meaning that a combined analysis is required. Here, we present the first such analysis in the solar wind, using the long-wavelength stability parameters to combine the anisotropies and drifts of all major species (core and beam protons, alphas, and electrons). At the threshold, the firehose parameter was found to be dominated by protons (67%), but also to have significant contributions from electrons (18%) and alphas (15%). Drifts were also found to be important, contributing 57% in the presence of a proton beam. A similar situation was found for the mirror, with contributions of 61%, 28%, and 11% for protons, electrons, and alphas, respectively. The parallel electric field contribution, however, was found to be small at 9%. Overall, the long-wavelength thresholds constrain the data well (<1% unstable), and the implications of this are discussed.

Magnetohydronamic turbulence is believed to play a crucial role in heating the laboratorial, space, and astrophysical plasmas. However, the precise connection between the turbulent fluctuations and the particle kinetics has not yet been established. Here we present clear evidence of plasma turbulence heating based on diagnosed wave features and proton velocity distributions from solar wind measurements by the Wind spacecraft. For the first time, we can report the simultaneous observation of counter-propagating magnetohydrodynamic waves in the solar wind turbulence. Different from the traditional paradigm with counter-propagating Alfv\'en waves, anti-sunward Alfv\'en waves (AWs) are encountered by sunward slow magnetosonic waves (SMWs) in this new type of solar wind compressible turbulence. The counter-propagating AWs and SWs correspond respectively to the dominant and sub-dominant populations of the imbalanced Els\"asser variables. Nonlinear interactions between the AWs and SMWs are inferred from the non-orthogonality between the possible oscillation direction of one wave and the possible propagation direction of the other. The associated protons are revealed to exhibit bi-directional asymmetric beams in their velocity distributions: sunward beams appearing in short and narrow patterns and anti-sunward broad extended tails. It is suggested that multiple types of wave-particle interactions, i.e., cyclotron and Landau resonances with AWs and SMWs at kinetic scales, are taking place to jointly heat the protons perpendicularly and parallel.

The collisionless dissipation of anisotropic Alfv\'enic turbulence is a promising candidate to solve the solar wind heating problem. Extensive studies examined the kinetic properties of Alfv\'en waves in simple Maxwellian or bi-Maxwellian plasmas. However, the observed electron velocity distribution functions in the solar wind are more complex. In this study, we analyze the properties of kinetic Alfv\'en waves in a plasma with two drifting electron populations. We numerically solve the linearized Maxwell-Vlasov equations and find that the damping rate and the proton-electron energy partition for kinetic Alfv\'en waves are significantly modified in such plasmas, compared to plasmas without electron drifts. We suggest that electron drift is an important factor to take into account when considering the dissipation of Alfv\'enic turbulence in the solar wind or other $\beta \sim 1$ astrophysical plasmas.

The goal of the present document is to present the idea of, and convince the community to participate in, Turbulent Dissipation Challenge. The idea was discussed in Solar Heliospheric and Interplanetary ENvironment (SHINE) 2012 meeting. The proponents of the idea Tulasi Parashar and Chadi Salem have prepared this document to circulate the idea in the community. The Turbulent Dissipation Challenge idea is to bring the community together and simulate the same set of problems and try to come to a common set of conclusions about the relative strengths of two different kinds of dissipative processes (current sheets/ reconnection sites vs. wave particle interactions). To take the challenge further, the simulators will provide artificial spacecraft data from the simulations for the observers to analyze.

We use a statistically significant set of measurements to show that the field-aligned electron heat flux $q_\parallel$ in the solar wind at 1 AU is consistent with the Spitzer-H\"{a}rm collisional heat flux $q_{sh}$ for temperature gradient scales larger than a few mean free paths $L_T \gtrsim 3.5 ~\lambda_{fp}$. This represents about 65% of the measured data and corresponds primarily to high $\beta$, weakly collisional plasma ('slow solar wind'). In the more collisionless regime $\lambda_{fp}/L_T \gtrsim 0.28$, the electron heat flux is limited to $q_\parallel/q_0 \sim 0.3$, independent of mean free path, where $q_0$ is the 'free-streaming' value; the measured $q_\parallel$ does not achieve the full $q_0$. This constraint $q_\parallel/q_0 \sim 0.3$ might be attributed to wave-particle interactions, an interplanetary electric potential, or inherent flux limitation. We also show a $\beta_e$ dependence to these results that is consistent with a local radial electron temperature profile $T_e \sim r^{-\alpha}$ that is a function of the thermal electron beta $\alpha = \alpha(\beta_e)$ and that the $\beta$ dependence of the collisionless regulation constraint is not obviously consistent with a whistler heat flux instability. We discuss the results in a broader astrophysical context.

It has long been known that the energy in velocity and magnetic field fluctuations in the solar wind is not in equipartition. In this paper, we present an analysis of 5 years of Wind data at 1 AU to investigate the reason for this. The residual energy (difference between energy in velocity and magnetic field fluctuations) was calculated using both the standard magnetohydrodynamic (MHD) normalization for the magnetic field and a kinetic version, which includes temperature anisotropies and drifts between particle species. It was found that with the kinetic normalization, the fluctuations are closer to equipartition, with a mean normalized residual energy of sigma_r = -0.19 and mean Alfven ratio of r_A = 0.71. The spectrum of residual energy, in the kinetic normalization, was found to be steeper than both the velocity and magnetic field spectra, consistent with some recent MHD turbulence predictions and numerical simulations, having a spectral index close to -1.9. The local properties of residual energy and cross helicity were also investigated, showing that globally balanced intervals with small residual energy contain local patches of larger imbalance and larger residual energy at all scales, as expected for non-linear turbulent interactions.

Kinetic plasma theory is used to generate synthetic spacecraft data to analyze and interpret the compressible fluctuations in the inertial range of solar wind turbulence. The kinetic counterparts of the three familiar linear MHD wave modes---the fast, Alfven, and slow waves---are identified and the properties of the density-parallel magnetic field correlation for these kinetic wave modes is presented. The construction of synthetic spacecraft data, based on the quasi-linear premise---that some characteristics of magnetized plasma turbulence can be usefully modeled as a collection of randomly phased, linear wave modes---is described in detail. Theoretical predictions of the density-parallel magnetic field correlation based on MHD and Vlasov-Maxwell linear eigenfunctions are presented and compared to the observational determination of this correlation based on 10 years of Wind spacecraft data. It is demonstrated that MHD theory is inadequate to describe the compressible turbulent fluctuations and that the observed density-parallel magnetic field correlation is consistent with a statistically negligible kinetic fast wave energy contribution for the large sample used in this study. A model of the solar wind inertial range fluctuations is proposed comprised of a mixture of a critically balanced distribution of incompressible Alfvenic fluctuations and a critically balanced or more anisotropic than critical balance distribution of compressible slow wave fluctuations. These results imply that there is little or no transfer of large scale turbulent energy through the inertial range down to whistler waves at small scales.

We present a measurement of the spectral index of density fluctuations between ion and electron scales in solar wind turbulence using the EFI instrument on the ARTEMIS spacecraft. The mean spectral index at 1 AU was found to be -2.75 +/- 0.06, steeper than predictions for pure whistler or kinetic Alfven wave turbulence, but consistent with previous magnetic field measurements. The steep spectra are also consistent with expectations of increased intermittency or damping of some of the turbulent energy over this range of scales. Neither the spectral index nor the flattening of the density spectra before ion scales were found to depend on the proximity to the pressure anisotropy instability thresholds, suggesting that they are features inherent to the turbulent cascade.

We present the first survey of electric field data using the ARTEMIS spacecraft in the solar wind to study inertial range turbulence. It was found that the average perpendicular spectral index of the electric field depends on the frame of measurement. In the spacecraft frame it is -5/3, which matches the magnetic field due to the large solar wind speed in Lorentz transformation. In the mean solar wind frame, the electric field is primarily due to the perpendicular velocity fluctuations and has a spectral index slightly shallower than -3/2, which is close to the scaling of the velocity. These results are an independent confirmation of the difference in scaling between the velocity and magnetic field, which is not currently well understood. The spectral index of the compressive fluctuations was also measured and found to be close to -5/3, suggesting that they are not only passive to the velocity but may also interact nonlinearly with the magnetic field.

We use a large, statistical set of measurements from the Wind spacecraft at 1 AU, and supporting synthetic spacecraft data based on kinetic plasma theory, to show that the compressible component of inertial range solar wind turbulence is primarily in the kinetic slow mode. The zero-lag cross correlation C(delta n, delta B_parallel) between proton density fluctuations delta n and the field-aligned (compressible) component of the magnetic field delta B_parallel is negative and close to -1. The typical dependence of C(delta n,delta B_parallel) on the ion plasma beta_i is consistent with a spectrum of compressible wave energy that is almost entirely in the kinetic slow mode. This has important implications for both the nature of the density fluctuation spectrum and for the cascade of kinetic turbulence to short wavelengths, favoring evolution to the kinetic Alfven wave mode rather than the (fast) whistler mode.

The proton temperature anisotropy in the solar wind is known to be constrained by the theoretical thresholds for pressure anisotropy-driven instabilities. Here we use approximately 1 million independent measurements of gyroscale magnetic fluctuations in the solar wind to show for the first time that these fluctuations are enhanced along the temperature anisotropy thresholds of the mirror, proton oblique firehose, and ion cyclotron instabilities. In addition, the measured magnetic compressibility is enhanced at high plasma beta ($\beta_\parallel \gtrsim 1$) along the mirror instability threshold but small elsewhere, consistent with expectations of the mirror mode. The power in this frequency (the 'dissipation') range is often considered to be driven by the solar wind turbulent cascade, an interpretation which should be qualified in light of the present results. In particular, we show that the short wavelength magnetic fluctuation power is a strong function of collisionality, which relaxes the temperature anisotropy away from the instability conditions and reduces correspondingly the fluctuation power.