No Arabic abstract
We present an analysis of the Counts-in-Cells (CiC) statistics of subhaloes in the publicly available IllustrisTNG cosmological simulations (TNG100-1, TNG100-3 and TNG300-3), considering their full and dark-on
We present the full public release of all data from the TNG50, TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~1.1 PB, including 2,000 full volume snapshots and ~110,000 high time-resolution subbox snapshots. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets.
We study the distribution of cold dark matter (CDM) in cosmological simulations from the FIRE (Feedback In Realistic Environments) project, for $M_{ast}sim10^{4-11},M_{odot}$ galaxies in $M_{rm h}sim10^{9-12},M_{odot}$ halos. FIRE incorporates explicit stellar feedback in the multi-phase ISM, with energetics from stellar population models. We find that stellar feedback, without fine-tuned parameters, greatly alleviates small-scale problems in CDM. Feedback causes bursts of star formation and outflows, altering the DM distribution. As a result, the inner slope of the DM halo profile ($alpha$) shows a strong mass dependence: profiles are shallow at $M_{rm h}sim10^{10}-10^{11},M_{odot}$ and steepen at higher/lower masses. The resulting core sizes and slopes are consistent with observations. This is broadly consistent with previous work using simpler feedback schemes, but we find steeper mass dependence of $alpha$, and relatively late growth of cores. Because the star formation efficiency $M_{ast}/M_{rm h}$ is strongly halo mass dependent, a rapid change in $alpha$ occurs around $M_{rm h}sim 10^{10},M_{odot}$ ($M_{ast}sim10^{6}-10^{7},M_{odot}$), as sufficient feedback energy becomes available to perturb the DM. Large cores are not established during the period of rapid growth of halos because of ongoing DM mass accumulation. Instead, cores require several bursts of star formation after the rapid buildup has completed. Stellar feedback dramatically reduces circular velocities in the inner kpc of massive dwarfs; this could be sufficient to explain the Too Big To Fail problem without invoking non-standard DM. Finally, feedback and baryonic contraction in Milky Way-mass halos produce DM profiles slightly shallower than the Navarro-Frenk-White profile, consistent with the normalization of the observed Tully-Fisher relation.
We use the IllustrisTNG simulations to show how the fractions of quenched galaxies vary across different environments and cosmic time, and to quantify the role AGN feedback and preprocessing play in quenching group and cluster satellites. At $z=0$, we select galaxies with $M_* = 10^{9-12} M_{odot}$ residing within ($leq R_{200c}$) groups and clusters of total host mass $M_{200c}=10^{13-15.2} M_{odot}$. TNG predicts a quenched fraction of $sim70-90%$ (on average) for centrals and satellites $gtrsim 10^{10.5} M_{odot}$, regardless of host mass, cosmic time ($0leq zleq0.5$), clustercentric distance and time since infall in the $z=0$ host. Low-mass centrals ($lesssim 10^{10} M_{odot}$), instead, are rarely quenched unless they become members of groups ($10^{13-14} M_{odot}$) or clusters ($geq10^{14} M_{odot}$), where the quenched fraction rises to $sim80%$. The fraction of low-mass passive galaxies is higher closer to the host center and for more massive hosts. The population of low-mass satellites accreted $gtrsim$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, $sim30%$ 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 to 10 billion years ago. For low-mass galaxies ($lesssim10^{10-10.5}M_{odot}$), this is due to preprocessing, 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.
Hydrodynamical simulations of galaxy formation have now reached sufficient volume to make precision predictions for clustering on cosmologically relevant scales. Here we use our new IllustrisTNG simulations to study the non-linear correlation functions and power spectra of baryons, dark matter, galaxies and haloes over an exceptionally large range of scales. We find that baryonic effects increase the clustering of dark matter on small scales and damp the total matter power spectrum on scales up to k ~ 10 h/Mpc by 20%. The non-linear two-point correlation function of the stellar mass is close to a power-law over a wide range of scales and approximately invariant in time from very high redshift to the present. The two-point correlation function of the simulated galaxies agrees well with SDSS at its mean redshift z ~ 0.1, both as a function of stellar mass and when split according to galaxy colour, apart from a mild excess in the clustering of red galaxies in the stellar mass range 10^9-10^10 Msun/h^2. Given this agreement, the TNG simulations can make valuable theoretical predictions for the clustering bias of different galaxy samples. We find that the clustering length of the galaxy auto-correlation function depends strongly on stellar mass and redshift. Its power-law slope gamma is nearly invariant with stellar mass, but declines from gamma ~ 1.8 at redshift z=0 to gamma ~ 1.6 at redshift z ~ 1, beyond which the slope steepens again. We detect significant scale-dependencies in the bias of different observational tracers of large-scale structure, extending well into the range of the baryonic acoustic oscillations and causing nominal (yet fortunately correctable) shifts of the acoustic peaks of around ~5%.
The IllustrisTNG project is a new suite of cosmological magneto-hydrodynamical simulations of galaxy formation performed with the Arepo code and updated models for feedback physics. Here we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters ($10^{13} leq M_{rm 200c}/M_{rm sun} leq 10^{15}$) at recent times ($z leq 1$). The richest clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intra-cluster light. The exact ICL fraction depends sensitively on the definition of a central galaxys mass and varies in our most massive clusters between 20 to 40% of the total stellar mass. Haloes of $5times 10^{14}M_{rm sun}$ and above have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial mass density distribution as shallow as the dark matters ( $-3.5 < alpha_{rm 3D} < -3$). Total halo mass is a very good predictor of stellar mass, and vice versa: at $z=0$, the 3D stellar mass measured within 30 kpc scales as $propto (M_{rm 500c})^{0.49}$ with a $sim 0.12$ dex scatter. This is possibly too steep in comparison to the available observational constraints, even though the abundance of TNG less massive galaxies ($< 10^{11}M_{rm sun}$ in stars) is in good agreement with the measured galaxy stellar mass functions at recent epochs. The 3D sizes of massive galaxies fall too on a tight ($sim$0.16 dex scatter) power-law relation with halo mass, with $r^{rm stars}_{rm 0.5} propto (M_{rm 500c})^{0.53}$. Even more fundamentally, halo mass alone is a good predictor for the whole stellar mass profiles beyond the inner few kpc, and we show how on average these can be precisely recovered given a single mass measurement of the galaxy or its halo.