Delayed and Displaced: The Impact of Binary Interactions on Core-collapse SN Feedback

Core-collapse supernova feedback models in hydrodynamical simulations typically assume that all stars evolve as single stars. However, the majority of massive stars are formed in binaries and multiple systems, where interactions with a companion can affect stars' subsequent evolution and kinematics. We assess the impact of binary interactions on the timing and spatial distribution of core-collapse supernovae, using `cogsworth` simulations to evolve binary star populations, and their subsequent galactic orbits, within state-of-the-art hydrodynamical zoom-in galaxy simulations. We show that binary interactions: (a) displace supernovae, with ~13% of all supernovae occurring more than 0.1 kpc from their parent cluster; and (b) produce delayed supernovae, such that ~25% of all supernovae occur after the final supernova from a single star population. Delays are largest for low-mass merger products, which can explode more than 200 Myr after a star formation event. We characterize our results as a function of: (1) initial binary population distributions, (2) binary physics parameters and evolutionary pathways, (3) birth cluster dissolution assumptions, and (4) galaxy models (which vary metallicity, star formation history, gravitational potential and simulation codes), and show that the overall timing and spatial distributions of supernovae are surprisingly insensitive to most of these variations. We provide metallicity-dependent analytic fits that can be substituted for single-star subgrid feedback prescriptions in hydrodynamical simulations, and discuss some of the possible implications for binary-driven feedback in galaxies, which may become particularly important at high redshift.