We examine with particle-in-cell (PIC) simulations how a parallel shock in pair plasma reacts to upstream waves, which are driven by escaping downstream particles. Initially, the shock is sustained in the two-dimensional simulation by a magnetic filamentation (beam-Weibel) instability. Escaping particles drive an electrostatic beam instability upstream. Modifications of the upstream plasma by these waves hardly affect the shock. In time, a decreasing density and increasing temperature of the escaping particles quench the beam instability. A larger thermal energy along than perpendicular to the magnetic field destabilizes the pair-Alfv\'en mode. In the rest frame of the upstream plasma, the group velocity of the growing pair-Alfv\'en waves is below that of the shock and the latter catches up with the waves. Accumulating pair-Alfv\'en waves gradually change the shock in the two-dimensional simulation from a Weibel-type shock into an Alfv\'enic shock with a Mach number that is about 6 for our initial conditions.

We study with a one-dimensional particle-in-cell (PIC) simulation the expansion of a pair cloud into a magnetized electron-proton plasma as well as the formation and subsequent propagation of a tangential discontinuity that separates both plasmas. Its propagation speed takes the value that balances the magnetic pressure of the discontinuity against the thermal pressure of the pair cloud and the ram pressure of the protons. Protons are accelerated by the discontinuity to a speed that exceeds the fast magnetosonic speed by the factor 10. A supercritical fast magnetosonic shock forms at the front of this beam. An increasing proton temperature downstream of the shock and ahead of the discontinuity leaves the latter intact. We create the discontinuity by injecting a pair cloud at a simulation boundary into a uniform electron-proton plasma, which is permeated by a perpendicular magnetic field. Collisionless tangential discontinuities in the relativistic pair jets of X-ray binaries (microquasars) are in permanent contact with the relativistic leptons of its inner cocoon and they become sources of radio synchrotron emissions.

We present the results from a particle-in-cell (PIC) simulation that models the interaction between a spatially localized electron-positron cloud and an electron-ion plasma. The latter is permeated by a magnetic field that is initially spatially uniform and aligned with the mean velocity vector of the pair cloud. The pair cloud expels the magnetic field and piles it up into an electromagnetic piston. Its electromagnetic field is strong enough to separate the pair cloud from the ambient plasma in the direction that is perpendicular to the cloud propagation direction. The piston propagates away from the spine of the injected pair cloud and it accelerates the protons to a high nonrelativistic speed. The accelerated protons form an outer cocoon that will eventually become separated from the unperturbed ambient plasma by a fast magnetosonic shock. No electromagnetic piston forms at the front of the cloud and a shock is mediated here by the filamentation instability. The final plasma distribution resembles that of a hydrodynamic jet. Collisionless plasma jets may form in the coronal plasma of accreting black holes and the interaction between the strong magnetic field of the piston and the hot pair cloud may contribute to radio emissions by such objects.

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.

We study the effect a guiding magnetic field has on the formation and structure of a pair jet that propagates through a collisionless electron-proton plasma at rest. We model with a PIC simulation a pair cloud with the temperature 400 keV and mean speed 0.9c. The cloud propagates through a spatially uniform, magnetized and cool ambient electron-proton plasma that is at rest. Its mean velocity vector is aligned with the background magnetic field. A jet forms in time. Its outer cocoon consists of jet-accelerated ambient plasma and is separated from the inner cocoon by an electromagnetic piston with a thickness that is comparable to the thermal gyroradius of jet particles. A supercritical fast magnetosonic shock will form between the pristine ambient plasma and the jet-accelerated one on a time scale that exceeds our simulation time by an order of magnitude. The inner cocoon is pair plasma that lost its directed flow energy while it swept out the background magnetic field. A beam of electrons and positrons moves along the jet spine at its initial speed. Its electrons are slowed down and some positrons are accelerated as they cross the jet's head. The latter escape upstream along the magnetic field, which yields an excess of MeV positrons ahead of the jet. Some of the protons, which were located behind the electromagnetic piston at the time it formed, are accelerated to MeV energies

By modelling the expansion of a cloud of electrons and positrons with the temperature 400 keV that propagates at the mean speed 0.9c ($c:$ speed of light) through an initially unmagnetized electron-proton plasma with a particle-in-cell (PIC) simulation, we find a mechanism that collimates the pair cloud into a jet. A filamentation instability develops between the protons at rest and the moving positrons. Its magnetic field collimates the positrons and drives an electrostatic shock into the electron-proton plasma. The magnetic field acts as a discontinuity that separates the protons of the shocked ambient plasma, known as the outer cocoon, from the jet's interior region. The outer cocoon expands at the speed 0.15c along the jet axis and at 0.03c perpendicularly to it. The filamentation instability converts the jet's directed flow energy into magnetic energy in the inner cocoon. The magnetic discontinuity cannot separate the ambient electrons from the jet electrons. Both species rapidly mix and become indistinguishable. The spatial distribution of the positive charge carriers is in agreement with the distributions of the ambient material and the jet material predicted by a hydrodynamic model apart from a dilute positronic outflow that is accelerated by the electromagnetic field at the jet's head.

The expansion of a charge-neutral cloud of electrons and positrons with the temperature 1 MeV into an unmagnetized ambient plasma is examined with a 2D particle-in-cell (PIC) simulation. The pair outflow drives solitary waves in the ambient protons. Their bipolar electric fields attract electrons of the outflowing pair cloud and repel positrons. These fields can reflect some of the protons thereby accelerating them to almost an MeV. Ion acoustic solitary waves are thus an efficient means to couple energy from the pair cloud to protons. The scattering of the electrons and positrons by the electric field slows down their expansion to a nonrelativistic speed. Only a dilute pair outflow reaches the expansion speed expected from the cloud's thermal speed. Its positrons are more energetic than its electrons. In time an instability grows at the front of the dense slow-moving part of the pair cloud, which magnetizes the plasma. The instability is driven by the interaction of the outflowing positrons with the protons. These results shed light on how magnetic fields are created and ions are accelerated in pair-loaded astrophysical jets and winds.

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.

We report on the first experimental observation of a current-driven instability developing in a quasi-neutral matter-antimatter beam. Strong magnetic fields ($\geq$ 1 T) are measured, via means of a proton radiography technique, after the propagation of a neutral electron-positron beam through a background electron-ion plasma.The experimentally determined equipartition parameter of $\epsilon_B \approx 10^{-3}$, is typical of values inferred from models of astrophysical gamma-ray bursts, in which the relativistic flows are also expected to be pair dominated. The data, supported by Particle-In-Cell simulations and simple analytical estimates, indicate that these magnetic fields persist in the background plasma for thousands of inverse plasma frequencies. The existence of such long-lived magnetic fields can be related to analog astrophysical systems, such as those prevalent in lepton-dominated jets.

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.

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.

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.

A pair of curved shocks in a collisionless plasma is examined with a two-dimensional particle-in-cell (PIC) simulation. The shocks are created by the collision of two electron-ion clouds at a speed that exceeds everywhere the threshold speed for shock formation. A variation of the collision speed along the initially planar collision boundary, which is comparable to the ion acoustic speed, yields a curvature of the shock that increases with time. The spatially varying Mach number of the shocks results in a variation of the downstream density in the direction along the shock boundary. This variation is eventually equilibrated by the thermal diffusion of ions. The pair of shocks is stable for tens of inverse ion plasma frequencies. The angle between the mean flow velocity vector of the inflowing upstream plasma and the shock's electrostatic field increases steadily during this time. The disalignment of both vectors gives rise to a rotational electron flow, which yields the growth of magnetic field patches that are coherent over tens of electron skin depths.

The jets of compact accreting objects are composed of electrons and a mixture of positrons and ions. These outflows impinge on the interstellar or intergalactic medium and both plasmas interact via collisionless processes. Filamentation (beam-Weibel) instabilities give rise to the growth of strong electromagnetic fields. These fields thermalize the interpenetrating plasmas. Hitherto, the effects imposed by a spatial non-uniformity on filamentation instabilities have remained unexplored. We examine the interaction between spatially uniform background electrons and a minuscule cloud of electrons and positrons. A square micro-cloud of equally dense electrons and positrons impinges in our particle-in-cell (PIC) simulation on a spatially uniform plasma at rest. The mean speed of the micro-cloud corresponds to a relativistic factor of 15, which is relevant for laboratory experiments and for relativistic astrophysical outflows. The spatial distributions of the leptons and of the electromagnetic fields are examined at several times. A filamentation instability develops between the magnetic field carried by the micro-cloud and the background electrons. The electromagnetic fields, which grow from noise levels, redistribute the electrons and positrons within the cloud, which boosts the peak magnetic field amplitude. The current density and the moduli of the electromagnetic fields grow aperiodically in time and steadily along the direction that is anti-parallel to the cloud's velocity vector. The micro-cloud remains conjoined during the simulation. The instability induces an electrostatic wakefield in the background plasma.

The collision of two plasma clouds at a speed that exceeds the ion acoustic speed can result in the formation of shocks. This phenomenon is observed not only in astrophysical scenarios such as the propagation of supernova remnant (SNR) blast shells into the interstellar medium, but also in laboratory-based laser-plasma experiments. These experiments and supporting simulations are thus seen as an attractive platform for the small-scale reproduction and study of astrophysical shocks in the laboratory. We model two plasma clouds, which consist of electrons and ions, with a 2D PIC simulation. The ion temperatures of both clouds differ by a factor of 10. Both clouds collide at a speed, which is realistic for laboratory studies and for SNR shocks in their late evolution phase like that of RCW86. A magnetic field, which is orthogonal to the simulation plane, has a strength that is comparable to that at SNR shocks. A forward shock forms between the overlap layer of both plasma clouds and the cloud with the cooler ions. A large-amplitude ion acoustic wave is observed between the overlap layer and the cloud with the hotter ions. It does not steepen into a reverse shock, because its speed is below the ion acoustic speed. A gradient of the magnetic field amplitude builds up close to the forward shock as it compresses the magnetic field. This gradient gives rise to an electron drift that is fast enough to trigger an instability. Electrostatic ion acoustic wave turbulence develops ahead of the shock. It widens its transition layer and thermalizes the ions, but the forward shock remains intact.

We report on the temporally and spatially resolved detection of the precursory stages that lead to the formation of an unmagnetized, supercritical collision-less shock in a laser-driven laboratory experiment. The measured evolution of the electrostatic potential associated with the shock unveils the transition from a current free double layer into a symmetric shock structure, stabilized by ion reflection at the shock front. Supported by a matching Particle-In-Cell simulation and theoretical considerations, we suggest that this process is analogeous to ion reflection at supercritical collisionless shocks in supernova remnants.

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.

Gamma ray bursts (GRBs) are thought to originate from highly relativistic jets. The fireball model predicts internal shocks in the jets, causing magnetic field to be amplified and particles to be accelerated. We model the effects of an asymmetric density configuration for an internal plasma collision in a quasi-parallel magnetic field. We measured electron acceleration and found that a tenuous population of electrons is accelerated to Lorentz factors of \sim 300 - close to energy equipartition with ions. We found that the filaments did not remain static, but were deflected by the Lorentz force and rolled up into small vortices, which themselves merge to form a larger current vortex.

The energetic electromagnetic eruptions observed during the prompt phase of gamma-ray bursts are attributed to synchrotron emissions. The internal shocks moving through the ultrarelativistic jet, which is ejected by an imploding supermassive star, are the likely source of this radiation. Synchrotron emissions at the observed strength require the simultaneous presence of powerful magnetic fields and highly relativistic electrons. We explore with one and three-dimensional relativistic particle-in-cell simulations the transition layer of a shock, that evolves out of the collision of two plasma clouds at a speed 0.9c and in the presence of a quasi-parallel magnetic field. The cloud densities vary by a factor of 10. The number densities of ions and electrons in each cloud, which have the mass ratio 250, are equal. The peak Lorentz factor of the electrons is determined in the 1D simulation, as well as the orientation and the strength of the magnetic field at the boundary of the two colliding clouds. The relativistic masses of the electrons and ions close to the shock transition layer are comparable as in previous work. The 3D simulation shows rapid and strong plasma filamentation behind the transient precursor. The magnetic field component orthogonal to the initial field direction is amplified in both simulations to values that exceed those expected from the shock compression by over an order of magnitude. The forming shock is quasi-perpendicular due to this amplification. The simultaneous presence of highly relativistic electrons and strong magnetic fields will give rise to significant synchrotron emissions.

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

An instability driven by the thermal anisotropy of a single electron species is investigated in a 2D particle-in-cell (PIC) simulation. This instability is the one considered by Weibel and it differs from the beam driven filamentation instability. A comparison of the simulation results with analytic theory provides similar exponential growth rates of the magnetic field during the linear growth phase of the instability. We observe in accordance with previous works the growth of electric fields during the saturation phase of the instability. Some components of this electric field are not accounted for by the linearized theory. A single-fluid-based theory is used to determine the source of this nonlinear electric field. It is demonstrated that the magnetic stress tensor, which vanishes in a 1D geometry, is more important in this 2-dimensional model used here. The electric field grows to an amplitude, which yields a force on the electrons that is comparable to the magnetic one. The peak energy density of each magnetic field component in the simulation plane agrees with previous estimates. Eddy currents develop, which let the amplitude of the third magnetic field component grow, which is not observed in a 1D simulation.

A 2D particle simulation models the collision of two electron-ion plasma clouds along a quasi-parallel magnetic field. The collision speed is 0.9c and the density ratio 10. A current sheet forms at the front of the dense cloud, in which the electrons and the magnetic field reach energy equipartition with the ions. A structure composed of a solenoidal and a toroidal magnetic field grows in this sheet. It resembles that in the cross-section of the torus of a force-free spheromak, which may provide the coherent magnetic fields in gamma-ray burst (GRB) jets needed for their prompt emissions.

Plasma processes close to SNR shocks result in the amplification of magnetic fields and in the acceleration of electrons, injecting them into the diffusive acceleration mechanism. The acceleration of electrons and the B field amplification by the collision of two plasma clouds, each consisting of electrons and ions, at a speed of 0.5c is investigated. A quasi-parallel guiding magnetic field, a cloud density ratio of 10 and a plasma temperature of 25 keV are considered. A quasi-planar shock forms at the front of the dense plasma cloud. It is mediated by a circularly left-hand polarized electromagnetic wave with an electric field component along the guiding magnetic field. Its propagation direction is close to that of the guiding field and orthogonal to the collision boundary. It has a low frequency and a wavelength that equals several times the ion inertial length, which would be indicative of a dispersive Alfven wave close to the ion cyclotron resonance frequency of the left-handed mode (ion whistler), provided that the frequency is appropriate. However, it moves with the super-alfvenic plasma collision speed, suggesting that it is an Alfven precursor or a nonlinear MHD wave such as a Short Large-Amplitude Magnetic Structure (SLAMS). The growth of the magnetic amplitude of this wave to values well in excess of those of the quasi-parallel guiding field and of the filamentation modes results in a quasi-perpendicular shock. We present evidence for the instability of this mode to a four wave interaction. The waves developing upstream of the dense cloud give rise to electron acceleration ahead of the collision boundary. Energy equipartition between the ions and the electrons is established at the shock and the electrons are accelerated to relativistic speeds.

The Weibel instability is investigated with PIC simulations of an initially unmagnetized and spatially uniform electron plasma. This instability, which is driven by the thermally anisotropic electron distribution, generates electromagnetic waves with wave vectors perpendicular to the direction of the higher temperature. Two simulations are performed: A 2D simulation, with a simulation plane that includes the direction of higher temperature, demonstrates that the wave spectrum is initially confined to one dimension. The electric field components in the simulation plane generated by the instability equalize at the end of the simulation through a secondary instability. A 1D PIC simulation with a high resolution, where the simulation box is aligned with the wave vectors of the growing waves, reveals details of the electron phase space distribution and permits a comparison of the magnetic and electric fields when the instability saturates. It is shown that the electrostatic field is driven by the magnetic pressure gradient and that it and the magnetic field redistribute the electrons in space.

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.

In continuation of previous work, numerical results are presented, concerning relativistically counter-streaming plasmas. Here, the relativistic mixed mode instability evolves through, and beyond, the linear saturation -- well into the nonlinear regime. Besides confirming earlier findings, that wave power initially peaks on the mixed mode branch, it is observed that, during late time evolution wave power is transferred to other wave numbers. It is argued that the isotropization of power in wavenumber space may be a consequence of weak turbulence. Further, some modifications to the ideal weak turbulence limit is observed. Development of almost isotropic predominantly electrostatic -- partially electromagnetic -- turbulent spectra holds relevance when considering the spectral emission signatures of the plasma, namely bremsstrahlung, respectively magneto-bremsstrahlung (synchrotron radiation and jitter radiation) from relativistic shocks in astrophysical jets and shocks from gamma-ray bursts and active galactic nuclei.

Supernova remnant (SNR) blast shells can reach the flow speed $v_s = 0.1 c$ and shocks form at its front. Instabilities driven by shock-reflected ion beams heat the plasma in the foreshock, which may inject particles into diffusive acceleration. The ion beams can have the speed $v_b \approx v_s$. For $v_b \ll v_s$ the Buneman or upper-hybrid instabilities dominate, while for $v_b \gg v_s$ the filamentation and mixed modes grow faster. Here the relevant waves for $v_b \approx v_s$ are examined and how they interact nonlinearly with the particles. The collision of two plasma clouds at the speed $v_s$ is modelled with particle-in-cell (PIC) simulations, which convect with them magnetic fields oriented perpendicular to their flow velocity vector. One simulation models equally dense clouds and the other one uses a density ratio of 2. Both simulations show upper-hybrid waves that are planar over large spatial intervals and that accelerate electrons to $\sim$ 10 keV. The symmetric collision yields only short oscillatory wave pulses, while the asymmetric collision also produces large-scale electric fields, probably through a magnetic pressure gradient. The large-scale fields destroy the electron phase space holes and they accelerate the ions, which facilitates the formation of a precursor shock.

We model relativistically colliding plasma by PIC simulations in one and two spatial dimensions, taking an ion-to-electron mass ratio of 400. Energy dissipation by a wave precursor of mixed polarity and different densities of the colliding plasma slabs results in a relativistic forward shock forming on millisecond timescales. The forward shock accelerates electrons to ultrarelativistic energies and reflects upstream ions, which drag the electrons along to preserve the plasma quasi-neutrality. No reverse shock forms. The shock may be representative for internal gamma ray burst shocks.

Particle acceleration by means of non-linear plasma wave interactions is of great topical interest. Accordingly, in this paper we focus on the electron surfing process. Self-consistent kinetic simulations, using both relativistic Vlasov and PIC (Particle In Cell) approaches, show here that electrons can be accelerated to highly relativistic energies (up to 100 m_e c^2) if the phase speed of the electrostatic wave is mildly relativistic (0.6c to 0.9c for the magnetic field strengths considered). The acceleration is strong because of relativistic stabilisation of the nonlinearly saturated electrostatic wave, seen in both relativistic Vlasov and PIC simulations. An inverse power law momentum distribution can arise for the most strongly accelerated electrons. These results are of relevance to observed rapid changes in the radio synchrotron emission intensities from microquasars, gamma ray bursts and other astrophysical objects that require rapid acceleration mechanisms for electrons.

Observations of synchrotron radiation across a wide range of wavelengths provide clear evidence that electrons are accelerated to relativistic energies in supernova remnants (SNRs). However, a viable mechanism for the pre-acceleration of such electrons to mildly relativistic energies has not yet been established. In this paper an electromagnetic particle-in-cell (PIC) code is used to simulate acceleration of electrons from background energies to tens of keV at perpendicular collisionless shocks associated with SNRs. Free energy for electron energization is provided by ions reflected from the shock front, with speeds greater than the upstream electron thermal speed. The PIC simulation results contain several new features, including: the acceleration, rather than heating, of electrons via the Buneman instability; the acceleration of electrons to speeds exceeding those of the shock-reflected ions producing the instability; and strong acceleration of electrons perpendicular to the magnetic field. Electron energization takes place through a variety of resonant and non-resonant processes, of which the strongest involves stochastic wave-particle interactions. In SNRs the diffusive shock process could then supply the final step required for the production of fully relativistic electrons. The mechanisms identified in this paper thus provide a possible solution to the electron pre-acceleration problem.