Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes.

We introduce gbpTrees: an algorithm for constructing merger trees from cosmological simulations, designed to identify and correct for pathological cases introduced by errors or ambiguities in the halo finding process. gbpTrees is built upon a halo matching method utilising pseudo-radial moments constructed from radially-sorted particle ID lists (no other information is required) and a scheme for classifying merger tree pathologies from networks of matches made to-and-from haloes across snapshots ranging forward-and-backward in time. Focusing on Subfind catalogs for this work, a sweep of parameters influencing our merger tree construction yields the optimal snapshot cadence and scanning range required for converged results. Pathologies proliferate when snapshots are spaced by $\lesssim{0.128}$ dynamical times; conveniently similar to that needed for convergence of semi-analytical modelling, as established by Benson etal Total merger counts are converged at the level of $\sim{5}$% for friends-of-friends (FoF) haloes of size $n_{\rm p}\gtrsim{75}$ across a factor of 512 in mass resolution, but substructure rates converge more slowly with mass resolution, reaching convergence of $\sim{10}$% for $n_{\rm p}\gtrsim{100}$ and particle mass $m_{\rm p}{\lesssim}10^{9}M_\odot$. We present analytic fits to FoF and substructure merger rates across nearly all observed galactic history ($z{\le}8.5$). While we find good agreement with the results presented by Fakhouri etal for FoF haloes, a slightly flatter dependance on merger ratio and increased major merger rates are found, reducing previously reported discrepancies with extended Press-Schechter estimates. When appropriately defined, substructure merger rates show a similar mass ratio dependance as FoF rates, but with stronger mass and redshift dependencies for their normalisation.

Annihilating dark matter (DM) models offer promising avenues for future DM detection, in particular via modification of astrophysical signals. However when modelling such potential signals at high redshift the emergence of both dark matter and baryonic structure, as well as the complexities of the energy transfer process, need to be taken into account. In the following paper we present a detailed energy deposition code and use this to examine the energy transfer efficiency of annihilating dark matter at high redshift, including the effects on baryonic structure. We employ the PYTHIA code to model neutralino-like DM candidates and their subsequent annihilation products for a range of masses and annihilation channels. We also compare different density profiles and mass-concentration relations for 10^5-10^7 M_sun haloes at redshifts 20 and 40. For these DM halo and particle models, we show radially dependent ionisation and heating curves and compare the deposited energy to the haloes' gravitational binding energy. We use the "filtered" annihilation spectra escaping the halo to calculate the heating of the circumgalactic medium and show that the mass of the minimal star forming object is increased by a factor of 2-3 at redshift 20 and 4-5 at redshift 40 for some DM models.

An interesting transient has been detected in one of our three Dark Energy Camera deep fields. Observations of these deep fields take advantage of the high red sensitivity of DECam on the Cerro Tololo Interamerican Observatory Blanco telescope. The survey includes the Y band with rest wavelength 1430{\AA} at z = 6. Survey fields (the Prime field 0555-6130, the 16hr field 1600-75 and the SUDSS New Southern Field) are deeper in Y than other infrared surveys. They are circumpolar, allowing all night to be used efficiently, exploiting the moon tolerance of 1 micron observations to minimize conflict with the Dark Energy Survey. As an i-band dropout (meaning that the flux decrement shortward of Lyman alpha is in the i bandpass), the transient we report here is a supernova candidate with z ~ 6, with a luminosity comparable to the brightest known current epoch superluminous supernova (i.e., ~ 2 x 10^11 solar luminosities).

We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis).

We present the Dark-ages Reionization and Galaxy-formation Observables from Numerical Simulations (DRAGONS) program and Tiamat, the collisionless N-body simulation program upon which DRAGONS is built. The primary trait distinguishing Tiamat from other large simulation programs is its density of outputs at high redshift (100 from z=35 to z=5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from "universal" behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20%) of galactic halos are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: i) halos need ~1.5 dynamical times following formation, and ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic halos and we suggest that care must be taken to ensure a physically meaningful match between the galaxy-formation physics of semi-analytic models and the halo finders supplying their input.

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $z\sim8$ UV luminosity function using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^\star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $\alpha = -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} \Psi^\star [\textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $z\sim8$, and cannot account for the apparent overdensity of bright galaxies compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $z\sim8$ source in our sample is in the range $\sim 3-15 \%$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $\mu>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.

An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-frequency signals from the tiles are amplified and band limited using analog signal conditioning units; sampled and channelized by digital-receivers. The signals from eight tiles are processed by a single digital-receiver, thus requiring 16 digital-receivers for the MWA. The main function of the digital-receivers is to digitize the broad-band signals from each tile, channelize them to form the sky-band, and transport it through optical fibers to a centrally located correlator for further processing. The digital-receiver firmware also implements functions to measure the signal power, perform power equalization across the band, detect interference-like events, and invoke diagnostic modes. The digital-receiver is controlled by high-level programs running on a single-board-computer. This paper presents the digital-receiver design, implementation, current status, and plans for future enhancements.

We predict Lyman-$\alpha$ (Ly$\alpha$) luminosity functions (LFs) of Ly$\alpha$-selected galaxies (Ly$\alpha$ emitters, or LAEs) at $z=3-6$ using the phenomenological model of Dijkstra & Wyithe (2012). This model combines observed UV-LFs of Lyman-break galaxies (LBGs, or drop out galaxies), with constraints on their distribution of Ly$\alpha$ line strengths as a function of UV-luminosity and redshift. Our analysis shows that while Ly$\alpha$ LFs of LAEs are generally not Schechter functions, these provide a good description over the luminosity range of $\log_{10}( L_{\alpha}/{\rm erg}\,{\rm s}^{-1})=41-44$. Motivated by this result, we predict Schechter function parameters at $z=3-6$. Our analysis further shows that (i) the faint end slope of the Ly$\alpha$ LF is steeper than that of the UV-LF of Lyman-break galaxies, (with a median $\alpha_{Ly\alpha} < -2.0$ at $z\gtrsim 4$), and (ii) a turn-over in the Ly$\alpha$ LF of LAEs at Ly$\alpha$ luminosities $10^{40}$ erg s$^{-1}<L_{\alpha}\lesssim 10^{41}$ erg s$^{-1}$ may signal a flattening of UV-LF of Lyman-break galaxies at $-12>M_{\rm UV}>-14$. We discuss the implications of these results - which can be tested directly with upcoming surveys - for the Epoch of Reionization.

The ionisation structure of the Intergalactic Medium (IGM) during reionisation is sensitive to the unknown galaxy formation physics that prevailed at that time. This structure introduces non-Gaussian statistics into the redshifted 21 cm fluctuation amplitudes that can only be studied through tomographic imaging, which will clearly discriminate between different galaxy formation scenarios. Imaging the ionisation structure and cosmological HII regions during reionisation is therefore a key goal for the SKA. For example, the SKA1-LOW baseline design with a 1 km diameter core will resolve HII regions expected from galaxy formation models which include strong feedback on low-mass galaxy formation. Imaging the smaller HII regions that result from galaxy formation in the absence of SNe feedback will also be possible for SKA1-LOW in the later stages of reionisation, but may require the greater sensitivity of SKA early in the reionisation era. In addition to having baselines long enough to resolve the HII regions, the field of view for SKA1-LOW reionisation experiments should be at least several degrees in order to image the largest HI structures towards the end of reionisation. The baseline design with 35 meter diameter stations has a field of view within a single primary pointing which is sufficient for this purpose.

We investigate the evolution of the star formation rate-stellar mass relation (SFR-M*) and Galaxy Stellar Mass Function (GSMF) of z ~ 4-7 galaxies, using cosmological simulations run with the smoothed particle hydrodynamics code P-GADGET3(XXL). We explore the effects of different feedback prescriptions (supernova driven galactic winds and AGN feedback), initial stellar mass functions and metal cooling. We show that our fiducial model, with strong energy-driven winds and early AGN feedback, is able to reproduce the observed stellar mass function obtained from Lyman-break selected samples of star forming galaxies at redshift 6 < z < 7. At z ~ 4, observed estimates of the GSMF vary according to how the sample was selected. Our simulations are more consistent with recent results from K-selected samples, which provide a better proxy of stellar masses and are more complete at the high mass end of the distribution. We find that in some cases simulated and observed SFR-M* relations are in tension, and this can lead to numerical predictions for the GSMF in excess of the GSMF observed. By combining the simulated SFR(M*) relationship with the observed star formation rate function at a given redshift, we argue that this disagreement may be the result of the uncertainty in the SFR-M* (Luv-M*) conversion. Our simulations predict a population of faint galaxies not seen by current observations.

The rapid decline in the number of strong Lyman Alpha (Lya) emitting galaxies at z > 6 provides evidence for neutral hydrogen in the IGM, but is difficult to explain with plausible models for reionization. We demonstrate that the observed reduction in Lya flux from galaxies at z > 6 can be explained by evolution in the escape fraction of ionizing photons, f_esc. We find that the median observed drop in the fraction of galaxies showing strong Lya emission, as well as the observed evolution of the Lya luminosity function both follow from a small increase in f_esc of Delta f_esc ~ 0.1 from f_esc ~ 0.6 at z ~ 6. This high escape fraction may be at odds with current constraints on the ionising photon escape fraction, which favor smaller values of f_esc < 20%. However, models that invoke a redshift evolution of f_ esc that is consistent with these constraints can suppress the z~7 Lya flux to the observed level, if they also include a small evolution in global neutral fraction of Delta x_HI ~ 0.2. Thus, an evolving escape fraction of ionising photons can be a plausible part of the explanation for evolution in the Lya emission of high redshift galaxies. More generally, our analysis also shows that the drop in the Lya fraction is quantitatively consistent with the observed evolution in the Lya luminosity functions of Lya Emitters.

We study the role of feedback from supernovae and black holes in the evolution of the star formation rate function (SFRF) of z~4-7 galaxies. We use a new set of cosmological hydrodynamic simulations, ANGUS (AustraliaN GADGET-3 early Universe Simulations), run with a modified and improved version of the parallel TreePM-smoothed particle hydrodynamics code GADGET-3 called P-GADGET3(XXL), that includes a self-consistent implementation of stellar evolution and metal enrichment. In our simulations both Supernova (SN) driven galactic winds and Active Galactic Nuclei (AGN) act simultaneously in a complex interplay. The SFRF is insensitive to feedback prescription at z>5, meaning that it cannot be used to discriminate between feedback models during reionisation. However, the SFRF is sensitive to the details of feedback prescription at lower redshift. By exploring different SN driven wind velocities and regimes for the AGN feedback, we find that the key factor for reproducing the observed SFRFs is a combination of "strong" SN winds and early AGN feedback in low mass galaxies. Conversely, we show that the choice of initial mass function and inclusion of metal cooling have less impact on the evolution of the SFRF. When variable winds are considered, we find that a non-aggressive wind scaling is needed to reproduce the SFRFs at z>4. Otherwise, the amount of objects with low SFRs is greatly suppressed and at the same time winds are not effective enough in the most massive systems.

We show that a bursty model of high redshift star formation explains several puzzling observations of the high redshift galaxy population. We begin by pointing out that the observed specific star formation rate requires a duty-cycle of ~10%, which is much lower than found in many hydro-dynamical simulations. This value follows directly from the fact that the observed star formation rate in galaxies integrated over a Hubble time exceeds the observed stellar mass by an order of magnitude. We use the large observed specific star formation rate to calibrate the efficiency of feedback in a model for the high redshift star formation rate which includes merger driven star formation regulated by SNe feedback. This model reproduces the star formation rate density function and the stellar mass function of galaxies at 4<z<7. A prediction of the model is that the specific star formation rate does not evolve very rapidly with either mass or redshift, in agreement with observation. This is in contrast to results from hydrodynamical simulations where the star formation closely follows the accretion rate, and so increases strongly towards high redshift. The bursty star formation model naturally explains the observation that there is not enough stellar mass at z~2-4 to account for all of the star-formation observed, without invoking properties like an evolving initial mass function of stars. The finding of a duty cycle that is ~10% implies that there should be ten times the number of known galaxies at fixed stellar mass that have not yet been detected through standard UV selection at high redshift. We therefore predict the existence of a large undetected population of UV-faint galaxies that accounts for most of the stellar mass density at z~4-8.

Motivated by recent observations of the star formation rate density function out to z~7, we describe a simple model for the star formation rate density function at high redshift based on the extended Press-Schechter formalism. This model postulates a starburst following each major merger, lasting for a time t_SF and converting at most f_star of galactic gas into stars. We include a simple physical prescription for supernovae feedback that suppresses star formation in low mass galaxies. Constraining t_SF and f_star to describe the observed star formation rate density at high redshifts, we find that individual starbursts were terminated after a time t_SF~10^7 years. This is comparable to the main-sequence lifetimes of supernova progenitors, indicating that high redshift starbursts are quenched once supernovae feedback had time to develop. High redshift galaxies convert ~10% of their mass into stars for galaxies with star formation rates above ~1 solar mass per year, but a smaller fraction for lower luminosity galaxies. Our best fit model successfully predicts the observed relation between star formation rate and stellar mass at z>~4, while our deduced relation between stellar mass and halo mass is also consistent with data on the dwarf satellites of the Milky Way. We find that supernovae feedback lowers the efficiency of star formation in the lowest mass galaxies and makes their contribution to reionization small. As a result, photo-ionization feedback on low mass galaxy formation does not significantly affect the reionization history. Using a semi-analytic model for the reionization history, we infer that approximately half of the ionizing photons needed to complete reionization have already been observed in star-forming galaxies.

We introduce the Fast Holographic Deconvolution method for analyzing interferometric radio data. Our new method is an extension of A-projection/software-holography/forward modeling analysis techniques and shares their precision deconvolution and widefield polarimetry, while being significantly faster than current implementations that use full direction-dependent antenna gains. Using data from the MWA 32 antenna prototype, we demonstrate the effectiveness and precision of our new algorithm. Fast Holographic Deconvolution may be particularly important for upcoming 21 cm cosmology observations of the Epoch of Reionization and Dark Energy where foreground subtraction is intimately related to the precision of the data reduction.

We study the visibility of the Lyman Alpha (Lya) emission line during the Epoch of Reionization (EoR). Combining galactic outflow models with large-scale semi-numeric simulations of reionization, we quantify the probability distribution function (PDF) of the fraction of Lya photons transmitted through the intergalactic medium (IGM), T_IGM. Our study focusses on galaxies populating dark matter halos with masses of M_halo=1e10 M_sun at z=8.6, which is inspired by the recent reported discovery of a galaxy at z=8.6 with strong Lya line emission. For reasonable model assumptions, we find that winds cause T_IGM>10% [50%], for the majority of galaxies, even when the Universe is ~80% [60%] neutral by volume. Thus, the observed strong Lya emission from the reported z=8.6 galaxy is consistent with a highly neutral IGM. We also investigate the implications of the recent tentative evidence for a observed decrease in the `LAE fraction' among drop-out galaxies between z=6 and z=7. If confirmed, we show that a rapid evolution in x_HI will be required to explain this observation via the effects of reionization.

Lyman Alpha Emitters (LAEs) are galaxies that have been selected on the basis of a strong Ly{\alpha} emission line in their spectra. Observational campaigns have dramatically increased the sample of known LAEs, which now extends out to z=7. These discoveries have motivated numerous theoretical studies on the subject, which usually define LAEs in their models based on sharp Ly{\alpha} luminosity and equivalent width (EW) cuts. While broadly representative, this procedure does not mimic the selection from observational programs in detail, which instead use cuts in various colour-spaces. We investigate what implications this disjoint may have for studies that aim to model LAEs. We construct an empirical model for the number density of star forming galaxies as a function of their UV and Ly{\alpha} luminosity, utilising measured constraints on the luminosity functions (LFs) of drop-out galaxies, and their luminosity dependent probability distribution function of Ly{\alpha} EW. In particular, we investigate whether the LAE LFs can be reproduced by defining LAEs using a (z-dependent) Ly{\alpha} luminosity and EW threshold. While we are able to reproduce the observed distribution of Ly{\alpha} EW among LAEs out to restframe EW 200 A, we find that our formalism over-predicts both the UV and Ly{\alpha} LFs of LAEs by a factor of 2-3, and is inconsistent with observations at the ~95% level. This tension is partially resolved if we assume the Ly{\alpha} EW-distribution of drop-out galaxies to be truncated at restframe EW>150 A. However the overprediction indicates that modeling LAEs with simple REW and luminosity cuts does not accurately mimic observed selection criteria, and can lead to uncertainties in the predicted number density of LAEs. On the other hand, the predicted z-evolution is not affected. We apply our formalism to drop-out galaxies at z>6, and predict the LFs of LAEs at z=7-9.

We identify a physical mechanism that would have resulted in rapid, obscured growth of seed super-massive black-holes in galaxies at z>6. Specifically, we find that the density at the centre of typical high redshift galaxies was at a level where the Bondi accretion rate implies a diffusion speed of photons that was slower than the gravitational infall velocity, resulting in photons being trapped within the accretion flow and advected into the black-hole. We show that there is a range of black-hole masses (M_bh ~ 10^3-10^5 solar masses) where the accretion flow traps radiation, corresponding to black-holes that were massive enough to generate a photon trapping accretion flow, but small enough that their Bondi radii did not exceed the isothermal scale height of self-gravitating gas. Under these conditions we find that the accretion reaches levels far in excess of the Eddington rate. A prediction of this scenario is that X-ray number counts of active galactic nuclei at z>6 would exhibit a cutoff at the low luminosities corresponding to black-hole masses below ~10^5 solar masses. At low redshifts we find photon trapping to be unimportant because it could only occur in rare low spin halos, and would require black-hole masses in excess of expectations from the observed black-hole - halo mass relation. The super-Eddington growth of ~10^5 solar mass seed black-holes at high redshift may have provided a natural acceleration towards the growth of super-massive black-holes at z~6-7, less than a billion years after the Big Bang.

It is commonly thought that stars are responsible for reionizing the Universe. However, deep searches for star-forming galaxies during the epoch of reionization have not yet found sufficient galaxies to provide the necessary ionising flux. Fast accretion shocks associated with gravitational infall of baryons during the formation of galaxies have recently been proposed as an alternative method of generating the required ionising photons. In this Letter we assess the contribution to hydrogen reionization from shocked gas associated with gravitational infall. We find that shocks can ionize no more than a few percent of the cosmic hydrogen by z~6. However, the small fraction of ionizing radiation produced by fast accretion shocks would be significantly more biased than that associated with stars, leading to a modification of the luminosity weighted source clustering by ~10%. This modification of the bias may be measurable with future precision experiments utilising the redshifted 21cm line to study the distribution of hydrogen during the reionization era.

We show that the dependence of Ly-alpha absorption on environment leads to significant non-gravitational features in the redshift space power-spectrum of Ly-alpha selected galaxies. We derive a physically motivated fitting formula that can be included in clustering analyses, and use this to discuss the predicted features in the Ly-alpha galaxy power-spectrum based on detailed models in which Ly-alpha absorption is influenced by gas infall and/or by strong galactic outflows. We show that power-spectrum measurements could be used to study the astrophysics of the galaxy-IGM connection, and to measure the properties of outflows from star-forming galaxies. Applying the modified redshift space power-spectrum to a Ly-alpha survey with parameters corresponding to the planned Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), we find that the dependence of observed Ly-alpha flux on velocity gradient and ionising background may compromise the ability of Ly-alpha selected galaxy redshift surveys to constrain cosmology using information from the full power-spectrum. This is because the effects of fluctuating ionizing background and velocity gradients effect the shape of the observed power-spectrum in ways that are similar to the shape of the primordial power-spectrum and redshift space distortions respectively. We use the Alcock-Paczynski test to show that without prior knowledge of the details of Ly-alpha absorption in the IGM, the precision of line-of-sight and transverse distance measurements for HETDEX will be ~1.3-1.7%, decreased by a factor of ~1.5-2 relative to the best case precision of ~0.8% available in a traditional galaxy redshift survey. We specify the precision with which modelling of Ly-alpha radiative transfer must be understood in order for HETDEX to achieve distance measurements that are better than 1%.

We use data on the high-redshift evolution of the size distribution and luminosity function of galaxies to constrain the relationship between their star formation efficiency and starburst lifetime. Based on the derived scaling relations, we predict the angular sizes and average surface brightnesses of faint galaxies that will be discovered with JWST. We find that JWST will be able to resolve galaxies at the magnitude limit m<31 out to a redshift of z~14. The next generation of large ground-based telescopes will resolve all galaxies discovered with JWST, provided they are sufficiently clumpy to enable detection above the bright thermal sky. We combine our constraints with simple models for self regulation of star formation, and show that feedback from supernovae at redshifts z>3 is likely mediated through momentum transfer, with the starburst timescale set by the lifetime of the massive stars rather than the dynamical time in the host galactic disk.

The discovery of quasars with redshifts higher than six has prompted a great deal of discussion in the literature regarding the role of quasars, both as sources of reionization, and as probes of the ionization state of the IGM. However the extreme ultra-violet (EUV) spectral index cannot be measured directly for high redshift quasars owing to absorption at frequencies above the Lyman limit, and as a result, studies of the impact of quasars on the intergalactic medium during reionization must assume a spectral energy distribution in the extreme ultra-violet based on observations at low redshift, z<1. In this paper we use regions of high Ly-alpha transmission (near-zones) around the highest redshift quasars to measure the quasar EUV spectral index at z~6. We jointly fit the available observations for variation of near-zone size with both redshift and luminosity, and propose that the observed relation provides evidence for an EUV spectral index that varies with absolute magnitude in the high redshift quasar sample, becoming softer at higher luminosity. Using a large suite of detailed numerical simulations we find that the typical value of spectral index for a luminous quasar at z~6 is constrained to be alpha=1.3+/-0.4 for a specific luminosity of the form L\propto\nu^{-alpha}. We find the scatter in spectral index among individual quasars to be in the range ~0.75-1.25. These values are in agreement with direct observations at low redshift, and indicate that there has been no significant evolution in the EUV spectral index of quasars over 90% of cosmic time.

The next generation of telescopes aim to directly observe the first generation of galaxies that initiated the reionization process in our Universe. The Lyman Alpha (Lya) emission line is robustly predicted to be the most prominent intrinsic spectral feature of these galaxies, making it an ideal target to search for and study high redshift galaxies. Unfortunately the large Gunn-Peterson optical depth of the surrounding neutral intergalactic medium (IGM) is thought to render this line extremely difficult to detect prior to reionization. In this paper we demonstrate that the radiative transfer effects in the interstellar medium (ISM), which cause Lya flux to emerge from galaxies at frequencies where the Gunn-Peterson optical depth is reduced, can substantially enhance the prospects for detection of the Lya line at high redshift. In particular, scattering off outflows of interstellar HI gas can modify the Lya spectral line shape such that >5% of the emitted Lya radiation is transmitted directly to the observer, even through a fully neutral IGM. It may therefore be possible to directly observe `strong' Lya emission lines (EW > 50 Angstrom rest frame) from the highest redshift galaxies that reside in the smallest HII `bubbles' early in the reionization era with JWST. In addition, we show that outflows can boost the fraction of Lya radiation that is transmitted through the IGM during the latter stages of reionization, and even post-reionization. Coupled with the fact that the first generation of galaxies are thought to have very large intrinsic equivalent Lya equivalent widths, our results suggest that the search for galaxies in their redshifted Lya emission line can be competitive with the drop-out technique out to the highest redshifts that can be probed in the JWST era.

GRB afterglows offer a probe of the intergalactic medium out to high redshift which complements observations along more abundant quasar lines-of-sight. Although both quasars and GRB afterglows should provide a-priori random sight-lines through the intervening IGM, it has been observed that strong Mg-II absorbers are twice as likely to be found along sight-lines toward GRBs. Several proposals to reconcile this discrepancy have been put forward, but none has been found sufficient to explain the magnitude of the effect. In this paper we estimate the effect of gravitational lensing by galaxies and their surrounding mass distributions on the statistics of Mg-II absorption. We find that the multi-band magnification bias could be very strong in the spectroscopic GRB afterglow population and that gravitational lensing can explain the discrepancy in density of absorbers, for plausibly steep luminosity functions. The model makes the prediction that approximately 20%-60% of the spectroscopic afterglow sample (i.e. ~ 5-15 of 26 sources) would have been multiply imaged, and hence result in repeating bursts. We show that despite this large lensing fraction it is likely that none would yet have been identified by chance owing to the finite sky coverage of GRB searches. We predict that continued optical monitoring of the bright GRB afterglow locations in the months and years following the initial decay would lead to identification of lensed GRB afterglows. A confirmation of the lensing hypothesis would allow us to constrain the GRB luminosity function down to otherwise inaccessibly faint levels, with potential consequences for GRB models.

The spatial clustering properties of HI galaxies can be studied using the formalism of the halo occupation distribution (HOD). The resulting parameter constraints describe properties like gas richness verses environment. Unfortunately, clustering studies based on individual HI galaxies will be restricted to the local Universe for the foreseeable future, even with the deepest HI surveys. Here we discuss how clustering studies of the HI HOD could be extended to moderate redshift, through observations of fluctuations in the combined 21cm intensity of unresolved galaxies. In particular we make an analytic estimate for the clustering of HI in the HOD. Our joint goals are to estimate i) the amplitude of the signal, and ii) the sensitivity of telescopes like the Australian SKA Pathfinder to HOD parameters. We find that the power spectrum of redshifted 21cm intensity could be used to study the distribution of HI within dark matter halos at z>0.5 where individual galaxies cannot be detected. In addition to the HOD of HI, the amplitude of the 21cm power spectrum would also yield estimates of the cosmic HI content at epochs between the local Universe, and redshifts probed by damped Ly-alpha absorbers.

21cm emission from residual neutral hydrogen after the epoch of reionization can be used to trace the cosmological power spectrum of density fluctuations. Using a Fisher matrix formulation, we provide a detailed forecast of the constraints on cosmological parameters that are achievable with this probe. We consider two designs: a scaled-up version of the MWA observatory as well as a Fast Fourier Transform Telescope. We find that 21cm observations dedicated to post-reionization redshifts may yield significantly better constraints than next generation Cosmic Microwave Background (CMB) experiments. We find the constraints on $\Omega_\Lambda$, $\Omega_{\rm m}h^2$, and $\Omega_\nu h^2$ to be the strongest, each improved by at least an order of magnitude over the Planck CMB satellite alone for both designs. Our results do not depend as strongly on uncertainties in the astrophysics associated with the ionization of hydrogen as similar 21cm surveys during the epoch of reionization. However, we find that modulation of the 21cm power spectrum from the ionizing background could potentially degrade constraints on the spectral index of the primordial power spectrum and its running by more than an order of magnitude. Our results also depend strongly on the maximum wavenumber of the power spectrum which can be used due to non-linearities.

We use the large catalogues of haloes available for the Millennium Simulation to test whether recently merged haloes exhibit stronger large-scale clustering than other haloes of the same mass. This effect could help to understand the very strong clustering of quasars at high redshift. However, we find no statistically significant excess bias for recently merged haloes over the redshift range 2 < z < 5, with the most massive haloes showing an excess of at most ~5%. We also consider galaxies extracted from a semianalytic model built on the Millennium Simulation. At fixed stellar mass, we find an excess bias of ~ 20-30% for recently merged objects, decreasing with increasing stellar mass. The fact that recently-merged galaxies are found in systematically more massive haloes than other galaxies of the same stellar mass accounts for about half of this signal, and perhaps more for high-mass galaxies. The weak merger bias of massive systems suggests that objects of merger-driven nature, such as quasars, do not cluster significantly differently than other objects of the same characteristic mass. We discuss the implications of these results for the interpretation of clustering data with respect to quasar duty cycles, visibility times, and evolution in the black hole-host mass relation.

We perform an analysis of the spatial clustering properties of HI selected galaxies from the HI Parkes All Sky Survey (HIPASS) using the formalism of the halo occupation distribution (HOD). The resulting parameter constraints show that the fraction of satellite galaxies (i.e. galaxies which are not the central member of their host dark matter halo) among HIPASS galaxies is <20%, and that satellite galaxies are therefore less common in HIPASS than in optically selected galaxy redshift surveys. Moreover the lack of fingers-of-god in the redshift space correlation function of HIPASS galaxies may indicate that the HI rich satellites which do exist are found in group mass rather than cluster mass dark matter halos. We find a minimum halo mass for HIPASS galaxies at the peak of the redshift distribution of M~10^11 solar masses, and show that less than 10% of baryons in HIPASS galaxies are in the form of HI. Quantitative constraints on HOD models from HIPASS galaxies are limited by uncertainties introduced through the small survey volume. However our results imply that future deeper surveys will allow the distribution of HI with environment to be studied in detail via clustering of HI galaxies.

The fraction of ionizing photons that escape their host galaxies and so are able to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era and early galaxy formation. Studies of the reionization history normally suffer from a degeneracy between the unknown values for the efficiency with which high redshift galaxies turn mass into stars and the escape fraction of ionizing photons. Recent gamma-ray burst (GRB) measurements of the star formation rate density during reionization provide the first opportunity to break this degeneracy. We confront a semi-analytic model for reionization with the GRB-derived star formation rate, as well as observations of the Ly-alpha forest and the CMB. Assuming that UV photons produced in star-forming galaxies dominate the reionization process, we show that the escape fraction of ionizing photons from high redshift galaxies is ~5% [log(f_esc)=-1.35+/-0.15 (68%)] for our fiducial model. This value is reasonably stable against uncertainties in the modeling, including the implementation of radiative feedback, the possibility of an evolving escape fraction, and the unknown shape of the IMF, which in sum contribute ~0.2 dex of additional systematic uncertainty on the value of escape fraction.

We discuss the 21cm power spectrum (PS) following the completion of reionization. In contrast to the reionization era, this PS is proportional to the PS of mass density fluctuations, with only a small modulation due to fluctuations in the ionization field on scales larger than the mean-free-path of ionizing photons. We derive the form of this modulation, and demonstrate that its effect on the 21cm PS will be smaller than 1% for physically plausible models of damped Ly-alpha systems. In contrast to the 21cm PS observed prior to reionization, in which HII regions dominate the ionization structure, the simplicity of the 21cm PS after reionization will enhance its utility as a cosmological probe by removing the need to separate the PS into physical and astrophysical components. As a demonstration, we consider the Alcock-Paczynski test and show that the next generation of low-frequency arrays could measure the angular distortion of the PS at the percent level for z~3-5.

We incorporate a contribution to reionization from X-rays within analytic and semi-numerical simulations of the 21-cm signal arising from neutral hydrogen during the epoch of reionization. We explore the impact that X-ray ionizations have on the power spectrum (PS) of 21-cm fluctuations by varying both the average X-ray MFP and the fractional contribution of X-rays to reionization. In general, prior to the epoch when the intergalactic medium is dominated by ionized regions (H {\sevensize II} regions), X-ray-induced ionization enhances fluctuations on spatial scales smaller than the X-ray MFP, provided that X-ray heating does not strongly supress galaxy formation. Conversely, at later times when \H2 regions dominate, small-scale fluctuations in the 21-cm signal are suppressed by X-ray ionization. Our modelling also shows that the modification of the 21-cm signal due to the presence of X-rays is sensitive to the relative scales of the X-ray MFP, and the characteristic size of \H2 regions. We therefore find that X-rays imprint an epoch and scale-dependent signature on the 21-cm PS, whose prominence depends on fractional X-ray contribution. The degree of X-ray heating of the IGM also determines the extent to which these features can be discerned. We show that the MWA will have sufficient sensitivity to detect this modification of the PS, so long as the X-ray photon MFP falls within the range of scales over which the array is most sensitive ($\sim0.1$ Mpc$^{-1}$). In cases in which this MFP takes a much smaller value, an array with larger collecting area would be required.

Using the 21 cm line, observed all-sky and across the redshift range from 0 to 5, the large scale structure of the Universe can be mapped in three dimensions. This can be accomplished by studying specific intensity with resolution ~ 10 Mpc, rather than via the usual galaxy redshift survey. The data set can be analyzed to determine Baryon Acoustic Oscillation wavelengths, in order to address the question: 'What is the nature of Dark Energy?' In addition, the study of Large Scale Structure across this range addresses the questions: 'How does Gravity effect very large objects?' and 'What is the composition our Universe?' The same data set can be used to search for and catalog time variable and transient radio sources.

Recently, a very large clustering length has been measured for quasars at a redshift of z~4. In combination with the observed quasar luminosity function we assess the implications of this clustering for the relationship between quasar luminosity and dark matter halo mass. Our analysis allows for non-linearity and finite scatter in the relation between quasar luminosity and halo mass, as well as a luminosity dependent quasar lifetime. The additional novel ingredient in our modelling is the allowance for an excess in the observed bias over the underlying halo bias owing to the merger driven nature of quasar activity. We find that the observations of clustering and luminosity function can be explained only if both of the following conditions hold: (i) The luminosity to halo mass ratio increases with halo mass; (ii) The observed clustering amplitude is in excess of that expected solely from halo bias. The latter result is statistically significant at the 99% level. Taken together, the observations provide compelling evidence for merger driven quasar activity, with a black-hole growth that is limited by feedback. In difference from previous analyses, we show that there could be scatter in the luminosity halo mass relation of up to 1 dex, and that quasar clustering can not be used to estimate the quasar lifetime.

The earliest generation of stars and black holes must have established an early 'Lyman-Werner' background (LWB) at high redshift, prior to the epoch of reionization. Because of the long mean free path of photons with energies E<13.6 eV, the LWB was nearly uniform. However, some variation in the LWB is expected due to the discrete nature of the sources, and their highly clustered spatial distribution. In this paper, we compute the probability distribution function (PDF) of the LW flux that irradiates dark matter (DM) halos collapsing at high-redshift (z~10). Our model accounts for (i) the clustering of DM halos, (ii) Poisson fluctuations in the number of corresponding star forming galaxies, and (iii) scatter in the LW luminosity produced by halos of a given mass (calibrated using local observations). We find that > 99% of the DM halos are illuminated by a LW flux within a factor of 2 of the global mean value. However, a small fraction, ~1e-8 to 1e-6, of DM halos with virial temperatures above 1e4 K have a close luminous neighbor within < 10 kpc, and are exposed to a LW flux exceeding the global mean by a factor of > 20, or to J_(21,LW)> 1e3 (in units of 1e-21 erg/s/Hz/sr/cm^2). This large LW flux can photo--dissociate H_2 molecules in the gas collapsing due to atomic cooling in these halos, and prevent its further cooling and fragmentation. Such close halo pairs therefore provide possible sites in which primordial gas clouds collapse directly into massive black holes (M_BH~ 1e4 - 1e6 M_sun), and subsequently grow into supermassive (M_BH > 1e9 M_sun) black holes by z~6.

We investigate the impact of neutral hydrogen (HI) in galaxies on the statistics of 21-cm fluctuations using analytic and semi-numerical modelling. Following the reionisation of hydrogen the HI content of the Universe is dominated by damped absorption systems (DLAs), with a cosmic density in HI that is observed to be constant at a level equal to ~2% of the cosmic baryon density from z~1 to z~5. We show that extrapolation of this constant fraction into the reionisation epoch results in a reduction of 10-20% in the amplitude of 21-cm fluctuations over a range of spatial scales. The assumption of a different percentage during the reionisation era results in a proportional change in the 21-cm fluctuation amplitude. We find that consideration of HI in galaxies/DLAs reduces the prominence of the HII region induced shoulder in the 21-cm power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. We also estimate the 21cm-galaxy cross PS, and show that the cross PS changes sign on scales corresponding to the HII regions. From consideration of the sensitivity for forthcoming low-frequency arrays we find that the effects of HI in galaxies/DLAs on the statistics of 21-cm fluctuations will be significant with respect to the precision of a PS or cross PS measurement. In addition, since overdense regions are reionised first we demonstrate that the cross-correlation between galaxies and 21-cm emission changes sign at the end of the reionisation era, providing an alternative avenue to pinpoint the end of reionisation. The sum of our analysis indicates that the HI content of the galaxies that reionise the universe will need to be considered in detailed modelling of the 21-cm intensity PS in order to correctly interpret measurements from forthcoming low-frequency arrays.

The fraction of ionizing photons which escape their host galaxy and so are able to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in studies of the reionization era and early galaxy formation. In this paper we combine observations of Lyman-alpha absorption towards high redshift quasars with the measured UV luminosity function of high redshift galaxies to constrain the escape fraction (f_esc) of ionizing photons from galaxies at z ~ 5.5-6. The observed Lyman-alpha transmission constrains the escape fraction to lie in the range f_esc ~ 10-25 % (at z ~ 5.5-6). Excluding halos with M< 10^10 M_sun (as might be expected if galaxy formation is suppressed due to the reionization of the IGM) implies a larger escape fraction of f_esc ~ 20-45 %. Using the numerical results to calibrate an analytic relation between the escape fraction and minimum galaxy halo mass we also extrapolate our results to a mass (M~10^8 M_sun) corresponding to the hydrogen cooling threshold. In this case we find f_esc ~ 5-10 %, consistent with observed estimates at lower redshift. We find that the escape fraction of high redshift galaxies must be greater than 5 % irrespepctive of galaxy mass. Based on these results we use a semi-analytic description to model the reionization history of the IGM, assuming ionizing sources with escape fractions suggested by our numerical simulations. We find that the IBG observed at z ~ 5.5-6 implies a sufficient number of ionizing photons to have reionized the Universe by z ~ 6. However, if the minimum mass for star-formation were greater than 10^9 M_sun, the IBG would be over-produced at redshifts less than z ~ 5. In summary, our results support a scenario in which the IGM was reionized by low mass galaxies.

We assess the impact of Galactic synchrotron foreground removal on the observation of high-redshift quasar HII regions in redshifted 21-cm emission. We consider the case where a quasar is observed in an intergalactic medium (IGM) whose ionisation structure evolves slowly relative to the light crossing time of the HII region, as well as the case where the evolution is rapid. The latter case is expected towards the end of the reionisation era where the highest redshift luminous quasars will be observed. In the absence of foregrounds the fraction of neutral hydrogen in the IGM could be measured directly from the contrast between the HII region and surrounding IGM. However, we find that foreground removal lowers the observed contrast between the HII region and the IGM. This indicates that measurement of the neutral fraction would require modelling to correct for this systematic effect. On the other hand, foreground removal does not modify the most prominent features of the 21-cm maps. Using a simple algorithm we demonstrate that measurements of the size and shape of observed HII regions will not be affected by continuum foreground removal. Moreover, measurements of these quantities will not be adversely affected by the presence of a rapidly evolving IGM.

Observations of damped Ly-alpha absorbers (DLA) indicate that the fraction of hydrogen in its neutral form (HI) is significant by mass at all redshifts. This gas represents the reservoir of material that is available for star formation at late times. As a result, observational identification of the systems in which this neutral hydrogen resides is an important missing ingredient in models of galaxy formation. Precise identification of DLA host mass via traditional clustering studies is not practical owing to the small numbers of known systems being spread across sparsely distributed sight lines. However following the completion of reionization, 21cm surface brightness fluctuations will be dominated by neutral hydrogen in DLAs. Observations of these fluctuations will measure the combined clustering signal from all DLAs within a large volume. We show that measurement of the spherically averaged power-spectrum of 21cm intensity fluctuations due to DLAs could be used to measure the galaxy bias for DLA host galaxies when combined with an independent measurement of the cosmological HI mass density from quasar absorption studies. Utilising this technique, the low frequency arrays now under construction could measure the characteristic DLA host mass with a statistical precision as low as 0.3 dex at z~4. In addition, high signal-to-noise observations of the peculiar-motion induced anisotropy of the power-spectrum would facilitate measurement of both the DLA host mass and the cosmic HI density directly from 21cm fluctuations. By exploiting this anisotropy, a second generation of low frequency arrays with an order of magnitude increase in collecting area could measure the values of cosmic HI density and DLA host mass, with uncertainties of a few percent and a few tens of percent respectively.

The introduction of low-frequency radio arrays is expected to revolutionize the study of the reionization epoch. Observation of the contrast in redshifted 21cm emission between a large HII region and the surrounding neutral inter-galactic medium (IGM) will be the simplest and most easily interpreted signature. However the highest redshift quasars known are thought to reside in an ionized IGM. Using a semi-analytic model we describe the redshifted 21cm signal from the IGM surrounding quasars discovered using the i-drop out technique (i.e. quasars at z~6). We argue that while quasars at z<6.5 seem to reside in the post overlap IGM, they will still provide valuable probes of the late stages of the overlap era because the light-travel time across a quasar proximity zone should be comparable to the duration of overlap. For redshifted 21cm observations within a 32MHz bandpass, we find that the subtraction of a spectrally smooth foreground will not remove spectral features due to the proximity zone. These features could be used to measure the neutral hydrogen content of the IGM during the late stages of reionization. The density of quasars at z~6 is now well constrained. We use the measured quasar luminosity function to estimate the prospects for discovery of high redshift quasars in fields that will be observed by the Murchison Widefield Array.

Measurements of the 21cm line emission by residual cosmic hydrogen after reionization can be used to trace the power spectrum of density perturbations through a significant fraction of the observable volume of the Universe. We show that a dedicated 21cm observatory coule probe a number of independent modes that is two orders of magnitude larger than currently available, and enable a cosmic-variance limited detection of the signature of a neutrino mass ~0.05eV. The evolution of the linear growth factor with redshift could also constrain exotic theories of gravity or dark energy to an unprecedented precision.

We assess the impact of inhomogeneous reionization on detection of HII regions surrounding luminous high redshift quasars using planned low frequency radio telescopes. Our approach is to implement a semi-numerical scheme to calculate the 3-dimensional structure of ionized regions surrounding a massive halo at high redshift, including the ionizing influence of a luminous quasar. As part of our analysis we briefly contrast our scheme with published semi-numerical models. We calculate mock 21cm spectra along the line of sight towards high redshift quasars, and estimate the ability of the planned Murchison Widefield Array to detect the presence of HII regions. The signal-to-noise for detection will drop as the characteristic bubble size grows during reionization because the quasar's influence becomes less prominent. However, quasars will imprint a detectable signature on observed 21cm spectra that is distinct from a region of typical IGM. At epochs where the mean hydrogen neutral fraction is ~30% or greater we find that neutral gas in the IGM surrounding a single quasar will be detectable (at a significance of 5 sigma) within 100 hour integrations in more than 50% of cases. 1000 hour integrations will be required to detect a smaller neutral fraction of 15% in more than 50% of cases. A highly significant detection will be possible in only 100 hours for a stack of 10 smaller 3 proper Mpc HII regions. The accurate measurement of the global average neutral fraction (<x_HI>) will be limited by systematic fluctuations between lines of sight for single HII regions. We estimate the accuracy with which the global neutral fraction could be measured from a single HII region to be 50%, 30% and 20% for <x_HI> ~ 0.15, 0.3 and 0.5 respectively.

Low-frequency observatories are currently being constructed with the goal of detecting redshifted 21cm emission from the epoch of reionization. These observatories will also be able to detect intensity fluctuations in the cumulative 21cm emission after reionization, from hydrogen in unresolved damped Ly-alpha absorbers (such as gas rich galaxies) down to a redshift z~3.5. The inferred power spectrum of 21cm fluctuations at all redshifts will show acoustic oscillations, whose co-moving scale can be used as a standard ruler to infer the evolution of the equation of state for the dark energy. We find that the first generation of low-frequency experiments (such as MWA or LOFAR) will be able to constrain the acoustic scale to within a few percent in a redshift window just prior to the end of the reionization era, provided that foregrounds can be removed over frequency band-passes of >8MHz. This sensitivity to the acoustic scale is comparable to the best current measurements from galaxy redshift surveys, but at much higher redshifts. Future extensions of the first generation experiments (involving an order of magnitude increase in the antennae number of the MWA) could reach sensitivities below one percent in several redshift windows and could be used to study the dark energy in the unexplored redshift regime of 3.5<z<12. Moreover, new experiments with antennae designed to operate at higher frequencies would allow precision measurements (<1%) of the acoustic peak to be made at more moderate redshifts (1.5<z<3.5), where they would be competitive with ambitious spectroscopic galaxy surveys covering more than 1000 square degrees. Together with other data sets, observations of 21cm fluctuations will allow full coverage of the acoustic scale from the present time out to z~12.

The fluctuations in the emission of redshifted 21cm photons from neutral inter-galactic hydrogen will provide an unprecedented probe of the reionization era. Conventional wisdom assumes that this 21cm signal disappears as soon as reionization is complete, when little atomic hydrogen is left through most of the volume of the inter-galactic medium (IGM). However observations of damped Ly-alpha absorbers indicate that the fraction of hydrogen in its neutral form is significant by mass at all redshifts. Here we use a physically-motivated model to show that residual neutral gas, confined to dense regions in the IGM with a high recombination rate, will generate a significant post-reionization 21cm signal. We show that the power-spectrum of fluctuations in this signal will be detectable by the first generation of low-frequency observatories at a signal-to-noise that is comparable to that achievable in observations of the reionization era. The statistics of 21cm fluctuations will therefore probe not only the pre-reionization IGM, but rather the entire process of HII region overlap, as well as the appearance of the diffuse ionized IGM.

Absorption spectra of high redshift quasars exhibit an increasingly thick Ly-alpha forest towards z~6. However, the interpretation of these spectra is complicated by the fact that the Ly-alpha optical depth is already large for neutral hydrogen fractions in excess of 10^-4, and also because quasars are expected to reside in dense regions of the IGM. We present a model for the evolution of the ionization state of the IGM which is applicable to the dense, biased regions around high-redshift quasars as well as more typical regions in the IGM, and combine this with numerical radiative transfer simulations. Our model is able to simultaneously reproduce the observed Ly-alpha forest opacity at 4<z<6, the ionizing photon mean-free-path at z~4 and the rapid evolution of highly ionized near-zone sizes around high-redshift quasars at 5.8<z<6.4. We find that within 5 physical Mpc of a high redshift quasar, the evolution of the ionization state of the IGM precedes that in more typical regions by around 0.3 redshift units. More importantly, when combined with the rapid increase in the ionizing photon mean-free-path expected shortly after overlap, this offset results in an ionizing background near the quasar which exceeds the value in the rest of the IGM by a factor of ~2-3. We further find that in the post-overlap phase of reionization the size of the observed quasar near-zones is not directly sensitive to the neutral hydrogen fraction of the IGM. Instead, these sizes probe the level of the background ionization rate and the temperature of the surrounding IGM. The observed rapid evolution of the quasar near-zone sizes at 5.8<z<6.4 can thus be explained by the rapid evolution of the ionizing background, which in our model is caused by the completion of overlap at the end of reionization by 6<z<7.

Several gravitationally lensed quasars are observed with anomalous magnifications in pairs of images that straddle a critical curve. Simple theoretical arguments suggest that the magnification of these images should be approximately equivalent, whereas one image is observed to be significantly demagnified. Microlensing provides a possible explanation for this discrepancy. There are two key parameters when modelling this effect. The first, the fraction of smooth matter in the lens at the image positions, has been explored by Schechter and Wambsganss (2002). They have shown that the anomalous flux ratio observed in the lensed quasar MG 0414+0534 is a priori a factor of 5 more likely if the assumed smooth matter content in the lens model is increased from 0% to 93%. The second parameter, the size of the emission region, is explored in this paper, and shown to be more significant. We find that the broadening of the magnification probability distributions due to smooth matter content is washed out for source sizes that are predicted by standard models for quasars. We apply our model to the anomalous lensed quasar MG 0414+0534, and find a 95% upper limit of 2.62 x 10^(16) h^(-1/2) (M/Msun)^(1/2) cm on the radius of the I-band emission region. The smooth matter percentage in the lens is unconstrained.

A significant fraction of Lyman Alpha (Lya) emitting galaxies (LAEs) at z> 5.7 have rest-frame equivalent widths (EW) greater than ~100 Angstrom. However only a small fraction of the Lya flux produced by a galaxy is transmitted through the IGM, which implies intrinsic Lya EWs that are in excess of the maximum allowed for a population-II stellar population having a Salpeter mass function. In this paper we study characteristics of the sources powering Lya emission in high redshift galaxies. We propose a simple model for Lya emitters in which galaxies undergo a burst of very massive star formation that results in a large intrinsic EW, followed by a phase of population-II star formation with a lower EW. We confront this model with a range of high redshift observations and find that the model is able to simultaneously describe the following eight properties of the high redshift galaxy population with plausible values for parameters like the efficiency and duration of star formation: i-iv) the UV and Lya luminosity functions of LAEs at z=5.7 and 6.5, v-vi) the mean and variance of the EW distribution of Lya selected galaxies at z=5.7, vii) the EW distribution of i-drop galaxies at z~6, and viii) the observed correlation of stellar age with EW. Our modeling suggests that the observed anomalously large intrinsic equivalent widths require a burst of very massive star formation lasting no more than a few to ten percent of the galaxies star forming lifetime. This very massive star formation may indicate the presence of population-III star formation in a few per cent of i-drop galaxies, and in about half of the Lya selected galaxies.

We consider the effect of reionization on the clustering properties of galaxy samples at intermediate redshifts (z~0.3-5.5). Current models for the reionization of intergalactic hydrogen predict that overdense regions will be reionized early, thus delaying the build up of stellar mass in the progenitors of massive lower-redshift galaxies. As a result, the stellar populations observed in intermediate redshift galaxies are somewhat younger and hence brighter in overdense regions of the Universe. Galaxy surveys would therefore be sensitive to galaxies with a somewhat lower dark matter mass in overdense regions. The corresponding increase in the observed number density of galaxies can be parameterized as a galaxy bias due to reionization. We model this process using merger trees combined with a stellar synthesis code. Our model demonstrates that reionization has a significant effect on the clustering properties of galaxy samples that are selected based on their star-formation properties. The bias correction in Lyman-break galaxies (including those in proposed baryonic oscillation surveys at z<1) is at the level of 10-20% for a halo mass of 10^12 solar masses, leading to corrections factors of 1.5-2 in the halo mass inferred from measurements of clustering length. The reionization of helium could also lead to a sharp increase in the amplitude of the galaxy correlation function at z~3. We find that the reionization bias is approximately independent of scale and halo mass. However since the traditional galaxy bias is mass dependent, the reionization bias becomes relatively more important for lower mass systems. The correction to the bias due to reionization is very small in surveys of luminous red galaxies at z<1.

Recent observations imply that the observed number counts of Lya Emitters (LAEs) evolved significantly between z=5.7 and z=6.5. It has been suggested that this was due to a rapid evolution in the ionisation state, and hence transmission of the IGM which caused Lya flux from z=6.5 galaxies to be more strongly suppressed. In this paper we consider the joint evolution of the Lya and UV luminosity functions (LFs) and show that the IGM transmission evolved between z=6.5 and z=5.7 by a factor 1.1 <R < 1.8 (95% CL). This result is insensitive to the underlying model of the Lya LF (as well as cosmic variance). Using a model for IGM transmission, we find that the evolution of the mean IGM density through cosmic expansion alone may result in a value for the ratio of transmissions as high as R=1.3. Thus, the existing LFs do not provide evidence for overlap. Furthermore, the constraint R<1.8 suggests that the Universe at z=6.5 was more than half ionised by volume, i.e. x_i,V>0.5.

Reionization is thought to be dominated by low mass galaxies, while direct observations of resolved galaxies probe only the most massive, rarest objects. The cross-correlation between fluctuations in the surface brightness of the cumulative Ly-alpha emission (which serves as a proxy for the star formation rate) and the redshifted 21cm signal from neutral hydrogen in the intergalactic medium (IGM), will directly probe the causal link between the production of ionizing photons in galaxies and the reionization of the IGM. We discuss the prospects for detecting this cross-correlation for unresolved galaxies. We find that on angular scales <10 arc-minutes, detection will be practical using widefield near-IR imaging from space in combination with the forthcoming Mileura Widefield Array - Low Frequency Demonstrator. When redshifted 21cm observations of the neutral IGM are combined with space-based near-IR imaging of Ly-alpha emission, the detection on angular scales <3 arc-minutes will be limited by the sensitivity of the 21cm signal, even when a small aperture optical telescope (~2m) and a moderate field of view (~10 square degrees) are used. On scales >3 arc-minutes, the measurement of cross-correlation will be limited by the accuracy of the foreground sky subtraction.

We calculate the impact of the intergalactic medium (IGM) on the observed Lyman alpha lines (hereafter Lya) emitted by galaxies in an ionised IGM at z>4. Our model accounts for gas clumping in the IGM and for the fact that high-redshift galaxies reside in overdense regions, which causes the velocity field of the IGM to depart from the Hubble flow. The observed shape of the Lya line varies widely, with dependence on the intrinsic width and systemic velocity of the line, a galaxies star formation rate and the local extra-galactic UV-background. For large star formation rates and levels of the UV-background, absorption in the IGM does not result in a Lya line that is asymmetric as is common among known high-redshift Lya emitters. For models in which the lines do show the observed strong asymmetry, the IGM typically transmits only 10-30% of the Lya flux. The increase in the ionising background that accompanied the completion of reionisation barely increased the IGM transmission, which suggests that LAEs of comparable luminosity should not appear to be significantly dimmer prior to overlap. In this light, we briefly discuss the potential of Lya emitters as a probe into the epoch of reionisation.

We have searched for HI 21cm absorption toward the two brightest radio AGN at high redshift, J0924--2201 at $z = 5.20$, and J0913+5919 at $z=5.11$, using the Giant Meter Wave Radio Telescope (GMRT). These data set a 3$\sigma$ upper limit to absorption of $< 30%$ at 40 km s$^{-1}$ resolution for the 30 mJy source J0913+5919, and $< 3%$ for the 0.55 Jy source J0924--2201 at 20 km s$^{-1}$ resolution. For J0924--2201, limits to broader lines at the few percent level are set by residual spectral baseline structure. For J0924--2201 the column density limit per 20 km s $^{-1}$ channel is: N(HI) $< 2.2\times 10^{18} \rm T_s$ cm$^{-2}$ over a velocity range of -700 km s$^{-1}$ to $+1180$ km s$^{-1}$ centered on the galaxy redshift determined through CO emission, assuming a covering factor of one. For J0913+5919 the column density limit per 40 km s$^{-1}$ channel is: N(HI) $< 2.2\times 10^{19} \rm T_s$ cm$^{-2}$ within $\pm 2400$ km s$^{-1}$ of the optical redshift. These data rule out any cool, high column density HI clouds within roughly $\pm 1000$ km s$^{-1}$ of the galaxies, as are often seen in Compact Steep Spectrum radio AGN, or clouds that might correspond to residual gas left over from cosmic reionization.

Spatial dependence in the statistics of redshifted 21cm fluctuations promises to provide the most powerful probe of the reionisation epoch. In this paper we consider the second and third moments of the redshifted 21cm intensity distribution using a simple model that accounts for galaxy bias during the reionisation process. We demonstrate that skewness in redshifted 21cm maps should be substantial throughout the reionisation epoch and on all angular scales, owing to the effects of galaxy bias which leads to early reionisation in over-dense regions of the IGM. The variance (or power-spectrum) of 21cm fluctuations will exhibit a minimum in redshift part way through the reionisation process, when the global ionisation fraction is around 50%. This minimum is generic, and is due to the transition from 21cm intensity being dominated by over-dense too under-dense regions as reionisation progresses. We show that the details of the reionisation history, including the presence of radiative feedback are encoded in the evolution of the auto-correlation and skewness functions with redshift and mean IGM neutral fraction. We discuss the detection of skewness by first generation instruments, and conclude that the Mileura Widefield Array - Low Frequency Demonstrator will have sufficient sensitivity to detect skewness on a range of angular scales at redshifts near the end of reionisation, while a subsequent instrument of 10 times the collecting area could map out the evolution of skewness in detail. The observation of a minimum in variance during the reionisation history, and the detection of skewness would both provide important confirmation of the cosmological origin of redshifted 21cm intensity fluctuations.

Absorption spectra of high redshift quasars suggest that the reionisation of cosmic hydrogen was complete near z~6. The dominant sources of ionising photons responsible for this reionisation are generally thought to be stars and quasars. In this paper we make a quantitative estimate of the relative contributions made by these sources. Our approach is to compute the evolution of the post overlap ionising background radiation by combining semi-analytic descriptions of reionisation in a clumpy medium with a model for the quasar luminosity function. Our overall model has two free parameters, the star formation efficiency and the minimum quasar luminosity. By adjusting these parameters, we constrain the relative contributions made by stars and quasars through comparison with reported observations (Fan et al. 2005). We find that the relative quasar contribution (at z=5.7) to the ionising background was between 1.4% and 14.5%. The range of uncertainty is dominated by the unknown minimum quasar luminosity.

Deep surveys have recently discovered galaxies at the tail end of the epoch of reionization. In the near future, these discoveries will be complemented by a new generation of low-frequency radio observatories that will map the distribution of neutral hydrogen in the intergalactic medium through its redshifted 21cm emission. In this paper we calculate the expected cross-correlation between the distribution of galaxies and intergalactic 21cm emission at high redshifts. We demonstrate using a simple model that overdense regions are expected to be ionized early as a result of their biased galaxy formation. This early phase leads to an anti-correlation between the 21cm emission and the overdensities in galaxies, matter, and neutral hydrogen. Existing Ly-alpha surveys probe galaxies that are highly clustered in overdense regions. By comparing 21cm emission from regions near observed galaxies to those away from observed galaxies, future observations will be able to test this generic prediction and calibrate the ionizing luminosity of high-redshift galaxies.

Recent surveys have detected Lya emission from z=4.5-6.5 at luminosities as low as 10^41 erg/s. There is good evidence that low numbers of AGN are among observed faint Lya emitters. Combining these observations with an empirical relation between the intrinsic Lya and B-band luminosities of AGN, we obtain an upper limit on the number density of AGN with absolute magnitudes M_B=[-16,-19] at z=4.5-6.5. These AGN are up to two orders of magnitude fainter than those discovered in the Chandra Deep Field, resulting in the faintest observational constraints to date at these redshifts. At z=4.5, the powerlaw slope of the very faint end of the luminosity function of AGN is shallower than the slope observed at lower redshifts, beta <1.6, at the 98% confidence level. In fact, we find marginal evidence that the luminosity function rises with luminosity, corresponding to a powerlaw slope beta <0, at magnitudes fainter than M_B~-20 (75% confidence level). These results suggest either that accretion onto lower mass black holes is less efficient than onto their more massive counterparts, or that the number of black holes powering AGN with M_B >-20 is lower than expected from the M_BH-sigma relation by one-two orders of magnitude. Extrapolating from reverberation-mapping studies suggests that these black holes would have M_BH=10^6-10^7 Msun. To facilitate the identification of AGN among observed Lya emitters, we derive observational properties of faint AGN in the Lya line, as well as in the X-ray and optical bands.

We demonstrate that a log-linear relation does not provide an adequate description of the correlation between the masses of Super-Massive Black-Holes (SMBH, M_bh) and the velocity dispersions of their host spheroid (sigma). An unknown relation between log(M_bh) and log(sigma) may be expanded to second order to obtain a log-quadratic relation of the form log(M_bh)=alpha+beta log(sigma/200) + beta_2[log(sigma/200)]^2. We perform a Bayesian analysis using the Nuker sample, and solve for beta, beta_2 and alpha, in addition to the intrinsic scatter (delta). We find unbiased parameter estimates of beta=4.2+/-0.37, beta_2=1.6+/-1.3 and delta=0.275+/-0.05. At the 80% level the M_bh-sigma relation does not follow a uniform power-law. Indeed, over the velocity range 70km/s<sigma<380km/s the logarithmic slope of the best fit relation varies between 2.7 and 5.1, which should be compared with a power-law estimate of 4.02+/-0.33. Assuming no systematic offset, single epoch virial SMBH masses estimated for AGN follow the same log-quadratic M_bh-sigma relation as the Nuker sample, but extend it downward in mass by an order of magnitude. The log-quadratic term in the M_bh-sigma relation has a significant effect on estimates of the local SMBH mass function at M_bh>10^9 solar masses, leading to densities of SMBHs with M_bh>10^10 solar masses that are several orders of magnitude larger than inferred from a log-linear relation. We also estimate unbiased parameters for the SMBH-bulge mass relation. With a parameterisation log(M_bh)=alpha_b + beta_b log(M_b/10^{11}) + beta_2b[log(M_b/10^{11})]^2, we find beta_b=1.15+/-0.18 and beta_2b=0.12+/-0.14. We determined an intrinsic scatter delta_b=0.41+/-0.07 which is ~50% larger than the scatter in the M_bh-sigma relation.

The HII regions around quasars and galaxies at redshifts beyond the epoch of reionisation will provide prime targets for upcoming 21cm campaigns using a new generation of low-frequency radio observatories. Here we show that the boundaries of these HII regions will not be sharp. Rather, the clustering of sources near massive galaxies results in a neutral fraction that rises gradually towards large radii from an interior value near zero. A neutral fraction corresponding to the global background value is typically reached at a distance of 2-5 times the radius of the HII region around the central massive galaxy.

A large number of faint galaxies, born less than a billion years after the big bang, have recently been discovered. The fluctuations in the distribution of these galaxies contributed to a scatter in the ionization fraction of cosmic hydrogen on scales of tens of Mpc, as observed along the lines of sight to the earliest known quasars. Theoretical simulations predict that the formation of dwarf galaxies should have been suppressed after cosmic hydrogen was reionized, leading to a drop in the cosmic star formation rate. Here we present evidence for this suppression. We show that the post-reionization galaxies which produced most of the ionizing radiation at a redshift z~5.5, must have had a mass in excess of ~10^{10.6+/-0.4} solar masses or else the aforementioned scatter would have been smaller than observed. This limiting mass is two orders of magnitude larger than the galaxy mass that is thought to have dominated the reionization of cosmic hydrogen (~10^8 solar masses). We predict that future surveys with space-based infrared telescopes will detect a population of smaller galaxies that reionized the Universe at an earlier time, prior to the epoch of dwarf galaxy suppression.

In the local universe, the masses of Super-Massive Black-Holes (SMBH) appear to correlate with the physical properties of their hosts, including the mass of the dark-matter halo. Using these clues as a starting point many studies have produced models that can explain phenomena like the quasar luminosity function. The shortcoming of this approach is that working models are not unique, and as a result it is not always clear what input physics is being constrained. Here we take a different approach. We identify critical parameters that describe the evolution of SMBHs at high redshift, and constrain their parameter space based on observations of high redshift quasars from the Sloan Digital Sky Survey. We find that the luminosity function taken in isolation is somewhat limited in its ability to constrain SMBH evolution due to some strong degeneracies. This explains the presence in the literature of a range of equally successful models based on different physical hypotheses. Including the constraint of the local SMBH to halo mass ratio breaks some of the degeneracies, and our results suggest halo masses at z~4.8 of 10^{12.5+/-0.3}M_solar (with 90% confidence), with a SMBH to halo mass ratio that decreases with time (>99%). We also find a quasar luminosity to halo mass ratio that increases with halo mass (>99%). These features need to be incorporated in all successful models of SMBH evolution. On the other hand current observations do not permit any conclusions regarding the evolution of quasar lifetime, or the SMBH occupation fraction in dark matter halos.

Population-III (Pop-III) starformation (SF) is thought to be quenched when the metallicity of the star-forming gas reaches a critical level. At high z, when the general intergalactic medium (IGM) was enriched with metals, the fraction of primordial gas already collapsed in minihalos was significantly larger than the fraction of primordial gas that had already been involved in Pop-III SF. We argue that this minihalo gas remained largely in a metal-free state, until these minihalos merged into large systems and formed stars. As a result, the era of Pop-III SF was significantly prolonged, leading to an integrated Pop-III SF an order of magnitude larger than expected for an abrupt transition redshift. The contribution of Pop-III SF to the reionization of hydrogen could have been significant until z~10 and may have extended to z~6. Our modeling allows for gradual enrichment of the IGM, feedback from photo-ionization and screening of reionization by minihalos. Nevertheless, extended Pop-III SF can result in complex, multi-peaked reionization histories. The contribution of Pop-III stars to reionization will be tested by the three-year WMAP results: (1) if Pop-III stars do not contribute to reionization, tau_es<0.05-0.06 and a rapid reionization at z~6 is expected; (2) if the product of star formation efficiency and escape fraction for Pop-III stars is significantly larger than for Pop-II stars, then a maximum tau_es=0.21 is achievable; (3) in a scenario where the product of star formation efficiency and escape fraction for Pop-III stars is comparable to that for Pop-II stars, tau_es=0.09-0.12 would be observed, with reionization histories characterized by an extended ionization plateau from z=7-12. This result holds regardless of the redshift where the IGM becomes enriched with metals.

Following the completion of cosmic reionization, the mean-free-path of ionizing photons was set by a population of Ly-limit absorbers. As the mean-free-path steadily grew, the intensity of the ionizing background also grew, thus lowering the residual neutral fraction of hydrogen in ionization equilibrium throughout the diffuse intergalactic medium (IGM). Ly-alpha photons provide a sensitive probe for tracing the distribution of this residual hydrogen at the end of reionization. Here we calculate the cosmic variance among different lines-of-sight in the distribution of the mean Ly-alpha optical depths. We find fractional variations in the effective post-reionization optical depth that are of order unity on a scale of ~100 co-moving Mpc, in agreement with observations towards high-redshift quasars. Significant contributions to these variations are provided by the cosmic variance in the density contrast on the scale of the mean-free-path for ionizing photons, and by fluctuations in the ionizing background induced by delayed or enhanced structure formation. Cosmic variance results in a highly asymmetric distribution of transmission through the IGM, with fractional fluctuations in Ly-alpha transmission that ar larger than in Ly-beta transmission.

We show that the observed scatter in the relations between the mass of supermassive black holes (SMBHs), M_bh, and the velocity dispersion sigma or mass M_sph of their host spheroid, place interesting constraints on the process that regulates SMBH growth in galaxies. When combined with the observed properties of early-type SDSS galaxies, the observed intrinsic scatters imply that SMBH growth is regulated by the spheroid velocity dispersion rather than its mass. The M_bh-M_sph relation is therefore a by-product of a more fundamental M_bh-sigma relation. We construct a theoretical model for the scatter among baryon modified dark matter halo profiles, out of which we generate a population of spheroid hosts and show that these naturally lead to a relation between effective radius and velocity dispersion of the form R_sph ~ sigma^1.5 with a scatter of ~0.2dex, in agreement with the corresponding projection of the fundamental plane for early type galaxies in SDSS. At the redshift of formation, our model predicts the minimum scatter that SMBHs can have at fixed velocity dispersion or spheroid mass under different formation scenarios. We also estimate the additional scatter that is introduced into these relations through collisionless mergers of purely stellar spheroids at z<1. We find that the observed scatter in the M_bh-sigma and M_bh-M_sph relations preclude the properties of dark matter halos from being the governing factor in SMBH growth. Finally, we show that SMBH growth governed by the properties of the host spheroid can lead to the observed values of scatter in the M_bh-sigma and M_bh-M_sph relations, only if the SMBH growth is limited by momentum or energy feedback over the dynamical time of the host spheroid.

The introduction of low-frequency radio arrays over the coming decade is expected to revolutionize the study of the reionization epoch. Observation of the contrast in redshifted 21cm emission between a large HII region and the surrounding neutral IGM will be the simplest and most easily interpreted signature. We find that an instrument like the planned Mileura Widefield Array Low-Frequency Demonstrator (LFD) will be able to obtain good signal to noise on HII regions around the most luminous quasars, and determine some gross geometric properties, e.g. whether the HII region is spherical or conical. A hypothetical follow-up instrument with 10 times the collecting area of the LFD (MWA-5000) will be capable of mapping the detailed geometry of HII regions, while SKA will be capable of detecting very narrow spectral features as well as the sharpness of the HII region boundary. The MWA-5000 will discover serendipitous HII regions in widefield observations. We estimate the number of HII regions which are expected to be generated by quasars. Assuming a late reionization at z~6 we find that there should be several tens of quasar HII regions larger than 4Mpc at z~6-8 per field of view. Identification of HII regions in forthcoming 21cm surveys can guide a search for bright galaxies in the middle of these regions. Most of the discovered galaxies would be the massive hosts of dormant quasars that left behind fossil HII cavities that persisted long after the quasar emission ended, owing to the long recombination time of intergalactic hydrogen. A snap-shot survey of candidate HII regions selected in redshifted 21cm image cubes may prove to be the most efficient method for finding very high redshift quasars and galaxies.

In the local universe, the masses of Super-Massive Black-Holes (SMBH) appear to correlate with physical properties of their hosts, including the mass of the dark-matter halos. At higher redshifts, we observe the growth of SMBHs indirectly through identification of high redshift quasars. However information on their hosts is difficult to obtain. In this paper we determine the masses of halos that host high redshift quasars (at z>4) by comparing the rate of growth of quasar density with that predicted by the Press-Schechter mass function. The host mass determined depends on how the ratio between SMBH and host halo mass evolves with redshift. Assuming the ratio between SMBH and halo mass does not evolve with redshift, we find a host halo mass of log(M)=11.7+/-0.3 solar masses. Even if the quasars shine at their Eddington limit, this mass is significantly smaller than that seen in the local universe. Indeed the null-hypothesis, of a constant SMBH to halo mass ratio at all redshifts, can be ruled out at greater than a 5-sigma level. SMBHs must therefore have contributed a larger fraction to the host mass in the past. Including redshift evolution of the SMBH to halo mass ratio, we find larger halo masses of log(M)=12.4+/-0.3 solar masses, and a ratio between SMBH and halo mass that increases with redshift as ~(1+z)^1.5 are required to be consistent with both local and high redshift observations. We investigate restrictions placed on the critical linear overdensity of quasar hosts at their epoch of virialisation and find that it cannot exceed the traditional value of delta_c=1.69 by more than a factor of two. Finally, we find that the high redshift quasars are hosted by fluctuations on scales that have a variance of (delta M/M) = 2-3, corresponding to (3-4.5)-sigma fluctuations in the density field.

We analyze the evolution of HII regions around the seven known SDSS quasars at z>6. The comparison between observed and model radii of the HII regions generated by these quasars individually, suggests that the surrounding intergalactic hydrogen is significantly neutral. When all constraints are combined, the existing quasar sample implies a volume averaged neutral fraction that is larger than 10% at z>6. This limited sample permits a preliminary analysis of the correlations between the quasar parameters, the sizes of their HII regions, and the associated constraints on the neutral hydrogen fraction. We find no evidence in these correlations to contradict the interpretation that the red side of the Gunn-Peterson trough corresponds to the boundary between an HII region and a partially neutral IGM.

The first galaxies to appear in the universe at redshifts z>20 created ionized bubbles in the intergalactic medium of neutral hydrogen left over from the Big-Bang. It is thought that the ionized bubbles grew with time, surrounded clusters of dwarf galaxies and eventually overlapped quickly throughout the universe over a narrow redshift interval near z~6. This event signaled the end of the reionization epoch when the universe was a billion years old. Measuring the hitherto unknown size distribution of the bubbles at their final overlap phase is a focus of forthcoming observational programs aimed at highly redshifted 21cm emission from atomic hydrogen. Here we show that the combined constraints of cosmic variance and causality imply an observed bubble size at the end of the overlap epoch of ~10 physical Mpc, and a scatter in the observed redshift of overlap along different lines-of-sight of ~0.15. This scatter is consistent with observational constraints from recent spectroscopic data on the farthest known quasars. Our novel result implies that future radio experiments should be tuned to a characteristic angular scale of ~0.5 degrees and have a minimum frequency band-width of ~8 MHz for an optimal detection of 21cm flux fluctuations near the end of reionization.

The discovery of Lyman-alpha emission from galaxies at redshifts beyond z~6.5 should not be naively interpreted as implying that the intergalactic medium (IGM) had been reionized at higher redsifts. We show that a cluster of faint undetected sources around each observed galaxy generates an HII region sufficiently large to allow transmission of the galaxy's Lyman-alpha line prior to reionization. We also show that quasars may contribute a significant fraction of the ionizing photons to HII regions around galaxies with a velocity dispersion larger than ~100km/s. These contributing quasars are not usually seen due to the small fraction of time they spend in a luminous phase.

The observed number counts of quasars may be explained either by long-lived activity within rare massive hosts, or by short-lived activity within smaller, more common hosts. It has been argued that quasar lifetimes may therefore be inferred from their clustering length, which determines the typical mass of the quasar host. Here we point out that the relationship between the mass of the black-hole and the circular velocity of its host dark-matter halo is more fundamental to the determination of the clustering length. In particular, the clustering length observed in the 2dF quasar redshift survey is consistent with the galactic halo - black-hole relation observed in local galaxies, provided that quasars shine at ~10-100% of their Eddington luminosity. The slow evolution of the clustering length with redshift inferred in the 2dF quasar survey favors a black-hole mass whose redshift-independent scaling is with halo circular velocity, rather than halo mass. These results are independent from observations of the number counts of bright quasars which may be used to determine the quasar lifetime and its dependence on redshift. We show that if quasar activity results from galaxy mergers, then the number counts of quasars imply an episodic quasar lifetime that is set by the dynamical time of the host galaxy rather than by the Salpeter time. Our results imply that as the redshift increases, the central black-holes comprise a larger fraction of their host galaxy mass and the quasar lifetime gets shorter.

The Ly-alpha absorption spectrum of the highest redshift quasars indicates that they are surrounded by giant HII regions, a few Mpc in size. The neutral gas around these HII regions should emit 21cm radiation in excess of the Cosmic Microwave Background, and enable future radio telescopes to measure the transverse extent of these HII regions. At early times, the HII regions expand with a relativistic speed. Consequently, their measured sizes along the line-of-sight (via Ly-alpha absorption) and transverse to it (via 21 cm emission) should have different observed values due to relativistic time-delay. We show that the combined measurement of these sizes would directly constrain the neutral fraction of the surrounding intergalactic medium (IGM) as well as the quasar lifetime. Based on current number counts of luminous quasars at z>6, an instrument like LOFAR should detect >2 redshifted 21cm shells per field (with a radius of 11 degrees) around active quasars as bright as those already discovered by SDSS, and >200 relic shells of inactive quasars per field. We show that Ly-alpha photons from the quasar are unable to heat the IGM or to couple the spin and kinetic temperatures of atomic hydrogen beyond the edge of the HII region. The detection of the IGM in 21cm emission around high redshift quasars would therefore gauge the presence of a cosmic Ly-alpha background during the reionization epoch.

The ionization fraction of cosmic hydrogen, left over from the big bang, provides crucial fossil evidence for when the first stars and quasar black holes formed in the infant universe. Spectra of the two most distant quasars known show nearly complete absorption of photons with wavelengths shorter than the Ly-alpha transition of neutral hydrogen, indicating that hydrogen in the intergalactic medium (IGM) had not been completely ionized at a redshift z~6.3, about a billion years after the big bang. Here we show that the radii of influence of ionizing radiation from these quasars imply that the surrounding IGM had a neutral hydrogen fraction of tens of percent prior to the quasar activity, much higher than previous lower limits of ~0.1%. When combined with the recent inference of a large cumulative optical depth to electron scattering after cosmological recombination from the WMAP data, our result suggests the existence of a second peak in the mean ionization history, potentially due to an early formation episode of the first stars.

The comoving mass density of massive black hole (MBH) remnants from pre-galactic star formation could have been similar in magnitude to the mass-density of supermassive black holes (SMBHs) in the present-day universe. We show that the fraction of MBHs that coalesce during the assembly of SMBHs can be extracted from the rate of ring-down gravitational waves that are detectable by Advanced LIGO. Based on the SMBH formation history inferred from the evolution of the quasar luminosity function, we show that an observed event rate of 1 per year will constrain the SMBH mass fraction that was contributed by MBHs coalescence down to a level of ~10^-6 for 20 solar mass MBH remnants (or ~10^-4 for 260 solar mass remnants).

We place a limit on the logarithmic slope of the luminous quasar luminosity function at z~6 of beta>-3.0 (90%) using gravitational lensing constraints to build on the limit of beta>-3.3 (90%) derived from an analysis of the luminosity distribution (Fan et al. 2003). This tight constraint is obtained by noting that of the two quasars which are lensed by foreground galaxies, neither are multiply imaged. These observations are surprising if the luminosity function is steep because magnification bias results in an overabundance of multiply imaged relative to singly imaged lensed quasars. Our Bayesian analysis uses the a-posteriori information regarding alignments with foreground galaxies of the two lensed quasars, and provides a constraint on beta that is nearly independent of the uncertain evolution in the lens population. The results suggest that the bright end of the quasar luminosity function continues to flatten out to z~6, as is observed between z~3 and z~5 (Fan et al. 2001). Provided that SDSS J1148-5251 at z=6.37 is magnified by an intervening lens galaxy at z~5 (White et al. 2003), we also show that the high lens redshift in this system implies a co-moving density of massive galaxies that is close to constant out to high redshift. This is in agreement with the lack of redshift evolution in the velocity function of dark-matter halos with velocity dispersions near 200 km/sec as predicted by the Press-Schechter formalism. The combination of constraints on the quasar luminosity function and lens galaxy evolution are used to compute an improved estimate for the z~6 multiple image lens fraction of ~1-3%.

The WMAP satellite has measured a large optical depth to electron scattering after cosmological recombination of 0.17+-0.04, implying significant reionization of the primordial gas only ~200 million years after the big bang. However, the most recent overlap of intergalactic HII regions must have occured at z<9 based on the Lyman-alpha forest constraint on the thermal history of the intergalactic medium. Here we argue that a first generation of metal-free stars with a heavy (rather than Salpeter) mass function is therefore required to account for much of the inferred optical depth. This conclusion holds if feedback regulates star formation in early dwarf galaxies as observed in present-day dwarfs.

We postulate that supermassive black-holes grow in the centers of galaxies until they unbind the galactic gas that feeds them. We show that the corresponding self-regulation condition yields a correlation between black-hole mass (Mbh) and galaxy velocity dispersion (sigma) as inferred in the local universe, and recovers the observed optical and X-ray luminosity functions of quasars at redshifts up to z~6 based on the hierarchical evolution of galaxy halos in a Lambda-CDM cosmology. With only one free parameter and a simple algorithm, our model yields the observed evolution in the number density of optically bright or X-ray faint quasars between 2<z<6 across 3 orders of magnitude in bolometric luminosity and 3 orders of magnitude in comoving density per logarithm of luminosity. The self-regulation condition identifies the dynamical time of galactic disks during the epoch of peak quasar activity (z~2.5) as the origin of the inferred characteristic quasar lifetime of ~10 million years. Since the lifetime becomes comparable to the Salpeter e-folding time at this epoch, the model also implies that the Mbh-sigma relation is a product of feedback regulated accretion during the peak of quasar activity. The mass-density in black-holes accreted by that time is consistent with the local black-hole mass density of ~(0.8-6.3) times 10^5 solar masses per cubic Mpc, which we have computed by combining the Mbh-sigma relation with the measured velocity dispersion function of SDSS galaxies (Sheth et al.~2003). Applying a similar self-regulation principle to supernova-driven winds from starbursts, we find that the ratio between the black hole mass and the stellar mass of galactic spheroids increases with redshift as (1+z)^1.5 although the Mbh-sigma relation is redshift-independent.

We compute the reionization histories of hydrogen and helium due to the ionizing radiation fields produced by stars and quasars. For the quasars we use a model based on halo-merger rates that reproduces all known properties of the quasar luminosity function at high redshifts. The less constrained properties of the ionizing radiation produced by stars are modeled with two free parameters: (i) a transition redshift, z_tran, above which the stellar population is dominated by massive, zero-metallicity stars and below which it is dominated by a Scalo mass function; (ii) the product of the escape fraction of stellar ionizing photons from their host galaxies and the star-formation efficiency, f_esc f_*. We constrain the allowed range of these free parameters at high redshifts based on the lack of the HI Gunn-Peterson trough at z<6 and the upper limit on the total intergalactic optical depth for electron scattering, tau_es<0.18, from recent cosmic microwave background (CMB) experiments. We find that quasars ionize helium by a redshift z~4, but cannot reionize hydrogen by themselves before z~6. A major fraction of the allowed combinations of f_esc f_* and z_tran lead to an early peak in the ionized fraction due to metal-free stars at high redshifts. This sometimes results in two reionization epochs, namely an early HII or HeIII overlap phase followed by recombination and a second overlap phase. Even if early overlap is not achieved, the peak in the visibility function for scattering of the CMB often coincides with the early ionization phase rather than with the actual reionization epoch. Consequently, tau_es does not correspond directly to the reionization redshift. We generically find values of tau_es>7%, that should be detectable by the MAP satellite.

We provide a simple theoretical model for the quasar luminosity function at high redshifts that naturally reproduces the statistical properties of the luminous SDSS quasar sample at redshifts z~4.3 and z>5.7. Our model is based on the assumptions that quasar emission is triggered by galaxy mergers, and that the black hole mass is proportional to a power-law in the circular velocity v of the host galactic halo. We assume that quasars shine at their Eddington luminosity over a time proportional to the mass ratio between the small and final galaxies in the merger. This simple model fits the quasar luminosity function at z~2-3, reproduces the normalization and logarithmic slope (beta -2.58) at z~4.3, explains the space density of bright SDSS quasars at z~6.0, reproduces the black hole - halo mass relation for dormant black holes in the local universe, and matches the estimated duty cycle of quasar activity (~10^7 years) in Lyman-break galaxies at z~3. Based on the derived luminosity function we predict the resulting gravitational lensing rates for high redshift quasars. The lens fractions in the SDSS samples are predicted to be ~2% at z~4.3 and ~10% at z>5.7. Interestingly, the limiting quasar luminosity in our best-fit relation of L proportional to v^5/G, scales as the binding energy of the host galaxy divided by its dynamical time, implying that feedback is the mechanism that regulates black hole growth in galactic potential wells.

We predict the effects of gravitational lensing on the color-selected flux-limited samples of z~4.3 and z>5.8 quasars, recently published by the Sloan Digital Sky Survey (SDSS). Our main findings are: (i) The lensing probability should be 1-2 orders of magnitude higher than for conventional surveys. The expected fraction of multiply-imaged quasars is highly sensitive to redshift and the uncertain slope of the bright end of the luminosity function, beta_h. For beta_h=2.58 (3.43) we find that at z~4.3 and i*<20.0 the fraction is ~4% (13%) while at z~6 and z*<20.2 the fraction is ~7% (30%). (ii) The distribution of magnifications is heavily skewed; sources having the redshift and luminosity of the SDSS z>5.8 quasars acquire median magnifications of med(mu_obs)~1.1-1.3 and mean magnifications of <mu_obs>~5-50. Estimates of the quasar luminosity density at high redshift must therefore filter out gravitationally-lensed sources. (iii) The flux in the Gunn-Peterson trough of the highest redshift (z=6.28) quasar is known to be f_lambda<3 10^-19 erg/sec/cm^2/Angstrom. Should this quasar be multiply imaged, we estimate a 40% chance that light from the lens galaxy would have contaminated the same part of the quasar spectrum with a higher flux. Hence, spectroscopic studies of the epoch of reionization need to account for the possibility that a lens galaxy, which boosts the quasar flux, also contaminates the Gunn-Peterson trough. (iv) Microlensing by stars should result in ~1/3 of multiply imaged quasars in the z>5.8 catalog varying by more than 0.5 magnitudes over the next decade. The median equivalent width would be lowered by ~20% with respect to the intrinsic value due to differential magnification of the continuum and emission-line regions.

Exceptionally bright quasars with redshifts up to z=6.28 have recently been discovered. Quasars are thought to be powered by the accretion of gas onto supermassive black holes at the centres of galaxies. Their maximum (Eddington) luminosity is proportional to the mass of the black hole, and so these bright quasars are inferred to have black holes with masses of more than a few billion solar masses. The existence of such massive black holes poses a challenge to models for the formation of structures in the early Universe, as it requires that the black holes would grow so massive in less than a billion years after the Big Bang. Here we show that up to a third of known quasars with z~6 will have their observed flux magnified by a factor of 10 or more through gravitational lensing by galaxies along the line of sight. The inferred abundance of quasar host galaxies, as well as the luminosity density provided by the quasars, are therefore substantially overestimated.

Intensive monitoring campaigns have recently attempted to measure the time delays between multiple images of gravitational lenses. Some of the resulting light-curves show puzzling low-level, rapid variability which is unique to individual images, superimposed on top of (and concurrent with) longer time-scale intrinsic quasar variations which repeat in all images. We demonstrate that both the amplitude and variability time-scale of the rapid light-curve anomalies, as well as the correlation observed between intrinsic and microlensed variability, are naturally explained by stellar microlensing of a smooth accretion disk which is occulted by optically-thick broad-line clouds. The rapid time-scale is caused by the high velocities of the clouds (~5x10^3 km/s), and the low amplitude results from the large number of clouds covering the magnified or demagnified parts of the disk. The observed amplitudes of variations in specific lenses implies that the number of broad-line clouds that cover ~10% of the quasar sky is ~10^5 per 4 pi steradian. This is comparable to the expected number of broad line clouds in models where the clouds originate from bloated stars.