Search

Shockwaves in plasma are usually dealt with using Magnetohydrodynamics (MHD). Yet, MHD entails the assumption of a short mean free path, which is not fulfilled in a collisionless plasma. Recently, for pair plasmas, we devised a model allowing to account for kinetic effects within an MHD-like formalism. Its relies on an estimate of the anisotropy generated when crossing the front, with a subsequent assessment of the stability of this anisotropy in the downstream. We solved our model for parallel, perpendicular and switch-on shocks. Here we bridge between all these cases by treating the problem of an arbitrarily, but coplanar, oriented magnetic field. Even though the formalism presented is valid for anisotropic upstream temperatures, only the case of a cold upstream is solved. We find extra solutions which are not part of the MHD catalog, and a density jump that is notably less in the quasi parallel, highly magnetized, regime. Given the complexity of the calculations, this work is mainly devoted to the presentation of the mathematical aspect of our model. A forthcoming article will be devoted to the physics of the shocks here defined., 16 pages, 8 figures, to appear in Journal of Plasma Physics

Shocks waves are a ubiquitous feature of many astrophysical plasma systems, and an important process for energy dissipation and transfer. The physics of these shock waves are frequently treated/modeled as a collisional, fluid MHD discontinuity, despite the fact that many shocks occur in the collisionless regime. In light of this, using fully kinetic, 3D simulations of non-relativistic, parallel propagating collisionless shocks comprised of electron-positron plasma, we detail the deviation of collisionless shocks form MHD predictions for varying magnetization/Alfv��nic Mach numbers, with particular focus on systems with Alf��nic Mach numbers much smaller than sonic Mach numbers. We show that the shock compression ratio decreases for sufficiently large upstream magnetic fields, in agreement with the predictions of Bret & Narayan (2018). Additionally, we examine the role of magnetic field strength on the shock front width. This work reinforces a growing body of work that suggest that modeling many astrophysical systems with only a fluid plasma description omits potentially important physics., 7 pages, 5 figures + a diagram, submitted to MNRAS

The properties of collisionless shocks, like the density jump, are usually derived from magnetohydrodynamics (MHD), where isotropic pressures are assumed. Yet, in a collisionless plasma, an external magnetic field can sustain a stable anisotropy. In \cite{BretJPP2018}, we devised a model for the kinetic history of the plasma through the shock front, allowing to self-consistently compute the downstream anisotropy, hence the density jump, in terms of the upstream parameters. This model dealt with the case of a parallel shock, where the magnetic field is normal to the front both in the upstream and the downstream. Yet, MHD also allows for shock solutions, the so-called switch-on solutions, where the field is normal to the front only in the upstream. This article consists in applying our model to these switch-on shocks. While MHD offers only 1 switch-on solution within a limited range of Alfv\'{e}n Mach numbers, our model offers 2 kinds of solutions within a slightly different range of Alfv\'{e}n Mach numbers. These 2 solutions are most likely the outcome of the intermediate and fast MHD shocks under our model. While the intermediate and fast shocks merge in MHD for the parallel case, they do not within our model. For simplicity, the formalism is restricted to non-relativistic shocks in pair plasmas where the upstream is cold., Comment: 18 pages. 6 figures. To appear in Journal of Plasma Physics

In a recent work [Bret, EPL \textbf{135} (2021) 35001], quantum electrodynamic (QED) effects were evaluated for the two-stream instability. It pertains to the growth of perturbations with a wave vector oriented along the flow in a collisionless counter-streaming system. Here, the analysis is extended to every possible orientation of the wave vector. The previous result for the two-stream instability is recovered, and it is proved that even within the framework of a 3D analysis, this instability remains fundamentally 1D even when accounting for QED effects. The filamentation instability, found for wave vectors normal to the flow, is weakly affected by QED corrections. As in the classical case, its growth rate saturates at large $k_\perp$. The saturation value is found independent of QED corrections. Also, the smallest unstable $k_\perp$ is independent of QED corrections. Surprisingly, unstable modes found for oblique wave vectors do \emph{not} follow the same pattern. For some, QED corrections do reduce the growth rate. But for others, the same corrections increase the growth rate instead. The possibility for QED effects to play a role in un-magnetized systems is evaluated. Pair production resulting from gamma emission by particles oscillating in the exponentially growing fields, is not accounting for., Comment: 19 pages, 7 figures, to appear in Physical Review E

Within the context of Magnetohydrodynamics (MHD), the properties of a parallel shock do not depend on the field strength, as the field and the fluid are disconnected for such a geometry. However, in the collisionless case, the field can sustain a stable anisotropy in the downstream, triggering a departure from the expected MHD behavior. In a recent work [A. Bret and R. Narayan, J. Plasma Phys. \textbf{84}, 905840604 (2018)], a theoretical model was presented allowing to derive the density ratio of a non-relativistic parallel collisionless shock in an electron/positron plasma, as a function of the field. Here we derive the entropy, pressure and temperature jumps stemming from this model. It is found to offer a transition between a 3D and a 1D downstream for the jumps in density, entropy, parallel temperature and parallel pressure., 6 pages, 4 figure. To appear in Physics of Plasmas

In shockwave theory, the density, velocity and pressure jumps are derived from the conservation equations. Here, we address the physics of a weak shock the other way around. We first show that the density profile of a weak shockwave in a fluid can be expressed as a sum of linear acoustic modes. The shock so built propagates at the speed of sound and matter is exactly conserved at the front crossing. Yet, momentum and energy are only conserved up to order 0 in powers of the shock amplitude. The density, velocity and pressure jumps are similar to those of a fluid shock, and an equivalent Mach number can be defined. A similar process is possible in magnetohydrodynamic. Yet, such a decomposition is found impossible for collisionless shocks due to the dispersive nature of ion acoustic waves. Weakly nonlinear corrections to their frequency do not solve the problem. Weak collisionless shocks could be inherently nonlinear, non-amenable to any linear superposition. Or they could be nonexistent, as hinted by recent works., Comment: 8 pages, 2 figures, to appear in Journal of Plasma Physics

When applied to compute the density jump of a shock, the standard magnetohydrodynamic (MHD) formalism assumes, 1) that all the upstream material passes downstream, together with the momentum and energy it carries, and 2) that pressures are isotropic. In a collisionless shock, shock accelerated particles going back and forth around the front can invalid the first assumption. In addition, an external magnetic field can sustain stable pressure anisotropies, invaliding the second assumption. It is therefore unclear whether the density jump of a collisionless shock fulfils the MHD jump or not. Here we try to clarify this issue. A literature review is conducted on 68 articles dealing with Particle-In-Cell simulations of collisionless shocks. We analyze the factors triggering departure from the MHD density jump and quantify their influence on $\Delta_{RH}$, the relative departure from the Rankine-Hugoniot jump. For small departures we propose $\Delta_{RH} = + \mathcal{O}(10^{-1-3.7\kappa})t^\kappa - \sigma \mathcal{O}(1)$ where $t$ is the timescale of the simulation, $\sigma$ the magnetization parameter and $\kappa$ a constant of order unity. The first term stems from the energy leakage into accelerated particle. The second term stems from the downstream anisotropy triggered by the field (assuming an isotropic upstream). This relation allows to assess to which extent a collisionless shock fulfils the RH density jump. In the strong field limit and for parallel shocks, the departure caused by the field saturates at a finite, negative, value. For perpendicular shocks, the departure goes to zero at small and high $\sigma$'s so that we find here a departure window. The results obtained have to be checked against full 3D simulations., Comment: To appear to ApJ. 13 pages, 5 figures

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., Comment: 10 pages, 6 figures, to appear in Journal of Plasma Physics

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., Comment: 17 pages, 9 figures, to appear in Physics of Plasmas

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., Comment: To appear in Journal of Plasma Physics, 16 pages, 7 figures

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., Comment: To appearin Journal of Plasma Physics - 35 pages, 11 figures

The expansion of a radial blast shell into an ambient plasma is modeled with a particle-in-cell (PIC) simulation. The unmagnetized plasma consists of electrons and protons. The formation and evolution of an electrostatic shock is observed, which is trailed by ion-acoustic solitary waves that grow on the beam of the blast shell ions in the post-shock plasma. In spite of the initially radially symmetric outflow, the solitary waves become twisted and entangled and, hence, they break the radial symmetry of the flow. The waves and their interaction with the shocked ambient ions slows down the blast shell protons and brings the post-shock plasma closer to an equilibrium., Comment: 5 pages, 4 figures, accepted by Physics of Plasmas

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., Comment: To appear in Physics of Plasmas

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., Comment: 17 pages, 9 figures, accepted for publication by the Journal of Plasma Physics

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., Comment: 11 pages, 6 figures. To appear in Journal of Plasma Physics

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., To Appear in ApJ

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., Comment: To appear in Physics of Plasmas

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, Comment: This is an author-created, un-copyedited version of an article accepted for publication/published in Reports on progress in physics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at 10.1088/0034-4885/79/4/046901

We present a \emph{Mathematica} notebook allowing for the symbolic calculation of the $3\times3$ dielectric tensor of a electron-beam plasma system in the fluid approximation. Calculation is detailed for a cold relativistic electron beam entering a cold magnetized plasma, and for arbitrarily oriented wave vectors. We show how one can elaborate on this example to account for temperatures, arbitrarily oriented magnetic field or a different kind of plasma., To appear in "Computer Physics Communications"

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. info:eu-repo/semantics/acceptedVersion

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., Comment: 5 pages, 6 figures

We analyze the full shock formation process in electron-ion plasmas in theory and simulations. It is accepted that electromagnetic shocks in initially unmagnetized relativistic plasmas are triggered by 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 filament merging in the nonlinear phase of the filamentation instability. This process is much slower than in electron-positron pair shocks, and so the shock formation is longer by a factor proportional to root m(i)/m(e) ln (m(i)/m(e)). info:eu-repo/semantics/publishedVersion

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., Comment: 20 pages, 4 figures, to appear in the Lab Astrophysics Special Issue of "Journal of Plasma Physics"

The expansion of a dense plasma into a dilute plasma across an initially uniform perpendicular magnetic field is followed with a one-dimensional particle-in-cell (PIC) simulation over MHD time scales. The dense plasma expands in the form of a fast rarefaction wave. The accelerated dilute plasma becomes separated from the dense plasma by a tangential discontinuity at its back. A fast magnetosonic shock with the Mach number 1.5 forms at its front. Our simulation demonstrates how wave dispersion widens the shock transition layer into a train of nonlinear fast magnetosonic waves., Comment: 5 pages, 4 figures, accepted for publication in Physics of Plasmas

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., Comment: To appear in Physics of Plasmas

Collisionless shocks are pervasive in astrophysics and they are critical to understand cosmic ray acceleration. Laboratory experiments with intense lasers are now opening the way to explore and characterise the underlying microphysics, which determine the acceleration process of collisionless shocks. We determine the shock character – electrostatic or electromagnetic – based on the stability of electrostatic shocks to transverse electromagnetic fluctuations as a function of the electron temperature and flow velocity of the plasma components and we compare the analytical model with particle-in-cell simulations. By making the connection with the laser parameters driving the plasma flows, we demonstrate that shocks with different and distinct underlying microphysics can be explored in the laboratory with state-of-the-art laser systems.

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., 3 pages, 2 figures, to appear in Physics of Plasmas

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., Comment: 10pages, 7 figures, accepted for publication in Physics of Plasmas

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., Comment: 22 pages, 18 figures

Cylindrical implosions driven by intense heavy ions beams should be instrumental in a near future to study High Energy Density Matter. By rotating the beam by means of a high frequency wobbler, it should be possible to deposit energy in the outer layers of a cylinder, compressing the material deposited in its core. The beam temporal profile should however generate an inevitable irradiation asymmetry likely to feed the Rayleigh-Taylor instability (RTI) during the implosion phase. In this paper, we compute the Fourier components of the target irradiation in order to make the junction with previous works on RTI performed in this setting. Implementing a 1D and 2D beam models, we find these components can be expressed exactly in terms of the Fourier transform of the temporal beam profile. If $T$ is the beam duration and $\Omega$ its rotation frequency, "magic products" $\Omega T$ can be identified which cancel the first harmonic of the deposited density, resulting in an improved irradiation symmetry., Comment: 19 pages, 8 figures, to appear in PRE

The Buneman instability occurring when an electron population is drifting with respect to the ions is analyzed in the quantum linear and nonlinear regimes. The one-dimensional low-frequency and collisional model of Shokri and Niknam [Phys. Plasmas, v. 12, p. 062110 (2005)] is revisited introducing the Bohm potential term in the momentum equation. The linear regime is investigated analytically, and quantum effects result in a reduction of the instability. The nonlinear regime is then assessed both numerically and analytically, and pure quantum density oscillations are found to appear during the late evolution of the instability.

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., Comment: To appear in Physics of Plasmas Letters

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)., 12 pages, 15 Figures, accepted to A&A

Particle-in-cell (PIC) simulations are widely used as a tool to investigate instabilities that develop between a collisionless plasma and beams of charged particles. However, even on contemporary supercomputers, it is not always possible to resolve the ion dynamics in more than one spatial dimension with such simulations. The ion mass is thus reduced below 1836 electron masses, which can affect the plasma dynamics during the initial exponential growth phase of the instability and during the subsequent nonlinear saturation. The goal of this article is to assess how far the electron to ion mass ratio can be increased, without changing qualitatively the physics. It is first demonstrated that there can be no exact similarity law, which balances a change of the mass ratio with that of another plasma parameter, leaving the physics unchanged. Restricting then the analysis to the linear phase, a criterion allowing to define a maximum ratio is explicated in terms of the hierarchy of the linear unstable modes. The criterion is applied to the case of a relativistic electron beam crossing an unmagnetized electron-ion plasma., Comment: To appear in Physics of Plasmas

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., 14 pages, 8 figures, accepted to ApJ

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., Comment: 4 pages, 3 figures, to appear in Phys. Lett. A

The filamentation instability occurring when a non relativistic electron beam passes through a quantum magnetized plasma is investigated by means of a cold quantum magnetohydrodynamic model. It is proved that the instability can be completely suppressed by quantum effects if and only if a finite magnetic field is present. A dimensionless parameter is identified which measures the strength of quantum effects. Strong quantum effects allow for a much smaller magnetic field to suppress the instability than in the classical regime., 9 pages, 5 figures, to appear in Phys. Plasmas

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., To appear in PRL

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., Comment: To appear in Physics of Plasmas

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)., Comment: 5 pages, 3 figures, to appear in EuroPhysics Letters

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., 4 page, 3 figures, to appear in PRL

We consider the collective neutrino plasma interactions and study the electron plasma instabilities produced by a nearly mono-energetic neutrino beam in a plasma. We describe the mutual interaction between neutrino flavor oscillations and electron plasma waves. We show that the neutrino flavor oscillations are not only perturbed by electron plasmas waves but also contribute to the dispersion relation and the growth rates of neutrino beam instabilities.

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.

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. The physics of the parallel and oblique modes is similar, and thus, we refer to the entire spectrum of electrostatic modes as 'two-stream-like'. 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 ultrarelativistic 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 electrostaticmore » modes are present in any type of plasma and for any shock parameters. We further find that two modes, one parallel (k{sub perpendicular} = 0) and the other one oblique (k{sub perpendicular} {approx} k{sub ||}), 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.« less

The quantum electromagnetic dielectric tensor for a multi species plasma is re-derived from the gauge invariant Wigner-Maxwell system and presented under a form very similar to the classical one. The resulting expression is then applied to a quantum kinetic theory of the electromagnetic filamentation instability. Comparison is made with the quantum fluid theory including a Bohm pressure term, and with the cold classical plasma result. A number of analytical expressions are derived for the cutoff wave vector, the largest growth rate and the most unstable wave vector.

The stability of two quantum counter-streaming electron beams is investigated within the quantum plasma fluid equations for arbitrarily oriented wave vectors. The analysis reveals that the two quantum two-stream unstable branches are indeed connected by a continuum of unstable modes with oblique wave vectors. Using the longitudinal approximation, the stability domain for any k is analytically explained, together with the growth rate.

The filamentation instability (FI) of counter-propagating beams of electrons is modelled with a particle-in-cell simulation in one spatial dimension and with a high statistical plasma representation. The simulation direction is orthogonal to the beam velocity vector. Both electron beams have initially equal densities, temperatures and moduli of their nonrelativistic mean velocities. The FI is electromagnetic in this case. A previous study of a small filament demonstrated, that the magnetic pressure gradient force (MPGF) results in a nonlinearly driven electrostatic field. The probably small contribution of the thermal pressure gradient to the force balance implied, that the electrostatic field performed undamped oscillations around a background electric field. Here we consider larger filaments, which reach a stronger electrostatic potential when they saturate. The electron heating is enhanced and electrostatic electron phase space holes form. The competition of several smaller filaments, which grow simultaneously with the large filament, also perturbs the balance between the electrostatic and magnetic fields. The oscillations are damped but the final electric field amplitude is still determined by the MPGF., Comment: 14 pages, 10 plots, accepted for publication in Physica Scripta

Proton beam transport in the context of proton driven Fast Ignition is usually assumed to be stable due to protons high inertia, but an analytical analysis of the process is still lacking. The stability of a charge and current neutralized proton beam passing through a plasma is therefore conducted here, for typical proton driven Fast Ignition parameters. In the cold regime, two fast growing Buneman-like modes are found, with an inverse growth-rate much smaller than the beam time-of-flight to the target core. The stability issue is thus not so obvious, and Kinetic effects are investigated. One unstable mode is found stabilized by the background plasma protons and electrons temperatures. The second mode is also damped, providing the proton beam thermal spread is larger than $\sim$ 10 keV. In Fusion conditions, the beam propagation should therefore be stable., Submitted to PoP

Quantum and collisional effects on the unstable spectrum of a relativistic electron beam-plasma system are investigated through a two-fluids model. Application is made to the near target center interaction of the relativistic electron beam in the Fast Ignition Scenario. Partial degeneracy effects are found negligible while the most influential factors are the beam temperature and the electron-ion collision frequency of the plasma. The introduction of the latter triggers some oblique unstable modes of much larger wave length than the collisionless ones. Transition from the collisionless regime to the resistive one is thus documented and found discontinuous., Comment: To appear in PPCF

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., Comment: To appear in Physics of Plasmas