skip to main content

Title: Star formation histories of UV-luminous galaxies at z ≃ 6.8: implications for stellar mass assembly at early cosmic times

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.

more » « less
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Medium: X Size: p. 5859-5881
["p. 5859-5881"]
Sponsoring Org:
National Science Foundation
More Like this

    With JWST, new opportunities to study the evolution of galaxies in the early Universe are emerging. Spitzer constraints on rest-optical properties of z ≳ 7 galaxies demonstrated the power of using galaxy stellar masses and star formation histories (SFHs) to indirectly infer the cosmic star formation history. However, only the brightest individual z ≳ 8 objects could be detected with Spitzer, making it difficult to robustly constrain activity at z ≳ 10. Here, we leverage the greatly improved rest-optical sensitivity of JWST at z ≳ 8 to constrain the ages of seven UV-bright ($M_{\rm uv}\lesssim -19.5$) galaxies selected to lie at z ∼ 8.5–11, then investigate implications for z ≳ 15 star formation. We infer the properties of individual objects with two spectral energy distribution modelling codes, then infer a distribution of ages for bright z ∼ 8.5–11 galaxies. We find a median age of ∼20 Myr, younger than that inferred at z ∼ 7 with a similar analysis, consistent with an evolution towards larger specific star formation rates at early times. The age distribution suggests that only ∼3 per cent of bright z ∼ 8.5–11 galaxies would be similarly luminous at z ≳ 15, implying that the number density of bright galaxies declines by at least an order of magnitude between z ∼ 8.5–11 and $z \sim 15$. This evolution is challenging to reconcile with some early JWST results suggesting the abundance of bright galaxies does not significantly decrease towards very early times, but we suggest this tension may be eased if young stellar populations form on top of older stellar components, or if bright z ∼ 15 galaxies are observed during a burst of star formation.

    more » « less
  2. Abstract Accurate models of the star formation histories (SFHs) of recently quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of nonparametric SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs), which have sharp changes in their recent SFH. In this paper, we create mock PSBs, then use the Prospector SED fitting software to test how well four different SFH models recover key properties. We find that a two-component parametric model performs well for our simple mock galaxies, but is sensitive to model mismatches. The fixed- and flexible-bin nonparametric models included in Prospector are able to rapidly quench a major burst of star formation, but systematically underestimate the post-burst age by up to 200 Myr. We develop a custom SFH model that allows for additional flexibility in the recent SFH. Our flexible nonparametric model is able to constrain post-burst ages with no significant offset and just ∼90 Myr of scatter. Our results suggest that while standard nonparametric models are able to recover first-order quantities of the SFH (mass, SFR, average age), accurately recovering higher-order quantities (burst fraction, quenching time) requires careful consideration of model flexibility. These mock recovery tests are a critical part of future SFH studies. Finally, we show that our new, public SFH model is able to accurately recover the properties of mock star-forming and quiescent galaxies and is suitable for broader use in the SED fitting community. 
    more » « less

    Spitzer/Infrared Array Camera (IRAC) imaging has revealed that the brightest z ∼ 7−8 galaxies often exhibit young ages and strong nebular line emission, hinting at high ionizing efficiency among early galaxies. However, IRAC’s limited sensitivity has long hindered efforts to study the fainter, more numerous population often thought largely responsible for reionization. Here, we use Cosmic Evolution Early Release Science (CEERS) JWST/NIRCam data to characterize 116 ultraviolet (UV)-faint (median MUV = −19.5) z ∼ 6.5−8 galaxies. The spectral energy distributions are typically dominated by young (∼10–50 Myr), low-mass (M* ∼ 108 M⊙) stellar populations, and we find no need for extremely high stellar masses (∼1011 M⊙). Considering previous studies of UV-bright (MUV ∼ −22) z ∼ 7−8 galaxies, we find evidence for a strong (5–10 times) increase in specific star formation rate (sSFR) toward lower luminosities (median sSFR = 103 Gyr−1 in CEERS). The larger sSFRs imply a more dominant contribution from OB stars in the relatively numerous UV-faint population, perhaps suggesting that these galaxies are very efficient ionizing agents (median ξion = 1025.7 erg−1 Hz). In spite of the much larger sSFRs, we find little increase in [O iii] + H  β equivalent widths towards fainter MUV (median ≈780 $\mathrm{\mathring{A}}$). If confirmed, this may indicate that a substantial fraction of our CEERS galaxies possess extremely low metallicities (≲3 per cent Z⊙) where [O iii] emission is suppressed. Alternatively, high ionizing photon escape fractions or bursty star formation histories can also weaken the nebular lines in a subset of our sample. While the majority of galaxies in our sample are very blue (median β = −2.0), we identify a significant tail of very dusty galaxies (β ∼ −1) at ≈0.5$L_\mathrm{UV}^\ast$ which may contribute significantly to the z ∼ 7−8 star formation rate density.

    more » « less

    The first deep-field observations of the JWST have immediately yielded a surprisingly large number of very high redshift candidates, pushing the frontier of observability well beyond z ≳ 10. We here present a detailed SED-fitting analysis of the 10 gravitationally lensed z ∼ 9–16 galaxy candidates detected behind the galaxy cluster SMACS J0723.3−7327 in a previous paper using the BEAGLE tool. Our analysis makes use of dynamical considerations to place limits on the ages of these galaxies and of all three published SL models of the cluster to account for lensing systematics. We find the majority of these galaxies to have relatively low stellar masses $M_{\star }\sim 10^7-10^8\, \mathrm{M}_{\odot }$ and young ages tage ∼ 10–100 Myr but with a few higher mass exceptions ($M_{\star }\sim 10^9\rm{-}10^{10}\, \mathrm{M}_{\odot }$) due to Balmer-break detections at z ∼ 9–10. Because of their very blue UV-slopes, down to β ∼ −3, all of the galaxies in our sample have extremely low dust attenuations AV ≲ 0.02. Placing the measured parameters into relation, we find a very shallow M⋆ − MUV-slope and high sSFRs above the main sequence of star formation with no significant redshift-evolution in either relation. This is in agreement with the bright UV luminosities measured for these objects and indicates that we are naturally selecting UV-bright galaxies that are undergoing intense star formation at the time they are observed. Finally, we discuss the robustness of our high-redshift galaxy sample regarding low-redshift interlopers and conclude that low-redshift solutions can safely be ruled out for roughly half of the sample, including the highest redshift galaxies at z ∼ 12–16. These objects represent compelling targets for spectroscopic follow-up observations with JWST and ALMA.

    more » « less
  5. ABSTRACT We present the star formation history (SFH) of the isolated (D ∼ 970 kpc) Local Group dwarf galaxy Wolf–Lundmark–Melotte (WLM) measured from colour–magnitude diagrams (CMDs) constructed from deep Hubble Space Telescope imaging. Our observations include a central ($0.5 \, r_h$) and outer field ($0.7 \, r_h$) that reach below the oldest main-sequence turn-off. WLM has no early dominant episode of star formation: 20 per cent of its stellar mass formed by ∼12.5 Gyr ago ($z$ ∼ 5). It also has an SFR that rises to the present with 50 per cent of the stellar mass within the most recent 5 Gyr ($z$ < 0.7). There is evidence of a strong age gradient: the mean age of the outer field is 5 Gyr older than the inner field despite being only 0.4 kpc apart. Some models suggest such steep gradients are associated with strong stellar feedback and dark-matter core creation. The SFHs of real isolated dwarf galaxies and those from the Feedback in Realistic Environment suite are in good agreement for M⋆($z$ = 0) ∼ 107–109M⊙, but in worse agreement at lower masses ($M_{\star }(z=0) \sim 10^5\!-\!10^7 \, \mathrm{M}_{\odot }$). These differences may be explainable by systematics in the models (e.g. reionization model) and/or observations (HST field placement). We suggest that a coordinated effort to get deep CMDs between HST/JWST (crowded central fields) and WFIRST (wide-area halo coverage) is the optimal path for measuring global SFHs of isolated dwarf galaxies. 
    more » « less