Accurate mass and velocity functions of dark matter haloes

N-body cosmological simulations are an essential tool to understand the observed distribution of galaxies. We use the MultiDark simulation suite, run with the Planck cosmological parameters, to revisit the mass and velocity functions. At redshift z = 0, the simulations cover four orders of magnitude in halo mass from ~1011M¿ with 8783 874 distinct haloes and 532 533 subhaloes. The total volume used is ~515 Gpc3, more than eight times larger than in previous studies. We measure and model the halo mass function, its covariance matrix w.r.t halo mass and the large-scale halo bias.With the formalism of the excursion-set mass function, we explicit the tight interconnection between the covariance matrix, bias and halo mass function. We obtain a very accurate ( < 2 per cent level) model of the distinct halo mass function. We also model the subhalo mass function and its relation to the distinct halo mass function. The set of models obtained provides a complete and precise framework for the description of haloes in the concordance Planck cosmology. Finally, we provide precise analytical fits of the Vmax maximum velocity function up to redshift z < 2.3 to push for the development of halo occupation distribution using Vmax. The data and the analysis code are made publicly available in the Skies and Universes data base. © 2017 The Authors.JC and FP acknowledge support from the Spanish MICINNs Consolider-Ingenio 2010 Programme under grant MD CSD2009-00064, MINECO Centro de Excelencia Severo Ochoa Programme under the grants SEV-20120249, FPA2012-34694, and the projects AYA2014-60641-C2-1-P and AYA2012-31101. GY acknowledges financial support from MINECO/FEDER (Spain) under project number AYA2012-31101 and AYA2015-63810-P.Peer Reviewe