Search

We study the implications of relaxing the requirement for ultralight axions to account for all dark matter in the Universe by examining mixed dark matter (MDM) cosmologies with axion fractions $f \leq 0.3$ within the fuzzy dark matter (FDM) window $10^{-25}$ eV $\lesssim m \lesssim 10^{-23}$ eV. Our simulations, using a new MDM gravity solver implemented in AxiREPO, capture wave dynamics across various scales with high accuracy down to redshifts $z\approx 1$. We identify halos with Rockstar using the CDM component and find good agreement of inferred halo mass functions (HMFs) and concentration-mass relations with theoretical models across redshifts $z=1-10$. This justifies our halo finder approach a posteriori as well as the assumptions underlying the MDM halo model AxionHMcode. Using the inferred axion halo mass - cold halo mass relation $M_{\text{a}}(M_{\text{c}})$ and calibrating a generalised smoothing parameter $\alpha$ to our MDM simulations, we present a new version of AxionHMcode. The code exhibits excellent agreement with simulations on scales $k< 20 \ h$ cMpc$^{-1}$ at redshifts $z=1-3.5$ for $f\leq 0.1$ around the fiducial axion mass $m = 10^{-24.5}$ eV $ = 3.16\times 10^{-25}$ eV, with maximum deviations remaining below 10%. For axion fractions $f\leq 0.3$, the model maintains accuracy with deviations under 20% at redshifts $z\approx 1$ and scales $k< 10 \ h$ cMpc$^{-1}$, though deviations can reach up to 30% for higher redshifts when $f=0.3$. Reducing the run-time for a single evaluation of AxionHMcode to below $1$ minute, these results highlight the potential of AxionHMcode to provide a robust framework for parameter sampling across MDM cosmologies in Bayesian constraint and forecast analyses., Comment: 20 pages, 10 figures, 5 Tables, comments welcome

The radial distribution of gas within galactic haloes is connected to the star formation rate and the nature of baryon-driven feedback processes. Using six variants of the hydrodynamic simulation Simba, we study the impact of different stellar/AGN feedback prescriptions on the gas density profiles of haloes in the total mass range $10^{11} \, \mathrm{M}_{\odot} < M_{\mathrm{200c}} < 10^{14} \, \mathrm{M}_{\odot}$ and redshift interval $0

According to the well established hierarchical framework for galaxy evolution, galaxies grow through mergers with other galaxies and the LambdaCDM cosmological model predicts that the stellar halos of galaxies are rich in remnants from minor mergers. The Stellar Streams Legacy Survey has provided a first release of a catalogue with a statistically significant sample of stellar streams in the Local Universe that can be used to study minor mergers and test the cosmological models. The main objective is to compare the results of the observations of stellar tidal streams with the predictions of state-of-the-art cosmological simulations regarding the formation of stellar streams up to a redshift z < 0.02, according to the LambdaCDM model. We use the predictions of the cosmological simulations Copernicus Complexio, TNG50 of the IllustrisTNG project and Auriga to generate 225 mock-images of nearby halos at a distance of 70 Mpc, and search for stellar streams. We compare the obtained stream frequency and characteristics with those obtained from the Stellar Streams Legacy Survey. We find good agreement between the results of analysing real images from the Dark Energy Survey and mock-images from cosmological simulations. We obtained predictions for the detection rate of stellar streams to a surface brightness limit of 35 mag arcsec^-2. The cosmological simulations predict that for a surface brightness limit of 32 mag arcsec^-2 a frequency of almost 70% in the detection of streams around galaxies can be achieved., Comment: 18 pages, 17 figures

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 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

A major aim of current and upcoming cosmological surveys is testing deviations from the standard $\Lambda$CDM model, but the full scientific value of these surveys will only be realised through efficient simulation methods that keep up with the increasing volume and precision of observational data. $N$-body simulations of modified gravity (MG) theories are computationally expensive since highly non-linear equations need to be solved to model the non-linear matter evolution; this represents a significant bottleneck in the path to reach the data volume and resolution attained by equivalent $\Lambda$CDM simulations. We develop a field-level, neural-network-based emulator that generates density and velocity divergence fields under the $f(R)$ gravity MG model from the corresponding $\Lambda$CDM simulated fields. Using attention mechanisms and a complementary frequency-based loss function, our model is able to learn this intricate mapping. We further use the idea of latent space extrapolation to generalise our emulator to $f(R)$ models with differing field strengths. The predictions of our emulator agree with the $f(R)$ simulations to within 5% for matter density and to within 10% for velocity divergence power spectra up to $k \sim 2\,h$ $\mathrm{Mpc}^{-1}$. But for a few select cases, higher-order statistics are reproduced with $\lesssim$10% agreement with the $f(R)$ simulations. Latent extrapolation allows our emulator to generalise to different $f(R)$ model variants without explicitly training on those variants. Given a $\Lambda$CDM simulation, the GPU-based emulator is able to reproduce the equivalent $f(R)$ realisation $\sim$600 times faster than full $N$-body simulations. This lays the foundations for a valuable tool for realistic yet rapid mock field generation and robust cosmological analyses., Comment: 20 pages, 10 figures. To be submitted to MNRAS

We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to define insitu stars. To uncover the underlying factors contributing to deviations between TNG50 and H3 distributions, we scrutinize correlation coefficients among internal drivers(e.g., virial radius, star formation rate [SFR]) and merger-related parameters (such as the effective mass-ratio, mean distance, average redshift, total number of mergers, average spin-ratio and maximum spin alignment between merging galaxies). Notably, we identify significant correlations between deviations from observational data and key parameters such as the median slope of virial radius, mean SFR values, and the rate of SFR change across different redshift scans. Furthermore, positive correlations emerge between deviations from observational data and parameters related to galaxy mergers. We validate these correlations using the Random Forest Regression method. Our findings underscore the invaluable insights provided by the H3 survey in unraveling the cosmic history of galaxies akin to the Milky Way, thereby advancing our understanding of galactic evolution and shedding light on the formation and evolution of Milky Way-like galaxies in cosmological simulations., Comment: 20 pages, 13 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

We present an in-depth investigation of galaxy clustering based on a new suite of realistic large-box galaxy-formation simulations in $f(R)$ gravity, with a subgrid physics model that has been recalibrated to reproduce various observed stellar and gas properties. We focus on the two-point correlation functions of the luminous red galaxies (LRGs) and emission line galaxies (ELGs), which are primary targets of ongoing and future galaxy surveys such as DESI. One surprising result is that, due to several nontrivial effects of modified gravity on matter clustering and the galaxy-halo connection, the clustering signal does not depend monotonically on the fifth-force strength. For LRGs this complicated behaviour poses a challenge to meaningfully constraining this model. For ELGs, in contrast, this can be straightforwardly explained by the time evolution of the fifth force, which means that weaker $f(R)$ models can display nearly the same -- up to $25\%$ -- deviations from $\Lambda$CDM as the strongest ones, albeit at lower redshifts. This implies that even very weak $f(R)$ models can be strongly constrained, unlike with most other observations. Our results show that galaxy formation acquires a significant environment dependence in $f(R)$ gravity which, if not properly accounted for, may lead to biased constraints on the model. This highlights the essential role of hydrodynamical simulations in future tests of gravity exploring precision galaxy-clustering data from the likes of DESI and Euclid., Comment: 18 pages, 8 figures + 2 figures in appendix

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

Observations with the James Webb Space Telescope have revealed a high abundance of bright galaxies at redshift, $z\gtrsim 12$, which has been widely interpreted as conflicting with the $\Lambda$CDM model. In Cowley et al. (2018) predictions were made - prior to the JWST observations - for the expected abundance of these galaxies using the Durham semi-analytic galaxy formation model, GALFORM, which is known to produce a realistic population of galaxies at lower redshifts including the present day. Key to this model is the assumption of a "top-heavy" initial mass function of stars formed in bursts (required to explain the number counts and redshift distribution of sub-millimetre galaxies). Here, we compare the rest-frame ultraviolet luminosity functions derived from JWST observations with those predicted by the Cowley et al. model up to $z=14$ and make further predictions for $z=16$. We find that below $z\sim 10$, the Cowley et al. predictions agree very well with observations, while agreement at $z\gtrsim12$ requires extending the model to take into account the timescale for the growth of obscuring dust grains and its dependence on gas metallicity. We trace the evolution of these galaxies from $z=14$ to $z=0$ and find that their descendants typically reside in halos with a median mass of $10^{13.6}\,h^{-1}\,\mathrm{M_{\odot}}$. The stellar masses of the descendants range from $10^{7}\,h^{-1}\,\mathrm{M_{\odot}}$ to $10^{11.5}\,h^{-1}\,\mathrm{M_{\odot}}$. Although these galaxies were all central galaxies at $z=14$, nearly half of their descendants end up as satellites in massive halos., Comment: 14 pages, 8 figures, submitted to MNRAS on 4 June, 2024

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

The goal of this study is to conduct a timely analysis of the high-redshift star-forming galaxy populations, which will be informative in designing next-generation experiments and their extragalactic targets. We use the hydrodynamical simulation MillenniumTNG (MTNG) to model Lyman-alpha Emitting (LAE) galaxies to extract key properties such as their clustering and occupation statistics. We define LAEs through an empirical relation between star formation rate (SFR) and Lyman-alpha flux. We also explore two other definitions, finding that imposing an additional cut on the maximum stellar mass of the galaxy sample, which approximates the effect of a low escape fraction at high halo mass, leads to a 5-10\% decrease of the linear bias of the population. As expected, we find that the HOD mass parameters rapidly decrease with increasing number density. Additionally, the HOD parameter $\sigma$ also decreases with number density, implying that the SFR-halo mass relationship becomes tighter for low-luminosity objects. Surprisingly, the non-linear clustering, estimated by the parameter $r_0$, is fixed at fixed number density, whereas the linear bias parameter varies with redshift as $b(z) \propto (1 + z)$, suggesting that our LAE samples are relatively stable and long-lived. Finally, we study the amount of galaxy assembly bias present at $z = 2, \ 3$ and find that while at $z = 2$ it is roughly $\lesssim$10\%, at $z = 3$ it decreases significantly to $\lesssim$5\%. This suggests that assembly bias effects become less important at high $z$ likely due to the lower number of cumulative two-halo interactions (mergers, splashback, stripping, etc.). While our study is based on a single full-physics simulation, we expect our results to reflect the properties of LAEs in the Universe. We demonstrate that our findings are in good agreement with previous results using both observations and simulations., Comment: 19 pages, 10 figures, 1 table

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