The Evolution of the Stellar Mass Function of Galaxies from z = 4.0 and the First Comprehensive Analysis of its Uncertainties: Evidence for Mass-Dependent Evolution

We present the evolution of the stellar mass function (SMF) of galaxies from z = 4.0 to z = 1.3 measured from a sample constructed from the deep near-infrared Multi-wavelength Survey by Yale-Chile, the Faint Infrared Extragalactic Survey, and the Great Observatories Origins Deep Survey-Chandra Deep Field South surveys, all having very high-quality optical to mid-infrared data. This sample, unique in that it combines data from surveys with a large range of depths and areas in a self-consistent way, allowed us to (1) minimize the uncertainty due to cosmic variance and empirically quantify its contribution to the total error budget; (2) simultaneously probe the high-mass end and the low-mass end (down to ~0.05 times the characteristic stellar mass) of the SMF with good statistics; and (3) empirically derive the redshift-dependent completeness limits in stellar mass. We provide, for the first time, a comprehensive analysis of random and systematic uncertainties affecting the derived SMFs, including the effect of metallicity, extinction law, stellar population synthesis model, and initial mass function. We find that the mass density evolves by a factor of ~17⁺⁷ _-10 since z = 4.0, mostly driven by a change in the normalization Φ^sstarf. If only random errors are taken into account, we find evidence for mass-dependent evolution, with the low-mass end evolving more rapidly than the high-mass end. However, we show that this result is no longer robust when systematic uncertainties due to the SED-modeling assumptions are taken into account. Another significant uncertainty is the contribution to the overall stellar mass density of galaxies below our mass limit; future studies with WFC3 will provide better constraints on the SMF at masses below 10¹⁰ M _sun at z>2. Taking our results at face value, we find that they are in conflict with semianalytic models of galaxy formation. The models predict SMFs that are in general too steep, with too many low-mass galaxies and too few high-mass galaxies. The discrepancy at the high-mass end is susceptible to uncertainties in the models and the data, but the discrepancy at the low-mass end may be more difficult to explain.

Based on observations with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory (JPL), California Institute of Technology under NASA contract 1407. Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555. Based on observations collected at the European Southern Observatories, Chile (ESO Programme LP164.O-0612, 168.A-0485, 170.A-0788, 074.A-0709, 275.A-5060, and 171.A-3045). Based on observations obtained at the Cerro Tololo Inter-American Observatory, a division of the National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.