We study the stability of current filaments produced by the Weibel, or current filamentation, instability in weakly magnetized counter-streaming plasmas. It is shown that a resonance exists between the current-carrying ions and a longitudinal drift-kink mode that strongly deforms and eventually breaks the current filaments. Analytical estimates of the wavelength, growth rate and saturation level of the resonant mode are derived and validated by three-dimensional particle-in-cell simulations. Furthermore, self-consistent simulations of counter-streaming plasmas indicate that this drift-kink mode is dominant in the slow down of the flows and in the isotropization of the magnetic field, playing an important role in the formation of collisionless shocks.

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