Both kinetic instabilities and strong turbulence have potential to impact the behavior of space plasmas. To assess effects of these two processes we compare results from a 3 dimensional particle-in-cell (PIC) simulation of collisionless plasma turbulence against observations by the MMS spacecraft in the terrestrial magnetosheath and by the Wind spacecraft in the solar wind. The simulation develops coherent structures and anisotropic ion velocity distributions that can drive micro-instabilities. Temperature-anisotropy driven instability growth rates are compared with inverse nonlinear turbulence time scales. Large growth rates occur near coherent structures; nevertheless linear growth rates are, on average, substantially less than the corresponding nonlinear rates. This result casts some doubt on the usual basis for employing linear instability theory, and raises questions as to why the linear theory appears to work in limiting plasma excursions in anisotropy and plasma beta.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

We study spectral features of ion velocity and magnetic field correlations in the solar wind and in the magnetosheath using data from the Magnetospheric Multi-Scale (MMS) spacecraft. High resolution MMS observations enable the study of transition of these correlations between their magnetofluid character at larger scales into the sub-proton kinetic range, previously unstudied in spacecraft data. Cross-helicity, angular alignment and energy partitioning is examined over a suit- able range of scales, employing measurements based on the Taylor frozen-in approximation as well as direct two-spacecraft correlation measurements. The results demonstrate signatures of alignment at large scales. As kinetic scales are approached, the alignment between v and b is destroyed by demagnetization of protons.

The heliosphere is formed due to interaction between the solar wind (SW) and local interstellar medium (LISM). The shape and position of the heliospheric boundary, the heliopause, in space depend on the parameters of interacting plasma flows. The interplay between the asymmetrizing effect of the interstellar magnetic field and charge exchange between ions and neutral atoms plays an important role in the SW-LISM interaction. By performing three-dimensional, MHD plasma / kinetic neutral atom simulations, we determine the width of the outer heliosheath - the LISM plasma region affected by the presence of the heliosphere - and analyze quantitatively the distributions in front of the heliopause. It is shown that charge exchange modifies the LISM plasma to such extent that the contribution of a shock transition to the total variation of plasma parameters becomes small even if the LISM velocity exceeds the fast magnetosonic speed in the unperturbed medium. By performing adaptive mesh refinement simulations, we show that a distinct boundary layer of decreased plasma density and enhanced magnetic field should be observed on the interstellar side of the heliopause. We show that this behavior is in agreement with the plasma oscillations of increasing frequency observed by the plasma wave instrument onboard Voyager 1. We also demonstrate that Voyager observations in the inner heliosheath between the heliospheric termination shock and the heliopause are consistent with dissipation of the heliospheric magnetic field. The choice of LISM parameters in this analysis is based on the simulations that fit observations of energetic neutral atoms performed by IBEX.

Kinetic plasma turbulence cascade spans multiple scales ranging from macroscopic fluid flow to sub-electron scales. Mechanisms that dissipate large scale energy, terminate the inertial range cascade and convert kinetic energy into heat are hotly debated. Here we revisit these puzzles using fully kinetic simulation. By performing scale-dependent spatial filtering on the Vlasov equation, we extract information at prescribed scales and introduce several energy transfer functions. This approach allows highly inhomogeneous energy cascade to be quantified as it proceeds down to kinetic scales. The pressure work, $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$, can trigger a channel of the energy conversion between fluid flow and random motions, which is a collision-free generalization of the viscous dissipation in collisional fluid. Both the energy transfer and the pressure work are strongly correlated with velocity gradients.