skip to main content


Title: Kinetic entropy-based measures of distribution function non-Maxwellianity: theory and simulations
We investigate kinetic entropy-based measures of the non-Maxwellianity of distribution functions in plasmas, i.e. entropy-based measures of the departure of a local distribution function from an associated Maxwellian distribution function with the same density, bulk flow and temperature as the local distribution. First, we consider a form previously employed by Kaufmann & Paterson ( J. Geophys. Res. , vol. 114, 2009, A00D04), assessing its properties and deriving equivalent forms. To provide a quantitative understanding of it, we derive analytical expressions for three common non-Maxwellian plasma distribution functions. We show that there are undesirable features of this non-Maxwellianity measure including that it can diverge in various physical limits and elucidate the reason for the divergence. We then introduce a new kinetic entropy-based non-Maxwellianity measure based on the velocity-space kinetic entropy density, which has a meaningful physical interpretation and does not diverge. We use collisionless particle-in-cell simulations of two-dimensional anti-parallel magnetic reconnection to assess the kinetic entropy-based non-Maxwellianity measures. We show that regions of non-zero non-Maxwellianity are linked to kinetic processes occurring during magnetic reconnection. We also show the simulated non-Maxwellianity agrees reasonably well with predictions for distributions resembling those calculated analytically. These results can be important for applications, as non-Maxwellianity can be used to identify regions of kinetic-scale physics or increased dissipation in plasmas.  more » « less
Award ID(s):
1804428
NSF-PAR ID:
10280917
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Journal of Plasma Physics
Volume:
86
Issue:
5
ISSN:
0022-3778
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed. 
    more » « less
  2. Abstract

    In the diffusion region of magnetotail reconnection, particle distributions are highly structured, exhibiting triangular shapes and multiple striations that deviate dramatically from the Maxwellian distribution. Fully kinetic simulations have been demonstrated to be capable of producing the essential structures of the observed distribution functions, yet are computationally not feasible for 3D global simulations. The fluid models used for large‐scale simulations, on the other hand, do not have the kinetic physics necessary for describing reconnection accurately. Our study aims to bridge fully kinetic and fluid simulations by quantifying the information required to capture the non‐Maxwellian features in the distributions underlying the closures used in the fluid code. We compare the results of fully kinetic simulations with observed electron velocity distributions in a magnetotail reconnection diffusion region and use the maximum entropy model to reconstruct electron and ion distributions using various numbers of moments obtained from the simulation. Our results indicate that using only local moments, the maximum entropy model can reproduce many of the features of the distributions: (1) the electron outflow distribution with a tilted triangular structure is reproduced with 21 or more moments in agreement with Ng et al. (2018,https://doi.org/10.1063/1.5041758) and (2) counterstreaming distributions can be captured with the 35‐moment model when the separation in velocity space between the populations is large.

     
    more » « less
  3. Abstract

    We investigate the local proton energization at magnetic discontinuities/intermittent structures and the corresponding kinetic signatures in velocity phase space in Alfvénic (high cross helicity) and non-Alfvénic (low cross helicity) wind streams observed by Parker Solar Probe. By means of the partial variance of increments method, we find that the hottest proton populations are localized around compressible, coherent magnetic structures in both types of wind. Analysis of parallel and perpendicular temperature distributions suggest that the Alfvénic wind undergoes preferential enhancements ofTat such structures, whereas the non-Alfvénic wind experiences preferentialTenhancements. Although proton beams are present in both types of wind, the proton velocity distribution function displays distinct features. Hot beams, i.e., beams with beam-to-core perpendicular temperatureT⊥,b/T⊥,cup to three times larger than the total distribution anisotropy, are found in the non-Alfvénic wind, whereas colder beams are in the Alfvénic wind. Our data analysis is complemented by 2.5D hybrid simulations in different geometrical setups, which support the idea that proton beams in Alfvénic and non-Alfvénic wind have different kinetic properties and different origins. The development of a perpendicular nonlinear cascade, favored in balanced turbulence, allows a preferential relative enhancement of the perpendicular plasma temperature and the formation of hot beams. Cold field-aligned beams are instead favored by Alfvén wave steepening. Non-Maxwellian distribution functions are found near discontinuities and intermittent structures, pointing to the fact that the nonlinear formation of small-scale structures is intrinsically related to the development of highly nonthermal features in collisionless plasmas. Our results contribute to understanding the role of different coherent structures in proton energization and their implication in collisionless energy dissipation processes in space plasmas.

     
    more » « less
  4. Optical Thomson scattering is now a mature diagnostic tool for precisely measuring local plasma density and temperature. These measurements typically take advantage of a simplified analytical model of the scattered spectrum, which is built upon the assumption that each plasma species is in thermal equilibrium. However, this assumption fails for most laboratory plasmas of interest, which are often produced through high field ionization of atoms via ultrashort laser pulses and vulnerable to several kinetic instabilities. While it is possible to analytically model the Thomson scattered spectrum for some non-Maxwellian distribution functions, it is often not practical to do so for laboratory plasmas with highly complex and unstable distribution functions. We present a new method for predicting the Thomson scattered spectrum from any plasma directly from fully kinetic particle-in-cell simulations. This approach allows us to model the contributions of kinetic instabilities to the Thomson spectrum that aren’t taken into account in Maxwellian theory. We demonstrate this method’s capability to capture nonthermal features in the Thomson spectrum by simulating a simple bumpon- tail plasma as well as a more complex laser-ionized plasma. The versatility of this approach makes it an effective aid in the experimental design of Thomson diagnostics to directly characterize kinetic instabilities in laboratory plasmas. Index Terms—plasma measurement, low-temperature plasmas, plasma diagnostics, plasma simulation, plasma stability, plasma density, plasma temperature 
    more » « less
  5. Abstract

    Plasma escape from the high‐latitude ionosphere (ion outflow) serves as a significant source of heavy plasma to the magnetospheric plasma sheet and ring current regions. Outflows alter mass density and reconnection rates, hence global responses of the magnetosphere. A new fully kinetic and semi‐kinetic model, KAOS (Kinetic model of Auroral ion OutflowS), is constructed from first principles which traces large numbers of individual O+ion macro‐particles along curved magnetic field lines, using a guiding‐center approximation, in order to facilitate calculation of ion distribution functions and moments. Particle forces include mirror and parallel electric field forces, a self‐consistent ambipolar electric field, and a parameterized source of ion cyclotron resonance wave heating, thought to be central to the transverse energization of ions. The model is initiated with a steady‐state ion density altitude profile and Maxwellian velocity distribution and particle trajectories are advanced via a direct simulation Monte Carlo scheme. This outlines the implementation of the kinetic outflow model, demonstrates the model's ability to achieve near‐hydrostatic equilibrium necessary for simulation spin‐up, and investigates L‐shell dependent wave heating and pressure cookers scenarios. This paper illustrates the model initialization process and numerical investigations of L‐shell dependent outflows and pressure cooker environments and serves to advance our understanding of the drivers and particle dynamics in the auroral ionosphere.

     
    more » « less