skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The synchrotron maser emission from relativistic shocks in Fast Radio Bursts: 1D PIC simulations of cold pair plasmas
ABSTRACT The emission process of Fast Radio Bursts (FRBs) remains unknown. We investigate whether the synchrotron maser emission from relativistic shocks in a magnetar wind can explain the observed FRB properties. We perform particle-in-cell (PIC) simulations of perpendicular shocks in cold pair plasmas, checking our results for consistency among three PIC codes. We confirm that a linearly polarized X-mode wave is self-consistently generated by the shock and propagates back upstream as a precursor wave. We find that at magnetizations σ ≳ 1 (i.e. ratio of Poynting flux to particle energy flux of the pre-shock flow) the shock converts a fraction $$f_\xi ^{\prime } \approx 7 \times 10^{-4}/\sigma ^2$$ of the total incoming energy into the precursor wave, as measured in the shock frame. The wave spectrum is narrow-band (fractional width ≲1−3), with apparent but not dominant line-like features as many resonances concurrently contribute. The peak frequency in the pre-shock (observer) frame is $$\omega ^{\prime \prime }_{\rm peak} \approx 3 \gamma _{\rm s | u} \omega _{\rm p}$$, where γs|u is the shock Lorentz factor in the upstream frame and ωp the plasma frequency. At σ ≳ 1, where our estimated $$\omega ^{\prime \prime }_{\rm peak}$$ differs from previous works, the shock structure presents two solitons separated by a cavity, and the peak frequency corresponds to an eigenmode of the cavity. Our results provide physically grounded inputs for FRB emission models within the magnetar scenario.  more » « less
Award ID(s):
1804048
PAR ID:
10127060
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
485
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
p. 3816-3833
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT Electromagnetic precursor waves generated by the synchrotron maser instability at relativistic magnetized shocks have been recently invoked to explain the coherent radio emission of fast radio bursts. By means of 2D particle-in-cell simulations, we explore the properties of the precursor waves in relativistic electron–positron perpendicular shocks as a function of the pre-shock magnetization σ ≳ 1 (i.e. the ratio of incoming Poynting flux to particle energy flux) and thermal spread Δγ ≡ kT/mc2 = 10−5−10−1. We measure the fraction fξ of total incoming energy that is converted into precursor waves, as computed in the post-shock frame. At fixed magnetization, we find that fξ is nearly independent of temperature as long as Δγ ≲ 10−1.5 (with only a modest decrease of a factor of 3 from Δγ = 10−5 to Δγ = 10−1.5), but it drops by nearly two orders of magnitude for Δγ ≳ 10−1. At fixed temperature, the scaling with magnetization $$f_\xi \sim 10^{-3}\, \sigma ^{-1}$$ is consistent with our earlier 1D results. For our reference σ = 1, the power spectrum of precursor waves is relatively broad (fractional width ∼1 − 3) for cold temperatures, whereas it shows pronounced line-like features with fractional width ∼0.2 for 10−3 ≲ Δγ ≲ 10−1.5. For σ ≳ 1, the precursor waves are beamed within an angle ≃σ−1/2 from the shock normal (as measured in the post-shock frame), as required so they can outrun the shock. Our results can provide physically grounded inputs for FRB emission models based on maser emission from relativistic shocks. 
    more » « less
  2. Abstract Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (vsh≈ (5500–22000) km s−1) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (vsh≈ 300 km s−1), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations. 
    more » « less
  3. Abstract Particle-in-cell (PIC) simulations have shown that relativistic collisionless shocks mediated by the Weibel instability accelerate ∼1% of incoming particles, while the majority are transmitted through the shock and become thermalized. The microphysical processes that determine whether an incoming particle will be transmitted or reflected are poorly understood. We study the microphysics of particle reflection in Weibel-mediated shocks by tracking a shell of test particles in a PIC simulation of a shock in pair plasma. We find that electrons in positron-dominated filaments and positrons in electron-dominated filaments efficiently reflect off of strong magnetic structures at the shock. To participate in diffusive shock acceleration, however, these reflected particles headed toward the upstream must avoid getting advected downstream. This is enabled by incoming filaments, which trap reflected particles carrying the same sign of current as the filaments. The final injection efficiency on the order of ∼1% thus results from the effectiveness of the initial reflection at the shock and the reflected particles’ probability of survival in the upstream postreflection. We develop a model that predicts the fraction of high-energy particles as a function of the properties of Weibel filamentation. 
    more » « less
  4. Abstract Collisionless shocks tend to send charged particles into the upstream, driving electric currents through the plasma. Using kinetic particle-in-cell simulations, we investigate how the background thermal plasma neutralizes such currents in the upstream of quasi-parallel non-relativistic electron–proton shocks. We observe distinct processes in different regions: the far upstream, the shock precursor, and the shock foot. In the far upstream, the current is carried by nonthermal protons, which drive electrostatic modes and produce suprathermal electrons that move toward upstream infinity. Closer to the shock (in the precursor), both the current density and the momentum flux of the beam increase, which leads to electromagnetic streaming instabilities that contribute to the thermalization of suprathermal electrons. At the shock foot, these electrons are exposed to shock-reflected protons, resulting in a two-stream type instability. We analyze these processes and the resulting heating through particle tracking and controlled simulations. In particular, we show that the instability at the shock foot can make the effective thermal speed of electrons comparable to the drift speed of the reflected protons. These findings are important for understanding both the magnetic field amplification and the processes that may lead to the injection of suprathermal electrons into diffusive shock acceleration. 
    more » « less
  5. Abstract The observed energy spectra of accelerated particles at interplanetary shocks often do not match the diffusive shock acceleration (DSA) theory predictions. In some cases, the particle flux forms a plateau over a wide range of energies, extendingupstream of the shockfor up to seven fluxe-folds before submerging into the background spectrum. Remarkably, at and downstream of the shock we have studied in detail, the flux falls off in energy asϵ−1, consistent with the DSA prediction for a strong shock. The upstream plateau suggests a particle transport mechanism different from those traditionally employed in DSA models. We show that a standard (linear) DSA solution based on a widely accepted diffusive particle transport with an underlying resonant wave–particle interaction is inconsistent with the plateau in the particle flux. To resolve this contradiction, we modify the DSA theory in two ways. First, we include a dependence of the particle diffusivityκon the particle fluxF(nonlinear particle transport). Second, we invoke short-scale magnetic perturbations that are self-consistently generated by, but not resonant with, accelerated particles. They lead to the particle diffusivity increasing with the particle energy as ∝ϵ3/2that simultaneously decreases with the particle flux as 1/F. The combination of these two trends results in the flat spectrum upstream. We speculate that nonmonotonic spatial variations of the upstream spectrum, apart from being time-dependent, may also result from non-DSA acceleration mechanisms at work upstream, such as stochastic Fermi or magnetic pumping acceleration. 
    more » « less