COSMOS2020: The galaxy stellar mass function
Authors: Weaver J. R. , Davidzon I. , Toft S. , Ilbert O. , McCracken H. J. , Gould K. M. L. , Jespersen C. K. , Steinhardt C. , Lagos C. D. P. , Capak P. L. , Casey C. M. , Chartab N. , Faisst A. L. , Hayward C. C. , Kartaltepe J. S. , Kauffmann O. B. , Koekemoer A. M. , Kokorev V. , Laigle C. , Liu D. , Long A. , Magdis G. E. , McPartland C. J. R. , Milvang-Jensen B. , Mobasher B. , Moneti A. , Peng Y. , Sanders D. B. , Shuntov M. , Sneppen A. , Valentino F. , Zalesky L. , Zamorani G. 
Published: May 2023

Context. How galaxies form, assemble, and cease their star formation is a central question within the modern landscape of galaxy evolution studies. These processes are indelibly imprinted on the galaxy stellar mass function (SMF), and its measurement and understanding is key to uncovering a unified theory of galaxy evolution.
Aims. We present constraints on the shape and evolution of the galaxy SMF, the quiescent galaxy fraction, and the cosmic stellar mass density across 90% of the history of the Universe from z = 7.5 → 0.2 as a means to study the physical processes that underpin galaxy evolution.
Methods. The COSMOS survey is an ideal laboratory for studying representative galaxy samples. Now equipped with deeper and more homogeneous near-infrared coverage exploited by the COSMOS2020 catalog, we leverage the large 1.27 deg2 effective area to improve sample statistics and understand spatial variations (cosmic variance) – particularly for rare, massive galaxies – and push to higher redshifts with greater confidence and mass completeness than previous studies. We divide the total stellar mass function into star-forming and quiescent subsamples through NUVrJ color-color selection. The measurements are then fit with single- and double-component Schechter functions to infer the intrinsic galaxy stellar mass function, the evolution of its key parameters, and the cosmic stellar mass density out to z = 7.5. Finally, we compare our measurements to predictions from state-of-the-art cosmological simulations and theoretical dark matter halo mass functions.
Results. We find a smooth, monotonic evolution in the galaxy stellar mass function since z = 7.5, in general agreement with previous studies. The number density of star-forming systems have undergone remarkably consistent growth spanning four decades in stellar mass from z = 7.5 → 2 whereupon high-mass systems become predominantly quiescent (“downsizing”). Meanwhile, the assembly and growth of low-mass quiescent systems only occurred recently, and rapidly. An excess of massive systems at z ≈ 2.5 − 5.5 with strikingly red colors, with some being newly identified, increase the observed number densities to the point where the SMF cannot be reconciled with a Schechter function.
Conclusions. Systematics including cosmic variance and/or active galactic nuclei contamination are unlikely to fully explain this excess, and so we speculate that they may be dust-obscured populations similar to those found in far infrared surveys. Furthermore, we find a sustained agreement from z ≈ 3 − 6 between the stellar and dark matter halo mass functions for the most massive systems, suggesting that star formation in massive halos may be more efficient at early times.