We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$ at z = 0) as part of the ‘Feedback in Realistic Environment’ (FIRE) project. Compact (r1/2 < 30 pc), relatively massive (0.5 × 105 ≲ M⋆/M⊙ ≲ 5 × 105), self-bound stellar clusters form at 11 ≳ z ≳ 5 in progenitors with $M_{\rm vir} \approx 10^9\, \mathrm{M}_{\odot }$. Cluster formation is triggered when at least $10^7\, \mathrm{M}_{\odot }$ of dense, turbulent gas reaches $\Sigma _{\rm gas} \approx 10^4\, {\rm M}_\odot \, {\rm pc}^{-2}$ as a result of the compressive effects of supernova feedback or from cloud–cloud collisions. The clusters can survive for $2-3\, {\rm Gyr}$; absent numerical effects, they could possibly survive substantially longer, perhaps to z = 0. The longest lived clusters are those that form at significant distance – several hundreds of pc – from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter haloes as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in Hubble Space Telescope lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with JWST. By contrast, the star clusters’ host galaxies will remain undetectable.
We use two independent galaxy-formation simulations, flares, a cosmological hydrodynamical simulation, and shark, a semi-analytic model, to explore how well the JWST will be able to uncover the existence and parameters of the star-forming main sequence (SFS) at z = 5 → 10, i.e. shape, scatter, normalization. Using two independent simulations allows us to isolate predictions (e.g. stellar mass, star formation rate, SFR, luminosity functions) that are robust to or highly dependent on the implementation of the physics of galaxy formation. Both simulations predict that JWST can observe ≥70–90 per cent (for shark and flares, respectively) of galaxies up to z ∼ 10 (down to stellar masses of ${\approx}10^{8.3}\rm M_{\odot }$ and SFRs of ${\approx}10^{0.5}{\rm M}_{\odot }\,{\rm yr}^{-1}$) in modest integration times and given current proposed survey areas (e.g. the Web COSMOS 0.6 deg2) to accurately constrain the parameters of the SFS. Although both simulations predict qualitatively similar distributions of stellar mass and SFR. There are important quantitative differences, such as the abundance of massive, star-forming galaxies with flares predicting a higher abundance than shark; the early onset of quenching as a result of black hole growth in flares (at z ≈ 8), not seen in shark until much lower redshifts; and the implementation of synthetic photometry with flares predicting more JWST-detected galaxies (∼90 per cent) than shark (∼70 per cent) at z = 10. JWST observations will distinguish between these models, leading to a significant improvement upon our understanding of the formation of the very first galaxies.
more » « less- NSF-PAR ID:
- 10380445
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 518
- Issue:
- 1
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 456-476
- Size(s):
- ["p. 456-476"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT Both observations and simulations have shown strong evidence for highly time-variable star formation in low-mass and/or high-redshift galaxies, which has important observational implications because high-redshift galaxy samples are rest-ultraviolet (rest-UV) selected and therefore particularly sensitive to the recent star formation. Using a suite of cosmological ‘zoom-in’ simulations at z > 5 from the Feedback in Realistic Environments project, we examine the implications of bursty star formation histories for observations of high-redshift galaxies with JWST. We characterize how the galaxy observability depends on the star formation history. We also investigate selection effects due to bursty star formation on the physical properties measured, such as the gas fraction, specific star formation rate, and metallicity. We find the observability to be highly time-dependent for galaxies near the survey’s limiting flux due to the star formation rate variability: as the star formation rate fluctuates, the same galaxy oscillates in and out of the observable sample. The observable fraction $f_\mathrm{obs} = 50~{{\ \rm per\ cent}}$ at z ∼ 7 and M⋆ ∼ 108.5–$10^{9}\, {\rm M}_{\odot }$ for a JWST/NIRCam survey reaching a limiting magnitude of $m^\mathrm{lim}_\mathrm{AB} \sim 29{\!-\!}30$, representative of surveys such as JADES and CEERS. JWST-detectable galaxies near the survey limit tend to have properties characteristic of galaxies in the bursty phase: on average, they show approximately 2.5 times higher cold, dense gas fractions and 20 times higher specific star formation rates at a given stellar mass than galaxies below the rest-UV detection threshold. Our study represents a first step in quantifying selection effects and the associated biases due to bursty star formation in studying high-redshift galaxy properties.
-
more » « less
ABSTRACT We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\lesssim$4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $\gt t_{\rm H}^{-1}$ to $\lt 0.2t_{\rm H}^{-1}$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\lesssim 1{{\ \rm per\ cent}}$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the H2 content, but for $M_*\gtrsim 10^{10.5} \,\mathrm{ M}_{\odot }$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift, ∝(1 + z)3.5, but the quenching rate evolves much more slowly, ∝(1 + z)0.9; there are insufficient mergers to explain the quenching rate at $z\lesssim 1.5$. simba first quenches galaxies at $z\gtrsim 3$, with a number density in good agreement with observations. The quenching time-scales τq are strongly bimodal, with ‘slow’ quenchings (τq ∼ 0.1tH) dominating overall, but ‘fast’ quenchings (τq ∼ 0.01tH) dominating in M* ∼ 1010–1010.5 M$\odot$ galaxies, likely induced by simba’s jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.
-
ABSTRACT We measure the rate of environmentally driven star formation quenching in galaxies at z ∼ 1, using eleven massive ($M\approx 2\times 10^{14}\, \mathrm{M}_\odot$) galaxy clusters spanning a redshift range 1.0 < z < 1.4 from the GOGREEN sample. We identify three different types of transition galaxies: ‘green valley’ (GV) galaxies identified from their rest-frame (NUV − V) and (V − J) colours; ‘blue quiescent’ (BQ) galaxies, found at the blue end of the quiescent sequence in (U − V) and (V − J) colour; and spectroscopic post-starburst (PSB) galaxies. We measure the abundance of these galaxies as a function of stellar mass and environment. For high-stellar mass galaxies (log M/M⊙ > 10.5) we do not find any significant excess of transition galaxies in clusters, relative to a comparison field sample at the same redshift. It is likely that such galaxies were quenched prior to their accretion in the cluster, in group, filament, or protocluster environments. For lower stellar mass galaxies (9.5 < log M/M⊙ < 10.5) there is a small but significant excess of transition galaxies in clusters, accounting for an additional ∼5–10 per cent of the population compared with the field. We show that our data are consistent with a scenario in which 20–30 per cent of low-mass, star-forming galaxies in clusters are environmentally quenched every Gyr, and that this rate slowly declines from z = 1 to z = 0. While environmental quenching of these galaxies may include a long delay time during which star formation declines slowly, in most cases this must end with a rapid (τ < 1 Gyr) decline in star formation rate.more » « less
-
more » « less
ABSTRACT 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.
