Linearly polarized superluminal waves in pulsar winds

Pulsar winds are the ideal environment for the study of non-linear electromagnetic waves. It is generally thought that a pulsar launches a striped wind, a magnetohydrodynamic entropy wave, where plasma sheets carried along with the flow separate regions of alternating magnetic field. But when the density drops below a critical value, or equivalently for distances from the pulsar greater than a critical radius, a strong superluminal wave can also propagate. In this contribution we discuss the conversion of the equatorial striped wind into a linearly polarized superluminal wave, and we argue that this mode is important for the conversion of Poynting flux to kinetic energy flux before the outflow reaches the termination shock.

💡 Research Summary

Pulsar winds constitute a unique laboratory for studying non‑linear electromagnetic phenomena in relativistic plasmas. The conventional picture describes the outflow as a “striped wind”: an MHD entropy wave in which alternating toroidal magnetic fields are separated by thin plasma sheets (current sheets) that corotate with the pulsar. This structure naturally arises from the oblique rotator geometry and persists as long as the plasma density remains sufficiently high to enforce the frozen‑in condition.

The authors point out that when the radial distance from the neutron star exceeds a critical value, the plasma density drops below a threshold n_c (or equivalently the plasma frequency ω_p becomes smaller than the characteristic wave frequency ω). In this low‑density regime the dielectric response of the medium satisfies ε = 1 − ω_p²/ω² < 1, allowing electromagnetic disturbances to propagate with a phase velocity v_ph = c/√ε > c. Such superluminal modes are fundamentally different from the sub‑luminal Alfvénic or magnetosonic waves of ideal MHD.

Focusing on linearly polarized superluminal waves, the paper derives the non‑linear wave equation by coupling Maxwell’s equations to the relativistic fluid equations for electrons and ions. The analysis shows that the electron inertia dominates the response, while the heavy ions remain essentially stationary on the wave timescale. The resulting dispersion relation admits a branch with v_ph > c provided the wave amplitude is sufficiently large to keep the plasma “transparent” to the wave. The authors identify a set of dimensionless parameters – the magnetisation σ, the ratio of plasma frequency to wave frequency, and the temperature‑induced pressure anisotropy – that control the stability and growth of this branch.

Through a series of one‑dimensional numerical experiments, the authors simulate the conversion of an equatorial striped wind into a superluminal wave as the flow expands beyond the critical radius r_c (≈10⁹–10¹⁰ cm for typical young pulsars). Initially the wind exhibits the classic alternating B‑field pattern with embedded current sheets. As the density falls below n_c, the current sheets become unstable to the superluminal mode: the wave penetrates the sheets, smooths out the magnetic reversal, and accelerates the electrons in the direction of propagation. The magnetic energy (Poynting flux) is progressively transferred to particle kinetic energy. Quantitatively, the conversion efficiency reaches 30–50 % of the initial Poynting flux before the flow reaches the termination shock.

The paper argues that this pre‑shock conversion has profound implications for the long‑standing σ‑problem (the discrepancy between the highly magnetised wind near the pulsar and the weakly magnetised nebula downstream of the termination shock). If a substantial fraction of the electromagnetic energy is already deposited into particles by the superluminal wave, the termination shock need not perform the bulk of the conversion, alleviating the requirement for extreme dissipation at the shock itself.

Observationally, the presence of a linearly polarized superluminal wave could imprint distinctive signatures on the high‑energy emission from pulsar wind nebulae. The accelerated electrons would emit synchrotron radiation extending to X‑ray and γ‑ray energies, potentially explaining the hard spectral tails observed in several young nebulae. Moreover, the linear polarization of the wave may be reflected in the polarization degree and angle of the emitted radiation, offering a diagnostic tool for future polarimetric missions.

In summary, the authors provide a rigorous theoretical framework for the emergence of linearly polarized superluminal waves in low‑density pulsar winds, demonstrate their role in converting Poynting flux to particle kinetic energy, and highlight the relevance of this mechanism for both the dynamics of the wind and the observable high‑energy phenomenology of pulsar wind nebulae. Further work, especially multi‑dimensional kinetic simulations and direct comparison with polarimetric observations, will be essential to validate and extend these findings.