Weakly magnetized, relativistic collisionless shock waves are not only the natural offsprings of relativistic jets in high-energy astrophysical sources, they are also associated with some of the most outstanding displays of energy dissipation through particle acceleration and radiation. Perhaps their most peculiar and exciting feature is that the magnetized turbulence that sustains the acceleration process, and (possibly) the secondary radiation itself, is self-excited by the accelerated particles themselves, so that the phenomenology of these shock waves hinges strongly on the microphysics of the shock. In this review, we draw a status report of this microphysics, benchmarking analytical arguments with particle-in-cell simulations, and extract consequences of direct interest to the phenomenology, regarding in particular the so-called microphysical parameters used in phenomenological studies.

We develop a comprehensive theoretical model of relativistic collisionless pair shocks mediated by the current filamentation instability. We notably characterize the noninertial frame in which this instability is of a mostly magnetic nature, and describe at a microscopic level the deceleration and heating of the incoming background plasma through its collisionless interaction with the electromagnetic turbulence. Our model compares well to large-scale 2D3V PIC simulations, and provides an important touchstone for the phenomenology of such plasma systems.

In this first paper of a series dedicated to the microphysics of unmagnetized, relativistic collisionless pair shocks, we discuss the physics of the Weibel-type transverse current filamentation instability (CFI) that develops in the shock precursor, through the interaction of an ultrarelativistic suprathermal particle beam with the background plasma. We introduce in particular the notion of "Weibel frame", or scattering center frame, in which the microturbulence is of mostly magnetic nature. We calculate the properties of this frame, using first a kinetic formulation of the linear phase of the instability, relying on Maxwell-J\"uttner distribution functions, then using a quasistatic model of the nonlinear stage of the instability. Both methods show that: (i) the Weibel frame moves at subrelativistic velocities relative to the background plasma, therefore at relativistic velocities relative to the shock front; (ii) the velocity of the Weibel frame relative to the background plasma scales with $\xi_{\rm b}$, i.e., the pressure of the suprathermal particle beam in units of the momentum flux density incoming into the shock; and (iii), the Weibel frame moves slightly less fast than the background plasma relative to the shock front. Our theoretical results are found to be in satisfactory agreement with the measurements carried out in dedicated large-scale 2D3V PIC simulations.

In this second paper of a series, we discuss the dynamics of a plasma entering the precursor of an unmagnetized, relativistic collisionless pair shock. We discuss how this background plasma is decelerated and heated through its interaction with a microturbulence that results from the growth of a current filamentation instability (CFI) in the shock precursor. We make use, in particular, of the reference frame $\mathcal R_{\rm w}$ in which the turbulence is mostly magnetic. This frame moves at relativistic velocities towards the shock front at rest, decelerating gradually from the far to the near precursor. In a first part, we construct a fluid model to derive the deceleration law of the background plasma expected from the scattering of suprathermal particles off the microturbulence. This law leads to the relationship $\gamma_{\rm p}\,\sim\,\xi_{\rm b}^{-1/2}$ between the background plasma Lorentz factor $\gamma_{\rm p}$ and the normalized pressure of the beam $\xi_{\rm b}$; it is found to match nicely the spatial profiles observed in large-scale 2D3V particle-in-cell simulations. In a second part, we model the dynamics of the background plasma at the kinetic level, incorporating the inertial effects associated with the deceleration of $\mathcal R_{\rm w}$ into a Vlasov-Fokker-Planck equation for pitch-angle diffusion. We show how the effective gravity in $\mathcal R_{\rm w}$ drives the background plasma particles through friction on the microturbulence, leading to efficient plasma heating. Finally, we compare a Monte Carlo simulation of our model with dedicated PIC simulations and conclude that it can satisfactorily reproduce both the heating and the deceleration of the background plasma in the shock precursor, thereby providing a successful 1D description of the shock transition at the microscopic level.

In this third paper of a series, we discuss the physics of the population of accelerated particles in the precursor of an unmagnetized, relativistic collisionless pair shock. In particular, we provide a theoretical estimate of their scattering length $l_{scatt}(p)$ in the self-generated electromagnetic turbulence, as well as an estimate of their distribution function. We obtain $l_{scatt}(p) \simeq (\gamma_p /\epsilon_B)(p/\gamma_{\infty} mc)^2 (c/\omega_p)$, with p the particle momentum in the rest frame of the shock front, $\epsilon_B$ the strength parameter of the microturbulence, $\gamma_p$ the Lorentz factor of the background plasma relative to the shock front and $\gamma_{\infty}$ its asymptotic value outside the precursor. We compare this scattering length to large-scale PIC simulations and find good agreement for the various dependencies.

The nonlinear evolution of current filaments generated by the Weibel-type filamentation instability is a topic of prime interest in space and laboratory plasma physics. In this paper, we investigate the stability of a stationary periodic chain of nonlinear current filaments in counterstreaming pair plasmas. We make use of a relativistic four-fluid model and apply the Floquet theory to compute the two-dimensional unstable eigenmodes of the spatially periodic system. We examine three different cases, characterized by various levels of nonlinearity and asymmetry between the plasma streams: a weakly nonlinear symmetric system, prone to purely transverse merging modes; a strongly nonlinear symmetric system, dominated by coherent drift-kink modes whose transverse periodicity is equal to, or an integer fraction of the unperturbed filaments; a moderately nonlinear asymmetric system, subject to a mix of kink and bunching-type perturbations. The growth rates and profiles of the numerically computed eigenmodes agree with particle-in-cell simulation results. In addition, we derive an analytic criterion for the transition between dominant filament-merging and drift-kink instabilites in symmetric two-beam systems.