We combine fast radiative transfer calculations with high resolution hydrodynamical simulations to study an epoch of early hydrogen reionisation by primordial stellar sources at redshifts 15<z<30. With relatively conservative assumptions, population III star formation proceeds in a self-regulated manner both locally and globally and, for a conventional LCDM cosmology, can significantly reionise the intergalactic medium between 15<z<20 as long as a large fraction of ionising photons can escape from these earliest galaxies. We then combine these results with our earlier work focusing on the role of population II stars in galaxies with virial temperatures >10^4K at redshifts 5<z< 20. Hence, we construct a complete reionisation history of the Universe which matches the Thomson optical depths as measured by the WMAP satellite as well as the evolution of the Gunn Peterson optical depth as seen in the asborption spectra of the higest redshift quasars. We find that even with conservative estimates for the impact of negative feedback mechanisms, primordial stellar sources contribute significantly to early reionisation. Future observations of a Thomson optical depth of tau_e>~0.13 would bolster the claim for the existence of population III stars similar to the ones studied here.

(abridged) We study high redshift structure formation and reionization in a LCDM universe under the assumption that the spectral power index of primordial density fluctuations is a function of length scale. We adopt a particular formulation of the running spectral index (RSI) model as suggested by the recent WMAP data. While early structure forms hierarchically in the RSI model, the reduced power on small scales causes a considerable delay in the formation epoch of low mass (~ 10^6 Msun) ``mini-halos'' compared to the LCDM model. The extremely small number of gas clouds in the RSI model indicates that reionization is initiated later than z<15, generally resulting in a smaller total Thomson optical depth than in the LCDM model. By carrying out radiative transfer calculations, we also study reionization by stellar populations formed in galaxies. Even with a top-heavy intial mass function representing an early population of massive stars and/or an extraordinarily high photon emission rate from galaxies, the total optical depth can only be as large as tau ~ 0.1 for reasonable models of early star-formation. The RSI model is thus in conflict with the large Thomson optical depth inferred by the WMAP satellite.

We study first structure formation in Lambda-dominated universes using large cosmological N-body/SPH simulations. We consider a standard LCDM model and a LWDM model in which the mass of the dark matter particles is taken to be m_X=10 keV. The linear power spectrum for the LWDM model has a characteristic cut-off at a wavenumber k=200 /Mpc, suppressing the formation of low mass (< 10^6 Msun) nonlinear objects early on. The absence of low mass halos in the WDM model makes the formation of primordial gas clouds with molecular hydrogen very inefficient at high redshifts. The first star-forming gas clouds form at z~21 in the WDM model, considerably later than in the CDM counterpart, and the abundance of these gas clouds differs by an order of magnitude between the two models. We carry out radiative transfer calculations by embedding massive Population III stars in the gas clouds. We show that the volume fraction of ionized gas rises up close to 100% by z=18 in the CDM case, whereas that of the WDM model remains extremely small at a level of a few percent. Thus the WDM model with m_X=10 keV is strongly inconsistent with the observed high optical depth by the WMAP satellite.

We study the observability of the reionization epoch through the 21 cm hyperfine transition of neutral hydrogen. We use a high-resolution cosmological simulation (including hydrodynamics) together with a fast radiative transfer algorithm to compute the evolution of 21 cm emission from the intergalactic medium (IGM) in several different models of reionization. We show that the mean brightness temperature of the IGM drops from dT_b~25 mK to dT_b~0.01 mK during overlap (over a frequency interval of ~25 MHz), while the root mean square brightness temperature fluctuations on small scales drop abruptly from ~10 mK before overlap to ~0.1 mK at the end of overlap. We show that 21 cm observations can efficiently discriminate models with a single early reionization epoch from models with two distinct reionization episodes.

We study the reionisation of the Universe by stellar sources using a numerical approach that combines fast 3D radiative transfer calculations with high resolution hydrodynamical simulations. Ionising fluxes for the sources are derived from intrinsic star formation rates computed in the underlying hydrodynamical simulations. Our mass resolution limit for sources is M~ 4.0 x 10^7 h^-1 M_sol, which is roughly an order of magnitude smaller than in previous studies of this kind. Our calculations reveal that the reionisation process is sensitive to the inclusion of dim sources with masses below ~10^9 h^-1 M_sol. We present the results of our reionisation simulation assuming a range of escape fractions for ionising photons and make statistical comparisons with observational constraints on the neutral fraction of hydrogen at z~6 derived from the z=6.28 SDSS quasar of Becker and coworkers. Our best fitting model has an escape fraction of ~20% and causes reionisation to occur by z~8, although the IGM remains fairly opaque until z~6. In order to simultaneously match the observations from the z=6.28 SDSS quasar and the optical depth measurement from WMAP with the sources modeled here, we require an evolving escape fraction that rises from f_esc=0.20 near z~6 to f_esc>~10 at z~18.

Using radiative transfer calculations and cosmological simulations of structure formation, we study constraints that can be placed on the nature of the cosmic ultraviolet (UV) background in the redshift interval 2.5< z <5. Our approach makes use of observational estimates of the opacities of hydrogen and singly ionised helium in the intergalactic medium during this epoch. In particular, we model the reionisation of He II by sources of hard ultraviolet radiation, i.e. quasars, and infer values for our parameterisation of this population from observational estimates of the opacity of the He II Lyman-alpha forest. Next, we estimate the photoionisation rate of H I from these sources and find that their contribution to the ionising background is insufficient to account for the measured opacity of the H I Lyman-alpha forest at a redshift z 3. This motivates us to include a soft, stellar component to the ionising background to boost the hydrogen photoionisation rate, but which has a negligible impact on the He II opacity. In order to simultaneously match observational estimates of the H I and He II opacities, we find that galaxies and quasars must contribute about equally to the ionising background in H I at z~ 3.

We study the reionization of Helium II by quasars using a numerical approach that combines 3D radiative transfer calculations with cosmological hydrodynamical simulations. Sources producing the ionizing radiation are selected according to an empirical quasar luminosity function and are assigned luminosities according to their intrinsic masses. We present models in which these parameters are varied and examine characteristics of the resultant reionization process that distinguish the various cases. In addition, we extract artificial spectra from the simulations and quantify statistical properties of the spectral features in each model. We find that the most important factor affecting the evolution of He II reionization is the cumulative number of ionizing photons that are produced by the sources. Comparisons between He II opacities measured observationally and those obtained by our analysis reveal that the available ranges in plausible values for the parameters provide enough leeway to provide a satisfactory match. However, one property common to all our calculations is that the epoch of Helium II reionization must have occurred at a redshift between 3 < z < 4. If so, future observational programs will be able to directly trace the details of the ionization history of helium and probe the low density phase of the intergalactic medium during this phase of the evolution of the Universe.

The incorporation of radiative transfer effects into cosmological hydrodynamical simulations is essential for understanding how the intergalactic medium (IGM) makes the transition from a neutral medium to one that is almost fully ionized. Here, we present an approximate numerical method designed to study in a statistical sense how a cosmological density field is ionized by a set of discrete point sources. A diffuse background radiation field is also computed self-consistently in our procedure. The method requires relatively few time steps and can be employed with simulations having high resolution. We describe the details of the algorithm and provide a description of how the method can be applied to the output from a pre-existing cosmological simulation to study the systematic reionization of a particular ionic species. As a first application, we compute the reionization of He II by quasars in the redshift range 3 to 6.

We estimate the contribution of extragalactic radio sources to fluctuations in sky temperature over the range of frequencies (10-300 GHz) used for Cosmic Microwave Background (CMB) anisotropy measurements. CMB anisotropy observations at high resolution and low frequencies are especially sensitive to this foreground. We have compiled a catalog of 2207 bright radio sources, including 758 sources with flux measurements at 90 GHz. We develop a method to extrapolate the source spectra and predict skymaps of extragalactic radio sources at instrument resolutions of 10 arcmin to 10 degrees FWHM. Our results indicate that the brightest sources will dominate microwave anisotropy for a wide range of resolutions and frequencies. Our skymaps predict the location and flux of the brightest radio sources at each frequency, making it straightforward to develop a template for masking the pixels containing them. This masking should be sufficient to protect high resolution CMB anisotropy observations from unacceptable radio source confusion.