skip to main content

Title: Simulating JWST deep extragalactic imaging surveys and physical parameter recovery
We present a new prospective analysis of deep multi-band imaging with the James Webb Space Telescope (JWST). In this work, we investigate the recovery of high-redshift 5 <   z  <  12 galaxies through extensive image simulations of accepted JWST programs, including the Early Release Science in the EGS field and the Guaranteed Time Observations in the HUDF. We introduced complete samples of ∼300 000 galaxies with stellar masses of log( M * / M ⊙ ) > 6 and redshifts of 0 <   z  <  15, as well as galactic stars, into realistic mock NIRCam, MIRI, and HST images to properly describe the impact of source blending. We extracted the photometry of the detected sources, as in real images, and estimated the physical properties of galaxies through spectral energy distribution fitting. We find that the photometric redshifts are primarily limited by the availability of blue-band and near-infrared medium-band imaging. The stellar masses and star formation rates are recovered within 0.25 and 0.3 dex, respectively, for galaxies with accurate photometric redshifts. Brown dwarfs contaminating the z  >  5 galaxy samples can be reduced to < 0.01 arcmin −2 with a limited impact on galaxy completeness. We investigate multiple high-redshift galaxy selection more » techniques and find that the best compromise between completeness and purity at 5 <   z  <  10 using the full redshift posterior probability distributions. In the EGS field, the galaxy completeness remains higher than 50% at magnitudes m UV  <  27.5 and at all redshifts, and the purity is maintained above 80 and 60% at z  ≤ 7 and 10, respectively. The faint-end slope of the galaxy UV luminosity function is recovered with a precision of 0.1–0.25, and the cosmic star formation rate density within 0.1 dex. We argue in favor of additional observing programs covering larger areas to better constrain the bright end. « less
; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Astronomy & Astrophysics
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Galaxy clustering measurements can be used to constrain many aspects of galaxy evolution, including galaxy host halo masses, satellite quenching efficiencies, and merger rates. We simulate JWST galaxy clustering measurements at z ∼ 4–10 by utilizing mock galaxy samples produced by an empirical model, the universemachine. We also adopt the survey footprints and typical depths of the planned joint NIRCam and NIRSpec Guaranteed Time Observation program planned for Cycle 1 to generate realistic JWST survey realizations and to model high-redshift galaxy selection completeness. We find that galaxy clustering will be measured with ≳5σ significance at z ∼ 4–10. Halo mass precisions resulting frommore »Cycle 1 angular clustering measurements will be ∼0.2 dex for faint (−18 ≳ $\mathit {M}_{\mathrm{UV}}^{ }$ ≳ −19) galaxies at z ∼ 4–10 as well as ∼0.3 dex for bright ($\mathit {M}_{\mathrm{UV}}^{ }$ ∼ −20) galaxies at z ∼ 4–7. Dedicated spectroscopic follow-up over ∼150 arcmin2 would improve these precisions by ∼0.1 dex by removing chance projections and low-redshift contaminants. Future JWST observations will therefore provide the first constraints on the stellar–halo mass relation in the epoch of reionization and substantially clarify how this relation evolves at z > 4. We also find that ∼1000 individual satellites will be identifiable at z ∼ 4–8 with JWST, enabling strong tests of satellite quenching evolution beyond currently available data (z ≲ 2). Finally, we find that JWST observations can measure the evolution of galaxy major merger pair fractions at z ∼ 4–8 with ∼0.1–0.2 dex uncertainties. Such measurements would help determine the relative role of mergers to the build-up of stellar mass into the epoch of reionization.« less
  2. We calibrate and validate different methods of rest-frame color-color selection to identify galaxies in active star-forming and quiescent stages of their evolution. Our method is similar to the widely-used UVJ color-color diagram, which is an effective way to distinguish between quiescent and star-forming galaxies using their rest-frame U-V and V-J colors. UVJ colors suffer known systematics, and at z > 4 the method must be extrapolated because the rest-frame J-band moves beyond the coverage of the deepest bandpasses (typically IRAC 4.5 µm). This leads to biases: for example, spectroscopic campaigns have shown that UVJ-quiescent samples include ~10-30% contamination from galaxiesmore »with significant amounts of star formation. Alternative selection methods will be important not just to mitigate these biases, but also in the JWST era where NIRCam coverage is also limited to ~5 µm . In this poster, we present calibrations of alternative rest-frame filter combinations that are applicable for galaxies at redshifts z = 4 - 6. We apply our method to a stellar mass-limited sample of galaxies at 4 < z < 6 from the FLAMINGOS-2 Extragalactic Near-Infrared K-Split (FENIKS) survey. FENIKS is a deep (23.1 - 24.5 AB mag) survey employing two novel filters which split the Ks band ( λc = 2.2 µm) K-blue and K-red filters ( λc = 1.9 and 2.3 µm, respectively), allowing for finer sampling of the Balmer/4000 Å break of galaxies with evolved populations. We quantify the improvement in the selection of quiescent and star-forming galaxies using the alternative color-color selection methods. Furthermore, we investigate correlations between galaxy properties and their rest-frame colors, in particular examining purity and completeness of these selection methods. Finally, we explore the above for a wide range of synthetic filter combinations to inform accurate selections of various galaxy populations and rule out unphysical areas of parameter space for these populations.« less
  3. ABSTRACT The James Webb Space Telescope (JWST) is expected to observe galaxies at z > 10 that are presently inaccessible. Here, we use a self-consistent empirical model, the universemachine, to generate mock galaxy catalogues and light-cones over the redshift range z = 0−15. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy–halo relationships, and galaxy–galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at z < 10. We predict that Cycle 1 JWST surveys will very likely detect galaxies with M* > 107 M⊙ and/or M1500 < −17 out tomore »at least z ∼ 13.5. Number density uncertainties at z > 12 expand dramatically, so efforts to detect z > 12 galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits (M* < 107 M⊙ or M1500 > −17). For reionization models, extrapolating luminosity functions with a constant faint-end slope from M1500 = −17 down to M1500 = −12 gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to z ∼ 12. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and light-cones for common fields are available online.« less
  4. 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)more »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.« less
  5. 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 Carlomore »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}$.« less