Search

We investigate the redshift evolution of the concentration-mass relationship of dark matter haloes in state-of-the-art cosmological hydrodynamic simulations and their dark-matter-only counterparts. By combining the IllustrisTNG suite and the novel MillenniumTNG simulation, our analysis encompasses a wide range of box size ($50 - 740 \: \rm cMpc$) and mass resolution ($8.5 \times 10^4 - 3.1 \times 10^7 \: \rm M_{\odot}$ per baryonic mass element). This enables us to study the impact of baryons on the concentration-mass relationship in the redshift interval $0

We present a comprehensive examination of the three most recent versions of the L-Galaxies semi-analytic galaxy formation model, focusing on the evolution of galaxy properties across a broad stellar mass range ($10^7\:{\rm M}_{\odot}\lesssim{M_\star}\lesssim10^{12}\:{\rm M}_{\odot}$) from $z=0$ to $z\simeq10$. We compare the predictions with the latest multiband data from key astronomical surveys, including SDSS, CANDELS, and COSMOS along with HST, JWST, and ALMA. We assess the models' ability to reproduce various time-dependent galaxy scaling relations for star-forming and quenched galaxies. Key focus areas include global galaxy properties such as stellar mass functions, cosmic star formation rate density, and the evolution of the main sequence of star-forming galaxies. Additionally, we examine resolved morphological properties such as the galaxy mass-size relation, alongside core $(R<1\,{\rm{kpc}})$ and effective $(R

We discuss the cold dark matter plus massive neutrinos simulations of the MillenniumTNG (MTNG) project, which aim to improve understanding of how well ongoing and future large-scale galaxy surveys will measure neutrino masses. Our largest simulations, $3000\,{\rm Mpc}$ on a side, use $10240^3$ particles of mass $m_{p} = 6.66\times 10^{8}\,h^{-1}{\rm M}_\odot$ to represent cold dark matter, and $2560^3$ to represent a population of neutrinos with summed mass $M_\nu = 100\,{\rm meV}$. Smaller volume runs with $\sim 630\,{\rm Mpc}$ also include cases with $M_\nu = 0\,\textrm{and}\, 300\,{\rm meV}$. All simulations are carried out twice using the paired-and-fixed technique for cosmic variance reduction. We evolve the neutrino component using the particle-based $\delta f$ importance sampling method, which greatly reduces shot noise in the neutrino density field. In addition, we modify the GADGET-4 code to account both for the influence of relativistic and mildly relativistic components on the expansion rate and for non-Newtonian effects on the largest represented simulation scales. This allows us to quantify accurately the impact of neutrinos on basic statistical measures of nonlinear structure formation, such as the matter power spectrum and the halo mass function. We use semi-analytic models of galaxy formation to predict the galaxy population and its clustering properties as a function of summed neutrino mass, finding significant ($\sim 10\%$) impacts on the cosmic star formation rate history, the galaxy mass function, and the clustering strength. This offers the prospect of identifying combinations of summary statistics that are optimally sensitive to the neutrino mass., Comment: 19 pages, 15 figures

The statistical analysis of cosmic large-scale structure is most often based on simple two-point summary statistics, like the power spectrum or the two-point correlation function of a sample of galaxies or other types of tracers. In contrast, topological measures of clustering are also sensitive to higher-order correlations, and thus offer the prospect to access additional information that may harbor important constraining power. We here revisit one such geometric measure of the cosmic web in the form of the so-called percolation analysis, using the recent MillenniumTNG simulation suite of the LCDM paradigm. We analyze continuum percolation statistics both for high resolution dark matter particle distributions, as well as for galaxy mock catalogues from a semi-analytic galaxy formation model within a periodic simulation volume of 3000 Mpc on a side. For comparison, we also investigate the percolation statistics of random particle sets and neutrino distributions with two different summed particle masses. We find that the percolation statistics of the dark matter distribution evolves strongly with redshift and thus clustering strength, yielding progressively lower percolation threshold towards later times. However, there is a sizable residual dependence on numerical resolution which we interpret as a residual influence of different levels of shot noise. This is corroborated by our analysis of galaxy mock catalogues whose results depend on sampling density more strongly than on galaxy selection criteria. While this limits the discriminative power of percolation statistics, our results suggest that it still remains useful as a complementary cosmological test when controlled for sampling density., Comment: 14 pages, 12 figures, accepted by ApJ, updated to match published version

Approximate N-body methods, such as FastPM and COLA, have been successful in modelling halo and galaxy clustering statistics, but their low resolution on small scales is a limitation for applications that require high precision. Full N-body simulations can provide better accuracy but are too computationally expensive for a quick exploration of cosmological parameters. This paper presents a method for correcting distinct haloes identified in fast N-body simulations, so that various halo statistics improve to a percent level accuracy. The scheme seeks to find empirical corrections to halo properties such that the virial mass is the same as that of a corresponding halo in a full N-body simulation. The modified outer density contour of the corrected halo is determined on the basis of the FastPM settings and the number of particles inside the halo. This method only changes some parameters of the halo finder, and does not require any extra CPU-cost. We demonstrate that the adjusted halo catalogues of FastPM simulations significantly improve the precision of halo mass-based statistics from redshifts $z=0.0$ to $1.0$, and that our calibration can be applied to different cosmologies without needing to be recalibrated., Comment: 11pages,15figures

Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed and the burstiness of the star formation history, are highly sensitive to numerical choices adopted in these subgrid models. Using the {\small SMUGGLE} stellar feedback model, we compute idealized simulations of a $M_{\rm vir} \sim 10^{10} \, \msun$ dwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of $\sim 3$. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies and a retention of cuspy dark matter density profiles. An improved understanding of the non-linear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment., Comment: Submitted ApJ; 15 pages, 12 figures; comments welcome

Weak gravitational lensing is a powerful tool for precision tests of cosmology. As the expected deflection angles are small, predictions based on non-linear N-body simulations are commonly computed with the Born approximation. Here we examine this assumption using ${\small DORIAN}$, a newly developed full-sky ray-tracing scheme applied to high-resolution mass-shell outputs of the two largest simulations in the MillenniumTNG suite, each with a 3000 Mpc box containing almost 1.1 trillion cold dark matter particles in addition to 16.7 billion particles representing massive neutrinos. We examine simple two-point statistics like the angular power spectrum of the convergence field, as well as statistics sensitive to higher order correlations such as peak and minimum statistics, void statistics, and Minkowski functionals of the convergence maps. Overall, we find only small differences between the Born approximation and a full ray-tracing treatment. While these are negligibly small at power-spectrum level, some higher order statistics show more sizable effects; ray-tracing is necessary to achieve percent level precision. At the resolution reached here, full-sky maps with 0.8 billion pixels and an angular resolution of 0.43 arcmin, we find that interpolation accuracy can introduce appreciable errors in ray-tracing results. We therefore implemented an interpolation method based on nonuniform fast Fourier transforms (NUFFT) along with more traditional methods. Bilinear interpolation introduces significant smoothing, while nearest grid point sampling agrees well with NUFFT, at least for our fiducial source redshift, $z_s=1.0$, and for the 1 arcmin smoothing we use for higher-order statistics., Comment: 13 pages, 7 figures, submitted to MNRAS

Galactic bars drive the internal evolution of spiral galaxies, while their formation is tightly coupled to the properties of their host galaxy and dark matter halo. To explore what drives bar formation in the cosmological context and how these structures evolve throughout cosmic history, we use the Auriga suite of magneto-hydrodynamical cosmological zoom-in simulations. We find that bars are robust and long-lived structures, and we recover a decreasing bar fraction with increasing redshift which plateaus around $\sim20\%$ at $z\sim3$. We find that bars which form at low and intermediate redshifts grow longer with time, while bars that form at high redshifts are born `saturated' in length, likely due to their merger-induced formation pathway. This leads to a larger bar-to-disc size ratio at high redshifts as compared to the local Universe. We subsequently examine the multi-dimensional parameter space thought to drive bar formation. We find that barred galaxies tend to have lower Toomre $Q$ values at the time of their formation, while we do not find a difference in the gas fraction of barred and unbarred populations when controlling for stellar mass. Barred galaxies tend to be more baryon-dominated at all redshifts, assembling their stellar mass earlier, while galaxies that are baryon-dominated but that do not host a bar, have a higher ex-situ bulge fraction. We explore the implications of the baryon-dominance of barred galaxies on the Tully-Fisher relation, finding an offset from the unbarred relation; confirming this in observations would serve as additional evidence for dark matter, as this behaviour is not readily explained in modified gravity scenarios., Comment: 21 pages including Appendices, 17 figures, submitted to MNRAS

Black holes are believed to be crucial in regulating star formation in massive galaxies, which makes it essential to faithfully represent the physics of these objects in cosmological hydrodynamics simulations. Limited spatial and mass resolution and the associated discreteness noise make following the dynamics of black holes especially challenging. In particular, dynamical friction, which is responsible for driving massive black holes towards the centres of galaxies, cannot be accurately modelled with softened $N$-body interactions. A number of subgrid models have been proposed to mimic dynamical friction or directly include its full effects in simulations. Each of these methods has its individual benefits and shortcomings, while all suffer from a common issue of being unable to represent black holes with masses below a few times the simulated dark matter particle mass. In this paper, we propose a correction for unresolved dynamical friction, which has been calibrated on simulations run with the code KETJU, in which gravitational interactions of black holes are not softened. We demonstrate that our correction is able to sink black holes with masses greater than the dark matter particle mass at the correct rate. We show that the impact of stochasticity is significant for low-mass black holes ($M_{\rm BH} \leq 5 M_{\rm DM}$) and propose a correction for stochastic heating. Combined, this approach is applicable to next generation cosmological hydrodynamics simulations that jointly track galaxy and black hole growth with realistic black hole orbits., Comment: 20 pages, 14 figures, accepted in MNRAS

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 HI 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 UVB 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 HI and MgII the column densities within $\sim 0.3\,R_{\rm vir}$ are reduced, while the OVI 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 AGN during phases of quasar activity and from soft X-ray sources, which have not yet been included in the simulations of the present study., Comment: 24 pages, 22 figures, accepted for the publication in MNRAS

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% 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., Comment: 19 pages, 13 figures

The Voids-within-Voids-within-Voids (VVV) project used dark-matter-only simulations to study the abundance and structure of dark matter haloes over the full mass range populated in the standard $\Lambda\mathrm{CDM}$ cosmology. Here we explore how baryonic effects modify these results for $z=0$ halo masses in the range $10^4$ to $10^7~\mathrm{M_\odot}$, below the threshold for galaxy formation. Our main study focuses on three simulations from identical initial conditions at $z=127$, one following dark matter only, one including non-radiative gas, and one additionally including the baryonic physics relevant in this halo mass range (cooling and photoheating). In the non-radiative simulation, above $10^{5.5}~\mathrm{M_\odot}$, halo abundance and internal structure are very similar to the dark-matter-only simulation, and the baryon to dark matter ratio is everywhere close to the cosmic value. At lower mass, this ratio drops and haloes are less concentrated and less massive in the non-radiative case. Test simulations at higher resolution show this to be mainly a resolution effect; the expected drop in baryon content due to residual pressure effects only becomes substantial for $z=0$ haloes below $\sim 10^{2.7}~\mathrm{M_\odot}$. However, gas is heated by reionization at $z=6$ in our "full physics" run, and this results in almost complete expulsion of gas from all haloes in our simulated mass range. This suppresses the halo mass function by $\sim 30 \%$, lowers halo concentration, and consequently weakens the dark matter annihilation signal by $\sim 40-60 \%$., Comment: 12+2 pages, 9 figures

We present a machine-learning-based model for the total density profiles of subhaloes with masses $M \gtrsim 7\times 10^8\,h^{-1}{\rm M}_\odot$ in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the density profiles themselves. The IVE returns accurate and unbiased predictions on all radial ranges, including the outer region profile where the subhaloes experience tidal stripping; here its fit accuracy exceeds that of the commonly used Einasto profile. The IVE discovers three independent degrees of freedom in the profiles, which can be interpreted in terms of the formation history of the subhaloes. In addition to the two parameters controlling the normalization and inner shape of the profile, the IVE discovers a third parameter that accounts for the impact of tidal stripping onto the subhalo outer profile; this parameter is sensitive to the mass loss experienced by the subhalo after its infall onto its parent halo. Baryonic physics in the IllustrisTNG galaxy formation model does not impact the number of degrees of freedom identified in the profile compared to the pure dark matter expectations, nor their physical interpretation. Our newly proposed profile fit can be used in strong lensing analyses or other observational studies which aim to constrain cosmology from small-scale structures., Comment: 13 pages, 7 figures. Minor changes to match version accepted for publication in MNRAS

Dwarf satellites in galaxy groups are distributed in an anisotropic and asymmetric manner, which is called the ``lopsided satellite distribution''. This lopsided signal has been observed not only in galaxy pairs but also in isolated systems. However, the physical origin of the lopsided signal in isolated systems is still unknown. In this work, we investigate this in the state-of-the-art hydrodynamical simulation of the MillenniumTNG Project by tracing each system back to high redshift. We find that the lopsided signal is dominated by satellites located in the outer regions of the halo and is also dominated by recently accreted satellites. The lopsided signal originates from the anisotropic accretion of galaxies from the surrounding large-scale structure and that, after accretion, the nonlinear evolution of satellites inside the dark-matter halo weakens the lopsidedness. The signal decreases as cosmic time passes because of a competition between anisotropic accretion and internal evolution within dark matter halos. Our findings provide a useful perspective for the study of galaxy evolution, especially for the origin of the spatial satellite galaxy distributions., Comment: 10 pages, 8 figures, 1 table. Submitted to MNRAS

For decades, the boundary of cosmic filaments have been a subject of debate. In this work, we determine the physically-motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at z = 0 in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from z = 4 to z = 1, and then slightly increases. Moreover, we find that the filament radius depends on the filament length, the distance from connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experience rapid contraction before z = 1, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves., Comment: 12 pages, 11 figures, accepted by mnras

The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above $10^{14}\rm{M_\odot}$, we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of $\sim1\rm{kpc}$ in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of $1.48\rm{kpc}$. We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using gas as a tracer for the potential can overestimate offsets by factors of $\sim10$ and $\sim30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section., Comment: 16 pages, 6 figures; v3: OJA published

We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical "zoom-in" simulations of 40 Milky Way-mass halos and 26 dwarf galaxy-mass halos run with the moving-mesh code Arepo. Auriga adopts the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) 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 main sequence, and HI gas fraction/disc thickness. Finally, we show that star-forming 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 utilise the released data., Comment: Matches final MNRAS version