The paucity of old, isolated accreting neutron stars in ROSAT observations is used to derive a lower limit on the mean velocity of neutron stars at birth. The secular evolution of the population is simulated following the paths of a statistical sample of stars for different values of the initial kick velocity, drawn from an isotropic, Gaussian distribution with mean velocity 0 ≤ 〈V〉 ≤ 550 km s−1. The spin-down, induced by dipole losses and the interaction with the ambient medium, is tracked together with the dynamical evolution in the Galactic potential, allowing for the determination of the fraction of stars which are, at present, in each of the four possible stages: Ejector, Propeller, Accretor, and Georotator. Taking from the ROSAT All-Sky Survey an upper limit of ~ 10 accreting neutron stars within ~ 140 pc from the Sun, we infer a lower bound for the mean kick velocity, 〈V〉 ≳ 200–300 km s−1. The same conclusion is reached for both a constant (B ~ 1012 G) and an exponentially decaying magnetic field with a timescale ~ 109 yr. Present results, moreover, constrain the fraction of low-velocity stars which could have escaped pulsar statistics to ≲ 1%.
Population Synthesis of Old Neutron Stars in the Galaxy