Galaxy luminosity functions from far-UV to submillimetre at $z=0$ in the COLIBRE simulations

TL;DR

This paper uses the COLIBRE simulations combined with SKIRT radiative transfer to accurately predict galaxy luminosity functions across a wide wavelength range at z=0, showing excellent agreement with observational data.

Contribution

It demonstrates that the COLIBRE simulations, when post-processed with dust modeling, can reliably reproduce galaxy luminosity functions from UV to submillimetre wavelengths without calibration.

Findings

Good convergence across different simulation resolutions.

Excellent match with observed luminosity functions from FUV to submillimetre.

Underprediction of MIR-bright and very luminous galaxies at the bright end.

Abstract

We present predictions from the recent COLIBRE cosmological hydrodynamical simulations of galaxy formation for the present-day galaxy luminosity functions (LFs) at wavelengths ranging from the far-ultraviolet (FUV) to the submillimetre. The simulations are post-processed with the radiative transfer code SKIRT, accounting for dust attenuation and emission using the distribution and properties of dust grains predicted directly by COLIBRE. Results from simulations varying in mass resolution by a factor of $\sim 10^2$ ($\sim 10^5 - 10^7\,\mathrm{M_{\odot}}$) show very good convergence over most luminosity ranges. The COLIBRE-SKIRT LFs match the data remarkably well from the FUV to the near-infrared ($3.4\,\mathrm{\mu m}$) and also in the far-infrared and submillimetre wavelength range ($70-850\,\mathrm{\mu m}$). In the mid-infrared (MIR; $8-24\,\mathrm{\mu m}$), COLIBRE-SKIRT matches the data well at low luminosities but significantly underpredicts the luminosities of MIR-bright galaxies, with the discrepancy increasing towards longer wavelengths. The total infrared LF, obtained by integrating the spectral energy distributions over $8-1000\,\mathrm{\mu m}$, also matches observations well at the faint end but underpredicts the number of very bright galaxies. The unprecedented agreement at all other wavelengths indicates that COLIBRE, coupled with this calibration-free SKIRT post-processing framework, successfully predicts the properties of stellar populations at the present day and the amount and distribution of interstellar dust.