Pratika Dayal

The earliest JWST observations have revealed an unexpected abundance of super-early (z > 10), massive ($M_*\, \approx 10^9 {\rm M}_{\odot }$ ) galaxies at the bright-end (MUV ≈ −21) of the ultraviolet luminosity function (UV LF). We present a minimal physical model that explains the observed galaxy abundance at z = 10–14. The model primarily combines (i) the halo mass function, with (ii) an obscured star formation fraction prescription that is consistent with findings of the ALMA REBELS dusty galaxy survey. It has been successfully tested on well-known UV LFs up to z = 7. We argue that the weak evolution from z = 7 to z ≈ 14 of the LF bright-end can arise from a conspiracy between a decreasing dust attenuation, making galaxies brighter, that almost exactly compensates for the increasing shortage of their host halos. Our minimal model naturally reproduces the z = 10–14 LF if galaxies at $z\lower.5ex\rm{\,\, \buildrel\gt \over \sim \,\,}11$ contain a negligible amounts of dust. We speculate that dust could have been efficiently ejected during the very first phases of galaxy build-up.

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman break galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at z ≃ 7. We find that the galaxy dust mass, Md, is regulated by the combination of Type II supernova dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard time-scale τ0 = 30 Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07\!-\!0.1{{\ \rm per\ cent}}$ and a UV-to-total star formation rate relation such that log(ψUV) = −0.05 [log(ψ)]2 + 0.86 log(ψ) − 0.05 (implying that 55–80 per cent of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9}\!-\!10^{10} \rm M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at z = 7. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94\!-\!1.1{{\ \rm per\ cent}}$) that are up to 18 times larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing τ0 to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor of ∼2. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions.