Star formation and stellar feedback are interlinked processes that redistribute energy, turbulence, and material throughout galaxies. Because young and massive stars form in spatially clustered environments, they create pockets of expanding gas termed superbubbles, which retain information about the physical processes that drive them. As these processes play a critical role in shaping galaxy discs and regulating the baryon cycle, measuring the properties of superbubbles provides important input for galaxy evolution models. With the wide coverage and high angular resolution (~50-150pc) of the PHANGS-ALMA CO (J=2-1) survey, we can now resolve, identify and characterise a statistically representative number of superbubbles using molecular gas in nearby galaxies. We identify superbubbles by requiring spatial correspondence between shells in CO with stellar populations identified in PHANGS-HST. Then, by combining the properties of the stellar populations with the CO, we quantify the energetics of the stars and constrain feedback models. We visually find 325 cavities across 18 PHANGS--ALMA galaxies, 88 of which have clear superbubble signatures (unbroken shells, central clusters, kinematic signatures of expansion). We measure their radii and expansion velocities using CO (2-1) to dynamically derive their ages and the mechanical power driving the bubbles, which we use to compute the expected properties of the parent stellar populations driving the bubbles. We find consistency between the predicted and derived stellar ages and masses of the stellar populations if we use a supernova (SN) model that injects energy with a coupling efficiency of ~10%. Not only does this confirm that molecular gas accurately traces superbubble properties, but it also provides key observational constraints for superbubble models. We also find evidence that the bubbles are sweeping up gas as they expand, and speculate that these sites have the potential to host new generations of stars. This work demonstrates that molecular superbubbles provide novel quantitative constraints on SNe feedback efficiencies and gas clearing times, and represent a promising environment to search for the propagation of star formation, all of which are needed to understand what sets the observed star formation rates in galaxies.

Cone search capability for table J/A+A/676/A67/table2 (Superbubble properties ordered by RA derived using CO)