Search

The star formation and metal enrichment histories of galaxies - at any epoch - constitute one of the key properties of galaxies, and their measurement is a core aim of observational extragalactic astronomy. The lack of deep rest-frame optical coverage at high redshift has made robust constraints elusive, but this is now changing thanks to JWST. In preparation for the constraints provided by JWST, we explore the star formation and metal enrichment histories of galaxies at z = 5-13 using the First Light And Reionization Epoch Simulations (FLARES) suite. Built on the EAGLE model, the unique strategy of FLARES allows us to simulate galaxies with a wide range of stellar masses (and luminosities) and environments. While we predict significant redshift evolution of average ages and specific star formation rates, our core result is mostly a flat relationship of age and specific star formation rate with stellar mass. We also find that galaxies in this epoch predominantly have strongly rising star formation histories, albeit with the normalization dropping with redshift and stellar mass. In terms of chemical enrichment, we predict a strong stellar mass-metallicity relation present at z = 10 and beyond alongside significant a-enhancement. Finally, we find no large-scale environmental dependence of the relationship between age, specific star formation rate, or metallicity with stellar mass.

The history of reionization is highly dependent on the ionizing properties of high-redshift galaxies. It is therefore important to have a solid understanding of how the ionizing properties of galaxies are linked to physical and observable quantities. In this paper, we use the First Light and Reionization Epoch Simulations (Flares) to study the Lyman-continuum (LyC, i.e. hydrogen-ionizing) emission of massive (⁠M∗>108M⊙⁠) galaxies at redshifts z = 5 − 10. We find that the specific ionizing emissivity (i.e. intrinsic ionizing emissivity per unit stellar mass) decreases as stellar mass increases, due to the combined effects of increasing age and metallicity.

In the Flares (First Light And Reionisation Epoch Simulations) suite of hydrodynamical simulations, we find the high redshift (z > 5) intrinsic size-luminosity relation is, surprisingly, negatively sloped. However, after including the effects of dust attenuation we find a positively sloped UV observed size-luminosity relation in good agreement with other simulated and observational studies. In this work, we extend this analysis to probe the underlying physical mechanisms driving the formation and evolution of the compact galaxies driving the negative size-mass/size-luminosity relation. We find the majority of compact galaxies (R1/2, ⋆ < 1 physical kpc, which drive the negative slope of the size-mass relation, have transitioned from extended to compact sizes via efficient centralised cooling, resulting in high specific star formation rates in their cores. These compact stellar systems are enshrouded by non-star forming gas distributions as much as 100 × larger than their stellar counterparts. By comparing with galaxies from the Eagle simulation suite, we find that these extended gas distributions ‘turn on’ and begin to form stars between z = 5 and z = 0 leading to increasing sizes, and thus the evolution of the size-mass relation from a negative to a positive slope. This explicitly demonstrates the process of inside-out galaxy formation in which compact bulges form earlier than the surrounding discs.

Passive galaxies are ubiquitous in the local universe, and various physical channels have been proposed that lead to this passivity. To date, robust passive galaxy candidates have been detected up to z ≤ 5, but it is still unknown if they exist at higher redshifts, what their relative abundances are, and what causes them to stop forming stars. We present predictions from the first light and reionisation epoch simulations (flares), a series of zoom simulations of a range of overdensities using the eagle code. Passive galaxies occur naturally in the eagle model at high redshift, and are in good agreement with number density estimates from Hubble Space Telescope (HST) and early JWST results at 3 ≤ z ≤ 5. Due to the unique flares approach, we extend these predictions to higher redshifts, finding passive galaxy populations up to z ∼ 8. Feedback from supermassive black holes is the main driver of passivity, leading to reduced gas fractions and star forming gas reservoirs. We find that passive galaxies at z ≥ 5 are not identified in the typical UVJ selection space due to their still relatively young stellar populations, and present new rest-frame selection regions. We also produce mock NIRCam and MIRI fluxes, and find that significant numbers of passive galaxies at z ≥ 5 should be detectable in upcoming wide surveys with JWST. Finally, we present JWST colour distributions, with new selection regions in the observer–frame for identifying these early passive populations.

We present results on the star cluster properties from a series of high resolution smoothed particles hydrodynamics (SPH) simulations of isolated dwarf galaxies as part of the GRIFFIN project. The simulations at sub-parsec spatial resolution and a minimum particle mass of 4 M-circle dot incorporate non-equilibrium heating, cooling, and chemistry processes, and realize individual massive stars. The simulations follow feedback channels of massive stars that include the interstellar-radiation field variable in space and time, the radiation input by photo-ionization and supernova explosions. Varying the star formation efficiency per free-fall time in the range epsilon(ff) = 0.2-50 per cent neither changes the star formation rates nor the outflow rates. While the environmental densities at star formation change significantly with epsilon(ff), the ambient densities of supernovae are independent of epsilon(ff) indicating a decoupling of the two processes. At low epsilon(ff), gas is allowed to collapse more before star formation, resulting in more massive, and increasingly more bound star clusters are formed, which are typically not destroyed. With increasing epsilon(ff), there is a trend for shallower cluster mass functions and the cluster formation efficiency Gamma for young bound clusters decreases from 50 per cent to similar to 1 per cent showing evidence for cluster disruption. However, none of our simulations form low mass (<10(3) M-circle dot) clusters with structural properties in perfect agreement with observations. Traditional star formation models used in galaxy formation simulations based on local free-fall times might therefore be unable to capture star cluster properties without significant fine tuning.

We present the intrinsic and observed sizes of galaxies at z ≥ 5 in the First Light And Reionisation Epoch Simulations (FLARES). We employ the large effective volume of FLARES to produce a sizeable sample of high redshift galaxies with intrinsic and observed luminosities and half light radii in a range of rest frame UV and visual photometric bands. This sample contains a significant number of intrinsically ultra-compact galaxies in the far-UV (1500 Å), leading to a negative intrinsic far-UV size-luminosity relation. However, after the inclusion of the effects of dust these same compact galaxies exhibit observed sizes that are as much as 50 times larger than those measured from the intrinsic emission, and broadly agree with a range of observational samples. This increase in size is driven by the concentration of dust in the core of galaxies, heavily attenuating the intrinsically brightest regions. At fixed luminosity we find a galaxy size redshift evolution with a slope of m = 1.21 − 1.87 depending on the luminosity sample in question, and we demonstrate the wavelength dependence of the size-luminosity relation which will soon be probed by the Webb Space Telescope.

We present the intrinsic and observed sizes of galaxies at z >= 5 in the First Light And Reionisation Epoch Simulations (flares). We employ the large effective volume of flares to produce a sizeable sample of high-redshift galaxies with intrinsic and observed luminosities and half-light radii in a range of rest-frame ultraviolet (UV) and visual photometric bands. This sample contains a significant number of intrinsically ultracompact galaxies in the far-UV (1500 angstrom), leading to a negative intrinsic far-UV size-luminosity relation. However, after the inclusion of the effects of dust these same compact galaxies exhibit observed sizes that are as much as 50 times larger than those measured from the intrinsic emission, and broadly agree with a range of observational samples. This increase in size is driven by the concentration of dust in the core of galaxies, heavily attenuating the intrinsically brightest regions. At fixed luminosity we find a galaxy size redshift evolution with a slope of m = 1.21-1.87 depending on the luminosity sample in question, and we demonstrate the wavelength dependence of the size-luminosity relation that will soon be probed by the James Webb Space Telescope.

With its exquisite sensitivity, wavelength coverage, and spatial and spectral resolution, the James Webb Space Telescope (JWST) is poised to revolutionize our view of the distant, high-redshift (z > 5) Universe. While Webb’s spectroscopic observations will be transformative for the field, photometric observations play a key role in identifying distant objects and providing more comprehensive samples than accessible to spectroscopy alone. In addition to identifying objects, photometric observations can also be used to infer physical properties and thus be used to constrain galaxy formation models. However, inferred physical properties from broad-band photometric observations, particularly in the absence of spectroscopic redshifts, often have large uncertainties. With the development of new tools for forward modelling simulations, it is now routinely possible to predict observational quantities, enabling a direct comparison with observations. With this in mind, in this work, we make predictions for the colour evolution of galaxies at z = 5–15 using the First Light And Reionisation Epoch Simulations (flares) cosmological hydrodynamical simulation suite. We predict a complex evolution with time, driven predominantly by strong nebular line emission passing through individual bands. These predictions are in good agreement with existing constraints from Hubble and Spitzer as well as some of the first results from Webb. We also contrast our predictions with other models in the literature: While the general trends are similar, we find key differences, particularly in the strength of features associated with strong nebular line emission. This suggests photometric observations alone should provide useful discriminating power between different models and physical states of galaxies.