When two collisionless plasma shells collide, they interpenetrate and the overlapping region may turn Weibel unstable for some values of the collision parameters. This instability grows magnetic filaments which, at saturation, have to block the incoming flow if a Weibel shock is to form. In a recent paper [J. Plasma Phys. (2016), vol. 82, 905820403], it was found implementing a toy model for the incoming particles trajectories in the filaments, that a strong enough external magnetic field $\mathbf{B}_0$ can prevent the filaments to block the flow if it is aligned with. Denoting $B_f$ the peak value of the field in the magnetic filaments, all test particles stream through them if $\alpha=B_0/B_f > 1/2$. Here, this result is extended to the case of an oblique external field $B_0$ making an angle $\theta$ with the flow. The result, numerically found, is simply $\alpha > \kappa(\theta)/\cos\theta$, where $\kappa(\theta)$ is of order unity. Noteworthily, test particles exhibit chaotic trajectories.

In the absence of frequent binary collisions to isotropize the plasma, the fulfillment of the magnetohydrodynamic (MHD) Rankine-Hugoniot jump conditions by collisionless shocks is not trivial. In particular, the presence of an external magnetic field can allow for stable anisotropies, implying some departures from the isotropic MHD jumps. The functional dependence of such anisotropies in terms of the field is yet to be determined. By hypothesizing a kinetic history of the plasma through the shock front, we recently devised a theory of the downstream anisotropy, hence of the density jump, in terms of the field strength for a parallel shock [J. Plasma Phys. (2018), vol. 84, 905840604]. Here, we extend the analysis to the case of a perpendicular shock. We still find that the field reduces the density jump, but the effect is less pronounced than in the parallel case.

Collisionless shocks follow the Rankine-Hugoniot jump conditions to a good approximation. However, for a shock propagating parallel to a magnetic field, magnetohydrodynamics states that the shock properties are independent of the field strength, whereas recent Particle-in-Cell simulations reveal a significant departure from magnetohydrodynamics behavior for such shocks in the collisionless regime. This departure is found to be caused by a field-driven anisotropy in the downstream pressure, but the functional dependence of this anisotropy on the field strength is yet to be determined. Here, we present a non-relativistic model of the plasma evolution through the shock front, allowing for a derivation of the downstream anisotropy in terms of the field strength. Our scenario assumes double adiabatic evolution of a pair plasma through the shock front. As a result, the perpendicular temperature is conserved. If the resulting downstream is firehose stable, then the plasma remains in this state. If unstable, it migrates towards the firehose stability threshold. In both cases, the conservation equations, together with the relevant hypothesis made on the temperature, allows a full determination of the downstream anisotropy in terms of the field strength.

When two plasmas collide, their interaction can be mediated by collisionless plasma instabilities or binary collisions between particles of each shell. By comparing the maximum growth rate of the collisionless instabilities with the collision frequency between particles of the shells, we determine the critical density separating the collisionless formation from the collisional formation of the resulting shock waves. This critical density is also the density beyond which the shock downstream is field free, as plasma instabilities do not have time to develop electromagnetic patterns. We further determine the conditions on the shells initial density and velocity for the downstream to be collisional. If these quantities fulfill the determined conditions, the collisionality of the downstream also prevents the shock from accelerating particles or generating strong magnetic fields. We compare the speed of sound with the relative speed of collision between the two shells, thus determining the portion of the parameters space where strong shock formation is possible for both classical and degenerate plasmas. Finally, we discuss the observational consequences in several astrophysical settings.

Energetic electromagnetic emissions by astrophysical jets like those that are launched during the collapse of a massive star and trigger gamma-ray bursts (GRBs) are partially attributed to relativistic internal shocks. The shocks are mediated in the collisionless plasma of such jets by the filamentation instability of counterstreaming particle beams. The filamentation instability grows fastest only if the beams move at a relativistic relative speed. We model here with a particle-in-cell (PIC) simulation the collision of two cold pair clouds at the speed c/2 (c: speed of light). We demonstrate that the two-stream instability outgrows the filamentation instability for this speed and is thus responsible for the shock formation. The incomplete thermalization of the upstream plasma by its quasi-electrostatic waves allows other instabilities to grow. A shock transition layer forms, in which a filamentation instability modulates the plasma far upstream of the shock. The inflowing upstream plasma is progressively heated by a two-stream instability closer to the shock and compressed to the expected downstream density by the Weibel instability. The strong magnetic field due to the latter is confined to a layer 10 electron skin depths wide.

The hierarchy of unstable modes when two counter-streaming pair plasmas interact over a flow-aligned magnetic field has been recently investigated [PoP \textbf{23}, 062122 (2016)]. The analysis is here extended to the case of an arbitrarily tilted magnetic field. The two plasma shells are initially cold and identical. For any angle $\theta \in [0,\pi/2]$ between the field and the initial flow, the hierarchy of unstable modes is numerically determined in terms of the initial Lorentz factor of the shells $\gamma_0$, and the field strength as measured by a parameter denoted $\sigma$. For $\theta=0$, four different kinds of mode are likely to lead the linear phase. The hierarchy simplifies for larger $\theta$'s, partly because the Weibel instability can no longer be cancelled in this regime. For $\theta>0.78$ (44$^\circ$) and in the relativistic regime, the Weibel instability always govern the interaction. In the non-relativistic regime, the hierarchy becomes $\theta$-independent because the interaction turns to be field-independent. As a result, the two-stream instability becomes the dominant one, regardless of the field obliquity.

According to magnetohydrodynamics (MHD), the encounter of two collisional magnetized plasmas at high velocity gives rise to shock waves. Investigations conducted so far have found that the same conclusion still holds in the case of collisionless plasmas. For the case of a flow-aligned field, MHD stipulates that the field and the fluid are disconnected, so that the shock produced is independent of the field. We present a violation of this MHD prediction when considering the encounter of two cold pair plasmas along a flow-aligned magnetic field. As the guiding magnetic field grows, isotropization is progressively suppressed, resulting in a strong influence of the field on the resulting structure. A micro-physics analysis allows to understand the mechanisms at work. Particle-in-cell simulations also support our conclusions and show that the results are not restricted to a strictly parallel field.

We examine with a particle-in-cell (PIC) simulation the collision of two equally dense clouds of cold pair plasma. The clouds interpenetrate until instabilities set in, which heat up the plasma and trigger the formation of a pair of shocks. The fastest-growing waves at the collision speed c/5 and low temperature are the electrostatic two-stream mode and the quasi-electrostatic oblique mode. Both waves grow and saturate via the formation of phase space vortices. The strong electric fields of these nonlinear plasma structures provide an efficient means of heating up and compressing the inflowing upstream leptons. The interaction of the hot leptons, which leak back into the upstream region, with the inflowing cool upstream leptons continuously drives electrostatic waves that mediate the shock. These waves heat up the inflowing upstream leptons primarily along the shock normal, which results in an anisotropic velocity distribution in the post-shock region. This distribution gives rise to the Weibel instability. Our simulation shows that even if the shock is mediated by quasi-electrostatic waves, strong magnetowaves will still develop in its downstream region.

For a Weibel shock to form, two plasma shells have to collide and trigger the Weibel instability. At saturation, this instability generates in the overlapping region magnetic filaments with peak field $B_f$. In the absence of an external guiding magnetic field, these filaments can block the incoming flow, initiating the shock formation, if their size is larger than the Larmor radius of the incoming particles in the peak field. Here we show that this results still holds in the presence of an external magnetic field $B_0$, provided it is not too high. Yet, for $B_0 \gtrsim B_f/2$, the filaments become unable to stop any particle, regardless of its initially velocity.

The instabilities triggered when two counter-streaming pair beams collide are analyzed. A guiding magnetic field is accounting for, while both beams are considered identical and cold. The instability analysis is conducted over the full \textbf{k}-spectrum, allowing to derive the hierarchy map of the dominant unstable modes, in terms of the initial beams energy $\gamma_0$ and a magnetic field strength parameter $\Omega_B$. Four different regions of the $(\Omega_B,\gamma_0)$ phase space are identified, each one governed by a different kind of mode. The analysis also unravels the existence of a "triple point", where 3 different modes grow exactly the same rate. A number of analytical expressions can be derived, either for the modes growth-rates, or for the frontiers between the 4 regions.

Collisionless shocks, that is shocks mediated by electromagnetic processes, are customary in space physics and in astrophysics. They are to be found in a great variety of objects and environments: magnetospheric and heliospheric shocks, supernova remnants, pulsar winds and their nebul\ae, active galactic nuclei, gamma-ray bursts and clusters of galaxies shock waves. Collisionless shock microphysics enters at different stages of shock formation, shock dynamics and particle energization and/or acceleration. It turns out that the shock phenomenon is a multi-scale non-linear problem in time and space. It is complexified by the impact due to high-energy cosmic rays in astrophysical environments. This review adresses the physics of shock formation, shock dynamics and particle acceleration based on a close examination of available multi-wavelength or in-situ observations, analytical and numerical developments. A particular emphasize is made on the different instabilities triggered during the shock formation and in association with particle acceleration processes with regards to the properties of the background upstream medium. It appears that among the most important parameters the background magnetic field through the magnetization and its obliquity is the dominant one. The shock velocity that can reach relativistic speeds has also a strong impact over the development of the micro-instabilities and the fate of particle acceleration. Recent developments of laboratory shock experiments has started to bring some new insights in the physics of space plasma and astrophysical shock waves. A special section is dedicated to new laser plasma experiments probing shock physics

Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical environments and they are mediated by the Weibel or filamentation instability. In such shocks, charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are characterized by a strong electrostatic field, which leads to electron trapping. Ions are accelerated by reflection from the electrostatic potential. Shock formation and particle acceleration will be discussed in theory and simulations.

The thin-shell instability has been named as one process, which can generate entangled structures in astrophysical plasma on collisional (fluid) scales. It is driven by a spatially varying imbalance between the ram pressure of the inflowing upstream plasma and the downstream's thermal pressure at a non-planar shock. Here we show by means of a particle-in-cell (PIC) simulation that an analogue process can destabilize a thin shell formed by two interpenetrating, unmagnetized and collisionless plasma clouds. The amplitude of the shell's spatial modulation grows and saturates after about ten inverse proton plasma frequencies, when the shell consists of connected piecewise linear patches.

The motion of a particle in a spatially harmonic magnetic field is a basic problem involved, for example, in the mechanism of formation of a collisionless shock. In such settings, it is generally reasoned that particles entering a Weibel generated turbulence are trapped inside it, provided their Larmor radius in the peak field is smaller than the field coherence length. The goal of this work is to put this heuristic conclusion on firm ground by studying, both analytically and numerically, such motion. A toy model is analyzed, consisting of a relativistic particle entering a region of space occupied by a spatially harmonic field. The particle penetrates the magnetic structure in a direction aligned with the magnetic filaments. Although the conclusions are not trivial, the main result is confirmed.

We analyse the full shock formation process in electron-ion plasmas in theory and simulations. It is accepted that electromagnetic shocks in initially unmagnetised relativistic plasmas are triggered by the filamentation instability. However, the transition from the first unstable phase to the quasi-steady shock is still missing. We derive a theoretical model for the shock formation time, taking into account the filament merging in the non-linear phase of the filamentation instability. This process is much slower than in electron-positron pair shocks, so that the shock formation is longer by a factor proportional to sqrt(m_i/m_e) ln(m_i/m_e).

In collisional fluids, a number of key processes rely on the frequency of binary collisions. Collisions seem necessary to generate a shock wave when two fluids collide fast enough, to fulfill the Rankine-Hugoniot relations, to establish an equation of state or a Maxwellian distribution. Yet, these seemingly collisional features are routinely either observed or assumed, in relation with collision\emph{less} astrophysical plasmas. This article will review our current answers to the following questions: How do colliding collisionless plasmas end-up generating a shock as if they were fluids? To which extent are the Rankine-Hugoniot relations fulfilled in this case? Do collisionless shocks propagate like fluid ones? Can we use an equation of state to describe collisionless plasmas, like MHD codes for astrophysics do? Why are Maxwellian distributions ubiquitous in Particle-In-Cell simulations of collisionless shocks? Time and length scales defining the border between the collisional and the collisionless behavior will be given when relevant. In general, when the time and length scales involved in the collisionless processes responsible for the fluid-like behavior may be neglected, the system may be treated like a fluid.

Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation of a collisionless shock. We unravel the physical mechanism at work and show that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter. The predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas. By predicting the shock formation time, experimental setups aiming at producing such shocks can be optimised to favourable conditions.

The filamentation (Weibel) instability plays a key role in the formation of collisionless shocks which are thought to produce Gamma-Ray-Bursts and High-Energy-Cosmic-Rays in astrophysical environments. While it has been known for long that a flow-aligned magnetic field can completely quench the instability, it was recently proved in 2D that in the cold regime, such cancelation is possible if and only if the field is perfectly aligned. Here, this result is finally extended to a 3D geometry. Calculations are conducted for symmetric and asymmetric counter-streaming relativistic plasma shells. 2D results are retrieved in 3D: the instability can never be completely canceled for an oblique magnetic field. In addition, the maximum growth-rate is always larger for wave vectors lying in the plan defined by the flow and the oblique field. On the one hand, this bears consequences on the orientation of the generated filaments. On the other hand, it certifies 2D simulations of the problem can be performed without missing the most unstable filamentation modes.

The filamentation instability triggered when two counter streaming plasma shells overlap appears to be the main mechanism by which collisionless shocks are generated. It has been known for long that a flow aligned magnetic field can completely suppress this instability. In a recent paper [PHYSICS OF PLASMAS 18, 080706 (2011)], it was demonstrated in two dimensions that for the case of two cold, symmetric, relativistically colliding shells, such cancellation cannot occur if the field is not perfectly aligned. Here, this result is extended to the case of two asymmetric shells. The filamentation instability appears therefore as an increasingly robust mechanism to generate shocks.

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.

We present results of 2D3V particle-in-cell simulations of non-relativistic plasma collisions with absent or parallel large-scale magnetic field for parameters applicable to the conditions at young supernova remnants. We study the collision of plasma slabs of different density, leading to two different shocks and a contact discontinuity. Electron dynamics play an important role in the development of the system. While non-relativistic shocks in both unmagnetized and magnetized plasmas can be mediated by Weibel-type instabilities, the efficiency of shock-formation processes is higher when a large-scale magnetic field is present. The electron distributions downstream of the forward and reverse shocks are generally isotropic, whereas that is not always the case for the ions. We do not see any significant evidence of pre-acceleration, neither in the electron population nor in the ion distribution.

Notes from a series of 13 one hour (or more) lectures on Plasma Physics given to Ramesh Narayan' research group at the Harvard-Smithsonian Center for Astrophysics, between January and July 2012. Lectures 1 to 5 cover various key Plasma Physics themes. Lectures 6 to 12 mainly go over the Review Paper on "Multidimensional electron beam-plasma instabilities in the relativistic regime" [\emph{Physics of Plasmas} \textbf{17}, 120501 (2010)]. Lectures 13 talks about the so-called Biermann battery and its ability to generate magnetic fields from scratch.

The growth and saturation of Buneman-type instabilities is examined with a particle-in-cell (PIC) simulation for parameters that are representative for the foreshock region of fast supernova remnant (SNR) shocks. A dense ion beam and the electrons correspond to the upstream plasma and a fast ion beam to the shock-reflected ions. The purpose of the 2D simulation is to identify the nonlinear saturation mechanisms, the electron heating and potential secondary instabilities that arise from anisotropic electron heating and result in the growth of magnetic fields. We confirm that the instabilities between both ion beams and the electrons saturate by the formation of phase space holes by the beam-aligned modes. The slower oblique modes accelerate some electrons, but they can not heat up the electrons significantly before they are trapped by the faster beam-aligned modes. Two circular electron velocity distributions develop, which are centred around the velocity of each ion beam. They develop due to the scattering of the electrons by the electrostatic wave potentials. The growth of magnetic fields is observed, but their amplitude remains low.

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here if the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here for that this electric field modification triggers a second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally.

There is an observational correlation between astrophysical shocks and non-thermal particle distributions extending to high energies. As a first step toward investigating the possible feedback of these particles on the shock at the microscopic level, we perform particle-in-cell (PIC) simulations of a simplified environment consisting of uniform, interpenetrating plasmas, both with and without an additional population of cosmic rays. We vary the relative density of the counterstreaming plasmas, the strength of a homogeneous parallel magnetic field, and the energy density in cosmic rays. We compare the early development of the unstable spectrum for selected configurations without cosmic rays to the growth rates predicted from linear theory, for assurance that the system is well represented by the PIC technique. Within the parameter space explored, we do not detect an unambiguous signature of any cosmic-ray-induced effects on the microscopic instabilities that govern the formation of a shock. We demonstrate that an overly coarse distribution of energetic particles can artificially alter the statistical noise that produces the perturbative seeds of instabilities, and that such effects can be mitigated by increasing the density of computational particles.

The filamentation instability (sometimes also referred to as "Weibel") is a key process in many astrophysical scenario. In the Fireball model for Gamma Ray Bursts, this instability is believed to mediate collisionless shock formation from the collision of two plasma shells. It has been known for long that a flow aligned magnetic field can completely cancel this instability. We show here that in the general case where there is an angle between the field and the flow, the filamentation instability can never be stabilized, regardless of the field strength. The presented model analyzes the stability of two symmetric counter-streaming cold electron/proton plasma shells. Relativistic effects are accounted for, and various exact analytical results are derived. This result guarantees the occurrence of the instability in realistic settings fulfilling the cold approximation.

Relativistic collisionless shocks are believed to be efficient particle accelerators. Nonlinear outcome of the interaction of accelerated particles that run ahead of the shock, the so-called "precursor", with the unperturbed plasma of the shock upstream, is thought to facilitate additional acceleration of these particles and to possibly modify the hydrodynamic structure of the shock. We explore here the linear growth of kinetic modes appearing in the precursor-upstream interaction in relativistic shocks propagating in non and weakly magnetized plasmas: electrostatic two-stream parallel mode and electrostatic oblique modes. These modes are of particular interest because they are the fastest growing modes known in this type of system. Using a simplified distribution function for a dilute ultra-relativistic beam that is relativistically hot in its own rest frame, yet has momenta that are narrowly collimated in the frame of the cold upstream plasma into which it propagates, we identify the fastest growing mode in the full $k$-space and calculate its growth rate. We consider all types of plasma (pairs and ions-electrons) and beam (charged and charge-neutral). We find that unstable electrostatic modes are present in any type of plasma and for any shock parameters. We further find that two modes, one parallel ($k_\perp=0$) and the other one oblique ($k_\perp \sim k_\|$), are competing for dominance and that either one may dominate the growth rate in different regions of the phase space. The dominant mode is determined mostly by the perpendicular spread of the accelerated particle momenta in the upstream frame, which reflects the shock Lorentz factor. The parallel mode becomes more dominant in shocks with lower Lorentz factors (i.e., with larger momentum spreads). We briefly discuss possible implications of our results for external shocks in gamma-ray burst sources.

The Rayleigh-Taylor instability is a key process in many fields of Physics ranging from astrophysics to inertial confinement fusion. It is usually analyzed deriving the linearized fluid equations, but the physics behind the instability is not always clear. Recent works on this instability allow for an very intuitive understanding of the phenomenon and for a straightforward calculation of the linear growth rate. In this Letter, it is shown that the same reasoning allows for a direct derivation of the relativistic expression of the linear growth rate for an incompressible fluid.

The prompt emissions of gamma-ray bursts are seeded by radiating ultrarelativistic electrons. Internal shocks propagating through a jet launched by a stellar implosion, are expected to amplify the magnetic field & accelerate electrons. We explore the effects of density asymmetry & a quasi-parallel magnetic field on the collision of plasma clouds. A 2D relativistic PIC simulation models the collision of two plasma clouds, in the presence of a quasi-parallel magnetic field. The cloud density ratio is 10. The densities of ions & electrons & the temperature of 131 keV are equal in each cloud. The mass ratio is 250. The peak Lorentz factor of the electrons is determined, along with the orientation & strength of the magnetic field at the cloud collision boundary. The magnetic field component orthogonal to the initial plasma flow direction is amplified to values that exceed those expected from shock compression by over an order of magnitude. The forming shock is quasi-perpendicular due to this amplification, caused by a current sheet which develops in response to the differing deflection of the incoming upstream electrons & ions. The electron deflection implies a charge separation of the upstream electrons & ions; the resulting electric field drags the electrons through the magnetic field, whereupon they acquire a relativistic mass comparable to the ions. We demonstrate how a magnetic field structure resembling the cross section of a flux tube grows in the current sheet of the shock transition layer. Plasma filamentation develops, as well as signatures of orthogonal magnetic field striping. Localized magnetic bubbles form. Energy equipartition between the ion, electron & magnetic energy is obtained at the shock transition layer. The electronic radiation can provide a seed photon population that can be energized by secondary processes (e.g. inverse Compton).

Magnetic-field generation by a relativistic ion beam propagating through an electron-ion plasma along a homogeneous magnetic field is investigated with 2.5D high-resolution particle-in-cell (PIC) simulations. The studies test predictions of a strong amplification of short-wavelength modes of magnetic turbulence upstream of nonrelativistic and relativistic parallel shocks associated with supernova remnants, jets of active galactic nuclei, and gamma-ray bursts. We find good agreement in the properties of the turbulence observed in our simulations compared with the dispersion relation calculated for linear waves with arbitrary orientation of ${\vec k}$. Depending on the parameters, the backreaction on the ion beam leads to filamentation of the ambient plasma and the beam, which in turn influences the properties of the magnetic turbulence. For mildly- and ultra-relativistic beams, the instability saturates at field amplitudes a few times larger than the homogeneous magnetic field strength. This result matches our recent studies of nonrelativistically drifting, hot cosmic-ray particles upstream of supernova-remnant shocks which indicated only a moderate magnetic-field amplification by nonresonant instabilities. We also demonstrate that the aperiodic turbulence generated by the beam can provide efficient particle scattering with a rate compatible with Bohm diffusion. Representing the ion beam as a constant external current, i.e. excluding a backreaction of the magnetic turbulence on the beam, we observe non-resonant parallel modes with wavelength and growth rate as predicted by analytic calculations. In this unrealistic setup the magnetic field is amplified to amplitudes far exceeding the homogeneous field, as observed in recent MHD and PIC simulations.

Many competing linear instabilities are likely to occur in astrophysical settings, and it is important to assess which one grows faster for a given situation. An analytical model including the main beam plasma instabilities is developed. The full 3D dielectric tensor is thus explained for a cold relativistic electron beam passing through a cold plasma, accounting for a guiding magnetic field, a return electronic current and moving protons. Considering any orientations of the wave vector allows to retrieve the most unstable mode for any parameters set. An unified description of the Filamentation (Weibel), Two-Stream, Buneman, Bell instabilities (and more) is thus provided, allowing for the exact determination of their hierarchy in terms of the system parameters. For relevance to both real situations and PIC simulations, the electron-to-proton mass ratio is treated as a parameter, and numerical calculations are conducted with two different values, namely 1/1836 and 1/100. In the system parameters phase space, the shape of the domains governed by each kind of instability is far from being trivial. For low density beams, the ultra-magnetized regime tends to be governed by either the Two-Stream or the Buneman instabilities. For beam densities equalling the plasma one, up to four kinds of modes are likely to play a role, depending of the beam Lorentz factor. In some regions of the system parameters phase space, the dominant mode may vary with the electron-to-proton mass ratio. Application is made to Solar Flares, Intergalactic Streams and Relativistic shocks physics.

Unstable modes growing when two plasma shells cross over a background plasma at arbitrary angle $\theta$, are investigated using a non-relativistic three cold fluids model. Parallel flows with $\theta=0$ are slightly more unstable than anti-parallel ones with $\theta=\pi$. The case $\theta=\pi/2$ is as unstable as the $\theta=0$ one, but the fastest growing modes are oblique. While the most unstable wave vector varies with orientation, its growth rate slightly evolves and there is no such thing as a stable configuration. A number of exact results can be derived, especially for the $\theta=\pi/2$ case.

Two charge- and current neutral plasma beams are modelled with a one-dimensional PIC simulation. The beams are uniform and unbounded. The relative speed between both beams is 0.4c. One beam is composed of electrons and protons and one out of protons and negatively charged oxygen (dust). All species have the temperature 9 keV. A Buneman instability develops between the electrons of the first beam and the protons of the second beam. The wave traps the electrons, which form plasmons. The plasmons couple energy into the ion acoustic waves, which trap the protons of the second beam. A proton phase space hole grows, which develops through its interaction with the oxygen and the heated electrons into a rarefaction pulse. This pulse drives a strong ion acoustic double layer, which accelerates a beam of electrons to about 50 MeV, which is comparable to the proton kinetic energy. The proton distribution eventually evolves into an electrostatic shock. Beams of charged particles moving at such speeds may occur in the foreshock of supernova remnant shocks. This double layer is thus potentially relevant for the electron acceleration (injection) into the diffusive shock acceleration by supernova remnants shocks.

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.

The growth rate of the filamentation instability triggered when a diluted cold electron beam passes through a cold plasma is evaluated using the quantum hydrodynamic equations. Compared with a cold fluid model, quantum effects reduce both the unstable wave vector domain and the maximum growth rate. Stabilization of large wave vector modes is always achieved, but significant reduction of the maximum growth rate depends on a dimensionless parameter that is provided. Although calculations are extended to the relativistic regime, they are mostly relevant to the non-relativistic one.

We present an investigation of the electromagnetic instabilities which are trig gered when an ultra relativistic electron beam passes through a plasma. The linear growth rate is computed for every direction of propagation of the unstable modes, and temperatures are modelled using simple waterbag distribution functions. The ultra relativistic unstable spectrum is located around a very narrow band centered on a critical angle which value is given analytically. The growth rate of modes propagating in this direction decreases like k^(-1/3).

The linear instability that induces a relativistic electron beam passing through a return plasma current to filament transversely is often related to some filamentation mode with wave vector normal to the beam or confused with Weibel modes. We show that these modes may not be relevant in this matter and identify the most unstable mode on the two-stream/filamentation branch as the main trigger for filamentation. This sets both the characteristic transverse and longitudinal filamentation scales in the non-resistive initial stage.