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:
10128596
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
485
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3816 to 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 We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations σ ≲ 10 −4 and pair-loading factors Z ± ≲ 10 are studied, where Z ± is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field when σ exceeds a critical value σ L that decreases with Z ± . At σ ≲ σ L the shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease with Z ± , leading to lower σ L . (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as E ¯ e ∝ Z ± + 1 − 1 . (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as σ ≈ 5 × 10 −6 . The ions then become essentially thermal with mean energy E ¯ i , while electrons form a nonthermal tail, extending from E ∼ Z ± + 1 − 1 E ¯ i to E ¯ i . When σ = 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows. 
    more » « less
  5. Zank, Gary (Ed.)
    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