Abstract Understanding the interplay of stellar feedback and turbulence in the interstellar medium (ISM) is essential to modeling the evolution of galaxies. To determine the timescales over which stellar feedback drives turbulence in the ISM, we performed a spatially resolved, multiwavelength study of the nearby star-forming dwarf galaxy UGC 4305. As indicators of turbulence on local scales (400 pc), we utilized ionized gas velocity dispersion derived from IFU Hαobservations and atomic gas velocity dispersion and energy surface densities derived from Hisynthesis observations with the Very Large Array. These indicators of turbulence were tested against star formation histories over the past 560 Myr derived from color–magnitude diagrams using Spearman’s rank correlation coefficient. The strongest correlation identified at the 400 pc scale is between measures of Hiturbulence and star formation 70–140 Myr ago. We repeated our analysis of UGC 4305's current turbulence and past star formation activity on multiple physical scales (∼560 and 800 pc) to determine whether there are indications of changes in the correlation timescale with changes to the physical scale. No notable correlations were found at larger physical scales, emphasizing the importance of analyzing star formation-driven turbulence as a local phenomenon.
more »
« less
Testing the Momentum-driven Supernova Feedback Paradigm in M31
Abstract Momentum feedback from isolated supernova remnants (SNRs) have been increasingly recognized by modern cosmological simulations as a resolution-independent means to implement the effects of feedback in galaxies, such as turbulence and winds. However, the integrated momentum yield from SNRs is uncertain due to the effects of SN clustering and interstellar medium (ISM) inhomogeneities. In this paper, we use spatially resolved observations of the prominent 10 kpc star-forming ring of M31 to test models of mass-weighted ISM turbulence driven by momentum feedback from isolated, nonoverlapping SNRs. We use a detailed stellar age distribution (SAD) map from the Panchromatic Hubble Andromeda Treasury survey, observationally constrained SN delay-time distributions, and maps of the atomic and molecular hydrogen to estimate the mass-weighted velocity dispersion using the Martizzi et al. ISM turbulence model. Our estimates are within a factor of two of the observed mass-weighted velocity dispersion in most of the ring, but exceed observations at densities ≲0.2 cm −3 and SN rates >2.1 × 10 −4 SN yr −1 kpc −2 , even after accounting for plausible variations in SAD models and ISM scale height assumptions. We conclude that at high SN rates the momentum deposited is most likely suppressed by the nonlinear effects of SN clustering, while at low densities, SNRs reach pressure equilibrium before the cooling phase. These corrections should be introduced in models of momentum-driven feedback and ISM turbulence.
more »
« less
- Award ID(s):
- 1907790
- PAR ID:
- 10376597
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 928
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 54
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We study star formation-driven outflows in a z ∼ 0.02 starbursting disc galaxy, IRAS08339+6517, using spatially resolved measurements from the Keck Cosmic Web Imager (KCWI). We develop a new method incorporating a multistep process to determine whether an outflow should be fit in each spaxel, and then subsequently decompose the emission line into multiple components. We detect outflows ranging in velocity, vout, from 100 to 600 km s−1 across a range of star formation rate surface densities, ΣSFR, from ∼0.01 to 10 M⊙ yr−1 kpc−2 in resolution elements of a few hundred parsec. Outflows are detected in ∼100 per cent of all spaxels within the half-light radius, and ∼70 per cent within r90, suggestive of a high covering fraction for this starbursting disc galaxy. Around 2/3 of the total outflowing mass originates from the star forming ring, which corresponds to $${\lt}10{{\ \rm per\ cent}}$$ of the total area of the galaxy. We find that the relationship between vout and the ΣSFR, as well as between the mass loading factor, η, and the ΣSFR, are consistent with trends expected from energy-driven feedback models. We study the resolution effects on this relationship and find stronger correlations above a re-binned size-scale of ∼500 pc. Conversely, we do not find statistically significant consistency with the prediction from momentum-driven winds.more » « less
-
Context. Supersonic disordered flows accompany the formation and evolution of molecular clouds (MCs). It has been argued that this is turbulence that can support against gravitational collapse and form hierarchical sub-structures. Aims. We examine the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Are they consistent with uniform turbulence or gravitational collapse? Do the simulated flows agree with observations? Methods. We analysed three MCs that have formed self-consistently within kiloparsec-scale numerical simulations of the interstellar medium (ISM). The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using velocity structure functions (VSFs) with and without density weighting or a density cutoff, and computed in one or three dimensions. However, we do not include optical depth effects that can hide motions in the densest gas, limiting comparison of our results with observations. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence, with a first-order power-law exponent consistent with Larson’s size-velocity relationship. Supernova blast wave impacts on MCs produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Removing the density weighting eliminates this effect as it emphasises the diffuse ISM. Conclusions. We conclude that two different effects coincidentally reproduce Larson’s size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs that are consistent with Larson’s relationship. Later, contraction dominates and the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson’s relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.more » « less
-
Abstract Massive elliptical galaxies harbor large amounts of hot gas (T≳ 106K) in their interstellar medium (ISM) but are typically quiescent in star formation. The jets of active galactic nuclei (AGNs) and Type Ia supernovae (SNe Ia) inject energy into the ISM, which offsets its radiative losses and keeps it hot. SNe Ia deposit their energy locally within the galaxy compared to the larger few ×10 kiloparsec-scale AGN jets. In this study, we perform high-resolution (5123) hydrodynamic simulations of a local (1 kpc3) density-stratified patch of the ISM of massive galaxies. We include radiative cooling and shell-averaged volume heating, as well as randomly exploding SN Ia. We study the effect of different fractions of supernova (SN) heating (with respect to the net cooling rate), different initial ISM density/entropy (which controls the growth timettiof the thermal instability), and different degrees of stratification (which affect the freefall timetff). We find that SNe Ia drive predominantly compressive turbulence in the ISM with a velocity dispersion ofσvup to 40 km s−1and logarithmic density dispersion ofσs∼ 0.2–0.4. These fluctuations trigger multiphase condensation in regions of the ISM, where , in agreement with theoretical expectations that large density fluctuations efficiently trigger multiphase gas formation. Since the SN Ia rate is not self-adjusting, when the net cooling drops below the net heating rate, SNe Ia drive a hot wind which sweeps out most of the mass in our local model. Global simulations are required to assess the ultimate fate of this gas.more » « less
-
ABSTRACT The physics of magnetic fields (B) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components affect galaxy evolution. To understand their common observational tracers, we analyse the magnetic fields in a set of high-resolution, magnetohydrodynamic, cosmological simulations of Milky-Way-like galaxies from the FIRE-2 project. We compare mock observables of magnetic field tracers for simulations with and without CRs to observations of Zeeman splitting and rotation/dispersion measures. We find reasonable agreement between simulations and observations in both the neutral and the ionized interstellar medium (ISM). We find that the simulated galaxies with CRs show weaker ISM |B| fields on average compared to their magnetic-field-only counterparts. This is a manifestation of the effects of CRs in the diffuse, low density inner circumgalactic medium (CGM). We find that equipartition between magnetic and cosmic ray energy densities may be valid at large (> 1 kpc) scales for typical ISM densities of Milky-Way-like galaxies, but not in their haloes. Within the ISM, the magnetic fields in our simulated galaxies follow a power-law scaling with gas density. The scaling extends down to neutral hydrogen number densities < 300 cm−3, in contrast to observationally derived models, but consistent with the observational measurements. Finally, we generate synthetic rotation measure (RM) profiles for projections of the simulated galaxies and compare to observational constraints in the CGM. While consistent with upper limits, improved data are needed to detect the predicted CGM RMs at 10–200 kpc and better constrain theoretical predictions.more » « less