skip to main content


Title: Simulating the interstellar medium and stellar feedback on a moving mesh: implementation and isolated galaxies
ABSTRACT We introduce the Stars and MUltiphase Gas in GaLaxiEs – SMUGGLE model, an explicit and comprehensive stellar feedback model for the moving-mesh code arepo. This novel sub-resolution model resolves the multiphase gas structure of the interstellar medium and self-consistently generates gaseous outflows. The model implements crucial aspects of stellar feedback including photoionization, radiation pressure, energy, and momentum injection from stellar winds and from supernovae. We explore this model in high-resolution isolated simulations of Milky Way like disc galaxies. Stellar feedback regulates star formation to the observed level and naturally captures the establishment of a Kennicutt–Schmidt relation. This result is achieved independent of the numerical mass and spatial resolution of the simulations. Gaseous outflows are generated with average mass loading factors of the order of unity. Strong outflow activity is correlated with peaks in the star formation history of the galaxy with evidence that most of the ejected gas eventually rains down on to the disc in a galactic fountain flow that sustains late-time star formation. Finally, the interstellar gas in the galaxy shows a distinct multiphase distribution with a coexistence of cold, warm, and hot phases.  more » « less
Award ID(s):
1814053 1814259 1817233
NSF-PAR ID:
10175243
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
489
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
4233 to 4260
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    The formation and evolution of galaxies have proved sensitive to the inclusion of stellar feedback, which is therefore crucial to any successful galaxy model. We present INFERNO, a new model for hydrodynamic simulations of galaxies, which incorporates resolved stellar objects with star-by-star calculations of when and where the injection of enriched material, momentum, and energy takes place. INFERNO treats early stellar kinematics to include phenomena such as walkaway and runaway stars. We employ this innovative model on simulations of a dwarf galaxy and demonstrate that our physically motivated stellar feedback model can drive vigorous galactic winds. This is quantified by mass and metal loading factors in the range of 10–100, and an energy loading factor close to unity. Outflows are established close to the disc, are highly multiphase, spanning almost 8 orders of magnitude in temperature, and with a clear dichotomy between mass ejected in cold, slow-moving (T ≲ 5 × 104 K, v < 100 km s−1) gas and energy ejected in hot, fast-moving (T > 106 K, v > 100 km s−1) gas. In contrast to massive disc galaxies, we find a surprisingly weak impact of the early stellar kinematics, with runaway stars having little to no effect on our results, despite exploding in diffuse gas outside the dense star-forming gas, as well as outside the galactic disc entirely. We demonstrate that this weak impact in dwarf galaxies stems from a combination of strong feedback and a porous interstellar medium, which obscure any unique signatures that runaway stars provide.

     
    more » « less
  2. ABSTRACT

    Star-forming galaxies like the Milky Way are surrounded by a hot gaseous halo at the virial temperature – the so-called galactic corona – that plays a fundamental role in their evolution. The interaction between the disc and the corona has been shown to have a direct impact on accretion of coronal gas onto the disc with major implications for galaxy evolution. In this work, we study the gas circulation between the disc and the corona of star-forming galaxies like the Milky Way. We use high-resolution hydrodynamical N-body simulations of a Milky Way-like galaxy with the inclusion of an observationally motivated galactic corona. In doing so, we use SMUGGLE, an explicit interstellar medium (ISM), and stellar feedback model coupled with the moving-mesh code arepo. We find that the reservoir of gas in the galactic corona is sustaining star formation: the gas accreted from the corona is the primary fuel for the formation of new stars, helping in maintaining a nearly constant level of cold gas mass in the galactic disc. Stellar feedback generates a gas circulation between the disc and the corona (the so-called galactic fountain) by ejecting different gas phases that are eventually re-accreted onto the disc. The accretion of coronal gas is promoted by its mixing with the galactic fountains at the disc–corona interface, causing the formation of intermediate temperature gas that enhances the cooling of the hot corona. We find that this process acts as a positive feedback mechanism, increasing the accretion rate of coronal gas onto the galaxy.

     
    more » « less
  3. ABSTRACT

    Arkenstone is a new model for multiphase, stellar feedback-driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealized non-cosmological simulations of a galaxy with a realistic circumgalactic medium (CGM), using the arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the transition region between the interstellar medium (ISM) and the CGM, while also providing refinement at the base of the wind. Without our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the star formation rate of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model.

     
    more » « less
  4. ABSTRACT

    Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionizes and heats the gas, and (2) the depletion of metals on to dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M⊙ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the fire galaxy formation models together with the chimes non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme; and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in H i and H2 mass with stellar mass, and in line luminosity versus star formation rate for [C ii]$_{158 \rm {\mu m}}$, [O i]$_{63 \rm {\mu m}}$, [O iii]$_{88 \rm {\mu m}}$, [N ii]$_{122 \rm {\mu m}}$, and H α6563Å. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [O iii]$_{88 \rm {\mu m}}$ and H α6563Å, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor ≈ 2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions.

     
    more » « less
  5. Abstract We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation (available at http://flathub.flatironinstitute.org/fire ) from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multiphase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Type Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from three suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way (MW)–mass galaxies, five SMC/LMC-mass galaxies, and four lower-mass galaxies including one ultrafaint; we release 39 snapshots across z = 0–10. The second comprises four massive galaxies, with 19 snapshots across z = 1–10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5–10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all MW-mass galaxies, we release the formation coordinates and an “ex situ” flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. 
    more » « less