Stellar ejection velocities from the binary supernova scenario: A comparison across population synthesis codes

The vast majority of binary systems are disrupted at the moment of the first supernova, resulting in an unbound compact object and companion star. These ejected companion stars contribute to the observed population of runaway stars. Therefore, an understanding of their ejection velocities is essential to interpreting observations, particularly in the Gaia era of high-precision astronomy. We present a comparison of the predicted ejection velocities of disrupted binary companions in three different population synthesis codes: COSMIC, COMPAS, and binary_c, which use two independent algorithms for the treatment of natal kicks. We confirm that, despite the codes producing different pre-supernova evolution from the same initial conditions, they each find the ejection velocities of secondary stars from disrupted binaries are narrowly distributed about their pre-supernova orbital velocity. We additionally include a correction to the derivation included in Kiel & Hurley 2009 that brings it into agreement with methods from other works for determining post-supernova binary orbital parameters. During this comparison, we identified and resolved bugs in the kick prescriptions of \textit{all three} codes we considered, highlighting how open-science practices and code comparisons are essential for addressing implementation issues.