Origin of Pulsed Radio Emission from Magnetars

Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron-positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag by the radiation field of the magnetar. This combination forces the plasma into a ``radiatively locked’’ state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power, and incorporate the local simulation results into a global magnetospheric model. The model explains key features of observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities $\sim 10^{30}{\rm{erg/s}}$, and a broad range of frequencies extending up to $\sim 100, \mathrm{GHz}$.

💡 Research Summary

The paper addresses the long‑standing puzzle of why a handful of magnetars emit bright, long‑lasting radio pulsations that differ markedly from ordinary radio pulsars. The authors propose and demonstrate, for the first time from first‑principles kinetic simulations, that the radio emission originates in a closed, twisted magnetic flux bundle anchored near the magnetic pole and extending to distances of order 100 km or more. Within this bundle, an electron‑positron (e±) plasma must carry the magnetospheric electric current j₀ required by the twisted field (j₀≈(c/4π)∇×B). Because the plasma is charge‑neutral (n₊≈n₋), the current is sustained by a small velocity difference δv between the two species. Simultaneously, the intense thermal X‑ray field from the hot neutron‑star surface exerts a strong radiative drag on the particles via resonant inverse‑Compton scattering. This drag drives both species toward a “radiatively‑locked” velocity v⋆ at which the drag vanishes. To maintain the required δv against the drag, a parallel electric field E₀ must be present.

In the two‑fluid description, the authors introduce the pair multiplicity M (the ratio of total charge flux to the required current). For M≫1 the required δv≈2c/M is tiny, and the radiative drag quickly pushes the particle momenta toward v⋆. However, the drag also compresses the distribution functions, and when the thermal spread of each species becomes comparable to the separation between the two streams, a two‑stream instability inevitably develops. The instability growth rate is set by the local plasma frequency ωₚ, while the cooling time t_cool (set by the drag) is orders of magnitude longer (ωₚ t_cool≈10³–10⁴ in the relevant region). Consequently, the instability can sustain a turbulent state for a long time, continuously feeding low‑frequency electromagnetic fluctuations.

To capture this subtle physics, the authors employ a relativistic Vlasov‑Maxwell solver based on the discontinuous Galerkin (DG) method (the Gkeyll framework). Standard particle‑in‑cell (PIC) techniques would be swamped by Poisson noise because the signal is weak; the DG approach resolves the distribution function directly in phase space. Simulations are performed in the radiatively‑locked (RL) frame, where the bulk flow is at v⋆ and the drag appears as a simple friction term F′_rad≈−p′/t′_cool.

1‑D (1X1V) electrostatic runs: The authors initialize a Maxwell‑Jüttner plasma with a prescribed multiplicity (M=100) and a small initial drift v_d. Uniform drag is applied with ωₚ t_cool=300, 1000, 3000. The simulations show rapid growth of electric field energy, saturating at ε_E≈2(ωₚ t_cool)⁻¹, in agreement with analytic expectations. The particle distribution evolves from a single thermal peak into two counter‑streaming beams separated by ±v_d, confirming that the drag alone cannot produce cold streams; instead, a sustained turbulent state persists.

2‑D (2X1V) electromagnetic runs: By allowing magnetic field perturbations (via Faraday’s law), the authors observe the emergence of genuine electromagnetic waves. Notably, superluminal modes (phase velocity > c) appear, which can escape the magnetosphere without being trapped by the plasma. The wave spectrum extends from hundreds of MHz up to tens of GHz, and the emitted power, when integrated over the simulated volume, matches the observed radio luminosities of magnetars (L_radio≈10³⁰ erg s⁻¹).

The authors then embed the local simulation results into a global magnetospheric model. They argue that after an X‑ray outburst, the twisted bundle is freshly energized; the radiative drag becomes strong in the region 10 R_* ≲ r ≲ 50 R_*, where the dimensionless drag parameter D≈(r/c)/t_cool lies between 1 and 100. In this zone the two‑stream instability is most vigorous, producing the radio‑emitting turbulence. As the twist untwists on a timescale of years, the instability weakens, naturally explaining the observed decay of radio emission over similar periods. The model also accounts for the unusually wide pulse profiles (the closed bundle subtends a large fraction of the rotational phase) and the hard radio spectra extending to ∼100 GHz.

Key contributions:

1. First kinetic demonstration that a radiatively‑locked e± flow inevitably drives a sustained two‑stream instability.
2. Identification of superluminal electromagnetic modes capable of escaping the magnetosphere, providing a concrete radio emission channel.
3. Quantitative connection between microphysical turbulence (from Vlasov‑DG simulations) and macroscopic observables (luminosity, spectrum, pulse width, temporal evolution).

The work thus bridges the gap between magnetar magnetospheric electrodynamics and observed radio phenomenology, offering a self‑consistent, long‑lived radio emission mechanism that does not require unrealistically high voltages or open‑field‑line dissipation. It opens avenues for future work, such as detailed polarization predictions, the role of higher multiplicities, and the impact of more complex field geometries on the escaping wave modes.