Search

We study weak gravitational lensing convergence maps produced from the MillenniumTNG (MTNG) 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., 14 pages, 10 figures, comments welcome

In the coming years, Sunyaev-Zel'dovich (SZ) measurements can dramatically improve our understanding of the Intergalactic Medium (IGM) and the role of feedback processes on 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, $\tau$, and the Compton-y parameter, $Y$, receive large line-of-sight contributions which can be removed effectively by applying a Compensated Aperture Photometry (CAP) filter. In contrast with other $\tau$ probes (e.g., X-rays and Fast Radio Bursts), the kSZ-inferred $\tau$ 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 $\tau-M$ scaling relations and report best-fit parameters adopting the smoothly broken power law (SBPL) formalism. We note that subgrid physics modeling can broaden the error bar on these by 30\% for intermediate-mass halos ($\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 halos in real observations, which can bias the inferred SZ profiles by $\sim$20\% for $L_\ast$-galaxies., Comment: 14 pages, 6 figures

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 \textsc{FORGE} and \textsc{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 $\Omega_{\mathrm{m0}}, H_0, \sigma_8$, and one MG parameter, either $\bar{f}_{R0}$ or $r_{\mathrm{c}}$. Our halo property emulators achieve an accuracy of $\lesssim 1\%$ 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., Comment: 18 pages, 15 figures

Wide-area deep imaging surveys have discovered large numbers of extremely low surface brightness dwarf galaxies, which challenge galaxy formation theory and, potentially, offer new constraints on the nature of dark matter. Here we discuss one as-yet unexplored formation mechanism that may account for a fraction of low surface brightness dwarfs. We call this the `ghost galaxy' scenario. In this scenario, inefficient radiative cooling prevents star formation in the `main branch' of the merger tree of a low mass dark matter halo, such that almost all its stellar mass is acquired through mergers with less massive (but nevertheless star-forming) progenitors. Present-day systems formed in this way would be `ghostly' isolated stellar halos with no central galaxy. We use merger trees based on the Extended Press-Schechter formalism and the COCO cosmological N-body simulation to demonstrate that mass assembly histories of this kind can occur for low-mass halos in Lambda-CDM, but they are rare. They are most probable in isolated halos of present-day mass ~4x10^9 M_sun, occurring for ~5 per cent of all halos of that mass under standard assumptions about the timing and effect of cosmic reionization. The stellar masses of star-forming progenitors in these systems are highly uncertain; abundance-matching arguments imply a bimodal present-day mass function having a brighter population (median M_star ~3x10^6 M_sun) consistent with the tail of the observed luminosity function of ultra-diffuse galaxies. This suggests observable analogues of these systems may await discovery. We find that a stronger ionizing background (globally or locally) produces brighter and more extended ghost galaxies., Comment: 19 pages, 13 figures, submitted to ApJ

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 \ {\rm Mpc}/h$, 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 modeling 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 subhalos and particles in the simulation. We find that ELGs are anisotropically distributed within halos, 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., Comment: 14 pages, 10 figures; submitted to MNRAS

Approximate methods to populate dark matter halos with galaxies are of great utility to large 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\sim 0.1{\rm Mpc}/h$, while here we focus on modeling galaxy assembly bias in the two-halo regime, $r\gtrsim 1{\rm Mpc}/h$. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales ($r\sim 20{\rm Mpc}/h$), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10-15\% discrepancy present in the standard halo model prescription 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-neighbor counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters., Comment: 16 pages, 8 figures; submitted to MNRAS

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5\mu$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $\Lambda$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $2\mu$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($\Delta z\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical., Comment: Submitted to ApJL. Fig. 1 shows the "superposed" state, Fig. 2 presents the potential z~5 protocluster, Fig. 3 shows some atypical z~5 SEDs, Fig. 4 studies the morphology, Fig. 5 summarizes the extraordinary situation at z~17, and Fig. 6 argues for why such interlopers may not be a major issue for z>10 searches. Comments warmly welcomed

We show that gravitational waves have the potential to unravel the microphysical properties of dark matter due to the dependence of the binary black hole merger rate on cosmic structure formation, which is itself highly dependent on the dark matter scenario. In particular, we demonstrate that suppression of small-scale structure -- such as that caused by interacting, warm, or fuzzy dark matter -- leads to a significant reduction in the rate of binary black hole mergers at redshifts $z\gtrsim5$. This shows that future gravitational-wave observations will provide a new probe of physics beyond the $\Lambda$CDM model., Comment: 15 pages, 12 figures, matches version published in PRD

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 $\approx15$ 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 $\approx10$ Mpc, there is significant scale-dependent bias on smaller scales. The amplitude of this effect varies between the two galaxy types, and also 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., Comment: 15 pages, 11 figures. Accepted for publication in MNRAS

The early release science results from $\textit{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\gtrsim8$) galaxy population and compare them to recent $\textit{JWST}$ observations. We show that the simulated galaxy population is broadly consistent with observations until $z\sim10$. From $z\approx10-12$, the observations indicate a preference for a galaxy population that is largely dust-free, but is still consistent with the simulations. Beyond $z\gtrsim12$, 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 $\textit{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 $\Lambda$CDM model. We emphasise that any new process invoked to explain this tension should only significantly influence the galaxy population beyond $z\gtrsim10$, while leaving the successful galaxy formation predictions of the fiducial model intact below this redshift., Comment: Accepted for publication in MNRAS -- Part of the initial set of papers introducing the MillenniumTNG project. Visit www.mtng-project.org for more details

We describe a method for generating halo catalogues on the light cone using the \Abacus{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 \Abacus{AbacusSummit} halo light cone catalogues. Our recommended use of this product is in the halo mass regime of $M_{\rm halo} > 2.1 \times 10^{11} \ M_\odot/h$ for the \texttt{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 \sim 1$ by adopting a modified version of the \Abacus{AbacusHOD} script, which builds on the standard halo occupation distribution (HOD) method by including various extensions. We find that the multipoles of the auto-correlation 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 subpercent level. Halo light cone catalogues for 25 \texttt{base} and 2 \texttt{huge} simulations at the fiducial cosmology is available at DOI:\href{https://www.doi.org/10.13139/OLCF/1825069}{10.13139/OLCF/1825069}, 16 pages, 11 figures, accepted by MNRAS

We perform a series of cosmological hydrodynamic simulations to investigate the effects of non-gravitational dark matter (DM) interactions on the intergalactic medium (IGM). In particular, we use the ETHOS framework (Cyr-Racine et al. 2016; Vogelsberger et al. 2016) to compare statistics of the Lyman-$��$ forest in cold dark matter (CDM) with an alternative model in which the DM couples strongly with a relativistic species in the early universe. These models are characterised by a cutoff in the linear power spectrum, followed by a series of 'dark acoustic oscillations' (DAOs) on sub-dwarf scales. While the primordial cutoff delays the formation of the first galaxies, structure builds-up more rapidly in the interacting DM model compared to CDM. We show that although DAOs are quickly washed away in the non-linear clustering of DM at $z\lesssim10$, their signature can be imprinted prominently in the Lyman-$��$ flux power spectrum at $z>5$. On scales larger than the cutoff ($k\sim0.08$ s/km for the specific model considered here), the relative difference to CDM is reminiscent of a warm dark matter (WDM) model with a similar initial cutoff; however, the redshift evolution on smaller scales is distinctly different. The appearance and disappearance of DAOs in the Lyman-$��$ flux spectrum provides a powerful way to distinguish interacting DM models from WDM and, indeed, variations in the thermal history of the IGM that may also induce a small-scale cutoff., 15 pages, 12 figures, 1 table. Updated version accepted for publication in MNRAS

We explore the connection between the UV luminosity functions (LFs) of high-$z$ galaxies and the distribution of stellar masses and star-formation histories (SFHs) in their host dark matter halos. We provide a baseline for a redshift-independent star-formation efficiency model to which observations and models can be compared. Our model assigns a star-formation rate (SFR) to each dark matter halo based on the growth rate of the halo and a redshift-independent star-formation efficiency. The dark matter halo accretion rate is obtained from a high-resolution $N$-body simulation in order to capture the stochasticity in accretion histories and to obtain spatial information for the distribution of galaxies. The halo mass dependence of the star-formation efficiency is calibrated at $z=4$ by requiring a match to the observed UV LF at this redshift. The model then correctly predicts the observed UV LF at $z=5-10$. We present predictions for the UV luminosity and stellar mass functions, \textit{JWST} number counts, and SFHs. In particular, we find a stellar-to-halo mass relation at $z=4-10$ that scales with halo mass at $M_{\rm h}, 28 pages with 17 figures + appendix, accepted for publication in ApJ