Neutron star kicks and their relationship to supernovae ejecta mass

We propose a simple model to explain the velocity of young neutron stars. We attempt to confirm a relationship between the amount of mass ejected in the formation of the neutron star and the 'kick' velocity imparted to the compact remnant resulting from the process. We assume that the velocity is given by v(kick) = alpha (M-ejecta/M-remnant) + beta. To test this simple relationship, we use the BPASS (Binary Population and Spectral Synthesis) code to create stellar population models from both single and binary star evolutionary pathways. We then use our Remnant Ejecta and Progenitor Explosion Relationship (REAPER) code to apply different alpha and beta values, and three different 'kick' orientations then record the resulting velocity probability distributions. We find that while a single star population provides a poor fit to the observational data, the binary population provides an excellent fit. Values of alpha = 70 km s(-1) and beta = 110 km s(-1) reproduce the Hobbs et al. observed two-dimensional velocities, and alpha = 70 km s(-1) and beta = 120 km s(-1) reproduce their inferred three-dimensional velocity distribution for nearby single neutron stars with ages less than 3 Myr. After testing isotropic, spin-axis aligned and orthogonal to spin-axis 'kick' orientations, we find no statistical preference for a 'kick' orientation. While ejecta mass cannot be the only factor that determines the velocity of supernova compact remnants, we suggest that it is a significant contributor and that the ejectabased 'kick' should replace the Maxwell-Boltzmann velocity distribution currently used in many population synthesis codes.