skip to main content

Title: A framework for simultaneously measuring field densities and the high-z luminosity function

Cosmic variance from large-scale structure will be a major source of uncertainty for galaxy surveys at $z \gtrsim 6$, but that same structure will also provide an opportunity to identify and study dense environments in the early Universe. Using a robust model for galaxy clustering, we directly incorporate large-scale densities into an inference framework that simultaneously measures the high-z ($z \gtrsim 6$) UV luminosity function and the average matter density of each distinct volume in a survey. Through this framework, we forecast the performance of several major upcoming James Webb Space Telescope (JWST) galaxy surveys. We find that they can constrain field matter densities down to the theoretical limit imposed by Poisson noise and unambiguously identify over-dense (and under-dense) regions on transverse scales of tens of comoving Mpc. We also predict JWST will measure the luminosity function with a precision at z = 12 comparable to existing Hubble Space Telescope’s constraints at z = 8 (and even better for the faint-end slope). We also find that wide-field surveys are especially important in distinguishing luminosity function models.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 4844-4856
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    One of the most exciting advances of the current generation of telescopes has been the detection of galaxies during the epoch of reionization, using deep fields that have pushed these instruments to their limits. It is essential to optimize our analyses of these fields in order to extract as much information as possible from them. In particular, standard methods of measuring the galaxy luminosity function discard information on large-scale dark matter density fluctuations, even though this large-scale structure drives galaxy formation and reionization during the Cosmic Dawn. Measuring these densities would provide a bedrock observable, connecting galaxy surveys to theoretical models of the reionization process and structure formation. Here, we use existing Hubble deep field data to simultaneously fit the universal luminosity function and measure large-scale densities for each Hubble deep field at z = 6–8 by directly incorporating priors on the large-scale density field and galaxy bias. Our fit of the universal luminosity function is consistent with previous methods but differs in the details. For the first time, we measure the underlying densities of the survey fields, including the most over/underdense Hubble fields. We show that the distribution of densities is consistent with current predictions for cosmic variance. This analysis on just 17 fields is a small sample of what will be possible with the James Webb Space Telescope, which will measure hundreds of fields at comparable (or better) depths and at higher redshifts.

    more » « less
  2. 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}$. 
    more » « less

    One of the key processes driving galaxy evolution during the Cosmic Dawn is supernova feedback. This likely helps regulate star formation inside of galaxies, but it can also drive winds that influence the large-scale intergalactic medium. Here, we present a simple semi-analytic model of supernova-driven galactic winds and explore the contributions of different phases of galaxy evolution to cosmic metal enrichment in the high-redshift (z ≳ 6) Universe. We show that models calibrated to the observed galaxy luminosity function at z ∼ 6–8 have filling factors $\sim 1{{\%}}$ at z ∼ 6 and $\sim 0.1{{\%}}$ at z ∼ 12, with different star formation prescriptions providing about an order of magnitude uncertainty. Despite the small fraction of space filled by winds, these scenarios predict an upper limit to the abundance of metal-line absorbers in quasar spectra at $z \gtrsim 5$ which is comfortably above that currently observed. We also consider enrichment through winds driven by Pop III star formation in minihalos. We find that these can dominate the total filling factor at $z \gtrsim 10$ and even compete with winds from normal galaxies at z ∼ 6, at least in terms of the total enriched volume. But these regions have much lower overall metallicities, because each one is generated by a small burst of star formation. Finally, we show that Compton cooling of these supernova-driven winds at $z \gtrsim 6$ has only a small effect on the cosmic microwave background.

    more » « less

    The high-redshift galaxy UV luminosity function (UVLF) has become essential for understanding the formation and evolution of the first galaxies. Yet, UVLFs only measure galaxy abundances, giving rise to a degeneracy between the mean galaxy luminosity and its stochasticity. Here, we show that upcoming clustering measurements with the JWST, as well as with Roman, will be able to break this degeneracy, even at redshifts z ≳ 10. First, we demonstrate that current Subaru Hyper Suprime-Cam (HSC) measurements of the galaxy bias at z ∼ 4–6 point to a relatively tight halo-galaxy connection, with low stochasticity. Then, we show that the larger UVLFs observed by JWST at z ≳ 10 can be explained with either a boosted average UV emission or an enhanced stochasticity. These two models, however, predict different galaxy biases, which are potentially distinguishable in JWST and Roman surveys. Galaxy-clustering measurements, therefore, will provide crucial insights into the connection between the first galaxies and their dark-matter haloes, and identify the root cause of the enhanced abundance of z ≳ 10 galaxies revealed with JWST during its first year of operations.

    more » « less

    We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z$ ≳ 6 quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in the IllustrisTNG simulations. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim 10^{12}\, {\rm M}_{\odot }\,h^{-1}$ haloes by $z$ = 7) using a  (9 Mpc h−1)3 simulated volume. BH growth is maximal within haloes that are compact and have a low tidal field. For these haloes, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4\!-\!10^6\, {\rm M}_{\odot }\,h^{-1}$ seeds are inserted in haloes above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of Mseed). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z$ ∼ 6 quasar with $\sim 10^9\, {\rm M}_{\odot }$ mass and ∼1047 erg s−1 luminosity. BH mergers play a crucial role at $z$ ≳ 9, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. With more stringent seeding conditions (e.g. $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z$ ≳ 6 quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors ≳10) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results demonstrate that progenitors of $z$ ∼ 6 quasars have distinct BH merger histories for different seeding models, which will be distinguishable with Laser Interferometer Space Antenna observations.

    more » « less