A Self-Consistent Model of Kinetic Alfven Solitons in Pulsar Wind Plasma: Linking Soliton Characteristics to Pulsar Observables

TL;DR

The paper develops a self-consistent kinetic Alfvén soliton model for the pulsar wind zone beyond the light cylinder by deriving a KdV equation from a two-fluid electron–positron–ion plasma with $\kappa$-distributed pairs and oblique KA propagation. Soliton amplitude and width are shown to depend on pulsar observables such as the spin period $P$, spin-down rate $\dot{P}$, and plasma composition (multiplicity $\delta$ and ion loading $\eta$), with the width scaling as $W \propto (P\dot{P})^{1/4}$. A population study of 1174 pulsars reveals a positive correlation between soliton width and spin period, implying distinct soliton flavors across pulsar classes and potential consequences for emission microphysics. The framework links microphysical KA dynamics to macroscopic observables and motivates future tests of soliton-driven processes in pulsar magnetospheres and striped wind regions, with implications for coherent curvature radiation and energy dissipation in low-$\beta$ plasmas.

Abstract

We present a self-consistent model for the formation and propagation of kinetic Alfven (KA) solitons in the pulsar wind zone, where a relativistic, magnetized electron positron ion plasma flows along open magnetic field lines beyond the light cylinder. Using a reductive perturbation approach, we derive a Korteweg de Vries (KdV) equation that governs the nonlinear evolution of KA solitons in this environment. The soliton amplitude and width are shown to depend sensitively on key pulsar observables, including spin period, spin-down rate, and pair multiplicity as well as plasma composition and suprathermal particle distributions. Our analysis reveals that soliton structures are strongly influenced by the presence of heavy ions, kappa-distributed pairs, and oblique propagation angles. Heavier ion species such as Fe26+ produce significantly broader solitons due to enhanced inertia and dispersion, while increasing pair multiplicity leads to smaller solitons through stronger screening. Oblique propagation (larger theta) results in wider but lower-amplitude solitons, and more thermalized pair plasmas (higher kappa) support taller and broader structures. A population-level analysis of 1174 pulsars shows a clear positive correlation between soliton width and spin period, with millisecond pulsars hosting the narrowest solitons. By linking soliton dynamics to measurable pulsar parameters, this work provides a framework for interpreting magnetospheric microphysics and its role in shaping pulsar emission signatures.

Figures (7)

• Figure 1: Schematic representation of the pulsar wind zone beyond the light cylinder, showing the magnetic axis, rotation axis, and the light cylinder at $R_{\rm LC} = c / \Omega$. The green field lines denote the toroidal magnetic field $\mathbf{B}_0$, which dominates in the wind zone. The outflowing electron–positron–ion plasma flows relativistically beyond $R_{\rm LC}$, where KA solitons can form in localized patches (shown in red) due to finite compressibility and weak dispersion. The region shown corresponds to the region, where conditions are favorable for nonlinear soliton formation.
• Figure 2: Local right-handed Cartesian coordinate system showing the magnetic field geometry and KA soliton (KAS) direction. The uniform background magnetic field $\bf B_0$ is aligned along the $z-$axis (toroidal direction), while the soliton propagates obliquely in the $x$-$z$ plane at angle $\theta$ to the field. This frame is co-moving with the relativistic wind, where two-fluid equations apply.
• Figure 3: Effect of ion species on the normalized KA soliton profile $\psi^{(1)}(\zeta)$. The plot compares soliton structures for Hydrogen (H$^+$, green-solid), Helium (He$^{2+}$, red-dotted), and Iron (Fe$^{26+}$, blue-dashed) ions, for the pulsar J1220-6318 with $P=0.7892$, $\dot{P}=7.760 \times10^{-17}$, and keeping all other plasma parameters fixed at $\kappa_e = \kappa_p = 16$, $\alpha = 0.99$, $\theta = 55^{\circ}$ and $u = -0.9$.
• Figure 4: Effect of pair multiplicity $\delta$ and ion loading factor $\eta$ on the normalized KA soliton profile $\psi^{(1)}(\zeta)$ for the same pulsar and baseline parameters as in Fig. \ref{['ions']}. Only $\delta$ and $\eta$ are varied to illustrate their influence on soliton amplitude and width in an ion-rich wind ($Z_i = 26$). Increasing $\delta$ reduces soliton amplitude and width due to enhanced pair screening, while increasing $\eta$ amplifies both via stronger inertial and dispersive effects.
• Figure 5: Effect of propagation angle $\theta$ on the normalized KA soliton profile $\psi^{(1)}(\zeta)$ for the same pulsar and baseline parameters as in Fig. \ref{['ions']}. Increasing $\theta$ from $10^\circ$ to $40^\circ$ leads to a reduction in soliton amplitude and an increase in width, reflecting the weakening of parallel electric fields and reduced charge bunching efficiency at larger obliquities.