Local kinetic effects in two-dimensional plasma turbulence

Quantifying Energy Conversion in Higher-Order Phase Space Density Moments in Plasmas

Weakly collisional and collisionless plasmas are typically far from local thermodynamic equilibrium (LTE), and understanding energy conversion in such systems is a forefront research problem. The standard approach is to investigate changes in internal (thermal) energy and density, but this omits energy conversion that changes any higher-order moments of the phase space density. In this Letter, we calculate from first principles the energy conversion associated with all higher moments of the phase space density for systems not in LTE. Particle-in-cell simulations of collisionless magnetic reconnection reveal that energy conversion associated with higher-order moments can be locally significant. The results may be useful in numerous plasma settings, such as reconnection, turbulence, shocks, and wave-particle interactions in heliospheric, planetary, and astrophysical plasmas.

Observations of cross scale energy transfer in the inner heliosphere by Parker Solar Probe

The solar wind, a continuous flow of plasma from the sun, not only shapes the near Earth space environment but also serves as a natural laboratory to study plasma turbulence in conditions that are not achievable in the lab. Starting with the Mariners, for more than five decades, multiple space missions have enabled in-depth studies of solar wind turbulence. Parker Solar Probe (PSP) was launched to explore the origins and evolution of the solar wind. With its state-of-the-art instrumentation and unprecedented close approaches to the sun, PSP is starting a new era of inner heliospheric exploration. In this review we discuss observations of turbulent energy flow across scales in the inner heliosphere as observed by PSP. After providing a quick theoretical overview and a quick recap of turbulence before PSP, we discuss in detail the observations of energy at various scales on its journey from the largest scales to the internal degrees of freedom of the plasma. We conclude with some open ended questions, many of which we hope that PSP will help answer.

Turbulence-Driven Ion Beams in the Magnetospheric Kelvin-Helmholtz Instability

The description of the local turbulent energy transfer and the high-resolution ion distributions measured by the Magnetospheric Multiscale mission together provide a formidable tool to explore the cross-scale connection between the fluid-scale energy cascade and plasma processes at subion scales. When the small-scale energy transfer is dominated by Alfvénic, correlated velocity, and magnetic field fluctuations, beams of accelerated particles are more likely observed. Here, for the first time, we report observations suggesting the nonlinear wave-particle interaction as one possible mechanism for the energy dissipation in space plasmas.

Particle Acceleration in Relativistic Plasma Turbulence

Due to its ubiquitous presence, turbulence is often invoked to explain the origin of nonthermal particles in astrophysical sources of high-energy emission. With particle-in-cell simulations, we study decaying turbulence in magnetically dominated (or, equivalently, "relativistic") pair plasmas. We find that the generation of a power-law particle energy spectrum is a generic by-product of relativistic turbulence. The power-law slope is harder for higher magnetizations and stronger turbulence levels. In large systems, the slope attains an asymptotic, system-size-independent value, while the high-energy spectral cutoff increases linearly with system size; both the slope and the cutoff do not depend on the dimensionality of our domain. By following a large sample of particles, we show that particle injection happens at reconnecting current sheets; the injected particles are then further accelerated by stochastic interactions with turbulent fluctuations. Our results have important implications for the origin of nonthermal particles in high-energy astrophysical sources.

Vlasov methods in space physics and astrophysics

This paper reviews Vlasov-based numerical methods used to model plasma in space physics and astrophysics. Plasma consists of collectively behaving charged particles that form the major part of baryonic matter in the Universe. Many concepts ranging from our own planetary environment to the Solar system and beyond can be understood in terms of kinetic plasma physics, represented by the Vlasov equation. We introduce the physical basis for the Vlasov system, and then outline the associated numerical methods that are typically used. A particular application of the Vlasov system is Vlasiator, the world’s first global hybrid-Vlasov simulation for the Earth’s magnetic domain, the magnetosphere. We introduce the design strategies for Vlasiator and outline its numerical concepts ranging from solvers to coupling schemes. We review Vlasiator’s parallelisation methods and introduce the used high-performance computing (HPC) techniques. A short review of verification, validation and physical results is included. The purpose of the paper is to present the Vlasov system and introduce an example implementation, and to illustrate that even with massive computational challenges, an accurate description of physics can be rewarding in itself and significantly advance our understanding. Upcoming supercomputing resources are making similar efforts feasible in other fields as well, making our design options relevant for others facing similar challenges.

Exact hybrid Vlasov equilibria for sheared plasmas with in-plane and out-of-plane magnetic field

The hybrid Vlasov-Maxwell system of equations is suitable to describe a magnetized plasma at scales on the order of or larger than proton kinetic scales. An exact stationary solution is presented by revisiting previous results with a uniform-density shear flow, directed either parallel or perpendicular to a uniform magnetic field, and by adapting the solution to the hybrid Vlasov-Maxwell model. A quantitative characterization of the equilibrium distribution function is provided by studying both analytically and numerically the temperature anisotropy and gyrotropy and the heat flux. In both cases, in the shear region, the velocity distribution significantly departs from local thermodynamical equilibrium. A comparison between the time behavior of the usual "fluidlike" equilibrium shifted Maxwellian and the exact stationary solutions is carried out by means of numerical simulations of the hybrid Vlasov-Maxwell equations. These hybrid equilibria can be employed as unperturbed states for numerous problems which involve sheared flows, such as the wave propagation in an inhomogeneous background and the onset of the Kelvin-Helmholtz instability.

Coherent Structures and Spectral Energy Transfer in Turbulent Plasma: A Space-Filter Approach

Plasma turbulence at scales of the order of the ion inertial length is mediated by several mechanisms, including linear wave damping, magnetic reconnection, the formation and dissipation of thin current sheets, and stochastic heating. It is now understood that the presence of localized coherent structures enhances the dissipation channels and the kinetic features of the plasma. However, no formal way of quantifying the relationship between scale-to-scale energy transfer and the presence of spatial structures has been presented so far. In the Letter we quantify such a relationship analyzing the results of a two-dimensional high-resolution Hall magnetohydrodynamic simulation. In particular, we employ the technique of space filtering to derive a spectral energy flux term which defines, in any point of the computational domain, the signed flux of spectral energy across a given wave number. The characterization of coherent structures is performed by means of a traditional two-dimensional wavelet transformation. By studying the correlation between the spectral energy flux and the wavelet amplitude, we demonstrate the strong relationship between scale-to-scale transfer and coherent structures. Furthermore, by conditioning one quantity with respect to the other, we are able for the first time to quantify the inhomogeneity of the turbulence cascade induced by topological structures in the magnetic field. Taking into account the low space-filling factor of coherent structures (i.e., they cover a small portion of space), it emerges that 80% of the spectral energy transfer (both in the direct and inverse cascade directions) is localized in about 50% of space, and 50% of the energy transfer is localized in only 25% of space.

Magnetospheric Multiscale Observation of Plasma Velocity-Space Cascade: Hermite Representation and Theory

Plasma turbulence is investigated using unprecedented high-resolution ion velocity distribution measurements by the Magnetospheric Multiscale mission (MMS) in the Earth's magnetosheath. This novel observation of a highly structured particle distribution suggests a cascadelike process in velocity space. Complex velocity space structure is investigated using a three-dimensional Hermite transform, revealing, for the first time in observational data, a power-law distribution of moments. In analogy to hydrodynamics, a Kolmogorov approach leads directly to a range of predictions for this phase-space transport. The scaling theory is found to be in agreement with observations. The combined use of state-of-the-art MMS data sets, novel implementation of a Hermite transform method, and scaling theory of the velocity cascade opens new pathways to the understanding of plasma turbulence and the crucial velocity space features that lead to dissipation in plasmas.

Turbulence generation during the head-on collision of Alfvénic wave packets

The description of the Moffatt and Parker problem recently revisited by O. Pezzi et al. [Astrophys. J. 834, 166 (2017)1538-435710.3847/1538-4357/834/2/166] is here extended by analyzing the features of the turbulence produced by the interaction of two colliding Alfvénic wave packets in a kinetic plasma. Although the approach based on the presence of linear modes features is still helpful in characterizing some low-energy fluctuations, other signatures, which go beyond the pure linear modes analysis, are recovered, such as the significant weakening of clear dispersion relations and the production of zero frequency fluctuations.

Energy transfer channels and turbulence cascade in Vlasov-Maxwell turbulence

Analysis of the Vlasov-Maxwell equations from the perspective of turbulence cascade clarifies the role of electromagnetic work, and reveals the importance of the pressure-strain relation in generating internal energy. Particle-in-cell simulation demonstrates the relative importance of the several energy exchange terms, indicating that the traceless pressure-strain interaction "Pi-D" is of particular importance for both electrons and protons. The Pi-D interaction and the second tensor invariants of the strain are highly localized in similar spatial regions, indicating that energy transfer occurs preferentially in coherent structures. The collisionless turbulence cascade may be fruitfully explored by study of these energy transfer channels, in addition to examining transfer across spatial scales.

A prospectus on kinetic heliophysics

Under the low density and high temperature conditions typical of heliospheric plasmas, the macroscopic evolution of the heliosphere is strongly affected by the kinetic plasma physics governing fundamental microphysical mechanisms. Kinetic turbulence, collisionless magnetic reconnection, particle acceleration, and kinetic instabilities are four poorly understood, grand-challenge problems that lie at the new frontier of kinetic heliophysics. The increasing availability of high cadence and high phase-space resolution measurements of particle velocity distributions by current and upcoming spacecraft missions and of massively parallel nonlinear kinetic simulations of weakly collisional heliospheric plasmas provides the opportunity to transform our understanding of these kinetic mechanisms through the full utilization of the information contained in the particle velocity distributions. Several major considerations for future investigations of kinetic heliophysics are examined. Turbulent dissipation followed by particle heating is highlighted as an inherently two-step process in weakly collisional plasmas, distinct from the more familiar case in fluid theory. Concerted efforts must be made to tackle the big-data challenge of visualizing the high-dimensional (3D-3V) phase space of kinetic plasma theory through physics-based reductions. Furthermore, the development of innovative analysis methods that utilize full velocity-space measurements, such as the field-particle correlation technique, will enable us to gain deeper insight into these four grand-challenge problems of kinetic heliophysics. A systems approach to tackle the multi-scale problem of heliophysics through a rigorous connection between the kinetic physics at microscales and the self-consistent evolution of the heliosphere at macroscales will propel the field of kinetic heliophysics into the future.

Explosive Particle Dispersion in Plasma Turbulence

Particle dynamics are investigated in plasma turbulence, using self-consistent kinetic simulations, in two dimensions. In the steady state, the trajectories of single protons and proton pairs are studied, at different values of plasma β (ratio between kinetic and magnetic pressure). For single-particle displacements, results are consistent with fluids and magnetic field line dynamics, where particles undergo normal diffusion for very long times, with higher β's being more diffusive. In an intermediate time range, with separations lying in the inertial range, particles experience an explosive dispersion in time, consistent with the Richardson prediction. These results, obtained for the first time with a self-consistent kinetic model, are relevant for astrophysical and laboratory plasmas, where turbulence is crucial for heating, mixing, and acceleration processes.

Pressure anisotropy and small spatial scales induced by velocity shear

By including the full pressure tensor dynamics in a fluid plasma model, we show that a sheared velocity field can provide an effective mechanism that makes the initial isotropic pressure nongyrotropic. This is distinct from the usual gyrotropic anisotropy related to the fluid compressibility and usually accounted for in double-adiabatic models. We determine the time evolution of the pressure agyrotropy and discuss how the propagation of "magnetoelastic perturbations" can affect the pressure tensor anisotropization and its spatial filamentation, which are due to the action of both the magnetic field and the flow strain tensor. We support this analysis with a numerical integration of the nonlinear equations describing the pressure tensor evolution.

Collisional Relaxation of Fine Velocity Structures in Plasmas

The existence of several characteristic times during the collisional relaxation of fine velocity structures is investigated by means of Eulerian numerical simulations of a spatially homogeneous force-free weakly collisional plasma. The effect of smoothing out velocity gradients on the evolution of global quantities, such as temperature and entropy, is discussed, suggesting that plasma collisionality can locally increase due to velocity space deformations of the particle velocity distribution function. These results support the idea that high-resolution measurements of the particle velocity distribution function are crucial for an accurate description of weakly collisional systems, such as the solar wind, in order to answer relevant scientific questions, related, for example, to particle heating and energization.

Intermittency, nonlinear dynamics and dissipation in the solar wind and astrophysical plasmas

An overview is given of important properties of spatial and temporal intermittency, including evidence of its appearance in fluids, magnetofluids and plasmas, and its implications for understanding of heliospheric plasmas. Spatial intermittency is generally associated with formation of sharp gradients and coherent structures. The basic physics of structure generation is ideal, but when dissipation is present it is usually concentrated in regions of strong gradients. This essential feature of spatial intermittency in fluids has been shown recently to carry over to the realm of kinetic plasma, where the dissipation function is not known from first principles. Spatial structures produced in intermittent plasma influence dissipation, heating, and transport and acceleration of charged particles. Temporal intermittency can give rise to very long time correlations or a delayed approach to steady-state conditions, and has been associated with inverse cascade or quasi-inverse cascade systems, with possible implications for heliospheric prediction.

Intermittent Dissipation and Heating in 3D Kinetic Plasma Turbulence

High resolution, fully kinetic, three dimensional (3D) simulation of collisionless plasma turbulence shows the development of turbulence characterized by sheetlike current density structures spanning a range of scales. The nonlinear evolution is initialized with a long wavelength isotropic spectrum of fluctuations having polarizations transverse to an imposed mean magnetic field. We present evidence that these current sheet structures are sites for heating and dissipation, and that stronger currents signify higher dissipation rates. The analyses focus on quantities such as J·E, electron, and proton temperatures, and conditional averages of these quantities based on local electric current density. Evidently, kinetic scale plasma, like magnetohydrodynamics, becomes intermittent due to current sheet formation, leading to the expectation that heating and dissipation in astrophysical and space plasmas may be highly nonuniform. Comparison with previous results from 2D kinetic simulations, as well as high frequency solar wind observational data, are discussed.

Proton kinetic effects and turbulent energy cascade rate in the solar wind

The first observed connection between kinetic instabilities driven by proton temperature anisotropy and estimated energy cascade rates in the turbulent solar wind is reported using measurements from the Wind spacecraft at 1 AU. We find enhanced cascade rates are concentrated along the boundaries of the (β∥, T⊥/T∥) plane, which includes regions theoretically unstable to the mirror and firehose instabilities. A strong correlation is observed between the estimated cascade rate and kinetic effects such as temperature anisotropy and plasma heating, resulting in protons 5-6 times hotter and 70%-90% more anisotropic than under typical isotropic plasma conditions. These results offer new insights into kinetic processes in a turbulent regime.

Von Kármán energy decay and heating of protons and electrons in a kinetic turbulent plasma

Decay in time of undriven weakly collisional kinetic plasma turbulence in systems large compared to the ion kinetic scales is investigated using fully electromagnetic particle-in-cell simulations initiated with transverse flow and magnetic disturbances, constant density, and a strong guide field. The observed energy decay is consistent with the von Kármán hypothesis of similarity decay, in a formulation adapted to magnetohydrodyamics. Kinetic dissipation occurs at small scales, but the overall rate is apparently controlled by large scale dynamics. At small turbulence amplitudes the electrons are preferentially heated. At larger amplitudes proton heating is the dominant effect. In the solar wind and corona the protons are typically hotter, suggesting that these natural systems are in the large amplitude turbulence regime.

Inhomogeneous kinetic effects related to intermittent magnetic discontinuities

A connection between kinetic processes and two-dimensional intermittent plasma turbulence is observed using direct numerical simulations of a hybrid Vlasov-Maxwell model, in which the Vlasov equation is solved for protons, while the electrons are described as a massless fluid. During the development of turbulence, the proton distribution functions depart from the typical configuration of local thermodynamic equilibrium, displaying statistically significant non-Maxwellian features. In particular, temperature anisotropy and distortions are concentrated near coherent structures, generated as the result of the turbulent cascade, such as current sheets, which are nonuniformly distributed in space. Here, the partial variance of increments (PVI) method has been employed to identify high magnetic stress regions within a two-dimensional turbulent pattern. A quantitative association between non-Maxwellian features and coherent structures is established.

Kinetic signatures and intermittent turbulence in the solar wind plasma

A connection between kinetic processes and intermittent turbulence is observed in the solar wind plasma using measurements from the Wind spacecraft at 1 A.U. In particular, kinetic effects such as temperature anisotropy and plasma heating are concentrated near coherent structures, such as current sheets, which are nonuniformly distributed in space. Furthermore, these coherent structures are preferentially found in plasma unstable to the mirror and firehose instabilities. The inhomogeneous heating in these regions, which is present in both the magnetic field parallel and perpendicular temperature components, results in protons at least 3-4 times hotter than under typical stable plasma conditions. These results offer a new understanding of kinetic processes in a turbulent regime, where linear Vlasov theory is not sufficient to explain the inhomogeneous plasma dynamics operating near non-Gaussian structures.