Search

The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the MBH − M⋆ relation and the BH mass function. All simulations predict tight MBH − M⋆ relations, and struggle to produce BHs of $M_{\rm BH}\leqslant 10^{7.5}\, \rm M_{\odot }$ in galaxies of $M_{\star }\sim 10^{10.5}\!-\!10^{11.5}\, \rm M_{\odot }$. While the time evolution of the mean MBH − M⋆ relation is mild ($\rm \Delta M_{\rm BH}\leqslant 1\, dex$ for 0 $\leqslant z \leqslant$ 5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for $M_{\star }\leqslant 10^{10.5}\, \rm M_{\odot }$. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with observations of the local Universe, we find an excess of BHs with $M_{\rm BH}\geqslant 10^{9}\, \rm M_{\odot }$ in most of the simulations. The BH mass function is dominated by efficiently accreting BHs ($\log _{10}\, f_{\rm Edd}\geqslant -2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.

© 2018 The Author(s). We use the IllustrisTNG simulations to investigate the evolution of the mass-metallicity relation (MZR) for star-forming cluster galaxies as a function of the formation history of their cluster host. The simulations predict an enhancement in the gas-phase metallicities of star-forming cluster galaxies (109 < M* < 1010M⊙ h-1) at z ≤ 1.0 in comparisons to field galaxies. This is qualitatively consistent with observations. We find that the metallicity enhancement of cluster galaxies appears prior to their infall into the central cluster potential, indicating for the first time a systematic 'chemical pre-processing' signature for infalling cluster galaxies. Namely, galaxies that will fall into a cluster by z = 0 show a ~0.05 dex enhancement in the MZR compared to field galaxies at z ≤ 0.5. Based on the inflow rate of gas into cluster galaxies and its metallicity, we identify that the accretion of pre-enriched gas is the key driver of the chemical evolution of such galaxies, particularly in the stellar mass range (109 < M* < 1010M⊙ h-1). We see signatures of an environmental dependence of the ambient/inflowing gas metallicity that extends well outside the nominal virial radius of clusters. Our results motivate future observations looking for pre-enrichment signatures in dense environments.

In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modelling of feedback from active galactic nuclei (AGNs). The amount of energy released by AGN feedback is proportional to the mass that has been accreted on to the black holes (BHs), but the exact subgrid modelling of AGN feedback differs in all simulations. While modern simulations reliably produce populations of quiescent massive galaxies at z ≤ 2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN populations of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z = 0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGNs of $L_{\rm x, 2\!-\!10 \, keV}\sim 10^{43\!-\!44}\, \rm erg\, s^{-1}$. AGN feedback in some simulations prevents the existence of any bright AGN with $L_{\rm x, 2\!-\!10 \, keV}\geqslant 10^{45}\rm \,erg\, s^{-1}$ (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z ≤ 2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{43-45}\, \rm erg\, s^{-1}$) and fainter ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{42-43}\, \rm erg\, s^{-1}$) AGNs and at both low and high redshifts. These differences can aid us in improving future BH and galaxy subgrid modelling in simulations. Upcoming X-ray missions (e.g. Athena, AXIS, and LynX) will bring faint AGNs to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGNs in future surveys spans at least one order of magnitude across the simul

In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modelling of feedback from active galactic nuclei (AGNs). The amount of energy released by AGN feedback is proportional to the mass that has been accreted on to the black holes (BHs), but the exact subgrid modelling of AGN feedback differs in all simulations. While modern simulations reliably produce populations of quiescent massive galaxies at z ≤ 2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN populations of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z = 0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGNs of $L_{\rm x, 2\!-\!10 \, keV}\sim 10^{43\!-\!44}\, \rm erg\, s^{-1}$. AGN feedback in some simulations prevents the existence of any bright AGN with $L_{\rm x, 2\!-\!10 \, keV}\geqslant 10^{45}\rm \,erg\, s^{-1}$ (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z ≤ 2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{43-45}\, \rm erg\, s^{-1}$) and fainter ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{42-43}\, \rm erg\, s^{-1}$) AGNs and at both low and high redshifts. These differences can aid us in improving future BH and galaxy subgrid modelling in simulations. Upcoming X-ray missions (e.g. Athena, AXIS, and LynX) will bring faint AGNs to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGNs in future surveys spans at least one order of magnitude across the simul

© 2019 The Author(s). We select galaxies from the IllustrisTNG hydrodynamical simulations (Mstars109Mo at 0 ≤ z ≤ 2) and characterize the shapes and evolutions of their UVJ and star formation rate Cstellar mass (SFR-CMstars) diagrams. We quantify the systematic uncertainties related to different criteria to classify star-forming versus quiescent galaxies, different SFR estimates, and by accounting for the star formation measured within different physical apertures. The TNG model returns the observed features of the UVJ diagram at z ≤ 2, with a clear separation between two classes of galaxies. It also returns a tight star-forming main sequence (MS) for Mstars < 1010.5 (Mo) with a ∼0.3 dex scatter at z ∼ 0 in our fiducial choices. If a UVJbased cut is adopted, the TNG MS exhibits a downwardly bending at stellar masses of about 1010.5-10.7 Mo. Moreover, the model predicts that ∼80 (50) per cent of 1010.5-11Mo galaxies at z = 0 (z = 2) are quiescent and the numbers of quenched galaxies at intermediate redshifts and high masses are in better agreement with observational estimates than previous models. However, shorter SFR-averaging time-scales imply higher normalizations and scatter of the MS, while smaller apertures lead to underestimating the galaxy SFRs: Overall we estimate the inspected systematic uncertainties to sum up to about 0.2.0.3 dex in the locus of the MS and to about 15 percentage points in the fraction of quenched galaxies. While TNG colour distributions are clearly bimodal, this is not the case for the SFR logarithmic distributions in bins of stellar mass (SFR ≥ 10.3 Moy-1). Finally, the slope and z = 0 normalization of the TNG MS are consistent with observational findings; however, the locus of the TNG MS remains lower by about 0.2.0.5 dex at 0.75 ≤z <2 than the available observational estimates taken at face value.

© 2018 The Author(s). Published by Oxford University Press on behalf of the Royal Astronomical Society. The distribution of elements in galaxies provides a wealth of information about their production sites and their subsequent mixing into the interstellar medium. Here we investigate the elemental distributions of stars in the IllustrisTNG simulations. We analyse the abundance ratios of magnesium and europium in Milky Way-like galaxies from the TNG100 simulation (stellar masses log (M*/M⊙) ~ 9.7-11.2). Comparison of observed magnesium and europium for individual stars in the Milky Way with the stellar abundances in our more than 850 Milky Way-like galaxies provides stringent constraints on our chemical evolutionary methods. Here, we use the magnesium-to-iron ratio as a proxy for the effects of our SNII (core-collapse supernovae) and SNIa (Type Ia supernovae) metal return prescription and as a comparison to a variety of galactic observations. The europium-to-iron ratio tracks the rare ejecta from neutron star-neutron star mergers, the assumed primary site of europium production in our models, and is a sensitive probe of the effects of metal diffusion within the gas in our simulations. We find that europium abundances in MilkyWay-like galaxies show no correlation with assembly history, present-day galactic properties, and average galactic stellar population age. We reproduce the europium-to-iron spread at low metallicities observed in the Milky Way, and find it is sensitive to gas properties during redshifts z ≈ 2-4. We show that while the overall normalization of [Eu/Fe] is susceptible to resolution and post-processing assumptions, the relatively large spread of [Eu/Fe] at low [Fe/H] when compared to that at high [Fe/H] is quite robust.

© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society. Using the IllustrisTNG simulations, we investigate the connection between galaxy morphology and star formation in central galaxies with stellar masses in the range 109-1011.5 M⊙. We quantify galaxy morphology by a kinematical decomposition of the stellar component into a spheroidal and a disc component (spheroid-to-total ratio, S/T) and by the concentration of the stellar mass density profile (C82). S/T is correlated with stellar mass and star formation activity, while C82 correlates only with stellar mass. Overall, we find good agreement with observational estimates for both S/T and C82. Low- and high-mass galaxies are dominated by random stellar motion, while only intermediate-mass galaxies (M∗ ≈ 1010-1010.5 M⊙) are dominated by ordered rotation. Whereas higher mass galaxies are typical spheroids with high concentrations, lower mass galaxies have low concentration, pointing to different formation channels. Although we find a correlation between S/T and star formation activity, in the TNG model galaxies do not necessarily change their morphology when they transition through the green valley or when they cease their star formation, this depending on galaxy stellar mass and morphological estimator. Instead, the morphology (S/T and C82) is generally set during the star-forming phase of galaxies. The apparent correlation between S/T and star formation arises because earlier forming galaxies had, on average, a higher S/T at a given stellar mass. Furthermore, we show that mergers drive in situ bulge formation in intermediate-mass galaxies and are responsible for the recent spheroidal mass assembly in the massive galaxies with M∗ > 1011 M⊙. In particular, these massive galaxies assemble about half of the spheroidal mass while star-forming and the other half through mergers while quiescent.

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (Mstar), and black hole mass (MBH) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and 10 model variations that alter the parameters of the black hole model. Focusing on central galaxies at z = 0 with Mstar > 1010 M⊙, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular MBH threshold above which galaxies are found to sharply transition from being mostly star forming to mostly quiescent. As a result of this behaviour, the fraction of quiescent galaxies as a function of Mstar is sensitive to both the normalization of the MBH–Mstar relation and the MBH threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical MBH measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in MBH at a given Mstar and show a smoother decline in sSFR with MBH.

© 2019 The Author(s) 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* ≳ 107 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.

© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society Motivated by the recently discovered kinematic 'Hubble sequence' shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range of M∗/M☉ ∈ [109.7, 1011.4] selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties. We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependences. In particular, the mean mass dependences of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at M∗ ≈ 1-2 × 1010 M☉ in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) M∗ ≈ 6 × 1010 M☉ are systematically higher than the observational data by ≾ 10 per cent (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed.

© 2016 The Authors. The inefficiency of star formation in massive elliptical galaxies is widely believed to be caused by the interactions of an active galactic nucleus (AGN) with the surrounding gas. Achieving a sufficiently rapid reddening of moderately massive galaxies without expelling too many baryons has however proven difficult for hydrodynamical simulations of galaxy formation, prompting us to explore a new model for the accretion and feedback effects of supermassive black holes. For high-accretion rates relative to the Eddington limit, we assume that a fraction of the accreted rest mass energy heats the surrounding gas thermally, similar to the 'quasar mode' in previouswork. For low-accretion rates, we invoke a new, pure kinetic feedback model that imparts momentum to the surrounding gas in a stochastic manner. These two modes of feedback are motivated both by theoretical conjectures for the existence of different types of accretion flows as well as recent observational evidence for the importance of kinetic AGN winds in quenching galaxies.We find that a large fraction of the injected kinetic energy in this mode thermalizes via shocks in the surrounding gas, thereby providing a distributed heating channel. In cosmological simulations, the resulting model produces red, non-star-forming massive elliptical galaxies, and achieves realistic gas fractions, black hole growth histories and thermodynamic profiles in large haloes.

© 2020 The Author(s). We use the IllustrisTNG hydrodynamical simulations to show how the fractions of quenched galaxies vary across different environments and cosmic time, and to quantify the role AGN feedback and pre-processing play in quenching group and cluster satellites. At z = 0, we select galaxies with Mstars = 109-12 M⊙ residing within (≤R200c) massive groups and clusters of total host mass M200c = 1013-15.2 M⊙ in TNG100 and TNG300. The model predicts a quenched fraction of ∼70-90 per cent (on average) for centrals and satellites of mass ≳1010.5 M⊙, regardless of host mass, cosmic time (0 ≤ z ≤ 0.5), cluster-centric distance, and time since infall in the z = 0 host. Low-mass central galaxies (≲ 1010 M⊙), on the other hand, are rarely quenched unless they become members of groups (1013-14 M⊙) or clusters (≥1014 M⊙), where the quenched fraction rises to ∼ 80 per cent. Typically, the fraction of low-mass passive galaxies is higher closer to the host centre and for progressively more massive hosts. The population of low-mass satellites accreted more than ∼ 4-6 Gyr ago in massive hosts is almost entirely passive, thus suggesting an upper limit for the time needed for environmental quenching to occur. In fact, about 30 per cent of group and cluster satellites that are quenched at z = 0 were already quenched before falling into their current host, and the bulk of them quenched as early as 4-10 billion years ago. For low-mass galaxies (Mstars ≲1010-10.5 M⊙), this is due to pre-processing, whereby current satellites may have been members of other hosts, and hence have undergone environmental processes, before falling into their final host, this mechanism being more common and more effective for the purposes of quenching for satellites found today in more massive hosts. On the other hand, massive galaxies quench on their own and because of AGN feedback, regardless of whether they are centrals or satellites.

© 2017 The Authors. We introduce an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO. The overall framework builds upon the successes of the Illustris galaxy formation model, and includes prescriptions for star formation, stellar evolution, chemical enrichment, primordial and metal-line cooling of the gas, stellar feedback with galactic outflows, and black hole formation, growth and multimode feedback. In this paper, we give a comprehensive description of the physical and numerical advances that form the core of the IllustrisTNG (The Next Generation) framework.We focus on the revised implementation of the galactic winds, of which we modify the directionality, velocity, thermal content and energy scalings, and explore its effects on the galaxy population. As described in earlier works, the model also includes a new black-hole-driven kinetic feedback at low accretion rates, magnetohydrodynamics and improvements to the numerical scheme. Using a suite of (25Mpc h-1)3 cosmological boxes, we assess the outcome of the new model at our fiducial resolution. The presence of a selfconsistently amplified magnetic field is shown to have an important impact on the stellar content of 1012M⊙ haloes and above. Finally, we demonstrate that the new galactic winds promise to solve key problems identified in Illustris in matching observational constraints and affecting the stellar content and sizes of the low-mass end of the galaxy population.

© 2019 The Author(s). We select galaxies from the IllustrisTNG hydrodynamical simulations (Mstars109Mo at 0 ≤ z ≤ 2) and characterize the shapes and evolutions of their UVJ and star formation rate Cstellar mass (SFR-CMstars) diagrams. We quantify the systematic uncertainties related to different criteria to classify star-forming versus quiescent galaxies, different SFR estimates, and by accounting for the star formation measured within different physical apertures. The TNG model returns the observed features of the UVJ diagram at z ≤ 2, with a clear separation between two classes of galaxies. It also returns a tight star-forming main sequence (MS) for Mstars < 1010.5 (Mo) with a ∼0.3 dex scatter at z ∼ 0 in our fiducial choices. If a UVJbased cut is adopted, the TNG MS exhibits a downwardly bending at stellar masses of about 1010.5-10.7 Mo. Moreover, the model predicts that ∼80 (50) per cent of 1010.5-11Mo galaxies at z = 0 (z = 2) are quiescent and the numbers of quenched galaxies at intermediate redshifts and high masses are in better agreement with observational estimates than previous models. However, shorter SFR-averaging time-scales imply higher normalizations and scatter of the MS, while smaller apertures lead to underestimating the galaxy SFRs: Overall we estimate the inspected systematic uncertainties to sum up to about 0.2.0.3 dex in the locus of the MS and to about 15 percentage points in the fraction of quenched galaxies. While TNG colour distributions are clearly bimodal, this is not the case for the SFR logarithmic distributions in bins of stellar mass (SFR ≥ 10.3 Moy-1). Finally, the slope and z = 0 normalization of the TNG MS are consistent with observational findings; however, the locus of the TNG MS remains lower by about 0.2.0.5 dex at 0.75 ≤z <2 than the available observational estimates taken at face value.

© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society. Using the IllustrisTNG simulations, we investigate the connection between galaxy morphology and star formation in central galaxies with stellar masses in the range 109-1011.5 M⊙. We quantify galaxy morphology by a kinematical decomposition of the stellar component into a spheroidal and a disc component (spheroid-to-total ratio, S/T) and by the concentration of the stellar mass density profile (C82). S/T is correlated with stellar mass and star formation activity, while C82 correlates only with stellar mass. Overall, we find good agreement with observational estimates for both S/T and C82. Low- and high-mass galaxies are dominated by random stellar motion, while only intermediate-mass galaxies (M∗ ≈ 1010-1010.5 M⊙) are dominated by ordered rotation. Whereas higher mass galaxies are typical spheroids with high concentrations, lower mass galaxies have low concentration, pointing to different formation channels. Although we find a correlation between S/T and star formation activity, in the TNG model galaxies do not necessarily change their morphology when they transition through the green valley or when they cease their star formation, this depending on galaxy stellar mass and morphological estimator. Instead, the morphology (S/T and C82) is generally set during the star-forming phase of galaxies. The apparent correlation between S/T and star formation arises because earlier forming galaxies had, on average, a higher S/T at a given stellar mass. Furthermore, we show that mergers drive in situ bulge formation in intermediate-mass galaxies and are responsible for the recent spheroidal mass assembly in the massive galaxies with M∗ > 1011 M⊙. In particular, these massive galaxies assemble about half of the spheroidal mass while star-forming and the other half through mergers while quiescent.

© 2018 The Author(s). Published by Oxford University Press on behalf of the Royal Astronomical Society. The distribution of elements in galaxies provides a wealth of information about their production sites and their subsequent mixing into the interstellar medium. Here we investigate the elemental distributions of stars in the IllustrisTNG simulations. We analyse the abundance ratios of magnesium and europium in Milky Way-like galaxies from the TNG100 simulation (stellar masses log (M*/M⊙) ~ 9.7-11.2). Comparison of observed magnesium and europium for individual stars in the Milky Way with the stellar abundances in our more than 850 Milky Way-like galaxies provides stringent constraints on our chemical evolutionary methods. Here, we use the magnesium-to-iron ratio as a proxy for the effects of our SNII (core-collapse supernovae) and SNIa (Type Ia supernovae) metal return prescription and as a comparison to a variety of galactic observations. The europium-to-iron ratio tracks the rare ejecta from neutron star-neutron star mergers, the assumed primary site of europium production in our models, and is a sensitive probe of the effects of metal diffusion within the gas in our simulations. We find that europium abundances in MilkyWay-like galaxies show no correlation with assembly history, present-day galactic properties, and average galactic stellar population age. We reproduce the europium-to-iron spread at low metallicities observed in the Milky Way, and find it is sensitive to gas properties during redshifts z ≈ 2-4. We show that while the overall normalization of [Eu/Fe] is susceptible to resolution and post-processing assumptions, the relatively large spread of [Eu/Fe] at low [Fe/H] when compared to that at high [Fe/H] is quite robust.

© 2019 The Author(s) 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* ≳ 107 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.

© 2017 The Author(s). The distribution of metals in the intra-cluster medium (ICM) encodes important information about the enrichment history and formation of galaxy clusters. Here, we explore the metal content of clusters in IllustrisTNG - a new suite of galaxy formation simulations building on the Illustris project. Our cluster sample contains 20 objects in TNG100-a ~(100 Mpc)3 volume simulation with 2 × 18203 resolution elements, and 370 objects in TNG300-a ~(300 Mpc)3 volume simulation with 2 × 25003 resolution elements. The z = 0 metallicity profiles agree with observations, and the enrichment history is consistent with observational data going beyond z~1, showing nearly no metallicity evolution. The abundance profiles vary only minimally within the cluster samples, especially in the outskirts with a relative scatter of ~15 per cent. The average metallicity profile flattens towards the centre, where we find a logarithmic slope of -0.1 compared to -0.5 in the outskirts. Cool core clusters have more centrally peaked metallicity profiles (~0.8 solar) compared to non-cool core systems (~0.5 solar), similar to observational trends. Si/Fe and O/Fe radial profiles follow positive gradients. The outer abundance profiles do not evolve below z ~ 2, whereas the inner profiles flatten towards z = 0. More than ~80 per cent of the metals in the ICM have been accreted from the proto-cluster environment, which has been enriched to ~0.1 solar already at z ~ 2. We conclude that the intra-cluster metal distribution is uniform among our cluster sample, nearly time-invariant in the outskirts for more than 10 Gyr, and forms through a universal enrichment history.

© 2018 The Author(s). We use the IllustrisTNG simulations to investigate the evolution of the mass-metallicity relation (MZR) for star-forming cluster galaxies as a function of the formation history of their cluster host. The simulations predict an enhancement in the gas-phase metallicities of star-forming cluster galaxies (109 < M* < 1010M⊙ h-1) at z ≤ 1.0 in comparisons to field galaxies. This is qualitatively consistent with observations. We find that the metallicity enhancement of cluster galaxies appears prior to their infall into the central cluster potential, indicating for the first time a systematic 'chemical pre-processing' signature for infalling cluster galaxies. Namely, galaxies that will fall into a cluster by z = 0 show a ~0.05 dex enhancement in the MZR compared to field galaxies at z ≤ 0.5. Based on the inflow rate of gas into cluster galaxies and its metallicity, we identify that the accretion of pre-enriched gas is the key driver of the chemical evolution of such galaxies, particularly in the stellar mass range (109 < M* < 1010M⊙ h-1). We see signatures of an environmental dependence of the ambient/inflowing gas metallicity that extends well outside the nominal virial radius of clusters. Our results motivate future observations looking for pre-enrichment signatures in dense environments.

The fundamental metallicity relation (FMR) is a postulated correlation between galaxy stellar mass, star formation rate (SFR), and gas-phase metallicity. At its core, this relation posits that offsets from the mass-metallicity relation (MZR) at a fixed stellar mass are correlated with galactic SFR. In this Letter, we use hydrodynamical simulations to quantify the time-scales over which populations of galaxies oscillate about the average SFR and metallicity values at fixed stellarmass.We find that Illustris and IllustrisTNG predict that galaxy offsets from the star formation main sequence and MZR oscillate over similar time-scales, are often anticorrelated in their evolution, evolve with the halo dynamical time, and produce a pronounced FMR. Our models indicate that galaxies oscillate about equilibrium SFR and metallicity values - set by the galaxy's stellar mass - and that SFR and metallicity offsets evolve in an anticorrelated fashion. This anticorrelated variability of the metallicity and SFR offsets drives the existence of the FMR in our models. In contrast to Illustris and IllustrisTNG, we speculate that the SFR and metallicity evolution tracks may become decoupled in galaxy formation models dominated by feedback-driven globally bursty SFR histories, which could weaken the FMR residual correlation strength. This opens the possibility of discriminating between bursty and non-bursty feedback models based on the strength and persistence of the FMR - especially at high redshift., United States. National Aeronautics and Space Administration (HST-HF2-51341.001-A), United States. Department of Energy. Computational Science Graduate Fellowship Program (Grant DEFG02-97ER25308)

© 2017 The Author(s). Published by Oxford University Press on behalf of the Royal Astronomical Society. We analyse scaling relations and evolution histories of galaxy sizes in TNG100, part of the IllustrisTNG simulation suite. Observational qualitative trends of size with stellar mass, star formation rate and redshift are reproduced, and a quantitative comparison of projected r band sizes at 0 ≲ z ≲ 2 shows agreement to much better than 0.25 dex. We follow populations of z = 0 galaxies with a range of masses backwards in time along their main progenitor branches, distinguishing between main-sequence and quenched galaxies. Our main findings are as follows. (i) At M*, z = 0 ≳ 109.5 M⊙, the evolution of the median main progenitor differs, with quenched galaxies hardly growing in median size before quenching, whereas main-sequence galaxies grow their median size continuously, thus opening a gap from the progenitors of quenched galaxies. This is partly because the main-sequence high-redshift progenitors of quenched z = 0 galaxies are drawn from the lower end of the size distribution of the overall population of main-sequence high-redshift galaxies. (ii) Quenched galaxies with M*, M*, z = 0 ≳ 109.5 M⊙ experience a steep size growth on the size-mass plane after their quenching time, but with the exception of galaxies withM*, z = 0 ≳ 1011 M⊙, the size growth after quenching is small in absolute terms, such that most of the size (and mass) growth of quenched galaxies (and its variation among them) occurs while they are still on the main sequence. After they become quenched, the size growth rate of quenched galaxies as a function of time, as opposed to versus mass, is similar to that of main-sequence galaxies. Hence, the size gap is retained down to z = 0.

© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society Motivated by the recently discovered kinematic 'Hubble sequence' shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range of M∗/M☉ ∈ [109.7, 1011.4] selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties. We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependences. In particular, the mean mass dependences of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at M∗ ≈ 1-2 × 1010 M☉ in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) M∗ ≈ 6 × 1010 M☉ are systematically higher than the observational data by ≾ 10 per cent (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed.

© 2018 The Author(s) Published by Oxford University Press on behalf of The Royal Astronomical Society. We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black holes cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about 108.5M⊙, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass - stellar bulge mass relation is in good agreement with observations, whereas the simulation tends to overpredict the high-redshift quasar luminosity function.

ABSTRACT Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (Mstar), and black hole mass (MBH) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and 10 model variations that alter the parameters of the black hole model. Focusing on central galaxies at z = 0 with Mstar > 1010 M⊙, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular MBH threshold above which galaxies are found to sharply transition from being mostly star forming to mostly quiescent. As a result of this behaviour, the fraction of quiescent galaxies as a function of Mstar is sensitive to both the normalization of the MBH–Mstar relation and the MBH threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical MBH measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in MBH at a given Mstar and show a smoother decline in sSFR with MBH.

© 2018 The Author(s). The galaxy Malin 1 contains the largest stellar disc known but the formation mechanism of this structure has been elusive. In this paper, we report a Malin 1 analogue in the 100 Mpc IllustrisTNG simulation and describe its formation history. At redshift zero, this massive galaxy, having a maximum circular velocity Vmax of 430 km s-1, contains a 100 kpc gas/stellar disc with morphology similar toMalin 1. The simulated galaxy reproduces well many observed features of Malin 1's vast disc, including its stellar ages, metallicities, and gas rotation curve. We trace the extended disc back in time and find that a large fraction of the cold gas at redshift zero originated from the cooling of hot halo gas, triggered by the merger of a pair of intruding galaxies. Our finding provides a novel way to form large galaxy discs as extreme as Malin 1 within the current galaxy formation framework.

The coevolution of galaxies and their metal content serves as an important test for galaxy feedback models. We analyze the distribution and evolution of metals within the IllustrisTNG simulation suite with a focus on the gas-phase mass-metallicity relation (MZR). We find that the IllustrisTNG model broadly reproduces the slope and normalization evolution of the MZR across the redshift range $02$) evolution fit by $\mathrm{d\;log(Z)}/\mathrm{dz} \approx - 0.064$. Our simulations indicate that the metal retention efficiency of the interstellar medium (ISM) is low: a majority of gas-phase metals ($\sim$ 85 per cent at $z=0$) live outside of the ISM, either in an extended gas disk, the circumgalactic medium, or outside the halo. Nevertheless, the redshift evolution in the simulated MZR normalization is driven by the higher gas fractions of high redshift galaxies, not by changes to the metal retention efficiency. The scatter in the simulated MZR contains a clear correlation with the gas-mass or star formation rate of the system, in agreement with the observed fundamental metallicity relation. The scatter in the MZR is driven by a competition between periods of enrichment- and accretion-dominated metallicity evolution. We expect that while the normalization of the MZR declines with redshift, the strength of the correlation between metallicity and gas-mass at fixed stellar mass is not a strong function of redshift. Our results indicate that the "regulator" style models are best suited for simultaneously explaining the shape, redshift evolution, and existence of correlated scatter with gas fraction about the MZR.

© 2016 The Authors. The inefficiency of star formation in massive elliptical galaxies is widely believed to be caused by the interactions of an active galactic nucleus (AGN) with the surrounding gas. Achieving a sufficiently rapid reddening of moderately massive galaxies without expelling too many baryons has however proven difficult for hydrodynamical simulations of galaxy formation, prompting us to explore a new model for the accretion and feedback effects of supermassive black holes. For high-accretion rates relative to the Eddington limit, we assume that a fraction of the accreted rest mass energy heats the surrounding gas thermally, similar to the 'quasar mode' in previouswork. For low-accretion rates, we invoke a new, pure kinetic feedback model that imparts momentum to the surrounding gas in a stochastic manner. These two modes of feedback are motivated both by theoretical conjectures for the existence of different types of accretion flows as well as recent observational evidence for the importance of kinetic AGN winds in quenching galaxies.We find that a large fraction of the injected kinetic energy in this mode thermalizes via shocks in the surrounding gas, thereby providing a distributed heating channel. In cosmological simulations, the resulting model produces red, non-star-forming massive elliptical galaxies, and achieves realistic gas fractions, black hole growth histories and thermodynamic profiles in large haloes.

© 2018 The Author(s). We have generated synthetic images of ∼27 000 galaxies from the IllustrisTNG and the original Illustris hydrodynamic cosmological simulations, designed to match Pan-STARRS observations of log10(M∗/M·)-∼9.8-11.3 galaxies at z ∼0.05. Most of our synthetic images were created with the SKIRT radiative transfer code, including the effects of dust attenuation and scattering, and performing the radiative transfer directly on the Voronoi mesh used by the simulations themselves. We have analysed both our synthetic and real Pan-STARRS images with the newly developed statmorph code, which calculates non-parametric morphological diagnostics-including the Gini-M20 and concentration-asymmetry-smoothness statistics-and performs 2D Sersic fits. Overall, we find that the optical morphologies of IllustrisTNG galaxies are in good agreement with observations, and represent a substantial improvement compared to the original Illustris simulation. In particular, the locus of the Gini-M20 diagram is consistent with that inferred from observations, while the median trends with stellar mass of all the morphological, size and shape parameters considered in this work lie within the ∼1s scatter of the observational trends. However, the IllustrisTNG model has some difficulty with more stringent tests, such as producing a strong morphology-colour relation. This results in a somewhat higher fraction of red discs and blue spheroids compared to observations. Similarly, the morphology-size relation is problematic: while observations show that discs tend to be larger than spheroids at a fixed stellar mass, such a trend is not present in IllustrisTNG.

© 2018 The Author(s). Published by Oxford University Press on behalf of the Royal Astronomical Society. We explore the abundance, spatial distribution, and physical properties of the OVI, OVII, and OVIII ions of oxygen in circumgalactic and intergalactic media (the CGM, IGM, and WHIM). We use the TNG100 and TNG300 large volume cosmological magnetohydrodynamical simulations. Modelling the ionization states of simulated oxygen, we find good agreement with observations of the low-redshiftOVI column density distribution function (CDDF), and present its evolution for all three ions from z = 0 to z = 4. Producing mock quasar absorption line spectral surveys, we show that the IllustrisTNG simulations are fully consistent with constraints on the OVI content of the CGM from COS-haloes and other low-redshift observations, producing columns as high as observed. We measure the total amount of mass and average column densities of each ion using hundreds of thousands of simulated galaxies spanning 1011 < Mhalo/M⊙ < 1015 corresponding to 109 < M/M⊙ < 1012 in stellar mass. Stacked radial profiles of OVI are computed in 3D number density and 2D projected column density, decomposing into 1-halo and 2-halo terms. Relating halo OVI to properties of the central galaxy, we find a correlation between the (g - r) colour of a galaxy and the total amount of OVI in its CGM. In comparison to the COS-Haloes finding, this leads to a dichotomy of columns around star-forming versus passive galaxies at fixed stellar (or halo) mass. We demonstrate that this correlation is a direct result of black hole feedback associated with quenching and represents a causal consequence of galactic-scale baryonic feedback impacting the physical state of the circumgalactic medium.

We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good agreement with observations, but we note an excess of faint AGN in hard X-ray (L[subscript x] ∼ 10[superscript 43-44] rm erg/s), and a lower number of bright AGN (L[subscript x] > 110[superscript 44] rm erg/s), a conclusion that varies quantitatively but not qualitatively with BH luminosity estimation method. The lower Eddington ratios of the 10[superscript 9] M[subscript ⊙] BHs compared to observations suggest that AGN feedback may be too efficient in this regime. We study galaxy star formation activity and structural properties, and design sample-dependent criteria to identify different galaxy types (star-forming/quiescent, extended/compact) that we apply both to the simulations and observations from the candels fields. We analyse how the simulated and observed galaxies populate the specific star formation rate - stellar mass surface density diagram. A large fraction of the z = 0 M{subscript *] > 10[superscript 11] M[subscript ⊙] quiescent galaxies first experienced a compaction phase (i.e. reduction of galaxy size) while still forming stars, and then a quenching event. We measure the dependence of AGN fraction on galaxies' locations in this diagram. After correcting the simulations with a redshift and AGN luminosity-dependent model for AGN obscuration, we find good qualitative and quantitative agreement with observations. The AGN fraction is the highest among compact star-forming galaxies (16-20% at z ∼1.5-2), and the lowest among compact quiescent galaxies (6-10 rm per cent at z ∼ 1.5-2). Key words: galaxies: formation - galaxies: evolution - methods: numerical

We use the IllustrisTNG simulations to study the demographics and properties of jellyfish galaxies in the full cosmological context. By jellyfish galaxies, we mean satellites orbiting in massive groups and clusters that exhibit highly asymmetric distributions of gas and gas tails. In particular, we select TNG100 galaxies at low redshifts (⁠z ≤ 0.6) with stellar mass exceeding 10[superscript 9.5] M[subscript ⊙] and with host halo masses in the range 10[superscript 13] ≤ M[subscript 200c]/M[subscript ⊙] ≤ 1014.6⁠. Among more than about 6000 (2600) galaxies with stars (and some gas), we identify 800 jellyfish galaxies by visually inspecting their gas and stellar mass maps in random projections. Namely, about 31 per cent of cluster satellites are found with signatures of ram-pressure stripping and gaseous tails stemming from their main luminous bodies. This is a lower limit: the random orientation entails a loss of about 30 per cent of galaxies that in an optimal projection would otherwise be identified as jellyfish. Furthermore, jellyfish galaxies are more frequent at intermediate and large cluster-centric distances (r/R[subscript 200c] ≳ 0.25), in more massive hosts and at smaller satellite masses, and they typically orbit supersonically. The gaseous tails usually extend in opposite directions to the galaxy trajectory, with no relation between tail orientation and position of the host’s centre. Finally, jellyfish galaxies are late infallers (<2.5–3 Gyr ago, at z = 0) and the emergence of gaseous tails correlates well with the presence of bow shocks in the intracluster medium. Keywords: methods:numerical; galaxies:clusters:general; galaxies:clusters:intracluster medium; galaxies:evolution; galaxies:groups:general; cosmology:theory