2023MNRAS.522.3926B

We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a ultraviolet (UV)-selected galaxy sample with the Atacama Large Millimeter Array (ALMA) spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at z = 6.4-7.7 with [C II] 158 $\rm{\mu m}$ line scans. 16 sources exhibit dust detection, 15 of which are also spectroscopically confirmed through the [C II] line. The infrared (IR) luminosities of the sample range from log L_IR/L_⊙ = 11.4 to 12.2. Using the UV luminosity function as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/\mathrm{ L}_\odot =11.6^{+0.2}_{-0.1}$ . Our observational results are in broad agreement with the average of predicted IRLFs from simulations at z ∼ 7. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from z ∼ 4 to z ∼ 7, into the reionization epoch. The IR obscured contribution to the cosmic star formation rate density at z ∼ 7 amounts to $\mathrm{log(SFRD/{\rm M}_{\odot }\,yr^{-1}\,Mpc^{-3}) = -2.66^{+0.17}_{-0.14} }$ that is at least ∼10 per cent of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations.