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.
- Award ID(s):
- 1716567
- NSF-PAR ID:
- 10291263
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 499
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2884 to 2895
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
ABSTRACT -
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
-
Context. Standing and moving shocks in relativistic astrophysical jets are very promising sites for particle acceleration to large Lorentz factors and for the emission from the radio up to the γ -ray band. They are thought to be responsible for at least part of the observed variability in radio-loud active galactic nuclei. Aims. We aim to simulate the interactions of moving shock waves with standing recollimation shocks in structured and magnetized relativistic jets and to characterize the profiles of connected flares in the radio light curve. Methods. Using the relativistic magneto-hydrodynamic code MPI-AMRVAC and a radiative transfer code in post-processing, we explore the influence of the magnetic-field configuration and transverse stratification of an over-pressured jet on its morphology, on the moving shock dynamics, and on the emitted radio light curve. First, we investigate different large-scale magnetic fields with their effects on the standing shocks and on the stratified jet morphology. Secondly, we study the interaction of a moving shock wave with the standing shocks. We calculated the synthetic synchrotron maps and radio light curves and analyze the variability at two frequencies 1 and 15.3 GHz and for several observation angles. Finally, we compare the characteristics of our simulated light curves with radio flares observed from the blazar 3C 273 with the Owens Valley Radio Observatory and Very Long Baseline Array in the MOJAVE survey between 2008 and 2019. Results. We find that in a structured over-pressured relativistic jet, the presence of the large-scale magnetic field structure changes the properties of the standing shock waves and leads to an opening in the jet. The interaction between waves from inner and outer jet components can produce strong standing shocks. When crossing such standing shocks, moving shock waves accompanying overdensities injected in the base of the jet cause very luminous radio flares. The observation of the temporal structure of these flares under different viewing angles probes the jet at different optical depths. At 1 GHz and for small angles, the self-absorption caused by the moving shock wave becomes more important and leads to a drop in the observed flux after it interacts with the brightest standing knot. A weak asymmetry is seen in the shape of the simulated flares, resulting from the remnant emission of the shocked standing shocks. The characteristics of the simulated flares and the correlation of peaks in the light curve with the crossing of moving and standing shocks favor this scenario as an explanation of the observed radio flares of 3C 273.more » « less
-
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 as . (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ ≈ 5 × 10−6. The ions then become essentially thermal with mean energy , while electrons form a nonthermal tail, extending from to . 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. -
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 ratiom i/m e= 200, fast Mach number –4, and upstream magnetic field angleθ Bn= 55°–85° from the shock normal . It is known that shock electron heating is described by an ambipolar, -parallel electric potential jump, ΔB ϕ ∥, that scales roughly linearly with the electron temperature jump. Our simulations have –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: a ( ) case with a long, 30d iprecursor of whistler waves along , and a ( ) case with a shorter, 5d iprecursor of whistlers oblique to both and ;B d iis 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 the ,θ Bn= 65° case,ϕ ∥shows a weak dependence on the electron plasma-to-cyclotron frequency ratioω pe/Ωce, andϕ ∥decreases by a factor of 2 asm i/m eis raised to the true proton–electron value of 1836.