Fully kinetic two-dimensional particle-in-cell simulations are used to study electron acceleration at high-Mach-number nonrelativistic perpendicular shocks. SNR shocks are mediated by the Weibel instability which is excited because of an interaction between shock-reflected and upstream ions. Nonlinear evolution of the Weibel instability leads to the formation of current sheets. At the turbulent shock ramp the current sheets decay through magnetic reconnection. The number of reconnection sites strongly depends on the ion-to-electron mass ratio and the Alfv\'enic Mach number of the simulated shock. Electron acceleration is observed at locations where magnetic reconnection operates. For the highest mass ratios almost all electrons are involved in magnetic reconnection, which makes the magnetic reconnection the dominant acceleration process for electrons at these shocks. We discuss the relevance of our results for 3D systems with realistic ion-to-electron mass ratio.

We revisit the effect of non-linear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in ${\bf k}$-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that non-linear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wavenumbers where wave driving is weak or absent. The relaxation time also strongly depends on the IGM temperature, $T_\mathrm{IGM}$, and for $T_\mathrm{IGM}\ll10$ eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small $k$, that non-linearly drains the wave energy at the resonant $\mathbf{k}$ of the pair-beam instability. Any other dissipation process operating at small $k$ would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other non-linear processes, such as the modulation instability, could provide additional damping of the non-resonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam.

This work revisits the electrostatic instability for blazar-induced pair beams propagating through IGM with the methods of linear analysis and PIC simulations. We study the impact of the realistic distribution function of pairs resulting from interaction of high-energy gamma-rays with the extragalactic background light. We present analytical and numerical calculations of the linear growth rate of the instability for arbitrary orientation of wave vectors. Our results explicitly demonstrate that the finite angular spread of the beam dramatically affects the growth rate of the waves, leading to fastest growth for wave vectors quasi-parallel to the beam direction and a growth rate at oblique directions that is only by a factor of 2-4 smaller compared to the maximum. To study the non-linear beam relaxation, we performed PIC simulations that take into account a realistic wide-energy distribution of beam particles. The parameters of the simulated beam-plasma system provide an adequate physical picture that can be extrapolated to realistic blazar-induced pairs. In our simulations the beam looses only 1\% percent of its energy, and we analytically estimate that the beam would lose its total energy over about $100$ simulation times. Analytical scaling is then used to extrapolate to the parameters of realistic blazar-induced pair beams. We find that they can dissipate their energy slightly faster by the electrostatic instability than through inverse-Compton scattering. The uncertainties arising from, e.g., details of the primary gamma-ray spectrum are too large to make firm statements for individual blazars, and an analysis based on their specific properties is required.

Pair beams produced by very high-energy radiation from TeV blazars emit gamma rays in the GeV band by inverse-Compton scattering of soft photons. The observed GeV-band signal is smaller than that expected from the full electromagnetic cascade. This means that the pair beams must be affected by other physical processes reducing their energy flux. One possible loss mechanism involves beam-plasma instabilities that we consider in the present work. For realistic parameters the pair beams can not be simulated by modern computers. Instead, we use a simple analytical model to find a range of the beam parameters that (i) provides a physical picture similar to that of realistic pair beams and (ii) at the same time can be handled by available computational resources. Afterwards, we performed corresponding 2D PIC simulations. We confirm that the beams experience only small changes in the relevant parameter regime, and other processes such as deflection in magnetic field must be at play.

In plasmas where the thermal energy density exceeds the magnetic energy density ($\beta_\parallel > 1$), the aperiodic ordinary mode (O-mode) instability is driven by an excess of parallel temperature $A = T_\perp /T_\parallel < 1$ (where $\parallel$ and $\perp$ denote directions relative to the uniform magnetic field). When stimulated by parallel plasma streams the instability conditions extend to low beta states, i.e., $\beta_\parallel <1$, and recent studies have proven the existence of a new regime, where the anisotropy threshold decreases steeply with lowering $\beta_\parallel \to 0$ if the streaming velocity is sufficiently high. However, the occurrence of this instability is questionable especially in the low-beta plasmas, where the electrostatic two-stream instabilities are expected to develop much faster in the process of relaxation of the counterstreams. It is therefore proposed here to identify the instability conditions for the O-mode below those required for the onset of the electrostatic instability. An hierarchy of these two instabilities is established for both the low $\beta_\parallel <1$ and large $\beta_\parallel > 1$ plasmas. The conditions where the O-mode instability can operate efficiently are markedly constrained by the electrostatic instabilities especially in the low-beta plasmas.