Abstract

Our current understanding of the cosmic star formation history at z > 3 is primarily based on UV-selected galaxies (Lyman-break galaxies, i.e., LBGs). Recent studies of H-dropouts (HST-dark galaxies) have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z > 3. In this work, we extend the H-dropout criterion to lower masses to select optically dark or faint galaxies (OFGs) at high redshifts in order to complete the census between LBGs and H-dropouts. Our criterion (H > 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies (typically E(B - V) > 0.4) with lower stellar masses at high redshifts. In addition, with this criterion, our sample is not contaminated by massive passive or old galaxies. In total, we identified 27 OFGs at (Zphot) > 3 (with a median of z(med) = 4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log(M-star/M-circle dot) = 9.4-11.1 (with a median of log(M-star med/M-circle dot) = 10.3). We find that up to 75% of the OFGs with log(M-star/M-circle dot) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing an optical-to-millimeter stacking analysis of the OFGs, we find that rather than being limited to a rare population of extreme starbursts, these OFGs represent a normal population of dusty star-forming galaxies at z > 3. The OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than the scaling relation of typical star-forming galaxies. Additionally, the total star formation rate (SFRtot = SFRIR + SFRUV) of the stacked OFGs is much higher than the SFRUVcorr (SFRUV corrected for dust extinction), with an average SFRtot/SFRUVcorr = 8 +/- 1, which lies above (similar to 0.3 dex) the 16-84th percentile range of typical star-forming galaxies at 3 <= z <= 6. All of the above suggests the presence of hidden dust regions in the OFGs that absorb all UV photons, which cannot be reproduced with dust extinction corrections. The effective radius of the average dust size measured by a circular Gaussian model fit in the uv plane is R-e(1.13 mm) = 1.01 +/- 0.05 kpc. After excluding the five LBGs in the OFG sample, we investigated their contributions to the cosmic star formation rate density (SFRD). We found that the SFRD at z > 3 contributed by massive OFGs (log(M-star/M-circle dot) > 10.3) is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculated the combined contribution of OFGs and LBGs to the cosmic SFRD at z = 4-5 to be 4 x 10(-2) M-circle dot yr(-1) Mpc(-3), which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift. This value could be even larger, as our calculations were performed in a very conservative way.