Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

TL;DR

This study uses particle-in-cell simulations to explore the microphysics of relativistic collisionless shocks loaded with electron-positron pairs, revealing how magnetization and pair-loading affect shock mediation, particle acceleration, and energy distribution.

Contribution

The paper provides new insights into the microphysics of relativistic shocks with pair-loading, including the transition between shock mediation mechanisms and effects on particle acceleration.

Findings

Shock mediation shifts from ion Larmor gyration to microturbulent scattering with increasing magnetization.

Post-shock pairs carry 20-50% of upstream ion energy, with electron energy inversely related to pair-loading.

Pair loading suppresses ion acceleration at low magnetizations, leading to thermal ions and nonthermal electrons.

Abstract

We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations $\sigma\lesssim 10^{-4}$ and pair-loading factors $Z_\pm \lesssim 10$ are studied, where $Z_\pm$ is the number of loaded electrons and positrons per ion. We find the following. (1) The shock becomes mediated by the ion Larmor gyration in the mean field when $\sigma$ exceeds a critical value $\sigma_{\rm L}$ that decreases with $Z_\pm$. At $\sigma\lesssim\sigma_{\rm L}$ the shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease with $Z_\pm$, leading to lower $\sigma_{\rm L}$. (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as $\overline{E}_{\rm e}\propto (Z_\pm+1)^{-1}$. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as $\sigma\approx 5\times 10^{-6}$. The ions then become essentially thermal with mean energy $\overline{E}_{\rm i}$, while electrons form a nonthermal tail, extending from $E\sim (Z_\pm + 1)^{-1}\overline{E}_{\rm i}$ to $\overline{E}_{\rm i}$. When $\sigma = 0$, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows.