skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

more » « less
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 3034-3045
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    The radiation drag in photon-rich environments of cosmic explosions can seed kinetic instabilities by inducing velocity spreads between relativistically streaming plasma components. Such microturbulence is likely imprinted on the breakout signals of radiation-mediated shocks. However, large-scale, transverse magnetic fields in the deceleration region of the shock transition can suppress the dominant kinetic instabilities by preventing the development of velocity separations between electron–positron pairs and a heavy ion species. We use a 1D five-fluid radiative transfer code to generate self-consistent profiles of the radiation drag force and plasma composition in the deceleration region. For increasing magnetization, our models predict rapidly growing pair multiplicities and a substantial radiative drag developing self-similarly throughout the deceleration region. We extract the critical magnetization parameter σc, determining the limiting magnetic field strength at which a three-species plasma can develop kinetic instabilities before reaching the isotropized downstream. For a relativistic, single ion plasma drifting with γu = 10 in the upstream of a relativistic radiation-mediated shock, we find the threshold σc ≈ 10−7 for the onset of microturbulence. Suppression of plasma instabilities in the case of multi-ion composition would likely require much higher values of σc. Identifying high-energy signatures of microturbulence in shock breakout signals and combining them with the magnetization limits provided in this work will allow a deeper understanding of the magnetic environment of cosmic explosions like supernovae, gamma-ray bursts, and neutron star binary mergers.

    more » « less
  2. Abstract

    Over the last few decades, different types of plasma waves (e.g., the ion acoustic waves (IAWs), electrostatic solitary waves, upper/lower hybrid waves, and Langmuir waves) have been observed in the upstream, downstream, and ramp regions of the collisionless interplanetary (IP) shocks. These waves may appear as short-duration (only a few milliseconds at 1 au) electric field signatures in the in-situ measurements, with typical frequencies of ∼1–10 kHz. A number of IAW features at the IP shocks seem to be unexplained by kinetic models and require a new modeling effort. Thus, this paper is dedicated to bridging this gap in understanding. In this paper, we model the linear IAWs inside the shock ramp by devising a novel linearization method for the two-fluid magnetohydrodynamic equations with spatially dependent shock parameters. It is found that, for parallel propagating waves, the linear dispersion relation leads to a finite growth rate, which is dependent on the shock density compression ratio, as Wind data suggest. Further analysis reveals that the wave frequency grows towards the downstream region within the shock ramp, and the wave growth rate is independent of the electron-to-ion temperature ratio, as Magnetospheric Multiscale (MMS) in-situ measurements suggest, and is uniform within the shock ramp. Thus, this study helps in understanding the characteristics of the IAWs at the collisionless IP shocks.

    more » « less
  3. Abstract

    We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (βp= 0.25) collisionless ion–electron shocks with mass ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ, that scales roughly linearly with the electron temperature jump. Our simulations haveΔϕ/(0.5miush2)0.1–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕin our low-βpshocks. We further focus on twoθBn= 65° shocks: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis the ion skin depth. Within the precursors,ϕhas a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

    more » « less
  4. Abstract

    The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.

    more » « less
  5. This paper develops a unified linear theory of cross field plasma instabilities, including the Farley–Buneman, electron thermal, and ion thermal instabilities, in spatially uniform collisional plasmas with partially unmagnetized multi-species ions. Collisional plasma instabilities in weakly ionized, highly dissipative, weakly magnetized plasmas play an important role in the lower Earth's ionosphere and may be of importance in other planetary ionospheres, stellar atmospheres, cometary tails, molecular clouds, accretion disks, etc. In the Earth's ionosphere, these collisional plasma instabilities cause intense electron heating. In the solar chromosphere, they can do the same—an effect originally suggested from spectroscopic observations and modeling. Based on a simplified 5-moment multi-fluid model, the theoretical analysis presented in this paper produces the linear dispersion relation for the combined Thermal Farley–Buneman Instability with an important long-wavelength limit analyzed in detail. This limit provides an easy interpretation of different instability drivers and wave dissipation. This analysis of instability, combined with simulations, will enable us to better understand plasma waves and turbulence in these commonly occurring collisional space plasmas.

    more » « less