We apply field-particle correlations -- a technique that tracks the time-averaged velocity-space structure of the energy density transfer rate between electromagnetic fields and plasma particles -- to data drawn from a hybrid Vlasov-Maxwell simulation of Alfv\'en Ion-Cyclotron turbulence. Energy transfer in this system is expected to include both Landau and cyclotron wave-particle resonances, unlike previous systems to which the field-particle correlation technique has been applied. In this simulation, the energy transfer rate mediated by the parallel electric field $E_\parallel$ comprises approximately $60\%$ of the total rate, with the remainder mediated by the perpendicular electric field $E_\perp$. The parallel electric field resonantly couples to protons, with the canonical bipolar velocity-space signature of Landau damping identified at many points throughout the simulation. The energy transfer mediated by $E_\perp$ preferentially couples to particles with $v_{tp} \lesssim v_\perp \lesssim 3 v_{tp}$ in agreement with the expected formation of a cyclotron diffusion plateau. Our results demonstrate clearly that the field-particle correlation technique can distinguish distinct channels of energy transfer using single-point measurements, even at points in which multiple channels act simultaneously, and can be used to determine quantitatively the rates of particle energization in each channel.

Monte Carlo methods are often employed to numerically integrate kinetic equations, such as the particle-in-cell method for the plasma kinetic equation, but these methods suffer from the introduction of counting noise to the solution. We report on a cautionary tale of counting noise modifying the nonlinear saturation of kinetic instabilities driven by unstable beams of plasma. We find a saturated magnetic field in under-resolved particle-in-cell simulations due to the sampling error in the current density. The noise-induced magnetic field is anomalous, as the magnetic field damps away in continuum kinetic and increased particle count particle-in-cell simulations. This modification of the saturated state has implications for a broad array of astrophysical phenomena beyond the simple plasma system considered here, and it stresses the care that must be taken when using particle methods for kinetic equations.

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio $Q_{\text{i}}/Q_{\text{e}}$ as a function of the compressive-to-Alfv\'enic driving power ratio $P_{\text{compr}}/P_{\text{AW}}$, of the ratio of ion thermal pressure to magnetic pressure $\beta_{\text{i}}$, and of the ratio of ion-to-electron background temperatures $T_{\text{i}}/T_{\text{e}}$. It is shown that $Q_{\text{i}}/Q_{\text{e}}$ is an increasing function of $P_{\text{compr}}/P_{\text{AW}}$. When the compressive driving is sufficiently large, $Q_{\text{i}}/Q_{\text{e}}$ approaches $\simeq P_{\text{compr}}/P_{\text{AW}}$. This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfv\'en waves is absent for both low and high $\beta_{\text{i}}$, meaning that the compressive driving is directly connected to the ion entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low $\beta_{\text{i}}$ and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

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

Turbulence is ubiquitously observed in nearly collisionless heliospheric plasmas, including the solar wind and corona and the Earth's magnetosphere. Understanding the collisionless mechanisms responsible for the energy transfer from the turbulent fluctuations to the particles is a frontier in kinetic turbulence research. Collisionless energy transfer from the turbulence to the particles can take place reversibly, resulting in non-thermal energy in the particle velocity distribution functions (VDFs) before eventual collisional thermalization is realized. Exploiting the information contained in the fluctuations in the VDFs is valuable. Here we apply a recently developed method based on VDFs, the field-particle correlation technique, to a $\beta=1$, solar-wind-like, low-frequency Alfv\'enic turbulence simulation with well resolved phase space to identify the field-particle energy transfer in velocity space. The field-particle correlations reveal that the energy transfer, mediated by the parallel electric field, results in significant structuring of the ion and electron VDFs in the direction parallel to the magnetic field. Fourier modes representing the length scales between the ion and electron gyroradii show that energy transfer is resonant in nature, localized in velocity space to the Landau resonances for each Fourier mode. The energy transfer closely follows the Landau resonant velocities with varying perpendicular wavenumber $k_\perp$ and plasma $\beta$. This resonant signature, consistent with Landau damping, is observed in all diagnosed Fourier modes that cover the dissipation range of the simulation.

A multi-institutional, multi-national science team will soon submit a NASA proposal to build a constellation of spacecraft to fly into the near-Earth solar wind in a swarm spanning a multitude of scales in order to obtain critically needed measurements that will reveal the underlying dynamics of magnetized turbulence. This white paper, submitted to the Plasma 2020 Decadal Survey Committee, provides a brief overview of turbulent systems that constitute an area of compelling plasma physics research, including why this mission is needed, and how this mission will achieve the goal of revealing how energy is transferred across scales and boundaries in plasmas throughout the universe.

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 Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multi-fluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. We focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations.

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.

We present the first 2X2V continuum Vlasov-Maxwell simulations of interpenetrating, unmagnetized plasmas to study the competition between two-stream, Oblique, and filamentation modes in the weakly relativistic regime. We find that after nonlinear saturation of the fastest-growing two-stream and Oblique modes, the effective temperature anisotropy, which drives current filament formation via the secular Weibel instability, has a strong dependence on the internal temperature of the counter-streaming plasmas. The effective temperature anisotropy is significantly more reduced in colder than in hotter plasmas, leading to orders of magnitude lower magnetization for colder plasmas. A strong dependence of the energy conversion efficiency of Weibel-type instabilities on internal beam temperature has implications for determining their contribution to the observed magnetization of many astrophysical and laboratory plasmas.

We present a new algorithm for the discretization of the Vlasov-Maxwell system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin finite element method for the spatial discretization, we obtain a high order accurate solution for the plasma's distribution function. Time stepping for the distribution function is done explicitly with a third order strong-stability preserving Runge-Kutta method. Since the Vlasov equation in the Vlasov-Maxwell system is a high dimensional transport equation, up to six dimensions plus time, we take special care to note various features we have implemented to reduce the cost while maintaining the integrity of the solution, including the use of a reduced high-order basis set. A series of benchmarks, from simple wave and shock calculations, to a five dimensional turbulence simulation, are presented to verify the efficacy of our set of numerical methods, as well as demonstrate the power of the implemented features.

Determining the physical mechanisms that extract energy from turbulent fluctuations in weakly collisional magnetized plasmas is necessary for a more complete characterization of the behavior of a variety of space and astrophysical plasmas. Such a determination is complicated by the complex nature of the turbulence as well as observational constraints, chiefly that in situ measurements of such plasmas are typically only available at a single point in space. Recent work has shown that correlations between electric fields and particle velocity distributions constructed from single-point measurements produce a velocity-dependent signature of the collisionless damping mechanism. We extend this work by constructing field-particle correlations using data sets drawn from single points in strongly driven, turbulent, electromagnetic gyrokinetic simulations to demonstrate that this technique can identify the collisionless mechanisms operating in such systems. The correlation's velocity-space structure agrees with expectations of resonant mechanisms transferring energy collisionlessly in turbulent systems. This work motivates the eventual application of field-particle correlations to spacecraft measurements in the solar wind, with the ultimate goal to determine the physical mechanisms that dissipate magnetized plasma turbulence.

Plasma turbulence is ubiquitous in space and astrophysical plasmas, playing an important role in plasma energization, but the physical mechanisms leading to dissipation of the turbulent energy remain to be definitively identified. Kinetic simulations in two dimensions (2D) have been extensively used to study the dissipation process. How the limitation to 2D affects energy dissipation remains unclear. This work provides a model of comparison between two- and three-dimensional (3D) plasma turbulence using gyrokinetic simulations; it also explores the dynamics of distribution functions during the dissipation process. It is found that both 2D and 3D nonlinear gyrokinetic simulations of a low-beta plasma generate electron velocity-space structures with the same characteristics as that of linear Landau damping of Alfv\'en waves in a 3D linear simulation. The continual occurrence of the velocity-space structures throughout the turbulence simulations suggests that the action of Landau damping may be responsible for the turbulent energy transfer to electrons in both 2D and 3D, and makes possible the subsequent irreversible heating of the plasma through collisional smoothing of the velocity-space fluctuations. Although, in the 2D case where variation along the equilibrium magnetic field is absent, it may be expected that Landau damping is not possible, a common trigonometric factor appears in the 2D resonant denominator, leaving the resonance condition unchanged from the 3D case. The evolution of the 2D and 3D cases is qualitatively similar. However, quantitatively the nonlinear energy cascade and subsequent dissipation is significantly slower in the 2D case.

Motivated by the upcoming Solar Orbiter and Solar Probe Plus missions, qualitative and quantitative predictions are made for the effects of the violation of the Taylor hypothesis on the magnetic energy frequency spectrum measured in the near-Sun environment. The synthetic spacecraft data method is used to predict observational signatures of the violation for critically balanced Alfv\'enic turbulence or parallel fast/whistler turbulence. The violation of the Taylor hypothesis can occur in the slow flow regime, leading to a shift of the entire spectrum to higher frequencies, or in the dispersive regime, in which the dissipation range spectrum flattens at high frequencies. It is found that Alfv\'enic turbulence will not significantly violate the Taylor hypothesis, but whistler turbulence will. The flattening of the frequency spectrum is therefore a key observational signature for fast/whistler turbulence.

The interpretation of single-point spacecraft measurements of solar wind turbulence is complicated by the fact that the measurements are made in a frame of reference in relative motion with respect to the turbulent plasma. The Taylor hypothesis---that temporal fluctuations measured by a stationary probe in a rapidly flowing fluid are dominated by the advection of spatial structures in the fluid rest frame---is often assumed to simplify the analysis. But measurements of turbulence in upcoming missions, such as Solar Probe Plus, threaten to violate the Taylor hypothesis, either due to slow flow of the plasma with respect to the spacecraft or to the dispersive nature of the plasma fluctuations at small scales. Assuming that the frequency of the turbulent fluctuations is characterized by the frequency of the linear waves supported by the plasma, we evaluate the validity of the Taylor hypothesis for the linear kinetic wave modes in the weakly collisional solar wind. The analysis predicts that a dissipation range of solar wind turbulence supported by whistler waves is likely to violate the Taylor hypothesis, while one supported by kinetic Alfven waves is not.

The quasilinear premise is a hypothesis for the modeling of plasma turbulence in which the turbulent fluctuations are represented by a superposition of randomly-phased linear wave modes, and energy is transferred among these wave modes via nonlinear interactions. We define specifically what constitutes the quasilinear premise, and present a range of theoretical arguments in support of the relevance of linear wave properties even in a strongly turbulent plasma. We review evidence both in support of and in conflict with the quasilinear premise from numerical simulations and measurements of plasma turbulence in the solar wind. Although the question of the validity of the quasilinear premise remains to be settled, we suggest that the evidence largely supports the value of the quasilinear premise in modeling plasma turbulence and that its usefulness may also be judged by the insights gained from such an approach, with the ultimate goal to develop the capability to predict the evolution of any turbulent plasma system, including the spectrum of turbulent fluctuations, their dissipation, and the resulting plasma heating.

Motivated by recent observations of distinct parallel and perpendicular signatures in magnetic helicity measurements segregated by wave period and angle between the local magnetic field and the solar wind velocity, this paper undertakes a comparison of three intervals of \emph{Ulysses} data with synthetic time series generated from a physically motivated turbulence model. From these comparisons, it is hypothesized that the observed signatures result from a perpendicular cascade of \Alfvenic fluctuations and a local, non-turbulent population of ion cyclotron or whistler waves generated by temperature anisotropy instabilities. By constraining the model's free parameters through comparison to \emph{in situ} data, it is found that, on average, $ \sim 95\%$ of the power near dissipative scales is contained in a perpendicular \Alfvenic cascade and that the parallel fluctuations are propagating nearly unidirectionally. The effects of aliasing on magnetic helicity measurements are considered and shown to be significant near the Nyquist frequency.

Results of the first validation of large guide field, $B_g / \delta B_0 \gg 1$, gyrokinetic simulations of magnetic reconnection at a fusion and solar corona relevant $\beta_i = 0.01$ and solar wind relevant $\beta_i = 1$ are presented, where $\delta B_0$ is the reconnecting field. Particle-in-cell (PIC) simulations scan a wide range of guide magnetic field strength to test for convergence to the gyrokinetic limit. The gyrokinetic simulations display a high degree of morphological symmetry, to which the PIC simulations converge when $\beta_i B_g / \delta B_0 \gtrsim 1$ and $B_g / \delta B_0 \gg 1$. In the regime of convergence, the reconnection rate, relative energy conversion, and overall magnitudes are found to match well between the PIC and gyrokinetic simulations, implying that gyrokinetics is capable of making accurate predictions well outside its regime of formal applicability. These results imply that in the large guide field limit many quantities resulting from the nonlinear evolution of reconnection scale linearly with the guide field.

A unique method of driving Alfvenic turbulence via an oscillating Langevin antenna is presented. This method of driving is motivated by a desire to inject energy into a finite domain numerical simulation in a manner that models the nonlinear transfer of energy from fluctuations in the turbulent cascade at scales larger than the simulation domain.. The oscillating Langevin antenna is shown to capture the essential features of the larger scale turbulence and efficiently couple to the plasma, generating steady-state turbulence within one characteristic turnaround time. The antenna is also sufficiently flexible to explore both strong and weak regimes of Alfvenic plasma turbulence.

We present the first study of the formation and dissipation of current sheets at electron scales in a wave-driven, weakly collisional, 3D kinetic turbulence simulation. We investigate the relative importance of dissipation associated with collisionless damping via resonant wave-particle interactions versus dissipation in small-scale current sheets in weakly collisional plasma turbulence. Current sheets form self-consistently from the wave-driven turbulence, and their filling fraction is well correlated to the electron heating rate. However, the weakly collisional nature of the simulation necessarily implies that the current sheets are not significantly dissipated via Ohmic dissipation. Rather, collisionless damping via the Landau resonance with the electrons is sufficient to account for the measured heating as a function of scale in the simulation, without the need for significant Ohmic dissipation. This finding suggests the possibility that the dissipation of the current sheets is governed by resonant wave-particle interactions and that the locations of current sheets correspond spatially to regions of enhanced heating.

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.

The magnetic variance anisotropy ($\mathcal{A}_m$) of the solar wind has been used widely as a method to identify the nature of solar wind turbulent fluctuations; however, a thorough discussion of the meaning and interpretation of the $\mathcal{A}_m$ has not appeared in the literature. This paper explores the implications and limitations of using the $\mathcal{A}_m$ as a method for constraining the solar wind fluctuation mode composition and presents a more informative method for interpreting spacecraft data. The paper also compares predictions of the $\mathcal{A}_m$ from linear theory to nonlinear turbulence simulations and solar wind measurements. In both cases, linear theory compares well and suggests the solar wind for the interval studied is dominantly Alfv\'{e}nic in the inertial and dissipation ranges to scales $k \rho_i \simeq 5$.

A numerical simulation of kinetic plasma turbulence is performed to assess the applicability of critical balance to kinetic, dissipation scale turbulence. The analysis is performed in the frequency domain to obviate complications inherent in performing a local analysis of turbulence. A theoretical model of dissipation scale critical balance is constructed and compared to simulation results, and excellent agreement is found. This result constitutes the first evidence of critical balance in a kinetic turbulence simulation and provides evidence of an anisotropic turbulence cascade extending into the dissipation range. We also perform an Eulerian frequency analysis of the simulation data and compare it to the results of a previous study of magnetohydrodynamic turbulence simulations.

A refined cascade model for kinetic turbulence in weakly collisional astrophysical plasmas is presented that includes both the transition between weak and strong turbulence and the effect of nonlocal interactions on the nonlinear transfer of energy. The model describes the transition between weak and strong MHD turbulence and the complementary transition from strong kinetic Alfven wave (KAW) turbulence to weak dissipating KAW turbulence, a new regime of weak turbulence in which the effects of shearing by large scale motions and kinetic dissipation play an important role. The inclusion of the effect of nonlocal motions on the nonlinear energy cascade rate in the dissipation range, specifically the shearing by large-scale motions, is proposed to explain the nearly power-law energy spectra observed in the dissipation range of both kinetic numerical simulations and solar wind observations.

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 first three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms. The resulting energy spectra are quantitatively consistent with a magnetic power spectrum scaling of $k^{-2.8}$ as observed in \emph{in situ} spacecraft measurements of the "dissipation range" of solar wind turbulence. Despite the strongly nonlinear nature of the turbulence, the linear kinetic \Alfven wave mode quantitatively describes the polarization of the turbulent fluctuations. The collisional ion heating is measured at sub-ion-Larmor radius scales, which provides the first evidence of the ion entropy cascade in an electromagnetic turbulence simulation.

Weak gravitational lensing due to the inhomogeneous matter distribution in the universe is an important systematic uncertainty in the use of standard candles in cosmology. There are two different descriptions of weak lensing amplification, one uses a local smoothness parameter $\tilde{\alpha}$, the other uses reduced convergence $\eta= 1+ \kappa/|\kappa_{min}|$ (where $\kappa$ is convergence). The $\tilde{\alpha}$ description involves Dyer-Roeder distance $D_A(\tilde{\alpha}|z)$ ($\tilde{\alpha}=1$ corresponds to a smooth universe); it is simple and convenient, and has been used by the community to illustrate the effect of weak lensing on point sources such as type Ia supernovae. Wang (1999) has shown that the $\tilde{\alpha}$ description can be made realistic by allowing $\tilde{\alpha}$ to be a local variable, the local smoothness parameter. The $\eta$ description has been used by Wang, Holz, & Munshi (2002) to derive a universal probability distribution (UPDF) for weak lensing amplification. In this paper, we bridge the two different descriptions of weak lensing amplification by relating the reduced convergence $\eta$ and the local smoothness parameter $\tilde{\alpha}$. We give the variance of $\tilde{\alpha}$ in terms of the matter power spectrum, thus providing a quantitative guidance to the use of Dyer-Roeder distances in illustrating the effect of weak lensing. The by-products of this work include a corrected definition of the reduced convergence, and simple and accurate analytical expressions for $D_A(\tilde{\alpha}|z)$. Our results should be very useful in studying the weak lensing of standard candles.