Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmological simulations. Using fast and versatile semi-numeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the HII morphology of reionization can drive the large-scale distribution of absorption systems. Specifically, we find that absorbers are highly correlated with the density field on small scales, and then become anti-correlated with the UVB on large scales. After reionization, the large-scale power spectrum of the absorbers traces the UVB power spectrum, which can be predicted with a simple analytic extension of the halo-model. During reionization, absorbers tend to preferentially lie inside overdensities (i.e. filaments) of the recently-ionized intergalactic medium (IGM). Absorbers may also dominate the small-scale (k > 1/Mpc) 21-cm power during and after reionization. Conversely, they smooth the contrast on moderate scales. Once the HII regions grow to surpass the mfp, the absorbers add to the large-scale 21-cm power. Our results should prove useful in interpreting future observations of the reionization epoch.

We investigate the effects of non-Gaussianity in the primordial density field on the reionization history. We rely on a semi-analytic method to describe the processes acting on the intergalactic medium (IGM), relating the distribution of the ionizing sources to that of dark matter haloes. Extending previous work in the literature, we consider models in which the primordial non-Gaussianity is measured by the dimensionless non-linearity parameter f_NL, using the constraints recently obtained from cosmic microwave background data. We predict the ionized fraction and the optical depth at different cosmological epochs assuming two different kinds of non-Gaussianity, characterized by a scale-independent and a scale-dependent f_NL and comparing the results to those for the standard Gaussian scenario. We find that a positive f_NL enhances the formation of high-mass haloes at early epochs, when reionization begins, and, as a consequence, the IGM ionized fraction can grow by a factor up to 5 with respect to the corresponding Gaussian model. The increase of the filling factor has a small impact on the reionization optical depth and is of order ~ 10 per cent if a scale-dependent non-Gaussianity is assumed. Our predictions for non-Gaussian models are in agreement with the latest WMAP results within the error bars, but a higher precision is required to constrain the scale dependence of non-Gaussianity.

In this paper we investigate the effects that a dynamic dark energy component dominant in the universe at late epochs has on reionization. We follow the evolution of HII regions with the analytic approach of Furlanetto and Oh (2005) in two different universes for which we assume the Peebles and Ratra (2003) and Brax and Martin (2000) quintessence models and we compare our results to the LCDM scenario. We show that, for a fixed ionization efficiency, at the same cosmological epoch the topology of bubbles is dominated by high-mass objects and the characteristic size of the ionized regions is slightly smaller than in the LCDM model, especially at the latest stages of reionization, due to the higher recombination efficiency. As a consequence, the bubbles' `epoch of overlap' happens earlier than in LCDM. Finally, we show how the different evolution of the HII regions affects the transmission of the high-z QSO spectra, reducing the Lyman flux absorption at small optical depths.

We investigate the ability of hierarchical models of QSO formation and evolution to match the observed luminosity, number counts and spatial clustering of quasars at redshift z<2. These models assume that the QSO emission is triggered by galaxy mergers, that the mass of the central black hole correlates with halo properties and that quasars shine at their Eddington luminosity except, perhaps, during the very early stages of evolution. We find that models based on simple analytic approximations successfully reproduce the observed B-band QSO luminosity function at all redshifts, provided that some mechanisms is advocated to quench mass accretion within haloes larger than about 1e13 Msun that host bright quasars. These models also match the observed strength of QSO clustering at z~0.8. At larger redshifts, however, they underpredict the QSO biasing which, instead, is correctly reproduced by semi-analytic models in which the halo merger history and associated BHs are followed by Monte Carlo realizations of the merger hierarchy. We show that the disagreement between the luminosity function predicted by semi-analytic models and observations can be ascribed to the use of B-band data, which are a biased tracer of the quasar population, due to obscuration.