2025MNRAS.543.1705N

High-mass X-ray binaries (HMXBs) are systems in which a neutron star or black hole accretes material from a massive companion. HMXBs are expected to have experienced a supernova in their evolution. The impulsive kick associated with this event should affect the space velocity of the system in a way that depends on the nature and state of the progenitor binary. Here, we test whether the different evolutionary histories of HMXBs have left a detectable imprint on their peculiar velocities ( $V_{\rm pec}$). Using data from Gaia Data Release 3 (Gaia DR3), we first calculate the $V_{\rm pec}$ values for 63 well-known HMXBs hosting a black hole or neutron star and estimate the associated uncertainties via Monte Carlo resampling. We then analyse their distribution and check for differences between classes. Overall, $V_{\rm pec}$ estimates extend up to 100 km s$^{-1}$, but with Be/X-ray binaries ( BeXRBs) favouring $V_{\rm pec}$ $\lesssim 40$ km s$^{-1}$and supergiant X-ray binaries (SgXRBs) favouring $V_{\rm pec}$ $\gtrsim 40$ km s$^{-1}$. Based on a Kolmogorov-Smirnov (K-S) test, the null hypothesis that the peculiar velocities of both classes are drawn from the same parent distribution can be robustly rejected, irrespective of the background stellar velocity dispersion. Tests with binary population synthesis demonstrate that SgXRBs typically have shorter orbital periods and higher fractional mass loss than BeXRBs at supernova. We argue that the magnitude of $V_{\rm pec}$ could be used as a complementary feature to distinguish between Be and supergiant systems. These findings extend previous inferences based on 2D kinematics from Hipparcos, and may be explained by the differing nature of the respective progenitors systems between the source classes at the instant of supernova.