skip to main content

Title: High-redshift JWST predictions from IllustrisTNG: dust modelling and galaxy luminosity functions
ABSTRACT The James Webb Space Telescope (JWST) promises to revolutionize our understanding of the early Universe, and contrasting its upcoming observations with predictions of the Λ cold dark matter model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multiband galaxy luminosity functions from z = 2 to z = 10. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalization and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalization compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust model calibrations. Furthermore, we also recover the observed high-redshift (z = 4–6) UV luminosity versus stellar mass relation, the H α versus star formation rate relation, and the H α luminosity function at z = 2. The bright end more » (MUV > −19.5) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect ∼80 (∼200) galaxies with a signal-to-noise ratio of 10 (5) within the NIRCam field of view, $2.2\times 2.2 \, {\rm arcmin}^{2}$, for a total exposure time of $10^5\, {\rm s}$ in the redshift range z = 8 ± 0.5. These numbers drop to ∼10 (∼40) for an exposure time of $10^4\, {\rm s}$. « less
; ; ; ; ; ; ; ; ; ;
Award ID(s):
1814053 1814259
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
5167 to 5201
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We present predictions for high redshift (z = 2−10) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the H α and H β + $[\rm O \,{\small III}]$ luminosity functions up to z = 8. The predicted H β + $[\rm O \,{\small III}]$ luminosity functions are consistent with present observations at z ≲ 3 with ${\lesssim} 0.1\, {\rm dex}$ differences in luminosities. However, the predicted H α luminosity function is ${\sim }0.3\, {\rm dex}$ dimmer than the observed one at z ≃ 2. Furthermore, we explore continuum spectral indices, the Balmer break at 4000 Å; (D4000) and the UV continuum slope β. The median D4000 versus specific star formation rate relation predicted at z = 2 is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed AUV versus β relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at z = 2−6 and investigate their dependence on galaxy colours andmore »stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colours, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high-redshift galaxies. Future James Webb Space Telecope observations will further test these predictions.« less
  2. ABSTRACT We post-process galaxies in the IllustrisTNG simulations with skirt radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at z ≥ 4. The rest-frame K- and z-band galaxy luminosity functions from TNG are overall consistent with observations, despite ${\sim}0.5\, \mathrm{dex}$ underprediction at z = 4 for MK ≲ −25 and Mz ≲ −24. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. Based on theoretical estimations, we show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $500\, {\rm arcmin}^{2}$ at z = 6 (z = 8). As opposed to the consistency in the UV, optical, and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshiftmore »are qualitatively consistent with observations, the peak dust temperature of z ≥ 6 galaxies are overestimated by about $20\, {\rm K}$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts.« less
  3. Abstract We investigate the relationship between dust attenuation and stellar mass ( M * ) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field survey at z ∼ 2.3 and the Sloan Digital Sky Survey (SDSS) at z ∼ 0, augmenting the latter optical data set with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of H α and H β emission lines, and photometric measurements of the rest-UV stellar continuum. The H α /H β ratio is used to determine the magnitude of attenuation at the wavelength of H α , A H α . Rest-UV colors and spectral energy distribution fitting are used to estimate A 1600 , the magnitude of attenuation at a rest wavelength of 1600 Å. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and M * over the redshift range z ∼ 0 to z ∼ 2.3. Folding in the latest estimates of the evolution of M dust , ( M dust / M gas ), and gasmore »surface density at fixed M * , we find that the expected M dust and dust mass surface density are both significantly higher at z ∼ 2.3 than at z ∼ 0. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation versus M * , it is essential to determine the relationship between metallicity and ( M dust / M gas ), the dust mass absorption coefficient and dust geometry, and the evolution of these relations and quantities from z ∼ 0 to z ∼ 2.3.« less
  4. Abstract We present Ly α and ultraviolet (UV)-continuum luminosity functions (LFs) of galaxies and active galactic nuclei (AGNs) at z = 2.0–3.5 determined by the untargeted optical spectroscopic survey of the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in 11.4 deg 2 of fiber spectra sky coverage, obtaining 18,320 galaxies spectroscopically identified with Ly α emission, 2126 of which host type 1 AGNs showing broad (FWHM > 1000 km s −1 ) Ly α emission lines. We derive the Ly α (UV) LF over 2 orders of magnitude covering bright galaxies and AGNs in log L Ly α / [ erg s − 1 ] = 43.3 – 45.5 (−27 < M UV < −20) by the 1/ V max estimator. Our results reveal that the bright-end hump of the Ly α LF is composed of type 1 AGNs. In conjunction with previous spectroscopic results at the faint end, we measure a slope of the best-fit Schechter function to be α Sch = − 1.70 − 0.14 + 0.13 , which indicates that α Sch steepens from z = 2–3 toward high redshift. Our UV LF agrees well with previous AGN UVmore »LFs and extends to faint-AGN and bright-galaxy regimes. The number fraction of Ly α -emitting objects ( X LAE ) increases from M UV * ∼ − 21 to bright magnitude due to the contribution of type 1 AGNs, while previous studies claim that X Ly α decreases from faint magnitudes to M UV * , suggesting a valley in the X Ly α –magnitude relation at M UV * . Comparing our UV LF of type 1 AGNs at z = 2–3 with those at z = 0, we find that the number density of faint ( M UV > −21) type 1 AGNs increases from z ∼ 2 to 0, as opposed to the evolution of bright ( M UV < −21) type 1 AGNs, suggesting AGN downsizing in the rest-frame UV luminosity.« less
  5. Star formation rate (SFR) measurements at z  > 4 have relied mostly on the rest-frame far-ultraviolet (FUV) observations. The corrections for dust attenuation based on the IRX- β relation are highly uncertain and are still debated in the literature. Hence, rest-frame far-infrared (FIR) observations are necessary to constrain the dust-obscured component of the SFR. In this paper, we exploit the rest-frame FIR continuum observations collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) to directly constrain the obscured SFR in galaxies at 4.4 <  z  < 5.9. We used stacks of continuum images to measure average infrared luminosities taking both detected and undetected sources into account. Based on these measurements, we measured the position of the main sequence of star-forming galaxies and the specific SFR (sSFR) at z  ∼ 4.5 and z  ∼ 5.5. We find that the main sequence and sSFR do not significantly evolve between z  ∼ 4.5 and z  ∼ 5.5, as opposed to lower redshifts. We developed a method to derive the obscured SFR density (SFRD) using the stellar masses or FUV-magnitudes as a proxy of FIR fluxes measured on the stacks and combining them with the galaxy stellar mass functions and FUV luminosity functions from the literature. Wemore »obtain consistent results independent of the chosen proxy. We find that the obscured fraction of SFRD is decreasing with increasing redshift, but even at z  ∼ 5.5 it constitutes around 61% of the total SFRD.« less