Mid- and Far-infrared Luminosity Functions and Galaxy Evolution from Multiwavelength Spitzer Observations up to z~2.5
[Abridged]We exploit a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent 24, 15, 12 and 8micron monochromatic and bolometric IR luminosity functions (LF) over the 0<z<2.5 redshift range. Our analysis is based on the combination of data from deep Spitzer surveys in the VVDS-SWIRE and GOODS areas. To our limiting flux of S(24)=400microJy our derived sample in VVDS-SWIRE includes 1494 sources, and 666 and 904 sources brighter than S(24)=80microJy are catalogued in GOODS-S and GOODS-N, respectively, for a total area of ~0.9 square degs. We obtain reliable optical identifications and redshifts, providing us a rich and robust dataset for our luminosity function determination. Based on the multi-wavelength information available, we constrain the LFs at 8, 12, 15 and 24micron. We also extrapolate total IR luminosities from our best-fit to the observed SEDs of each source, and use this to derive the bolometric LF and comoving volume emissivity up to z~2.5. In the 0<z<1 interval, the bolometric IR luminosity density evolves as (1+z)^3.8+/-0.4. Although more uncertain at higher-z, our results show a flattening of the IR luminosity density at z>1. The mean redshift of the peak in the source number density shifts with luminosity: the brighest IR galaxies appear to be forming stars earlier in cosmic time (z>1.5), while the less luminous ones keep doing it at more recent epochs (z~1 for L(IR)<10^11L_sun). Our results suggest a rapid increase of the galaxy IR comoving volume emissivity back to z~1 and a constant average emissivity at z>1. We also seem to find a difference in the evolution rate of the source number densities as a function of luminosity, a downsizing evolutionary pattern similar to that reported from other samples of cosmic sources.