skip to main content

Title: A better way to define dark matter haloes

Dark matter haloes have long been recognized as one of the fundamental building blocks of large-scale structure formation models. Despite their importance – or perhaps because of it! – halo definitions continue to evolve towards more physically motivated criteria. Here, we propose a new definition that is physically motivated, effectively unique, and parameter-free: ‘A dark matter halo is comprised of the collection of particles orbiting in their own self-generated potential’. This definition is enabled by the fact that, even with as few as ≈300 particles per halo, nearly every particle in the vicinity of a halo can be uniquely classified as either orbiting or infalling based on its dynamical history. For brevity, we refer to haloes selected in this way as physical haloes. We demonstrate that (1) the mass function of physical haloes is Press–Schechter, provided the critical threshold for collapse is allowed to vary slowly with peak height; and (2) the peak-background split prediction of the clustering amplitude of physical haloes is statistically consistent with the simulation data, with accuracy no worse than ≈5 per cent.

; ; ; ; ; ; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 2464-2476
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The splashback radius, Rsp, is a physically motivated halo boundary that separates infalling and collapsed matter of haloes. We study Rsp in the hydrodynamic and dark matter-only IllustrisTNG simulations. The most commonly adopted signature of Rsp is the radius at which the radial density profiles are steepest. Therefore, we explicitly optimize our density profile fit to the profile slope and find that this leads to a $\sim 5{{\ \rm per\ cent}}$ larger radius compared to other optimizations. We calculate Rsp for haloes with masses between 1013 and 15 M⊙ as a function of halo mass, accretion rate, and redshift. Rsp decreases with mass and with redshift for haloes of similar M200 m in agreement with previous work. We also find that Rsp/R200 m decreases with halo accretion rate. We apply our analysis to dark matter, gas, and satellite galaxies associated with haloes to investigate the observational potential of Rsp. The radius of steepest slope in gas profiles is consistently smaller than the value calculated from dark matter profiles. The steepest slope in galaxy profiles, which are often used in observations, tends to agree with dark matter profiles but is lower for less massive haloes. We compare Rsp in hydrodynamic and N-body dark matter-only simulationsmore »and do not find a significant difference caused by the addition of baryonic physics. Thus, results from dark matter-only simulations should be applicable to realistic haloes.« less

    We present a per cent-level accurate model of the line-of-sight velocity distribution of galaxies around dark matter haloes as a function of projected radius and halo mass. The model is developed and tested using synthetic galaxy catalogues generated with the UniverseMachine run on the Multi-Dark Planck 2 N-body simulations. The model decomposes the galaxies around a cluster into three kinematically distinct classes: orbiting, infalling, and interloping galaxies. We demonstrate that: (1) we can statistically distinguish between these three types of galaxies using only projected line-of-sight velocity information; (2) the halo edge radius inferred from the line-of-sight velocity dispersion is an excellent proxy for the three-dimensional halo edge radius; and (3) we can accurately recover the full velocity dispersion profile for each of the three populations of galaxies. Importantly, the velocity dispersion profiles of the orbiting and infalling galaxies contain five independent parameters – three distinct radial scales and two velocity dispersion amplitudes – each of which is correlated with mass. Thus, the velocity dispersion profile of galaxy clusters has inherent redundancies that allow us to perform non-trivial systematics checks from a single data set. We discuss several potential applications of our new model for detecting the edge radius and constraining cosmology andmore »astrophysics using upcoming spectroscopic surveys.

    « less

    We analyse strongly lensed images in eight galaxy clusters to measure their dark matter density profiles in the radial region between 10 kpc and 150 kpc, and use this to constrain the self-interaction cross-section of dark matter (DM) particles. We infer the mass profiles of the central DM haloes, bright central galaxies, key member galaxies, and DM subhaloes for the member galaxies for all eight clusters using the qlens code. The inferred DM halo surface densities are fit to a self-interacting dark matter model, which allows us to constrain the self-interaction cross-section over mass σ/m. When our full method is applied to mock data generated from two clusters in the Illustris-TNG simulation, we find results consistent with no dark matter self-interactions as expected. For the eight observed clusters with average relative velocities of $1458_{-81}^{+80}$ km s−1, we infer $\sigma /m = 0.082_{-0.021}^{+0.027} \rm cm^2\, g^{ -1}$ and $\sigma /m \lt 0.13~ \rm cm^2\, g^{ -1}$ at the 95 per cent confidence level.


    Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on amore »modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

    « less

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

    « less