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 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.

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 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.

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.

As a shock front interacts with turbulence, it develops corrugation which induces outgoing wave modes in the downstream plasma. For a fast shock wave, the incoming wave modes can either be fast magnetosonic waves originating from downstream, outrunning the shock, or eigenmodes of the upstream plasma drifting through the shock. Using linear perturbation theory in relativistic MHD, this paper provides a general analysis of the corrugation of relativistic magnetized fast shock waves resulting from their interaction with small amplitude disturbances. Transfer functions characterizing the linear response for each of the outgoing modes are calculated as a function of the magnetization of the upstream medium and as a function of the nature of the incoming wave. Interestingly, if the latter is an eigenmode of the upstream plasma, we find that there exists a resonance at which the (linear) response of the shock becomes large or even diverges. This result may have profound consequences on the phenomenology of astrophysical relativistic magnetized shock waves.

We present a predictive model of the nonlinear phase of the Weibel instability induced by two symmetric, counter-streaming ion beams in the non-relativistic regime. This self-consistent model combines the quasilinear kinetic theory of Davidson et al. [Phys. Fluids 15, 317 (1972)] with a simple description of current filament coalescence. It allows us to follow the evolution of the ion parameters up to a stage close to complete isotropization, and is thus of prime interest to understand the dynamics of collisionless shock formation. Its predictions are supported by 2-D and 3-D particle-in-cell simulations of the ion Weibel instability. The derived approximate analytical solutions reveal the various dependencies of the ion relaxation to isotropy. In particular, it is found that the influence of the electron screening can affect the results of simulations using an unphysical electron mass.

We present a particle-in-cell simulation of the generation of a collisionless turbulent shock in a dense plasma driven by an ultra-high-intensity laser pulse. From the linear analysis, we highlight the crucial role of the laser-heated and return-current electrons in triggering a strong Weibel-like instability, giving rise to a magnetic turbulence able to isotropize the target ions.

We report here on a fast current-driven instability at relativistic collisionless shocks, triggered by the perpendicular current carried by the supra-thermal particles as they gyrate around the background magnetic field in the shock precursor. We show that this instability grows faster than any other instability studied so far in this context, and we argue that it is likely to shape the physics of the shock and of particle acceleration in a broad parameter range.

The physics of collisionless relativistic shocks with a moderate magnetization is presented. Micro-physics is relevant to explain the most energetic radiative phenomena of Nature, namely that of the termination shock of Gamma Ray Bursts. A transition towards Fermi process occurs for decreasing magnetization around a critical value which turns out to be the condition for the scattering to break the mean field inhibition. Scattering is produced by magnetic micro-turbulence driven by the current carried by returning particles, which had not been considered till now, but turns out to be more intense than Weibel's one around the transition. The current is also responsible for a buffer effect on the motion of the incoming flow, on which the threshold for the onset of turbulence depends.

The physics of instabilities in the precursor of relativistic collisionless shocks is of broad importance in high energy astrophysics, because these instabilities build up the shock, control the particle acceleration process and generate the magnetic fields in which the accelerated particles radiate. Two crucial parameters control the micro-physics of these shocks: the magnetization of the ambient medium and the Lorentz factor of the shock front; as of today, much of this parameter space remains to be explored. In the present paper, we report on a new instability upstream of electron-positron relativistic shocks and we argue that this instability shapes the micro-physics at moderate magnetization levels and/or large Lorentz factors. This instability is seeded by the electric current carried by the accelerated particles in the shock precursor as they gyrate around the background magnetic field. The compensation current induced in the background plasma leads to an unstable configuration, with the appearance of charge neutral filaments carrying a current of the same polarity, oriented along the perpendicular current. This ``current-driven filamentation'' instability grows faster than any other instability studied so far upstream of relativistic shocks, with a growth rate comparable to the plasma frequency. Furthermore, the compensation of the current is associated with a slow-down of the ambient plasma as it penetrates the shock precursor (as viewed in the shock rest frame). This slow-down of the plasma implies that the ``current driven filamentation'' instability can grow for any value of the shock Lorentz factor, provided the magnetization \sigma <~ 10^{-2}. We argue that this instability explains the results of recent particle-in-cell simulations in the mildly magnetized regime.

We present an exact calculation of the power spectrum of the electromagnetic fluctuations in a relativistic equilibrium plasma described by Maxwell-J\"uttner distribution functions. We consider the cases of wave vectors parallel or normal to the plasma mean velocity. The relative contributions of the subluminal and supraluminal fluctuations are evaluated. Analytical expressions of the spatial fluctuation spectra are derived in each case. These theoretical results are compared to particle-in-cell simulations, showing a good reproduction of the subluminal fluctuation spectra.

Collisionless shocks are ubiquitous in astrophysics and in the lab. Recent numerical simulations and experiments have shown how they can arise from the encounter of two collisionless plasma shells. When the shells interpenetrate, the overlapping region turns unstable, triggering the shock formation. As a first step towards a microscopic understanding of the process, we analyze here in detail the initial instability phase. On the one hand, 2D relativistic PIC simulations are performed where two symmetric initially cold pair plasmas collide. On the other hand, the instabilities at work are analyzed, as well as the field at saturation and the seed field which gets amplified. For mildly relativistic motions and onward, Weibel modes govern the linear phase. We derive an expression for the duration of the linear phase in good agreement with the simulations. This saturation time constitutes indeed a lower-bound for the shock formation time.

Besides being one of the most fundamental basic issues of plasma physics, the stability analysis of an electron beam-plasma system is of critical relevance in many areas of physics. Surprisingly, decades of extensive investigation had not yet resulted in a realistic unified picture of the multidimensional unstable spectrum within a fully relativistic and kinetic framework. All attempts made so far in this direction were indeed restricted to simplistic distribution functions and/or did not aim at a complete mapping of the beam-plasma parameter space. The present paper comprehensively tackles this problem by implementing an exact linear model. We show that three kinds of modes compete in the linear phase, which can be classified according to the direction of their wavenumber with respect to the beam. We then determine their respective domain of preponderance in a three-dimensional parameter space. All these results are supported by multidimensional particle-in-cell simulations.