skip to main content

Title: Gas infall and radial transport in cosmological simulations of milky way-mass discs

Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, discy galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies at low redshift, we analyse the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (Mhalo ∼ 1012M⊙), focusing on simulations with cosmic ray physics. We find that at z ∼ 0, gas approaches the disc with angular momentum similar to the gaseous disc edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominately parallel to the disc and joins largely in the outskirts. Immediately prior to joining the disc, trajectories briefly become more vertical on average. Within the disc, gas motion is complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show slow net radial infall with transport speeds of 1–3 km s−1 and net mass fluxes through the disc of ∼M⊙ yr−1, comparable to the galaxies’ star more » formation rates and decreasing towards galactic centre as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1–7 km s−1, ∼4–5 M⊙ yr−1). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies.

« less
; ; ; ; ; ; ; ; ;
Award ID(s):
1715101 2108230 1715216 2045928 2108318 1911233 2108314
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 4149-4170
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Observations of ultraviolet (UV) metal absorption lines have provided insight into the structure and composition of the circumgalactic medium (CGM) around galaxies. We compare these observations with the low-redshift (z ≤ 0.3) CGM around dwarf galaxies in high-resolution cosmological zoom-in runs in the FIRE-2 (Feedback In Realistic Environments) simulation suite. We select simulated galaxies that match the halo mass, stellar mass, and redshift of the observed samples. We produce absorption measurements using trident for UV transitions of C iv, O vi, Mg ii, and Si iii. The FIRE equivalent width (EW) distributions and covering fractions for the C iv ion are broadly consistent with observations inside 0.5Rvir, but are underpredicted for O vi, Mg ii, and Si iii. The absorption strengths of the ions in the CGM are moderately correlated with the masses and star formation activity of the galaxies. The correlation strengths increase with the ionization potential of the ions. The structure and composition of the gas from the simulations exhibit three zones around dwarf galaxies characterized by distinct ion column densities: the discy interstellar medium, the inner CGM (the wind-dominated regime), and the outer CGM (the IGM accretion-dominated regime). We find that the outer CGM in the simulations is nearly but not quite supported bymore »thermal pressure, so it is not in hydrostatic equilibrium, resulting in halo-scale bulk inflow and outflow motions. The net gas inflow rates are comparable to the star formation rate of the galaxy, but the bulk inflow and outflow rates are greater by an order of magnitude, with velocities comparable to the virial velocity of the halo. These roughly virial velocities (${\sim } 100 \, \rm km\, s^{-1}$) produce large EWs in the simulations. This supports a picture for dwarf galaxies in which the dynamics of the CGM at large scales are coupled to the small-scale star formation activity near the centre of their haloes.« less

    The globular cluster (GC) systems of low-mass late-type galaxies, such as NGC 2403, have been poorly studied to date. As a low mass galaxy (M*  = 7 × 109 M⊙), cosmological simulations predict NGC 2403 to contain few, if any, accreted GCs. It is also isolated, with a remarkably undisturbed HI disc. Based on candidates from the literature, Sloan Digital Sky Survey and Hyper Suprime-Cam imaging, we selected several GCs for follow-up spectroscopy using the Keck Cosmic Web Imager. From their radial velocities and other properties, we identify eight bona-fide GCs associated with either the inner halo or the disc of this bulgeless galaxy. A stellar population analysis suggests a wide range of GC ages from shortly after the big bang until the present day. We find all of the old GCs to be metal-poor with [Fe/H] ≤ −1. The age–metallicity relation for the observed GCs suggests that they were formed over many Gyr from gas with a low effective yield, similar to that observed in the SMC. Outflows of enriched material may have contributed to the low yield. With a total system of ∼50 GCs expected, our study is the first step in fully mapping the star cluster history of NGC 2403more »in both space and time.

    « less

    We present a detailed study of a galaxy merger taking place at z = 1.89 in the GOODS-S field. Here, we analyse Keck/MOSFIRE spectroscopic observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey along with multiwavelength photometry assembled by the 3D-HST survey. The combined data set is modelled to infer the past star formation histories (SFHs) of both merging galaxies. They are found to be massive, with log10(M*/M⊙) > 11, with a close mass ratio satisfying the typical major-merger definition. Additionally, in the context of delayed-τ models, GOODS-S 43114, and GOODS-S 43683 have similar SFHs and low star formation rates (log10(SFR(SED)/${\rm M}_{\odot }\,\rm {yr}^{-1}$) < 1.0) compared to their past averages. The best-fitting model SEDs show elevated H δA values for both galaxies, indicating that their stellar spectra are dominated by A-type stars, and that star formation peaked ∼0.5−1 Gyr ago and has recently declined. Additionally, based on SED fitting both merging galaxies turned on and shut off star formation within a few hundred Myr of each other, suggesting that their bursts of star formation may be linked. Combining the SFHs and H δA results with recent galaxy merger simulations, we infer that these galaxies have recently completed their first pericentric passage andmore »are moving apart. Finally, the relatively low second velocity moment of GOODS-S 43114, given its stellar mass suggests a disc-like structure. However, including the geometry of the galaxy in the modelling does not completely resolve the discrepancy between the dynamical and stellar masses. Future work is needed to resolve this inconsistency in mass.

    « less

    We study the gas distribution and kinematics of the inner kpc of six moderately luminous (43.43 ≤ log Lbol ≤ 44.83) nearby (0.004 ≤ z ≤ 0.014) Seyfert galaxies observed with the Near-infrared Integral Field Spectrograph (NIFS) in the J ($1.25\,\mu$m) and K ($2.2\,\mu$m) bands. We analyse the most intense emission lines detected on these spectral wavebands: [Fe ii] $1.2570\, \mu$m and Paβ, which trace the ionized gas in the partially and fully ionized regions, and $\mathrm{ H}_2 \ 2.1218\, \mu$m, which traces the hot (∼2000 K) molecular gas. The dominant kinematic component is rotation in the disc of the galaxies, except for the ionized gas in NGC 5899 that shows only weak signatures of a disc component. We find ionized gas outflow in four galaxies, while signatures of H2 outflows are seen in three galaxies. The ionized gas outflows display velocities of a few hundred km s−1, and their mass outflow rates are in the range 0.005–12.49 M⊙ yr−1. Their kinetic powers correspond to 0.005–0.7 per cent of the active galactic nuclei (AGN) bolometric luminosities. Besides rotation and outflows signatures in some cases, the H2 kinematics also reveals inflows in three galaxies. The inflow velocities are 50–80 km s−1 and the mass inflow rates are in the range 1–9 × 10−4 M⊙ yr−1 formore »hot molecular gas. These inflows might be only the hot skin of the total inflowing gas, which is expected to be dominated by colder gas. The mass inflow rates are lower than the current accretion rates to the AGN, and the ionized outflows are apparently disturbing the gas in the inner kpc.

    « less
  5. Of the almost 40 star-forming galaxies at z≳ 5 (not counting quasi-stellar objects) observed in [{{C}} {{II}}] to date, nearly half are either very faint in [{{C}} {{II}}] or not detected at all, and fall well below expectations based on locally derived relations between star formation rate and [{{C}} {{II}}] luminosity. This has raised questions as to how reliable [{{C}} {{II}}] is as a tracer of star formation activity at these epochs and how factors such as metallicity might affect the [{{C}} {{II}}] emission. Combining cosmological zoom simulations of galaxies with SÍGAME (SImulator of GAlaxy Millimeter/submillimeter Emission), we modeled the multiphased interstellar medium (ISM) and its emission in [{{C}} {{II}}], as well as in [O I] and [O III], from 30 main-sequence galaxies at z≃ 6 with star formation rates ˜3-23 {M}⊙ {yr}}-1, stellar masses ˜ (0.7{--}8)× {10}9 {M}⊙ , and metallicities ˜ (0.1{--}0.4)× {Z}⊙ . The simulations are able to reproduce the aforementioned [{{C}} {{II}}] faintness of some normal star-forming galaxy sources at z≥slant 5. In terms of [O I] and [O III], very few observations are available at z≳ 5, but our simulations match two of the three existing z≳ 5 detections of [O III] and are furthermoremore »roughly consistent with the [O I] and [O III] luminosity relations with star formation rate observed for local starburst galaxies. We find that the [{{C}} {{II}}] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ˜66% and ˜27%, respectively. The molecular gas, which constitutes only ˜ 10 % of the total gas mass, is thus a more efficient emitter of [{{C}} {{II}}] than the ionized gas, which makes up ˜85% of the total gas mass. A principal component analysis shows that the [{{C}} {{II}}] luminosity correlates with the star formation activity of a galaxy as well as its average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [{{C}} {{II}}] faintness, and we suggest that these factors may also be responsible for the [{{C}} {{II}}]-faint normal galaxies observed at these early epochs.« less