Waves in pulsar winds

The radio, optical, X-ray and gamma-ray nebulae that surround many pulsars are thought to arise from synchrotron and inverse Compton emission. The energy powering this emission, as well as the magnetic fields and relativistic particles, are supplied by a “wind” driven by the central object. The inner parts of the wind can be described using the equations of MHD, but these break down in the outer parts, when the density of charge carriers drops below a critical value. This paper reviews the wave properties of the inner part (striped wind), and uses a relativistic two-fluid model (cold electrons and positrons) to re-examine the nonlinear electromagnetic modes that propagate in the outer parts. It is shown that in a radial wind, two solutions exist for circularly polarised electromagnetic modes. At large distances one of them turns into a freely expanding flow containing a vacuum wave, whereas the other decelerates, corresponding to a confined flow.

💡 Research Summary

The nebular emission that surrounds many rotation‑powered pulsars – from radio through optical and X‑ray to γ‑ray – is powered by a relativistic wind launched by the neutron star. In the inner region of this wind the plasma density is high enough that the charged particles tightly couple the electric and magnetic fields, allowing the flow to be described by ideal magnetohydrodynamics (MHD). In this regime the wind adopts the well‑known “striped” structure: alternating current sheets of electrons and positrons, set up by the oblique rotation of the magnetic dipole, propagate outward with the bulk flow. The MHD description, however, breaks down once the wind expands and the particle number density falls below a critical value n_c ≈ Ω B/(2π e c), where Ω is the pulsar spin frequency and B the magnetic field strength. Below this threshold the plasma can no longer screen the electromagnetic fields, and the dynamics must be treated with a kinetic or multi‑fluid approach.

The authors adopt a cold two‑fluid model, treating electrons and positrons as separate, pressure‑less fluids that interact with the full set of Maxwell equations. By assuming a circularly polarised electromagnetic wave propagating radially, they reduce the problem to a nonlinear dispersion relation that couples the wave number k, angular frequency ω, the plasma dielectric response ε(ω,γ) and magnetic permeability μ(ω,γ), and the Lorentz factor γ of the bulk flow. The resulting quadratic equation in k yields two distinct real solutions.

The first solution satisfies k ≈ ω/c, i.e., the wave behaves essentially as a vacuum electromagnetic wave. In this “free‑expanding” mode the plasma inertia is negligible; the wave carries almost all of the wind’s energy in the fields, and the particles are only weakly perturbed. Energy is transmitted to large radii with little conversion into particle kinetic energy, providing a natural explanation for the high‑energy photons observed far from the pulsar.

The second solution has k < ω/c, indicating a slower phase velocity and a strong coupling between the wave and the plasma. Here the electromagnetic fields do work on the electrons and positrons, transferring a substantial fraction of the Poynting flux into particle kinetic energy. The flow decelerates, and the wave becomes a “confined” mode in which the plasma and fields remain intertwined. The group velocity of this mode can be much smaller than c, leading to a gradual buildup of particle pressure and a modification of the wind’s magnetic structure. The degree of deceleration depends sensitively on the local density and on the initial Lorentz factor of the wind; in the extreme low‑density limit the wave can become almost stationary while still carrying a large field amplitude.

These two modes coexist in the outer wind, and the transition between them is governed by the radial decline of the particle density. In the transition zone non‑linear effects such as wave steepening, parametric decay, and even pair creation may become important, potentially producing the rapid variability and spectral changes observed in pulsar wind nebulae. The authors argue that the free‑expanding mode accounts for the extended high‑energy emission, whereas the confined mode explains regions where the wind appears magnetically dominated and particle acceleration is suppressed.

By quantifying the energy budget of each mode, the paper shows that the free‑expanding solution transports Poynting flux efficiently to large distances, while the confined solution converts a sizable fraction of that flux into particle kinetic energy, thereby influencing the pressure balance and magnetic topology of the nebula. The relative importance of the two solutions depends on the pulsar’s spin‑down power, magnetic field strength, and the ambient medium density.

In summary, the work extends the classic striped‑wind picture by incorporating a relativistic two‑fluid description that remains valid beyond the MHD breakdown radius. It demonstrates that radial, circularly polarised electromagnetic waves admit two physically distinct solutions: a vacuum‑like, freely expanding wave and a decelerating, confined wave. This dual‑mode framework provides a coherent theoretical basis for interpreting the diverse observational signatures of pulsar wind nebulae, and it suggests concrete avenues for future numerical simulations and multi‑wavelength observations to discriminate between the two regimes.