skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: Galactic coronae in Milky Way-like galaxies: the role of stellar feedback in gas accretion

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
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 4091-4108
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We use FIRE simulations to study disc formation in z ∼ 0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar discs is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation prior to joining the galaxy. Among galaxies with a high fraction ($\gt 70{{\ \rm per\ cent}}$) of their young stars in a thin disc (h/R ∼ 0.1), we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (≳20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to $\lesssim 10^4\, {\rm K}$; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disc, and while cooling transitions from a quasi-spherical spatial configuration to a more-flattened, disc-like configuration. We show that the existence of this ‘rotating cooling flow’ accretion mode is strongly correlated with the fraction of stars forming in a thin disc, using a sample of 17 z ∼ 0 galaxies spanning a halo mass range of 1010.5 M⊙ ≲ Mh ≲ 1012 M⊙ and stellar mass range of 108 M⊙ ≲ M⋆ ≲ 1011 M⊙. Notably, galaxies with a thick disc or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion on to the galaxy are likely a necessary condition for the formation of thin, star-forming disc galaxies in a ΛCDM universe.

    more » « less
  2. 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
  3. Abstract

    Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy’s star-forming efficiency. Mass, energy, and metal outflows (ηM,ηE, andηZ, here normalized by the star formation rate, the SNe energy, and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas, which enters the interstellar medium (ISM), although whether such high mass loadings agree with observations is still unclear. To better understand how the relative importance of ejective (i.e., high mass loading) versus preventative (i.e., high energy loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from 1010to 1012M, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low-mass-loaded (ηM∼ 0.1–10) but high-energy-loaded (ηE∼ 0.1–1) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass loading of outflows but are quite sensitive to our choice of the energy loading, preferringηE∼ 1 for the lowest-mass halos and ∼0.1 for Milky Way–like halos.

    more » « less

    We use accurate estimates of aluminium abundance from the APOGEE Data Release 17 and Gaia Early Data Release 3 astrometry to select a highly pure sample of stars with metallicity −1.5 ≲ [Fe/H] ≲ 0.5 born in-situ in the Milky Way proper. The low-metallicity ([Fe/H]  ≲ −1.3) in-situ component we dub Aurora is kinematically hot with an approximately isotropic velocity ellipsoid and a modest net rotation. Aurora stars exhibit large scatter in metallicity and in many element abundance ratios. The median tangential velocity of the in-situ stars increases sharply with metallicity between [Fe/H] = −1.3 and −0.9, the transition that we call the spin-up. The observed and theoretically expected age–metallicity correlations imply that this increase reflects a rapid formation of the MW disc over ≈1–2 Gyr. The transformation of the stellar kinematics as a function of [Fe/H] is accompanied by a qualitative change in chemical abundances: the scatter drops sharply once the Galaxy builds up a disc during later epochs corresponding to [Fe/H] > −0.9. Results of galaxy formation models presented in this and other recent studies strongly indicate that the trends observed in the MW reflect generic processes during the early evolution of progenitors of MW-sized galaxies: a period of chaotic pre-disc evolution, when gas is accreted along cold narrow filaments and when stars are born in irregular configurations, and subsequent rapid disc formation. The latter signals formation of a stable hot gaseous halo around the MW progenitor, which changes the mode of gas accretion and allows development of coherently rotating disc.

    more » « less
  5. null (Ed.)
    ABSTRACT We examine the 1.4 GHz radio luminosities of galaxies arising from star formation and active galactic nuclei (AGNs) within the state-of-the-art cosmological hydrodynamic simulation Simba. Simba grows black holes via gravitational torque limited accretion from cold gas and Bondi accretion from hot gas, and employs AGN feedback including jets at low Eddington ratios. We define a population of radio loud AGNs (RLAGNs) based on the presence of ongoing jet feedback. Within RLAGN, we define high and low excitation radio galaxies (HERGs and LERGs) based on their dominant mode of black hole accretion: torque limited accretion representing feeding from a cold disc, or Bondi representing advection-dominated accretion from a hot medium. Simba predicts good agreement with the observed radio luminosity function (RLF) and its evolution, overall as well as separately for HERGs and LERGs. Quiescent galaxies with AGN-dominated radio flux dominate the RLF at $\gtrsim 10^{22-23}$ W Hz−1, while star formation dominates at lower radio powers. Overall, RLAGNs have higher black hole accretion rates and lower star formation rates than non-RLAGN at a given stellar mass or velocity dispersion, but have similar black hole masses. Simba predicts an LERG number density of 8.53 Mpc−3, ∼10× higher than for HERGs, broadly as observed. While LERGs dominate among most massive galaxies with the largest black holes and HERGs dominate at high specific star formation rates, they otherwise largely populate similar-sized dark matter haloes and have similar host galaxy properties. Simba thus predicts that deeper radio surveys will reveal an increasing overlap between the host galaxy demographics of HERGs and LERGs. 
    more » « less