Unveiling the population of massive quenched galaxies at zge2 in the COLIBRE simulations -- I. Galaxy demographics

JWST has uncovered a substantial population of Massive ($M_{\star} \gtrsim 10^{10}\,\mathrm{M_{\odot}}$), Quenched Galaxies (MQGs) in the early Universe ($z \gtrsim 2$), whose properties challenge current galaxy formation models. In this series, we examine this population of MQGs within the new COLIBRE cosmological hydrodynamical simulations, which introduce key innovations in their sub-grid physics. In this first paper, we find a dependence of MQG number densities on both mass resolution and the Active Galactic Nucleus (AGN) feedback implementation, as well as a significant impact from potential observational uncertainties. Using the fiducial $(200\,\rm cMpc)^3$ volume L200m6 simulation, which provides adequate volume, mass and spatial resolution to study these systems, we report number densities and stellar mass functions in broad agreement with the latest observations. The predicted quenching and formation timescales are qualitatively consistent with observational inferences, indicating extended formation (medians $t_{50}\approx0.5-1.5\,\mathrm{Gyr}$) followed by rapid quenching (medians $t_{\mathrm{q}}\lesssim0.6\,\mathrm{Gyr}$) with strong starburst episodes. Leveraging the state-of-the-art physics in COLIBRE, the model predicts that MQGs have dust and $\rm H_{2}$ fractions more than $1$ dex lower than their massive star-forming counterparts; generally consistent with the (scarce) observational estimates. MQGs and massive star-forming systems show broadly similar sizes and kinematics, suggesting that size or morphological transformations occur after quenching in COLIBRE. Our results provide robust predictions for MQGs and show that tensions with observations are reduced when an effective observational uncertainty is forward-modelled.