The detection of an unexpectedly large number of faint candidate quasars (QSO) at $z \gtrsim 4$ has motivated Madau and Haardt to investigate a reionization scenario in which active galactic nuclei dominate the photon budget at all times. Their analytical study reveals that this picture is in agreement with the evolution of the HI volume fraction in the intergalactic medium (IGM) and the optical depth to Thomson scattering of the cosmic microwave background. We employ a suite of hydrodynamical simulations post-processed with a radiative transfer code to further investigate the properties of such a reionization history. In particular, we generate synthetic absorption-line spectra for HI and HeII that we compare with observational data. Although we confirm the analytical results mentioned above, we also find that the evolution of the IGM temperature and of the HeII optical depth are at odds with current observational constraints. Nevertheless, the QSO-dominated scenario presents the attractive feature that it naturally generates an inhomogeneous IGM and thus produces an extended tail at high values in the distribution function of the HI effective optical depth, in agreement with recent observations. We find tentative evidence that considering some QSO emission at redshift $z \gtrsim 5$ could reconcile such observations with numerical simulations that have so far failed to reproduce this feature. Finally, we show that collecting more HeII-absorption spectra at $z \gtrsim 3$ and studying their distinctive characteristics will be key to precisely constraining the QSO contribution to reionization.

We present a detailed analysis of the ionization and thermal structure of the intergalactic medium (IGM) around a high-redshift QSO using a large suite of cosmological, multi-frequency radiative transfer (RT) simulations, exploring the contribution from galaxies as well as the QSO, and the effect of X-rays and secondary ionization. We show that in high-z QSO environments both the central QSO and the surrounding galaxies concertedly control the reionization morphology of hydrogen and helium and have a non-linear impact on the thermal structure of the IGM. A QSO imprints a distinctive morphology on H II regions if its total ionizing photon budget exceeds that of the surrounding galaxies since the onset of hydrogen reionization; otherwise, the morphology shows little difference from that of H II regions produced only by galaxies. In addition, the spectral shape of the collective radiation field from galaxies and QSOs controls the thickness of the I-fronts. While a UV-obscured QSO can broaden the I-front, the contribution from other UV sources, either galaxies or unobscured QSO, is sufficient to maintain a sharp I-front. X-rays photons from the QSO are responsible for a prominent extended tail of partial ionization ahead of the I-front. QSOs leave a unique imprint on the morphology of He II / He III regions. We suggest that, while the physical state of the IGM is modified by QSOs, the most direct test to understand the role of galaxies and QSOs during reionization is to perform galaxy surveys in a region of sky imaged by 21 cm tomography.

We use hydrodynamic simulations post-processed with the radiative-transfer code RADAMESH to assess recent claims that the low HeII opacity observed in z>3 quasar spectra may be incompatible with models of HeII reionization driven by the observed population of active galactic nuclei (AGNs). In particular, building upon our previous work, we consider an early population of sources and start the radiative-transfer calculation at redshifts z>=5. Our model faithfully reproduces the emissivity of optically selected AGNs as inferred from measurements of their luminosity function. We find that HeII reionization is very extended in redshift ({\Delta} z>=2) and highly spatially inhomogeneous. In fact, mock spectra extracted from the simulations show a large variability in the evolution of the HeII effective optical depth within chunks of size {\Delta} z=0.04. Regions with low opacity ({\tau}^{eff}_{HeII}<3) can be found at high redshift, in agreement with the most recent observations of UV-transmitting quasars. At the highest redshift currently probed by observations (z~3.4), our updated model predicts a much lower HeII effective optical depth than previous simulations in the literature relieving most of the tension with the current data, that, however, still persists at about the (Gaussian) 1{\sigma} to 2{\sigma} level. Given the very small number of observed lines of sight, our analysis indicates that current data cannot rule out a purely AGN-driven scenario with high statistical significance.

We use a set of AMR hydrodynamic simulations post-processed with the radiative-transfer code RADAMESH to study how inhomogeneous HeII reionization affects the intergalactic medium (IGM). We propagate radiation from active galactic nuclei (AGNs) considering two scenarios for the time evolution of the ionizing sources. We find that HeII reionization takes place in a very inhomogeneous fashion, through the production of well separated bubbles of the ionized phase that subsequently percolate. Overall, the reionization process is extended in time and lasts for a redshift interval Delta z>1. At fixed gas density, the temperature distribution is bimodal during the early phases of HeII reionization and cannot be described by a simple effective equation of state. When HeII reionization is complete, the IGM is characterized by a polytropic equation of state with index gamma~1.20. This relation is appreciably flatter than at the onset of the reionization process (gamma=1.56) and also presents a much wider dispersion around the mean. We extract HI and HeII Ly-alpha absorption spectra from the simulations and fit Voigt profiles to them. We find that the regions where helium is doubly ionized are characterized by different probability density functions of the curvature and of the Doppler b parameters of the HI Ly-alpha forest as a consequence of the bimodal temperature distribution during the early phases of HeII reionization. The column-density ratio in HeII and HI shows a strong spatial variability. Its probability density function rapidly evolves with time reflecting the increasing volume fraction in which ionizing radiation is harder due to the AGN contribution. Finally we show that the number density of the flux-transmission windows per unit redshift and the mean size of the dark gaps in the HeII spectra have the potential to distinguish between different reionization scenarios. (abridged)