DUCA: Dynamic Universe Cosmological Analysis. I. The halo mass function in dynamical dark energy cosmologies

The halo mass function (HMF) is fundamental for interpreting the number counts of galaxy clusters, serving as a pivotal theoretical tool in cosmology. With the advent of high-precision surveys such as LSST, eROSITA, DESI, and Euclid, accurate HMF modeling becomes indispensable to avoid systematic biases in cosmological parameter estimation from cluster cosmology. Moreover, these surveys aim to shed light on the dark sector and uncover dark energy's puzzling nature, necessitating models that faithfully capture its features to ensure robust parameter inference. We aim to construct a model for the HMF in dynamical dark energy cosmologies that preserves the accuracy achieved for the standard $\Lambda (\nu)$CDM model of cosmology, while meeting the precision requirements necessary for future cosmological surveys. Our approach models the HMF parameters as functions of the deceleration parameter at the turnaround, a quantity shown to encapsulate essential information regarding the impact of dynamical dark energy on structure formation. We calibrate the model using results from a comprehensive suite of $N$-body simulations spanning various cosmological scenarios, ensuring sub-percent systematic accuracy. We present an HMF model tailored for dynamical dark energy cosmologies. The model is calibrated following a Bayesian approach, and its uncertainty is characterized by a single parameter controlling its systematic error, which remains at the sub-percent level. This ensures that theoretical uncertainties from our model are subdominant relative to other error sources in future cluster number counts analyses.