New stellar bow shocks and bubbles found around runaway stars

Runaway stars with peculiar high velocities can generate stellar bow shocks. Only a few bow shocks show clear radio emission. Our goal is to identify and characterize new stellar bow shocks around O and Be runaway stars in the infrared (IR), and to study their possible radio emission and nature. Our input data is a catalog of O and Be runaways compiled using Gaia DR3. We used WISE IR images to search for bow shocks around these runaways, Gaia DR3 data to determine the actual motion of the runaway stars corrected for interstellar medium (ISM) motion caused by Galactic rotation, and archival radio data to search for emission signatures. We finally explored the radio detectability of these sources under thermal and nonthermal scenarios. We found 9 new stellar bow shock candidates, 3 new bubble candidates, and 1 intermediate structure candidate. One of them is an in situ bow shock candidate. We also found 17 already known bow shocks in our sample, though we discarded one, and 62 miscellaneous sources showing some IR emission around the runaways. We geometrically characterized the sources in IR using the WISE-4 band and estimated the ISM density at the bow shock positions, obtaining median values of ~6 and ~4 cm$^{-3}$ using 2D and 3D peculiar velocities. Most of the new discovered bow shocks come from new runaway discoveries. Within our samples we found that ~24% of the O-type runaway stars show bow shocks, while this decreases to ~3% for Be-type runaway stars. Two bow shocks present radio emission but not as clear counterparts, and two others show hints of radio emission. The physical scenarios indicate that two sources could still be compatible with nonthermal radio emission. The new sample of O and Be runaway stars allowed us to discover both new stellar bow shocks and bubbles. Their geometrical characterization can be used to assess the physical scenario of the radio emission. (Abridged)