Halo Distribution Research Articles

Charged fuzzy dark matter black holes

We investigate the impact of fuzzy dark matter (FDM) on supermassive black holes (SMBHs) characterized by a spherical charge distribution. This work introduces a new class of spherically symmetric, self-gravitational relativistic charged models for FDM haloes, using the Einasto density model. This study enables the dark matter (DM) to appear as the matter ingredient, which constructs the black hole and extends the non-commutative mini black hole stellar solutions. By considering the charged anisotropic energy–momentum tensor with an equation of state (EoS) pr=−ρ, we explore various black hole solutions for different values of the Einasto index and mass parameter. Our approach suggests that the central density of the resulting black hole model mimics the usual de Sitter core. Furthermore, we discuss the possibility of constructing a charged self-gravitational droplet by replacing the above-mentioned EoS with a non-local one. However, under these circumstances, the radial pressure is observed to be negative. Ultimately, we consider various possibilities of constructing DM black holes, featuring intermediate masses that could evolve into galaxies. Consequently, some of these theoretical models have the potential to replace the usual black hole solutions of the galactic core. Simultaneously, these models are physically beneficial for being comprised of the fundamental matter component of the cosmos. Due to the outcomes of this paper, we would be able to study the connection between BH and DM by formulating stable stellar structures featuring fuzzy mass distributions derived from the Einasto distribution of DM halos.

How Do Primordial Black Holes Change the Halo Mass Function and Structure?

We examine the effects of massive primordial black holes (PBHs) on cosmic structure formation, employing both a semianalytical approach and cosmological simulations. Our simulations incorporate PBHs with a monochromatic mass distribution centered around 106 M ⊙, constituting a fraction of 10−2 to 10−4 of the dark matter (DM) in the Universe, with the remainder being collisionless particle DM. Additionally, we conduct a ΛCDM simulation for comparative analysis with runs that include PBHs. At smaller scales, halos containing PBHs exhibit similar density and velocity dispersion profiles to those without PBHs. Conversely, at larger scales, PBHs can expedite the formation of massive halos and reside at their centers owing to the “seed effect.” To analyze the relative distribution of PBH host halos compared to non-PBH halos, we apply nearest neighbor statistics. Our results suggest that PBH host halos, through gravitational influence, significantly impact the structure formation process, compared to the ΛCDM case, by attracting and engulfing nearby newly formed minihalos. Should PBHs constitute a fraction of DM significantly larger than ∼10−3, almost all newly formed halos will be absorbed by PBH-seeded halos. Consequently, our simulations predict a bimodal feature in the halo mass function, with most of the massive halos containing at least one PBH at their core and the rest being less massive non-PBH halos.

Secondary halo bias through cosmic time II. Reconstructing halo properties using clustering information

When constructing mock galaxy catalogs based on suites of dark matter halo catalogs generated with approximated, calibrated, or machine-learning approaches, assigning intrinsic properties for these tracers is a step of paramount importance, given that they can shape the abundance and spatial distribution of mock galaxies and galaxy clusters. We explore the possibility of assigning properties of dark matter halos within the context of calibrated or learning approaches, explicitly using clustering information. The goal is to retrieve the correct signal of primary and secondary large-scale effective bias as a function of properties reconstructed solely based on phase-space properties of the halo distribution and dark matter density field. The algorithm reconstructs a set of halo properties (such as virial mass, maximum circular velocity, concentration, and spin) constrained to reproduce both primary and secondary (or assembly) bias. The key ingredients of the algorithm are the implementation of individually-assigned large-scale effective bias, a multi-scale approach to account for halo exclusion, and a hierarchical assignment of halo properties. The method facilitates the assignment of halo properties, aiming to replicate the large-scale effective bias, both primary and secondary. This constitutes an improvement over previous methods in the literature, especially for the high-mass end population. We have designed a strategy for reconstructing the main properties of dark matter halos obtained using calibrated or learning algorithms, such that the one- and two-point statistics (on large scales) replicate the signal from detailed $N$-body simulations. We encourage the application of this strategy (or the implementation of our algorithm) for the generation of mock catalogs of dark matter halos based on approximated methods.

Fast simulation mapping: From standard to modified gravity cosmologies using the bias assignment method

Context. We assess the effectiveness of a non-parametric bias model in generating mock halo catalogues for modified gravity (MG) cosmologies, relying on the distribution of dark matter from either MG or Λ cold dark matter (ΛCDM) simulations. Aims. We aim to generate halo catalogues that effectively capture the distinct impact of MG, ensuring high accuracy in both two- and three-point statistics for a comprehensive analysis of large-scale structures. We investigated the inclusion of MG in non-local bias to directly map the tracers onto ΛCDM fields, which would significantly reduce computational costs. Methods. We employed the bias assignment method (BAM) to model halo distribution statistics by leveraging seven high-resolution COLA simulations of MG cosmologies. Taking cosmic-web dependences into account when learning the bias relations, we designed two experiments to map the MG effects: one utilising the consistent MG density fields and the other employing the benchmark ΛCDM density field. Results. BAM generates MG halo catalogues from both calibration experiments with excellent summary statistics, achieving a ~1% accuracy in the power spectrum across a wide range of k modes, with minimal differences well below 10% for modes subject to cosmic variance, particularly below k < 0.07 h Mpc−1. The reduced bispectrum remains consistent with the reference catalogues within 10% for the studied configuration. Our results demonstrate that a non-linear and non-local bias description can model the effects of MG starting from a ΛCDM field.

The intrinsic distribution of Lyman-alpha halos

The emission and escape of Lyman-alpha photons from star-forming galaxies is determined through complex interactions between the emitted photons and a galaxy's interstellar and circumgalactic gas.\ This causes Lyman-alpha emitters (LAEs) to commonly appear not as point sources but in spatially extended halos with complex spectral profiles. We developed a 3D spatial-spectral model of Lyman-alpha halos (LAHs) to replicate LAH observations in integral field spectroscopic studies, such as those made with VLT/MUSE. The profile of this model is a function of six key halo properties: the halo- and compact-source exponential scale lengths sH $ and $r_ sC $), the halo flux fraction ($f_H$), the compact component ellipticity ($q$), the spectral line width (sigma ), and the spectral line skewness parameter (gamma ). Placing a series of model LAHs into datacubes that reflect observing conditions in the MUSE UDF-Mosaic survey, we tested their detection recoverability and determine that sigma , $r_ sH $, and $f_H$ are expected to have the most significant effect on the detectability of the overall LAH at a given central wavelength and intrinsic line luminosity. We developed a general selection function model that spans a grid of these halo parameters. Using it with a sample of 145 LAHs with measured halo properties observed in the UDF-Mosaic survey, we derived completeness-corrected, intrinsic distributions of the values of sigma , $r_ sH $, and $f_H$ for $3<z<5$ LAHs. We present the best-fit functional forms of the distributions as well as a sigma distribution corrected for instrumental line-spread function broadening, and thereby show the physical line-spread distribution of the intrinsic population. Finally, we discuss possible implications for these distributions for the nature of Lyalpha emission through the circumgalactic medium, finding that observations may undercount LAHs with extended halo scale lengths compared to the intrinsic population.

Identifying Halos in Cosmological Simulations with Continuous Wavelet Analysis: The 2D Case

Continuous wavelet analysis is gaining popularity in science and engineering for its ability to analyze data across both spatial and scale domains simultaneously. In this study, we introduce a wavelet-based method for halo identification and assess its feasibility in two-dimensional (2D) scenarios. We begin by generating four pseudo-2D data sets from the SIMBA dark matter simulation by compressing thin slices of three-dimensional (3D) data into 2D. Subsequently, we calculate the continuous wavelet transform (CWT) directly from the particle distributions, identify local maxima that represent actual halos, and segment the CWT to delineate halo boundaries. A comparison with the traditional friends-of-friends (FOF) method shows that while our CWT-identified halos contain slightly fewer particles, they have smoother boundaries and are more compact in dense regions. In contrast, the CWT method can link particles over greater distances to form halos in sparse regions due to its spatial segmentation scheme. The spatial distribution and halo power spectrum of both CWT and FOF halos demonstrate substantial consistency, validating the 2D applicability of CWT for halo detection. While our identification scheme has been tested solely in a limited 2D context and has shown some performance limitations, its linear time complexity of O(N) and consistency with the FOF method suggest its suitability for analyzing significantly larger data sets in the future and its potential to be extended to the 3D case.

Integrating multispectral remote sensing and geological investigation for gold prospecting in the Borongo-Mborguene gold field, eastern Cameroon

Multispectral remote sensing data has great promise for mineral prospecting via regional-scale mapping of alteration zones and geological structures, once the method is calibrated using in-situ geological observations. The Borongo-Mborguene area, located in eastern Cameroon at the transition from sub-tropical to arid climates, hosts several alluvial gold deposits. In this investigation, maps developed by processing of Landsat-8 OLI images, ASTER satellite images, and SRTM DEM data were compared to results of a comprehensive field survey to locate zones of hydrothermal alteration and regional structural lineaments, and to verify the association of these features with gold mineralization. Specific Landsat-8 and ASTER band ratios, as well as principal component analysis of the multispectral images, were selected to isolate spectral signals that characterize particular groups of alteration minerals, such as clays, hydrated carbonates, hydrated sulphates, iron oxides, and iron hydroxides. Lineaments and their preferred orientations were extracted by filtering of the digital elevation data. A fuzzy logic approach was subsequently applied to the ensemble of resulting maps to highlight five areas of concentrated hydrothermal alteration and structural complexity as potential gold prospects (Mborguene, Boutila, Nandoungue, Borongo and Ouaden villages). Petrographic study of samples collected from the remotely-identified alteration zones revealed a hydrothermal ore assemblage of gold, galena, pyrite and hematite in mineralized quartz veins. Structural control of the alteration zones in the Borongo-Mborguene gold field was confirmed by field mapping along a shear zone. Hence, our petrographic and field studies indicate that the multispectral remote sensing data and data processing methods developed herein are appropriate for mineral exploration in remote sub-tropical areas of Africa. More generally, the computationally efficient and simple procedure developed here should yield accurate preliminary maps of the spatial distribution of hydrothermal halos and related mineral deposits in challenging mixed-type vegetated terrain around the globe.

Clustering of dark matter in the cosmic web as a probe of massive neutrinos

ABSTRACT The large-scale structure of the Universe is distributed in a cosmic web. Studying the distribution and clustering of dark matter particles and haloes may open up a new horizon for studying the physics of the dark Universe. In this work, we investigate the nearest neighbour statistics and spherical contact function in cosmological models with massive neutrinos. For this task, we use the relativistic N-body code, gevolution, and study particle snapshots at three different redshifts. In each snapshot, we find the haloes and evaluate the letter functions for them. We show that a generic behaviour can be found in the nearest neighbour, G(r), and spherical contact functions, F(r), which makes these statistics promising tools to constrain the total neutrino mass.

Secondary halo bias through cosmic time

Context. The spatial distribution of dark matter halos carries cosmological and astrophysical information. Cosmological information can be considered to be contained in the connection between halo main properties and the large-scale halo bias, while the astrophysical information would be encoded in the scaling relations between halo properties. The combination of these two contributions leads to the effect of secondary halo bias. Aims. Our goal is to measure the signal of secondary halo bias as a function of a variety of intrinsic and environmental halo properties and to characterize its statistical significance as a function of cosmological redshift. Methods. Using fixed and paired N-body simulations of dark-matter halos – the UNIT simulation – with masses above ∼1011 M⊙h−1 identified over a wide range of cosmological redshifts (0 < z < 5), we explored the behavior of the scaling relations among different halo properties. We included novel environmental properties based on the halo distribution as well as the underlying dark-matter field. We implemented an object-by-object estimator of large-scale effective bias and tested its validity against standard approaches. With a bias assigned to each tracer, we performed a statistical analysis aimed at characterizing the distribution of the bias and the signal of the secondary halo bias. Results. We show how the halo scaling relations linking direct probes of the halo potential well do not depend on the environment. On the contrary, links between the halo mass and the so-called set of secondary halo properties are sensitive to the cosmological environment, mainly to under-dense regions. We show that the signal of secondary bias is derived statistically from secondary correlations beyond the standard link to the halo mass. Conclusions. We show that the secondary bias arises through nonlocal and/or environmental properties related either to the halo distribution or to the properties of the underlying dark-matter field. In particular, properties such as the tidal field (a measure of the anisotropy of the density field) and the local Mach number (a measure of the local kinetic temperature of the halo distribution) generate the signals of the secondary bias with the highest significance. We propose applications of the assignment of individual bias for the generation of mock catalogs containing the signal of secondary bias, as well as a series of cosmological analyses aimed at mining large galaxy datasets.

Stellar halo density with LAMOST K and M giants

Aims. We derive the morphology of the stellar component in the outer halo volume, and search for possible overdensities due to substructures therein. Methods. We made use of some of the data releases of the spectroscopic survey LAMOST DR8-DR9 in tandem with distance determinations for two subsamples, that is, of K-giants and M-giants, respectively, making up 60 000 stars. These distance are obtained through Bayesian techniques that derive absolute magnitudes as a function of measured spectroscopic parameters. Our calculation of the density from these catalogues requires: (1) derivation of the selection function; and (2) a correction for the convolution of the distance errors, which we carried out with Lucy’s inversion of the corresponding integral equation. Results. The stellar density distribution of the outer halo (distance to the Galactic centre, rG, of between 25 and 90 kpc) is a smooth monotonously decreasing function with a dependence of approximately ρ ∝ rG−n, with n = 4.6 ± 0.4 for K-giants and n = 4.5 ± 0.2 for M-giants, and with a insignificant oblateness. The value of n is independent of the angular distance to the Sagittarius tidal stream plane, which is what would be expected if such a stream did not exist in the anticenter positions or had a negligible imprint in the density distribution in the outer halo. Apart from random fluctuations or minor anomalies in some lines of sight, we do not see substructures superimposed in the outer halo volume within the resolution that we are using and limited by the error bars. This constrains the mass of over- and under-densities in the outer halo to be of ≲103 M⊙ deg−2, whereas the total mass of the stellar halo, including inner and outer parts, is ∼7 × 108 M⊙.

Topological bias: how haloes trace structural patterns in the cosmic web

ABSTRACT We trace the connectivity of the cosmic web as defined by haloes in the Planck-Millennium simulation using a persistence and Betti curve analysis. We normalize clustering up to the second-order correlation function and use our systematic topological analysis to correlate local information and properties of haloes with their multiscale geometrical environment of the cosmic web (elongated filamentary bridges and sheetlike walls). We capture the multiscale topology traced by the halo distribution through filtrations of the corresponding Delaunay tessellation. The resulting nested alpha shapes are sensitive to the local density, perfectly outline the local geometry, and contain the complete information on the multiscale topology. We find a remarkable linear relationship between halo masses and topology: haloes of different mass trace environments with different topological signature. This is a topological bias, an environmental structure bias independent of the halo clustering bias associated with the two-point correlation function. This mass-dependent linear scaling relation allows us to take clustering into account and determine the overall connectivity from a limited sample of galaxies. The presence of topological bias has major implications for the study of voids and filaments in the observed distribution of galaxies. The (infra)structure and shape of these key cosmic web components will strongly depend on the underlying galaxy sample. Their use as cosmological probes, with their properties influenced by cosmological parameters, will have to account for the subtleties of topological bias. This is of particular relevance with the large upcoming galaxy surveys such as DESI, Euclid, and the Vera Rubin telescope surveys.

Revisiting the extreme clustering of z ≈ 4 quasars with large volume cosmological simulations

ABSTRACT Observations from wide-field quasar surveys indicate that the quasar autocorrelation length increases dramatically from z ≈ 2.5 to ≈ 4. This large clustering amplitude at z ≈ 4 has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter haloes residing in the most extreme environments at that redshift. In this work, we present a model that simultaneously reproduces both the observed quasar autocorrelation and quasar luminosity functions. The spatial distribution of haloes and their relative abundance are obtained via a novel method that computes the halo mass and halo cross-correlation functions by combining multiple large-volume dark-matter-only cosmological simulations with different box sizes and resolutions. Armed with these halo properties, our model exploits the conditional luminosity function framework to describe the stochastic relationship between quasar luminosity, L, and halo mass, M. Assuming a simple power-law relation L ∝ Mγ with lognormal scatter, σ, we are able to reproduce observations at z ∼ 4 and find that: (i) the quasar luminosity–halo mass relation is highly non-linear (γ ≳ 2), with very little scatter (σ ≲ 0.3 dex); (ii) luminous quasars ($\log _{10} L/{\rm erg}\, {\rm s}^{-1}\gtrsim 46.5-47$) are hosted by haloes with mass log10M/M⊙ ≳ 13–13.5; and (iii) the implied duty cycle for quasar activity approaches unity ($\varepsilon _{\rm DC}\approx 10\,\mathrm{ per}\,\mathrm{ cent}-60~{{\ \rm per\ cent}}$). We also consider observations at z ≈ 2.5 and find that the quasar luminosity–halo mass relation evolves significantly with cosmic time, implying a rapid change in quasar host halo masses and duty cycles, which in turn suggests concurrent evolution in black hole scaling relations and/or accretion efficiency.

Accurate halo mass functions from the simplest excursion set theory

ABSTRACT Excursion set theory is a powerful and widely used tool for describing the distribution of dark matter haloes, but it is normally applied with simplifying approximations. We use numerical sampling methods to study the mass functions predicted by the theory without approximations. With a spherical top-hat window and a constant δ = 1.5 threshold, the theory accurately predicts mass functions with the M200 mass definition, both unconditional and conditional, in simulations of a range of matter-dominated cosmologies. For Λ cold dark matter at the present epoch, predictions lie between the M200m and M200c mass functions. In contrast, with the same window function, a non-constant threshold based on ellipsoidal collapse predicts uniformly too few haloes. This work indicates a new way to simply and accurately evaluate halo mass functions, clustering bias, and assembly histories for a range of cosmologies. We provide a fitting function that accurately represents the predictions of the theory for a wide range of parameters.

An almost dark galaxy with the mass of the Small Magellanic Cloud

Almost dark galaxies are objects that have eluded detection by traditional surveys such as the Sloan Digital Sky Survey (SDSS). The low surface brightness of these galaxies (μr(0) > 26 mag arcsec−2), and hence their low surface stellar mass density (a few solar masses per pc2 or less), suggest that the energy density released by baryonic feedback mechanisms is inefficient in modifying the distribution of the dark matter halos they inhabit. For this reason, almost dark galaxies are particularly promising for probing the microphysical nature of dark matter. In this paper, we present the serendipitous discovery of Nube, an almost dark galaxy with ⟨μV⟩e ∼ 26.7 mag arcsec−2. The galaxy was identified using deep optical imaging from the IAC Stripe82 Legacy Project. Follow-up observations with the 100 m Green Bank Telescope strongly suggest that the galaxy is at a distance of 107 Mpc. Ultra-deep multi-band observations with the 10.4 m Gran Telescopio Canarias favour an age of ∼10 Gyr and a metallicity of [Fe/H] ∼ −1.1. With a stellar mass of ∼4 × 108 M⊙ and a half-mass radius of Re = 6.9 kpc (corresponding to an effective surface density of ⟨Σ⟩e ∼ 0.9 M⊙ pc−2), Nube is the most massive and extended object of its kind discovered so far. The galaxy is ten times fainter and has an effective radius three times larger than typical ultradiffuse galaxies with similar stellar masses. Galaxies with comparable effective surface brightness within the Local Group have very low mass (tens of 105 M⊙) and compact structures (effective radius Re < 1 kpc). Current cosmological simulations within the cold dark matter scenario, including baryonic feedback, do not reproduce the structural properties of Nube. However, its highly extended and flattened structure is consistent with a scenario where the dark matter particles are ultralight axions with a mass of mB = (0.8−0.2+0.4) × 10−23 eV.

A Comparison of Void-finding Algorithms Using Crossing Numbers

We study how well void-finding algorithms identify cosmic void regions and whether we can quantitatively and qualitatively compare the voids they find with dynamical information from the underlying matter distribution. Using the ORIGAMI algorithm to determine the number of dimensions along which dark matter particles have undergone shell crossing (crossing number) in N-body simulations from the AbacusSummit simulation suite, we identify dark matter particles that have undergone no shell crossing as belonging to voids. We then find voids in the corresponding halo distribution using two different void-finding algorithms: VoidFinder and V2, a ZOBOV-based algorithm. The resulting void catalogs are compared to the distribution of dark matter particles to examine how their crossing numbers depend on void proximity. While both algorithms’ voids have a similar distribution of crossing numbers near their centers, we find that beyond 0.25 times the effective void radius, voids found by VoidFinder exhibit a stronger preference for particles with low crossing numbers than those found by V2. We examine two possible methods of mitigating this difference in efficacy between the algorithms. While we are able to partially mitigate the ineffectiveness of V2 by using the distance from the void edge as a measure of centrality, we conclude that VoidFinder more reliably identifies dynamically distinct regions of low crossing number.

A physical and concise halo model based on the depletion radius

ABSTRACT We develop a self-consistent and accurate halo model by partitioning matter according to the depletion radii of haloes. Unlike conventional models that define haloes with the virial radius while relying on a separate exclusion radius or ad hoc fixes to account for halo exclusion, our model distributes mass across all scales self-consistently and accounts for both the virialized and non-virialized matter distribution around each halo. Using a cosmological simulation, we show that our halo definition leads to very simple and intuitive model components, with the one-halo term given by the Einasto profile with no truncation needed, and the halo–halo correlation function following a universal power-law form down to the halo boundary. The universal halo–halo correlation also allows us to easily model the distribution of unresolved haloes as well as diffuse matter. Convolving the halo profile with the halo–halo correlation function, we obtain a complete description of the halo–matter correlation across all scales, which self-consistently accounts for halo exclusion at the transition scale. Mass conservation is explicitly maintained in our model, and the scale dependence of the classical halo bias is easily reproduced. Our model can successfully reconstruct the halo–matter correlation function within an accuracy of 9 per cent for halo virial masses in the range of 1011.5h−1 M⊙ < Mvir < 1015.35h−1 M⊙ at z = 0, and covers the radial range of 0.01 h−1 Mpc < r < 20 h−1 Mpc. We also show that our model profile can accurately predict the characteristic depletion radius at the minimum bias and the splash-back radius at the steepest density slope locations.

Mapping the milky way’s stellar halo with 2D data

ABSTRACT We propose a new method for measuring the spatial density distribution of the stellar halo of the Milky Way. Our method is based on a pairwise statistic of the distribution of stars in the sky, the angular two-point correlation function (ATPCF). The ATPCF utilizes two-dimensional data of stars only and is therefore immune to the large uncertainties in the determination of distances to stars. We test our method using mock stellar data coming from various models including the single power-law (SPL) and the broken power-law (BPL) density profiles. We also test the influence of axisymmetric flattening factors using both constant and varying values. We find that the ATPCF is a powerful tool for recovering the spatial distributions of the stellar haloes in our models. We apply our method to observational data from the type ab RR Lyrae catalogue in the Catalina Survey Data Release 1. In the 3-parameter BPL model, we find that $s_{1}=2.46_{-0.20}^{+0.18}, s_{2}=3.99_{-1.33}^{+0.75}$, and $r_{0}=31.11_{-5.88}^{+7.61}$, which are in good agreement with previous results. We also find that introducing an extra parameter, the radially varying flattening factor, greatly improves our ability to model accurately the observed data distribution. This implies perhaps that the stellar halo of the Milky Way should be regarded as oblate.

The MillenniumTNG Project: high-precision predictions for matter clustering and halo statistics

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

The MillenniumTNG Project: the large-scale clustering of galaxies

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

Topological data analysis reveals differences between simulated galaxies and dark matter haloes

ABSTRACT We use topological summaries based on Betti curves to characterize the large-scale spatial distribution of simulated dark matter haloes and galaxies. Using the IllustrisTNG and CAMELS-SAM simulations, we show that the topology of the galaxy distribution is significantly different from the topology of the dark matter halo distribution. Further, there are significant differences between the distributions of star-forming and quiescent galaxies. These topological differences are broadly consistent across all simulations, while at the same time there are noticeable differences when comparing between different models. Finally, using the CAMELS-SAM simulations, we show that the topology of the quiescent galaxies, in particular, depends strongly on the amount of supernova feedback. These results suggest that topological summary statistics could be used to help better understand the processes of galaxy formation and evolution.