skip to main content


Title: Efficient early stellar feedback can suppress galactic outflows by reducing supernova clustering
ABSTRACT We present a novel set of stellar feedback models, implemented in the moving-mesh code arepo, designed for galaxy formation simulations with near-parsec (or better) resolution. These include explicit sampling of stars from the IMF, allowing feedback to be linked to individual massive stars, an improved method for the modelling of H ii regions, photoelectric (PE) heating from a spatially varying FUV field and supernova feedback. We perform a suite of 32 simulations of isolated $M_\mathrm{vir} = 10^{10}\, \mathrm{M_\odot }$ galaxies with a baryonic mass resolution of $20\, \mathrm{M_\odot }$ in order to study the non-linear coupling of the different feedback channels. We find that photoionization (PI) and supernova feedback are both independently capable of regulating star formation to the same level, while PE heating is inefficient. PI produces a considerably smoother star formation history than supernovae. When all feedback channels are combined, the additional suppression of star formation rates is minor. However, outflow rates are substantially reduced relative to the supernova only simulations. We show that this is directly caused by a suppression of supernova clustering by the PI feedback, disrupting star-forming clouds prior to the first supernovae. We demonstrate that our results are robust to variations of our star formation prescription, feedback models and the baryon fraction of the galaxy. Our results also imply that the burstiness of star formation and the mass loading of outflows may be overestimated if the adopted star particle mass is considerably larger than the mass of individual stars because this imposes a minimum cluster size.  more » « less
Award ID(s):
1835509
NSF-PAR ID:
10355199
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
506
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3882 to 3915
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$), run to z = 0 with $30\, \mathrm{M}_{\odot }$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs ($M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M*/Mhalo > 10−4 to form a dark matter core ${\gt}200\rm \, pc$, while lower mass UFDs exhibit cusps down to ${\lesssim}100\rm \, pc$, as expected from energetic arguments. Our dwarfs with $M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).

     
    more » « less
  2. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant. 
    more » « less
  3. ABSTRACT

    Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odot }\, {\rm kpc^{-2}}$ (where $\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band $\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. $\Sigma _{1\, {\rm kpc}}$), not ‘total’ galaxy properties.

     
    more » « less
  4. Abstract Type Ia supernovae are critical for feedback and elemental enrichment in galaxies. Recent surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific supernova Ia rate at z ∼ 0 may be ≲ 20 − 50 × higher in lower-mass galaxies than at Milky Way-mass. Independently, observations show that the close-binary fraction of solar-type Milky Way stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of metallicity-dependent rate models on galaxies of $M_* \sim 10^7\, \rm {M}_{\odot }-10^{11}\, \rm {M}_{\odot }$. First, we benchmark our simulated star-formation histories (SFHs) against observations, and show that assumed stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent rate models, potentially causing ASAS-SN and DES observations to agree more than might appear. Models in which the supernova Ia rate increases with decreasing metallicity ($\propto Z^{-0.5 \; \rm {to} \; -1}$) provide significantly better agreement with observations. Encouragingly, these rate increases (≳ 10 × in low-mass galaxies) do not significantly impact galaxy masses and morphologies, which remain largely unaffected except for our most extreme models. We explore implications for both [Fe/H] and [$\alpha /\rm {Fe}$] enrichment; metallicity-dependent rate models can improve agreement with the observed stellar mass-metallicity relations in low-mass galaxies. Our results demonstrate that a range of metallicity-dependent rate models are viable for galaxy formation and motivate future work. 
    more » « less
  5. null (Ed.)
    ABSTRACT Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the gizmo code that follow the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbers found in observed GMCs ($\mathcal {M} \sim 10\!-\!50$). As in past works, the mean and median stellar masses are sensitive to numerical resolution, because they are sensitive to low-mass stars that contribute a vanishing fraction of the overall stellar mass. The mass-weighted median stellar mass M50 becomes insensitive to resolution once turbulent fragmentation is well resolved. Without imposing Larson-like scaling laws, our simulations find $M_\mathrm{50} \,\, \buildrel\propto \over \sim \,\,M_\mathrm{0} \mathcal {M}^{-3} \alpha _\mathrm{turb}\, \mathrm{SFE}^{1/3}$ for GMC mass M0, sonic Mach number $\mathcal {M}$, virial parameter αturb, and star formation efficiency SFE = M⋆/M0. This fit agrees well with previous IMF results from the ramses, orion2, and sphng codes. Although M50 has no significant dependence on the magnetic field strength at the cloud scale, MHD is necessary to prevent a fragmentation cascade that results in non-convergent stellar masses. For initial conditions and SFE similar to star-forming GMCs in our Galaxy, we predict M50 to be $\gt 20 \, \mathrm{M}_{\odot }$, an order of magnitude larger than observed ($\sim 2 \, \mathrm{M}_\odot$), together with an excess of brown dwarfs. Moreover, M50 is sensitive to initial cloud properties and evolves strongly in time within a given cloud, predicting much larger IMF variations than are observationally allowed. We conclude that physics beyond MHD turbulence and gravity are necessary ingredients for the IMF. 
    more » « less