2017A&A...608A..57V

We argue that comparison with observations of theoretical models for the velocity distribution of pulsars must be done directly with the observed quantities, that is parallax and the two components of proper motion. We have developed a formalism to do so, and applied it to pulsars with accurate VLBI measurements. For computational convenience, we model the data with Maxwellians. We find that a distribution with two Maxwellians improves significantly on a single Maxwellian. The "mixed" model takes into account that pulsars move away from their place of birth, a narrow region around the Galactic plane. The best model has 42% of the pulsars in a Maxwellian with average velocity {sigma(sqrt(8/π))=120km/s, and 58% in a Maxwellian with average velocity 540km/s. About 5% of the pulsars has a velocity at birth less than 60km/s. For the youngest pulsars (τ_c<10Myr), these numbers are 32% with 130km/s, 68% with 520km/s, and 3%, with appreciable uncertainties. Our analysis shows that the velocity distribution is wider than can be described with a single Maxwellian; it does not prove that two Maxwellians provide a better description than other wide models.