No Arabic abstract
We investigate the effect of helium on hydrogen reionisation using a hydrodynamical simulation combined with the cosmological radiative transfer code CRASH. The simulations are run in a 35.12/h comoving Mpc box using a variety of assumptions for the amplitude and power-law extreme-UV (EUV) spectral index, alpha, of the ionising emissivity. We use an empirically motivated prescription for ionising sources which ensures all of the models are consistent with constraints on the Thomson scattering optical depth and the hydrogen photo-ionisation rate at z=6. The inclusion of helium slightly delays reionisation due to the small number of ionising photons which reionise neutral helium instead of hydrogen. However, helium has a significant impact on the thermal state of the IGM. Models with alpha=3 produce IGM temperatures at the mean density at z=6 which are about 20 % higher compared to models without helium photo-heating. Harder EUV indices produce even larger IGM temperature boosts. A comparison to recent observational estimates of the IGM temperature at z=5 - 6 suggests that hydrogen reionisation was primarily driven by pop-II stellar sources with a soft EUV index, alpha<3. We also find that faint, as yet undetected galaxies, characterised by a luminosity function with a steepening faint-end slope and an increasing Lyman continuum escape fraction (fesc=0.5), are required to reproduce the ionising emissivity used in our simulations at z>6. Finally, we note there is some tension between recent observational constraints which indicate the IGM is > 10% neutral by volume z=7, and estimates of the ionising emissivity at z=6 which indicate only between 1 and 3 ionising photons are emitted per hydrogen atom over a Hubble time. This tension may be alleviated by either a lower neutral fraction at z=7 or an IGM which still remains a few % neutral by volume at z=6.
We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the Lyman-$alpha$ forest transmitted flux at redshifts $4.2leq z leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k sim 0.1rm,km^{-1},s$, by $sim 10$ per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, $k<0.03rm,km^{-1},s$, the power spectrum is instead enhanced by $10$-$50$ per cent by large scale spatial variations in the neutral hydrogen fraction. The effect of inhomogeneous reionisation must therefore be accounted for in analyses of forthcoming high precision measurements. We provide a correction for the Lyman-$alpha$ forest power spectrum at $4.1leq z leq 5.4$ in a form that can be easily applied within other parameter inference frameworks. We perform a Bayesian analysis of mock data to assess the extent of systematic biases that may arise in measurements of the intergalactic medium if ignoring this correction. At the scales probed by current high resolution Lyman-$alpha$ forest data at $z>4$, $0.006 rm ,km^{-1},sleq k leq 0.2 rm, km^{-1},s$, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of $sim 10$ per cent. However, for $5$ per cent uncertainties, $sim 1sigma$ shifts between the estimated and true parameters occur.
Modelling reionization often requires significant assumptions about the properties of ionizing sources. Here, we infer the total output of hydrogen-ionizing photons (the ionizing emissivity, $dot{N}_textrm{ion}$) at $z=4-14$ from current reionization constraints, being maximally agnostic to the properties of ionizing sources. We use a Bayesian analysis to fit for a non-parametric form of $dot{N}_textrm{ion}$, allowing us to flexibly explore the entire prior volume. We infer a declining $dot{N}_textrm{ion}$ with redshift at $z>6$, which can be used as a benchmark for reionization models. Model-independent reionization constraints from the CMB optical depth and Ly$alpha$ and Ly$beta$ forest dark pixel fraction produce $dot{N}_textrm{ion}$ evolution ($dlog_{10}dot{N}_textrm{ion}/dz|_{z=6rightarrow8} = -0.31pm0.35$ dex) consistent with the declining UV luminosity density of galaxies, assuming constant ionizing photon escape fraction and efficiency. Including measurements from Ly$alpha$ damping of galaxies and quasars produces a more rapid decline: $dlog_{10}dot{N}_textrm{ion}/dz|_{z=6rightarrow8} =-0.44pm0.22$ dex, steeper than the declining galaxy luminosity density (if extrapolated beyond $M_mathrm{UV} lesssim -13$), and constrains the mid-point of reionization to $z = 6.93pm0.14$.
The effect of the primeval sources of radiation and particles on the thermodynamical state of the intergalactic medium during the Epoch of Reionisation is still unclear. In this work, we explore the contribution of electrons accelerated in the jets of high-redshift microquasars to heating and ionising the intergalactic medium. We develop Monte Carlo simulations of the propagation and energy deposition of these electrons as they travel away from their sources. We find that microquasars contribute significantly to heating the intergalactic medium and are effective ionisers only near the galaxies. Their effect on heating is of the same order of magnitude than that of CRs from SNe.
The recent measurement of an ionising mean free path $lambda_{text{mfp}}<1$ pMpc at $z=6$ challenges our understanding of the small-scale structure of the intergalactic medium (IGM) at the end of reionisation. We introduce a new method to constrain mfp at $z=6$ by using lower limits on the individual free paths of ionisation around quasars. Lyman-limit absorbers with a density sufficient to halt ionising photons produce strong absorption in the 6 lowest-energy Lyman transitions, in the absence of which a robust lower limit can be placed on the individual free path. Applying this method to a set of $26$ quasars at $5.5<z<6.5$, we find that $80%$ of bright quasars ($M_{1450}<-26.5$) require individual free paths larger than $2$ pMpc. We model the relation between opacity $kappa$ and photo-ionisation rate $Gamma$ via the parameter $xi$ such that $kappaproptoGamma^{-xi}$, and pose joint limits on mfp and $xi$. For the nominal value of $xi=2/3$, we constrain $lambda_{text{mfp}}>0.31 (0.18)$ pMpc at $2sigma (3sigma)$: a much tighter lower bound than obtained through traditional stacking methods. Our constraints get significantly stronger for lower values of $xi$. New constraints on mfp and $xi$ are crucial to our understanding of the reionisation-era IGM.
We present an updated determination of the z ~ 4 QSO luminosity function (QLF), improving the quality of the determination of the faint end of the QLF presented in Glikman et al. (2010). We have observed an additional 43 candidates from our survey sample, yielding one additional QSO at z = 4.23 and increasing the completeness of our spectroscopic follow-up to 48% for candidates brighter than R = 24 over our survey area of 3.76 deg2. We study the effect of using K-corrections to compute the rest-frame absolute magnitude at 1450A compared with measuring M1450 directly from the object spectra. We find a luminosity-dependent bias: template-based K-corrections overestimate the luminosity of low-luminosity QSOs, likely due to their reliance on templates derived from higher luminosity QSOs. Combining our sample with bright quasars from the Sloan Digital Sky Survey and using spectrum-based M1450 for all the quasars, we fit a double-power-law to the binned QLF. Our best fit has a bright-end slope, {alpha} = 3.3pm0.2, and faint-end slope, {beta} = 1.6(+0.8/-0.6). Our new data revise the faint-end slope of the QLF down to flatter values similar to those measured at z ~ 3. The break luminosity, though poorly constrained, is at M* = -24.1(+0.7/-1.9), approximately 1 - 1.5 mag fainter than at z ~ 3. This QLF implies that QSOs account for about half the radiation needed to ionize the IGM at these redshifts.