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

    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

    The ionizing photon escape fraction [Lyman continuum (LyC) fesc] of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshifts. Here we circumvent this uncertainty by inferring fesc using resolved (R > 4000) Lyman α (Lyα) profiles from the X-SHOOTER Lyα survey at z = 2 (XLS-z2). With empirically motivated criteria, we use Lyα profiles to select leakers ($f_{\mathrm{ esc}} > 20{{\ \rm per\ cent}}$) and non-leakers ($f_{\mathrm{ esc}} < 5{{\ \rm per\ cent}}$) from a representative sample of >0.2L* Lyman α emitters (LAEs). We use median stacked spectra of these subsets over λrest ≈ 1000–8000 Å to investigate the conditions for LyC fesc. Our stacks show similar mass, metallicity, MUV, and βUV. We find the following differences between leakers versus non-leakers: (i) strong nebular C iv and He ii emission versus non-detections; (ii) [O iii]/[O ii] ≈ 8.5 versus ≈3; (iii) Hα/Hβ indicating no dust versus E(B − V) ≈ 0.3; (iv) Mg ii emission close to the systemic velocity versus redshifted, optically thick Mg ii; and (v) Lyα fesc of ${\approx} 50{{\ \rm per\ cent}}$ versus ${\approx} 10{{\ \rm per\ cent}}$. The extreme equivalent widths (EWs) in leakers ([O iii]+$\mathrm{ H}\beta \approx 1100$ Å rest frame)more »constrain the characteristic time-scale of LyC escape to ≈3–10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers – extremely ionizing stellar populations, low column densities, a dust-free, high-ionization state interstellar medium (ISM) – occur simultaneously in the $f_{\rm esc} > 20{{\ \rm per\ cent}}$ stack, suggesting they are causally connected, and motivating why indicators like [O iii]/[O ii] may suffice to constrain fesc at z > 6 with the James Webb Space Telescope (JWST). The leakers comprise half of our sample, have a median LyC$f_{\rm esc} \approx 50{{\ \rm per\ cent}}$ (conservative range: $20\!-\!55{{\ \rm per\ cent}}$), and an ionizing production efficiency $\log ({\xi _{\rm {ion}}/\rm {Hz\ erg^{-1}}})\approx 25.9$ (conservative range: 25.7–25.9). These results show LAEs – the type of galaxies rare at z ≈ 2, but that become the norm at higher redshift – are highly efficient ionizers, with extreme ξion and prolific fesc occurring in sync.

    « less
  4. 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 atz∼ 2.3 and the Sloan Digital Sky Survey (SDSS) atz∼ 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α,AHα. Rest-UV colors and spectral energy distribution fitting are used to estimateA1600, 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 andM*over the redshift rangez∼ 0 toz∼ 2.3. Folding in the latest estimates of the evolution ofMdust, (Mdust/Mgas), and gas surface density at fixedM*, we find that the expectedMdustand dust mass surface density are both significantly higher atz∼ 2.3 than atz∼ 0. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancymore »in attenuation versusM*, it is essential to determine the relationship between metallicity and (Mdust/Mgas), the dust mass absorption coefficient and dust geometry, and the evolution of these relations and quantities fromz∼ 0 toz∼ 2.3.

    « less
  5. 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