Halo Mass Function Research Articles

Euclid preparation. Simulations and non-linearities beyond Lambda cold dark matter. I. Numerical methods and validation

To constrain cosmological models beyond Lambda CDM, the development of the analysis pipeline requires simulations that capture the non-linear phenomenology of such models. We present an overview of numerical methods and $N$-body simulation codes developed to study the non-linear regime of structure formation in alternative dark energy and modified gravity theories. We review a variety of numerical techniques and approximations employed in cosmological $N$-body simulations to model the complex phenomenology of scenarios beyond Lambda CDM. This includes discussions on solving non-linear field equations, accounting for fifth forces, and implementing screening mechanisms. Furthermore, we conduct a code comparison exercise to assess the reliability and convergence of different simulation codes across a range of models. Our analysis demonstrates a high degree of agreement among the outputs of different simulation codes, typically within 2<!PCT!> for the predicted modification of the matter power spectrum and within 4<!PCT!> for the predicted modification of the halo mass function, although some approximations degrade accuracy a bit further. This provides confidence in current numerical methods of modelling cosmic structure formation beyond Lambda CDM. We highlight recent advances made in simulating the non-linear scales of structure formation, which are essential for leveraging the full scientific potential of the forthcoming observational data from the mission.

DiffSph: a Python tool to compute diffuse signals from dwarf spheroidal galaxies

So far no diffuse emissions in dwarf spheroidal satellites of the Milky Way have ever been observed. Given that dwarf galaxies are predominantly composed of Dark Matter, the discovery of these signals could offer valuable insights into understanding the nature of Dark Matter. We present “diffSph”, a Python tool which in its present version provides fast predictions of such diffuse signals in radio frequencies. It also features a very comprehensive module for the computation of “J” and “D” factors that are relevant for indirect Dark Matter detection using gamma rays. Routines are coupled to parton-shower algorithms and Dark Matter halo mass functions from state-of-the-art kinematic fits. This code is also useful for testing generic hypotheses (not necessarily associated with any Dark Matter candidate) about the cosmic-ray electron/positron sources in the dwarf galaxies. The diffSph tool has already been employed in searches for diffuse signals from dwarf spheroidal galaxies using the LOw Frequency ARray (LOFAR).

The importance of stochasticity in determining galaxy emissivities and UV LFs during cosmic dawn and reionization

The stochastic nature of star formation and photon propagation in high-redshift galaxies can result in sizable galaxy-to-galaxy scatter in their properties. Ignoring this scatter by assuming mean quantities can bias estimates of their emissivity and corresponding observables. We constructed a flexible, semi-empirical model, sampling scatter around the following mean relations: (i) the conditional halo mass function (CHMF); (ii) the stellar-to-halo mass relation (SHMR); (iii) the galaxy star formation main sequence (SFMS); (iv) the fundamental metallicity relation (FMR); (v) the conditional intrinsic luminosity; and (vi) the photon escape fraction. In our fiducial model, ignoring scatter in these galaxy properties overestimates the duration of the Epoch of Reionization (EoR), delaying its completion by $ z 1--2. We quantified the relative importance of each of the above sources of scatter in determining the ionizing, soft-band X-ray, and Lyman Werner (LW) emissivities as a function of scale and redshift. We find that scatter around the SFMS is important for all bands, especially at the highest redshifts where the emissivity is dominated by the faintest, most "bursty" galaxies. Ignoring this scatter would underestimate the mean emissivity and its standard deviation computed over 5 cMpc regions by factors of up to sim 2--10 at $5 z 15$. The scatter around the X-ray luminosity to star formation rate and metallicity relation is important for determining X-ray emissivity, accounting for roughly half of its mean and standard deviation. The importance of scatter in the ionizing escape fraction depends on its functional form, while scatter around the SHMR contributes at the level of sim 10--20<!PCT!>. Other sources of scatter have a negligible contribution to the emissivities. Although scatter does flatten the UV luminosity functions, shifting the bright end by 1--2 magnitudes, the level of scatter in our fiducial model is insufficient to fully explain recent estimates from JWST photometry (consistent with previous studies). We conclude that models of the EoR should account for the burstiness of star formation, while models for the cosmic 21cm signal should additionally account for scatter in intrinsic X-ray production.

The baryon census and the mass-density of stars, neutral gas, and hot gas as a function of halo mass

Abstract We study the stellar, neutral gas content within halos over a halo mass range 1010to1015.5M⊙ and hot X-ray gas content over a halo mass range 1012.8to1015.5M⊙ in the local universe. We combine various empirical datasets of stellar, H i and X-ray observations of galaxies, groups and clusters to establish fundamental baryonic mass vs halo mass scaling relations. These scaling relations are combined with halo mass function to obtain the baryon densities of stars, neutral gas and hot gas (T > 106K), as a function of halo mass. We calculate the contributions of the individual baryonic components to the cosmic baryon fraction. Cosmic stellar mass density ($\Omega _\text{star}=2.09^{+0.21}_{-0.18} \times 10^{-3}$), cosmic H i mass density ($\Omega _\rm{H\,{\small I}}=0.49^{+0.25}_{-0.12} \times 10^{-3}$) and cosmic neutral gas mass density ($\Omega _\text{neutral gas}=0.71^{+0.39}_{-0.18} \times 10^{-3}$) estimates are consistent with previous more direct method measurements of these values, thereby establishing the veracity of our method. We also give an estimate of the cosmic hot plasma density ($\Omega _\text{hot gas}=2.58^{+2.1}_{-0.66} \times 10^{-3}$).

The VST ATLAS Quasar Survey – III. Halo mass function via quasar clustering and quasar-CMB lensing cross-clustering

ABSTRACT We exploit the VST ATLAS quasar (QSO) catalogue to perform three measurements of the quasar halo mass profile. First, we make a new estimate of the angular autocorrelation function of ≈230 000 ATLAS quasars with $z_{\rm photo}\lesssim 2.5$ and $17 < g < 22$. By comparing with the $\Lambda$CDM mass clustering correlation function, we measure the quasar bias to be $b_{\rm Q}\approx 2.1$, implying a quasar halo mass of $M_{\rm halo} \approx 8.5\times 10^{11}\,h^{-1}\,{\rm M}_\odot$. Second, we cross-correlate these $z\approx 1.7$ ATLAS quasars with the Planck cosmic microwave background (CMB) lensing maps, detecting a somewhat stronger signal at $4\,{\rm arcmin} < \theta < 60\,{\rm arcmin}$ than previous authors. Scaling these authors’ model fit to our data, we estimate a quasar host halo mass of $M_{\rm halo}\approx 8.3\times 10^{11}\,h^{-1}\,{\rm M}_\odot$. Third, we fit halo occupation sistribution (HOD) model parameters to our quasar autocorrelation function and from the derived halo mass function, we estimate a quasar halo mass of $M_{\rm halo}\approx 2.5\times 10^{12}\,h^{-1}\,{\rm M}_\odot$. We then compare our HOD model prediction to our quasar-CMB lensing result, confirming their consistency. We find that most (≈2/3) QSOs have halo masses within a factor of ≈3 of this average mass. An analysis based on the probability of X-ray detections of AGN in galaxies and the galaxy stellar mass function gives a similarly small mass range. Finally, we compare the quasar halo mass and luminosity functions and suggest that gravitational growth may produce the constant space density with redshift seen in the quasar luminosity function.

JWST lensed quasar dark matter survey – II. Strongest gravitational lensing limit on the dark matter free streaming length to date

ABSTRACT This is the second in a series of papers in which we use JWST Mid Infrared Instrument multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of nine systems. The warm dust region is compact and sensitive to perturbations by populations of haloes down to masses $\sim 10^6$ M$_{\odot }$. Using these warm dust flux-ratio measurements in combination with five previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third- and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhaloes. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6}\, {\rm M}_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly $\alpha$ forest and Milky Way satellite galaxies.

Quijote-PNG: Optimizing the Summary Statistics to Measure Primordial Non-Gaussianity

We apply a suite of different estimators to the Quijote-png halo catalogs to find the best approach to constrain Primordial non-Gaussianity (PNG) at nonlinear cosmological scales, up to kmax=0.5hMpc−1 . The set of summary statistics considered in our analysis includes the power spectrum, bispectrum, halo mass function, marked power spectrum, and marked modal bispectrum. Marked statistics are used here for the first time in the context of the PNG study. We perform a Fisher analysis to estimate their cosmological information content, showing substantial improvements when marked observables are added to the analysis. Starting from these summaries, we train deep neural networks to perform likelihood-free inference of cosmological and PNG parameters. We assess the performance of different subsets of summary statistics; in the case of fNLequil , we find that a combination of the power spectrum and a suitable marked power spectrum outperforms the combination of power spectrum and bispectrum, the baseline statistics usually employed in PNG analysis. A minimal pipeline to analyze the statistics we identified can be implemented either with our ML algorithm or via more traditional estimators, if these are deemed more reliable.

Cosmological simulations of scale-dependent primordial non-Gaussianity

We present the results of a set of cosmological N-body simulations with standard ΛCDM cosmology but characterized by a scale-dependent primordial non-Gaussianity of the local type featuring a power-law dependence of the f NL loc(k) at large scales followed by a saturation to a constant value at smaller scales where non-linear growth leads to the formation of collapsed cosmic structures. Such models are built to ensure consistency with current Cosmic Microwave Background bounds on primordial non-Gaussianity yet allowing for large effects of the non-Gaussian statistics on the properties of non-linear structure formation. We show the impact of such scale-dependent non-Gaussian scenarios on a wide range of properties of the resulting cosmic structures, such as the non-linear matter power spectrum, the halo and sub-halo mass functions, the concentration-mass relation, the halo and void density profiles, and we highlight for the first time that some of these models might mimic the effects of Warm Dark Matter for several of such observables.

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.

Euclid preparation

The Euclid mission, designed to map the geometry of the dark Universe, presents an unprecedented opportunity for advancing our understanding of the cosmos through its photometric galaxy cluster survey. Central to this endeavor is the accurate calibration of the mass- and redshift-dependent halo bias (HB), which is the focus of this paper. Our aim is to enhance the precision of HB predictions, which is crucial for deriving cosmological constraints from the clustering of galaxy clusters. Our study is based on the peak-background split (PBS) model linked to the halo mass function (HMF), and it extends it with a parametric correction to precisely align with results from an extended set of N-body simulations carried out with the OpenGADGET3 code. Employing simulations with fixed and paired initial conditions, we meticulously analyzed the matter-halo cross-spectrum and modeled its covariance using a large number of mock catalogs generated with Lagrangian perturbation theory simulations with the PINOCCHIO code. This ensures a comprehensive understanding of the uncertainties in our HB calibration. Our findings indicate that the calibrated HB model is remarkably resilient against changes in cosmological parameters, including those involving massive neutrinos. The robustness and adaptability of our calibrated HB model provide an important contribution to the cosmological exploitation of the cluster surveys to be provided by the Euclid mission. This study highlights the necessity of continuously refining the calibration of cosmological tools such as the HB to match the advancing quality of observational data. As we project the impact of our calibrated model on cosmological constraints, we find that given the sensitivity of the Euclid survey, a miscalibration of the HB could introduce biases in cluster cosmology analysis. Our work fills this critical gap, ensuring the HB calibration matches the expected precision of the Euclid survey.

Cosmic stability of dark matter from Pauli blocking

Why does dark matter (DM) live longer than the age of the Universe? Here we study a novel sub-eV scalar DM candidate whose stability is due to the Pauli exclusion of its fermionic decay products. We analyze the stability of the DM condensate against decays, scatterings (i.e., evaporation), and parametric resonance, delineating the viable parameter regions in which DM is cosmologically stable. In a minimal scenario in which the scalar DM decays to a pair of new exotic fermions, we find that scattering can populate an interacting thermal dark sector component to energies far above the DM mass. This self-interacting dark radiation may potentially alleviate the Hubble tensions. Furthermore, our scenario can be probed through precise measurements of the halo-mass function or the masses of dwarf spheroidal galaxies since scattering prevents the DM from becoming too dense. On the other hand, if the lightest neutrino stabilizes the DM, the cosmic neutrino background (CνB) can be significantly altered from the Lambda cold dark matter prediction and thus be probed in the future by CνB detection experiments. Published by the American Physical Society 2024

Bayesian inference of initial conditions from non-linear cosmic structures using field-level emulators

ABSTRACT Analysing next-generation cosmological data requires balancing accurate modelling of non-linear gravitational structure formation and computational demands. We propose a solution by introducing a machine learning-based field-level emulator, within the Hamiltonian Monte Carlo-based Bayesian Origin Reconstruction from Galaxies (BORG) inference algorithm. Built on a V-net neural network architecture, the emulator enhances the predictions by first-order Lagrangian perturbation theory to be accurately aligned with full N-body simulations while significantly reducing evaluation time. We test its incorporation in BORG for sampling cosmic initial conditions using mock data based on non-linear large-scale structures from N-body simulations and Gaussian noise. The method efficiently and accurately explores the high-dimensional parameter space of initial conditions, fully extracting the cross-correlation information of the data field binned at a resolution of $1.95\,h^{-1}$ Mpc. Percent-level agreement with the ground truth in the power spectrum and bispectrum is achieved up to the Nyquist frequency $k_\mathrm{N} \approx 2.79h \,\, \mathrm{Mpc}^{-1}$. Posterior resimulations – using the inferred initial conditions for N-body simulations – show that the recovery of information in the initial conditions is sufficient to accurately reproduce halo properties. In particular, we show highly accurate $M_{200\mathrm{c}}$ halo mass function and stacked density profiles of haloes in different mass bins $[0.853,16]\times 10^{14}\,{\rm M}_{\odot }\,h^{-1}$. As all available cross-correlation information is extracted, we acknowledge that limitations in recovering the initial conditions stem from the noise level and data grid resolution. This is promising as it underscores the significance of accurate non-linear modelling, indicating the potential for extracting additional information at smaller scales.

The e-MANTIS emulator: Fast and accurate predictions of the halo mass function in f(R)CDM and wCDM cosmologies

In this work, we present a novel emulator of the halo mass function (HMF), which we implemented in the framework of the e-MANTIS emulator of $f(R)$ gravity models. We also extended e-MANTIS to cover a larger cosmological parameter space and to include models of dark energy with a constant equation of state $w$CDM. We used a Latin hypercube sampling of the $w$CDM and $f(R)$CDM cosmological parameter spaces, over a wide range, and carried out a large suite of more than $10000$ $N$-body simulations with a different volume, mass resolution, and random phase for the initial conditions. For each simulation in the suite, we generated halo catalogues using the friends-of-friends (FoF) halo finder, as well as the spherical overdensity (SO) algorithm for different overdensity thresholds (200, 500, and 1000 times the critical density). We decomposed the corresponding HMFs on a B-spline basis, while adopting a minimal set of assumptions on their shape. We used this decomposition to train an emulator based on Gaussian processes. The resulting emulator is able to predict the HMF for redshifts $ 1.5$ and for halo masses $M_h M_ The typical HMF errors for SO haloes with $ c $ at $z=0$ in $w$CDM (respectively $f(R)$CDM) are of the order of $ up to a transition mass $M_t M_ ($M_t M_ For larger masses, the errors are dominated by the shot noise and scale as $ with $ up to $M_h M_ Independently of this general trend, the emulator is able to provide an estimation of its own error as a function of the cosmological parameters, halo mass, and redshift. We have performed an extensive comparison against analytical parametrizations and shown that e-MANTIS is able to better capture the cosmological dependence of the HMF, while being complementary to other existing emulators. The e-MANTIS emulator, which is publicly available, can be used to obtain fast and accurate predictions of the HMF in the $f(R)$CDM and $w$CDM non-standard cosmological models. As such, it represents a useful theoretical tool to constrain the nature of dark energy using data from galaxy cluster surveys.

The JWST EXCELS survey: too much, too young, too fast? Ultra-massive quiescent galaxies at 3 < z < 5

ABSTRACT We report ultra-deep, medium-resolution spectroscopic observations for four quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot })\gt 11$ at $3 \lt z \lt 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first two galaxies are newly selected from PRIMER UDS imaging, both at $z=4.62$ and separated by 860 pkpc on the sky, within a larger structure for which we confirm several other members. Both formed at $z\simeq 8-10$. These systems could plausibly merge by the present day to produce a local massive elliptical galaxy. The other two ultra-massive quiescent galaxies are previously known at $z=3.99$ and 3.19, with the latter (ZF-UDS-7329) having been the subject of debate as potentially too old and too massive to be accommodated by the $\Lambda$-CDM halo-mass function. Both exhibit high stellar metallicities, and for ZF-UDS-7329 we are able to measure the $\alpha -$enhancement, obtaining [Mg/Fe] = $0.42^{+0.19}_{-0.17}$. We finally evaluate whether these four galaxies are consistent with the $\Lambda$-CDM halo-mass function using an extreme value statistics approach. We find that the $z=4.62$ objects and the $z=3.19$ object are unlikely within our area under the assumption of standard stellar fractions ($f_*\simeq 0.1-0.2$). However, these objects roughly align with the most massive galaxies expected under the assumption of 100 per cent conversion of baryons to stars ($f_*$=1). Our results suggest extreme galaxy formation physics during the first billion years, but no conflict with $\Lambda$-CDM cosmology.

The circular velocity and halo mass functions of galaxies in the nearby Universe

ABSTRACT The circular velocity function (CVF) of galaxies is a fundamental test of the Lambda cold dark matter ($\Lambda$CDM) paradigm as it traces the variation of galaxy number densities with circular velocity ($v_{\rm {circ}}$), a proxy for dynamical mass. Previous observational studies of the CVF have either been based on H i-rich galaxies, or encompassed low-number statistics and probed narrow ranges in $v_{\rm {circ}}$. We present a benchmark computation of the CVF between $100\,{\text{and}}\,350\ \rm {km\ s^{-1}}$ using a sample of 3527 nearby Universe galaxies, representative for stellar masses between $10^{9.2}\,{\text{and}}\,10^{11.9} \rm {{\rm M}_{\odot }}$. We find significantly larger number densities above 150 $\rm {km\ s^{-1}}$ compared to results from H i surveys, pertaining to the morphological diversity of our sample. Leveraging the fact that circular velocities are tracing the gravitational potential of haloes, we compute the halo mass function (HMF), covering $\sim$1 dex of previously unprobed halo masses ($10^{11.7}{\!-\!}10^{12.7} \rm {{\rm M}_{\odot }}$). The HMF for our sample, representative of the galaxy population with $M_{200}\geqslant 10^{11.35} \rm {{\rm M}_{\odot }}$, shows that spiral morphologies contribute 67 per cent of the matter density in the nearby Universe, while early types account for the rest. We combine our HMF data with literature measurements based on H i kinematics and group/cluster velocity dispersions. We constrain the functional form of the HMF between $10^{10.5}-10^{15.5} \rm {{\rm M}_{\odot }}$, finding a good agreement with $\Lambda$CDM predictions. The halo mass range probed encompasses 72$\substack{+5 -6}$ per cent ($\Omega _{\rm {M,10.5-15.5}} = 0.227 \pm 0.018$) of the matter density in the nearby Universe; 31$\substack{+5 -6}$ per cent is accounted for by haloes below $10^{12.7}\rm {{\rm M}_{\odot }}$ occupied by a single galaxy.

Not-so-little Red Dots: Two Massive and Dusty Starbursts at z ∼ 5–7 Pushing the Limits of Star Formation Discovered by JWST in the COSMOS-Web Survey

We present the properties of two candidate massive (M ⋆ ∼ 1011 M ⊙) and dusty (A v > 2.5 mag) galaxies at z = 5–7 in the first 0.28 deg2 of the COSMOS-Web survey. One object is spectroscopically confirmed at z spec = 5.051, while the other has a robust z phot = 6.7 ± 0.3. Thanks to their extremely red colors (F277W–F444W ∼ 1.7 mag), these galaxies satisfy the nominal color selection for the widely studied “little red dot” (LRD) population with the exception of their spatially resolved morphologies. The morphology of our targets allows us to conclude that their red continuum is dominated by highly obscured stellar emission and not by reddened nuclear activity. Using a variety of spectral energy distribution fitting tools and star formation histories, we estimate the stellar masses to be log(M⋆)=11.32−0.15+0.07M⊙ and log(M⋆)=11.2−0.2+0.1M⊙ , respectively, with a red continuum emission dominated by a recent episode of star formation. We then compare their number density to the halo mass function to infer stellar baryon fractions of ϵ ⋆ ∼ 0.25 and ϵ ⋆ ∼ 0.5. Both are significantly higher than what is commonly observed in lower-z galaxies or more dust-obscured galaxies at similar redshifts. With very bright ultra-high-z Lyman-Break Galaxies and some non-AGN-dominated LRDs, such “extended” LRDs represent another population that may require very efficient star formation at early times.

JWST PRIMER: a new multifield determination of the evolving galaxy UV luminosity function at redshifts z ≃ 9 – 15

ABSTRACT We present a new determination of the evolving galaxy ultraviolet (UV) luminosity function (LF) over the redshift range $8.5< z< 15.5$ using a combination of several major Cycle-1 JWST imaging programmes – Public Release IMaging for Extragalactic Research, JWST Advanced Deep Extragalactic Survey, and Next Generation Deep Extragalactic Exploratory Public Survey. This multifield approach yields a total of $\simeq 370$ arcmin2 of JWST/NIRCam imaging, reaching (5-$\sigma$) depths of $\simeq 30$ AB mag in the deepest regions. We select a sample of 2548 galaxies with a significant probability of lying at high redshift ($p(z> 8.5)> 0.05$) to undertake a statistical calculation of the UV LF. Our new measurements span $\simeq 4$ mag in UV luminosity at $z=9-12.5$, placing new constraints on both the shape and evolution of the LF at early times. Our measurements yield a new estimate of the early evolution of cosmic star-formation rate density ($\rho _{\rm {SFR}}$) confirming the gradual decline deduced from early JWST studies, at least out to $z \simeq 12$. Finally we show that the observed early evolution of the galaxy UV LF (and $\rho _{\rm {SFR}}$) can be reproduced in a ${\rm \Lambda }$cold dark matter Universe, with no change in dust properties or star-formation efficiency required out to $z \simeq 12$. Instead, a progressive trend towards younger stellar population ages can reproduce the observations, and the typical ages required at $z \simeq$ 8, 9, 10, and 11 all converge on $\simeq 380-330$ Myr after the big bang, indicative of a rapid emergence of early galaxies at $z \simeq 12 - 13$. This is consistent with the first indications of a steeper drop-off in $\rho _{\rm {SFR}}$ we find beyond $z \simeq 13$, possibly reflecting the rapid evolution of the halo mass function at earlier times.

Exploring the role of the halo mass function for inferring astrophysical parameters during reionisation

Abstract Detecting the 21-cm signal at z ≳ 6 will reveal insights into the properties of the first galaxies responsible for driving reionisation. To extract this information, we perform parameter inference with 3D simulations of the 21-cm signal embedded within a Bayesian inference pipeline. Presently, when performing inference, we must choose which sources of uncertainty to sample and which to hold fixed. Since the astrophysics of galaxies is much more uncertain than that of the underlying halo-mass function (HMF), we typically parameterise and model the former while fixing the latter. However, doing so may bias our inference of the galaxy properties. In this work, we explore the consequences of assuming an incorrect HMF and quantify the relative biases on our inferred astrophysical model parameters when considering the wrong HMF. We then relax this assumption by constructing a generalised five parameter HMF model and simultaneously recover it with our underlying astrophysical model. For this, we use 21CMFAST and perform Simulation-Based Inference using marginal neural ratio estimation to learn the likelihood-to-evidence ratio with Swyft. Using a mock 1000-hour observation of the 21-cm power spectrum from the forthcoming Square Kilometre Array, conservatively assuming foreground wedge avoidance, we find that assuming the incorrect HMF can bias the recovered astrophysical parameters by up to ∼3 − 4σ even when including independent information from observed luminosity functions. Using our generalised HMF model, although we recover our astrophysical parameters with a factor of ∼2 − 4 larger marginalised uncertainties, the constraints are unbiased, agnostic to the underlying HMF and therefore more conservative.

Measuring the Conditional Luminosity and Stellar Mass Functions of Galaxies by Combining the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys Data Release 9, Survey Validation 3, and Year 1 Data

In this investigation, we leverage the combination of the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 9, Survey Validation 3, and Year 1 data sets to estimate the conditional luminosity functions and conditional stellar mass functions (CLFs and CSMFs) of galaxies across various halo mass bins and redshift ranges. To support our analysis, we utilize a realistic DESI mock galaxy redshift survey (MGRS) generated from a high-resolution Jiutian simulation. An extended halo-based group finder is applied to both MGRS catalogs and DESI observation. By comparing the r- and z-band luminosity functions (LFs) and stellar mass functions (SMFs) derived using both photometric and spectroscopic data, we quantified the impact of photometric redshift (photo-z) errors on the galaxy LFs and SMFs, especially in the low-redshift bin at the low-luminosity/mass end. By conducting prior evaluations of the group finder using MGRS, we successfully obtain a set of CLF and CSMF measurements from observational data. We find that at low redshift, the faint-end slopes of CLFs and CSMFs below ∼109 h −2 L ⊙ (or h −2 M ⊙) evince a compelling concordance with the subhalo mass functions. After correcting the cosmic variance effect of our local Universe following Chen et al., the faint-end slopes of the LFs/SMFs turn out to also be in good agreement with the slope of the halo mass function.

The hot circumgalactic medium in the eROSITA All-Sky Survey II. Scaling relations between X-ray luminosity and galaxies' mass

Understanding how the properties of galaxies relate to the properties of the hot circum-galactic medium (CGM) around them can constrain galaxy evolution models. We aim to measure the scaling relations between the X-ray luminosity of the hot CGM and the fundamental properties (stellar mass and halo mass) of a galaxy. We measured the X-ray luminosity of the hot CGM based on the surface brightness profiles of central galaxy samples measured from Spectrum Roentgen Gamma (SRG)/eROSITA all-sky survey data. We related the X-ray luminosity to the galaxies' stellar and halo mass, and we compared the observed relations to the self-similar model and intrinsic (i.e., not forward-modeled) output of the IllustrisTNG, EAGLE, and SIMBA simulations. The average hot CGM X-ray luminosity ($L_ X,CGM $) correlates with the galaxy's stellar mass ($M_*$). It increases from $(1.6 $ to $(3.4 $, when $ increases from 10.0 to 11.5. A power law describes the correlation as $ X,CGM The hot CGM X-ray luminosity as a function of halo mass is measured within $ 500c )=11.3-13.7$, extending our knowledge of the scaling relation by more than two orders of magnitude. $L_ X,CGM $ increases with $M_ 500c $ from $(3.0 $ at $ 500c )=11.3$ to $(1.3 $ at $ 500c )=13.7$. The relation follows a power law of $ X,CGM 500c Our observations highlight the necessity of non-gravitational processes at the galaxy group scale while suggesting these processes are sub-dominant at the galaxy scale. We show that the outputs of current cosmological galaxy simulations generally align with the observational results uncovered here but with possibly important deviations in selected mass ranges. We explore, at the low mass end, the average scaling relations between the CGM X-ray luminosity and the galaxy's stellar mass or halo mass, which constitutes a new benchmark for galaxy evolution models and feedback processes.