Radio pulsars are neutron stars that receive natal kicks of hundreds of km s−1^{-1} at birth, but whether the kick direction is tied to the stellar spin axis — aligned, orthogonal, isotropic, or mixed — remains unsettled. Biryukov & Beskin (2025, PASA 42, e106) reported the first population-level kinematic evidence for spin–kick alignment: strongly oblique pulsars show larger, more dispersed transverse velocities than weakly oblique ones (25 vs 13 objects; KS p=0.016p=0.016, AD p=0.007p=0.007), as expected if kicks follow spin axes. The result rests on a small subsample with unquantified dispersion-measure distance and proper-motion systematics; a survey-forwarded population synthesis was run only for aligned and isotropic kicks (orthogonal and mixed geometries were explored with a toy Monte Carlo alone); the axial sign of the kick (parallel vs anti-parallel to the spin axis) cannot be fixed by sky-plane projections; and the intrinsic three-dimensional kick-magnitude distribution, including possible multimodality, remains undetermined. Follow-up work through mid-2026 (Li et al. 2026; Lu & Blackman 2026) either assumes the alignment rather than testing it or leaves its facets unresolved, so the geometry, its axial sign, and the 3D kick distribution remain open.

Background

Radio pulsars — spinning, strongly magnetized neutron stars observed through their clock-like radio pulses — are born in core-collapse supernovae and recoil with 'natal kick' velocities of order 100–1000 km s−1^{-1}, thought to arise from asymmetric mass ejection or anisotropic neutrino emission. Three observables enter the present question. The magnetic obliquity α\alpha is the angle between the pulsar's magnetic dipole axis and its rotation axis; it is estimated by fitting the rotating vector model (RVM) — the standard geometric model in which the linear-polarization position angle swings as the rotating beam sweeps past the observer — to polarimetric data. The sky-projected (transverse) velocity vtv_t follows from the pulsar's proper motion and its distance, the latter obtained from parallax measurements or, less reliably, from the dispersion measure (DM), the frequency-dependent pulse delay proportional to the free-electron column density toward the source. Finally, population synthesis — Monte Carlo simulation of the birth, evolution, beam geometry, luminosity, and survey detection of synthetic pulsar populations, in the tradition of Faucher-Giguère & Kaspi — is the standard tool for connecting kick physics to what surveys actually observe.

Whether the natal kick is correlated with the neutron-star spin axis is a decades-old open question. Direct three-dimensional evidence exists for a single object, PSR J0538+2817, whose velocity and spin axis were found to be aligned (best-fit three-dimensional angle near 6°, with a 68% confidence interval extending to about 28°) by combining FAST polarimetry, scintillation analysis, and supernova-remnant geometry (Yao et al. 2021). Statistically, the projected directions of spin axes (from RVM fits) and proper motions are correlated (Johnston et al. 2005; Noutsos et al. 2012, based on 54 pulsars), but the RVM suffers a 90° degeneracy between the two orthogonal polarization modes, so polarization statistics alone cannot distinguish an aligned from an orthogonal kick geometry. The kick-magnitude distribution is itself debated: Hobbs et al. (2005) fitted a single Maxwellian (mean 3D speed near 400 km s−1^{-1}) to 233 proper motions, while several analyses favor a two-component or bimodal distribution. These two questions are coupled: Mandel & Igoshev (2023) showed that assuming spin–kick alignment rather than isotropy substantially changes the kick distribution inferred from observed transverse velocities. Large recent polarimetric surveys — MeerKAT Thousand-Pulsar-Array RVM fits (Johnston et al. 2023) and the FAST polarization database (Wang et al. 2023) — have now supplied obliquity estimates for hundreds of pulsars, making a kinematic test feasible.

Biryukov & Beskin (2025) supplied that test. Under the spin-aligned-kick hypothesis, a pulsar classified as weakly oblique (α\alpha near 0) is seen nearly pole-on, so its spin axis — and hence its kick velocity — points close to the line of sight, producing small, weakly scattered transverse velocities; magnetically orthogonal pulsars instead move nearly in the plane of the sky, so their transverse velocities statistically reflect their full 3D speeds. Comparing 13 weakly and 25 strongly oblique pulsars with known distances and proper motions, they found exactly this split (Kolmogorov–Smirnov p=0.016p=0.016; Anderson–Darling p=0.007p=0.007), reproduced its amplitude in a survey-aware population synthesis (with Rankin's beam-opening-angle formula, a single-normalization radio pseudo-luminosity law, and kick directions drawn within σkick=15∘\sigma_{\rm kick}=15^\circ of the spin axis, following Noutsos et al. 2012) but not under isotropic kicks; a simple Monte Carlo shows an orthogonal kick would produce the same split with the opposite sign. A mixed-alignment population is a population-level mixture rather than a single geometry: a fraction of the natal kicks follow the spin axis within σkick\sigma_{\rm kick} (the aligned case) while the remainder are drawn isotropically, so that the aligned and isotropic limits are its endpoints and the mixing fraction — equivalently, the spread of the spin–kick alignment-angle distribution between those endpoints — is its parameter. The authors themselves flag what remains unresolved: the small kinematic subsample; unquantified DM-distance and proper-motion systematics; unexplored variations of kick prescription, obliquity-evolution model, and survey selection function; the lack of a survey-forwarded synthesis for orthogonal or mixed-alignment populations; the undetermined axial sign of the kick (parallel vs anti-parallel to the spin axis, which projection alone cannot fix); and the undeveloped program of using transverse velocities of young orthogonal rotators — whose vtv_t nearly equals their 3D speed under alignment — to constrain the natal kick-magnitude distribution and its possible multimodality. Theory, for its part, is divided: fallback accretion can generate spin–kick alignment in core-collapse simulations (Janka et al. 2022; Müller 2023), while the gravitational tug-boat mechanism produces alignment only under additional conditions (Lu & Blackman 2026), so the observational geometry is a live discriminator among kick mechanisms.

Problem Statement

Determine the geometric relationship between pulsar natal kicks and spin axes, and exploit that geometry to recover the birth-kick velocity distribution, in three coupled parts that jointly define completion:

(1) Establish whether the obliquity-dependent split in transverse-velocity distributions between weakly and strongly oblique rotators (Biryukov & Beskin 2025: 13 vs 25 objects, KS p=0.016p=0.016, AD p=0.007p=0.007) persists and tightens in enlarged samples — combining large-scale polarimetric obliquity measurements with parallax distances and proper motions, and examining how the signal depends on pulsar age, period, and Galactic location — after the selection and measurement biases of DM-based distances and proper-motion estimates are quantified and corrected and the split is demonstrably not an artifact of the kick-magnitude prescription, the obliquity-evolution model, or the beam-geometry and luminosity assumptions, with orthogonal and mixed-alignment populations examined as alternatives to both the aligned and isotropic cases.

(2) Determine the axial sign of the kick: whether natal kicks are parallel or anti-parallel to the spin axis. Sky-plane projections alone cannot fix the axial sense; resolving it requires genuinely three-dimensional velocity constraints for a sample of pulsars with known spin-axis orientation.

(3) Constrain the intrinsic three-dimensional natal kick-magnitude distribution — in particular whether it is unimodal or multimodal — by exploiting the spin–kick geometry, with quantified projection corrections.

These parts form one objective: the kick geometry controls the projection corrections needed to turn observed transverse velocities into 3D kicks, and the obliquity-split statistics are the observational test of that geometry, so the sign, the geometry class (aligned, orthogonal, isotropic, or mixed), and the recovered kick-magnitude distribution must be mutually consistent. A non-confirming outcome of part (1) — the split does not persist once the selection and measurement biases are corrected — is itself a valid resolution of part (1), not a failure of the task; whether the split persists and whether its significance tightens are separate outcomes to be reported, and tightening is not a precondition for resolution. If the enlarged, bias-corrected sample does not support spin–kick alignment, parts (2) and (3) are to be reported conditionally on the geometry the data actually support, with the axial sign and the projection corrections defined relative to that geometry.

References

  1. Biryukov, A., & Beskin, G. 2025, 'Evidence for the spin-kick alignment of pulsars from the statistics of their magnetic inclinations', Publications of the Astronomical Society of Australia 42, e106, doi:10.1017/pasa.2025.10069, arXiv:2412.12017, https://arxiv.org/abs/2412.12017
  2. Li, Z., Liu, X., You, Z.-Q., Yao, J., & Zhu, X.-J. 2026, 'Constraining the Pulsar 3D Velocity Distribution: The Impact of Spin–Velocity Alignment', The Astrophysical Journal 1000, 288, doi:10.3847/1538-4357/ae4ebf, arXiv:2603.00390, https://arxiv.org/abs/2603.00390
  3. Lu, Y., & Blackman, E. G. 2026, 'Effects of Rotation on the Gravitational Tug-Boat Mechanism for Neutron-Star Kicks and Implications for Spin-Kick Alignment', arXiv:2606.17496, https://arxiv.org/abs/2606.17496
  4. Yao, J. M., Zhu, W. W., Manchester, R. N., Coles, W. A., Li, D., Wang, N., Kramer, M., Stinebring, D. R., Feng, Y., Yan, W. M., Miao, C. C., Yuan, M., Wang, P., & Lu, J. G. 2021, 'Evidence for three-dimensional spin-velocity alignment in a pulsar', Nature Astronomy 5, 788–795, doi:10.1038/s41550-021-01360-w, arXiv:2103.01839, https://arxiv.org/abs/2103.01839
  5. Noutsos, A., Kramer, M., Carr, P., & Johnston, S. 2012, 'Pulsar spin-velocity alignment: further results and discussion', Monthly Notices of the Royal Astronomical Society 423, 2736–2752, doi:10.1111/j.1365-2966.2012.21083.x, arXiv:1205.2305, https://arxiv.org/abs/1205.2305
  6. Johnston, S., Hobbs, G., Vigeland, S., Kramer, M., Weisberg, J. M., & Lyne, A. G. 2005, 'Evidence for alignment of the rotation and velocity vectors in pulsars', Monthly Notices of the Royal Astronomical Society 364, 1397–1412, doi:10.1111/j.1365-2966.2005.09669.x, arXiv:astro-ph/0510260, https://arxiv.org/abs/astro-ph/0510260
  7. Johnston, S., Kramer, M., Karastergiou, A., Keith, M. J., Oswald, L. S., Parthasarathy, A., & Weltevrede, P. 2023, 'The Thousand-Pulsar-Array programme on MeerKAT XI: Application of the rotating vector model', Monthly Notices of the Royal Astronomical Society 520, 4801–4814, doi:10.1093/mnras/stac3636, arXiv:2212.03988, https://arxiv.org/abs/2212.03988
  8. Wang, P. F., Han, J. L., Xu, J., Wang, C., Yan, Y., Jing, W. C., Su, W. Q., Zhou, D. J., & Wang, T. 2023, 'FAST pulsar database: I. Polarization profiles of 682 pulsars', Research in Astronomy and Astrophysics 23, 104002, doi:10.1088/1674-4527/acea1f, arXiv:2307.10340, https://arxiv.org/abs/2307.10340
  9. Mandel, I., & Igoshev, A. P. 2023, 'The impact of spin-kick alignment on the inferred velocity distribution of isolated pulsars', The Astrophysical Journal 944, 153, doi:10.3847/1538-4357/acb3c3, arXiv:2210.12305, https://arxiv.org/abs/2210.12305
  10. Hobbs, G., Lorimer, D. R., Lyne, A. G., & Kramer, M. 2005, 'A statistical study of 233 pulsar proper motions', Monthly Notices of the Royal Astronomical Society 360, 974–992, doi:10.1111/j.1365-2966.2005.09087.x, arXiv:astro-ph/0504584, https://arxiv.org/abs/astro-ph/0504584
  11. Janka, H.-T., Wongwathanarat, A., & Kramer, M. 2022, 'Supernova Fallback as Origin of Neutron Star Spins and Spin-kick Alignment', The Astrophysical Journal 926, 9, doi:10.3847/1538-4357/ac403c, arXiv:2104.07493, https://arxiv.org/abs/2104.07493
  12. Müller, B. 2023, 'Fallback onto Kicked Neutron Stars and its Effect on Spin-Kick Alignment', Monthly Notices of the Royal Astronomical Society 526, 2880–2888, doi:10.1093/mnras/stad2881, arXiv:2308.08312, https://arxiv.org/abs/2308.08312
  13. Wang, X. I., Gao, S.-J., & Li, X.-D. 2025, 'A Two-stage Kick Scenario for the Peculiar Low-mass X-Ray Binary GX 1+4', The Astrophysical Journal Letters 991, L37, doi:10.3847/2041-8213/ae060a, arXiv:2509.21177, https://arxiv.org/abs/2509.21177
  14. Shishkin, D., Bear, E., & Soker, N. 2025, 'Natal Kick by Early-asymmetrical Pairs of Jets to the Neutron Star of Supernova Remnant S147', The Astrophysical Journal 992, 190, doi:10.3847/1538-4357/ae0332, arXiv:2506.21548, https://arxiv.org/abs/2506.21548