A$^3$COSMOS: the infrared luminosity function and dust-obscured star formation rate density at $0.5<z<6$

TL;DR

This study uses the extensive A$^3$COSMOS ALMA survey to analyze the infrared luminosity function and dust-obscured star formation rate density of galaxies from redshift 0.5 to 6, revealing evolutionary trends and extending high-redshift constraints.

Contribution

The paper introduces a new method to correct for inhomogeneous survey sampling, enabling accurate derivation of luminosity functions and star formation rates across a wide redshift range.

Findings

Galaxy population is mostly massive, IR-bright, and highly star-forming.

The IR luminosity function shows evolution with decreasing density and increasing luminosity at higher redshifts.

Star formation rate density peaks between redshift 1 and 3, declining at higher redshifts.

Abstract

Aims: We leverage the largest available Atacama Large Millimetre/submillimetre Array (ALMA) survey from the archive (A$^3$COSMOS) to study to study infrared luminosity function and dust-obscured star formation rate density of sub-millimeter/millimeter (sub-mm/mm) galaxies from $z=0.5\,-\,6$. Methods: The A$^3$COSMOS survey utilizes all publicly available ALMA data in the COSMOS field, therefore having inhomogeneous coverage in terms of observing wavelength and depth. In order to derive the luminosity functions and star formation rate densities, we apply a newly developed method that corrects the statistics of an inhomogeously sampled survey of individual pointings to those representing an unbiased blind survey. Results: We find our sample to mostly consist of massive ($M_{\star} \sim 10^{10} - 10^{12}$ $\rm M_{\odot}$), IR-bright ($L_* \sim 10^{11}-10^{13.5} \rm L_{\odot}$), highly star-forming (SFR $\sim 100-1000$ $\rm M_{\odot}$ $\rm yr^{-1}$) galaxies. We find an evolutionary trend in the typical density ($\Phi^*$) and luminosity ($L^*$) of the galaxy population, which decrease and increase with redshift, respectively. Our IR LF is in agreement with previous literature results and we are able to extend to high redshift ($z > 3$) the constraints on the knee and bright-end of the LF, derived by using the Herschel data. Finally, we obtain the SFRD up to $z\sim 6$ by integrating the IR LF, finding a broad peak from $z \sim 1$ to $z \sim 3$ and a decline towards higher redshifts, in agreement with recent IR/mm-based studies, within the uncertainties, thus implying the presence of larger quantities of dust than what is expected by optical/UV studies.