R Endsley

The variety of star formation histories (SFHs) of z ≳ 6 galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations suggest that many z ∼ 6–9 galaxies are dominated by ≳200 Myr stellar populations, implying significant star formation at z ≳ 9, while others find that most reionization era galaxies are ≲10 Myr, consistent with little z ≳ 9 star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{\rm \small UV}\lesssim -21$) galaxies colour-selected to lie at z ≃ 6.6–6.9, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of ∼30–70 Myr depending on SED model. We quantify the fractions of ≤10 and ≥250 Myr galaxies, finding that these systems comprise ∼15–30 per cent and ∼20–25 per cent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied by the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for z ≳ 9 stellar mass assembly and discuss improvements expected from JWST.

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.