Search

We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical 'zoom-in' simulations of 40 Milky Way-mass haloes and 26 dwarf galaxy–mass haloes run with the moving-mesh code arepo. Auriga adopts the Lambda cold dark matter cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully–Fisher relation, the star-forming (SF) main sequence, and H  i gas fraction/disc thickness. Finally, we show that SF gas discs appear to build rotation and velocity dispersion rapidly for |$z\gtrsim 3$| before they 'settle' into ever-increasing rotation-dispersion ratios (⁠|$V/\sigma$|⁠). This evolution appears to be in rough agreement with some kinematic measurements from H |$\alpha$| observations, and demonstrates an application of how to utilize the released data. [ABSTRACT FROM AUTHOR]

In this study, we investigate the impact of local stellar radiation in cosmological zoom simulations of the formation of Milky Way-sized galaxies. We include the radiation field as an additional feedback component that is computed alongside gravity with a tree code in an optically thin approximation. We resimulate the initial conditions of five Milk Way-like systems taken from the Auriga project with and without stellar radiation, and study the effects of local stellar radiation on several properties of the galaxies and the circumgalactic medium (CGM). Similar to previous findings, we observe with our current model that local stellar radiation can modify gas cooling in the CGM and thus suppress star formation and the surface densities of young stars and H  i gas, while having little impact on the total gas content. In particular, it also suppresses the peak of the rotation curve and reduces the mass of the stellar bulge. In the CGM region, the young stellar radiation exceeds the external ultraviolet background and dominates the radiation field within the virial halo at all redshifts. Nevertheless, we find that the local stellar radiation, as implemented in the current study, has overall little impact on the radial density and temperature profile of the CGM gas. However, for the ion species H  i and Mg  ii , the column densities within |$\sim 0.3\, R_{\rm vir}$| are reduced, while the O  vi column density is hardly impacted by the radiation field due to a lack of soft X-ray components in our current model. Additional effects can be expected from the radiation of the central active galactic nucleus during phases of quasar activity and from soft X-ray sources, which have not yet been included in the simulations of the present study. [ABSTRACT FROM AUTHOR]

The properties of galaxies in low-density regions of the universe suggest an interplay between galaxy formation and environment. However, the specific reason why this particular large-scale environment influences the evolution of galaxies remains unclear. This paper examines the properties and evolutionary paths of galaxies within cosmic voids using the Illustris TNG300 simulation. The population of void galaxies at z = 0 has a higher star formation rate, a smaller stellar-to-halo-mass ratio, higher gas metallicity, and lower stellar metallicity in comparison with non-void galaxies at fixed stellar mass. Our analysis shows that these differences are mainly due to the characteristics of galaxies classified as satellites, for which the largest differences between void and non-void samples are found. Although the mean number of mergers is similar between void and non-void samples at a fixed stellar mass, void galaxies tend to experience mergers at later times, resulting in a more recent accumulation of accreted stellar mass. While the mean net accreted mass is comparable for high-mass galaxies, low-mass void galaxies tend to exhibit higher fractions of accreted stars than non-void galaxies. This finding challenges the common notion that void galaxies predominantly experience growth with infrequent mergers or interactions. With this paper, we also publicly release our void catalogue as part of the IllustrisTNG project. [ABSTRACT FROM AUTHOR]

We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar distributions of around 70 per cent of the simulated galaxies exhibit an 'inside–out' pattern, with older (younger) stellar populations preferentially located in the inner (outer) disc regions. In all cases, we find a very tight correlation between the infall, outflow, and net accretion rates, as well as between these three quantities and the star formation rate. This is because the amount of gas which is ultimately available for star formation in each radial ring depends not only on the infall rates, but also on the amount of gas leaving the disc in outflows, which directly relates to the local star formation level. Therefore, any of these rates can be used to identify galaxies with inside–out growth. For these galaxies, the correlation between the dominant times of accretion/star formation and disc radius is well fitted by a linear function. We also find that, when averaged over galaxies with formation histories similar to the Milky Way, the simulated accretion rates show a similar evolution (both temporally and radially integrated) to the usual accretion prescriptions used in chemical evolution models, although some major differences arise at early times and in the inner disc regions. [ABSTRACT FROM AUTHOR]

Cosmic rays (CRs) are dynamically important in the evolution of galaxies by regulating star formation and powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magnetohydrodynamical simulations of entire galaxies with masses ranging from 1010 to |$10^{12}\, \mathrm{M}_\odot$| and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from |$0.1\, \mathrm{GeV} c^{-1}$| to |$100\, \mathrm{TeV} c^{-1}$|⁠. We find that hadronic cooling of CRs dominates over Alfvén cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of 4 in dwarf galaxies. All simulated low-mass haloes (M  = 1010, 1011, and |$3\times 10^{11}\, \mathrm{M}_\odot$|⁠) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with |$M=10^{12}\, \mathrm{M}_\odot$| does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of |$\langle {D}\rangle _{e_\mathrm{cr}}\sim 10^{28}\, \mathrm{cm^2\, s^{-1}}$| in the disc up to |$3\times 10^{29}\, \mathrm{cm^2\, s^{-1}}$| in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models. [ABSTRACT FROM AUTHOR]

Upcoming large galaxy surveys will subject the standard cosmological model, Lambda Cold Dark Matter, to new precision tests. These can be tightened considerably if theoretical models of galaxy formation are available that can predict galaxy clustering and galaxy–galaxy lensing on the full range of measurable scales, throughout volumes as large as those of the surveys, and with sufficient flexibility that uncertain aspects of the underlying astrophysics can be marginalized over. This, in particular, requires mock galaxy catalogues in large cosmological volumes that can be directly compared to observation, and can be optimized empirically by Monte Carlo Markov Chains or other similar schemes, thus eliminating or estimating parameters related to galaxy formation when constraining cosmology. Semi-analytic galaxy formation methods implemented on top of cosmological dark matter simulations offer a computationally efficient approach to construct physically based and flexibly parametrized galaxy formation models, and as such they are more potent than still faster, but purely empirical models. Here, we introduce an updated methodology for the semi-analytic L-Galaxies code, allowing it to be applied to simulations of the new MillenniumTNG project, producing galaxies directly on fully continuous past lightcones, potentially over the full sky, out to high redshift, and for all galaxies more massive than |$\sim 10^8\, {\rm M}_\odot$|⁠. We investigate the numerical convergence of the resulting predictions, and study the projected galaxy clustering signals of different samples. The new methodology can be viewed as an important step towards more faithful forward-modelling of observational data, helping to reduce systematic distortions in the comparison of theory to observations. [ABSTRACT FROM AUTHOR]

The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of m b ≈ 104 M⊙), to additionally enable the study of 'mini-haloes' with virial temperatures T vir < 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M 1500 =   − 10) by z  = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties. [ABSTRACT FROM AUTHOR]

The early release science results from JWST have yielded an unexpected abundance of high-redshift luminous galaxies that seems to be in tension with current theories of galaxy formation. However, it is currently difficult to draw definitive conclusions form these results as the sources have not yet been spectroscopically confirmed. It is in any case important to establish baseline predictions from current state-of-the-art galaxy formation models that can be compared and contrasted with these new measurements. In this work, we use the new large-volume (⁠|$L_\mathrm{box}\sim 740 \, \mathrm{cMpc}$|⁠) hydrodynamic simulation of the MillenniumTNG project, suitably scaled to match results from higher resolution – smaller volume simulations, to make predictions for the high-redshift (z ≳ 8) galaxy population and compare them to recent JWST observations. We show that the simulated galaxy population is broadly consistent with observations until z ∼ 10. From z ≈ 10–12, the observations indicate a preference for a galaxy population that is largely dust-free, but is still consistent with the simulations. Beyond z ≳ 12, however, our simulation results underpredict the abundance of luminous galaxies and their star-formation rates by almost an order of magnitude. This indicates either an incomplete understanding of the new JWST data or a need for more sophisticated galaxy formation models that account for additional physical processes such as Population III stars, variable stellar initial mass functions, or even deviations from the standard ΛCDM model. We emphasize that any new process invoked to explain this tension should only significantly influence the galaxy population beyond z ≳ 10, while leaving the successful galaxy formation predictions of the fiducial model intact below this redshift. [ABSTRACT FROM AUTHOR]

Momentum-dependent cosmic ray diffusion drives powerful galactic winds The largest difference between the non-cooling and the correct model is the fraction of low energy CRs to the total CR energy. In the paper we present simulations with a new spectral cosmic ray (CR) solver, that includes energy dependent spatial diffusion and spectral cooling as introduced in Girichidis et al. ([1]). [Extracted from the article]

Dwarf galaxies are thought to quench primarily due to environmental processes most typically occurring in galaxy groups and clusters or around single, massive galaxies. However, at earlier epochs, (5 < z < 2), the collapse of large-scale structure (forming Zel'dovich sheets and subsequently filaments of the cosmic web) can produce volume-filling accretion shocks which elevate large swaths of the intergalactic medium (IGM) in these structures to a hot (T > 106 K) phase. We study the impact of such an event on the evolution of central dwarf galaxies (5.5 < log  M * < 8.5) in the field using a spatially large, high resolution cosmological zoom simulation which covers the cosmic web environment between two protoclusters. We find that the shock-heated sheet acts as an environmental quencher much like clusters and filaments at lower redshift, creating a population of quenched, central dwarf galaxies. Even massive dwarfs that do not quench are affected by the shock, with reductions to their sSFR and gas accretion. This process can potentially explain the presence of isolated quenched dwarf galaxies, and represents an avenue of pre-processing, via which quenched satellites of bound systems quench before infall. [ABSTRACT FROM AUTHOR]

We study the presence and importance of stellar migration in the evolution of 17 Milky-Way like disc galaxies with stellar mass 10 < log(M */M⊙) < 11 from the Auriga suite of zoom-in cosmological hydrodynamical simulations. We compare the birth radii of the stars to their radii at z  = 0 for each system and present mean values of the strength of stellar migration as a function of radius and stellar age which vary between 1–4 kpc. We also investigate the effect of migration on age and metallicity radial profiles in the discs. We find several cases of age gradient flattening due to migration, but significant changes to metallicity profiles only for older stellar populations and discs that develop a strong bar. Furthermore, we study stellar migration from the perspective of the change of the galactocentric radius (Δ R) and orbital guiding centre radius (Δ Rg) of stellar particles between given time intervals. We find that stars migrate approximately as a diffusion process only in the outer parts of the discs and for particular galaxies that have a weak bar. Strongly barred galaxies in our sample show larger stellar migration but its timestep evolution is slower-than-diffusion. Finally, we give parametrizations that encapsulate the dependence of the strength of the radial migration as a function of time and radius, for incorporation into (semi-)analytic models of galaxy evolution. [ABSTRACT FROM AUTHOR]

Modern hydrodynamical simulations reproduce many properties of the real Universe. These simulations model various physical processes, but many of these are included using 'subgrid models' due to resolution limits. Although different subgrid models have been successful in modelling the effects of supernovae (SNe) and active galactic nuclei (AGNs) feedback on galactic properties, it remains unclear if, and by how much, these differing implementations affect observable halo gas properties. In this work, we use 'zoom-in' cosmological initial conditions of two volumes selected to resemble the Local Group (LG) evolved with both the Auriga and Apostle galaxy formation models. While the subgrid physics models in both simulations reproduce realistic stellar components of L ⋆ galaxies, they exhibit different gas properties. Namely, Auriga predicts that the Milky Way is almost baryonically closed, whereas Apostle suggests that only half of the expected baryons reside within the halo. Furthermore, Apostle predicts that this baryon deficiency extends to the LG (r ≤ 1 Mpc). Some of the baryon deficiency in Apostle is due to SNe feedback at high redshift, which generates halo-wide outflows, with high covering fractions and radial velocities, which both eject baryons and significantly impede cosmic gas accretion. Conversely, in Auriga , gas accretion into the halo appears to be almost unaffected by feedback. These differences appear to result from the different energy injection methods from feedback (both SNe and AGNs) to gas. Our results suggest that both quasar absorption lines and fast radio burst dispersion measures could constrain these two regimes with future observations. [ABSTRACT FROM AUTHOR]

Early-type galaxies (ETGs) possess total density profiles close to isothermal, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at 1.5 R eff (effective radius) that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on these kinematic features through mock ETGs in IllustrisTNG. We create mock Integral-Field-Unit observations to extract projected stellar kinematic features for 207 ETGs with stellar mass |$M_{\ast }\geqslant 10^{11} \, \mathrm{M_{\odot}}$| in TNG100-1. The mock observations reproduce the key outer (1.5 R eff) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope γouter and the kurtosis h 4's gradient. We find that h 4 is uncorrelated with stellar orbital anisotropy beyond R eff; instead, we find that the variations in γouter and outer h 4 (a good proxy for h 4 gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the origin of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates the two quantities. [ABSTRACT FROM AUTHOR]

Feedback from active galactic nuclei (AGNs) has become established as a fundamental process in the evolution of the most massive galaxies. Its impact on Milky Way (MW)-mass systems, however, remains comparatively unexplored. In this work, we use the auriga simulations to probe the impact of AGN feedback on the dynamical and structural properties of galaxies, focusing on the bar, bulge, and disc. We analyse three galaxies – two strongly and one unbarred/weakly barred – using three setups: (i) the fiducial auriga model, which includes both radio and quasar mode feedback, (ii) a setup with no radio mode, and (iii) one with neither the radio nor the quasar mode. When removing the radio mode, gas in the circumgalactic medium cools more efficiently and subsequently settles in an extended disc, with little effect on the inner disc. Contrary to previous studies, we find that although the removal of the quasar mode results in more massive central components, these are in the form of compact discs, rather than spheroidal bulges. Therefore, galaxies without quasar mode feedback are more baryon-dominated and thus prone to forming stronger and shorter bars, which reveals an anticorrelation between the ejective nature of AGN feedback and bar strength. Hence, we report that the effect of AGN feedback (i.e. ejective or preventive) can significantly alter the dynamical properties of MW-like galaxies. Therefore, the observed dynamical and structural properties of MW-mass galaxies can be used as additional constraints for calibrating the efficiency of AGN feedback models. [ABSTRACT FROM AUTHOR]

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as z  = 4, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scalelengths. The bar fraction mildly evolves with redshift for galaxies with |$M_{*}\ge 10^{10}\rm M_{\odot }$|⁠ , being greater than |$\sim 40{{\ \rm per\ cent}}$| at 0.5 < z < 3 and |$\sim 30{{\ \rm per\ cent}}$| at z  = 0. When bars larger than a given physical size (⁠|$\ge 2\, \rm kpc$|⁠) or the angular resolution limit of twice the I -band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions, and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at six redshifts between z  = 4 and 0. [ABSTRACT FROM AUTHOR]

We post-process galaxies in the IllustrisTNG simulations with skirt radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at z ≥ 4. The rest-frame K - and z -band galaxy luminosity functions from TNG are overall consistent with observations, despite |${\sim}0.5\, \mathrm{dex}$| underprediction at z = 4 for M K ≲ −25 and M z ≲ −24. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. Based on theoretical estimations, we show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of |$500\, {\rm arcmin}^{2}$| at z = 6 (z = 8). As opposed to the consistency in the UV, optical, and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshift are qualitatively consistent with observations, the peak dust temperature of z ≥ 6 galaxies are overestimated by about |$20\, {\rm K}$|⁠. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts. [ABSTRACT FROM AUTHOR]

We present the cumulative star formation histories (SFHs) of >15 000 dwarf galaxies (⁠|$M_{\rm *}=10^{7-10}\, {\rm M}_{\odot }$|⁠) simulated with the TNG50 run of the IllustrisTNG suite across a vast range of environments. The key factors that determine the dwarfs' SFHs are their central/satellite status and stellar mass, with centrals and more massive dwarfs assembling their stellar mass at later times, on average, compared to satellites and lower mass dwarfs. Satellites (in hosts of mass |$M_{\rm 200c, host}=10^{12-14.3}\, {\rm M}_{\odot }$|⁠) assembled 90 per cent of their stellar mass |${\sim}7.0_{-5.5}^{+3.3}$| Gyr ago, on average and within the 10th to 90th percentiles, while the centrals did so only |${\sim}1.0_{-0.5}^{+4.0}$| Gyr ago. TNG50 predicts a large diversity in SFHs, so that individual dwarfs can have significantly different cumulative SFHs compared to the stacked median SFHs. Satellite dwarfs with the highest stellar mass to host cluster mass ratios have the latest stellar mass assembly. Conversely, satellites at fixed stellar and host halo mass found closer to the cluster centre or accreted at earlier times show significantly earlier stellar mass assembly. These trends and the shapes of the SFHs themselves are a manifestation of the varying proportions within a given subsample of quenched versus star-forming galaxies, which exhibit markedly distinct SFH shapes. Finally, satellite dwarfs in the most massive hosts have higher SFRs at early times, well before accretion into their z = 0 host, compared to a control sample of centrals mass-matched at the time of accretion. This is the result of the satellites being preprocessed in smaller hosts prior to accretion. Our findings are useful theoretical predictions for comparison to future resolved stellar population observations. [ABSTRACT FROM AUTHOR]

Submillimetre galaxies (SMGs) have long posed a challenge for theorists, and self-consistently reproducing the properties of the SMG population in a large-volume cosmological hydrodynamical simulation has not yet been achieved. We use a scaling relation derived from previous simulations plus radiative transfer calculations to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects based on the simulated galaxies' star formation rates (SFRs) and dust masses, and compare the predicted number counts with observations. We find that the predicted SMG number counts based on IllustrisTNG are significantly less than observed (more than 1 dex at S 850 ≳ 4 mJy). The simulation from the original Illustris project yields more SMGs than IllustrisTNG : the predicted counts are consistent with those observed at both S 850 ≲ 5 mJy and S 850 ≳ 9 mJy and only a factor of ∼2 lower than those observed at intermediate flux densities. The redshift distribution of SMGs with S 850 > 3 mJy in IllustrisTNG is consistent with the observed distribution, whereas the Illustris redshift distribution peaks at significantly lower redshift (1.5 versus 2.8). We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass (⁠|$M_{\star }\sim 10^{11} \, \text{M}_{\odot }$|⁠) z ∼ 2–3 galaxies have lower dust masses and SFRs than in Illustris owing to differences in the subgrid models for stellar and/or active galactic nucleus feedback between the two simulations (we unfortunately cannot isolate the specific cause(s) post hoc). Our results demonstrate that because our method enables predicting SMG number counts in post-processing with a negligible computational expense, SMGs can provide useful constraints for tuning subgrid models in future large-volume cosmological simulations. [ABSTRACT FROM AUTHOR]

We compare predictions for galaxy–galaxy lensing profiles and clustering from the Henriques et al. public version of the Munich semi-analytical model (SAM) of galaxy formation and the IllustrisTNG suite, primarily TNG300, with observations from KiDS + GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, and stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with |$M_\ast \gt 10^{11}\, \mathrm{ M}_\odot$|⁠. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode active galactic nucleus accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high-mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only (χ2 = 1.81 compared to χ2 = 7.79 for the SAM) and locally brightest galaxy samples (χ2 = 3.80 compared to χ2 = 5.01). With added dust corrections to the colours it matches the SDSS clustering signal well for red low-mass galaxies. We find that both the SAMs and TNG300 predict |$\sim 50\, {{\ \rm per\ cent}}$| excessive lensing signals for intermediate-mass red galaxies with 10.2 < log10 M *[M⊙] < 11.2 at |$r \approx 0.6\, h^{-1}\, \text{Mpc}$|⁠ , which require further theoretical development. [ABSTRACT FROM AUTHOR]

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 t cool/ t ff < 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 m gas ≃ 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 FROM AUTHOR]

We present a new implementation for active galactic nucleus (AGN) feedback through small-scale, ultrafast winds in the moving-mesh hydrodynamic code arepo. The wind is injected by prescribing mass, momentum, and energy fluxes across a spherical boundary centred on a supermassive black hole according to available constraints for accretion disc winds. After sweeping-up a mass equal to their own, small-scale winds thermalize, powering energy-driven outflows with dynamics, structure, and cooling properties in excellent agreement with those of analytic wind solutions. Momentum-driven solutions do not easily occur, because the Compton cooling radius is usually much smaller than the free-expansion radius of the small-scale winds. Through various convergence tests, we demonstrate that our implementation yields wind solutions, which are well converged down to the typical resolution achieved in cosmological simulations. We test our model in hydrodynamic simulations of isolated Milky Way – mass galaxies. Above a critical AGN luminosity, initially spherical, small-scale winds power bipolar, energy-driven superwinds that break out of the galactic nucleus, flowing at speeds |$\gt 1000 \rm \, km \, s^{-1}$| out to |$\sim 10 \, \rm kpc$|⁠. These energy-driven outflows result in moderate, but long-term, reduction in star formation, which becomes more pronounced for higher AGN luminosities and faster small-scale winds. Suppression of star formation proceeds through a rapid mode that involves the removal of the highest density, nuclear gas, and through a slower mode that effectively halts halo gas accretion. Our new implementation makes it possible to model AGN-driven winds in a physically meaningful and validated way in simulations of galaxy evolution, the interstellar medium and black hole accretion flows. [ABSTRACT FROM AUTHOR]

We present predictions for high redshift (z  = 2−10) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the H α and H β +  |$[\rm O \,{\small III}]$| luminosity functions up to z  = 8. The predicted H β +  |$[\rm O \,{\small III}]$| luminosity functions are consistent with present observations at z ≲ 3 with |${\lesssim} 0.1\, {\rm dex}$| differences in luminosities. However, the predicted H α luminosity function is |${\sim }0.3\, {\rm dex}$| dimmer than the observed one at z ≃ 2. Furthermore, we explore continuum spectral indices, the Balmer break at 4000 Å; (D4000) and the UV continuum slope β. The median D4000 versus specific star formation rate relation predicted at z  = 2 is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed A UV versus β relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at z  = 2−6 and investigate their dependence on galaxy colours and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colours, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high-redshift galaxies. Future James Webb Space Telecope observations will further test these predictions. [ABSTRACT FROM AUTHOR]

The James Webb Space Telescope (JWST) promises to revolutionize our understanding of the early Universe, and contrasting its upcoming observations with predictions of the Λ cold dark matter model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multiband galaxy luminosity functions from z  = 2 to z  = 10. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalization and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalization compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust model calibrations. Furthermore, we also recover the observed high-redshift (z  = 4–6) UV luminosity versus stellar mass relation, the H α versus star formation rate relation, and the H α luminosity function at z  = 2. The bright end (M UV > −19.5) cumulative galaxy number densities are consistent with observational data. For the F 200 W NIRCam band, we predict that JWST will detect ∼80 (∼200) galaxies with a signal-to-noise ratio of 10 (5) within the NIRCam field of view, |$2.2\times 2.2 \, {\rm arcmin}^{2}$|⁠ , for a total exposure time of |$10^5\, {\rm s}$| in the redshift range z  = 8 ± 0.5. These numbers drop to ∼10 (∼40) for an exposure time of |$10^4\, {\rm s}$|⁠. [ABSTRACT FROM AUTHOR]

Using a set of 15 high-resolution magnetohydrodynamic cosmological simulations of Milky Way formation, we investigate the origin of the baryonic material found in stars at redshift zero. We find that roughly half of this material originates from subhalo/satellite systems and half is smoothly accreted from the intergalactic medium. About |$90 {{\ \rm per\ cent}}$| of all material has been ejected and re-accreted in galactic winds at least once. The vast majority of smoothly accreted gas enters into a galactic fountain that extends to a median galactocentric distance of ∼20 kpc with a median recycling time-scale of ∼500 Myr. We demonstrate that, in most cases, galactic fountains acquire angular momentum via mixing of low angular momentum, wind-recycled gas with high angular momentum gas in the circumgalactic medium (CGM). Prograde mergers boost this activity by helping to align the disc and CGM rotation axes, whereas retrograde mergers cause the fountain to lose angular momentum. Fountain flows that promote angular momentum growth are conducive to smooth evolution on tracks quasi-parallel to the disc sequence of the stellar mass-specific angular momentum plane, whereas retrograde minor mergers, major mergers, and bar-driven secular evolution move galaxies towards the bulge sequence. Finally, we demonstrate that fountain flows act to flatten and narrow the radial metallicity gradient and metallicity dispersion of disc stars, respectively. Thus, the evolution of galactic fountains depends strongly on the cosmological merger history and is crucial for the chemodynamical evolution of Milky-Way-sized disc galaxies. [ABSTRACT FROM AUTHOR]

We present shape measurements of Milky Way–sized dark matter haloes at redshift z  = 0 in a suite of 30 zoom simulations from the Auriga project. We compare the results in full magnetohydrodynamics against dark matter–only simulations and find a strong influence of baryons in making dark matter haloes rounder at all radii compared to their dark matter–only counterparts. At distances ≲30 kpc, rounder dark matter distributions correlate with extended massive stellar discs and low-core gas densities. We measure the alignment between the halo and the disc shapes at different radii and find a high degree of alignment at all radii for most of the galaxies. In some cases, the alignment significantly changes as a function of radius implying that the halo shape twists; this effect correlates with recently formed bulges and is almost absent in the dark matter–only simulations. In a comparison against observational constraints, we find that |$20{{\ \rm per\ cent}}$| of haloes in our sample are consistent with observational results derived from the Pal 5 stream that favours an almost spherical shape. Including baryons is a required element to achieve this level of agreement. In contrast, none of the simulations (neither dark matter only nor with baryons) match the constraints derived from the Sagittarius stream that favour an oblate dark matter halo. [ABSTRACT FROM AUTHOR]

We present a new cosmological, magnetohydrodynamical simulation for galaxy formation: TNG50, the third and final instalment of the IllustrisTNG project. TNG50 evolves 2 × 21603 dark matter particles and gas cells in a volume 50 comoving Mpc across. It hence reaches a numerical resolution typical of zoom-in simulations, with a baryonic element mass of |$8.5\times 10^4\, {\rm M}_{\odot }$| and an average cell size of 70–140 pc in the star-forming regions of galaxies. Simultaneously, TNG50 samples ∼700 (6500) galaxies with stellar masses above |$10^{10} \, (10^8)\, {\rm M}_{\odot }$| at |$z$|  = 1. Here we investigate the structural and kinematical evolution of star-forming galaxies across cosmic time (0 ≲  |$z$|  ≲ 6). We quantify their sizes, disc heights, 3D shapes, and degree of rotational versus dispersion-supported motions as traced by rest-frame V -band light (i.e. roughly stellar mass) and by |$\rm H\,\alpha$| light (i.e. star-forming and dense gas). The unprecedented resolution of TNG50 enables us to model galaxies with sub-kpc half-light radii and with ≲300-pc disc heights. Coupled with the large-volume statistics, we characterize a diverse, redshift- and mass-dependent structural and kinematical morphological mix of galaxies all the way to early epochs. Our model predicts that for star-forming galaxies the fraction of disc-like morphologies, based on 3D stellar shapes, increases with both cosmic time and galaxy stellar mass. Gas kinematics reveal that the vast majority of |$10^{9-11.5}\, {\rm M}_{\odot }$| star-forming galaxies are rotationally supported discs for most cosmic epochs (V rot/σ > 2–3, |$z$|  ≲ 5), being dynamically hotter at earlier epochs (⁠|$z$|  ≳ 1.5). Despite large velocity dispersion at high redshift, cold and dense gas in galaxies predominantly arranges in disky or elongated shapes at all times and masses; these gaseous components exhibit rotationally dominated motions far exceeding the collisionless stellar bodies. [ABSTRACT FROM AUTHOR]

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 FROM AUTHOR]

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 FROM AUTHOR]

We use the Illustris simulation to study the distributions of ellipticities and kinematic misalignments of galactic projections, as well as the intrinsic shapes and rotation of the simulated galaxies. Our results for the projections of galaxies display clear trends of an overall increase of kinematic misalignment and a slight decrease of ellipticity for fast rotators with increasing masses, while revealing no clear distinction between slow rotators of different mass. It is also found that the number of very slow rotators with large ellipticities is much larger than found in observations. The intrinsic properties of the galaxies are then analysed. The results for the intrinsic shapes of the galaxies are mostly consistent with previous results inferred from observational data. The distributions of intrinsic misalignment of the galaxies suggest that some of the galaxies produced by Illustris have significant rotation around their medium axes. Further analysis reveals that most of these galaxies display signs of non-equilibrium. We then study the evolution of the intrinsic misalignments and shapes of three specific Illustris galaxies, which we consider as typical ones, along the main progenitor line of their merger trees, revealing how mergers influence the intrinsic shapes and kinematics: the spin axis in general stays close to the shortest axis, and tends to quickly relax to such an equilibrium state within a few dynamical times of the galaxy after major perturbations; triaxiality and intrinsic flatness in general decrease with time, however, sometimes increases occur that are clearly seen to correlate with major merger events. [ABSTRACT FROM AUTHOR]

Galaxies are surrounded by massive gas reservoirs (i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated Milky Way-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L⋆ galaxies is extremely diverse: column densities of commonly observed species span ∼3 − 4 dex and their covering fractions range from |${\sim } 5$| to |$90{{\ \rm per\ cent}}$|⁠. Despite this diversity, we identify the following correlations: 1) the covering fractions (CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H  i , C  iv , and Si  ii anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H  i , C  iv , and Si  ii positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L⋆ galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM. [ABSTRACT FROM AUTHOR]

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 (C 82). S/T is correlated with stellar mass and star formation activity, while C 82 correlates only with stellar mass. Overall, we find good agreement with observational estimates for both S/T and C 82. 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 C 82) 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. [ABSTRACT FROM AUTHOR]

We select galaxies from the IllustrisTNG hydrodynamical simulations (⁠|$M_{\rm stars}\gt 10^9 \, {\rm M}_\odot$| at 0 ≤ z ≤ 2) and characterize the shapes and evolutions of their UVJ and star formation rate–stellar mass (SFR– M stars) 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 |$M_{\rm stars}\lt 10^{10.5} \, ({\rm M}_\odot)$| with a ∼0.3 dex scatter at z ∼ 0 in our fiducial choices. If a UVJ -based cut is adopted, the TNG MS exhibits a downwardly bending at stellar masses of about 1010.5−10.7 M⊙. Moreover, the model predicts that |${\sim }80\, (50)$|  per cent of 1010.5−11 M⊙ 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 M⊙yr−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. [ABSTRACT FROM AUTHOR]

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_{\rm x}\sim 10^{43-44}\, \rm erg/s$|⁠), and a lower number of bright AGN (⁠|$L_{\rm x}\gt 10^{44} \, \rm erg/s$|⁠), a conclusion that varies quantitatively but not qualitatively with BH luminosity estimation method. The lower Eddington ratios of the |$10^{9}\, \rm M_{\odot }$| 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_{\star }\geqslant 10^{11}\, \rm M_{\odot }$| 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{{\ \rm per\ cent}}$| at |$z$|  ∼ 1.5–2), and the lowest among compact quiescent galaxies (⁠|$6-10{{\ \rm per\ cent}}$| at |$z$|  ∼ 1.5–2). [ABSTRACT FROM AUTHOR]

The coevolution of galaxies and their metal content serves as an important test for galaxy feedback models. We analyse 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 0 < z < 2 and mass range 109 < M */M⊙ < 1010.5. We make predictions for the high-redshift (2 < z < 10) metal content of galaxies which is described by a gradual decline in the normalization of the metallicity with an average high-redshift (z > 2) evolution fit by d log(Z)/d z ≈ −0.064. Our simulations indicate that the metal retention efficiency of the interstellar medium (ISM) is low: a majority of gas-phase metals (∼85 per cent at z  = 0) live outside of the ISM, either in an extended gas disc, 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 slope of the correlation between metallicity and gas-mass at fixed stellar mass is not a strong function of redshift. Our results indicate that the gas fraction dependence of 'regulator' style models allows them to simultaneously explaining the shape, redshift evolution, and existence of correlated scatter with gas fraction about the MZR. [ABSTRACT FROM AUTHOR]

We explore orbit properties of 35 prolate-triaxial galaxies selected from the Illustris cosmological hydrodynamic simulation. We present a detailed study of their orbit families, and also analyse relations between the relative abundance of the orbit families and the spin parameter, triaxiality, the ratio of the angular momentum, and the baryon fraction. We find that box orbits dominate the orbit structure for most prolate-triaxial galaxies, especially in the central region. The fraction of irregular orbits in the prolate-triaxial galaxies is small, and it increases with galaxy radius. Both the x -tube and z -tube orbits are important in prolate-triaxial galaxies, especially in the outer regions. The fraction of box orbits for prolate-triaxial galaxies (0.7 < T < 1) decreases with the triaxiality of the stars, while the fraction of x -tube orbits increases with the triaxiality for a given axis ratio. The fraction of box orbits increases and the fractions of x -tube and z -tube orbits weakly decrease with increasing baryon fraction. These results help to understand the structure of prolate-triaxial galaxies and provide cross-checks for constructing dynamical models of prolate-triaxial galaxies by using the Schwarzschild or the made-to-measure methods. [ABSTRACT FROM AUTHOR]

The circumgalactic medium (CGM), i.e. the gaseous haloes around galaxies, is both the reservoir of gas that fuels galaxy growth and the repository of gas expelled by galactic winds. Most cosmological, hydrodynamical simulations focus their computational effort on the galaxies themselves and treat the CGM more coarsely, which means small-scale structure cannot be resolved. We get around this issue by running zoom-in simulations of a Milky Way-mass galaxy with standard mass refinement and additional uniform spatial refinement within the virial radius. This results in a detailed view of its gaseous halo at unprecedented (1 kpc) uniform resolution with only a moderate increase in computational time. The improved spatial resolution does not impact the central galaxy or the average density of the CGM. However, it drastically changes the radial profile of the neutral hydrogen column density, which is enhanced at galactocentric radii larger than 40 kpc. The covering fraction of Lyman-Limit Systems within 150 kpc is almost doubled. We therefore conclude that some of the observational properties of the CGM are strongly resolution dependent. Increasing the resolution in the CGM, without increasing the resolution of the galaxies, is a promising and computationally efficient method to push the boundaries of state-of-the-art simulations. [ABSTRACT FROM AUTHOR]

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 V max of 430 km s−1, contains a 100 kpc gas/stellar disc with morphology similar to Malin 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. [ABSTRACT FROM AUTHOR]

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 coincidentally 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 $$10^{8.5}\, {\rm M}_{\odot}$$, 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 FROM AUTHOR]

Observational studies have revealed that galaxy pairs tend to have lower gas-phase metallicity than isolated galaxies. This metallicity deficiency can be caused by inflows of low-metallicity gas due to the tidal forces and gravitational torques associated with galaxy mergers, diluting the metal content of the central region. In this work we demonstrate that such metallicity dilution occurs in state-of-the-art cosmological simulations of galaxy formation. We find that the dilution is typically 0.1 dex for major mergers, and is noticeable at projected separations smaller than 40 kpc. For minor mergers the metallicity dilution is still present, even though the amplitude is significantly smaller. Consistent with previous analysis of observed galaxies we find that mergers are outliers from the fundamental metallicity relation, with deviations being larger than expected for a Gaussian distribution of residuals. Our large sample of mergers within full cosmological simulations also makes it possible to estimate how the star formation rate enhancement and gas consumption timescale behave as a function of the merger mass ratio. We confirm that strong starbursts are likely to occur in major mergers, but they can also arise in minor mergers if more than two galaxies are participating in the interaction, a scenario that has largely been ignored in previous work based on idealised isolated merger simulations. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 stellar mass. 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. [ABSTRACT FROM AUTHOR]

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* < 1010 M⊙ 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* < 1010 M⊙ 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. [ABSTRACT FROM AUTHOR]

We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colours of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2 × 25003 resolution elements in a volume 20 times larger. Here, we present first results on the galaxy colour bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g - r) colours of 109 12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey.We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement with the observations, with a sharp transition in median colour from blue to red at a characteristic M* ∼ 1010.5☉. Investigating the build-up of the colour-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy colour transition is supermassive black hole feedback in its low accretion state. Across the entire population the median colour transition time-scale Δtgreen is ∼1.6Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, redder galaxies have lower star formation rates, gas fractions, and gas metallicities; their stellar populations are also older and their large-scale interstellar magnetic fields weaker than in bluer galaxies. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* >1011☉ which redden at z < 1 accumulate on average ∼25 per cent of their final z = 0 mass post-reddening; at the same time, ∼18 per cent of such massive galaxies acquire half or more of their final stellar mass while on the red sequence. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Galaxy mergers are believed to play a key role in transforming star-forming disc galaxies into quenched ellipticals. Most of our theoretical knowledge about such morphological transformations does, however, rely on idealized simulations where processes such as cooling of hot halo gas into the disc and gas accretion in the post-merger phase are not treated in a self-consistent cosmological fashion. In this paper, we study the morphological evolution of the stellar components of four major mergers occurring at z = 0.5 in cosmological hydrodynamical zoom simulations. In all simulations, the merger reduces the disc mass fraction, but all galaxies simulated at our highest resolution regrow a significant disc by z = 0 (with a disc fraction larger than 24 per cent). For runs with our default physics model, which includes galactic winds from star formation and black hole feedback, none of the merger remnants are quenched, but in a set of simulations with stronger black hole feedback, we find that major mergers can indeed quench galaxies. We conclude that major merger remnants commonly evolve into star-forming disc galaxies, unless sufficiently strong active galactic nucleus feedback assists in the quenching of the remnant. [ABSTRACT FROM AUTHOR]