2020MNRAS.494.3663I

Understanding the natal kicks, or birth velocities, of neutron stars is essential for understanding the evolution of massive binaries and double neutron star formation. We use maximum likelihood methods as published in Verbunt et al. to analyse a new large data set of parallaxes and proper motions measured by Deller et al. This sample is roughly three times larger than number of measurements available before. For both the complete sample and its younger part (spin-down ages τ < 3 Myr), we find that a bimodal Maxwellian distribution describes the measured parallaxes and proper motions better than a single Maxwellian with probability of 99.3 and 95.0 per cent, respectively. The bimodal Maxwellian distribution has three parameters: fraction of low-velocity pulsars and distribution parameters σ₁ and σ₂ for low- and high-velocity modes. For a complete sample, these parameters are as follows: 42_–15⁺¹⁷ per cent, σ₁=128_–18⁺²² km s^–1, and σ₂ = 298 ± 28 km s^–1. For younger pulsars, which are assumed to represent the natal kick, these parameters are as follows: 20_–10⁺¹¹ per cent, σ₁=56_–15⁺²⁵ km s^–1, and σ₂ = 336 ± 45 km s^–1. In the young population, 5 ± 3 per cent of pulsars have velocities less than 60 km s^–1. We perform multiple Monte Carlo tests for the method taking into account realistic observational selection. We find that the method reliably estimates all parameters of the natal kick distribution. Results of the velocity analysis are weakly sensitive to the exact values of scale lengths of the Galactic pulsar distribution.