2011ApJ...742...24L

Exploiting the Herschel Astrophysical Terahertz Large Area Survey Science Demonstration Phase survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 µm and at several redshifts z ≳ 1, for bright submillimeter galaxies with star formation rates (SFRs) ≳ 100 M_☉/yr. We find that the evolution of the comoving LF is strong up to z ~ 2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H ≳ 3x10¹² M_☉) amounts to ∼7x10⁸ yr. Given the SFRs, which are in the range of 10²-10³ M_☉/yr, this timescale implies final stellar masses of the order of 10¹¹-10¹² M_☉. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z ≳ 1. The comparison of the statistics for submillimeter and UV-selected galaxies suggests that the dust-free, UV bright phase is ≳ 10² times shorter than the submillimeter bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference spectral energy distribution (SED; the one of the z ~ 2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts >1 but also the counts at wavelengths ranging from 250 µm to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of submillimeter SEDs of z > 1 galaxies around the reference one is rather small.