skip to main content


Title: Decomposition of Plasma Kinetic Entropy into Position and Velocity Space and the Use of Kinetic Entropy in Particle-in-Cell Simulations
We describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems. First, we calculate kinetic entropy in two forms – the “combinatorial” form related to the logarithm of the number of microstates per macrostate and the “continuous” form related to f ln f, where f is the particle distribution function. We discuss the advantages and disadvantages of each and discuss subtleties about implementing them in PIC codes. Using collisionless PIC simulations that are two-dimensional in position space and three-dimensional in velocity space, we verify the implementation of the kinetic entropy diagnostics and discuss how to optimize numerical parameters to ensure accurate results. We show the total kinetic entropy is conserved to three percent in an optimized simulation of anti-parallel magnetic reconnection. Kinetic entropy can be decomposed into a sum of a position space entropy and a velocity space entropy, and we use this to investigate the nature of kinetic entropy transport during collisionless reconnection. We find the velocity space entropy of both electrons and ions increases in time due to plasma heating during magnetic reconnection, while the position space entropy decreases due to plasma compression. This project uses collisionless simulations, so it cannot address physical dissipation mechanisms; nonetheless, the infrastructure developed here should be useful for studies of collisional or weakly collisional heliospheric, planetary, and astrophysical systems. Beyond reconnection, the diagnostic is expected to be applicable to plasma  more » « less
Award ID(s):
1804428 1602769
NSF-PAR ID:
10107473
Author(s) / Creator(s):
Date Published:
Journal Name:
Physics of plasmas
ISSN:
1070-664X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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
  2. Magnetic reconnection plays an important role in the release of magnetic energy and consequent energization of particles in collisionless plasmas. Energy transfer in collisionless magnetic reconnection is inherently a two-step process: reversible, collisionless energization of particles by the electric field, followed by collisional thermalization of that energy, leading to irreversible plasma heating. Gyrokinetic numerical simulations are used to explore the first step of electron energization, and we generate the first examples of field–particle correlation signatures of electron energization in 2D strong-guide-field collisionless magnetic reconnection. We determine these velocity space signatures at the x-point and in the exhaust, the regions of the reconnection geometry in which the electron energization primarily occurs. Modeling of these velocity–space signatures shows that, in the strong-guide-field limit, the energization of electrons occurs through bulk acceleration of the out-of-plane electron flow by the parallel electric field that drives the reconnection, a non-resonant mechanism of energization. We explore the variation of these velocity–space signatures over the plasma beta range 0.01≤βi≤1. Our analysis goes beyond the fluid picture of the plasma dynamics and exploits the kinetic features of electron energization in the exhaust region to propose a single-point diagnostic, which can potentially identify a reconnection exhaust region using spacecraft observations. 
    more » « less
  3. Magnetic reconnection plays an important role in the release of magnetic energy and consequent energization of particles in collisionless plasmas. Energy transfer in collisionless magnetic reconnection is inherently a two-step process: reversible, collisionless energization of particles by the electric field, followed by collisional thermalization of that energy, leading to irreversible plasma heating. Gyrokinetic numerical simulations are used to explore the first step of electron energization, and we generate the first examples of field-particle correlation (FPC) signatures of electron energization in 2D strong-guide-field collisionless magnetic reconnection. We determine these velocity space signatures at the x-point and in the exhaust, the regions of the reconnection geometry in which the electron energization primarily occurs. Modeling of these velocity-space signatures shows that, in the strong-guide-field limit, the energization of electrons occurs through bulk acceleration of the out-of-plane electron flow by parallel electric field that drives the reconnection, a non-resonant mechanism of energization. We explore the variation of these velocity-space signatures over the plasma beta range 0.01 < beta_i < 1. Our analysis goes beyond the fluid picture of the plasma dynamics and exploits the kinetic features of electron energization in the exhaust region to propose a single-point diagnostic which can potentially identify a reconnection exhaust region using spacecraft observations. 
    more » « less
  4. 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
  5. Abstract Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing nonthermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations, we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without a guide field. In X-points, where diverging flows result in a nondiagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain, we identify the mechanisms driving the nonideal electric field using a full Ohm law by means of a statistical analysis based on our PIC simulations. We show that the nonideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated nonuniform effective resistivity model that is negligible on global scales and becomes significant only locally in X-points. It captures the properties of collisionless reconnection with the aim of mimicking its essentials in nonideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission. 
    more » « less