We combine high resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large scale distribution and redshift evolution of the visibi
lity of Lyman-alpha emission in 6<=z<=8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyman-alpha transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyman-alpha emission line due to resonant scattering. We reproduce a decline in the space density of Lyman-alpha emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z=6-8, and a typical Lyman-alpha line velocity offset of 100 km/s redward of systemic at z=6 which decreases toward higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.
The quasar ULAS J1120+0641 at redshift z=7.085 has a highly ionised near zone which is smaller than those around quasars of similar luminosity at z~6. The spectrum also exhibits evidence for a damping wing extending redward of the systemic Lya redshi
ft. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionisation state of the inhomogeneous IGM surrounding this quasar. Our simulations show that the transmission profile is consistent with an IGM in the vicinity of the quasar with a volume averaged HI fraction of f_HI>0.1 and that ULAS J1120+0641 has been bright for 10^6--10^7 yr. The observed spectrum is also consistent with smaller IGM neutral fractions, f_HI ~ 10^-3--10-4, if a damped Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ~5 per cent of our simulated sight-lines for a bright phase of 10^6--10^7 yr. Unless ULAS J1120+0641 grows during a previous optically obscured phase, the low age inferred for the quasar adds to the theoretical challenge of forming a 2x10^9 M_sol black hole at this high redshift.
We present a suite of full hydrodynamical cosmological simulations that quantitatively address the impact of neutrinos on the (mildly non-linear) spatial distribution of matter and in particular on the neutral hydrogen distribution in the Intergalact
ic Medium (IGM), which is responsible for the intervening Lyman-alpha absorption in quasar spectra. The free-streaming of neutrinos results in a (non-linear) scale-dependent suppression of power spectrum of the total matter distribution at scales probed by Lyman-alpha forest data which is larger than the linear theory prediction by about 25% and strongly redshift dependent. By extracting a set of realistic mock quasar spectra, we quantify the effect of neutrinos on the flux probability distribution function and flux power spectrum. The differences in the matter power spectra translate into a ~2.5% (5%) difference in the flux power spectrum for neutrino masses with Sigma m_{ u} = 0.3 eV (0.6 eV). This rather small effect is difficult to detect from present Lyman-alpha forest data and nearly perfectly degenerate with the overall amplitude of the matter power spectrum as characterised by sigma_8. If the results of the numerical simulations are normalized to have the same sigma_8 in the initial conditions, then neutrinos produce a smaller suppression in the flux power of about 3% (5%) for Sigma m_{ u} = 0.6$ eV (1.2 eV) when compared to a simulation without neutrinos. We present constraints on neutrino masses using the Sloan Digital Sky Survey flux power spectrum alone and find an upper limit of Sigma m_{ u} < 0.9$ eV (2 sigma C.L.), comparable to constraints obtained from the cosmic microwave background data or other large scale structure probes.
The next generation of X-ray telescopes have the potential to detect faint quasars at very high redshift and probe the early growth of massive black holes (BHs). We present modelling of the evolution of the optical and X-ray AGN luminosity function a
t 2 < z < 6 based on a CDM merger-driven model for the triggering of nuclear activity combined with a variety of fading laws. We extrapolate the merger-driven models to z > 6 for a range of BH growth scenarios. We predict significant numbers of sources at z ~ 6 with fluxes just an order of magnitude below the current detection limits and thus detectable with XEUS and Constellation-X, relatively independently of the fading law chosen. The predicted number of sources at even higher redshift depends sensitively on the early growth history of BHs. For passive evolution models in which BHs grow constantly at their Eddington limit, detectable BHs may be rare beyond z ~ 10 even with Generation-X. However, in the more probable scenario that BH growth at z > 6 can be described by passive evolution with a small duty cycle, or by our merger driven accretion model, then we predict that XEUS and Generation-X will detect significant numbers of black holes out to z ~ 10 and perhaps beyond.
