The "bubbly" interstellar medium as origin for the inhomogeneous internal metallicity distributions in large disk galaxies

Resolved metallicity studies of local disk galaxies have revealed that their interstellar media (ISMs) are far from chemically homogeneous, displaying significant ($\sim 0.05$ dex) variations in the metallicity on characteristic scales of a few hundred parsecs. Such data is at odds with most analytical models, where the ISM is predicted to be more well-mixed. Here, we suggest that the observed small-scale features seen in galaxies may be superbubbles of metal-enriched gas created by a collection of core collapse supernovae with tight spatial (and temporal) correlation. In this scenario, the size of the metallicity fluctuations (superbubble radius, $ϕ$) is set by the disk scale height of the galaxy in question (after which point shock breakout favours preferential expansion along directions perpendicular to the dense disc), and the amount of additional metals contained within a fluctuation is proportional to the star formation efficiency in superbubble regions ($ε$). To test this theory, we analysed metallicity maps from the PHANGS-MUSE sample of galaxies using a geostatistical forward-modelling approach. We find $ϕ\simeq 300$ pc and $ε= 0.1-0.2$, in good agreement with our theoretical model. Further, these small-scale parameters are found to be related to the global galaxy properties, suggesting that the local structure of the interstellar medium of galaxies is not universal. Such a model of star formation paints a new picture of galaxy evolution in the modern universe: in large local galaxies, star formation appears steady and regular when averaged over large scales. However, on small scales, these large galaxies remain intrinsically bursty like their smaller, high-redshift counterparts.

💡 Research Summary

The paper addresses a long‑standing tension between observations of local disk galaxies and analytical models of interstellar medium (ISM) mixing. Integral‑field spectroscopy (e.g., PHANGS‑MUSE) now resolves metallicity maps at ~20 pc scales and reveals that the ISM is not chemically homogeneous: typical variations of ~0.05 dex occur on spatial scales of a few hundred parsecs. Classical models that assume a smoothly distributed star formation history predict far smaller scatter, implying that some additional physical process must maintain these inhomogeneities.

The authors propose that the observed small‑scale metallicity fluctuations are the face‑on signatures of superbubbles—large, hot, metal‑rich cavities created by tightly clustered core‑collapse supernovae (SNe). In their picture, a burst of massive star formation (10–100 high‑mass stars) occurs within a ~10 pc region; over ~1 Myr these stars explode, injecting ~10⁵⁴ erg of energy. The resulting superbubble expands until it reaches the vertical scale height of the cold gas disc (≈200–300 pc). At that point the bubble “breaks out” vertically (the classic “chimney” model), venting hot gas into the circumgalactic medium while leaving behind a metal‑enriched pocket in the disc plane. The radius of the enriched region therefore reflects the disc scale height, and the amount of metals added scales with the star‑formation efficiency ε within the superbubble region.

A simple analytic framework is constructed. Using a Kroupa IMF, the fraction of stars that end as core‑collapse SNe is 0.69 %. To generate a superbubble one needs a star‑forming event converting 1.5–15 × 10³ M⊙ of gas into stars, corresponding to a local SFR of 0.1–1 M⊙ yr⁻¹ in molecular clouds—roughly 10–100 % of a galaxy’s total SFR. The metal yield per unit gas mass is y≈0.015, so the metallicity increase in a superbubble is ΔZ≈ε y. For ε=0.1–0.2 this yields ΔZ≈0.0015, i.e. an observable O/H enhancement of ~0.04 dex, matching the measured fluctuations. The superbubble lifetime (~10⁷ yr) and formation rate (γ_SB≈7×10⁻⁵ yr⁻¹) imply ~700 active superbubbles per Milky‑like galaxy at any time, covering ~20–30 % of the disc area—consistent with observed H II region covering fractions.

To test these predictions, the authors analyse metallicity maps of 19 PHANGS‑MUSE galaxies (stellar masses ~10⁹–10¹¹ M⊙, SFRs 0.1–5 M⊙ yr⁻¹). After standard data reduction (extinction correction, BPT classification, S/N>10), they apply a geostatistical forward‑modelling approach. By fitting the spatial covariance of metallicity residuals with a parametric model that includes a correlation length (the superbubble radius ϕ) and an amplitude (related to ε), they infer ϕ≈300 pc and ε=0.1–0.2 across the sample. These values are in excellent agreement with the analytic expectations. Moreover, the inferred ϕ and ε correlate with global galaxy properties such as stellar mass, star‑formation main‑sequence offset, and disc thickness, indicating that the small‑scale physics of superbubble formation is not universal but modulated by the host galaxy’s macro‑environment.

The paper discusses limitations: (i) strong‑line metallicity diagnostics can be contaminated by variations in ionisation parameter; the authors mitigate this by using three independent diagnostics. (ii) The model treats the superbubble as a simple spherical enrichment zone, ignoring possible anisotropies and mixing processes; (iii) The analysis requires sub‑100 pc resolution, limiting the sample to nearby, well‑observed galaxies. Future work is suggested: high‑resolution ALMA CO observations to map molecular gas within superbubbles, JWST NIRSpec to probe temperature and ionisation structure, and high‑resolution hydrodynamic simulations (e.g., with FLASH or AREPO) to capture the non‑linear interaction of clustered SNe and the resulting breakout dynamics.

In the discussion, the authors argue that their “bubbly ISM” model reconciles the apparent paradox of globally steady star formation (as seen in the star‑formation main sequence) with locally bursty behaviour. Large, present‑day disks appear smooth on kiloparsec scales, yet on sub‑kiloparsec scales they host frequent, energetic superbubble events that locally enrich the gas and drive vertical outflows. This picture aligns high‑redshift, clumpy, bursty star formation with the more quiescent local universe, suggesting a continuity of feedback‑driven processes across cosmic time.

In conclusion, the study provides a compelling, physically motivated explanation for the observed ~0.05 dex metallicity fluctuations in nearby disk galaxies. By linking the characteristic fluctuation scale to the disc scale height and the enrichment amplitude to the local star‑formation efficiency, and by confirming these relations with geostatistical analysis of PHANGS‑MUSE data, the authors establish superbubbles as a dominant driver of ISM chemical inhomogeneity. This framework offers a new lens through which to interpret metallicity maps, informs models of galactic chemical evolution, and highlights the importance of small‑scale feedback in shaping galaxy evolution.