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.
ABSTRACT
We use the high-resolution TNG50 cosmological magnetohydrodynamical simulation to explore the properties and origin of cold circumgalactic medium (CGM) gas around massive galaxies (M⋆ > 1011 ...M⊙ ) at intermediate redshift ($z \sim 0.5$). We discover a significant abundance of small-scale, cold gas structure in the CGM of ‘red and dead’ elliptical systems, as traced by neutral H i and Mg ii. Halos can host tens of thousands of discrete absorbing cloudlets, with sizes of order a kpc or smaller. With a Lagrangian tracer analysis, we show that cold clouds form due to strong $\delta \rho / \bar{\rho } \gg 1$ gas density perturbations that stimulate thermal instability. These local overdensities trigger rapid cooling from the hot virialized background medium at ∼107 K to radiatively inefficient ∼104 K clouds, which act as cosmologically long-lived, ‘stimulated cooling’ seeds in a regime where the global halo does not satisfy the classic tcool/tff < 10 criterion. Furthermore, these small clouds are dominated by magnetic rather than thermal pressure, with plasma β ≪ 1, suggesting that magnetic fields may play an important role. The number and total mass of cold clouds both increase with resolution, and the mgas ≃ 8 × 104 M⊙ cell mass of TNG50 enables the ∼ few hundred pc, small-scale CGM structure we observe to form. Finally, we make a preliminary comparison against observations from the COS-LRG, LRG-RDR, COS-Halos, and SDSS LRG surveys. We broadly find that our recent, high-resolution cosmological simulations produce sufficiently high covering fractions of extended, cold gas as observed to surround massive galaxies.
Abstract
We present full-volume cosmological simulations, using the moving-mesh code arepo to study the coevolution of dust and galaxies. We extend the dust model in arepo to include thermal ...sputtering of grains and investigate the evolution of the dust mass function, the cosmic distribution of dust beyond the interstellar medium and the dependence of dust-to-stellar mass ratio on galactic properties. The simulated dust mass function is well described by a Schechter fit and lies closest to observations at z = 0. The radial scaling of projected dust surface density out to distances of 10 Mpc around galaxies with magnitudes 17 < i < 21 is similar to that seen in Sloan Digital Sky Survey data, albeit with a lower normalization. At z = 0, the predicted dust density of Ωdust ≈ 1.3 × 10−6 lies in the range of Ωdust values seen in low-redshift observations. We find that the dust-to-stellar mass ratio anticorrelates with stellar mass for galaxies living along the star formation main sequence. Moreover, we estimate the 850 μm number density functions for simulated galaxies and analyse the relation between dust-to-stellar flux and mass ratios at z = 0. At high redshift, our model fails to produce enough dust-rich galaxies, and this tension is not alleviated by adopting a top-heavy initial mass function. We do not capture a decline in Ωdust from z = 2 to 0, which suggests that dust production mechanisms more strongly dependent on star formation may help to produce the observed number of dusty galaxies near the peak of cosmic star formation.
We introduce the Illustris Project, a series of large-scale hydrodynamical simulations of galaxy formation. The highest resolution simulation, Illustris-1, covers a volume of (106.5 Mpc)3, has a dark ...mass resolution of 6.26 × 106 M⊙, and an initial baryonic matter mass resolution of 1.26 × 106 M⊙. At z = 0 gravitational forces are softened on scales of 710 pc, and the smallest hydrodynamical gas cells have an extent of 48 pc. We follow the dynamical evolution of 2 × 18203 resolution elements and in addition passively evolve 18203 Monte Carlo tracer particles reaching a total particle count of more than 18 billion. The galaxy formation model includes: primordial and metal-line cooling with self-shielding corrections, stellar evolution, stellar feedback, gas recycling, chemical enrichment, supermassive black hole growth, and feedback from active galactic nuclei. Here we describe the simulation suite, and contrast basic predictions of our model for the present-day galaxy population with observations of the local universe. At z = 0 our simulation volume contains about 40 000 well-resolved galaxies covering a diverse range of morphologies and colours including early-type, late-type and irregular galaxies. The simulation reproduces reasonably well the cosmic star formation rate density, the galaxy luminosity function, and baryon conversion efficiency at z = 0. It also qualitatively captures the impact of galaxy environment on the red fractions of galaxies. The internal velocity structure of selected well-resolved disc galaxies obeys the stellar and baryonic Tully–Fisher relation together with flat circular velocity curves. In the well-resolved regime, the simulation reproduces the observed mix of early-type and late-type galaxies. Our model predicts a halo mass dependent impact of baryonic effects on the halo mass function and the masses of haloes caused by feedback from supernova and active galactic nuclei.
Abstract
We present the new TNG50 cosmological, magnetohydrodynamical simulation – the third and final volume of the IllustrisTNG project. This simulation occupies a unique combination of large ...volume and high resolution, with a 50 Mpc box sampled by 21603 gas cells (baryon mass of 8 × 104 M⊙). The median spatial resolution of star-forming interstellar medium gas is ∼100−140 pc. This resolution approaches or exceeds that of modern ‘zoom’ simulations of individual massive galaxies, while the volume contains ∼20 000 resolved galaxies with $M_\star \gtrsim 10^7$ M⊙. Herein we show first results from TNG50, focusing on galactic outflows driven by supernovae as well as supermassive black hole feedback. We find that the outflow mass loading is a non-monotonic function of galaxy stellar mass, turning over and rising rapidly above 1010.5 M⊙ due to the action of the central black hole (BH). The outflow velocity increases with stellar mass, and at fixed mass it is faster at higher redshift. The TNG model can produce high-velocity, multiphase outflows that include cool, dense components. These outflows reach speeds in excess of 3000 km s−1 out to 20 kpc with an ejective, BH-driven origin. Critically, we show how the relative simplicity of model inputs (and scalings) at the injection scale produces complex behaviour at galactic and halo scales. For example, despite isotropic wind launching, outflows exhibit natural collimation and an emergent bipolarity. Furthermore, galaxies above the star-forming main sequence drive faster outflows, although this correlation inverts at high mass with the onset of quenching, whereby low-luminosity, slowly accreting, massive BHs drive the strongest outflows.
ABSTRACT
Supermassive black holes (SMBHs) that reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star ...formation in massive galaxies. Here, we study the $10^{9-12.5}\, \mathrm{M_\odot }$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at z = 0, and show how the three components – SMBH, galaxy, and circumgalactic medium (CGM) – are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how the onset of the low-accretion black hole (BH) feedback mode, realized in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star formation quenching at stellar masses $\gtrsim 10^{10.5}\, \mathrm{M_\odot }$ but also to a change in thermodynamic properties of the (non-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10\!-\!100\, \mathrm{Myr}$ to $1\!-\!10\, \mathrm{Gyr}$, effectively ceasing ongoing star formation and inhibiting radiative cooling and future gas accretion. In practice, the same active galactic nucleus (AGN) feedback channel is simultaneously ‘ejective’ and ‘preventative’ and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times.