The Space Density of Ultra-luminous QSOs at the End of Reionization Epoch by the QUBRICS Survey and the AGN Contribution to the Hydrogen Ionizing Background

Motivated by evidences favoring a rapid and late hydrogen reionization process completing at z ∼ 5.2–5.5 and mainly driven by rare and luminous sources, we have reassessed the estimate of the space density of ultra-luminous QSOs at z ∼ 5 in the framework of the QUBRICS survey. A ∼ 90% complete sample of 14 spectroscopically confirmed QSOs at M 1450 ≤ −28.3 and 4.5 ≤ z ≤ 5.0 has been derived in an area of 12,400 deg2, thanks to multiwavelength selection and Gaia astrometry. The space density of z ∼ 5 QSOs within −29.3 ≤ M 1450 ≤ −28.3 is three times higher than previous determinations. Our results suggest a steep bright-end slope for the QSO luminosity function at z ∼ 5 and a mild redshift evolution of the space density of ultrabright QSOs (M 1450 ∼ −28.5) at 3 < z < 5.5, in agreement with the redshift evolution of the much fainter active galactic nucleus (AGN) population at M 1450 ∼ −23. These findings are consistent with a pure density evolution for the AGN population at z > 3. Adopting our z ∼ 4 QSO luminosity function and applying a mild density evolution in redshift, a photoionization rate of Γ HI = 0.46 − 0.09 + 0.17 × 10 − 12 s − 1 has been obtained at z = 4.75, assuming an escape fraction of ∼70% and a steep faint-end slope of the AGN luminosity function. The derived photoionization rate is ∼50–100% of the ionizing background measured at the end of the reionization epoch, suggesting that AGNs could play an important role in the cosmological reionization process.

Constraints on the temperature-density relation of the intergalactic medium with non-negligible absorber spatial structure

We investigate evolution of physical parameters of the intergalactic medium using an analysis of Lya forest lines detected towards distant quasars. We used the enlarged sample of 98 quasars obtained with Keck/HIRES and VLT/UVES. We show that taking into account a finite spatial size of absorbers, regulated by pressure smoothing, significantly affects the inferred thermal parameters of the intergalactic gas, such as the hydrogen photoionization rate and parameters of the temperature-density relation. Using Bayesian framework we constrained for the first time the scale parameter between the Jeans length and characteristic size of the absorbers. We also discuss limitations of the method based on the analysis of the minimal broadending of Lya lines, which stem from the patchy nature of He II reionization.

Exploring reionization and high-z galaxy observables with recent multiredshift MWA upper limits on the 21-cm signal

ABSTRACT We use the latest multiredshift (z = 6.5−8.7) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. We explore these limits in the context of reionization astrophysics by using 21CMMC to connect the disfavoured regions of parameter space to existing observational constraints on reionization such as high-z galaxy ultraviolet (UV) luminosity functions, the background UV photoionization rate, the intergalactic medium (IGM) neutral fraction, the electron scattering optical depth and the soft-band X-ray emissivity. We find the vast majority of disfavoured models to already be inconsistent with existing observational constraints. These can be broadly classified into two types of models: (i) ‘cold’ reionization and (ii) pure matter density fluctuations in a cold, neutral IGM (i.e. no reionization). Interestingly, a small subsample of models inconsistent with the MWA is consistent with the aforementioned constraints (excluding the X-ray emissivity). This implies that the current MWA limits are already providing unique information to disfavour models of reionization, albeit extremely weakly. We also provide the first limits on the soft-band X-ray emissivity from galaxies at high redshifts, finding 1σ lower limits of ϵX, 0.5−2 keV ≳ 1034.5 erg s−1 Mpc−3. Finally, we recover 95 per cent disfavoured limits on the IGM spin temperature of $\bar{T}_{\rm S}\lesssim$ 1.3, 1.4, 1.5, 1.8, 2.1, and 2.4 K at z = 6.5, 6.8, 7.1, 7.8, 8.2, and 8.7. With this, we infer the IGM must have undergone, at the very least, a small amount of X-ray heating. Note, the limits on ϵX, 0.5–2 keV and $\bar{T}_{\rm S}$ are conditional on the IGM neutral fraction.

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

Measuring the temperature and profiles of Ly α absorbers

ABSTRACT The distribution of the absorption line broadening observed in the Ly α forest carries information about the temperature, T, and widths, λF, of the filaments in the intergalactic medium (IGM), and the background hydrogen photoionization rate, $\Gamma _{\rm H\, \small {I}}$. In this work, we present and test a new method for inferring T and λF and $\Gamma _{\rm H\, \small {I}}$ from combining the distribution of the absorption line broadening and the median flux. The method accounts for any underlying degeneracies. We apply our method to mock spectra from the reference model of the EAGLE cosmological simulation, and we demonstrate that we are able to reconstruct the IGM properties.

Radio observations of evaporating objects in the Cygnus OB2 region

We present observations of the Cygnus OB2 region obtained with the Giant Metrewave Radio Telescope (GMRT) at frequencies of 325 and 610 MHz. In this contribution we focus on the study of proplyd-like objects (also known as free-floating evaporating gas globules or frEGGs) that typically show an extended cometary morphology. We identify eight objects previously studied at other wavelengths and derive their physical properties by obtaining their optical depth at radio-wavelengths. Using their geometry and the photoionization rate needed to produce their radio-continuum emission, we find that these sources are possibly ionized by a contribution of the stars Cyg OB2 #9 and Cyg OB2 #22. Spectral index maps of the eight frEGGs were constructed, showing a flat spectrum in radio frequencies in general. We interpret these as produced by optically thin ionized gas, although it is possible that a combination of thermal emission, not necessarily optically thin, produced by a diffuse gas component and the instrument response (which detects more diffuse emission at low frequencies) can artificially generate negative spectral indices. In particular, for the case of the Tadpole we suggest that the observed emission is not of non-thermal origin despite the presence of regions with negative spectral indices in our maps.

Evolution of the AGN UV luminosity function from redshift 7.5

ABSTRACTDeterminations of the ultraviolet (UV) luminosity function of active galactic nuclei (AGN) at high redshifts are important for constraining the AGN contribution to reionization and understanding the growth of supermassive black holes. Recent inferences of the luminosity function suffer from inconsistencies arising from inhomogeneous selection and analysis of data. We address this problem by constructing a sample of more than 80 000 colour-selected AGN from redshift $z$ = 0 to 7.5 using multiple data sets homogenized to identical cosmologies, intrinsic AGN spectra, and magnitude systems. Using this sample, we derive the AGN UV luminosity function from redshift $z$ = 0 to 7.5. The luminosity function has a double power-law form at all redshifts. The break magnitude M* shows a steep brightening from M* ∼ −24 at $z$ = 0.7 to M* ∼ −29 at $z$ = 6. The faint-end slope β significantly steepens from −1.9 at $z$ < 2.2 to −2.4 at $z$ ≃ 6. In spite of this steepening, the contribution of AGN to the hydrogen photoionization rate at $z$ ∼ 6 is subdominant (<3 per cent), although it can be non-negligible (∼10 per cent) if these luminosity functions hold down to M1450 = −18. Under reasonable assumptions, AGN can reionize He ii by redshift $z$ = 2.9. At low redshifts ($z$ < 0.5), AGN can produce about half of the hydrogen photoionization rate inferred from the statistics of H i absorption lines in the intergalactic medium. Our analysis also reveals important systematic errors in the data, which need to be addressed and incorporated in the AGN selection function in future in order to improve our results. We make various fitting functions, codes, and data publicly available.