Marchetti, L. and Vaccari, M. and Franceschini, A. and Arumugam, V. and Aussel, H. and Béthermin, M. and Bock, J. and Boselli, A. and Buat, V. and Burgarella, D. and Clements, D. L. and Conley, A. and Conversi, L. and Cooray, A. and Dowell, C. D. and Farrah, D. and Feltre, A. and Glenn, J. and Griffin, M. and Hatziminaoglou, E. and Heinis, S. and Ibar, E. and Ivison, R. J. and Nguyen, H. T. and O'Halloran, B. and Oliver, S. J. and Page, M. J. and Papageorgiou, A. and Pearson, C. P. and Pérez-Fournon, I. and Pohlen, M. and Rigopoulou, D. and Roseboom, I. G. and Rowan-Robinson, M. and Schulz, B. and Scott, Douglas and Seymour, N. and Shupe, D. L. and Smith, A. J. and Symeonidis, M. and Valtchanov, I. and Viero, M. and Wang, L. and Wardlow, J. and Xu, C. K. and Zemcov, M. (2016) The HerMES submillimetre local and low-redshift luminosity functions. Monthly Notices of the Royal Astronomical Society, 456 (2). pp. 1999-2023. ISSN 0035-8711
Abstract
We used wide-area surveys over 39 deg2 by the HerMES (Herschel Multi-tiered Extragalactic Survey) collaboration, performed with the Herschel Observatory SPIRE multiwavelength camera, to estimate the low-redshift, 0.02 <z <0.5, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500 μm. Within this redshift interval, we detected 7087 sources in five independent sky areas, ˜40 per cent of which have spectroscopic redshifts, while for the remaining objects photometric redshifts were used. The SPIRE LFs in different fields did not show any field-to-field variations beyond the small differences to be expected from cosmic variance. SPIRE flux densities were also combined with Spitzer photometry and multiwavelength archival data to perform a complete spectral energy distribution fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (IR; 8-1000 μm) LFs and their low-z evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The LFs show significant and rapid luminosity evolution already at low redshifts, 0.02 <z <0.2, with L_{IR}^{*} ∝ (1+z)^{6.0± 0.4} and Φ _{IR}^{*} ∝ (1+z)^{-2.1± 0.4}, L_{250}^{*} ∝ (1+z)^{5.3± 0.2} and Φ _{250}^{*} ∝ (1+z)^{-0.6± 0.4} estimated using the IR bolometric and the 250 μm LFs, respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter initial mass function and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate an SFRD scaling of SFRD0 + 0.08z, where SFRD0 ≃ (1.9 ± 0.03) × 10-2 [M⊙ Mpc-3] is our total SFRD estimate at z ˜ 0.02.