Halo Mass Function Research Articles

Earliest Galaxies in the JADES Origins Field: Luminosity Function and Cosmic Star Formation Rate Density 300 Myr after the Big Bang

We characterize the earliest galaxy population in the JADES Origins Field, the deepest imaging field observed with JWST. We make use of ancillary Hubble Space Telescope optical images (five filters spanning 0.4–0.9 μm) and novel JWST images with 14 filters spanning 0.8−5 μm, including seven medium-band filters, and reaching total exposure times of up to 46 hr per filter. We combine all our data at >2.3 μm to construct an ultradeep image, reaching as deep as ≈31.4 AB mag in the stack and 30.3–31.0 AB mag (5σ, r = 0.″1 circular aperture) in individual filters. We measure photometric redshifts and use robust selection criteria to identify a sample of eight galaxy candidates at redshifts z = 11.5−15. These objects show compact half-light radii of R 1/2 ∼ 50−200 pc, stellar masses of M ⋆ ∼ 107−108 M ☉, and star formation rates ∼ 0.1−1 M ☉ yr−1. Our search finds no candidates at 15 < z < 20, placing upper limits at these redshifts. We develop a forward-modeling approach to infer the properties of the evolving luminosity function without binning in redshift or luminosity that marginalizes over the photometric redshift uncertainty of our candidate galaxies and incorporates the impact of nondetections. We find a z = 12 luminosity function in good agreement with prior results, and that the luminosity function normalization and UV luminosity density decline by a factor of ∼2.5 from z = 12 to z = 14. We discuss the possible implications of our results in the context of theoretical models for evolution of the dark matter halo mass function.

Primordial magnetic fields: consistent initial conditions and impact on high-z structures

Primordial magnetic fields (PMFs) can enhance matter power spectrum on small scales (≲ Mpc) and still agree with bounds from cosmic microwave background (CMB) and Faraday rotation measurements. As modes on scales smaller than Mpc have already become non-linear today, exploring PMFs' impact on small-scale structures requires dedicated cosmological simulations. Here, for the first time, we perform a suite of hydrodynamical simulations that take into account the different impacts of PMFs on baryons and dark matter. Specifically, in the initial conditions we displace particles according to the Lorentz force from PMFs. We also highlight the large theoretical uncertainty in the peak enhancement of the matter power spectrum due to PMFs, which was not considered in previous studies. We present halo mass functions and show that they can be accurately reproduced using Sheth-Tormen formalism. Moreover, we show that PMFs can generate galaxies with baryon fraction several times larger than the cosmic average at high redshifts. This is simply a consequence of the fact that PMFs enhance baryon perturbations, causing them to be larger than dark matter perturbations. We argue that this scenario could be tested soon by obtaining accurate estimates of the baryon fraction in high redshift galaxies.

Constraining f(R) gravity using future galaxy cluster abundance and weak-lensing mass calibration datasets

We present forecasts for constraints on the Hu and Sawicki f(R) modified gravity model using realistic mock data representative of future cluster and weak lensing surveys. We create mock thermal Sunyaev-Zel’dovich effect selected cluster samples for SPT-3G and CMB-S4 and the corresponding weak gravitational lensing data from next-generation weak-lensing (ngWL) surveys like Euclid and Rubin. We employ a state-of-the-art Bayesian likelihood approach that includes all observational effects and systematic uncertainties to obtain constraints on the f(R) gravity parameter log10|fR0|. In this analysis we vary the cosmological parameters [Ωm,Ωνh2,h,As,ns,log10|fR0|], which allows us to account for possible degeneracies between cosmological parameters and f(R) modified gravity. The analysis accounts for f(R) gravity via its effect on the halo mass function which is enhanced on cluster mass scales compared to the expectations within general relativity (GR). Assuming a fiducial GR model, the upcoming cluster dataset SPT−3G×ngWL is expected to obtain an upper limit of log10|fR0|&lt;−5.95 at 95% credibility, which significantly improves upon the current best bounds. The CMB−S4×ngWL dataset is expected to improve this even further to log10|fR0|&lt;−6.23. Furthermore, f(R) gravity models with log10|fR0|≥−6, which have larger numbers of clusters, would be distinguishable from GR with both datasets. We also report degeneracies between log10|fR0| and Ωm as well as σ8 for log10|fR0|&gt;−6 and log10|fR0|&gt;−5 respectively. Our forecasts indicate that future cluster abundance studies of f(R) gravity will enable substantially improved constraints that are competitive with other cosmological probes. Published by the American Physical Society 2024

Cluster cosmology redux: a compact representation for the halo mass function

ABSTRACT Groups and clusters of galaxies imprint coherent, arcminute-scale features across the spectrophotometric sky, especially optical-IR clusters of galaxies, spectral distortions in the cosmic microwave background, and extended sources of X-ray emission. The space–time density of the host dark matter halo population – the halo mass function (HMF) – is a common theoretical basis for modelling such observable features. We explore a compact representation – a dual-quadratic (DQ-HMF) form – that features readily interpretable parameters representing polynomial expansions of the space–time number density surface, first in terms of log-mass, then in redshift. The DQ-HMF form fits Mira-Titan N-body emulator expectations for halo masses $10^{13.7-14.5} \, h^{-1}\, {\rm M}_\odot$ over redshifts, 0.1 &amp;lt; z &amp;lt; 1.5 to within $\sim \! 5~{{\ \rm per\ cent}}$. We provide best-fitting parameters for a Planck 2018 cosmology and demonstrate model self-similarity in the Ωm−$\, \sigma _8$ plane. Convolving with a minimal mass–observable relation (MOR) yields closed-form expressions for counts, mean mass, and mass variance of cluster samples characterized by an observable property. Performing information-matrix forecasts of potential parameter constraints from existing and future surveys, we demonstrate the potential for percent-level constraints on model parameters by an LSST-like optical cluster survey of 300 000 clusters and a richness–mass variance of 0.32. Even better constraints could potentially be achieved by a survey with one-tenth the sample size but with a reduced selection property variance of 0.12. Potential benefits and extensions to the basic MOR parametrization are discussed.

Characterizing ultra-high-redshift dark matter halo demographics and assembly histories with the gureft simulations

ABSTRACT Dark matter halo demographics and assembly histories are a manifestation of cosmological structure formation and have profound implications for the formation and evolution of galaxies. In particular, merger trees provide fundamental input for several modelling techniques, such as semi-analytic models (SAMs), sub-halo abundance matching (SHAM), and decorated halo occupation distribution models. Motivated by the new ultra-high-redshift (z ≳ 10) frontier enabled by JWST, we present a new suite of Gadget at Ultrahigh Redshift with Extra-Fine Timesteps (gureft) dark matter-only cosmological simulations that are carefully designed to capture halo merger histories and structural properties in the ultra-z universe. The simulation suite consists of four 10243-particle simulations with box sizes of 5, 15, 35, and 90 Mpc h−1, each with 170 snapshots stored between 40 ≥ z ≥ 6. With the unprecedented number of available snapshots and strategically chosen dynamic range covered by these boxes, gureft uncovers the emerging dark matter halo populations and their assembly histories in the earliest epochs of cosmic history. In this work, we present the halo mass functions (HMF) between z ∼ 20 and 6 down to log (Mvir/M⊙) ∼ 5, and show that at high redshift, these robust HMFs can differ substantially from commonly used analytic approximations or older fitting functions in the literature. We also present key physical properties of the ultra-high z halo population, such as concentration and spin, as well as their mass growth and merger rates, and again provide updated fitting functions.

Characterizing structure formation through instance segmentation

Dark matter haloes form from small perturbations to the almost homogeneous density field of the early universe. Although it is known how large these initial perturbations must be to form haloes, it is rather poorly understood how to predict which particles will end up belonging to which halo. However, it is this process that determines the Lagrangian shape of proto-haloes and it is therefore essential to understand their mass, spin, and formation history. We present a machine learning framework to learn how the proto-halo regions of different haloes emerge from the initial density field. We developed one neural network to distinguish semantically which particles become part of any halo and a second neural network that groups these particles by halo membership into different instances. This instance segmentation is done through the Weinberger method, in which the network maps particles into a pseudo-space representation where different instances can easily be distinguished through a simple clustering algorithm. Our model reliably predicts the masses and Lagrangian shapes of haloes object by object, as well as other properties such as the halo-mass function. We find that our model extracts information close to optimally by comparing it to the degree of agreement between two N-body simulations with slight differences in their initial conditions. We publish our model open source and suggest that it can be used to inform analytical methods of structure formation by studying the effect of systematic manipulations of the initial conditions.

Primordial Black Hole–Neutron Star Merger Rate in Modified Gravity

In this work, we investigate the merger rate of primordial black hole–neutron star (PBH-NS) binaries in two widely studied modified gravity (MG) models: Hu–Sawicki f(R) gravity and the normal branch of Dvali–Gabadadze–Porrati gravity. In our analysis, we take into account the effects of MG on the halo properties, including halo mass function, halo concentration parameter, halo density profile, and velocity dispersion of dark matter particles. We find that these MG models, due to their stronger gravitational field induced by an effective fifth force, predict enhanced merger rates compared to general relativity. This enhancement is found to be redshift-dependent and sensitive to model parameters and PBH mass and fraction. Assuming a PBH mass range of 5–50 M ⊙, we compare the predicted merger rate of PBH-NS binaries with those inferred from LIGO–Virgo–KAGRA observations of gravitational waves (GWs). We find that the merger rates obtained from MG models will be consistent with the GW observations if the abundance of PBHs is relatively large, with the exact amount depending on the MG model and its parameter values, as well as PBH mass. We also establish upper limits on the abundance of PBHs in these MG frameworks while comparing them with the existing non-GW constraints, which can potentially impose even more stringent constraints.

Cosmological Test of an Ultraviolet Origin of Dark Energy

The accelerated expansion of the Universe is impressively well described by a cosmological constant. However, the observed value of the cosmological constant is much smaller than expected based on quantum field theories. Recent efforts to achieve consistency in these theories have proposed a relationship between Dark Energy and the most compact objects, such as black holes (BHs). However, experimental tests are very challenging to devise and perform. In this article, we present a testable model with no cosmological constant in which the accelerated expansion can be driven by black holes. The model couples the expansion of the Universe (the Friedmann equation) with the mass function of cosmological halos (using the Press–Schechter formalism). Through the observed link between halo masses and BH masses, one thus gets a coupling between the expansion rate of the Universe and the BHs. We compare the predictions of this simple BH model with SN1a data and find poor agreement with observations. Our method is sufficiently general to allow us to also test a fundamentally different model, also without a cosmological constant, where the accelerated expansion is driven by a new force proportional to the internal velocity dispersion of galaxies. Surprisingly enough, this model cannot be excluded using the SN1a data.

Baryonic Imprints on DM Halos: the concentration-mass relation and its dependence on halo and galaxy properties

The halo concentration-mass relation has ubiquitous use in modeling the matter field for cosmological and astrophysical analyses, and including the imprints from galaxy formation physics is tantamount to its robust usage. Many analyses, however, probe the matter around halos selected by a given halo/galaxy property – rather than by halo mass – and the imprints under each selection choice can be different. We employ the CAMELS simulation suite to quantify the astrophysics and cosmology dependence of the concentration-mass relation, , when selected on five properties: (i) velocity dispersion, (ii) formation time, (iii) halo spin, (iv) stellar mass, and (v) gas mass. We construct simulation-informed nonlinear models for all properties as a function of halo mass, redshift, and six cosmological/astrophysical parameters, with a mass range . There are many mass-dependent imprints in all halo properties, with clear connections across different properties and non-linear couplings between the parameters. Finally, we extract the relation for subsamples of halos that have scattered above/below the mean property- relation for a chosen property. Selections on gas mass or stellar mass have a significant impact on the astrophysics/cosmology dependence of , while those on any of the other three properties have a significant (mild) impact on the cosmology (astrophysics) dependence. We show that ignoring such selection effects can lead to errors of ≈25% in baryon imprint modelling of . Our nonlinear model for all properties is made publicly available.

Compatibility of JWST results with exotic halos

The James Webb Space Telescope (JWST) is unveiling astounding results about the first few hundred million years of life of the Universe, delivering images of galaxies at very high redshifts. Here, we develop a UV luminosity function model for high-redshift galaxies, considering parameters such as the stellar formation rate, dust extinction, and halo mass function. Calibration of this luminosity function model using UV luminosity data at redshifts z=4-7 yields optimal parameter values. Testing the model against data at higher redshifts reveals successful accommodation of the data at z=8-9, but challenges emerge at z≃13. Our findings suggest a negligible role of dust extinction at the highest redshifts, prompting a modification of the stellar formation rate to incorporate a larger fraction of luminous objects per massive halo, consistently with similar recent studies. This effect could be attributed to mundane explanations such as unknown evolution of standard astrophysics at high redshift or to the existence of exotic objects at high redshift. We comment on this latter possibility.

The FLAMINGO project: baryonic impact on weak gravitational lensing convergence peak counts

ABSTRACT Weak gravitational lensing convergence peaks, the local maxima in weak lensing convergence maps, have been shown to contain valuable cosmological information complementary to commonly used two-point statistics. To exploit the full power of weak lensing for cosmology, we must model baryonic feedback processes because these reshape the matter distribution on non-linear and mildly non-linear scales. We study the impact of baryonic physics on the number density of weak lensing peaks using the FLAMINGO cosmological hydrodynamical simulation suite. We generate ray-traced full-sky convergence maps mimicking the characteristics of a Stage IV weak lensing survey. We compare the number densities of peaks in simulations that have been calibrated to reproduce the observed galaxy mass function and cluster gas fraction or to match a shifted version of these, and that use either thermally driven or jet active galactic nucleus feedback. We show that the differences induced by realistic baryonic feedback prescriptions (typically 5–30 per cent for κ = 0.1–0.4) are smaller than those induced by reasonable variations in cosmological parameters (20–60 per cent for κ = 0.1–0.4) but must be modelled carefully to obtain unbiased results. The reasons behind these differences can be understood by considering the impact of feedback on halo masses, or by considering the impact of different cosmological parameters on the halo mass function. Our analysis demonstrates that, for the range of models we investigated, the baryonic suppression is insensitive to changes in cosmology up to κ ≈ 0.4 and that the higher κ regime is dominated by Poisson noise and cosmic variance.

Constraining primordial non-Gaussianity from DESI quasar targets and Planck CMB lensing

We detect the cross-correlation between 2.7 million DESI quasar targets across 14,700 deg2 (180 quasars deg-2) and Planck 2018 CMB lensing at ∼30σ. We use the cross-correlation on very large scales to constrain local primordial non-Gaussianity via the scale dependence of quasar bias. The DESI quasar targets lie at an effective redshift of 1.51 and are separated into four imaging regions of varying depth and image quality. We select quasar targets from Legacy Survey DR9 imaging, apply additional flux and photometric redshift cuts to improve the purity and reduce the fraction of unclassified redshifts, and use early DESI spectroscopy of 194,000 quasar targets to determine their redshift distribution and stellar contamination fraction (2.6%). Due to significant excess large-scale power in the quasar autocorrelation, we apply weights to mitigate contamination from imaging systematics such as depth, extinction, and stellar density. We use realistic contaminated mocks to determine the greatest number of systematic modes that we can fit, before we are biased by overfitting and spuriously remove real power. We find that linear regression with one to seven imaging templates removed per region accurately recovers the input cross-power, f NL and linear bias. As in previous analyses, our f NL constraint depends on the linear primordial non-Gaussianity bias parameter, bϕ = 2(b - p)δc assuming universality of the halo mass function. We measure f NL = -26+45 -40 with p = 1.6 (f NL = -18+29 -27 with p = 1.0), and find that this result is robust under several systematics tests. Future spectroscopic quasar cross-correlations with Planck lensing can tighten the f NL constraints by a factor of 2 if they can remove the excess power on large scales in the quasar auto power spectrum.

Zero Metallicity with Zero CPU Hours: Masses of the First Stars on the Laptop

We develop an analytic model for the mass of the first stars forming in the centers of primordial gas clouds as a function of host halo mass, redshift, and degree of rotation. The model is based on the estimation of key timescales determining the following three processes: the collapse of the gas cloud, the accretion onto the protostellar core, and the radiative feedback of the protostellar core. The final stellar mass is determined by the total mass accreted until the radiative feedback halts the accretion. The analytic estimation, motivated by the result of the full numerical simulations, leads to algebraic expressions allowing an extremely fast execution. Despite its simplicity, the model reproduces the stellar mass scale and its parameter dependencies observed in state-of-the-art cosmological zoom-in simulations. This work clarifies the basic physical principles undergirding such numerical treatments and provides a path to efficiently calibrating numerical predictions against eventual observations of the first stars.

JWST lensed quasar dark matter survey – I. Description and first results

ABSTRACT The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar’s light is emitted from a physical region of the quasar with an angular scale of milliarcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10 μm. The JWST Mid-Infrared Instrument is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3 per cent precision. Our simulations suggest that this precision makes it feasible to detect the presence of 107 M⊙ dark matter haloes at cosmological distances. Such haloes are expected to be completely dark in cold dark matter models.

An environment-dependent halo mass function as a driver for the early quenching of z ≥ 1.5 cluster galaxies

ABSTRACT Many z ≈1.5 galaxies with a stellar mass (M⋆) $\ge 10^{10}\, \mathrm{M}_\odot$ are already quenched in both galaxy clusters (&amp;gt;50 per cent) and the field (&amp;gt;20 per cent), with clusters having a higher quenched fraction at all stellar masses compared to the field. A puzzling issue is that these massive quenched galaxies have stellar populations of similar age in both clusters and the field. This suggests that, despite the higher quenched fraction in clusters, the dominant quenching mechanism for massive galaxies is similar in both environments. In this work, we use data from the cosmological hydrodynamic simulations Hydrangea and EAGLE to test whether the excess quenched fraction of massive galaxies in z=1.5 clusters results from fundamental differences in their halo properties compared to the field. We find that (i) at $10^{10}\le \, M_{\star }/\mathrm{M}_\odot \, \le 10^{11}$, quenched fractions at 1.5&amp;lt;z&amp;lt;3.5 are consistently higher for galaxies with higher peak maximum circular velocity of the dark matter halo (vmax, peak), and (ii) the distribution of vmax, peak is strongly biased towards higher values for cluster satellites compared to the field centrals. Due to this difference in the halo properties of cluster and field galaxies, secular processes alone may account for (most of) the environmental excess of massive quenched galaxies in high-redshift (proto-)clusters. Taken at face value, our results challenge a fundamental assumption of popular quenching models that clusters are assembled from an unbiased subset of infalling field galaxies. If confirmed, this would imply that such models must necessarily fail at high redshift, as indicated by recent observations.

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.

The cosmological analysis of X-ray cluster surveys

Context. Cosmological studies have now entered Stage IV according to the Dark Energy Task Force (DETF) prescription. New missions (Euclid, Rubin Observatory, SRG/eROSITA) will cover very large fractions of the sky with unprecedented depth. These are expected to provide the required ultimate accuracy in the dark energy (DE) equation of state (EoS), which is required for the elucidation of the origin of the acceleration of cosmic expansion. However, none of these projects have the power to systematically unveil the galaxy cluster population in the 1 &lt; z &lt; 2 range. There therefore remains the need for an Athena-like mission to run independent cosmological investigations and scrutinise the consistency between the results from the 0 &lt; z &lt; 1 and 1 &lt; z &lt; 2 epochs. Aims. We study the constraints on the DE EoS and on primordial non-gaussanities for typical X-ray cluster surveys executed by a generic Athena-like Wide Field Imager. We focus on the impact of cluster number counts in the 1 &lt; z &lt; 2 range. Methods. We consider two survey designs: 50 deg2 at 80 ks (survey A) and 200 deg2 at 20 ks (survey B). We analytically derive cluster number counts and predict the cosmological potential of the corresponding samples, A and B, by means of a Fisher analysis. We adopt an approach that forward models the observed properties of the cluster population in the redshift–count rate–hardness ratio parameter space. Results. The achieved depth allows us to unveil the halo mass function down to the group scale out to z = 2. We predict the detection of thousands of clusters down to a few 1013h−1 M⊙, in particular 940 and 1400 clusters for surveys A and B, respectively, at z &gt; 1. Such samples will allow a detailed modelling of the evolution of cluster physics along with a standalone cosmological analysis. Our results suggest that survey B has the optimal design as it provides greater statistics. Remarkably, high-redshift clusters represent 15% or less of the full samples but contribute at a much higher level to the cosmological accuracy: by alleviating various degeneracies, these objects allow a significant reduction of the uncertainty on the cosmological parameters: Δwa is reduced by a factor of ∼2.3 and Δ fNLloc by a factor of ∼3. Conclusions. Inventorying the deep high-z X-ray cluster population can play a crucial role in ensuring overall cosmological consistency. This will be the major aim of future new-generation Athena-like missions.

The abundance of dark matter haloes down to Earth mass

ABSTRACT 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 (δ &amp;lt; −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.

Warm dark matter constraints from the JWST

ABSTRACT Warm dark matter (WDM) particles with masses (∼ kilo electronvolt) offer an attractive solution to the small-scale issues faced by the cold dark matter (CDM) paradigm. The delay of structure formation in WDM models and the associated dearth of low-mass systems at high redshifts makes this an ideal time to revisit WDM constraints in light of the unprecedented datasets from the JWST. Developing a phenomenological model based on the halo mass functions in CDM and WDM models, we calculate high redshift ($z \gtrsim 6$), the stellar mass functions (SMF) and the associated stellar mass density (SMD) and the maximum stellar mass allowed in a given volume. We find that: (i) WDM as light as 1.5 keV is already disfavoured by the low-mass end of the SMF (stellar mass $M_* \sim 10^7\,\rm M_\odot$) although caution must be exerted given the impact of lensing uncertainties; (ii) 1.5 keV WDM models predict SMD values that show a steep decrease from 108.8 to $10^{2}\,\rm M_\odot \rm cMpc^{-3}$ from z ∼ 4 to 17 for $M_* \gtrsim 10^8\,\rm M_\odot$; and (iii) the 1.5 keV WDM model predicts a sharp and earlier cut-off in the maximum stellar masses for a given number density (or volume) as compared to CDM or heavier WDM models. For example, with a number density of $10^{-3}\,\rm {cMpc^{-3}}$, 1.5 (3) KeV WDM models do not predict bound objects at $z \gtrsim 12$ (18). Forthcoming JWST observations of multiple blank fields can therefore be used as a strong probe of WDM at an epoch inaccessible by other means.

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.