2021MNRAS.508.3345I

Supernova explosion and the associated neutron star (NS) natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary, or a millisecond radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain NS natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass, and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries (BeXs). Our best model for natal kicks is consistent with both measurements and contains a fraction of w = 0.2 ± 0.1 weak natal kicks with σ₁=₄₅+25^_–15 km s^–1, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: σ₂ = 336 km s^–1. We found that currently used models for natal kicks of NSs produced by electron capture supernova (ecSN; combination of Maxwellian σ = 265 km s^–1 and σ = 30 km s^–1 for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and BeXs.