The galaxy ultraviolet luminosity function from $z=7$ to $15$ in the COLIBRE simulations

Abstract: JWST has enabled the detection of galaxies in the earliest stages of cosmic history. We compare the ultraviolet luminosity functions (UVLFs) at redshifts $z=7-15$ predicted by the new cosmological hydrodynamics simulations, COLIBRE with observations, including those from JWST. The UV luminosities of COLIBRE galaxies are derived using the radiative transfer code SKIRT, which tracks stellar emission and its processing through the multi-phase interstellar medium and dust distribution predicted by COLIBRE. We find that although COLIBRE is consistent with the observed evolution of the stellar mass function up to $z=12$, its dust-attenuated UVLFs fall systematically below the observations at the bright end: at the number density of $10^{-6}\,\mathrm{Mpc^{-3}\,mag^{-1}}$, the brightest galaxies are underluminous by $\approx 1\,\rm mag$ at $z=7$, increasing to $\approx 2.5\,\rm mag$ at $z=15$. Accounting for observational uncertainties brings the COLIBRE UVLFs closer to the observational data, but does not fully resolve the discrepancy. Ignoring dust attenuation allows COLIBRE to produce sufficiently bright galaxies at $7\lesssim z \lesssim 12$, while at $z=15$, COLIBRE still underpredicts the luminosities of the brightest galaxies, indicating the need for additional physical mechanisms to boost the UV luminosities at the earliest cosmic epochs, such as a ''top-heavy'' stellar initial mass function. We fit the COLIBRE UVLFs with Schechter functions and calculate the evolution of the best-fit parameters. We find that the galaxy number density decreases, the characteristic luminosity becomes fainter and the faint-end slope becomes steeper towards higher redshifts. The UV luminosity density decreases by a factor of $\approx 300$ from $z = 7$ to $z = 15$.