Importance of accurate consideration of the electron inertia in hybrid-kinetic simulations of collisionless plasma turbulence: 1. The 2D limit

Neeraj Jain,Patricio Muñoz,Meisam Farzalipour Tabriz,Markus Rampp,Jörg Büchner

DOI: https://doi.org/10.1063/5.0087103

2022-02-02

Abstract:The dissipation mechanism of the magnetic energy in turbulent collisionless space and astrophysical plasmas is still not well understood. Its investigation requires efficient kinetic simulations of the energy transfer in collisionless plasma turbulence. In this respect, hybrid-kinetic simulations, in which ions are treated as particles and electrons as an inertial fluid, have begun to attract a significant interest recently. Hybrid-kinetic models describe both ion- and electron scale processes by ignoring electron kinetic effects so that they are computationally much less demanding compared to fully kinetic plasma models. Hybrid-kinetic codes solve either the Vlasov equation for the ions (Eulerian Vlasov-hybrid codes) or the equations of motion of the ions as macro-particles (Lagrangian Particle-in-Cell (PIC)-hybrid codes). They consider the inertia of the electron fluid using different approximations. We check the validity of these approximations by employing our recently massively parallelized three-dimensional PIC-hybrid code CHIEF which considers the electron inertia without any of the common approximations. In particular we report the results of simulations of two-dimensional collisionless plasma turbulence. We conclude that the simulation results obtained using hybrid-kinetic codes which use approximations to describe the electron inertia need to be interpreted with caution. We also discuss the parallel scalability of CHIEF, to the best of our knowledge, the first PIC-hybrid code which without approximations describes the inertial electron fluid.

Plasma Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: In the hybrid kinetic simulation of collisionless plasma turbulence, the accurate consideration of electron inertial effects is crucial for the understanding of magnetic energy dissipation mechanisms. Specifically, the paper evaluates the effectiveness of two commonly - used approximation methods (A1 and A2) in describing electron inertial effects and explores the impact of these approximation methods on simulation results at different scales.

### Background and Importance of the Problem

1. **Collisionless Plasma Turbulence**:
- Collisionless plasma turbulence is a key unsolved problem in plasma physics, involving scales from macroscopic fluid scales to microscopic kinetic scales.
- This type of turbulence is very important in laboratory high - temperature plasmas and astrophysical plasmas.

2. **Energy Dissipation Mechanisms**:
- The dissipation mechanism of magnetic energy in collisionless plasma turbulence has not been fully understood.
- In order to study this mechanism, efficient kinetic simulations are required.

3. **Hybrid Kinetic Simulation**:
- Hybrid kinetic simulation is a relatively low - cost computational method, in which ions are treated as particles and electrons are treated as fluids with inertia.
- This method ignores the kinetic effects of electrons but can cover processes from ion scales to electron scales.

### Main Contributions of the Paper

1. **Evaluating the Effectiveness of Approximation Methods**:
- The authors used their recently developed three - dimensional parallelized PIC - hybrid code CHIEF, which does not use any common electron inertia approximations.
- By comparing the results of CHIEF with the results using approximation methods, the accuracy of these approximation methods was evaluated.

2. **Approximation Methods A1 and A2**:
- **A1**: Ignore the unsteady electron inertia term \(\frac{m_e}{e} \frac{\partial \vec{u}_e}{\partial t}\) in the generalized Ohm's law (equation (3)).
- **A2**: Ignore the \(\nabla (\nabla \cdot \vec{E})\) term when solving the elliptic partial differential equation (equation (8)) for the electric field.

3. **Experimental Results**:
- **A1**: Introduces errors at the electron scale, especially in the current sheet structure.
- **A2**: Is invalid from the ion scale to the electron scale because \(\nabla (\nabla \cdot \vec{E})\) and \(\nabla^2 \vec{E}\) almost completely cancel each other.

### Conclusions

- **Main Conclusion**: The results of hybrid kinetic simulations using approximation methods A1 and A2 need to be carefully interpreted, especially for physical processes at the electron scale.
- **Physical Significance**: Approximation method A1 may miss key physical processes of collisionless dissipation, which mainly occur in current sheets at the electron scale. Approximation method A2 has effectiveness problems across the entire scale range.

### Formula Summary

1. **Generalized Ohm's Law**:
\[
\vec{E} = -\vec{u}_e \times \vec{B} - \frac{1}{en} \nabla p_e - \frac{m_e}{e} \left( \frac{\partial \vec{u}_e}{\partial t} + (\vec{u}_e \cdot \nabla) \vec{u}_e \right) + \eta \vec{j}
\]

2. **Elliptic Partial Differential Equation**:
\[
\nabla \times \nabla \times \vec{E} + \frac{1}{d_e^2} \left( 1 + \frac{m_e}{m_i} \right) \vec{E} = -\frac{1}{d_e^2} \left( \vec{u}_e + \frac{m_e}{m_i} \vec{u}_i \right) \times \vec{B} + \frac{1}{n e d_e^2} \left(