Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

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.