Search for: All records

A common occurrence in nearly collisionless plasmas is the unequal electron–ion temperatures. The pressure–strain interaction provides a mechanism-agnostic pathway for increasing plasma internal energy through spatiotemporally local isotropic compression and volume-preserving deformation, yet its behavior under thermal disequilibrium remains unexplored. We investigate this using five fully kinetic 2.5-dimensional particle-in-cell simulations of undriven decaying turbulence by varying the initial electron-to-ion temperature ratios. By analyzing the species' internal energy density alongside a decomposition of the pressure–strain interaction, with a focus on the volume-preserving deformation that contains normal and shear contributions, we quantify how the initial temperature imbalance modifies the channels through which turbulence increases each species' internal energy density. The cumulative pressure–strain interaction tracks the change in average internal energy density for both electrons and ions, where the total deformation channel dominates energy evolution. We discover that changes to electron internal energy density are governed primarily by the shear deformation power density, concentrated in electron-scale current sheets, while the ion shear and normal deformation components cancel, yielding a much smaller net deformation power density that peaks around, rather than within, those electron-scale current structures. By varying the initial temperature ratio, we find that the amplitudes and localization of deformation change, but preserve these qualitative trends. Together, these results show how thermal disequilibrium could shape species-dependent turbulent “heating rate,” measured via pressure–strain interaction and now approximated via only its shear deformation part, and provide a framework for interpreting energy evolution and conversion in turbulent space plasmas where unequal species temperature is the norm. 

Abstract Collisionless plasma shocks are a common feature of many space and astrophysical systems. They are sources of high-energy particles and nonthermal emission, channeling as much as 20% of the shock’s energy into nonthermal particles. The generation and acceleration of these nonthermal particles have been previously studied and shown to affect shock hydrodynamics to the zeroth order. In this work, we use self-consistent hybrid particle-in-cell simulations to examine the effect of self-generated nonthermal ion populations on the nature of collisionless, quasi-parallel shocks. Accelerated nonthermal particles downstream of the shock diffuse into the upstream region, taking energy away from the shock, which increases the compression ratio, slows the shock down, and flattens the nonthermal population’s spectral index for lower-Mach-number shocks. We show that this enhances shock compressibility when the heat flux is included in the Rankine–Hugoniot jump conditions, results that are roughly consistent with previous theories of “cosmic-ray-modified shocks.” Additionally, the simulation data show that heat flux and enthalpy flux cancels out in the upstream region, yielding a relatively simple, alternative closure for the jump conditions which accurately predict for the shock speed and compression ratio. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies in predictions of the arrival times of coronal mass ejections and the conflicting radio and X-ray observations of intracluster shocks. These effects will likely need to be included in fluid modeling to predict shock evolution accurately. 

We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection. 

Anisotropic electron heating during electron-only magnetic reconnection with a large guide magnetic field is directly measured in a laboratory plasma through in situ measurements of electron velocity distribution functions. Electron heating preferentially parallel to the magnetic field is localized to one separatrix, and anisotropies of 1.5 are measured. The mechanism for electron energization is identified as the parallel reconnection electric field because of the anisotropic nature of the heating and spatial localization. These characteristics are reproduced in a 2D particle-in-cell simulation and are also consistent with numerous magnetosheath observations. A measured increase in the perpendicular temperature along both separatrices is not reproduced by our 2D simulations. This work has implications for energy partition studies in magnetosheath and laboratory reconnection. 

Magnetic reconnection often initiates abruptly and then rapidly progresses to a nonlinear quasi-steady state. While satellites frequently detect reconnection events, ascertaining whether the system has achieved steady-state or is still evolving in time remains challenging. Here, we propose that the relatively rapid opening of the reconnection separatrices within the electron diffusion region serves as an indicator of the growth phase of reconnection. The opening of the separatrices is produced by electron flows diverging away from the neutral line downstream of the X-line and flowing around a dipolarization front. This flow pattern leads to characteristic spatial structures in the electron pressure-strain interaction that could be a useful indicator for the growth phase of a reconnection event. We employ two-dimensional particle-in-cell numerical simulations of anti-parallel magnetic reconnection to validate this prediction. We find that the signature discussed here, alongside traditional reconnection indicators, can serve as a marker of the growth phase. This signature is potentially accessible using multi-spacecraft single-point measurements, such as with NASA's Magnetospheric Multiscale satellites in Earth's magnetotail. Applications to other settings where reconnection occurs are also discussed. 

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.