No Arabic abstract
Plasma outflows from gamma-ray bursts, supernovae, and relativistic jets, in general, interact with the surrounding medium through collisionless shocks. The microphysics of such shocks are still poorly understood, which, potentially, can introduce uncertainties in the interpretation of observations. It is now well established that the Weibel two-stream instability is capable of generating strong electromagnetic fields in the transition region between the jet and the ambient plasma. However, the parameter space of collisionless shocks is vast and still remains unexplored. In this Letter, we focus on how an ambient magnetic field affects the evolution of the electron Weibel instability and the associated shock. Using a particle-in-cell code, we have performed three-dimensional numerical experiments on such shocks. We compare simulations in which a jet is injected into an unmagnetized plasma with simulations in which the jet is injected into a plasma with an ambient magnetic field both parallel and perpendicular to the jet flow. We find that there exists a threshold of the magnetic field strength below which the Weibel two-stream instability dominates, and we note that the interstellar medium magnetic field strength lies well below this value. In the case of a strong magnetic field parallel to the jet, the Weibel instability is quenched. In the strong perpendicular case, ambient and jet electrons are strongly accelerated because of the charge separation between deflected jet electrons and less deflected jet ions. Also, the electromagnetic topologies become highly non-linear and complex with the appearance of anti-parallel field configurations.
In this second paper of a series, we discuss the dynamics of a plasma entering the precursor of an unmagnetized, relativistic collisionless pair shock. We discuss how this background plasma is decelerated and heated through its interaction with a microturbulence that results from the growth of a current filamentation instability (CFI) in the shock precursor. We make use, in particular, of the reference frame $mathcal R_{rm w}$ in which the turbulence is mostly magnetic. This frame moves at relativistic velocities towards the shock front at rest, decelerating gradually from the far to the near precursor. In a first part, we construct a fluid model to derive the deceleration law of the background plasma expected from the scattering of suprathermal particles off the microturbulence. This law leads to the relationship $gamma_{rm p},sim,xi_{rm b}^{-1/2}$ between the background plasma Lorentz factor $gamma_{rm p}$ and the normalized pressure of the beam $xi_{rm b}$; it is found to match nicely the spatial profiles observed in large-scale 2D3V particle-in-cell simulations. In a second part, we model the dynamics of the background plasma at the kinetic level, incorporating the inertial effects associated with the deceleration of $mathcal R_{rm w}$ into a Vlasov-Fokker-Planck equation for pitch-angle diffusion. We show how the effective gravity in $mathcal R_{rm w}$ drives the background plasma particles through friction on the microturbulence, leading to efficient plasma heating. Finally, we compare a Monte Carlo simulation of our model with dedicated PIC simulations and conclude that it can satisfactorily reproduce both the heating and the deceleration of the background plasma in the shock precursor, thereby providing a successful 1D description of the shock transition at the microscopic level.
Gamma ray bursts are among the most energetic events in the known universe. A highly relativistic fireball is ejected. In most cases the burst itself is followed by an afterglow, emitted under deceleration as the fireball plunges through the circum-stellar media. To interpret the observations of the afterglow emission, two physical aspects need to be understood: 1) The origin and nature of the magnetic field in the fireball and 2) the particle velocity distribution function behind the shock. Both are necessary in existing afterglow models to account for what is believed to be synchrotron radiation. To answer these questions, we need to understand the microphysics at play in collisionless shocks. Using 3D particle-in-cell simulations we can gain insight in the microphysical processes that take place in such shocks. We discuss the results of such computer experiments. It is shown how a Weibel-like two-stream plasma instability is able to create a strong transverse intermittent magnetic field and points to a connected mechanism for in situ particle acceleration in the shock region.
We study the long term evolution of magnetic fields generated by an initially unmagnetized collisionless relativistic $e^+e^-$ shock. Our 2D particle-in-cell numerical simulations show that downstream of such a Weibel-mediated shock, particle distributions are approximately isotropic, relativistic Maxwellians, and the magnetic turbulence is highly intermittent spatially, nonpropagating, and decaying. Using linear kinetic theory, we find a simple analytic form for these damping rates. Our theory predicts that overall magnetic energy decays like $(omega_p t)^{-q}$ with $q sim 1$, which compares favorably with simulations, but predicts overly rapid damping of short wavelength modes. Magnetic trapping of particles within the magnetic structures may be the origin of this discrepancy. We conclude that initially unmagnetized relativistic shocks in electron-positron plasmas are unable to form persistent downstream magnetic fields. These results put interesting constraints on synchrotron models for the prompt and afterglow emission from GRBs.
Relativistic astrophysical collisionless shocks represent outstanding dissipation agents of the huge power of relativistic outflows produced by accreting black holes, core collapsed supernovae and other objects into multi-messenger radiation (cosmic rays, neutrinos, electromagnetic radiation). This article provides a theoretical discussion of the fundamental physical ingredients of these extreme phenomena. In the context of weakly magnetized shocks, in particular, it is shown how the filamentation type instabilities, which develop in the precursor of pair dominated or electron-ion shocks, provide the seeds for the scattering of high energy particles as well as the agent which preheats and slows down the incoming precursor plasma. This analytical discussion is completed with a mesoscopic, non-linear model of particle acceleration in relativistic shocks based on Monte Carlo techniques. This Monte Carlo model uses a semi-phenomenological description of particle scattering which allows it to calculate the back-reaction of accelerated particles on the shock structure on length and momentum scales which are currently beyond the range of microscopic particle-in-cell (PIC) simulations.
During magnetically dominated relativistic reconnection, inflowing plasma depletes the initial relativistic pressure at the x-line and collisionless plasma heating inside the diffusion region is insufficient to overcome this loss. The resulting pressure drop causes a collapse at the x-line, essentially a localization mechanism of the diffusion region necessary for fast reconnection. The extension of this low-pressure region further explains the bursty nature of anti-parallel reconnection because a once opened outflow exhaust can also collapse, which repeatedly triggers secondary tearing islands. However, a stable single x-line reconnection can be achieved when an external guide field exists, since the reconnecting magnetic field component rotates out of the reconnection plane at outflows, providing additional magnetic pressure to sustain the opened exhausts.