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.

The slow solar wind is typically characterized as having low Alfv\'enicity. However, Parker Solar Probe (PSP) observed predominately Alfv\'enic slow solar wind during several of its initial encounters. From its first encounter observations, about 55.3\% of the slow solar wind inside 0.25 au is highly Alfv\'enic ($|\sigma_C| > 0.7$) at current solar minimum, which is much higher than the fraction of quiet-Sun-associated highly Alfv\'enic slow wind observed at solar maximum at 1 au. Intervals of slow solar wind with different Alfv\'enicities seem to show similar plasma characteristics and temperature anisotropy distributions. Some low Alfv\'enicity slow wind intervals even show high temperature anisotropies, because the slow wind may experience perpendicular heating as fast wind does when close to the Sun. This signature is confirmed by Wind spacecraft measurements as we track PSP observations to 1 au. Further, with nearly 15 years of Wind measurements, we find that the distributions of plasma characteristics, temperature anisotropy and helium abundance ratio ($N_\alpha/N_p$) are similar in slow winds with different Alfv\'enicities, but the distributions are different from those in the fast solar wind. Highly Alfv\'enic slow solar wind contains both helium-rich ($N_\alpha/N_p\sim0.045$) and helium-poor ($N_\alpha/N_p\sim0.015$) populations, implying it may originate from multiple source regions. These results suggest that highly Alfv\'enic slow solar wind shares similar temperature anisotropy and helium abundance properties with regular slow solar winds, and they thus should have multiple origins.

A familiar problem in space and astrophysical plasmas is to understand how dissipation and heating occurs. These effects are often attributed to the cascade of broadband turbulence which transports energy from large scale reservoirs to small scale kinetic degrees of freedom. When collisions are infrequent, local thermodynamic equilibrium is not established. In this case the final stage of energy conversion becomes more complex than in the fluid case, and both pressure-dilatation and pressure strain interactions (Pi-D $\equiv -\Pi_{ij} D_{ij}$) become relevant and potentially important. Pi-D in plasma turbulence has been studied so far primarily using simulations. The present study provides a statistical analysis of Pi-D in the Earth's magnetosheath using the unique measurement capabilities of the Magnetospheric Multiscale (MMS) mission. We find that the statistics of Pi-D in this naturally occurring plasma environment exhibit strong resemblance to previously established fully kinetic simulations results. The conversion of energy is concentrated in space and occurs near intense current sheets, but not within them. This supports recent suggestions that the chain of energy transfer channels involves regional, rather than pointwise, correlations.

When magnetohydrodynamic turbulence evolves in the presence of a large-scale mean magnetic field, an anisotropy develops relative to that preferred direction. The well-known tendency is to develop stronger gradients perpendicular to the magnetic field, relative to the direction along the field. This anisotropy of the spectrum is deeply connected with anisotropy of estimated timescales for dynamical processes, and requires reconsideration of basic issues such as scale locality and spectral transfer. Here analysis of high-resolution three-dimensional simulations of unforced magnetohydrodynamic turbulence permits quantitative assessment of the behavior of theoretically relevant timescales in Fourier wavevector space. We discuss the distribution of nonlinear times, Alfv\'en times, and estimated spectral transfer rates. Attention is called to the potential significance of special regions of the spectrum, such as the two-dimensional limit and the "critical balance" region. A formulation of estimated spectral transfer in terms of a suppression factor supports a conclusion that the quasi two-dimensional fluctuations (characterized by strong nonlinearities) are not a singular limit, but may be in general expected to make important contributions.

The magnetization $|\Omega_{\mathrm e}|/\omega_{\mathrm{e}}$ is an important parameter in plasma astrophysics, where $\Omega_{\mathrm e}$ and $\omega_{\mathrm{e}}$ are the electron gyro-frequency and electron plasma frequency, respectively. It only depends on the mass ratio $m_{\mathrm i}/m_{\mathrm e}$ and the light-to-Alfv\'en speed ratio $c/v_{\mathrm{Ai}}$, where $m_{\mathrm i}$ ($m_{\mathrm e}$) is the ion (electron) mass, $c$ is the speed of light, and $v_{\mathrm{Ai}}$ is the ion Alfv\'en speed. Nonlinear numerical plasma models such as particle-in-cell simulations must often assume unrealistic values for $m_{\mathrm i}/m_{\mathrm e}$ and for $c/v_{\mathrm{Ai}}$. Because linear theory yields exact results for parametric scalings of wave properties at small amplitudes, we use linear theory to investigate the dispersion relations of Alfv\'en/ion-cyclotron and fast-magnetosonic/whistler waves as prime examples for collective plasma behaviour depending on $m_{\mathrm i}/m_{\mathrm e}$ and $c/v_{\mathrm{Ai}}$. We analyse their dependence on $m_{\mathrm i}/m_{\mathrm e}$ and $c/v_{\mathrm{Ai}}$ in quasi-parallel and quasi-perpendicular directions of propagation with respect to the background magnetic field for a plasma with $\beta_j\sim1$, where $\beta_j$ is the ratio of the thermal to magnetic pressure for species $j$. Although their dispersion relations are largely independent of $c/v_{\mathrm{Ai}}$ for $c/v_{\mathrm{Ai}}\gtrsim 10$, the mass ratio $m_{\mathrm i}/m_{\mathrm e}$ has a strong effect at scales smaller than the ion inertial length. Moreover, we study the impact of relativistic electron effects on the dispersion relations. Based on our results, we recommend aiming for a more realistic value of $m_{\mathrm i}/m_{\mathrm e}$ than for a more realistic value of $c/v_{\mathrm{Ai}}$ in non-relativistic plasma simulations if such a choice is necessary, although $\dots$

Using in situ data, accumulated in the turbulent magnetosheath by the Magnetospheric Multiscale (MMS) Mission, we report a statistical study of magnetic field curvature and discuss its role in the turbulent space plasmas. Consistent with previous simulation results, the Probability Distribution Function (PDF) of the curvature is shown to have distinct power-law tails for both high and low value limits. We find that the magnetic-field-line curvature is intermittently distributed in space. High curvature values reside near weak magnetic-field regions, while low curvature values are correlated with small magnitude of the force acting normal to the field lines. A simple statistical treatment provides an explanation for the observed curvature distribution. This novel statistical characterization of magnetic curvature in space plasma provides a starting point for assessing, in a turbulence context, the applicability and impact of particle energization processes, such as curvature drift, that rely on this fundamental quantity.

Observations at 1 au have confirmed that enhancements in measured energetic particle fluxes are statistically associated with "rough" magnetic fields, i.e., fields having atypically large spatial derivatives or increments, as measured by the Partial Variance of Increments (PVI) method. One way to interpret this observation is as an association of the energetic particles with trapping or channeling within magnetic flux tubes, possibly near their boundaries. However, it remains unclear whether this association is a transport or local effect; i.e., the particles might have been energized at a distant location, perhaps by shocks or reconnection, or they might experience local energization or re-acceleration. The Parker Solar Probe (PSP), even in its first two orbits, offers a unique opportunity to study this statistical correlation closer to the corona. As a first step, we analyze the separate correlation properties of the energetic particles measured by the \isois instruments during the first solar encounter. The distribution of time intervals between a specific type of event, i.e., the waiting time, can indicate the nature of the underlying process. We find that the \isois observations show a power-law distribution of waiting times, indicating a correlated (non-Poisson) distribution. Analysis of low-energy \isois data suggests that the results are consistent with the 1 au studies, although we find hints of some unexpected behavior. A more complete understanding of these statistical distributions will provide valuable insights into the origin and propagation of solar energetic particles, a picture that should become clear with future PSP orbits.

Direct evidence of an inertial-range turbulent energy cascade has been provided by spacecraft observations in heliospheric plasmas. In the solar wind, the average value of the derived heating rate near 1 au is $\sim 10^{3}\, \mathrm{J\,kg^{-1}\,s^{-1}}$, an amount sufficient to account for observed departures from adiabatic expansion. Parker Solar Probe (PSP), even during its first solar encounter, offers the first opportunity to compute, in a similar fashion, a fluid-scale energy decay rate, much closer to the solar corona than any prior in-situ observations. Using the Politano-Pouquet third-order law and the von K\'arm\'an decay law, we estimate the fluid-range energy transfer rate in the inner heliosphere, at heliocentric distance $R$ ranging from $54\,R_{\odot}$ (0.25 au) to $36\,R_{\odot}$ (0.17 au). The energy transfer rate obtained near the first perihelion is about 100 times higher than the average value at 1 au. This dramatic increase in the heating rate is unprecedented in previous solar wind observations, including those from Helios, and the values are close to those obtained in the shocked plasma inside the terrestrial magnetosheath.

The solar wind shows periods of highly Alfv\'enic activity, where velocity fluctuations and magnetic fluctuations are aligned or anti-aligned with each other. It is generally agreed that solar wind plasma velocity and magnetic field fluctuations observed by Parker Solar Probe (PSP) during the first encounter are mostly highly Alfv\'enic. However, quantitative measures of Alfv\'enicity are needed to understand how the characterization of these fluctuations compares with standard measures from prior missions in the inner and outer heliosphere, in fast wind and slow wind, and at high and low latitudes. To investigate this issue, we employ several measures to quantify the extent of Alfv\'enicity -- the Alfv\'en ratio $r_A$, {normalized} cross helicity $\sigma_c$, {normalized} residual energy $\sigma_r$, and the cosine of angle between velocity and magnetic fluctuations $\cos\theta_{vb}$. We show that despite the overall impression that the Alfv\'enicity is large in the solar wind sampled by PSP during the first encounter, during some intervals the cross helicity starts decreasing at very large scales. These length-scales (often $> 1000 d_i$) are well inside inertial range, and therefore, the suppression of cross helicity at these scales cannot be attributed to kinetic physics. This drop at large scales could potentially be explained by large-scale shears present in the inner heliosphere sampled by PSP. In some cases, despite the cross helicity being constant down to the noise floor, the residual energy decreases with scale in the inertial range. These results suggest that it is important to consider all these measures to quantify Alfv\'enicity.

During the Parker Solar Probe's (PSP) first perihelion pass, the spacecraft reached within a heliocentric distance of \(\sim 37~R_\odot\) and observed numerous magnetic and flow structures characterized by sharp gradients. To better understand these intermittent structures in the young solar wind, an important property to examine is their degree of correlation in time and space. To this end, we use the well-tested Partial Variance of Increments (PVI) technique to identify intermittent events in FIELDS and SWEAP observations of magnetic and proton-velocity fields (respectively) during PSP's first solar encounter, when the spacecraft was within 0.25 au from the Sun. We then examine distributions of waiting times between events with varying separation and PVI thresholds. We find power-law distributions for waiting times shorter than a characteristic scale comparable to the correlation time, suggesting a high degree of correlation that may originate in a clustering process. Waiting times longer than this characteristic time are better described by an exponential, suggesting a random memory-less Poisson process at play. These findings are consistent with near-Earth observations of solar wind turbulence. The present study complements the one by Dudok de Wit et al. (2020, present volume), which focuses on waiting times between observed "switchbacks" in the radial magnetic field.

Magnetic field are transported and tangled by turbulence, even as they lose identity due to nonideal or resistive effects. On balance field lines undergo stretch-twist-fold processes. The curvature field, a scalar that measures the tangling of the magnetic field lines, is studied in detail here, in the context of magnetohydrodynamic turbulence. A central finding is that the magnitudes of the curvature and the magnetic field are anti-correlated. High curvature co-locates with low magnetic field, which gives rise to power-law tails of the probability density function of the curvature field. The curvature drift term that converts magnetic energy into flow and thermal energy, largely depends on the curvature field behavior, a relationship that helps to explain particle acceleration due to curvature drift. This adds as well to evidence that turbulent effects most likely play an essential role in particle energization since turbulence drives stronger tangled field configurations, and therefore curvature.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

The Parker Solar Probe (PSP) primary mission extends seven years and consists of 24 orbits of the Sun with descending perihelia culminating in a closest approach of ($\sim 9.8~R_\odot$). In the course of these orbits PSP will pass through widely varying conditions, including anticipated large variations of turbulence properties such as energy density, correlation scales and cross helicities. Here we employ global magnetohydrodynamics simulations with self-consistent turbulence transport and heating \citep{usmanov2018} to preview likely conditions that will be encountered by PSP, by assuming suitable boundary conditions at the coronal base. The code evolves large-scale parameters -- such as velocity, magnetic field, and temperature -- as well as turbulent energy density, cross helicity, and correlation scale. These computed quantities provide the basis for evaluating additional useful parameters that are derivable from the primary model outputs. Here we illustrate one such possibility in which computed turbulence and large-scale parameters are used to evaluate the accuracy of the Taylor "frozen-in" hypothesis along the PSP trajectory. Apart from the immediate purpose of anticipating turbulence conditions that PSP will encounter, as experience is gained in comparisons of observations with simulated data, this approach will be increasingly useful for planning and interpretation of subsequent observations.

In the context of space and astrophysical plasma turbulence and particle heating, several vocabularies emerge for estimating turbulent energy dissipation rate, including Kolmogorov-Yaglom third-order law and, in its various forms, $\boldsymbol{j}\cdot\boldsymbol{E}$ (work done by the electromagnetic field on particles), and $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ (pressure-strain interaction), to name a couple. It is now understood that these energy transfer channels, to some extent, are correlated with coherent structures. In particular, we find that different energy dissipation proxies, although not point-wise correlated, are concentrated in proximity to each other, for which they decorrelate in a few $d_i$(s). However, the energy dissipation proxies dominate at different scales. For example, there is an inertial range over which the third-order law is meaningful. Contributions from scale bands stemming from scale-dependent spatial filtering show that, the energy exchange through $\boldsymbol{j}\cdot\boldsymbol{E}$ mainly results from large scales, while the energy conversion from fluid flow to internal through $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ dominates at small scales.

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.

We study the effects of plasma \b{eta} (ratio of plasma pressure to magnetic pressure) on the evolution of kinetic plasma turbulence using fully kinetic particle-in-cell simulations of decaying turbulence. We find that the plasma \b{eta} systematically affects spectra, measures of intermittency, decay rates of turbulence fluctuations, and partitioning over different channels of energy exchange More specifically, an increase in plasma \b{eta} leads to greater total heating, with proton heating preferentially more than electrons. Implications for achieving magnetosheath like temperature ratios are discussed.

Protons (ionized hydrogen) in the solar wind frequently exhibit distinct temperatures ($T_{\perp p}$ and $T_{\parallel p}$) perpendicular and parallel to the plasma's background magnetic-field. Numerous prior studies of the interplanetary solar-wind have shown that, as plasma beta ($\beta_{\parallel p}$) increases, a narrower range of temperature-anisotropy ($R_p\equiv T_{\perp p}\,/\,T_{\parallel p}$) values is observed. Conventionally, this effect has been ascribed to the actions of kinetic microinstabilities. This study is the first to use data from the Magnetospheric Multiscale Mission (MMS) to explore such $\beta_{\parallel p}$-dependent limits on $R_p$ in Earth's magnetosheath. The distribution of these data across the $(\beta_{\parallel p},R_p)$-plane reveals limits on both $R_p>1$ and $R_p<1$. Linear Vlasov theory is used to compute contours of constant growth-rate for the ion-cyclotron, mirror, parallel-firehose, and oblique-firehose instabilities. These instability thresholds closely align with the contours of the data distribution, which suggests a strong association of instabilities with extremes of ion temperature anisotropy in the magnetosheath. The potential for instabilities to regulate temperature anisotropy is discussed.

Plasma turbulence is investigated using high-resolution ion velocity distributions measured by the Magnetospheric Multiscale Mission (MMS) in the Earth's magnetosheath. The particle distribution is highly structured, suggesting a cascade-like process in velocity space. This complex velocity space structure is investigated using a three-dimensional Hermite transform that reveals a power law distribution of moments. In analogy to hydrodynamics, a Kolmogorov approach leads directly to a range of predictions for this phase-space cascade. The scaling theory is in agreement with observations, suggesting a new path for the study of plasma turbulence in weakly collisional space and astrophysical plasmas.

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.

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of the Larmor radius of the charged particle to the characteristic outer length scale of the turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared with those from computer simulations, showing good agreement.

Intermittency of heating in weakly collisional plasma turbulence is an active subject of research, with significant potential impact on understanding of the solar wind, solar corona and astrophysical plasmas. Recent studies suggest a role of vorticity in plasma heating. In magnetohydrodynamics small scale vorticity is generated near current sheets and this effect persists in kinetic plasma, as demonstrated here with hybrid and fully kinetic Particle-In-Cell (PIC) simulations. Furthermore, vorticity enhances local kinetic effects, with a generalized resonance condition selecting sign-dependent enhancements or reductions of proton heating and thermal anisotropy. In such plasmas heating is correlated with vorticity and current density, but more strongly with vorticity. These results help explain several prior results that find kinetic effects and energization near to, but not centered on, current sheets. Evidently intermittency in kinetic plasma involves multiple physical quantities, and the associated coherent structures and nonthermal effects are closely related.

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.