scholarly journals High-redshift JWST predictions from IllustrisTNG: dust modelling and galaxy luminosity functions

2020 ◽  
Vol 492 (4) ◽  
pp. 5167-5201 ◽  
Author(s):  
Mark Vogelsberger ◽  
Dylan Nelson ◽  
Annalisa Pillepich ◽  
Xuejian Shen ◽  
Federico Marinacci ◽  
...  

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 (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}$.

2020 ◽  
Vol 495 (4) ◽  
pp. 4747-4768 ◽  
Author(s):  
Xuejian Shen ◽  
Mark Vogelsberger ◽  
Dylan Nelson ◽  
Annalisa Pillepich ◽  
Sandro Tacchella ◽  
...  

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 and 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.


2020 ◽  
Vol 493 (3) ◽  
pp. 3966-3984 ◽  
Author(s):  
A A Khostovan ◽  
S Malhotra ◽  
J E Rhoads ◽  
C Jiang ◽  
J Wang ◽  
...  

ABSTRACT We present our measurements of the H α, [O iii], and [O ii] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 z = 0.47 H α-, 3933 z = 0.93 [O iii]-, and 5367 z = 1.59 [O ii]-selected emission line galaxies in a 3 deg2 single, CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5σ depths of 8.2 × 10−18 erg s−1 cm−2 and comoving volumes of (1−7) × 105 Mpc3 making our survey one of the deepest narrow-band surveys. We select our emission line galaxies via spectroscopic confirmation, photometric redshifts, and colour–colour selections. We measure the observed luminosity functions for each sample and find best fits of $\phi ^\star = 10^{-3.16^{+0.09}_{-0.09}}$ Mpc−3 and $L^\star = 10^{41.72^{+0.09}_{-0.09}}$ erg s−1 for H α, $\phi ^\star = 10^{-2.16^{+0.10}_{-0.12}}$ Mpc−3 and $L^\star = 10^{41.38^{+0.07}_{-0.06}}$ erg s−1 for [O iii], and $\phi ^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc−3 and $L^\star = 10^{41.66^{+0.03}_{-0.03}}$ erg s−1 for [O ii], with α fixed to −1.75, −1.6, and −1.3, respectively. An excess of bright >1042 erg s−1 [O iii] emitters is observed and may be due to active galactic nucleus (AGN) contamination. Corrections for dust attenuation are applied assuming AHα = 1 mag. We also design our own empirical rest-frame g − r calibration using SDSS DR12 data, test it against our z = 0.47 H α emitters with zCOSMOS 1D spectra, and calibrate it for (g − r) between −0.8 and 1.3 mag. Dust and AGN-corrected star formation rate densities (SFRDs) are measured as log10ρSFR/(M⊙ yr−1 Mpc−3) = −1.63 ± 0.04, −1.07 ± 0.06, and −0.90 ± 0.10 for H α, [O iii], and [O ii], respectively. We find our [O iii] and [O ii] samples fully trace cosmic star formation activity at their respective redshifts in comparison to multiwavelength SFRDs, while the H α sample traces ∼70 per cent of the total z = 0.47 SFRD.


2019 ◽  
Vol 15 (S352) ◽  
pp. 64-68
Author(s):  
Xiangcheng Ma

AbstractWe present a suite of high-resolution cosmological zoom-in simulations of galaxies at z⩾ 5using the state-of-the-art models for the multi-phase ISM, star formation, and stellar feedback from the FIRE project. We present a series of key results from these simulations, including the stellar mass–halo mass relation, the ultraviolet luminosity functions, dust attenuation and dust temperatures, the ubiquitous formation of bound star clusters, morphology and clumpiness, and the escape fractions of ionizing photons from high-redshift galaxies. We discuss how different simulations in the literature agree and disagree and what observations are most useful for testing the models in the era of ALMA and JWST.


2019 ◽  
Vol 489 (1) ◽  
pp. 1357-1372 ◽  
Author(s):  
Yisheng Qiu ◽  
Simon J Mutch ◽  
Elisabete da Cunha ◽  
Gregory B Poole ◽  
J Stuart B Wyithe

ABSTRACT We present a new analysis of high-redshift UV observations using a semi-analytic galaxy formation model, and provide self-consistent predictions of the infrared excess (IRX)–β relations and cosmic star formation rate density. We combine the Charlot & Fall dust attenuation model with the meraxes semi-analytic model, and explore three different parametrizations for the dust optical depths, linked to star formation rate, dust-to-gas ratio, and gas column density, respectively. A Bayesian approach is employed to statistically calibrate model-free parameters including star formation efficiency, mass loading factor, dust optical depths, and reddening slope directly against UV luminosity functions and colour–magnitude relations at $z$ ∼ 4–7. The best-fitting models show excellent agreement with the observations. We calculate IRX using energy balance arguments and find that the large intrinsic scatter in the IRX–β plane correlates with specific star formation rate. Additionally, the difference among the three dust models suggests at least a factor of 2 systematic uncertainty in the dust-corrected star formation rate when using the Meurer IRX–β relation at $z$ ≳ 4.


2020 ◽  
Vol 500 (3) ◽  
pp. 3394-3412
Author(s):  
Steven R Furlanetto

ABSTRACT In recent years, simple models of galaxy formation have been shown to provide reasonably good matches to available data on high-redshift luminosity functions. However, these prescriptions are primarily phenomenological, with only crude connections to the physics of galaxy evolution. Here, we introduce a set of galaxy models that are based on a simple physical framework but incorporate more sophisticated models of feedback, star formation, and other processes. We apply these models to the high-redshift regime, showing that most of the generic predictions of the simplest models remain valid. In particular, the stellar mass–halo mass relation depends almost entirely on the physics of feedback (and is thus independent of the details of small-scale star formation) and the specific star formation rate is a simple multiple of the cosmological accretion rate. We also show that, in contrast, the galaxy’s gas mass is sensitive to the physics of star formation, although the inclusion of feedback-driven star formation laws significantly changes the naive expectations. While these models are far from detailed enough to describe every aspect of galaxy formation, they inform our understanding of galaxy formation by illustrating several generic aspects of that process, and they provide a physically grounded basis for extrapolating predictions to faint galaxies and high redshifts currently out of reach of observations. If observations show violations from these simple trends, they would indicate new physics occurring inside the earliest generations of galaxies.


2020 ◽  
Vol 500 (2) ◽  
pp. 2127-2145
Author(s):  
Christopher C Lovell ◽  
Aswin P Vijayan ◽  
Peter A Thomas ◽  
Stephen M Wilkins ◽  
David J Barnes ◽  
...  

ABSTRACT We introduce the First Light And Reionisation Epoch Simulations (FLARES), a suite of zoom simulations using the EAGLE model. We resimulate a range of overdensities during the Epoch of Reionization (EoR) in order to build composite distribution functions, as well as explore the environmental dependence of galaxy formation and evolution during this critical period of galaxy assembly. The regions are selected from a large $(3.2 \, \mathrm{cGpc})^{3}$ parent volume, based on their overdensity within a sphere of radius 14 h−1 cMpc. We then resimulate with full hydrodynamics, and employ a novel weighting scheme that allows the construction of composite distribution functions that are representative of the full parent volume. This significantly extends the dynamic range compared to smaller volume periodic simulations. We present an analysis of the galaxy stellar mass function (GSMF), the star formation rate distribution function (SFRF), and the star-forming sequence (SFS) predicted by FLARES, and compare to a number of observational and model constraints. We also analyse the environmental dependence over an unprecedented range of overdensity. Both the GSMF and the SFRF exhibit a clear double-Schechter form, up to the highest redshifts (z = 10). We also find no environmental dependence of the SFS normalization. The increased dynamic range probed by FLARES will allow us to make predictions for a number of large area surveys that will probe the EoR in coming years, carried out on new observatories such as Roman and Euclid.


2019 ◽  
Vol 15 (S341) ◽  
pp. 226-230
Author(s):  
Christian Binggeli ◽  
Erik Zackrisson ◽  
Xiangcheng Ma ◽  
Akio K. Inoue ◽  
Anton Vikaeus ◽  
...  

AbstractRecently, spectroscopic detections of O[III] 88 μm and Ly-α emission lines from the z ≍ 9.1 galaxy MACS1149-JD1 have been presented, and with these, some interesting properties of this galaxy were uncovered. One such property is that MACS1149-JD1 exhibits a significant Balmer break at around rest-frame 4000 Å, which may indicate that the galaxy has experienced large variations in star formation rate prior to z ∼ 9, with a rather long period of low star formation activity. While some simulations predict large variations in star formation activity in high-redshift galaxies, it is unclear whether the simulations can reproduce the kind of variations seen in MACS1149-JD1. Here, we utilize synthetic spectra of simulated galaxies from two simulation suites in order to study to what extent these can accurately reproduce the spectral features (specifically the Balmer break) observed in MACS1149-JD1. We show that while the simulations used in this study produce galaxies with varying star formation histories, galaxies such as MACS1149-JD1 would be very rare in the simulations. In principle, future observations with the James Webb Space Telescope may tell us if MACS1149-JD1 represents something rare, or if such galaxies are more common than predicted by current simulations.


2012 ◽  
Vol 10 (H16) ◽  
pp. 128-128
Author(s):  
Jamie R. Ownsworth ◽  
Christopher J. Conselice ◽  
Alice Mortlock ◽  
William G. Hartley ◽  
Fernando Buitrago

We investigate the resolved star formation properties of a sample of 45 massive galaxies (M* > 1011 M⊙) within a redshift range of 1.5 ⩽ z ⩽ 3 detected in the GOODS NICMOS Survey (Conselice et al. 2011), a HST H160-band imaging program. We derive the star formation rate as a function of radius using rest frame UV data from deep z850 ACS imaging. The star formation present at high redshift is then extrapolated to z = 0, and we examine the stellar mass produced in individual regions within each galaxy. We also construct new stellar mass profiles of the in situ stellar mass at high redshift from Sérsic fits to rest-frame optical, H160-band, data. We combine the two stellar mass profiles to produce an evolved stellar mass profile. We then fit a new Sérsic profile to the evolved profile, from which we examine what effect the resulting stellar mass distribution added via star formation has on the structure and size of each individual galaxy.


2019 ◽  
Vol 491 (3) ◽  
pp. 3891-3899 ◽  
Author(s):  
Jaehong Park ◽  
Nicolas Gillet ◽  
Andrei Mesinger ◽  
Bradley Greig

ABSTRACT Upcoming observations will probe the first billion years of our Universe in unprecedented detail. Foremost among these are 21-cm interferometry with the Hydrogen Epoch of Reionization Arrays (HERA) and the Square Kilometre Array (SKA), and high-z galaxy observations with the James Webb Space Telescope (JWST). Here, we quantify how observations from these instruments can be used to constrain the astrophysics of high-z galaxies. We generate several mock JWST luminosity functions (LFs) and SKA1 21-cm power spectra, which are consistent with current observations, but assume different properties for the unseen, ultrafaint galaxies driving the epoch of reionization (EoR). Using only JWST data, we predict up to a factor of 2–3 improvement (compared with Hubble Space Telescope, HST) in the fractional uncertainty of the star formation rate to halo mass relation and the turnover magnitude. Most parameters regulating the ultraviolet (UV) galaxy properties can be constrained at the level of ∼10 per cent or better, if either (i) we are able to better characterize systematic lensing uncertainties than currently possible; or (ii) the intrinsic LFs peak at magnitudes brighter than MUV ≲ −13. Otherwise, improvement over HST-based inference is modest. When combining with upcoming 21-cm observations, we are able to significantly mitigate degeneracies, and constrain all of our astrophysical parameters, even for our most pessimistic assumptions about upcoming JWST LFs. The 21-cm observations also result in an order of magnitude improvement in constraints on the EoR history.


2020 ◽  
Vol 634 ◽  
pp. A97 ◽  
Author(s):  
Y. Khusanova ◽  
O. Le Fèvre ◽  
P. Cassata ◽  
O. Cucciati ◽  
B. C. Lemaux ◽  
...  

Context. The star formation rate density (SFRD) evolution presents an area of great interest in the studies of galaxy evolution and reionization. The current constraints of SFRD at z >  5 are based on the rest-frame UV luminosity functions with the data from photometric surveys. The VIMOS UltraDeep Survey (VUDS) was designed to observe galaxies at redshifts up to ∼6 and opened a window for measuring SFRD at z >  5 from a spectroscopic sample with a well-controlled selection function. Aims. We establish a robust statistical description of the star-forming galaxy population at the end of cosmic HI reionization (5.0 ≤ z ≤ 6.6) from a large sample of 49 galaxies with spectroscopically confirmed redshifts. We determine the rest-frame UV and Lyα luminosity functions and use them to calculate SFRD at the median redshift of our sample z = 5.6. Methods. We selected a sample of galaxies at 5.0 ≤ zspec ≤ 6.6 from the VUDS. We cleaned our sample from low redshift interlopers using ancillary photometric data. We identified galaxies with Lyα either in absorption or in emission, at variance with most spectroscopic samples in the literature where Lyα emitters (LAE) dominate. We determined luminosity functions using the 1/Vmax method. Results. The galaxies in this redshift range exhibit a large range in their properties. A fraction of our sample shows strong Lyα emission, while another fraction shows Lyα in absorption. UV-continuum slopes vary with luminosity, with a large dispersion. We find that star-forming galaxies at these redshifts are distributed along the main sequence in the stellar mass vs. SFR plane, described with a slope α = 0.85 ± 0.05. We report a flat evolution of the specific SFR compared to lower redshift measurements. We find that the UV luminosity function is best reproduced by a double power law, while a fit with a Schechter function is only marginally inferior. The Lyα luminosity function is best fitted with a Schechter function. We derive a logSFRDUV(M⊙ yr−1 Mpc−3) = −1.45+0.06−0.08 and logSFRDLyα(M⊙ yr−1 Mpc−3) = −1.40+0.07−0.08. The SFRD derived from the Lyα luminosity function is in excellent agreement with the UV-derived SFRD after correcting for IGM absorption. Conclusions. Our new SFRD measurements at a mean redshift of z = 5.6 are ∼0.2 dex above the mean SFRD reported in Madau & Dickinson (2014, ARA&A, 52, 415), but in excellent agreement with results from Bouwens et al. (2015a, ApJ, 803, 34). These measurements confirm the steep decline of the SFRD at z >  2. The bright end of the Lyα luminosity function has a high number density, indicating a significant star formation activity concentrated in the brightest LAE at these redshifts. LAE with equivalent width EW > 25 Å contribute to about 75% of the total UV-derived SFRD. While our analysis favors low dust content in 5.0 <  z <  6.6, uncertainties on the dust extinction correction and associated degeneracy in spectral fitting will remain an issue, when estimating the total SFRD until future surveys extending spectroscopy to the NIR rest-frame spectral domain, such as with JWST.


Sign in / Sign up

Export Citation Format

Share Document