The mass function dependence on the dynamical state of dark matter haloes

¹ Max-Planck-Institut für extraterrestrische Physik (MPE), Giessenbachstrasse 1, 85748 Garching bei München, Germany
e-mail: rseppi@mpe.mpg.de
² Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía, 18080 Granada, Spain
³ Astronomy Department, New Mexico State University, Las Cruces, NM, USA
⁴ Department of Astronomy, University of Virginia, Charlottesville, VA, USA

Received: 7 August 2020
Accepted: 25 May 2021

Abstract

Context. Galaxy clusters are luminous tracers of the most massive dark matter haloes in the Universe. To use them as a cosmological probe, a detailed description of the properties of dark matter haloes is required.

Aims. We characterize how the dynamical state of haloes impacts the dark matter halo mass function at the high-mass end (i.e., for haloes hosting clusters of galaxies).

Methods. We used the dark matter-only MultiDark suite of simulations and the high-mass objects M > 2.7 × 10¹³ M_⊙ h⁻¹ therein. We measured the mean relations of concentration, offset, and spin as a function of dark matter halo mass and redshift. We investigated the distributions around the mean relations. We measured the dark matter halo mass function as a function of offset, spin, and redshift. We formulated a generalized mass function framework that accounts for the dynamical state of the dark matter haloes.

Results. We confirm the recent discovery of the concentration upturn at high masses and provide a model that predicts the concentration for different values of mass and redshift with one single equation. We model the distributions around the mean values of concentration, offset, and spin with modified Schechter functions. We find that the concentration of low-mass haloes shows a faster redshift evolution compared to high-mass haloes, especially in the high-concentration regime. We find that the offset parameter is systematically smaller at low redshift, in agreement with the relaxation of structures at recent times. The peak of its distribution shifts by a factor of ∼1.5 from z = 1.4 to z = 0. The individual models are combined into a comprehensive mass function model, which predicts the mass function as a function of spin and offset. Our model recovers the fiducial mass function with ∼3% accuracy at redshift 0 and accounts for redshift evolution up to z ∼ 1.5.

Results. This new approach accounts for the dynamical state of the halo when measuring the halo mass function. It offers a connection with dynamical selection effects in galaxy cluster observations. This is key toward precision cosmology using cluster counts as a probe.