Origin of radio polarization in pulsar polar caps

It is crucial to know the polarization properties of coherent radio waves that escape from pulsar polar caps to calculate the radiative transfer through the magnetosphere and to predict observable radio properties. We describe pair cascades in the pulsar polar cap, and we determine for the first time the Stokes parameters of the escaping radio waves from first-principle kinetic simulations for a pulsar with a magnetic obliquity of $60^{\circ}$. We present 3D particle-in-cell kinetic simulations that include quantum-electrodynamic pair cascades in a charge-limited flow from the stellar surface. Our model quantitatively and qualitatively explains the observed pulsar radio powers and spectra, the pulse profiles, polarization curves, their temporal variability, the strong Stokes-$L$ and weak Stokes-$V$ polarization components, the decline in the linear polarization with frequency, and the nonexistence of a radius-to-frequency relation. The observable properties of radio emission from the polar cap can vary and include single- or double-peaked profiles. Most of the Stokes~$V$ curves from our simulations appear to be antisymmetric, but symmetric curves are also present at some viewing angles. Although the polarization-angle (PA) swing of the radiation from the polar cap fits the rotating vector model (RVM) for most viewing angles, the angles obtained from the RVM do not correspond to the dipole geometry of the magnetic field. Instead, the PA is directly related to the plasma flows in the polar cap. Our simulations demonstrate that pair discharges close to the surface of the polar cap cause the radio emission of pulsars and determine the majority of their typically observed properties. The merits of RVM for estimations of the magnetic field geometry from observations need to be reevaluated.

💡 Research Summary

This paper presents the first comprehensive calculation of the Stokes parameters of pulsar radio emission directly from three‑dimensional particle‑in‑cell (PIC) simulations that incorporate quantum electrodynamic (QED) pair cascades in a charge‑limited flow above the polar cap. The authors model a neutron star with a 10 km radius, 1.5 M⊙ mass, a rotation period of 0.25 s, and a dipolar magnetic field of 10¹² G inclined by 60° to the rotation axis. The simulation domain spans roughly 1.2 km × 2.5 km × 1.2 km with a grid resolution of 1.67 m, allowing the electric gap height (~100 m) to be well resolved. Magnetospheric currents are imposed using analytic fits derived from global force‑free simulations (Gralla et al. 2017), ensuring that the Goldreich‑Julian charge density is reproduced at the stellar surface.

Pair creation proceeds via acceleration of primary particles in the gap, curvature‑radiated γ‑photons, and subsequent magnetic conversion into electron‑positron pairs. The resulting non‑stationary plasma exhibits strong spatial inhomogeneities: low‑density “Poynting‑flux channels” form along magnetic field lines where the current density is small, while dense plasma bunches occupy the surrounding regions. Electromagnetic waves generated by the intermittent charge separation are initially superluminal O‑mode disturbances. When these waves encounter the low‑density channels, their frequencies exceed the local plasma frequency, allowing them to decouple and propagate as vacuum‑like radio waves.

The authors extract the electric field at a virtual plane near the top of the simulation box, Fourier‑transform it, and compute Stokes I, Q, U, and V. The key observational signatures emerging from the simulations are:

1. Strong linear polarization – the linear component L = √(Q² + U²) carries 60–80 % of the total intensity across the simulated frequency band (≈100 MHz–few GHz).
2. Weak circular polarization – Stokes V remains at the 5–10 % level, with most V‑profiles showing an antisymmetric shape about the pulse centre; symmetric V‑curves appear for certain viewing geometries.
3. Frequency dependence – L systematically declines with increasing frequency, while V is essentially flat, reproducing the observed trend of decreasing linear polarization at higher radio frequencies.
4. Pulse morphology – depending on the observer’s line‑of‑sight relative to the current channels, the simulated pulse profile can be single‑peaked or double‑peaked, with the separation set by the spacing of the channels.
5. Polarization‑angle (PA) swing – the PA follows a rotating‑vector‑model (RVM)‑like S‑shaped curve for most angles, yet the fitted magnetic inclination and viewing angles derived from the RVM do not match the actual dipole geometry. Instead, the PA is directly linked to the direction of plasma flows within the channels.

These results challenge several long‑standing assumptions in pulsar radio theory. The absence of a radius‑to‑frequency mapping follows naturally because the waves are generated near the surface and then travel through low‑density channels without significant refraction or absorption; thus the emission height does not vary with observing frequency. The traditional use of the RVM to infer magnetic geometry is called into question, as the PA is governed by plasma dynamics rather than the static dipole field.

The study also clarifies why observed pulsar radio emission often shows strong linear but weak circular polarization, and why the linear fraction drops at higher frequencies: the low‑density channels become increasingly transparent at higher frequencies, reducing the differential propagation effects that sustain high linear polarization.

Limitations are acknowledged. The plasma density and skin depth are down‑scaled to make the problem computationally tractable, potentially under‑estimating wave‑plasma coupling. The simulation box does not extend to the full magnetosphere, so propagation effects beyond the polar cap (e.g., birefringence, mode coupling in the outer magnetosphere) are not captured. General relativistic light‑bending and aberration are omitted, which could modify the observed pulse shapes for fast rotators.

In conclusion, the paper demonstrates that non‑stationary pair discharges in the polar cap, together with the formation of low‑density Poynting‑flux channels, can account for the majority of observed radio properties of pulsars. The work provides the first first‑principles prediction of the full set of Stokes parameters from a realistic 3‑D kinetic model, and it urges a re‑evaluation of the rotating‑vector model as a tool for magnetic geometry diagnostics. Future work should aim at larger domains, higher‑fidelity plasma scaling, and inclusion of full general‑relativistic radiative transfer to bridge the gap between the polar‑cap emission and the observed radio signal at Earth.