skip to main content

Title: Bursting Bubbles: Clustered Supernova Feedback in Local and High-redshift Galaxies

We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact Hishells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows nearz∼ 2. We find that high-redshift galaxies may “capture” 20%–50% of their feedback momentum in the dense ISM (with the remainder escaping into the nearby CGM), whereas local galaxies may contain ≲10% of their feedback momentum from the central starburst. Using azimuthally averaged galaxy properties, we predict that most superbubbles stall and fragmentwithinthe ISM, and that this occurs at, or near, the gas scale height. We find a consistent interpretation in the observed Hibubble radii and velocities, and predict that most will fragment within the ISM, and that those able to break out originate from short dynamical time regions (where the dynamical time is shorter than feedback timescales). Additionally, we demonstrate that models with constant star cluster formation efficiency per Toomre mass are inconsistent with the occurrence of outflows from high-zstarbursts and local circumnuclear regions.

; ; ;
Publication Date:
Journal Name:
The Astrophysical Journal Letters
Page Range or eLocation-ID:
Article No. L28
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We use medium-resolution Keck/Echellette Spectrograph and Imager spectroscopy of bright quasars to study cool gas traced by Caiiλλ3934, 3969 and Naiλλ5891, 5897 absorption in the interstellar/circumgalactic media of 21 foreground star-forming galaxies at redshifts 0.03 <z< 0.20 with stellar masses 7.4 ≤ logM*/M≤ 10.6. The quasar–galaxy pairs were drawn from a unique sample of Sloan Digital Sky Survey quasar spectra with intervening nebular emission, and thus have exceptionally close impact parameters (R< 13 kpc). The strength of this line emission implies that the galaxies’ star formation rates (SFRs) span a broad range, with several lying well above the star-forming sequence. We use Voigt profile modeling to derive column densities and component velocities for each absorber, finding that column densitiesN(Caii) > 1012.5cm−2(N(Nai) > 1012.0cm−2) occur with an incidencefC(Caii) = 0.63+0.10−0.11(fC(Nai) = 0.57+0.10−0.11). We find no evidence for a dependence offCor the rest-frame equivalent widthsWr(CaiiK) orWr(Nai5891) onRorM*. Instead,Wr(CaiiK) is correlated with local SFR at >3σsignificance, suggesting that Caiitraces star formation-driven outflows. While most of the absorbers have velocities within ±50 km s−1of the host redshift, their velocity widths (characterized by Δv90) are universally 30–177 km s−1larger than that implied by tilted-ring modeling of the velocities of interstellar material. These kinematics mustmore »trace galactic fountain flows and demonstrate that they persist atR> 5 kpc. Finally, we assess the relationship between dust reddening andWr(CaiiK) (Wr(Nai5891)), finding that 33% (24%) of the absorbers are inconsistent with the best-fit Milky WayE(B−V)-Wrrelations at >3σsignificance.

    « less
  2. Abstract

    Stellar feedback is fundamental to the modeling of galaxy evolution, as it drives turbulence and outflows in galaxies. Understanding the timescales involved are critical for constraining the impact of stellar feedback on the interstellar medium. We analyzed the resolved star formation histories along with the spatial distribution and kinematics of the atomic and ionized gas of four nearby star-forming dwarf galaxies (NGC 4068, NGC 4163, NGC 6789, and UGC 9128) to determine the timescales over which stellar feedback drives turbulence. The four galaxies are within 5 Mpc and have a range of properties including current star formation rates of 0.0005–0.01Myr−1, log(M*/M) between 7.2 and 8.2, and log(MHi/M) between 7.2 and 8.3. Their color–magnitude diagram derived star formation histories over the past 500 Myr were compared to their atomic and ionized gas velocity dispersion and Hienergy surface densities as indicators of turbulence. The Spearman’s rank correlation coefficient was used to identify any correlations between their current turbulence and their past star formation activity on local scales (∼400 pc). The strongest correlation found was between the Hiturbulence measures and the star formation rate 100–200 Myr ago. This suggests a coupling between the star formation activity and atomic gas on this timescale.more »No strong correlation between the ionized gas velocity dispersion and the star formation activity between 5 and 500 Myr ago was found. The sample and analysis are the foundation of a larger program aimed at understanding the timescales over which stellar feedback drives turbulence.

    « less
  3. Context. The Lyman- α line in the ultraviolet (UV) and the [CII] line in the far-infrared (FIR) are widely used tools to identify galaxies in the early Universe and to obtain insights into interstellar medium (ISM) properties in high-redshift galaxies. By combining data obtained with ALMA in band 7 at ∼320 GHz as part of the ALMA Large Program to INvestigate [CII] at Early Times (ALPINE) with spectroscopic data from DEIMOS at the Keck Observatory, VIMOS and FORS2 at the Very Large Telescope, we assembled a unique sample of 53 main-sequence star-forming galaxies at 4.4 <   z  <  6 in which we detect both the Lyman- α line in the UV and the [CII] line in the FIR. Aims. The goal of this paper is to constrain the properties of the Ly α emission in these galaxies in relation to other properties of the ISM. Methods. We used [CII], observed with ALMA, as a tracer of the systemic velocity of the galaxies, and we exploited the available optical spectroscopy to obtain the Ly α -[CII] and ISM-[CII] velocity offsets. Results. We find that 90% of the selected objects have Ly α -[CII] velocity offsets in the range 0 <  Δ vmore »Ly α  − [CII]  <  400 km s −1 , in line with the few measurements available so far in the early Universe, and significantly smaller than those observed at lower redshifts. At the same time, we observe ISM-[CII] offsets in the range −500 <  Δ v ISM−[CII]  <  0 km s −1 , in line with values at all redshifts, which we interpret as evidence for outflows in these galaxies. We find significant anticorrelations between Δ v Ly α −[CII] and the Ly α rest-frame equivalent width EW 0 (Ly α ) (or equivalently, the Ly α escape fraction f esc (Ly α )): galaxies that show smaller Δ v Ly α −[CII] have larger EW 0 (Ly α ) and f esc (Ly α ). Conclusions. We interpret these results in the framework of available models for the radiative transfer of Ly α photons. According to the models, the escape of Ly α photons would be favored in galaxies with high outflow velocities, producing large EW 0 (Ly α ) and small Δ v Ly α -[CII] , in agreement with our observations. The uniform shell model would also predict that the Ly α escape in galaxies with slow outflows (0 <   v out  <  300 km s −1 ) is mainly determined by the neutral hydrogen column density (NHI) along the line of sight, while the alternative model by Steidel et al. (2010, ApJ, 717, 289) would more highly favor a combination of NHI at the systemic velocity and covering fraction as driver of the Ly α escape. We suggest that the increase in Ly α escape that is observed in the literature between z  ∼ 2 and z  ∼ 6 is not due to a higher incidence of fast outflows at high redshift, but rather to a decrease in average NHI along the line of sight, or alternatively, a decrease in HI covering fraction.« less
  4. Abstract

    We present the discovery of neutral gas detected in both damped Lyαabsorption (DLA) and Hi21 cm emission outside of the stellar body of a galaxy, the first such detection in the literature. A joint analysis between the Cosmic Ultraviolet Baryon Survey and the MeerKAT Absorption Line Survey reveals an Hibridge connecting two interacting dwarf galaxies (log (Mstar/M) = 8.5 ± 0.2) that host az= 0.026 DLA with log[N(Hi)/cm−2] = 20.60 ± 0.05 toward the QSO J2339−5523 (zQSO= 1.35). At impact parameters ofd= 6 and 33 kpc, the dwarf galaxies have no companions more luminous than ≈0.05L*within at least Δv= ±300 km s−1andd≈ 350 kpc. The Hi21 cm emission is spatially coincident with the DLA at the 2σ–3σlevel per spectral channel over several adjacent beams. However, Hi21 cm absorption is not detected against the radio-bright QSO; if the background UV and radio sources are spatially aligned, the gas is either warm or clumpy (with a spin temperature to covering factor ratioTs/fc> 1880 K). Observations with VLT-MUSE demonstrate that theα-element abundance of the ionized interstellar medium (ISM) is consistent with the DLA (≈10% solar), suggesting that the neutral gas envelope is perturbed ISM gas. This study showcases the impact of dwarf–dwarfmore »interactions on the physical and chemical state of neutral gas outside of star-forming regions. In the SKA era, joint UV and Hi21 cm analyses will be critical for connecting the cosmic neutral gas content to galaxy environments.

    « less
  5. Abstract

    We use a sample of 27 gamma-ray bursts (GRBs) at redshiftz= 2–6 to probe the outflows in their respective host galaxies (log(M*/M) ∼ 9–11) and search for possible relations between the outflow properties and those of the host galaxies, such asM*, the star formation rate (SFR), and the specific SFR (sSFR). First, we consider three outflow properties: outflow column density (Nout), maximum outflow velocity (Vmax), and normalized maximum velocity (Vnorm=Vmax/Vcirc,halo, whereVcirc,halois the halo circular velocity). We observe clear trends ofNoutandVmaxwith increasing SFR in high-ion-traced outflows, with a stronger (>3σ)Vmax–SFR correlation. We find that the estimated mass outflow rate and momentum flux of the high-ion outflows scale with SFR and can be supported by the momentum imparted by star formation (supernovae and stellar winds). The kinematic correlations of high-ion-traced outflows with SFR are similar to those observed for star-forming galaxies at low redshifts. The correlations with SFR are weaker in low-ion outflows. This, along with the lower detection fraction in low-ion outflows, indicates that the outflow is primarily high-ion dominated. We also observe a strong (>3σ) trend of normalized velocity (Vnorm) decreasing with halo mass and increasing with sSFR, suggesting that outflows from low-mass halos and high-sSFR galaxies are mostmore »likely to escape and enrich the outer circumgalactic medium (CGM) and intergalactic medium with metals. By comparing the CGM–GRB stacks with those of starbursts atz∼ 2 andz∼ 0.1, we find that over a broad redshift range, the outflow strength strongly depends on the main-sequence offset at the respective redshifts, rather than simply the SFR.

    « less