Search

We use the voids-within-voids-within-voids (VVV) simulations, a suite of successive nested N -body simulations with extremely high resolution (denoted, from low to high resolution, by L0 to L7), to test the Press–Schechter (PS), Sheth–Tormen (ST), and extended Press–Schechter (EPS) formulae for the halo abundance over the entire mass range, from minihaloes of 10−6 M⊙, to cluster haloes of 1015 M⊙, at different redshifts, from z = 30 to the present. We find that at z = 0 and z = 2, ST best reproduces the results of L0, which has the mean cosmic density (overdensity δ = 0), at 1011−15 M⊙. The higher resolution levels (L1–L7) are biased underdense regions (δ < −0.6). The EPS formalism takes this into account since it gives the mass function of a region conditioned, in this case, on having a given underdensity. EPS provides good matches to these higher levels, with deviations ≲20 per cent, at 10−6−12.5 M⊙. At z ∼ 7−15, the ST predictions for L0 and the EPS for L1–L7 show somewhat larger deviations from the simulation results. However, at even higher redshifts, z ∼ 30, EPS fits the simulations well again. We confirm our results by picking more subvolumes from the L0 simulation, finding that our conclusions depend only weakly on the size and overdensity of the region. The good agreement of EPS with the higher level simulations implies that PS (or ST) gives an accurate description of the total halo mass function in representative regions of the universe. [ABSTRACT FROM AUTHOR]

Modern redshift surveys are tasked with mapping out the galaxy distribution over enormous distance scales. Existing hydrodynamical simulations, however, do not reach the volumes needed to match upcoming surveys. We present results for the clustering of galaxies using a new, large volume hydrodynamical simulation as part of the MillenniumTNG (MTNG) project. With a computational volume that is ≈15 times larger than the next largest such simulation currently available, we show that MTNG is able to accurately reproduce the observed clustering of galaxies as a function of stellar mass. When separated by colour, there are some discrepancies with respect to the observed population, which can be attributed to the quenching of satellite galaxies in our model. We combine MTNG galaxies with those generated using a semi-analytic model to emulate the sample selection of luminous red galaxies (LRGs) and emission-line galaxies (ELGs) and show that, although the bias of these populations is approximately (but not exactly) constant on scales larger than ≈10 Mpc, there is significant scale-dependent bias on smaller scales. The amplitude of this effect varies between the two galaxy types and between the semi-analytic model and MTNG. We show that this is related to the distribution of haloes hosting LRGs and ELGs. Using mock SDSS-like catalogues generated on MTNG lightcones, we demonstrate the existence of prominent baryonic acoustic features in the large-scale galaxy clustering. We also demonstrate the presence of realistic redshift space distortions in our mocks, finding excellent agreement with the multipoles of the redshift-space clustering measured in SDSS data. [ABSTRACT FROM AUTHOR]

The stellar halo of the Milky Way records the history of its interactions with dwarf galaxies, whose subsequent destruction results in the formation of an extended stellar component. Recent works have suggested that galaxies with masses comparable to the Large Magellanic Cloud (LMC, |$M_\star \sim 10^9$| M☉) may be the primary building blocks of the stellar halo of our Galaxy. We use cosmological simulations of the Lambda cold dark matter model to investigate LMC-mass galaxies at |$z$| = 1–2 using a semi-analytical model of galaxy formation. We find that LMC analogues at |$z$| = 2 evolve until the present day along three distinct pathways: (i) those that are destroyed in Milky Way-mass hosts; (ii) those that are themselves the main progenitors of Milky Way-mass galaxies; and (iii) those that survive until |$z$| = 0, with stellar mass ∼1.0 dex lower than typical Milky Ways. We predict that the properties of these galaxies at |$z$| = 2 (stellar metallicities, sizes, gas content, etc.) are largely indistinguishable, irrespective of which of these pathways is eventually taken; a survey targeting such galaxies in this redshift range would struggle to tell apart a 'destroyed' stellar halo progenitor from a 'surviving' LMC analogue. The only factor that determines the eventual fate of these galaxies is their proximity to a neighbouring Milky Way main progenitor at |$z$| = 2: while the mean separation to a 'surviving' galaxy is around 7 Mpc, it is only 670 kpc to a 'destroyed' galaxy. This suggests that old stellar populations in the Milky Way may share intrinsic (i.e. non-dynamical) properties that are essentially indistinguishable from progenitors of its stellar halo. [ABSTRACT FROM AUTHOR]

Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in cold dark matter (CDM) and alternative DM models which have varying relativistic coupling and self-interaction strengths. The interacting DM scenarios produce a cutoff in the linear power spectrum on small-scales, followed by a series of 'dark acoustic oscillations'. We find that DM interactions suppress the abundance of galaxies below |$M_\star \sim 10^8\, {\rm M}_\odot$| for the models considered. The cutoff in the power spectrum delays structure formation relative to CDM. Objects in ETHOS that end up at the same final masses as their CDM counterparts are characterized by a more vigorous phase of early star formation. While galaxies with |$M_\star \lesssim 10^6\, {\rm M_\odot }$| make up more than 60 per cent of star formation in CDM at z ≈ 10, they contribute only about half the star formation density in ETHOS. These differences diminish with decreasing redshift. We find that the effects of DM self-interactions are negligible compared to effects of relativistic coupling (i.e. the effective initial conditions for galaxy formation) in all properties of the galaxy population we examine. Finally, we show that the clustering strength of galaxies at high redshifts depends sensitively on DM physics, although these differences are manifest on scales that may be too small to be measurable by JWST. [ABSTRACT FROM AUTHOR]

Tracking the formation and evolution of dark matter haloes is a critical aspect of any analysis of cosmological N -body simulations. In particular, the mass assembly of a halo and its progenitors, encapsulated in the form of its merger tree, serves as a fundamental input for constructing semi-analytic models of galaxy formation and, more generally, for building mock catalogues that emulate galaxy surveys. We present an algorithm for constructing halo merger trees from abacussummit , the largest suite of cosmological N -body simulations performed to date consisting of nearly 60 trillion particles, and which has been designed to meet the Cosmological Simulation Requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. Our method tracks the cores of haloes to determine associations between objects across multiple time slices, yielding lists of halo progenitors and descendants for the several tens of billions of haloes identified across the entire suite. We present an application of these merger trees as a means to enhance the fidelity of abacussummit halo catalogues by flagging and 'merging' haloes deemed to exhibit non-monotonic past merger histories. We show that this cleaning technique identifies portions of the halo population that have been deblended due to choices made by the halo finder, but which could have feasibly been part of larger aggregate systems. We demonstrate that by cleaning halo catalogues in this post-processing step, we remove potentially unphysical features in the default halo catalogues, leaving behind a more robust halo population that can be used to create highly accurate mock galaxy realizations from abacussummit. [ABSTRACT FROM AUTHOR]

We use N -body cosmological simulations and empirical galaxy models to study the merger history of dwarf-mass galaxies (with |$M_{\rm halo}\sim 10^{10}\, \mathrm{M}_\odot$|⁠). Our input galaxy models describe the stellar mass–halo mass relation, and the galaxy occupation fraction. The number of major and minor mergers depends on the type of dark matter; in particular, minor mergers are greatly suppressed in warm dark matter models. In addition, the number of mergers that bring in stars is strongly dependent on the galaxy occupation model. For example, minor mergers are negligible for stellar halo growth in models with a high mass threshold for galaxy formation (i.e. |$10^{9.3}\, \mathrm{M}_\odot$| at z = 0). Moreover, this threshold for galaxy formation can also determine the relative difference (if any) between the stellar haloes of satellite and field dwarfs. Using isolated simulations of dwarf–dwarf mergers, we show that the relative frequency of major and minor mergers predict very different stellar haloes: Typically, 'intermediate' dark matter merger ratios (∼1:5) maximize the growth of distant stellar haloes. We discuss the observability of dwarf stellar haloes and find that the surface brightness of these features are incredibly faint. However, when several dwarfs are stacked together, models that form particularly rich stellar haloes could be detectable. Finally, we show that stellar streams in the Galactic halo overlapping in phase space with known dwarf satellites are likely remnants of their stripped stellar haloes. The mere existence of dwarf stellar haloes can already put constraints on some small-scale models, and thus observational probes should be a high priority. [ABSTRACT FROM AUTHOR]

The James Webb Space Telescope (JWST) will revolutionize our understanding of early galaxy formation, and could potentially set stringent constraints on the nature of dark matter. We use a semi-empirical model of galaxy formation to investigate the extent to which uncertainties in the implementation of baryonic physics may be degenerate with the predictions of two different models of dark matter – cold dark matter (CDM) and a 7 keV sterile neutrino, which behaves as warm dark matter (WDM). Our models are calibrated to the observed UV luminosity function at z  = 4 using two separate dust attenuation prescriptions, which manifest as high and low star formation efficiency in low-mass haloes. These efficiencies capture the net effect of processes that regulate star formation. We find that while at fixed star formation efficiency, ε, there are marked differences in the abundance of faint galaxies in the two dark matter models at high- z ; these differences are mimicked easily by varying ε in the same dark matter model. We find that a high ε WDM and a low ε CDM model – which provide equally good fits to the z  = 4 UV luminosity function – exhibit nearly identical evolution in the cosmic stellar mass and star formation rate densities. We show that differences in the star formation rate at fixed stellar mass are larger for variations in ε in a given dark matter model than they are between dark matter models; however, the scatter in star formation rates is larger between the two models than they are when varying ε. Our results suggest that JWST will likely be more informative in constraining baryonic processes operating in high- z galaxies than it will be in constraining the nature of dark matter. [ABSTRACT FROM AUTHOR]

We take advantage of the unprecedented dynamical range provided by the 'Cosmic-Zoom' project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the Lambda cold dark matter paradigm when the dark matter is made of weakly interacting massive particles of mass 100 GeV. In particular, we complement previous studies by exploring the MAHs of haloes with mass from |$10^8\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$| down to Earth mass, |$10^{-6}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠. The formation redshift of low-mass haloes anticorrelates weakly with mass, peaking at z = 3 for haloes of mass |$10^{-4}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠. Even lower masses are affected by the free-streaming cut-off in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press–Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological N -body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, |$M\gt 10^{10}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠ , the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than 20 per cent at high redshift. If this result can be extrapolated to lower halo masses, we conclude that EPS theory predicts the hierarchical build up of dark matter haloes quite well over the entire halo mass range. [ABSTRACT FROM AUTHOR]

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 non-linear 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 haloes. Our findings provide a useful perspective for the study of galaxy evolution, especially for the origin of the spatial satellite galaxy distributions. [ABSTRACT FROM AUTHOR]

We employ the hydrodynamical simulation illustrisTNG to inform the galaxy–halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ∼ 0.8. Specifically, we model the galaxy colours of illustrisTNG and apply sliding DESI colour–magnitude cuts, matching the DESI target densities. We study the halo occupation distribution (HOD) model of the selected samples by matching them to their corresponding dark matter haloes in the illustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as 'galaxy assembly bias'. The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modelling galaxy–halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10–15  |${{\ \rm per\ cent}}$| excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures. [ABSTRACT FROM AUTHOR]

We describe a new method (compaso) for identifying groups of particles in cosmological N -body simulations. compaso builds upon existing spherical overdensity (SO) algorithms by taking into consideration the tidal radius around a smaller halo before competitively assigning halo membership to the particles. In this way, the compaso finder allows for more effective deblending of haloes in close proximity as well as the formation of new haloes on the outskirts of larger ones. This halo-finding algorithm is used in the abacussummit suite of N -body simulations, designed to meet the cosmological simulation requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. compaso is developed as a highly efficient on-the-fly group finder, which is crucial for enabling good load-balancing between the GPU and CPU and the creation of high-resolution merger trees. In this paper, we describe the halo-finding procedure and its particular implementation in abacus , accompanying it with a qualitative analysis of the finder. We test the robustness of the compaso catalogues before and after applying the cleaning method described in an accompanying paper and demonstrate its effectiveness by comparing it with other validation techniques. We then visualize the haloes and their density profiles, finding that they are well fit by the NFW formalism. Finally, we compare other properties such as radius–mass relationships and two-point correlation functions with that of another widely used halo finder, rockstar. [ABSTRACT FROM AUTHOR]

In this series of papers, we present an emulator-based halo model for the non-linear clustering of galaxies in modified gravity cosmologies. In the first paper, we present emulators for the following halo properties: the halo mass function, concentration–mass relation and halo-matter cross-correlation function. The emulators are trained on data extracted from the forge and bridge suites of N -body simulations, respectively, for two modified gravity (MG) theories: f (R) gravity, and the DGP model, varying three standard cosmological parameters Ωm0, H 0, σ8, and one MG parameter, either |$\bar{f}_{R0}$| or r c. Our halo property emulators achieve an accuracy of |${\lesssim}1\ \hbox{per cent}$| on independent test data sets. We demonstrate that the emulators can be combined with a galaxy–halo connection prescription to accurately predict the galaxy–galaxy and galaxy–matter correlation functions using the halo model framework. [ABSTRACT FROM AUTHOR]

Sunyaev–Zel'dovich (SZ) measurements can dramatically improve our understanding of the intergalactic medium and the role of feedback processes in galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, τ, and the Compton- y parameter, Y , receive large line-of-sight contributions that can be removed effectively by applying a compensated aperture photometry filter. In contrast with other τ probes (e.g. X-rays and fast radio bursts), the kinematic SZ-inferred τ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line of sight. Additionally, we perform fits to the Y – M and τ– M scaling relations and report best-fitting parameters adopting the smoothly broken power law formalism. We note that subgrid physics modelling can broaden the error bar on these by 30 per cent for intermediate-mass haloes (⁠|${\sim }10^{13} \, {\rm M}_{\odot }$|⁠). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than haloes in observations, which can bias the inferred profiles by ∼20 per cent for L * galaxies. [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]

We study weak gravitational lensing convergence maps produced from the MillenniumTNG simulations by direct projection of the mass distribution on the past backwards lightcone of a fiducial observer. We explore the lensing maps over a large dynamic range in simulation mass and angular resolution, allowing us to establish a clear assessment of numerical convergence. By comparing full physics hydrodynamical simulations with corresponding dark-matter-only runs, we quantify the impact of baryonic physics on the most important weak lensing statistics. Likewise, we predict the impact of massive neutrinos reliably far into the non-linear regime. We also demonstrate that the 'fixed & paired' variance suppression technique increases the statistical robustness of the simulation predictions on large scales not only for time slices but also for continuously output lightcone data. We find that both baryonic and neutrino effects substantially impact weak lensing shear measurements, with the latter dominating over the former on large angular scales. Thus, both effects must explicitly be included to obtain sufficiently accurate predictions for stage IV lensing surveys. Reassuringly, our results agree accurately with other simulation results where available, supporting the promise of simulation modelling for precision cosmology far into the non-linear regime. [ABSTRACT FROM AUTHOR]

Approximate methods to populate dark-matter haloes with galaxies are of great utility to galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs, z  = 1 and z  = 0, in the large-volume, high-resolution hydrodynamical simulation of the MillenniumTNG project. In a partner study we concentrated on the small-scale, one-halo regime down to r ∼ 0.1  h −1 Mpc, while here we focus on modelling galaxy assembly bias in the two-halo regime, r ≳ 1  h −1 Mpc. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales (r ∼ 20  h −1 Mpc), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10–15 per cent discrepancy is only reconciled when we augment our halo occupation model with a dependence on extrinsic halo properties ('shear' being the best-performing one) rather than intrinsic ones (e.g. concentration, peak mass). We argue that this fact constitutes evidence for two-halo galaxy conformity. Including tertiary assembly bias (i.e. a property beyond mass and 'shear') is not an essential requirement for reconciling the galaxy assembly bias signal of LRGs, but the combination of external and internal properties is beneficial for recovering ELG the clustering. We find that centrals in low-mass haloes dominate the assembly bias signal of both populations. Finally, we explore the predictions of our model for higher order statistics such as nearest neighbour counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters. [ABSTRACT FROM AUTHOR]

Luminous red galaxies (LRGs) and blue star-forming emission-line galaxies (ELGs) are key tracers of large-scale structure used by cosmological surveys. Theoretical predictions for such data are often done via simplistic models for the galaxy–halo connection. In this work, we use the large, high-fidelity hydrodynamical simulation of the MillenniumTNG project (MTNG) to inform a new phenomenological approach for obtaining an accurate and flexible galaxy-halo model on small scales. Our aim is to study LRGs and ELGs at two distinct epochs, z  = 1 and z  = 0, and recover their clustering down to very small scales, |$r \sim 0.1 \ h^{-1}\, {\rm Mpc}$|⁠ , i.e. the one-halo regime, while a companion paper extends this to a two-halo model for larger distances. The occupation statistics of ELGs in MTNG inform us that (1) the satellite occupations exhibit a slightly super-Poisson distribution, contrary to commonly made assumptions, and (2) that haloes containing at least one ELG satellite are twice as likely to host a central ELG. We propose simple recipes for modelling these effects, each of which calls for the addition of a single free parameter to simpler halo occupation models. To construct a reliable satellite population model, we explore the LRG and ELG satellite radial and velocity distributions and compare them with those of subhaloes and particles in the simulation. We find that ELGs are anisotropically distributed within haloes, which together with our occupation results provides strong evidence for cooperative galaxy formation (manifesting itself as one-halo galaxy conformity); i.e. galaxies with similar properties form in close proximity to each other. Our refined galaxy-halo model represents a useful improvement of commonly used analysis tools and thus can be of help to increase the constraining power of large-scale structure surveys. [ABSTRACT FROM AUTHOR]

Cosmological inference with large galaxy surveys requires theoretical models that combine precise predictions for large-scale structure with robust and flexible galaxy formation modelling throughout a sufficiently large cosmic volume. Here, we introduce the millenniumTNG (MTNG) project which combines the hydrodynamical galaxy formation model of illustrisTNG with the large volume of the millennium simulation. Our largest hydrodynamic simulation, covering |$(500 \, h^{-1}{\rm Mpc})^3 \simeq (740\, {\rm Mpc})^3$|⁠ , is complemented by a suite of dark-matter-only simulations with up to 43203 dark matter particles (a mass resolution of |$1.32\times 10^8 \, h^{-1}{\rm M}_\odot$|⁠) using the fixed-and-paired technique to reduce large-scale cosmic variance. The hydro simulation adds 43203 gas cells, achieving a baryonic mass resolution of |$2\times 10^7 \, h^{-1}{\rm M}_\odot$|⁠. High time-resolution merger trees and direct light-cone outputs facilitate the construction of a new generation of semi-analytic galaxy formation models that can be calibrated against both the hydro simulation and observation, and then applied to even larger volumes – MTNG includes a flagship simulation with 1.1 trillion dark matter particles and massive neutrinos in a volume of |$(3000\, {\rm Mpc})^3$|⁠. In this introductory analysis we carry out convergence tests on basic measures of non-linear clustering such as the matter power spectrum, the halo mass function and halo clustering, and we compare simulation predictions to those from current cosmological emulators. We also use our simulations to study matter and halo statistics, such as halo bias and clustering at the baryonic acoustic oscillation scale. Finally we measure the impact of baryonic physics on the matter and halo distributions. [ABSTRACT FROM AUTHOR]

We introduce a novel technique for constraining cosmological parameters and galaxy assembly bias using non-linear redshift-space clustering of galaxies. We scale cosmological N -body simulations and insert galaxies with the SubHalo Abundance Matching extended (SHAMe) empirical model to generate over 175 000 clustering measurements spanning all relevant cosmological and SHAMe parameter values. We then build an emulator capable of reproducing the projected galaxy correlation function at the monopole, quadrupole, and hexadecapole level for separations between |$0.1\, h^{-1}\, {\rm Mpc}$| and |$25\, h^{-1}\, {\rm Mpc}$|⁠. We test this approach by using the emulator and Monte Carlo Markov Chain (MCMC) inference to jointly estimate cosmology and assembly bias parameters both for the MTNG740 hydrodynamic simulation and for a semi-analytical model (SAM) galaxy formation built on the MTNG740-DM dark matter-only simulation, obtaining unbiased results for all cosmological parameters. For instance, for MTNG740 and a galaxy number density of |$n\sim 0.01 h^{3}\, {\rm Mpc}^{-3}$|⁠ , we obtain |$\sigma _{8}=0.799^{+0.039}_{-0.044}$| and |$\Omega _\mathrm{M}h^2= 0.138^{+ 0.025}_{- 0.018}$| (which are within 0.4 and 0.2σ of the MTNG cosmology). For fixed Hubble parameter (h), the constraint becomes |$\Omega _\mathrm{M}h^2= 0.137^{+ 0.011}_{- 0.012}$|⁠. Our method performs similarly well for the SAM and for other tested sample densities. We almost always recover the true amount of galaxy assembly bias within 1σ. The best constraints are obtained when scales smaller than |$2\, h^{-1}\, {\rm Mpc}$| are included, as well as when at least the projected correlation function and the monopole are incorporated. These methods offer a powerful way to constrain cosmological parameters using galaxy surveys. [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]

Cosmological simulations are an important theoretical pillar for understanding non-linear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large-volume simulations of cosmic structure formation aiming to better understand physical processes on large scales and to help interpret upcoming large-scale galaxy surveys. We here focus on the full physics box MTNG740 that computes a volume of |$740\, \mathrm{Mpc}^3$| with a baryonic mass resolution of |$3.1\times ~10^7\, \mathrm{M_\odot }$| using arepo with 80.6 billion cells and the IllustrisTNG galaxy formation model. We verify that the galaxy properties produced by MTNG740 are consistent with the TNG simulations, including more recent observations. We focus on galaxy clusters and analyse cluster scaling relations and radial profiles. We show that both are broadly consistent with various observational constraints. We demonstrate that the SZ-signal on a deep light-cone is consistent with Planck limits. Finally, we compare MTNG740 clusters with galaxy clusters found in Planck and the SDSS-8 RedMaPPer richness catalogue in observational space, finding very good agreement as well. However, simultaneously matching cluster masses, richness, and Compton- y requires us to assume that the SZ mass estimates for Planck clusters are underestimated by 0.2 dex on average. Due to its unprecedented volume for a high-resolution hydrodynamical calculation, the MTNG740 simulation offers rich possibilities to study baryons in galaxies, galaxy clusters, and in large-scale structure, and in particular their impact on upcoming large cosmological surveys. [ABSTRACT FROM AUTHOR]

In this series of papers, we present a simulation-based model for the non-linear clustering of galaxies based on separate modelling of clustering in real space and velocity statistics. In the first paper, we present an emulator for the real-space correlation function of galaxies, whereas the emulator of the real-to-redshift space mapping based on velocity statistics is presented in the second paper. Here, we show that a neural network emulator for real-space galaxy clustering trained on data extracted from the dark quest suite of N-body simulations achieves sub-per cent accuracies on scales 1 < r < 30 |$h^{-1} \, \mathrm{Mpc}$|⁠ , and better than 3 per cent on scales r < 1 |$h^{-1}\, \mathrm{Mpc}$| in predicting the clustering of dark-matter haloes with number density 10−3.5 |$(h^{-1}\, \mathrm{Mpc})^{-3}$|⁠ , close to that of SDSS LOWZ-like galaxies. The halo emulator can be combined with a galaxy–halo connection model to predict the galaxy correlation function through the halo model. We demonstrate that we accurately recover the cosmological and galaxy–halo connection parameters when galaxy clustering depends only on the mass of the galaxies' host halos. Furthermore, the constraining power in σ8 increases by about a factor of 2 when including scales smaller than 5 |$h^{-1} \, \mathrm{Mpc}$|⁠. However, when mass is not the only property responsible for galaxy clustering, as observed in hydrodynamical or semi-analytic models of galaxy formation, our emulator gives biased constraints on σ8. This bias disappears when small scales (r < 10 |$h^{-1}\, \mathrm{Mpc}$|⁠) are excluded from the analysis. This shows that a vanilla halo model could introduce biases into the analysis of future data sets. [ABSTRACT FROM AUTHOR]

We investigate cosmological structure formation in fuzzy dark matter (FDM) with the attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons. [ABSTRACT FROM AUTHOR]

To extract information from the clustering of galaxies on non-linear scales, we need to model the connection between galaxies and haloes accurately and in a flexible manner. Standard halo occupation distribution (HOD) models make the assumption that the galaxy occupation in a halo is a function of only its mass, however, in reality; the occupation can depend on various other parameters including halo concentration, assembly history, environment, and spin. Using the IllustrisTNG hydrodynamical simulation as our target, we show that machine learning tools can be used to capture this high-dimensional dependence and provide more accurate galaxy occupation models. Specifically, we use a random forest regressor to identify which secondary halo parameters best model the galaxy–halo connection and symbolic regression to augment the standard HOD model with simple equations capturing the dependence on those parameters, namely the local environmental overdensity and shear, at the location of a halo. This not only provides insights into the galaxy formation relationship but also, more importantly, improves the clustering statistics of the modelled galaxies significantly. Our approach demonstrates that machine learning tools can help us better understand and model the galaxy–halo connection, and are therefore useful for galaxy formation and cosmology studies from upcoming galaxy surveys. [ABSTRACT FROM AUTHOR]

We introduce the AbacusHOD model and present two applications of AbacusHOD and the AbacusSummit simulations to observations. AbacusHOD is a Halo Occupation Distribution (HOD) framework written in Python that is particle-based, multitracer, highly generalized, and highly efficient. It is designed specifically with multitracer/cosmology analyses for next-generation large-scale structure surveys in mind, and takes advantage of the volume and precision offered by the new state-of-the-art AbacusSummit cosmological simulations. The model is also highly customizable and should be broadly applicable to any upcoming surveys and a diverse range of cosmological analyses. In this paper, we demonstrate the capabilities of the AbacusHOD framework through two example applications. The first example demonstrates the high efficiency and the large HOD extension feature set through an analysis of full-shape redshift-space clustering of BOSS galaxies at intermediate to small scales (⁠|$\lt 30\, h^{-1}$|  Mpc), assessing the necessity of introducing secondary galaxy biases (assembly bias). We find strong evidence for using halo environment instead of concentration to trace secondary galaxy bias, a result which also leads to a moderate reduction in the 'lensing is low' tension. The second example demonstrates the multitracer capabilities of the AbacusHOD package through an analysis of the extended Baryon Oscillation Spectroscopic Survey cross-correlation measurements between three different galaxy tracers: luminous red galaxies, emission-line galaxies, and quasi-stellar objects. We expect the AbacusHOD framework, in combination with the AbacusSummit simulation suite, to play an important role in a simulation-based analysis of the upcoming Dark Energy Spectroscopic Instrument data sets. [ABSTRACT FROM AUTHOR]

We describe a method for generating halo catalogues on the light-cone using the AbacusSummit suite of N -body simulations. The main application of these catalogues is the construction of realistic mock galaxy catalogues and weak lensing maps on the sky. Our algorithm associates the haloes from a set of coarsely spaced snapshots with their positions at the time of light-cone crossing by matching halo particles to on-the-fly light-cone particles. It then records the halo and particle information into an easily accessible product, which we call the AbacusSummit halo light-cone catalogues. Our recommended use of this product is in the halo mass regime of M halo > 2.1 × 1011  M ⊙  h −1 for the base resolution simulations, i.e. haloes containing at least 100 particles, where the interpolated halo properties are most reliable. To test the validity of the obtained catalogues, we perform various visual inspections and consistency checks. In particular, we construct galaxy mock catalogues of emission-line galaxies (ELGs) at z ∼ 1 by adopting a modified version of the AbacusHOD script, which builds on the standard halo occupation distribution (HOD) method by including various extensions. We find that the multipoles of the autocorrelation function are consistent with the predictions from the full-box snapshot, implicitly validating our algorithm. In addition, we compute and output CMB convergence maps and find that the auto- and cross-power spectrum agrees with the theoretical prediction at the sub-per-cent level. Halo light-cone catalogues for 25 base and 2 huge simulations at the fiducial cosmology are available at DOI:10.13139/OLCF/1825069 [ABSTRACT FROM AUTHOR]

We present the public data release of the AbacusSummit cosmological N -body simulation suite, produced with the Abacus   N -body code on the Summit supercomputer of the Oak Ridge Leadership Computing Facility. Abacus achieves |$\mathcal {O}(10^{-5})$| median fractional force error at superlative speeds, calculating 70M particle updates per second per node at early times, and 45M particle updates per second per node at late times. The simulation suite totals roughly 60 trillion particles, the core of which is a set of 139 simulations with particle mass |$2\times 10^{9}\, h^{-1}\, \mathrm{M}_\odot$| in box size |$2\, h^{-1}\, \mathrm{Gpc}$|⁠. The suite spans 97 cosmological models, including Planck 2018, previous flagship simulation cosmologies, and a linear derivative and cosmic emulator grid. A subsuite of 1883 boxes of size |$500\, h^{-1}\, \mathrm{Mpc}$| is available for covariance estimation. AbacusSummit data products span 33 epochs from z  = 8 to 0.1 and include light cones, full particle snapshots, halo catalogues, and particle subsets sampled consistently across redshift. AbacusSummit is the largest high-accuracy cosmological N -body data set produced to date. [ABSTRACT FROM AUTHOR]