skip to main content


Title: Bursty Star Formation Naturally Explains the Abundance of Bright Galaxies at Cosmic Dawn
Abstract

Recent discoveries of a significant population of bright galaxies at cosmic dawnz10have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxys’ UV luminosity functions (UVLFs) appear higher than predicted by many models. Using approximately 25,000 galaxy snapshots at 8 ≤z≤ 12 in a suite of FIRE-2 cosmological “zoom-in” simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multichannel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications, such as a nonstandard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright (MUV≲ −20) galaxies atz≳ 10. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models.

 
more » « less
Award ID(s):
1652522 2108230 2307327
NSF-PAR ID:
10484963
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
955
Issue:
2
ISSN:
2041-8205
Format(s):
Medium: X Size: Article No. L35
Size(s):
["Article No. L35"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We study the kinematics of the interstellar medium (ISM) viewed “down the barrel” in 20 gravitationally lensed galaxies during cosmic noon (z= 1.5–3.5). We use moderate-resolution spectra (R∼ 4000) from Keck’s Echellette Spectrograph and Imager and Magellan/MagE to spectrally resolve the ISM absorption in these galaxies into ∼10 independent elements and use double Gaussian fits to quantify the velocity structure of the gas. We find that the bulk motion of gas in this galaxy sample is outflowing, with average velocity centroidvcent=141km s−1(±111 km s−1scatter) measured with respect to the systemic redshift. A total of 16 out of the 20 galaxies exhibit a clear positive skewness, with a blueshifted tail extending to ∼ −500 km s−1. We examine scaling relations in outflow velocities with galaxy stellar mass and star formation rate, finding correlations consistent with a momentum-driven wind scenario. Our measured outflow velocities are also comparable to those reported for FIRE-2 and TNG50 cosmological simulations at similar redshift and galaxy properties. We also consider implications for interpreting results from lower-resolution spectra. We demonstrate that while velocity centroids are accurately recovered, the skewness, velocity width, and probes of high-velocity gas (e.g.,v95) are subject to large scatter and biases at lower resolution. We find thatR≳ 1700 is required for accurate results for the gas kinematics of our sample. This work represents the largest available sample of well-resolved outflow velocity structure atz> 2 and highlights the need for good spectral resolution to recover accurate properties.

     
    more » « less
  2. Abstract

    We explore how the presence of detectable molecular gas depends on the inferred star formation histories (SFHs) in eight massive, quiescent galaxies atz∼ 0.7. Half of the sample have clear detections of molecular gas, traced by CO(2–1). We find that the molecular gas content is unrelated to the rate of star formation decline prior to the most recent 1 Gyr, suggesting that the gas reservoirs are not left over from their primary star formation epoch. However, the recent SFHs of CO-detected galaxies demonstrate evidence for secondary bursts of star formation in their last Gyr. The fraction of stellar mass formed in these secondary bursts ranges fromfburst≈ 0.3%–6% and ended betweentend-burst≈ 0–330 Myr ago. The CO-detected galaxies form a higher fraction of mass in the last Gyr (fM1Gyr=2.6%±1.8%) compared to the CO-undetected galaxies (fM1Gyr=0.2%±0.1%). The galaxies with gas reservoirs have enhanced late-time star formation, highlighting this as a contributing factor to the observed heterogeneity in the gas reservoirs in high-redshift quiescent galaxies. We find that the amount of gas and star formation driven by these secondary bursts are inconsistent with that expected from dry minor mergers, and instead are likely driven by recently accreted gas, i.e., gas-rich minor mergers. This conclusion would not have been made based on SFRUV+IRmeasurements alone, highlighting the power of detailed SFH modeling in the interpretation of gas reservoirs. Larger samples are needed to understand the frequency of low-level rejuvenation among quiescent galaxies at intermediate redshifts, and to what extent this drives the diversity of molecular gas reservoirs.

     
    more » « less
  3. Abstract

    We investigate how cosmic web structures affect galaxy quenching in the IllustrisTNG (TNG100) cosmological simulations by reconstructing the cosmic web within each snapshot using the DisPerSE framework. We measure the comoving distance from each galaxy with stellar masslog(M*/M)8to the nearest node (dnode) and the nearest filament spine (dfil) to study the dependence of both the median specific star formation rate (〈sSFR〉) and the median gas fraction (〈fgas〉) on these distances. We find that the 〈sSFR〉 of galaxies is only dependent on the cosmic web environment atz< 2, with the dependence increasing with time. Atz≤ 0.5,8log(M*/M)<9galaxies are quenched atdnode≲ 1 Mpc, and have significantly suppressed star formation atdfil≲ 1 Mpc, trends driven mostly by satellite galaxies. Atz≤ 1, in contrast to the monotonic drop in 〈sSFR〉 oflog(M*/M)<10galaxies with decreasingdnodeanddfil,log(M*/M)10galaxies—both centrals and satellites—experience an upturn in 〈sSFR〉 atdnode≲ 0.2 Mpc. Much of this cosmic web dependence of star formation activity can be explained by an evolution in 〈fgas〉. Our results suggest that in the past ∼10 Gyr, low-mass satellites are quenched by rapid gas stripping in dense environments near nodes and gradual gas starvation in intermediate-density environments near filaments. At earlier times, cosmic web structures efficiently channeled cold gas into most galaxies. State-of-the-art ongoing spectroscopic surveys such as the Sloan Digital Sky Survey and DESI, as well as those planned with the Subaru Prime Focus Spectrograph, JWST, and Roman, are required to test our predictions against observations.

     
    more » « less
  4. Abstract

    We present a search for extremely red, dust-obscured,z> 7 galaxies with JWST/NIRCam+MIRI imaging over the first 20 arcmin2of publicly available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W−F444W (∼2.5 mag) and detection in MIRI/F770W (∼25 mag), we identify two galaxies, COS-z8M1 and CEERS-z7M1, that have best-fit photometric redshifts ofz=8.40.4+0.3and7.60.1+0.1, respectively. We perform spectral energy distribution fitting with a variety of codes (includingbagpipes,prospector,beagle, andcigale) and find a >95% probability that these indeed lie atz> 7. Both sources are compact (Reff≲ 200 pc) and highly obscured (AV∼ 1.5–2.5) and, at our best-fit redshift estimates, likely have strong [Oiii]+Hβemission contributing to their 4.4μm photometry. We estimate stellar masses of ∼1010Mfor both sources; by virtue of detection in MIRI at 7.7μm, these measurements are robust to the inclusion of bright emission lines, for example, from an active galactic nucleus. We identify a marginal (2.9σ) Atacama Large Millimeter/submillimeter Array detection at 2 mm within 0.″5 of COS-z8M1, which, if real, would suggest a remarkably high IR luminosity of ∼1012L. These two galaxies, if confirmed atz∼ 8, would be extreme in their stellar and dust masses and may be representative of a substantial population of highly dust-obscured galaxies at cosmic dawn.

     
    more » « less
  5. Abstract

    We present a Keck/MOSFIRE rest-optical composite spectrum of 16 typical gravitationally lensed star-forming dwarf galaxies at 1.7 ≲z≲ 2.6 (zmean= 2.30), all chosen independent of emission-line strength. These galaxies have a median stellar mass oflog(M*/M)med=8.290.43+0.51and a median star formation rate ofSFRHαmed=2.251.26+2.15Myr1. We measure the faint electron-temperature-sensitive [Oiii]λ4363 emission line at 2.5σ(4.1σ) significance when considering a bootstrapped (statistical-only) uncertainty spectrum. This yields a direct-method oxygen abundance of12+log(O/H)direct=7.880.22+0.25(0.150.06+0.12Z). We investigate the applicability at highzof locally calibrated oxygen-based strong-line metallicity relations, finding that the local reference calibrations of Bian et al. best reproduce (≲0.12 dex) our composite metallicity at fixed strong-line ratio. At fixedM*, our composite is well represented by thez∼ 2.3 direct-method stellar mass—gas-phase metallicity relation (MZR) of Sanders et al. When comparing to predicted MZRs from the IllustrisTNG and FIRE simulations, having recalculated our stellar masses with more realistic nonparametric star formation histories(log(M*/M)med=8.920.22+0.31), we find excellent agreement with the FIRE MZR. Our composite is consistent with no metallicity evolution, at fixedM*and SFR, of the locally defined fundamental metallicity relation. We measure the doublet ratio [Oii]λ3729/[Oii]λ3726 = 1.56 ± 0.32 (1.51 ± 0.12) and a corresponding electron density ofne=10+215cm3(ne=10+74cm3) when considering the bootstrapped (statistical-only) error spectrum. This result suggests that lower-mass galaxies have lower densities than higher-mass galaxies atz∼ 2.

     
    more » « less