skip to main content

Title: The bolometric quasar luminosity function at z = 0–7
ABSTRACT In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from z = 0 to z = 7. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft, and hard X-ray in past decades. Our method follows Hopkins et al. with an updated quasar SED model and bolometric and extinction corrections. The new best-fitting bolometric quasar luminosity function behaves qualitatively different from the old Hopkins model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at z ≳ 2. Due to the paucity of measurements at the faint end, the faint end slope at z ≳ 5 is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at z ≳ 5. Further multiband observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this more » model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density, and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at z ≲ 3. The predicted CXB spectrum, cosmic SMBH mass density, and SMBH mass function are generally consistent with existing observations. « less
; ; ; ; ; ;
Award ID(s):
1652522 1715847
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
3252 to 3275
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT X-ray observations provide a unique probe of the accretion disc corona of supermassive black holes (SMBHs). In this paper, we present a uniform Chandra X-ray data analysis of a sample of 152 z ≥ 4.5 quasars. We firmly detect 46 quasars of this sample in 0.5–2 keV above 3σ and calculate the upper limits of the X-ray flux of the remaining. We also estimate the power-law photon index of the X-ray spectrum of 31 quasars. 24 of our sample quasars are detected in the FIRST or NVSS radio surveys; all of them are radio-loud. We statistically compare the X-ray propertiesmore »of our z ≥ 4.5 quasars to other X-ray samples of active galactic nuclei (AGNs) at different redshifts. The relation between the rest-frame X-ray luminosity and other quasar parameters, such as the bolometric luminosity, UV luminosity, or SMBH mass, shows large scatters. These large scatters can be attributed to the narrow luminosity range at the highest redshift, the large measurement error based on relatively poor X-ray data, and the inclusion of radio-loud quasars in the sample. The LX–LUV relationship is significantly sublinear. We do not find a significant redshift evolution of the LX–LUV relation, expressed either in the slope of this relation, or the departure of individual AGNs from the best-fitting αOX–LUV relation (ΔαOX). The median value of the X-ray photon index is Γ ≈ 1.79, which does not show redshift evolution from z = 0 to z ∼ 7. The X-ray and UV properties of the most distant quasars could potentially be used as a standard candle to constrain cosmological models. The large scatter of our sample on the Hubble diagram highlights the importance of future large unbiased deep X-ray and radio surveys in using quasars in cosmological studies.« less
  2. Abstract The luminosity function of active galactic nuclei (AGN) probes the history of supermassive black hole assembly and growth across cosmic time. To mitigate selection biases, we present a consistent analysis of the AGN luminosity functions (LFs) derived for both X-ray and mid-infrared (MIR) selected AGN in the XMM-Large Scale Structure field. There are 4268 AGN used to construct the MIR luminosity function (IRLF) and 3427 AGN used to construct the X-ray luminosity function (XLF), providing the largest census of the AGN population out to z = 4 in both bands with significant reduction in uncertainties. We are able formore »the first time to see the knee of the IRLF at z > 2 and observe a flattening of the faint-end slope as redshift increases. The bolometric luminosity density, a proxy for the cosmic black hole accretion history, computed from our LFs, shows a peak at z ≈ 2.25, consistent with recent estimates of the peak in the star formation rate density (SFRD). However, at earlier epochs, the AGN luminosity density is flatter than the SFRD. If confirmed, this result suggests that the build up of black hole mass outpaces the growth of stellar mass in high-mass systems at z ≳ 2.5. This is consistent with observations of redshift z ∼ 6 quasars that lie above the local M − σ relationship. The luminosity density derived from the IRLF is higher than that from the XLF at all redshifts. This is consistent with the dominant role of obscured AGN activity in the cosmic growth of supermassive black holes.« 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 The James Webb Space Telescope will have the power to characterize high-redshift quasars at z ≥ 6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e. with bolometric luminosity $L_{\rm bol}\geqslant 10^{46}\, \rm erg/s$) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with $L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$ opens a path to understand the build-up of more normal BHs at z ≥ 6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE,more »and SIMBA large-scale cosmological simulations do not agree on whether BHs at z ≥ 4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the MBH − M⋆ scaling relation at z = 0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z ≥ 4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with $L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$ observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z ∼ 6, with the same diversity of subgrid physics, in order to gain statistics on the most extreme objects at high redshift.« less
  5. ABSTRACT Observations of massive galaxies at low redshift have revealed approximately linear scaling relations between the mass of a supermassive black hole (SMBH) and properties of its host galaxy. How these scaling relations evolve with redshift and whether they extend to lower-mass galaxies, however, remain open questions. Recent galaxy formation simulations predict a delayed, or ‘two-phase,’ growth of SMBHs: slow, highly intermittent BH growth due to repeated gas ejection by stellar feedback in low-mass galaxies, followed by more sustained gas accretion that eventually brings BHs on to the local scaling relations. The predicted two-phase growth implies a steep increase, ormore »‘kink,’ in BH-galaxy scaling relations at a stellar mass $\rm {M}_{*}\sim 5\times 10^{10}$ M⊙. We develop a parametric, semi-analytic model to compare different SMBH growth models against observations of the quasar luminosity function (QLF) at z ∼ 0.5−4. We compare models in which the relation between SMBH mass and galaxy mass is purely linear versus two-phase models. The models are anchored to the observed galaxy stellar mass function, and the BH mass functions at different redshifts are consistently connected by the accretion rates contributing to the QLF. The best fits suggest that two-phase evolution is significantly preferred by the QLF data over a purely linear scaling relation. Moreover, when the model parameters are left free, the two-phase model fits imply a transition mass consistent with that predicted by simulations. Our analysis motivates further observational tests, including measurements of BH masses and active galactic nuclei activity at the low-mass end, which could more directly test two-phase SMBH growth.« less