An observational study of rotation and binarity of Galactic O-type runaway stars

Abstract: Gaia DR3 data have revealed new massive runaway stars, while spectroscopic surveys enable detailed characterization. The relative contributions of binary supernova (BSS) and dynamical ejection (DES) scenarios to explain their runaway origin remain poorly constrained, particularly in the Milky Way. We aim to characterize the largest sample of Galactic O-type runaway stars ever investigated through their kinematics, rotation, and binarity to shed light into their origins. We use the GOSC-Gaia DR3 catalog, and IACOB spectroscopic information to build a sample with 214 O-type stars with projected rotational velocities ($v \sin{i}$), and a subsample of 168 O-type stars with additional information about their likely single (LS) or single-lined (SB1) spectroscopic binary nature. We also consider an additional sample of 65 double-lined (SB2) spectroscopic binaries. We find that among our sample of Galactic O-type runaways, most (74%) have $v \sin{i}<200$ km/s, whereas for normal stars this fraction is slightly higher (82%). There are no fast-moving runaways being fast rotators, except for HD 124 979. Runaways show lower SB1 fractions than normal stars, with no runaway SB1 fast-rotating systems; on average, runaways rotate faster than normal stars; and their runaway fraction is higher among fast rotators (44%) vs. the slow rotators (34%). This is consistent with BSS dominance for fast rotators. We also found that SB2 systems hardly reach runaway velocities with a low runaway fraction (10%). Runaways with 2D velocities > 60 km/s are mostly single and interpreted as DES products, while runaways with 2D velocities > 85 km/s are also interpreted as two-step products. Three of 12 runaway SB1 systems are HMXBs. Our study reveals that most Galactic O-type runaways are slow rotators, suggests a dominance of BSS among fast-rotating runaways, and of DES and two-step among the high-velocity ones. (Abridged)