Abstract:The possibility of slow diffusion regions as the origin for extended TeV emission halos around some pulsars (such as PSR J0633+1746 and PSR B0656+14) challenges the standard scaling of the electron diffusion coefficient in the interstellar medium. Self-generated turbulence by electron-positron pairs streaming out of the pulsar wind nebula was proposed as a possible mechanism to produce the enhanced turbulence required to explain the morphology and brightness of these TeV halos. We perform fully kinetic 1D3V particle-in-cell simulations of this instability, considering the case where streaming electrons and positrons have the same density. This implies purely resonant instability as the beam does not carry any current. We compare the linear phase of the instability with analytical theory and find very reasonable agreement. The non-linear phase of the instability is also studied, which reveals that the intensity of saturated waves is consistent with a momentum exchange criterion between a decelerating beam and growing magnetic waves. With the adopted parameters, the instability-driven wavemodes cover both the Alfv\'enic (fluid) and kinetic scales. The spectrum of the produced waves is non-symmetric, with left-handed circular polarisation waves being strongly damped when entering the ion-cyclotron branch, while right-handed waves are suppressed at smaller wavelength when entering the Whistler branch. The low-wavenumber part of the spectrum remains symmetric when in the Alfv\'enic branch. As a result, positrons behave dynamically differently compared to electrons. The final drift velocity of positrons can maintain a larger value than the ambient Alfv\'en speed $V_A$ while the drift of electrons can drop below $V_A$. We also observed a second harmonic plasma emission in the wave spectrum. An MHD-PIC approach is warranted to probe hotter beams and investigate the Alfv\'en branch physics. We provide a few such test simulations to support this assertion. This work confirms that the self-confinement scenario develops essentially according to analytical expectations, but some of the adopted approximations (like the distribution of non-thermal particles in the beam) need to be revised and other complementary numerical techniques should be used to get closer to more realistic configuration.

3D cylindrical BGK model of electron phase-space holes with finite velocity and polarization drift

Theory of ion holes in plasmas with flat-topped electron distributions

Three-Dimensional Particle-In-Cell Simulations of Two-Dimensional Bernstein-Greene-Kruskal Modes

The dynamic of ion Bernstein-Greene-Kruskal holes in plasmas with regularized κ-distributed electrons

Magnetic Perturbations in Electron Phase-Space Holes: Contribution of Electron Polarization Drift

A Gyrokinetic Simulation Model for Low Frequency Electromagnetic Fluctuations in Magnetized Plasmas

Theory of electron and ion holes as vortices in the phase-space of collision-less plasmas

On Gyrokinetic-Fluid Model for Electromagnetic Fluctuations in Magnetized Plasmas

Magnetic field amplification in electron phase-space holes and related effects

A "Galactic Disk"-Model for Three-Dimensional Bernstein-Greene-Kruskal Modes in a Finite Magnetic Field

Electron vortex magnetic holes: a nonlinear coherent plasma structure

High-Resolution Particle-In-Cell Simulations of Two-Dimensional Bernstein-Greene-Kruskal Modes

An Unstructured Body-of-Revolution Electromagnetic Particle-in-Cell Algorithm with Radial Perfectly Matched Layers and Dual Polarizations

Electron Holes in a Regularized Kappa Background

Non-Maxwellian Ion Distribution in the Equatorial and Auroral Electrojets

Kinetic models of equilibrium structures in space plasmas

A gyrokinetic model for the plasma periphery of tokamak devices

Electron Acceleration by High Current-Density Relativistic Electron Bunch in Plasmas

Theory for the anomalous electron transport in Hall effect thrusters. I. Insights from particle-in-cell simulations

Nonlinear Velocity Redistribution Caused by Energetic-Particle-driven Geodesic Acoustic Modes, Mapped with the Beam-Plasma System

Kinetic simulations of electron-positron induced streaming instability in the context of gamma-ray halos around pulsars