skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Tracing the History of Obscured Star Formation with the SIMBA Cosmological Galaxy Evolution Simulation
Abstract We explore the cosmic evolution of the fraction of dust-obscured star formation predicted by thesimbacosmological hydrodynamic simulations featuring an on-the-fly model for dust formation, evolution, and destruction. We find that up toz= 3, our results are broadly consistent with previous observational results of little to no evolution in obscured star formation. However, atz> 3 we find strong evolution at fixed galaxy stellar mass toward greater amounts of obscured star formation, in tension with high-redshift observations. We explain the trend of increasing obscuration at higher redshifts by evolving star-dust geometry, as the dust-to-stellar mass ratios remain relatively constant across cosmic time. We additionally see that at a fixed redshift, more massive galaxies have a higher fraction of their star formation obscured, which is explained by increased dust-to-stellar mass ratios at higher stellar masses. Finally, we estimate the contribution of dust-obscured star formation to the total star formation rate budget and find that the dust-obscured star formation history peaks aroundz∼ 2−3, and becomes subdominant atz≳ 5. The dominance of obscured star formation at redshiftsz≲ 4 is consistent with our results for the evolution of the obscured star formation fraction at fixed stellar mass to higher values at higher redshift because there exist fewer massive, heavily obscured galaxies at high redshift.  more » « less
Award ID(s):
1909153
PAR ID:
10544678
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
973
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 146
Size(s):
Article No. 146
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Using stellar population synthesis models to infer star formation histories (SFHs), we analyze photometry and spectroscopy of a large sample of quiescent galaxies that are members of Sunyaev–Zel’dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 < z < 1.4 using rest-frame optical spectra and the Python-based Prospector framework, from 61 clusters in the SPT-GMOS Spectroscopic Survey (0.3 < z < 0.9) and three clusters in the SPT Hi-z cluster sample (1.25 < z < 1.4). We analyze spectra of subpopulations divided into bins of redshift, stellar mass, cluster mass, and velocity-radius phase-space location, as well as by creating composite spectra of quiescent member galaxies. We find that quiescent galaxies in our data set sample a diversity of SFHs, with a median formation redshift (corresponding to the lookback time from the redshift of observation to when a galaxy forms 50% of its mass, t 50 ) of z = 2.8 ± 0.5, which is similar to or marginally higher than that of massive quiescent field and cluster galaxy studies. We also report median age–stellar mass relations for the full sample (age of the universe at t 50 (Gyr) = 2.52 (±0.04)–1.66 (±0.12) log 10 ( M /10 11 M ⊙ )) and recover downsizing trends across stellar mass; we find that massive galaxies in our cluster sample form on aggregate ∼0.75 Gyr earlier than lower-mass galaxies. We also find marginally steeper age–mass relations at high redshifts, and report a bigger difference in formation redshifts across stellar mass for fixed environment, relative to formation redshifts across environment for fixed stellar mass. 
    more » « less
  2. Abstract We analyze the evolution of massive (log10[M/M] > 10) galaxies atz∼ 1–4 selected from JWST Cosmic Evolution Early Release Survey (CEERS). We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting withdense basisto select a sample of high-redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam WFSS. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei. We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos forz∼ 1–4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model, resulting in a shallower decline of observed volume densities of massive galaxies. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology. 
    more » « less
  3. Abstract We present results from Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm continuum observations of a sample of 27 star-forming galaxies at z = 2.1–2.5 from the MOSFIRE Deep Evolution Field survey with metallicity and star formation rate measurements from optical emission lines. Using stacks of Spitzer, Herschel, and ALMA photometry (rest frame ∼8–400 μ m), we examine the infrared (IR) spectral energy distributions (SED) of z ∼ 2.3 subsolar-metallicity (∼0.5 Z ⊙ ) luminous infrared galaxies (LIRGs). We find that the data agree well with an average template of higher-luminosity local low-metallicity dwarf galaxies (reduced χ 2 = 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ultraluminous IR galaxies, even in the most favorable case (with reduced χ 2 = 2.8), the templates are rejected at >98% confidence. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar-metallicity galaxies may result from different grain properties or a harder/more intense ionizing radiation field that increases the dust temperature. The obscured star formation rate (SFR) indicated by the far-IR emission of the subsolar-metallicity galaxies is only ∼60% of the total SFR, considerably lower than that of the local LIRGs with ∼96%–97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data overestimate the Rayleigh–Jeans tail measurements by a factor of 2–20. These templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4 dex. At a given stellar mass or metallicity, dust masses at z ∼ 2.3 are an order of magnitude higher than z ∼ 0. Given the predicted molecular gas fractions, the observed z ∼ 2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies with similar metallicities. 
    more » « less
  4. Abstract We make use of sensitive (9.3 μ Jy beam −1 rms) 1.2 mm continuum observations from the Atacama Large Millimeter/submillimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra-Deep Field (ASPECS) large program to probe dust-enshrouded star formation from 1362 Lyman-break galaxies spanning the redshift range z  = 1.5–10 (to ∼7–28 M ⊙ yr −1 at 4 σ over the entire range). We find that the fraction of ALMA-detected galaxies in our z  = 1.5–10 samples increases steeply with stellar mass, with the detection fraction rising from 0% at 10 9.0 M ⊙ to % at >10 10 M ⊙ . Moreover, on stacking all 1253 low-mass (<10 9.25 M ⊙ ) galaxies over the ASPECS footprint, we find a mean continuum flux of −0.1 ± 0.4 μ Jy beam −1 , implying a hard upper limit on the obscured star formation rate of <0.6 M ⊙ yr −1 (4 σ ) in a typical low-mass galaxy. The correlation between the infrared excess (IRX) of UV-selected galaxies ( L IR / L UV ) and the UV-continuum slope is also seen in our ASPECS data and shows consistency with a Calzetti-like relation at > and an SMC-like relation at lower masses. Using stellar mass and β measurements for z  ∼ 2 galaxies over the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, we derive a new empirical relation between β and stellar mass and then use this correlation to show that our IRX– β and IRX–stellar mass relations are consistent with each other. We then use these constraints to express the IRX as a bivariate function of β and stellar mass. Finally, we present updated estimates of star formation rate density determinations at z  > 3, leveraging present improvements in the measured IRX and recent probes of ultraluminous far-IR galaxies at z  > 2. 
    more » « less
  5. Abstract Over the past decade, rest-frame color–color diagrams have become popular tools for selecting quiescent galaxies at high redshift, breaking the color degeneracy between quiescent and dust-reddened star-forming galaxies. In this work, we study one such color–color selection tool—the rest-frameU−VversusV−Jdiagram—by employing mock observations of cosmological galaxy formation simulations. In particular, we conduct numerical experiments assessing both trends in galaxy properties inUVJspace and the color–color evolution of massive galaxies as they quench at redshiftsz∼ 1–2. We find that our models broadly reproduce the observedUVJdiagram atz= 1–2, including (for the first time in a cosmological simulation) reproducing the population of extremely dust-reddened galaxies in the top right of theUVJdiagram. However, our models primarily populate this region with low-mass galaxies and do not produce as clear a bimodality between star-forming and quiescent galaxies as is seen in observations. The former issue is due to an excess of dust in low-mass galaxies and relatively gray attenuation curves in high-mass galaxies, while the latter is due to the overpopulation of the green valley insimba. When investigating the time evolution of galaxies on theUVJdiagram, we find that the quenching pathway on theUVJdiagram is independent of the quenching timescale, and instead dependent primarily on the average specific star formation rate in the 1 Gyr prior to the onset of quenching. Our results support the interpretation of different quenching pathways as corresponding to the divergent evolution of post-starburst and green valley galaxies. 
    more » « less