skip to main content


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
NSF-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

    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
  3. 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−4and pair-loading factorsZ±≲ 10 are studied, whereZ±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σLthat decreases withZ±. AtσσLthe shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease withZ±, 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 asE¯eZ±+11. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ≈ 5 × 10−6. The ions then become essentially thermal with mean energyE¯i, while electrons form a nonthermal tail, extending fromEZ±+11E¯itoE¯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
  4. ABSTRACT During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS’s open magnetic flux – and hence increases its wind luminosity – relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral wind may (depending on pair loading) secularly increase, leading to self-interaction and internal shocks in the outflow beyond the binary orbit. The magnetized forward shock can generate coherent radio emission via the synchrotron maser process, resulting in an observable radio precursor to binary NS merger. We perform 1D relativistic hydrodynamical simulations of shock interaction in the accelerating binary NS wind, assuming that the inspiral wind efficiently converts its Poynting flux into bulk kinetic energy prior to the shock radius. This is combined with the shock maser spectrum from particle-in-cell simulations, to generate synthetic radio light curves. The precursor burst with a fluence of ∼1 Jy·ms at ∼GHz frequencies lasts ∼1–500 ms following the merger for a source at ∼3 Gpc (Bd/1012 G)8/9, where Bd is the dipole field strength of the more strongly magnetized star. Given an outflow geometry concentrated along the binary equatorial plane, the signal may be preferentially observable for high-inclination systems, that is, those least likely to produce a detectable gamma-ray burst. 
    more » « less
  5. ABSTRACT

    At least some fast radio bursts (FRBs) are produced by magnetars. Even though mounting observational evidence points towards a magnetospheric origin of FRB emission, the question of the location for FRB generation continues to be debated. One argument suggested against the magnetospheric origin of bright FRBs is that the radio waves associated with an FRB may lose most of their energy before escaping the magnetosphere because the cross-section for e± to scatter large-amplitude electromagnetic waves in the presence of a strong magnetic field is much larger than the Thompson cross-section. We have investigated this suggestion and find that FRB radiation travelling through the open field line region of a magnetar’s magnetosphere does not suffer much loss due to two previously ignored factors. First, the plasma in the outer magnetosphere ($r \gtrsim 10^9$ cm), where the losses are potentially most severe, is likely to be flowing outwards at a high Lorentz factor γp ≥ 103. Secondly, the angle between the wave vector and the magnetic field vector, θB, in the outer magnetosphere is likely of the order of 0.1 radian or smaller due in part to the intense FRB pulse that tilts open magnetic field lines so that they get aligned with the pulse propagation direction. Both these effects reduce the interaction between the FRB pulse and the plasma substantially. We find that a bright FRB with an isotropic luminosity $L_{\rm frb} \gtrsim 10^{42} \, {\rm erg \ s^{-1}}$ can escape the magnetosphere unscathed for a large section of the γp − θB parameter space, and therefore conclude that the generation of FRBs in magnetar magnetosphere passes this test.

     
    more » « less